def addDiMuonSeparationVariables(process, sequence, treeProducer):
from MuonAnalysis.TagAndProbe.nearbyMuonsInfo_cfi import nearbyMuonsInfo;
tpp = treeProducer.tagProbePairs.moduleLabel
if not hasattr(process, tpp+"DiMuonInfo"):
setattr(process, tpp+"DiMuonInfo", nearbyMuonsInfo.clone(src = treeProducer.tagProbePairs))
sequence.replace(getattr(process,tpp), getattr(process,tpp) + getattr(process, tpp+"DiMuonInfo"))
if not hasattr(treeProducer, 'pairVariables'):
treeProducer.pairVariables = cms.PSet()
treeProducer.pairFlags = cms.PSet()
treeProducer.pairVariables.dphiVtxTimesQ = cms.InputTag(tpp+"DiMuonInfo", "dphiVtxTimesQ")
treeProducer.pairVariables.drM2 = cms.InputTag(tpp+"DiMuonInfo", "drM2")
treeProducer.pairVariables.dphiM2 = cms.InputTag(tpp+"DiMuonInfo", "dphiM2")
treeProducer.pairVariables.distM2 = cms.InputTag(tpp+"DiMuonInfo", "distM2")
treeProducer.pairVariables.drStaIn = cms.InputTag(tpp+"DiMuonInfo", "drStaIn")
treeProducer.pairVariables.dphiStaIn = cms.InputTag(tpp+"DiMuonInfo", "dphiStaIn")
def addDiMuonSeparationVariables(process, sequence, treeProducer):
from MuonAnalysis.TagAndProbe.nearbyMuonsInfo_cfi import nearbyMuonsInfo
tpp = treeProducer.tagProbePairs.moduleLabel
if not hasattr(process, tpp + "DiMuonInfo"):
setattr(process, tpp + "DiMuonInfo",
nearbyMuonsInfo.clone(src=treeProducer.tagProbePairs))
sequence.replace(
getattr(process, tpp),
getattr(process, tpp) + getattr(process, tpp + "DiMuonInfo"))
if not hasattr(treeProducer, 'pairVariables'):
treeProducer.pairVariables = cms.PSet()
treeProducer.pairFlags = cms.PSet()
treeProducer.pairVariables.dphiVtxTimesQ = cms.InputTag(
tpp + "DiMuonInfo", "dphiVtxTimesQ")
treeProducer.pairVariables.drM2 = cms.InputTag(tpp + "DiMuonInfo", "drM2")
treeProducer.pairVariables.dphiM2 = cms.InputTag(tpp + "DiMuonInfo",
"dphiM2")
treeProducer.pairVariables.distM2 = cms.InputTag(tpp + "DiMuonInfo",
"distM2")
treeProducer.pairVariables.drStaIn = cms.InputTag(tpp + "DiMuonInfo",
"drStaIn")
treeProducer.pairVariables.dphiStaIn = cms.InputTag(
tpp + "DiMuonInfo", "dphiStaIn")
process.load(
"Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.GlobalTag.globaltag = 'GR_R_44_V12::All'
from MuonAnalysis.MuonAssociators.TnPTreeProducer_separation_variables_cfi import *
process.tnpTreeJPsi = tnpTreeJPsiProducer.clone()
###################
treeProducer = process.tnpTreeJPsi
from MuonAnalysis.TagAndProbe.nearbyMuonsInfo_cfi import nearbyMuonsInfo
nearbyMuonsInfo.filename = cms.string("Debug_JPsi_MuOnia_13_2.root")
tpp = treeProducer.tagProbePairs.moduleLabel
setattr(process, tpp + "DiMuonInfo",
nearbyMuonsInfo.clone(src=treeProducer.tagProbePairs))
if not hasattr(treeProducer, 'pairVariables'):
treeProducer.pairVariables = cms.PSet()
treeProducer.pairFlags = cms.PSet()
treeProducer.pairVariables.dphiVtxTimesQ = cms.InputTag(
tpp + "DiMuonInfo", "dphiVtxTimesQ")
treeProducer.pairVariables.drM2 = cms.InputTag(tpp + "DiMuonInfo", "drM2")
treeProducer.pairVariables.dphiM2 = cms.InputTag(tpp + "DiMuonInfo", "dphiM2")
treeProducer.pairVariables.distM2 = cms.InputTag(tpp + "DiMuonInfo", "distM2")
treeProducer.pairVariables.drStaIn = cms.InputTag(tpp + "DiMuonInfo",
"drStaIn")
treeProducer.pairVariables.dphiStaIn = cms.InputTag(tpp + "DiMuonInfo",
"dphiStaIn")
#####################
import FWCore.ParameterSet.Config as cms
from MuonAnalysis.TagAndProbe.nearbyMuonsInfo_cfi import nearbyMuonsInfo as tagProbeSeparation
tagProbeStaSeparation = tagProbeSeparation.clone(
src=cms.InputTag("tpPairsSta"), )
#########################################################################################
## Object counting modules ##
#########################################################################################
nverticesModule = cms.EDProducer(
"VertexMultiplicityCounter",
probes=cms.InputTag("tagMuons"),
objects=cms.InputTag("offlinePrimaryVertices"),
objectSelection=cms.string(
"!isFake && ndof > 4 && abs(z) <= 25 && position.Rho <= 2"),
)
njets30Module = cms.EDProducer(
"CandCleanedMultiplicityCounter",
pairs=cms.InputTag("tpPairs"),
objects=cms.InputTag("ak4PFJetsCHS"),
objectSelection=cms.string("abs(eta) < 5 && pt > 30"),
minTagObjDR=cms.double(0.3),
minProbeObjDR=cms.double(0.3),
)
#########################################################################################
## Tracking-related modules ##
#########################################################################################
process.load("Configuration.StandardSequences.Reconstruction_cff")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
#process.GlobalTag.globaltag = 'GR_R_42_V14::All'
#process.GlobalTag.globaltag = 'GR_R_42_V21::All'
process.GlobalTag.globaltag = 'GR_R_44_V12::All'
from MuonAnalysis.MuonAssociators.TnPTreeProducer_separation_variables_Z_cfi import *
process.tnpTreeZ = tnpTreeZProducer.clone()
###################
treeProducer = process.tnpTreeZ
from MuonAnalysis.TagAndProbe.nearbyMuonsInfo_cfi import nearbyMuonsInfo
nearbyMuonsInfo.filename = cms.string("Debug_Z_IsoMu24_2p1_1.root")
tpp = treeProducer.tagProbePairs.moduleLabel
setattr(process, tpp+"DiMuonInfo", nearbyMuonsInfo.clone(src = treeProducer.tagProbePairs))
if not hasattr(treeProducer, 'pairVariables'):
treeProducer.pairVariables = cms.PSet()
treeProducer.pairFlags = cms.PSet()
treeProducer.pairVariables.dphiVtxTimesQ = cms.InputTag(tpp+"DiMuonInfo", "dphiVtxTimesQ")
treeProducer.pairVariables.drM2 = cms.InputTag(tpp+"DiMuonInfo", "drM2")
treeProducer.pairVariables.dphiM2 = cms.InputTag(tpp+"DiMuonInfo", "dphiM2")
treeProducer.pairVariables.distM2 = cms.InputTag(tpp+"DiMuonInfo", "distM2")
treeProducer.pairVariables.drStaIn = cms.InputTag(tpp+"DiMuonInfo", "drStaIn")
treeProducer.pairVariables.dphiStaIn = cms.InputTag(tpp+"DiMuonInfo", "dphiStaIn")
#####################
process.TFileService = cms.Service("TFileService", fileName = cms.string('rfio:/castor/cern.ch/user/r/rybinska/Feb/Z_IsoMu24_2p1_1.root'))
process.tnpSequence = cms.Sequence(