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 customizeSignal(self, process): process.flashggGenPhotonsExtra.defaultType = 1 from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff import runMETs runMETs(process, True) #isMC # Default should be the right name for all signals process.load("flashgg/MicroAOD/flashggPDFWeightObject_cfi") process.p *= process.flashggPDFWeightObject
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 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'))
def customizeBackground(self, process): from flashgg.MicroAOD.flashggMet_RunCorrectionAndUncertainties_cff.py import runMETs runMETs(process, True) #isMC if "sherpa" in self.datasetName: process.flashggGenPhotonsExtra.defaultType = 1