def customizeSignal(self,process):
process.flashggGenPhotonsExtra.defaultType = 1
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,True) #isMC
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
setMetCorr(process,multPhiCorr_MC_DY_80X)
process.p *=process.flashggMetSequence
process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
process.rivetProducerHTXS = cms.EDProducer('HTXSRivetProducer',
HepMCCollection = cms.InputTag('myGenerator','unsmeared'),
ProductionMode = cms.string('PRODUCTIONMODENOTSET'),
)
process.mergedGenParticles = cms.EDProducer("MergedGenParticleProducer",
inputPruned = cms.InputTag("prunedGenParticles"),
inputPacked = cms.InputTag("packedGenParticles"),
)
process.myGenerator = cms.EDProducer("GenParticles2HepMCConverterHTXS",
genParticles = cms.InputTag("mergedGenParticles"),
genEventInfo = cms.InputTag("generator"),
)
process.p *= process.mergedGenParticles
process.p *= process.myGenerator
process.p *= process.rivetProducerHTXS
process.out.outputCommands.append("keep *_rivetProducerHTXS_*_*")
self.customizePDFs(process)
def customizeData(self,process):
## remove MC-specific modules
modules = process.flashggMicroAODGenSequence.moduleNames()
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,False) #!isMC
if "2016G" in customize.datasetName or "2016H" in customize.datasetName:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_G_80X
setMetCorr(process,multPhiCorr_Data_G_80X)
elif "2016" in customize.datasetName:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_B_80X
setMetCorr(process,multPhiCorr_Data_B_80X)
else:
pass
process.p *=process.flashggMetSequence
for pathName in process.paths:
path = getattr(process,pathName)
for mod in modules:
path.remove( getattr(process,mod))
print getattr(process,pathName)
process.out.outputCommands.append("drop *_*Gen*_*_*")
process.out.outputCommands.append("keep *_reducedEgamma_*RecHit*_*") # for bad events
delattr(process,"flashggPrunedGenParticles") # will be run due to unscheduled mode unless deleted
self.customizeHighMassIsolations(process)
process.load("flashgg/MicroAOD/flashggDiPhotonFilter_cfi")
process.p1 = cms.Path(process.diPhotonFilter) # Do not save events with 0 diphotons
process.out.SelectEvents = cms.untracked.PSet(SelectEvents=cms.vstring('p1'))
def customizeSignal(self,process):
process.flashggGenPhotonsExtra.defaultType = 1
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,True) #isMC
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
setMetCorr(process,multPhiCorr_MC_DY_80X)
process.p *=process.flashggMetSequence
process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
process.rivetProducerHTXS = cms.EDProducer('HTXSRivetProducer',
HepMCCollection = cms.InputTag('myGenerator','unsmeared'),
ProductionMode = cms.string('PRODUCTIONMODENOTSET'),
)
process.mergedGenParticles = cms.EDProducer("MergedGenParticleProducer",
inputPruned = cms.InputTag("prunedGenParticles"),
inputPacked = cms.InputTag("packedGenParticles"),
)
process.myGenerator = cms.EDProducer("GenParticles2HepMCConverterHTXS",
genParticles = cms.InputTag("mergedGenParticles"),
genEventInfo = cms.InputTag("generator"),
)
process.p *= process.mergedGenParticles
process.p *= process.myGenerator
process.p *= process.rivetProducerHTXS
process.out.outputCommands.append("keep *_rivetProducerHTXS_*_*")
self.customizePDFs(process)
def customizeBackground(self,process):
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,True) #isMC
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
setMetCorr(process,multPhiCorr_MC_DY_80X)
if "sherpa" in self.datasetName:
process.flashggGenPhotonsExtra.defaultType = 1
def customizeBackground(self, process):
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs, setMetCorr
runMETs(process, True) #isMC
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
setMetCorr(process, multPhiCorr_MC_DY_80X)
if "sherpa" in self.datasetName:
process.flashggGenPhotonsExtra.defaultType = 1
def customizeBackground(self,process):
if not os.environ["CMSSW_VERSION"].count("CMSSW_9"):
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,True) #isMC
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
setMetCorr(process,multPhiCorr_MC_DY_80X)
process.p *=process.flashggMetSequence
if "sherpa" in self.datasetName:
process.flashggGenPhotonsExtra.defaultType = 1
def customizeBackground(self,process):
if not os.environ["CMSSW_VERSION"].count("CMSSW_9"):
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,True) #isMC
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
if not os.environ["CMSSW_VERSION"].count("CMSSW_9_2"):
setMetCorr(process,multPhiCorr_MC_DY_80X)
process.p *=process.flashggMetSequence
if "sherpa" in self.datasetName:
process.flashggGenPhotonsExtra.defaultType = 1
def customizeSignal(self,process):
process.flashggGenPhotonsExtra.defaultType = 1
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,True) #isMC
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
if not (os.environ["CMSSW_VERSION"].count("CMSSW_9_2") or os.environ["CMSSW_VERSION"].count("CMSSW_9_4")):
setMetCorr(process,multPhiCorr_MC_DY_80X)
process.p *=process.flashggMetSequence
if os.environ["CMSSW_VERSION"].count("CMSSW_8_0"):
process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
process.rivetProducerHTXS = cms.EDProducer('HTXSRivetProducer',
HepMCCollection = cms.InputTag('myGenerator','unsmeared'),
ProductionMode = cms.string('PRODUCTIONMODENOTSET'),
)
process.mergedGenParticles = cms.EDProducer("MergedGenParticleProducer",
inputPruned = cms.InputTag("prunedGenParticles"),
inputPacked = cms.InputTag("packedGenParticles"),
)
process.myGenerator = cms.EDProducer("GenParticles2HepMCConverterHTXS",
genParticles = cms.InputTag("mergedGenParticles"),
genEventInfo = cms.InputTag("generator"),
)
process.p *= process.mergedGenParticles
process.p *= process.myGenerator
process.p *= process.rivetProducerHTXS
process.out.outputCommands.append("keep *_rivetProducerHTXS_*_*")
if os.environ["CMSSW_VERSION"].count("CMSSW_9_4"):
#raise Exception,"Debugging ongoing for HTXS in CMSSW 9"
process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
process.rivetProducerHTXS = cms.EDProducer('HTXSRivetProducer',
HepMCCollection = cms.InputTag('myGenerator','unsmeared'),
LHERunInfo = cms.InputTag('externalLHEProducer'),
ProductionMode = cms.string('AUTO'),
)
process.mergedGenParticles = cms.EDProducer("MergedGenParticleProducer",
inputPruned = cms.InputTag("prunedGenParticles"),
inputPacked = cms.InputTag("packedGenParticles"),
)
process.myGenerator = cms.EDProducer("GenParticles2HepMCConverter",
genParticles = cms.InputTag("mergedGenParticles"),
genEventInfo = cms.InputTag("generator"),
signalParticlePdgIds = cms.vint32(25), ## for the Higgs analysis
)
process.p *= process.mergedGenParticles
process.p *= process.myGenerator
process.p *= process.rivetProducerHTXS
process.out.outputCommands.append('keep *_*_*_runRivetAnalysis')
self.customizePDFs(process)
def customizeSignal(self,process):
process.flashggGenPhotonsExtra.defaultType = 1
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,True) #isMC
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
if not (os.environ["CMSSW_VERSION"].count("CMSSW_9_2") or os.environ["CMSSW_VERSION"].count("CMSSW_9_4")):
setMetCorr(process,multPhiCorr_MC_DY_80X)
process.p *=process.flashggMetSequence
if os.environ["CMSSW_VERSION"].count("CMSSW_8_0"):
process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
process.rivetProducerHTXS = cms.EDProducer('HTXSRivetProducer',
HepMCCollection = cms.InputTag('myGenerator','unsmeared'),
ProductionMode = cms.string('PRODUCTIONMODENOTSET'),
)
process.mergedGenParticles = cms.EDProducer("MergedGenParticleProducer",
inputPruned = cms.InputTag("prunedGenParticles"),
inputPacked = cms.InputTag("packedGenParticles"),
)
process.myGenerator = cms.EDProducer("GenParticles2HepMCConverterHTXS",
genParticles = cms.InputTag("mergedGenParticles"),
genEventInfo = cms.InputTag("generator"),
)
process.p *= process.mergedGenParticles
process.p *= process.myGenerator
process.p *= process.rivetProducerHTXS
process.out.outputCommands.append("keep *_rivetProducerHTXS_*_*")
if os.environ["CMSSW_VERSION"].count("CMSSW_9_4"):
raise Exception,"Debugging ongoing for HTXS in CMSSW 9"
process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
process.rivetProducerHTXS = cms.EDProducer('HTXSRivetProducer',
HepMCCollection = cms.InputTag('myGenerator','unsmeared'),
LHERunInfo = cms.InputTag('externalLHEProducer'),
ProductionMode = cms.string('AUTO'),
)
process.mergedGenParticles = cms.EDProducer("MergedGenParticleProducer",
inputPruned = cms.InputTag("prunedGenParticles"),
inputPacked = cms.InputTag("packedGenParticles"),
)
process.myGenerator = cms.EDProducer("GenParticles2HepMCConverter",
genParticles = cms.InputTag("mergedGenParticles"),
genEventInfo = cms.InputTag("generator"),
signalParticlePdgIds = cms.vint32(25), ## for the Higgs analysis
)
process.p *= process.mergedGenParticles
process.p *= process.myGenerator
process.p *= process.rivetProducerHTXS
process.out.outputCommands.append('keep *_*_*_runRivetAnalysis')
self.customizePDFs(process)
def customizeBackground(self,process):
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
if os.environ["CMSSW_VERSION"].count("CMSSW_8_0"):
era="Summer16_23Sep2016V4_MC"
elif os.environ["CMSSW_VERSION"].count("CMSSW_9_4"):
era="Fall17_17Nov2017_V6_MC"
runMETs(process,era)
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
if not (os.environ["CMSSW_VERSION"].count("CMSSW_9_2") or os.environ["CMSSW_VERSION"].count("CMSSW_9_4")):
setMetCorr(process,multPhiCorr_MC_DY_80X)
process.p *=process.flashggMetSequence
if "sherpa" in self.datasetName:
process.flashggGenPhotonsExtra.defaultType = 1
def customizeData(self, process):
print "CUSTOMIZE DATA"
## remove MC-specific modules
modules = process.flashggMicroAODGenSequence.moduleNames()
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs, setMetCorr
if os.environ["CMSSW_VERSION"].count("CMSSW_8_0"):
era = "Summer16_23Sep2016AllV4_DATA"
elif os.environ["CMSSW_VERSION"].count("CMSSW_9_4"):
era = "Fall17_17Nov2017BCDEF_V6_DATA"
runMETs(process, era)
if "2016G" in customize.datasetName or "2016H" in customize.datasetName:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_G_80X
setMetCorr(process, multPhiCorr_Data_G_80X)
elif "2016" in customize.datasetName:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_B_80X
setMetCorr(process, multPhiCorr_Data_B_80X)
else:
pass
process.p *= process.flashggMetSequence
for pathName in process.paths:
path = getattr(process, pathName)
for mod in modules:
path.remove(getattr(process, mod))
print getattr(process, pathName)
process.out.outputCommands.append("drop *_*Gen*_*_*")
process.out.outputCommands.append(
"keep *_reducedEgamma_*RecHit*_*") # for bad events
delattr(process, "flashggPrunedGenParticles"
) # will be run due to unscheduled mode unless deleted
delattr(process, "flashggGenPhotons"
) # will be run due to unscheduled mode unless deleted
delattr(process, "flashggGenPhotonsExtra"
) # will be run due to unscheduled mode unless deleted
delattr(process, "flashggGenNeutrinos"
) # will be run due to unscheduled mode unless deleted
from flashgg.MicroAOD.flashggJets_cfi import maxJetCollections
for vtx in range(0, maxJetCollections):
# getattr(process,"flashggPFCHSJets%i"%vtx).Debug = True
delattr(process, "patJetGenJetMatchAK4PFCHSLeg%i" % vtx)
delattr(process, "patJetFlavourAssociationAK4PFCHSLeg%i" % vtx)
delattr(process, "patJetPartons%i" % vtx)
delattr(process, "patJetPartonMatchAK4PFCHSLeg%i" % vtx)
self.customizeHighMassIsolations(process)
process.load("flashgg/MicroAOD/flashggDiPhotonFilter_cfi")
process.flashggDiPhotonFilterSequence += process.diPhotonSelector
process.flashggDiPhotonFilterSequence += process.diPhotonFilter # Do not continue running events with 0 diphotons passing pt cuts
process.p1 = cms.Path(
process.diPhotonFilter
) # Do not save events with 0 diphotons passing pt cuts
process.out.SelectEvents = cms.untracked.PSet(
SelectEvents=cms.vstring('p1'))
def customizeData(self,process):
print "CUSTOMIZE DATA"
## remove MC-specific modules
modules = process.flashggMicroAODGenSequence.moduleNames()
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
if os.environ["CMSSW_VERSION"].count("CMSSW_8_0"):
era="Summer16_23Sep2016AllV4_DATA"
elif os.environ["CMSSW_VERSION"].count("CMSSW_9_4"):
era="Fall17_17Nov2017BCDEF_V6_DATA"
runMETs(process,era)
if "2016G" in customize.datasetName or "2016H" in customize.datasetName:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_G_80X
setMetCorr(process,multPhiCorr_Data_G_80X)
elif "2016" in customize.datasetName:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_B_80X
setMetCorr(process,multPhiCorr_Data_B_80X)
else:
pass
process.p *= process.flashggMetSequence
for pathName in process.paths:
path = getattr(process,pathName)
for mod in modules:
path.remove( getattr(process,mod))
print getattr(process,pathName)
process.out.outputCommands.append("drop *_*Gen*_*_*")
process.out.outputCommands.append("keep *_reducedEgamma_*RecHit*_*") # for bad events
delattr(process,"flashggPrunedGenParticles") # will be run due to unscheduled mode unless deleted
delattr(process,"flashggGenPhotons") # will be run due to unscheduled mode unless deleted
delattr(process,"flashggGenPhotonsExtra") # will be run due to unscheduled mode unless deleted
delattr(process,"flashggGenNeutrinos") # will be run due to unscheduled mode unless deleted
from flashgg.MicroAOD.flashggJets_cfi import maxJetCollections
for vtx in range(0,maxJetCollections):
# getattr(process,"flashggPFCHSJets%i"%vtx).Debug = True
delattr(process,"patJetGenJetMatchAK4PFCHSLeg%i"%vtx)
delattr(process,"patJetFlavourAssociationAK4PFCHSLeg%i"%vtx)
delattr(process,"patJetPartons%i"%vtx)
delattr(process,"patJetPartonMatchAK4PFCHSLeg%i"%vtx)
self.customizeHighMassIsolations(process)
process.load("flashgg/MicroAOD/flashggDiPhotonFilter_cfi")
process.flashggDiPhotonFilterSequence += process.diPhotonSelector
process.flashggDiPhotonFilterSequence += process.diPhotonFilter # Do not continue running events with 0 diphotons passing pt cuts
process.p1 = cms.Path(process.diPhotonFilter) # Do not save events with 0 diphotons passing pt cuts
process.out.SelectEvents = cms.untracked.PSet(SelectEvents=cms.vstring('p1'))
###---Add HLT filter as first step of MicroAOD sequence
if self.addMicroAODHLTFilter:
process.triggerFilterModule = getMicroAODHLTFilter(customize.datasetName)
if process.triggerFilterModule:
process.p = cms.Path(process.triggerFilterModule*process.p._seq)
process.p1 = cms.Path(process.triggerFilterModule*process.p1._seq)
def customizeBackground(self, process):
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs, setMetCorr
if os.environ["CMSSW_VERSION"].count("CMSSW_8_0"):
era = "Summer16_23Sep2016V4_MC"
elif os.environ["CMSSW_VERSION"].count("CMSSW_9_4"):
era = "Fall17_17Nov2017_V6_MC"
runMETs(process, era)
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_MC_DY_80X
if not (os.environ["CMSSW_VERSION"].count("CMSSW_9_2")
or os.environ["CMSSW_VERSION"].count("CMSSW_9_4")):
setMetCorr(process, multPhiCorr_MC_DY_80X)
process.p *= process.flashggMetSequence
if "sherpa" in self.datasetName:
process.flashggGenPhotonsExtra.defaultType = 1
def customizeData(self,process):
print "CUSTOMIZE DATA"
## remove MC-specific modules
modules = process.flashggMicroAODGenSequence.moduleNames()
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,False) #!isMC
if "2016G" in customize.datasetName or "2016H" in customize.datasetName:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_G_80X
setMetCorr(process,multPhiCorr_Data_G_80X)
elif "2016" in customize.datasetName:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_B_80X
setMetCorr(process,multPhiCorr_Data_B_80X)
else:
pass
process.p *=process.flashggMetSequence
for pathName in process.paths:
path = getattr(process,pathName)
for mod in modules:
path.remove( getattr(process,mod))
print getattr(process,pathName)
process.out.outputCommands.append("drop *_*Gen*_*_*")
process.out.outputCommands.append("keep *_reducedEgamma_*RecHit*_*") # for bad events
delattr(process,"flashggPrunedGenParticles") # will be run due to unscheduled mode unless deleted
delattr(process,"flashggGenPhotons") # will be run due to unscheduled mode unless deleted
delattr(process,"flashggGenPhotonsExtra") # will be run due to unscheduled mode unless deleted
delattr(process,"flashggGenNeutrinos") # will be run due to unscheduled mode unless deleted
from flashgg.MicroAOD.flashggJets_cfi import maxJetCollections
for vtx in range(0,maxJetCollections):
# getattr(process,"flashggPFCHSJets%i"%vtx).Debug = True
delattr(process,"patJetGenJetMatchAK4PFCHSLeg%i"%vtx)
delattr(process,"patJetFlavourAssociationAK4PFCHSLeg%i"%vtx)
delattr(process,"patJetPartons%i"%vtx)
delattr(process,"patJetPartonMatchAK4PFCHSLeg%i"%vtx)
self.customizeHighMassIsolations(process)
process.load("flashgg/MicroAOD/flashggDiPhotonFilter_cfi")
process.p1 = cms.Path(process.diPhotonFilter) # Do not save events with 0 diphotons
process.out.SelectEvents = cms.untracked.PSet(SelectEvents=cms.vstring('p1'))
def customizeData(self,process):
## remove MC-specific modules
modules = process.flashggMicroAODGenSequence.moduleNames()
from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs,setMetCorr
runMETs(process,False) #!isMC
if "2016G" in customize.datasetName:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_G_80X
setMetCorr(process,multPhiCorr_Data_G_80X)
else:
from flashgg.MicroAOD.METcorr_multPhiCorr_80X_sumPt_cfi import multPhiCorr_Data_B_80X
setMetCorr(process,multPhiCorr_Data_B_80X)
for pathName in process.paths:
path = getattr(process,pathName)
for mod in modules:
path.remove( getattr(process,mod))
print getattr(process,pathName)
process.out.outputCommands.append("drop *_*Gen*_*_*")
process.out.outputCommands.append("keep *_*_*RecHit*_*") # for bad events
delattr(process,"flashggPrunedGenParticles") # will be run due to unscheduled mode unless deleted
self.customizeHighMassIsolations(process)
process.load("flashgg/MicroAOD/flashggDiPhotonFilter_cfi")
process.p1 = cms.Path(process.diPhotonFilter) # Do not save events with 0 diphotons
process.out.SelectEvents = cms.untracked.PSet(SelectEvents=cms.vstring('p1'))