process.goodOfflinePrimaryVertices = cms.EDFilter( "PrimaryVertexObjectFilter", filterParams = pvSelector.clone( minNdof = cms.double(4.0), maxZ = cms.double(24.0) ), src=cms.InputTag('offlinePrimaryVertices') ) #input for softLeptonTagInfos process.softPFElectronsTagInfos.primaryVertex = cms.InputTag('goodOfflinePrimaryVertices') process.softPFMuonsTagInfos.primaryVertex = cms.InputTag('goodOfflinePrimaryVertices') #for the flavour matching from PhysicsTools.JetMCAlgos.HadronAndPartonSelector_cfi import selectedHadronsAndPartons process.selectedHadronsAndPartons = selectedHadronsAndPartons.clone() from PhysicsTools.JetMCAlgos.AK5PFJetsMCFlavourInfos_cfi import ak5JetFlavourInfos process.jetFlavourInfosAK5PFJets = ak5JetFlavourInfos.clone() process.jetFlavourInfosAK5PFJets.jets = cms.InputTag("ak5PFJets") process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(11) ) process.source = cms.Source("PoolSource", skipEvents=cms.untracked.uint32(0), fileNames = cms.untracked.vstring( #'/store/mc/Fall11/QCD_Pt-50to80_Tune4C_7TeV_pythia8/AODSIM/PU_S6_START44_V9B-v1/0000/00225C27-4744-E111-8583-003048C692E4.root', #'/store/mc/Fall11/QCD_Pt-80to120_Tune4C_7TeV_pythia8/AODSIM/PU_S6_START44_V9B-v1/0000/00D668C4-5244-E111-A8F6-0025901D4936.root', #'/store/mc/Fall11/QCD_Pt-120to170_Tune4C_7TeV_pythia8/AODSIM/PU_S6_START44_V9B-v1/0000/023CE3AD-CF46-E111-90CB-003048C690A0.root', #'/store/mc/Fall11/QCD_Pt-170to300_Tune4C_7TeV_pythia8/AODSIM/PU_S6_START44_V9B-v1/0000/0433579C-7644-E111-9D93-00215AD4D6C8.root', #'/store/mc/Fall11/QCD_Pt-300to470_Tune4C_7TeV_pythia8/AODSIM/PU_S6_START44_V9B-v1/0000/00C128DE-923F-E111-AA76-003048C69312.root',
#fileNames = cms.untracked.vstring('root://eospublic.cern.ch//eos/opendata/cms/MonteCarlo2012/Summer12_DR53X/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball/AODSIM/PU_RD1_START53_V7N-v1/20000/003063B7-4CCF-E211-9FED-003048D46124.root')) # Set the maximum number of events to be processed (-1 processes all events) process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1)) # Number of events to be skipped (0 by default) process.source.skipEvents = cms.untracked.uint32(0) #------------------------------------------------------------------------- # AK5 jets #------------------------------------------------------------------------- from PhysicsTools.JetMCAlgos.HadronAndPartonSelector_cfi import selectedHadronsAndPartons process.selectedHadronsAndPartons = selectedHadronsAndPartons.clone() from PhysicsTools.JetMCAlgos.AK5PFJetsMCFlavourInfos_cfi import ak5JetFlavourInfos process.jetFlavourInfosAK5PFJets = ak5JetFlavourInfos.clone() ################################ ## MC Weights Analyzer ################################ process.mcweightanalyzer = cms.EDAnalyzer( "WeightAnalyzerBEff", JetsTag=cms.InputTag("ak5PFJets"), discriminator=cms.string("combinedSecondaryVertexBJetTags"), DiscriminatorValueTight=cms.double(0.898), DiscriminatorValueMedium=cms.double(0.679), DiscriminatorValueLoose=cms.double(0.244), ) process.TFileService = cms.Service( "TFileService", fileName=cms.string("flavortagefficiencies.root"))
import FWCore.ParameterSet.Config as cms # Flavour info: jet collection with all associated ghosts from PhysicsTools.JetMCAlgos.AK5PFJetsMCFlavourInfos_cfi import ak5JetFlavourInfos genJetFlavourPlusLeptonInfos = ak5JetFlavourInfos.clone( leptons = cms.InputTag("selectedHadronsAndPartons","leptons") ) from PhysicsTools.JetMCAlgos.GenHFHadronMatcher_cfi import matchGenHFHadron # Configuration for matching B-hadrons ================================================================ matchGenBHadron = matchGenHFHadron.clone( jetFlavourInfos = cms.InputTag("genJetFlavourPlusLeptonInfos"), flavour = 5 ) # Configuration for matching C-hadrons ================================================================= matchGenCHadron = matchGenHFHadron.clone( jetFlavourInfos = cms.InputTag("genJetFlavourPlusLeptonInfos"), flavour = 4 )
) process.flavourByValGenJet = cms.EDProducer("JetFlavourIdentifier", srcByReference = cms.InputTag("flavourByRefGenJet"), physicsDefinition = cms.bool(False) ) from PhysicsTools.JetMCAlgos.HadronAndPartonSelector_cfi import selectedHadronsAndPartons process.selectedHadronsAndPartons = selectedHadronsAndPartons.clone() from PhysicsTools.JetMCAlgos.AK4PFJetsMCFlavourInfos_cfi import ak4JetFlavourInfos process.myak4PFJetFlavourInfos = ak4JetFlavourInfos.clone(jets = cms.InputTag("ak4PFJets")) #ak5 jets from PhysicsTools.JetMCAlgos.AK5PFJetsMCFlavourInfos_cfi import ak5JetFlavourInfos process.myak5PFJetFlavourInfos = ak5JetFlavourInfos.clone(jets = cms.InputTag("ak5PFJets")) from PhysicsTools.JetMCAlgos.AK4PFJetsMCFlavourInfos_cfi import ak4JetFlavourInfos process.myak4PFJetCHSFlavourInfos = ak4JetFlavourInfos.clone(jets = cms.InputTag("ak4PFJetsCHS")) #from PhysicsTools.JetMCAlgos.AK4PFJetsMCFlavourInfos_cfi import ak4JetFlavourInfos #process.myak4CaloJetFlavourInfos = ak4JetFlavourInfos.clone(jets = cms.InputTag("ak4CaloJets")) from PhysicsTools.JetMCAlgos.AK4PFJetsMCFlavourInfos_cfi import ak4JetFlavourInfos process.myak4GenJetFlavourInfos = ak4JetFlavourInfos.clone(jets = cms.InputTag("ak4GenJets"))
# Herwig6: # /TT_CT10_AUET2_8TeV-powheg-herwig/Summer12_DR53X-PU_S10_START53_V19-v1/AODSIM #'/store/mc/Summer12_DR53X/TT_CT10_AUET2_8TeV-powheg-herwig/AODSIM/PU_S10_START53_V19-v1/10000/0039166F-C9CA-E211-9B10-0026189438E0.root' )) process.TFileService = cms.Service("TFileService", fileName=cms.string('histos.root')) #------------------------------------------------------------------------- # AK5 jets #------------------------------------------------------------------------- from PhysicsTools.JetMCAlgos.HadronAndPartonSelector_cfi import selectedHadronsAndPartons process.selectedHadronsAndPartons = selectedHadronsAndPartons.clone( particles=cms.InputTag("genParticles")) from PhysicsTools.JetMCAlgos.AK5PFJetsMCFlavourInfos_cfi import ak5JetFlavourInfos process.jetFlavourInfosAK5PFJets = ak5JetFlavourInfos.clone( jets=cms.InputTag("ak5PFJets")) #------------------------------------------------------------------------- process.jetFlavorAnalyzer = cms.EDAnalyzer( 'JetFlavorAnalyzer', src=cms.InputTag("genParticles"), jetFlavourInfos=cms.InputTag("jetFlavourInfosAK5PFJets"), jetPtMin=cms.double(30.)) process.p = cms.Path(process.selectedHadronsAndPartons * process.jetFlavourInfosAK5PFJets * process.jetFlavorAnalyzer)
import FWCore.ParameterSet.Config as cms # Flavour info: jet collection with all associated ghosts from PhysicsTools.JetMCAlgos.AK5PFJetsMCFlavourInfos_cfi import ak5JetFlavourInfos genJetFlavourPlusLeptonInfos = ak5JetFlavourInfos.clone( leptons=cms.InputTag("selectedHadronsAndPartons", "leptons")) from PhysicsTools.JetMCAlgos.GenHFHadronMatcher_cfi import matchGenHFHadron # Configuration for matching B-hadrons ================================================================ matchGenBHadron = matchGenHFHadron.clone( jetFlavourInfos=cms.InputTag("genJetFlavourPlusLeptonInfos"), flavour=5) # Configuration for matching C-hadrons ================================================================= matchGenCHadron = matchGenHFHadron.clone( jetFlavourInfos=cms.InputTag("genJetFlavourPlusLeptonInfos"), flavour=4)