"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,
