def defineAnalysis(process,castorDir="",reRunRho=False) :
# our MET producer
process.ClusteredPFMetProducer = ClusteredPFMetProducer.clone()
#
# configure the analyzer (cf. base values are in the StandardSelections_cfi)
#
process.evAnalyzer = cms.EDAnalyzer("DileptonPlusMETEventAnalyzer",
dtag=cms.string('llvv'),
Trigger = BaseTriggerSelection.clone(),
Generator = BaseGeneratorSelection.clone(),
Vertices = BaseVertexSelection.clone(),
Photons = BasePhotonsSelection.clone(),
LooseMuons = BaseLooseMuonsSelection.clone(),
Muons = BaseMuonsSelection.clone(),
LooseElectrons = BaseLooseElectronsSelection.clone(),
Electrons = BaseElectronsSelection.clone(),
Dileptons = BaseDileptonSelection.clone(),
Jets = BaseJetSelection.clone(),
AssocJets = AssocJetSelection.clone(),
MET = BaseMetSelection.clone()
)
if(castorDir.find('12_03_13_HZZ2l2v_pat')>0):
print ' Warning: applying PF muon patch to standard selection'
process.evAnalyzer.Muons.source=cms.InputTag("selectedPatMuonsPFlow")
process.evAnalyzer.Muons.sourceIsPF=cms.bool(True)
process.evAnalyzer.LooseMuons.source=process.evAnalyzer.Muons.source
process.evAnalyzer.LooseMuons.sourceIsPF=process.evAnalyzer.Muons.sourceIsPF
#MVAs for IDs
#process.load('EGamma.EGammaAnalysisTools.electronIdMVAProducer_cfi')
#process.mvaIDs = cms.Sequence( process.mvaTrigV0 + process.mvaNonTrigV0 )
from JetMETCorrections.Configuration.JetCorrectionProducers_cff import ak5PFJetsL1L2L3
process.ak5PFJetsL1L2L3ForMVAMET=ak5PFJetsL1L2L3.clone()
#rho for muon isolation
if(reRunRho) :
process.pfAllNeutralHadronsAndPhotons = pfAllNeutralHadronsAndPhotons.clone( src=cms.InputTag("particleFlow") )
process.kt6PFJetsCentralNeutral = kt4PFJets.clone( src = cms.InputTag("pfAllNeutralHadronsAndPhotons"),
rParam = cms.double(0.6),
doAreaFastjet = cms.bool(True),
doRhoFastjet = cms.bool(True),
Ghost_EtaMax = cms.double(3.1),
Rho_EtaMax = cms.double(2.5),
inputEtMin = cms.double(0.5)
)
process.rhoForIsolationSequence = cms.Sequence(process.pfAllNeutralHadronsAndPhotons*process.kt6PFJetsCentralNeutral)
process.analysis = cms.Path( #process.mvaIDs +
process.rhoForIsolationSequence +
process.ak5PFJetsL1L2L3ForMVAMET +
process.ClusteredPFMetProducer +
process.evAnalyzer)
else :
process.analysis = cms.Path( #process.mvaIDs +
process.ak5PFJetsL1L2L3ForMVAMET +
process.ClusteredPFMetProducer +
process.evAnalyzer)
def defineAnalysis(process,castorDir="",reRunRho=False) :
# our MET producer
process.ClusteredPFMetProducer = ClusteredPFMetProducer.clone()
#
# configure the analyzer (cf. base values are in the StandardSelections_cfi)
#
process.dataAnalyzer = cms.EDAnalyzer("DileptonPlusMETEventAnalyzer",
dtag=cms.string('llvv'),
Trigger = BaseTriggerSelection.clone(),
MetFilter = BaseMetFilters.clone(),
Generator = BaseGeneratorSelection.clone(),
Vertices = BaseVertexSelection.clone(), #checked
Photons = BasePhotonsSelection.clone(),
SoftMuons = BaseSoftMuonsSelection.clone(),
Muons = BaseMuonsSelection.clone(),
LooseElectrons = BaseLooseElectronsSelection.clone(),
Electrons = BaseElectronsSelection.clone(),
Dileptons = BaseDileptonSelection.clone(),
Jets = BaseJetSelection.clone(),
AssocJets = AssocJetSelection.clone(),
MET = BaseMetSelection.clone()
)
#if(castorDir.find('12_03_13_HZZ2l2v_pat')>0):
# print ' Warning: applying PF muon patch to standard selection'
# process.evAnalyzer.Muons.source=cms.InputTag("selectedPatMuonsPFlow")
# process.evAnalyzer.Muons.sourceIsPF=cms.bool(True)
# process.evAnalyzer.LooseMuons.source=process.evAnalyzer.Muons.source
# process.evAnalyzer.LooseMuons.sourceIsPF=process.evAnalyzer.Muons.sourceIsPF
#MVAs for IDs
#process.load('EGamma.EGammaAnalysisTools.electronIdMVAProducer_cfi')
#process.mvaIDs = cms.Sequence( process.mvaTrigV0 + process.mvaNonTrigV0 )
from JetMETCorrections.Configuration.JetCorrectionProducers_cff import ak5PFJetsL1L2L3
process.ak5PFJetsL1L2L3ForMVAMET=ak5PFJetsL1L2L3.clone()
