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