/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/Validation/EventGenerator/python/TTbar_Validation_cfi.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

TTbarAnalyzeSpinCorr = cms.EDAnalyzer("TTbarSpinCorrHepMCAnalyzer",
                                      genEventInfoProductTag = cms.InputTag("generator"),
                                      genParticlesTag = cms.InputTag("genParticles")
                                      )

from GeneratorInterface.LHEInterface.lheCOMWeightProducer import *
lheCOMWeightProducer.NewECMS = cms.double(8000)

## produce TtGenEvt
from TopQuarkAnalysis.TopEventProducers.sequences.ttGenEvent_cff import *
decaySubset.fillMode = "kME"
decaySubset.addRadiation = False

## get lorentzvectors
analyzeTopKinematics = cms.EDAnalyzer('TTbar_Kinematics',
                                      SaveTree = cms.untracked.bool(False),
                                      genEventInfoProductTag = cms.InputTag("generator")
                                      )

## analyze genjets
analyzeGenJets = cms.EDAnalyzer("TTbar_GenJetAnalyzer",
                                jets = cms.InputTag('ak5GenJets' ),
                                genEventInfoProductTag = cms.InputTag("generator")
                                )

# --- Create list of interesting genParticles ---
from SimGeneral.HepPDTESSource.pythiapdt_cfi import *

# --- ShortList is to prevent running multiple times over full particle list ---
#see https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideGenParticlePruner

genParticlesShortList = cms.EDProducer("GenParticlePruner",
                                       src = cms.InputTag("genParticles",""),
                                       select = cms.vstring("drop  *  ",
                                                            "keep pdgId = {mu+}      & status = 1",
                                                            "keep pdgId = {mu-}      & status = 1",
                                                            "keep pdgId = {e+}       & status = 1",
                                                            "keep pdgId = {e-}       & status = 1",
                                                            "keep pdgId = {nu_e}     & status = 1",
                                                            "keep pdgId = {nu_ebar}  & status = 1",
                                                            "keep pdgId = {nu_mu}    & status = 1",
                                                            "keep pdgId = {nu_mubar} & status = 1"
                                                            )
                                       )

#see https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideGenParticlePruner
genParticlesMuons     = cms.EDProducer("GenParticlePruner",
                                       src = cms.InputTag("genParticlesShortList"),
                                       select = cms.vstring("drop  *  ",
                                                            "keep pdgId = {mu+} & status = 1 & pt > 30 & abs(eta) < 2.5",
                                                            "keep pdgId = {mu-} & status = 1 & pt > 30 & abs(eta) < 2.5"
                                                            )
                                       )

#see https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideGenParticlePruner
genParticlesElectrons = cms.EDProducer("GenParticlePruner",
                                       src = cms.InputTag("genParticlesShortList"),
                                       select = cms.vstring("drop  *  ",
                                                            "keep pdgId = {e+} & status = 1 & pt > 30 & abs(eta) < 2.5",
                                                            "keep pdgId = {e-} & status = 1 & pt > 30 & abs(eta) < 2.5"
                                                            )
                                       )

#see https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideGenParticlePruner
genParticlesNeutrinos = cms.EDProducer("GenParticlePruner",
                                       src = cms.InputTag("genParticlesShortList"),
                                       select = cms.vstring("drop  *  ",
                                                            "keep pdgId = {nu_e}     & status = 1",
                                                            "keep pdgId = {nu_ebar}  & status = 1",
                                                            "keep pdgId = {nu_mu}    & status = 1",
                                                            "keep pdgId = {nu_mubar} & status = 1"
                                                            )
                                       )

## analyze gen leptons
analyzeGenMuons = cms.EDAnalyzer("TTbar_GenLepAnalyzer", leptons = cms.InputTag('genParticlesMuons' ))
analyzeGenElecs = cms.EDAnalyzer("TTbar_GenLepAnalyzer", leptons = cms.InputTag('genParticlesElectrons' ))
analyzeGenNtrns = cms.EDAnalyzer("TTbar_GenLepAnalyzer", leptons = cms.InputTag('genParticlesNeutrinos' ))