"leadChIso","subleadChIso", "leadSigmaIeIe","subleadSigmaIeIe", "leadHoE","subleadHoE", ]) # set up event selection(s) from diphotons.Analysis.DiPhotonAnalysis import DiPhotonAnalysis analysis = DiPhotonAnalysis(diphotonDumper, massCut=massCut,ptLead=customize.ptLead,ptSublead=customize.ptSublead,scaling=customize.scaling, ## kinematic cuts computeMVA=False, genIsoDefinition=("genIso",10.), dataTriggers=dataTriggers, mcTriggers=mcTriggers, askTriggerOnMc=askTriggerOnMc, ## if mcTriggers is not empty will still compute efficiencies singlePhoDumperTemplate=photonDumper, applyDiphotonCorrections=customize.applyDiphotonCorrections, diphotonCorrectionsVersion=customize.diphotonCorrectionsVersion, sourceDiphotons=sourceDiphotons,sourceSinglePhotons=sourceSinglePhotons, extraSysModules=extraSysModules ) dumpKinTree=False if customize.datasetName(): # drop samples overlap if "GJet-HT" in customize.datasetName() or "GJets_DR-0p4_HT" in customize.datasetName(): analysis.keepPFOnly = True elif "QCD" in customize.datasetName(): analysis.keepFFOnly = True elif "DiPhotonJetsBox_MGG-80toInf_13TeV-Sherpa" in customize.datasetName():
("EE","1",0) ], variables=["genMass := mass","pt", "genLeadPt := leadingPhoton.pt", "genSubeadPt := subLeadingPhoton.pt", ], histograms=["genMass>>genmass(1500,0,15000)", ] ) from diphotons.Analysis.DiPhotonAnalysis import DiPhotonAnalysis analysis = DiPhotonAnalysis(None, massCut=customize.massCut,ptLead=customize.ptLead,ptSublead=customize.ptSublead,scaling=customize.scaling, ## kinematic cuts genIsoDefinition=("genIso",10.), dataTriggers=[], mcTriggers=[], removeEEEE=customize.removeEEEE ) analysis.addGenOnlySelection(process,genDiphotonDumper) if not customize.lastAttempt: # make sure process doesn't get stuck due to slow I/O process.watchDog = cms.EDAnalyzer("IdleWatchdog", minIdleFraction=cms.untracked.double(0.1), tolerance=cms.untracked.int32(10), checkEvery=cms.untracked.int32(100), ) process.watch = cms.Path( process.watchDog
"/store/group/phys_higgs/cmshgg/musella/flashgg/ExoPhys14ANv1/diphotonsPhys14AnV1/GGJets_M-1000To2000_Pt-50_13TeV-sherpa/ExoPhys14ANv1-diphotonsPhys14AnV1-v0-Phys14DR-PU20bx25_PHYS14_25_V1-v1/150330_192709/0000/diphotonsMicroAOD_1.root" )) process.TFileService = cms.Service("TFileService", fileName=cms.string("test.root")) # # analysis configuration # from diphotons.Analysis.DiPhotonAnalysis import DiPhotonAnalysis analysis = DiPhotonAnalysis( diphotonDumper, massCut=200., ptLead=100., ptSublead=100., scaling=False, ## kinematic cuts computeMVA=False, genIsoDefinition=("userFloat('genIso')", 10.), dataTriggers=["HLT_DoublePhoton85*", "HLT_Photon250_NoHE*"], mcTriggers=[], ## ["HLT_DoublePhoton85*","HLT_Photon250_NoHE*"], askTriggerOnMc= False, ## if mcTriggers is not empty will still compute efficiencies singlePhoDumperTemplate=photonDumper) ## kinematic selection analysis.addKinematicSelection(process, dumpTrees=True, splitByIso=True) ## analysis selections # CiC from diphotons.Analysis.highMassCiCDiPhotons_cfi import highMassCiCDiPhotons analysis.addAnalysisSelection( process,
variables=[ "genMass := mass", "pt", "genLeadPt := leadingPhoton.pt", "genSubeadPt := subLeadingPhoton.pt", ], histograms=[ "genMass>>genmass(1500,0,15000)", ]) from diphotons.Analysis.DiPhotonAnalysis import DiPhotonAnalysis analysis = DiPhotonAnalysis( None, massCut=customize.massCut, ptLead=customize.ptLead, ptSublead=customize.ptSublead, scaling=customize.scaling, ## kinematic cuts genIsoDefinition=("genIso", 10.), dataTriggers=[], mcTriggers=[], removeEEEE=customize.removeEEEE) analysis.addGenOnlySelection(process, genDiphotonDumper) if not customize.lastAttempt: # make sure process doesn't get stuck due to slow I/O process.watchDog = cms.EDAnalyzer( "IdleWatchdog", minIdleFraction=cms.untracked.double(0.1), tolerance=cms.untracked.int32(10), checkEvery=cms.untracked.int32(100), )
"leadMatchType","leadGenIso", "subleadMatchType","subleadGenIso", "leadPhoIsoEA","subleadPhoIsoEA", "leadPhoIso","subleadPhoIso", "leadChIso","subleadChIso", "leadSigmaIeIe","subleadSigmaIeIe", "leadHoE","subleadHoE", ]) from diphotons.Analysis.DiPhotonAnalysis import DiPhotonAnalysis analysis = DiPhotonAnalysis(diphotonDumper, massCut=massCut,ptLead=customize.ptLead,ptSublead=customize.ptSublead,scaling=customize.scaling, ## kinematic cuts computeMVA=True, genIsoDefinition=("genIso",10.), dataTriggers=dataTriggers, mcTriggers=mcTriggers, askTriggerOnMc=askTriggerOnMc, ## if mcTriggers is not empty will still compute efficiencies singlePhoDumperTemplate=photonDumper, ) dumpKinTree=False if customize.datasetName(): # drop samples overlap if "GJet-HT" in customize.datasetName(): analysis.keepPFOnly = True elif "QCD" in customize.datasetName(): analysis.keepFFOnly = True elif "DiPhotonJetsBox_MGG-80toInf_13TeV-Sherpa" in customize.datasetName(): analysis.vetoGenDiphotons = 200.
fileNames=cms.untracked.vstring("/store/group/phys_higgs/cmshgg/musella/flashgg/ExoPhys14ANv1/diphotonsPhys14AnV1/GGJets_M-1000To2000_Pt-50_13TeV-sherpa/ExoPhys14ANv1-diphotonsPhys14AnV1-v0-Phys14DR-PU20bx25_PHYS14_25_V1-v1/150330_192709/0000/diphotonsMicroAOD_1.root") ) process.TFileService = cms.Service("TFileService", fileName = cms.string("test.root") ) # # analysis configuration # from diphotons.Analysis.DiPhotonAnalysis import DiPhotonAnalysis analysis = DiPhotonAnalysis(diphotonDumper, massCut=200.,ptLead=100.,ptSublead=100.,scaling=False, ## kinematic cuts computeMVA=False, genIsoDefinition=("userFloat('genIso')",10.), dataTriggers=["HLT_DoublePhoton85*","HLT_Photon250_NoHE*"], mcTriggers=[],## ["HLT_DoublePhoton85*","HLT_Photon250_NoHE*"], askTriggerOnMc=False, ## if mcTriggers is not empty will still compute efficiencies singlePhoDumperTemplate=photonDumper ) ## kinematic selection analysis.addKinematicSelection(process,dumpTrees=True,splitByIso=True ) ## analysis selections # CiC from diphotons.Analysis.highMassCiCDiPhotons_cfi import highMassCiCDiPhotons analysis.addAnalysisSelection(process,"cic",highMassCiCDiPhotons,dumpTrees=True,dumpWorkspace=False,dumpHistos=True,splitByIso=True, dumperTemplate=minimalDumper,