CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
TTbar_Validation_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 TTbarAnalyzeSpinCorr = cms.EDAnalyzer("TTbarSpinCorrHepMCAnalyzer",
4  genEventInfoProductTag = cms.InputTag("generator"),
5  genParticlesTag = cms.InputTag("genParticles")
6  )
7 
9 lheCOMWeightProducer.NewECMS = cms.double(8000)
10 
11 
12 ## get lorentzvectors
13 analyzeTopKinematics = cms.EDAnalyzer('TTbar_Kinematics',
14  SaveTree = cms.untracked.bool(False),
15  hepmcCollection = cms.InputTag("generatorSmeared"),
16  genEventInfoProductTag = cms.InputTag("generator")
17  )
18 
19 ## analyze genjets
20 analyzeGenJets = cms.EDAnalyzer("TTbar_GenJetAnalyzer",
21  jets = cms.InputTag('ak4GenJets' ),
22  genEventInfoProductTag = cms.InputTag("generator")
23  )
24 
25 # --- Create list of interesting genParticles ---
27 
28 # --- ShortList is to prevent running multiple times over full particle list ---
29 #see https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideGenParticlePruner
30 
31 genParticlesShortList = cms.EDProducer("GenParticlePruner",
32  src = cms.InputTag("genParticles",""),
33  select = cms.vstring("drop * ",
34  "keep pdgId = {mu+} & status = 1",
35  "keep pdgId = {mu-} & status = 1",
36  "keep pdgId = {e+} & status = 1",
37  "keep pdgId = {e-} & status = 1",
38  "keep pdgId = {nu_e} & status = 1",
39  "keep pdgId = {nu_ebar} & status = 1",
40  "keep pdgId = {nu_mu} & status = 1",
41  "keep pdgId = {nu_mubar} & status = 1"
42  )
43  )
44 
45 #see https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideGenParticlePruner
46 genParticlesMuons = cms.EDProducer("GenParticlePruner",
47  src = cms.InputTag("genParticlesShortList"),
48  select = cms.vstring("drop * ",
49  "keep pdgId = {mu+} & status = 1 & pt > 30 & abs(eta) < 2.5",
50  "keep pdgId = {mu-} & status = 1 & pt > 30 & abs(eta) < 2.5"
51  )
52  )
53 
54 #see https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideGenParticlePruner
55 genParticlesElectrons = cms.EDProducer("GenParticlePruner",
56  src = cms.InputTag("genParticlesShortList"),
57  select = cms.vstring("drop * ",
58  "keep pdgId = {e+} & status = 1 & pt > 30 & abs(eta) < 2.5",
59  "keep pdgId = {e-} & status = 1 & pt > 30 & abs(eta) < 2.5"
60  )
61  )
62 
63 #see https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideGenParticlePruner
64 genParticlesNeutrinos = cms.EDProducer("GenParticlePruner",
65  src = cms.InputTag("genParticlesShortList"),
66  select = cms.vstring("drop * ",
67  "keep pdgId = {nu_e} & status = 1",
68  "keep pdgId = {nu_ebar} & status = 1",
69  "keep pdgId = {nu_mu} & status = 1",
70  "keep pdgId = {nu_mubar} & status = 1"
71  )
72  )
73 
74 ## analyze gen leptons
75 analyzeGenMuons = cms.EDAnalyzer("TTbar_GenLepAnalyzer", leptons = cms.InputTag('genParticlesMuons' ))
76 analyzeGenElecs = cms.EDAnalyzer("TTbar_GenLepAnalyzer", leptons = cms.InputTag('genParticlesElectrons' ))
77 analyzeGenNtrns = cms.EDAnalyzer("TTbar_GenLepAnalyzer", leptons = cms.InputTag('genParticlesNeutrinos' ))
78 
79 
80