| 1 |
loizides |
1.3 |
# $Id:$
|
| 2 |
|
|
|
| 3 |
bendavid |
1.1 |
process Rec =
|
| 4 |
|
|
{
|
| 5 |
paus |
1.2 |
untracked PSet maxEvents = { untracked int32 input = 10 }
|
| 6 |
loizides |
1.3 |
|
| 7 |
bendavid |
1.1 |
include "Configuration/StandardSequences/data/Services.cff"
|
| 8 |
|
|
include "Configuration/StandardSequences/data/Geometry.cff"
|
| 9 |
|
|
include "Configuration/StandardSequences/data/MagneticField.cff"
|
| 10 |
|
|
include "Configuration/StandardSequences/data/FakeConditions.cff"
|
| 11 |
|
|
include "FWCore/MessageService/data/MessageLogger.cfi"
|
| 12 |
|
|
|
| 13 |
paus |
1.2 |
# Sample generator configuration
|
| 14 |
loizides |
1.3 |
source = PythiaSource
|
| 15 |
|
|
{
|
| 16 |
|
|
untracked int32 pythiaPylistVerbosity = 1
|
| 17 |
|
|
untracked bool pythiaHepMCVerbosity = false
|
| 18 |
|
|
untracked int32 maxEventsToPrint = 3
|
| 19 |
|
|
# put here the cross section of your process (in pb)
|
| 20 |
|
|
untracked double crossSection = 0.00212
|
| 21 |
|
|
# put here the efficiency of your filter (1. if no filter)
|
| 22 |
|
|
untracked double filterEfficiency = 1.
|
| 23 |
|
|
PSet PythiaParameters = {
|
| 24 |
|
|
vstring parameterSets = {
|
| 25 |
|
|
"processParameters"
|
| 26 |
|
|
}
|
| 27 |
bendavid |
1.1 |
|
| 28 |
loizides |
1.3 |
include "Configuration/Generator/data/PythiaUESettings.cfi"
|
| 29 |
bendavid |
1.1 |
|
| 30 |
loizides |
1.3 |
vstring processParameters = {
|
| 31 |
|
|
"PMAS(25,1)=180.0 !mass of Higgs",
|
| 32 |
|
|
"MSEL=0 ! user selection for process",
|
| 33 |
|
|
# "MSTJ(41)=1 !Switch off Pythia QED bremsshtrahlung",
|
| 34 |
|
|
"MSUB(102)=1 !ggH",
|
| 35 |
|
|
"MSUB(123)=1 !ZZ fusion to H",
|
| 36 |
|
|
"MSUB(124)=1 !WW fusion to H",
|
| 37 |
|
|
"CKIN(45)=5. !high mass cut on m2 in 2 to 2 process",
|
| 38 |
|
|
"CKIN(46)=150. !high mass cut on secondary resonance m1 in 2->1->2 process",
|
| 39 |
|
|
"CKIN(47)=5. !low mass cut on secondary resonance m2 in 2->1->2 process",
|
| 40 |
|
|
"CKIN(48)=150. !high mass cut on secondary resonance m2 in 2->1->2 process",
|
| 41 |
|
|
# Z decays
|
| 42 |
|
|
"MDME(174,1)=0 !Z decay into d dbar",
|
| 43 |
|
|
"MDME(175,1)=0 !Z decay into u ubar",
|
| 44 |
|
|
"MDME(176,1)=0 !Z decay into s sbar",
|
| 45 |
|
|
"MDME(177,1)=0 !Z decay into c cbar",
|
| 46 |
|
|
"MDME(178,1)=0 !Z decay into b bbar",
|
| 47 |
|
|
"MDME(179,1)=0 !Z decay into t tbar",
|
| 48 |
|
|
"MDME(182,1)=1 !Z decay into e- e+",
|
| 49 |
|
|
"MDME(183,1)=0 !Z decay into nu_e nu_ebar",
|
| 50 |
|
|
"MDME(184,1)=1 !Z decay into mu- mu+",
|
| 51 |
|
|
"MDME(185,1)=0 !Z decay into nu_mu nu_mubar",
|
| 52 |
|
|
"MDME(186,1)=1 !Z decay into tau- tau+",
|
| 53 |
|
|
"MDME(187,1)=0 !Z decay into nu_tau nu_taubar",
|
| 54 |
|
|
# Higgs decays
|
| 55 |
|
|
"MDME(210,1)=0 !Higgs decay into dd",
|
| 56 |
|
|
"MDME(211,1)=0 !Higgs decay into uu",
|
| 57 |
|
|
"MDME(212,1)=0 !Higgs decay into ss",
|
| 58 |
|
|
"MDME(213,1)=0 !Higgs decay into cc",
|
| 59 |
|
|
"MDME(214,1)=0 !Higgs decay into bb",
|
| 60 |
|
|
"MDME(215,1)=0 !Higgs decay into tt",
|
| 61 |
|
|
"MDME(216,1)=0 !Higgs decay into",
|
| 62 |
|
|
"MDME(217,1)=0 !Higgs decay into Higgs decay",
|
| 63 |
|
|
"MDME(218,1)=0 !Higgs decay into e nu e",
|
| 64 |
|
|
"MDME(219,1)=0 !Higgs decay into mu nu mu",
|
| 65 |
|
|
"MDME(220,1)=0 !Higgs decay into tau nu tau",
|
| 66 |
|
|
"MDME(221,1)=0 !Higgs decay into Higgs decay",
|
| 67 |
|
|
"MDME(222,1)=0 !Higgs decay into g g",
|
| 68 |
|
|
"MDME(223,1)=0 !Higgs decay into gam gam",
|
| 69 |
|
|
"MDME(224,1)=0 !Higgs decay into gam Z",
|
| 70 |
|
|
"MDME(225,1)=1 !Higgs decay into Z Z",
|
| 71 |
|
|
"MDME(226,1)=0 !Higgs decay into W W"
|
| 72 |
|
|
}
|
| 73 |
|
|
}
|
| 74 |
bendavid |
1.1 |
}
|
| 75 |
|
|
|
| 76 |
|
|
# load generator sequence (VtxSmeared is needed inside, missing dependence)
|
| 77 |
|
|
include "Configuration/StandardSequences/data/VtxSmearedBetafuncEarlyCollision.cff"
|
| 78 |
|
|
include "Configuration/StandardSequences/data/Generator.cff"
|
| 79 |
|
|
|
| 80 |
|
|
include "Configuration/StandardSequences/data/Simulation.cff"
|
| 81 |
|
|
include "Configuration/StandardSequences/data/MixingNoPileUp.cff"
|
| 82 |
|
|
include "Configuration/StandardSequences/data/L1Emulator.cff"
|
| 83 |
|
|
include "Configuration/StandardSequences/data/DigiToRaw.cff"
|
| 84 |
|
|
include "Configuration/StandardSequences/data/VtxSmearedBetafuncEarlyCollision.cff"
|
| 85 |
|
|
include "Configuration/StandardSequences/data/Reconstruction.cff"
|
| 86 |
|
|
#include "Configuration/StandardSequences/data/UseZeroBField.cff"
|
| 87 |
|
|
|
| 88 |
|
|
# define the tree service
|
| 89 |
|
|
service = TreeService {
|
| 90 |
|
|
untracked vstring treeNames = { "MitTree" }
|
| 91 |
|
|
untracked vstring fileNames = { "mit-full" }
|
| 92 |
|
|
untracked vstring pathNames = { "." }
|
| 93 |
|
|
untracked vuint32 maxSizes = { 1024 }
|
| 94 |
|
|
untracked vuint32 compLevels = { 9 }
|
| 95 |
|
|
untracked vuint32 splitLevels = { 99 }
|
| 96 |
|
|
untracked vuint32 brSizes = { 32000 }
|
| 97 |
|
|
}
|
| 98 |
|
|
|
| 99 |
|
|
# customize the MIT filler
|
| 100 |
|
|
module MitGenFiller = FillGenParts {
|
| 101 |
|
|
}
|
| 102 |
|
|
|
| 103 |
|
|
module MitTrackFiller = FillTracks {
|
| 104 |
|
|
}
|
| 105 |
|
|
|
| 106 |
|
|
module MitMuonFiller = FillMuons {
|
| 107 |
|
|
}
|
| 108 |
|
|
|
| 109 |
|
|
module MitElectronFiller = FillElectrons {
|
| 110 |
|
|
}
|
| 111 |
|
|
|
| 112 |
|
|
# also make Edm output for the events we generate
|
| 113 |
|
|
include "Configuration/EventContent/data/EventContent.cff"
|
| 114 |
|
|
module FEVT = PoolOutputModule
|
| 115 |
|
|
{
|
| 116 |
|
|
using FEVTSIMEventContent
|
| 117 |
|
|
untracked string fileName = "edm-full.root"
|
| 118 |
|
|
}
|
| 119 |
|
|
|
| 120 |
|
|
# output path for the Edm file
|
| 121 |
|
|
endpath outpath = { FEVT }
|
| 122 |
|
|
|
| 123 |
|
|
|
| 124 |
paus |
1.2 |
path p0 = { pgen} # generator
|
| 125 |
|
|
path p1 = { psim} # simulation
|
| 126 |
|
|
path p2 = { pdigi} # digitization
|
| 127 |
|
|
path p3 = { L1Emulator } # trigger: level-1 emulation
|
| 128 |
|
|
path p4 = { DigiToRaw } # converting digitized to raw data
|
| 129 |
|
|
path p5 = { reconstruction } # reconstruction
|
| 130 |
|
|
path p6 = { MitGenFiller,MitTrackFiller,MitMuonFiller,MitElectronFiller } # MitProd tree filler
|
| 131 |
bendavid |
1.1 |
|
| 132 |
|
|
schedule = {p0,p1,p2,p3,p4,p5,p6,outpath}
|
| 133 |
|
|
}
|
