TreeFiller/test/FullChainExample.cfg

process Rec = 
{
    untracked PSet maxEvents = {untracked int32 input = 10}
    # this example configuration offers some minimum 
    # annotation, to help users get through; please
    # don't hesitate to read through the comments
    
    
    # use MessageLogger to redirect/suppress multiple
    # service messages coming from the system
    #
    # in this config below, we use the replace option to make
    # the logger let out messages of severity ERROR (INFO level
    # will be suppressed), and we want to limit the number to 10
    #
    include "Configuration/StandardSequences/data/Services.cff"
    include "Configuration/StandardSequences/data/Geometry.cff"
    include "Configuration/StandardSequences/data/MagneticField.cff"
    include "Configuration/StandardSequences/data/FakeConditions.cff" 
    include "FWCore/MessageService/data/MessageLogger.cfi"
    
    # in you wish to quiet the logger even more and leave in 
    # only info from the Framework, use config example below :
    #
    #service = MessageLogger
    #{
    #   untracked vstring destinations = {"cout"} 
    #   untracked vstring categories = { "FwkJob" } 
    #   untracked PSet cout =
    #   {
    #      untracked PSet default = { untracked int32 limit = 0 }  # kill all messages in the log
    #      untracked PSet FwkJob  = { untracked int32 limit = -1 } # but FwkJob category - those unlimitted
    #   } 
    #}
    
    include "Configuration/StandardSequences/data/SimulationRandomNumberGeneratorSeeds.cff"
    include "SimGeneral/HepPDTESSource/data/pythiapdt.cfi"

    
  #Sample generator configuration  
  source = PythiaSource 
  { 
  untracked int32 pythiaPylistVerbosity = 1
  untracked bool pythiaHepMCVerbosity = false
  untracked int32 maxEventsToPrint = 3
  
  # put here the cross section of your process (in pb)
  untracked double crossSection = 0.00212
  # put here the efficiency of your filter (1. if no filter)
  untracked double filterEfficiency = 1.


  PSet PythiaParameters = {

    # This is a vector of ParameterSet names to be read, in this order
    vstring parameterSets = {
      "processParameters"
    }

   include "Configuration/Generator/data/PythiaUESettings.cfi"
 
       vstring processParameters = {
           "PMAS(25,1)=180.0        !mass of Higgs",
           "MSEL=0                  ! user selection for process",
#          "MSTJ(41)=1              !Switch off Pythia QED bremsshtrahlung",
           "MSUB(102)=1             !ggH",
           "MSUB(123)=1             !ZZ fusion to H",
           "MSUB(124)=1             !WW fusion to H",
           "CKIN(45)=5.             !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch",
           "CKIN(46)=150.           !high mass cut on secondary resonance m1 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch",
           "CKIN(47)=5.             !low mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch",
           "CKIN(48)=150.           !high mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch",
# Z decays
           "MDME(174,1)=0           !Z decay into d dbar",        
           "MDME(175,1)=0           !Z decay into u ubar",
           "MDME(176,1)=0           !Z decay into s sbar",
           "MDME(177,1)=0           !Z decay into c cbar",
           "MDME(178,1)=0           !Z decay into b bbar",
           "MDME(179,1)=0           !Z decay into t tbar",
           "MDME(182,1)=1           !Z decay into e- e+",
           "MDME(183,1)=0           !Z decay into nu_e nu_ebar",
           "MDME(184,1)=1           !Z decay into mu- mu+",
           "MDME(185,1)=0           !Z decay into nu_mu nu_mubar",
           "MDME(186,1)=1           !Z decay into tau- tau+",
           "MDME(187,1)=0           !Z decay into nu_tau nu_taubar",
# Higgs decays
           "MDME(210,1)=0           !Higgs decay into dd",
           "MDME(211,1)=0           !Higgs decay into uu",
           "MDME(212,1)=0           !Higgs decay into ss",
           "MDME(213,1)=0           !Higgs decay into cc",
           "MDME(214,1)=0           !Higgs decay into bb",
           "MDME(215,1)=0           !Higgs decay into tt",
           "MDME(216,1)=0           !Higgs decay into",
           "MDME(217,1)=0           !Higgs decay into Higgs decay",
           "MDME(218,1)=0           !Higgs decay into e nu e",
           "MDME(219,1)=0           !Higgs decay into mu nu mu",
           "MDME(220,1)=0           !Higgs decay into tau nu tau",
           "MDME(221,1)=0           !Higgs decay into Higgs decay",
           "MDME(222,1)=0           !Higgs decay into g g",
           "MDME(223,1)=0           !Higgs decay into gam gam",
           "MDME(224,1)=0           !Higgs decay into gam Z",
           "MDME(225,1)=1           !Higgs decay into Z Z",
           "MDME(226,1)=0           !Higgs decay into W W"
       }
    }
   }
    
    # load generator sequence (VtxSmeared is needed inside, missing dependence)
    include "Configuration/StandardSequences/data/VtxSmearedBetafuncEarlyCollision.cff"
    include "Configuration/StandardSequences/data/Generator.cff"
    
    # this config frament brings you 3 steps of the detector simulation:
    # -- vertex smearing (IR modeling)
    # -- G4-based hit level detector simulation
    # -- digitization (electronics readout modeling)
    # it returns 2 sequences : 
    # -- psim (vtx smearing + G4 sim)
    # -- pdigi (digitization in all subsystems, i.e. tracker=pix+sistrips,
    #           cal=ecal+ecal-0-suppression+hcal), muon=csc+dt+rpc)
    #
    include "Configuration/StandardSequences/data/Simulation.cff"
    
    # please note the IMPORTANT: 
    # in order to operate Digis, one needs to include Mixing module 
    # (pileup modeling), at least in the 0-pileup mode
    #
    # There're 3 possible configurations of the Mixing module :
    # no-pileup, low luminosity pileup, and high luminosity pileup
    #
    # they come, respectively, through the 3 config fragments below
    #
    # *each* config returns label "mix"; thus you canNOT have them
    # all together in the same configuration, but only one !!!
    #
    
    
    include "Configuration/StandardSequences/data/MixingNoPileUp.cff" 
    #include "Configuration/StandardSequences/data/MixingLowLumiPileUp.cff" 
    #include "Configuration/StandardSequences/data/MixingHighLumiPileUp.cff" 
    
    include "Configuration/StandardSequences/data/L1Emulator.cff" 

    include "Configuration/StandardSequences/data/DigiToRaw.cff"

    include "Configuration/StandardSequences/data/VtxSmearedBetafuncEarlyCollision.cff"
    
    include "Configuration/StandardSequences/data/Reconstruction.cff"
    #
    # You can set Mag Field to 0
    #
    #include "Configuration/StandardSequences/data/UseZeroBField.cff"
    
    # define the tree service
    service = TreeService { 
        untracked vstring treeNames   = { "MitTree" }
        untracked vstring fileNames   = { "mit-full" }
        untracked vstring pathNames   = { "." }
        untracked vuint32 maxSizes    = { 1024 }
        untracked vuint32 compLevels  = { 9 }
        untracked vuint32 splitLevels = { 99 }
        untracked vuint32 brSizes     = { 32000 }
    }
    
    # customize the MIT filler
    module MitGenFiller = FillGenParts {
    }
    
    module MitTrackFiller = FillTracks {
    }
    
    module MitMuonFiller = FillMuons {
    }
    
    module MitElectronFiller = FillElectrons {
    }
    
    
    # also make Edm output for the events we generate
    include "Configuration/EventContent/data/EventContent.cff"
    module FEVT = PoolOutputModule
    {
        using FEVTSIMEventContent
        untracked string fileName = "edm-full.root"
    }
    
    # output path for the Edm file
    endpath outpath = { FEVT }
    
    
    path p0 = {pgen} # generator
    path p1 = {psim} # simulation
    path p2 = {pdigi} # digitization
    path p3 = {L1Emulator}
    path p4 = {DigiToRaw}
    path p5 = {reconstruction} # reconstruction
    path p6 = { MitGenFiller, MitTrackFiller, MitMuonFiller, MitElectronFiller } #MitProd tree filler
    
    schedule = {p0,p1,p2,p3,p4,p5,p6,outpath}
    
    
}               
Revision:	1.1
Committed:	Mon Jun 9 00:06:36 2008 UTC (16 years, 11 months ago) by bendavid
Branch:	MAIN
Log Message:	adherance to coding conventions, reverted ParticleExample, new full chain example
#	User	Rev	Content
1	bendavid	1.1	process Rec =
2			{
3			untracked PSet maxEvents = {untracked int32 input = 10}
4			# this example configuration offers some minimum
5			# annotation, to help users get through; please
6			# don't hesitate to read through the comments
7
8
9			# use MessageLogger to redirect/suppress multiple
10			# service messages coming from the system
11			#
12			# in this config below, we use the replace option to make
13			# the logger let out messages of severity ERROR (INFO level
14			# will be suppressed), and we want to limit the number to 10
15			#
16			include "Configuration/StandardSequences/data/Services.cff"
17			include "Configuration/StandardSequences/data/Geometry.cff"
18			include "Configuration/StandardSequences/data/MagneticField.cff"
19			include "Configuration/StandardSequences/data/FakeConditions.cff"
20			include "FWCore/MessageService/data/MessageLogger.cfi"
21
22			# in you wish to quiet the logger even more and leave in
23			# only info from the Framework, use config example below :
24			#
25			#service = MessageLogger
26			#{
27			# untracked vstring destinations = {"cout"}
28			# untracked vstring categories = { "FwkJob" }
29			# untracked PSet cout =
30			# {
31			# untracked PSet default = { untracked int32 limit = 0 } # kill all messages in the log
32			# untracked PSet FwkJob = { untracked int32 limit = -1 } # but FwkJob category - those unlimitted
33			# }
34			#}
35
36			include "Configuration/StandardSequences/data/SimulationRandomNumberGeneratorSeeds.cff"
37			include "SimGeneral/HepPDTESSource/data/pythiapdt.cfi"
38
39
40			#Sample generator configuration
41			source = PythiaSource
42			{
43			untracked int32 pythiaPylistVerbosity = 1
44			untracked bool pythiaHepMCVerbosity = false
45			untracked int32 maxEventsToPrint = 3
46
47			# put here the cross section of your process (in pb)
48			untracked double crossSection = 0.00212
49			# put here the efficiency of your filter (1. if no filter)
50			untracked double filterEfficiency = 1.
51
52
53
54			PSet PythiaParameters = {
55
56			# This is a vector of ParameterSet names to be read, in this order
57			vstring parameterSets = {
58			"processParameters"
59			}
60
61			include "Configuration/Generator/data/PythiaUESettings.cfi"
62
63			vstring processParameters = {
64			"PMAS(25,1)=180.0 !mass of Higgs",
65			"MSEL=0 ! user selection for process",
66			# "MSTJ(41)=1 !Switch off Pythia QED bremsshtrahlung",
67			"MSUB(102)=1 !ggH",
68			"MSUB(123)=1 !ZZ fusion to H",
69			"MSUB(124)=1 !WW fusion to H",
70			"CKIN(45)=5. !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch",
71			"CKIN(46)=150. !high mass cut on secondary resonance m1 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch",
72			"CKIN(47)=5. !low mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch",
73			"CKIN(48)=150. !high mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch",
74			# Z decays
75			"MDME(174,1)=0 !Z decay into d dbar",
76			"MDME(175,1)=0 !Z decay into u ubar",
77			"MDME(176,1)=0 !Z decay into s sbar",
78			"MDME(177,1)=0 !Z decay into c cbar",
79			"MDME(178,1)=0 !Z decay into b bbar",
80			"MDME(179,1)=0 !Z decay into t tbar",
81			"MDME(182,1)=1 !Z decay into e- e+",
82			"MDME(183,1)=0 !Z decay into nu_e nu_ebar",
83			"MDME(184,1)=1 !Z decay into mu- mu+",
84			"MDME(185,1)=0 !Z decay into nu_mu nu_mubar",
85			"MDME(186,1)=1 !Z decay into tau- tau+",
86			"MDME(187,1)=0 !Z decay into nu_tau nu_taubar",
87			# Higgs decays
88			"MDME(210,1)=0 !Higgs decay into dd",
89			"MDME(211,1)=0 !Higgs decay into uu",
90			"MDME(212,1)=0 !Higgs decay into ss",
91			"MDME(213,1)=0 !Higgs decay into cc",
92			"MDME(214,1)=0 !Higgs decay into bb",
93			"MDME(215,1)=0 !Higgs decay into tt",
94			"MDME(216,1)=0 !Higgs decay into",
95			"MDME(217,1)=0 !Higgs decay into Higgs decay",
96			"MDME(218,1)=0 !Higgs decay into e nu e",
97			"MDME(219,1)=0 !Higgs decay into mu nu mu",
98			"MDME(220,1)=0 !Higgs decay into tau nu tau",
99			"MDME(221,1)=0 !Higgs decay into Higgs decay",
100			"MDME(222,1)=0 !Higgs decay into g g",
101			"MDME(223,1)=0 !Higgs decay into gam gam",
102			"MDME(224,1)=0 !Higgs decay into gam Z",
103			"MDME(225,1)=1 !Higgs decay into Z Z",
104			"MDME(226,1)=0 !Higgs decay into W W"
105			}
106			}
107			}
108
109			# load generator sequence (VtxSmeared is needed inside, missing dependence)
110			include "Configuration/StandardSequences/data/VtxSmearedBetafuncEarlyCollision.cff"
111			include "Configuration/StandardSequences/data/Generator.cff"
112
113			# this config frament brings you 3 steps of the detector simulation:
114			# -- vertex smearing (IR modeling)
115			# -- G4-based hit level detector simulation
116			# -- digitization (electronics readout modeling)
117			# it returns 2 sequences :
118			# -- psim (vtx smearing + G4 sim)
119			# -- pdigi (digitization in all subsystems, i.e. tracker=pix+sistrips,
120			# cal=ecal+ecal-0-suppression+hcal), muon=csc+dt+rpc)
121			#
122			include "Configuration/StandardSequences/data/Simulation.cff"
123
124			# please note the IMPORTANT:
125			# in order to operate Digis, one needs to include Mixing module
126			# (pileup modeling), at least in the 0-pileup mode
127			#
128			# There're 3 possible configurations of the Mixing module :
129			# no-pileup, low luminosity pileup, and high luminosity pileup
130			#
131			# they come, respectively, through the 3 config fragments below
132			#
133			# each config returns label "mix"; thus you canNOT have them
134			# all together in the same configuration, but only one !!!
135			#
136
137
138			include "Configuration/StandardSequences/data/MixingNoPileUp.cff"
139			#include "Configuration/StandardSequences/data/MixingLowLumiPileUp.cff"
140			#include "Configuration/StandardSequences/data/MixingHighLumiPileUp.cff"
141
142			include "Configuration/StandardSequences/data/L1Emulator.cff"
143
144			include "Configuration/StandardSequences/data/DigiToRaw.cff"
145
146			include "Configuration/StandardSequences/data/VtxSmearedBetafuncEarlyCollision.cff"
147
148			include "Configuration/StandardSequences/data/Reconstruction.cff"
149			#
150			# You can set Mag Field to 0
151			#
152			#include "Configuration/StandardSequences/data/UseZeroBField.cff"
153
154			# define the tree service
155			service = TreeService {
156			untracked vstring treeNames = { "MitTree" }
157			untracked vstring fileNames = { "mit-full" }
158			untracked vstring pathNames = { "." }
159			untracked vuint32 maxSizes = { 1024 }
160			untracked vuint32 compLevels = { 9 }
161			untracked vuint32 splitLevels = { 99 }
162			untracked vuint32 brSizes = { 32000 }
163			}
164
165			# customize the MIT filler
166			module MitGenFiller = FillGenParts {
167			}
168
169			module MitTrackFiller = FillTracks {
170			}
171
172			module MitMuonFiller = FillMuons {
173			}
174
175			module MitElectronFiller = FillElectrons {
176			}
177
178
179			# also make Edm output for the events we generate
180			include "Configuration/EventContent/data/EventContent.cff"
181			module FEVT = PoolOutputModule
182			{
183			using FEVTSIMEventContent
184			untracked string fileName = "edm-full.root"
185			}
186
187			# output path for the Edm file
188			endpath outpath = { FEVT }
189
190
191			path p0 = {pgen} # generator
192			path p1 = {psim} # simulation
193			path p2 = {pdigi} # digitization
194			path p3 = {L1Emulator}
195			path p4 = {DigiToRaw}
196			path p5 = {reconstruction} # reconstruction
197			path p6 = { MitGenFiller, MitTrackFiller, MitMuonFiller, MitElectronFiller } #MitProd tree filler
198
199			schedule = {p0,p1,p2,p3,p4,p5,p6,outpath}
200
201
202			}