Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-11-03 19:40:29

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0003 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0004 
0005 generator = cms.EDFilter("Pythia8GeneratorFilter",
0006     pythiaPylistVerbosity = cms.untracked.int32(0),
0007     pythiaHepMCVerbosity = cms.untracked.bool(False),
0008     maxEventsToPrint = cms.untracked.int32(0),
0009     comEnergy = cms.double(14000.0),
0010     ExternalDecays = cms.PSet(
0011         EvtGen130 = cms.untracked.PSet(
0012             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2014_NOLONGLIFE.DEC'),
0013             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_2014.pdl'),
0014             operates_on_particles = cms.vint32(20443,445),                       # we care just about our signal particles
0015             convertPythiaCodes = cms.untracked.bool(False),
0016             user_decay_file = cms.vstring('GeneratorInterface/ExternalDecays/data/Onia_chic_jpsigamma.dec'),
0017             list_forced_decays = cms.vstring('Mychi_c1','Mychi_c2'),
0018         ),
0019         parameterSets = cms.vstring('EvtGen130')
0020     ),
0021     PythiaParameters = cms.PSet(
0022         pythia8CommonSettingsBlock,
0023         pythia8CP5SettingsBlock,
0024         processParameters = cms.vstring(
0025     # generate just the needed and nothing else
0026             'Charmonium:states(3PJ) = 20443,445',
0027             'Charmonium:O(3PJ)[3P0(1)] = 0.05,0.05',
0028             'Charmonium:O(3PJ)[3S1(8)] = 0.0031,0.0031',
0029             'Charmonium:gg2ccbar(3PJ)[3PJ(1)]g = on,on',
0030             'Charmonium:qg2ccbar(3PJ)[3PJ(1)]q = on,on',
0031             'Charmonium:qqbar2ccbar(3PJ)[3PJ(1)]g = on,on',
0032             'Charmonium:gg2ccbar(3PJ)[3S1(8)]g = on,on',
0033             'Charmonium:qg2ccbar(3PJ)[3S1(8)]q = on,on',
0034             'Charmonium:qqbar2ccbar(3PJ)[3S1(8)]g = on,on',
0035 #
0036             'PhaseSpace:pTHatMin = 10.'                   # (filter efficiency 4.671e-02) be aware of this ckin(3) equivalent
0037             ),
0038         parameterSets = cms.vstring('pythia8CommonSettings',
0039                                     'pythia8CP5Settings',
0040                                     'processParameters',
0041                                     )
0042     )
0043 )
0044 
0045 # Next two muon filter are derived from muon reconstruction
0046 oniafilter = cms.EDFilter("MCParticlePairFilter",
0047     Status = cms.untracked.vint32(2, 1),
0048     MinPt = cms.untracked.vdouble(16., 0.2),
0049     MaxEta = cms.untracked.vdouble(1.6, 1.6),
0050     MinEta = cms.untracked.vdouble(-1.6, -1.6),
0051     ParticleCharge = cms.untracked.int32(0),
0052     MinP = cms.untracked.vdouble(0.,0.),
0053     ParticleID1 = cms.untracked.vint32(443),
0054     ParticleID2 = cms.untracked.vint32(22)
0055 )
0056 
0057 muminusfilter = cms.EDFilter("PythiaDauVFilter",
0058     MotherID = cms.untracked.int32(0),
0059     MinPt = cms.untracked.vdouble(2.5, 2.5, 3.5),
0060     ParticleID = cms.untracked.int32(443),
0061     ChargeConjugation = cms.untracked.bool(False),
0062     MinEta = cms.untracked.vdouble(1.2, -1.6, -1.2),
0063     MaxEta = cms.untracked.vdouble(1.6, -1.2, 1.2),
0064     NumberDaughters = cms.untracked.int32(1),
0065     DaughterIDs = cms.untracked.vint32(-13, -13, -13)
0066 )
0067 
0068 muplusfilter = cms.EDFilter("PythiaDauVFilter",
0069     MotherID = cms.untracked.int32(0),
0070     MinPt = cms.untracked.vdouble(2.5, 2.5, 3.5),
0071     ParticleID = cms.untracked.int32(443),
0072     ChargeConjugation = cms.untracked.bool(False),
0073     MinEta = cms.untracked.vdouble(1.2, -1.6, -1.2),
0074     MaxEta = cms.untracked.vdouble(1.6, -1.2, 1.2),
0075     NumberDaughters = cms.untracked.int32(1),
0076     DaughterIDs = cms.untracked.vint32(13, 13, 13)
0077 )
0078 
0079 ProductionFilterSequence = cms.Sequence(generator*oniafilter*muminusfilter*muplusfilter)