)
# restricted set of variables to be dumped for ancillary event selections
minimalDumper = diphotonDumper.clone()
cfgTools.dumpOnly(minimalDumper,
["mass","pt","genMass","satRegressedMass","regressedMass",
"leadEnergy","subLeadEnergy",
"leadSatRegressedEnergy","subLeadSatRegressedEnergy",
"leadRegressedEnergy","subLeadRegressedEnergy",
"leadIsSat","subLeadIsSat","leadIsWeird","subLeadIsWeird",
"leadPt","leadEta","leadScEta","leadPhi",
"subleadPt","subleadEta","subleadScEta","subleadPhi",
"leadBlockPhoIso","subleadBlockPhoIso",
"leadBlockChIso","subleadBlockChIso",
"leadRndConePhoIso","leadRndConeChIso",
"subleadRndConePhoIso","subleadRndConeChIso",
"leadRndConeChIso0","leadRndConeChIso1","leadRndConeChIso2","leadRndConeChIso3",
"leadRndConeChIso4","leadRndConeChIso5","leadRndConeChIso6","leadRndConeChIso7","leadRndConeChIso8",
"subleadRndConeChIso0","subleadRndConeChIso1","subleadRndConeChIso2","subleadRndConeChIso3",
"subleadRndConeChIso4","subleadRndConeChIso5","subleadRndConeChIso6","subleadRndConeChIso7","subleadRndConeChIso8",
"leadMatchType","leadGenIso",
"subleadMatchType","subleadGenIso",
"leadPhoIsoEA","subleadPhoIsoEA",
"leadPhoIso","subleadPhoIso",
"leadChIso","subleadChIso",
"leadSigmaIeIe","subleadSigmaIeIe",
"leadHoE","subleadHoE",
])
# set up event selection(s)
genDiphotonDumper.globalVariables = globalVariables
cfgTools.addCategories(genDiphotonDumper,
[("EB","max(abs(leadingPhoton.eta),abs(subLeadingPhoton.eta))<1.4442",0),
("EE","1",0)
],
variables=["genMass := mass","pt",
"genLeadPt := leadingPhoton.pt",
"genSubeadPt := subLeadingPhoton.pt",
],
histograms=["genMass>>genmass(1500,0,15000)",
]
)
# restricted set of variables to be dumped for ancillary event selections
minimalDumper = diphotonDumper.clone()
cfgTools.dumpOnly(minimalDumper,
minimalVariables)
# 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,
cfgTools.dumpOnly(minimalDumper, [
"mass",
"pt",
"leadPt",
"leadEta",
"leadScEta",
"leadPhi",
"subleadPt",
"subleadEta",
"subleadScEta",
"subleadPhi",
"leadBlockPhoIso",
"subleadBlockPhoIso",
"leadBlockChIso",
"subleadBlockChIso",
"leadRndConePhoIso",
"leadRndConeChIso",
"subleadRndConePhoIso",
"subleadRndConeChIso",
"leadMatchType",
"leadGenIso",
"subleadMatchType",
"subleadGenIso",
"leadPhoIsoEA",
"subleadPhoIsoEA",
"leadPhoIso",
"subleadPhoIso",
"leadChIso",
"subleadChIso",
"leadSigmaIeIe",
"subleadSigmaIeIe",
"leadHoE",
"subleadHoE",
])
]
)
minimalDumper = diphotonDumper.clone()
cfgTools.dumpOnly(minimalDumper,
["mass","pt","genMass","satRegressedMass","regressedMass",
"leadEnergy","subLeadEnergy",
"leadSatRegressedEnergy","subLeadSatRegressedEnergy",
"leadRegressedEnergy","subLeadRegressedEnergy",
"leadIsSat","subLeadIsSat","leadIsWeird","subLeadIsWeird",
"leadPt","leadEta","leadScEta","leadPhi",
"subleadPt","subleadEta","subleadScEta","subleadPhi",
"leadBlockPhoIso","subleadBlockPhoIso",
"leadBlockChIso","subleadBlockChIso",
"leadRndConePhoIso","leadRndConeChIso",
"subleadRndConePhoIso","subleadRndConeChIso",
"leadRndConeChIso0","leadRndConeChIso1","leadRndConeChIso2","leadRndConeChIso3",
"leadRndConeChIso4","leadRndConeChIso5","leadRndConeChIso6","leadRndConeChIso7","leadRndConeChIso8",
"subleadRndConeChIso0","subleadRndConeChIso1","subleadRndConeChIso2","subleadRndConeChIso3",
"subleadRndConeChIso4","subleadRndConeChIso5","subleadRndConeChIso6","subleadRndConeChIso7","subleadRndConeChIso8",
"leadMatchType","leadGenIso",
"subleadMatchType","subleadGenIso",
"leadPhoIsoEA","subleadPhoIsoEA",
"leadPhoIso","subleadPhoIso",
"leadChIso","subleadChIso",
"leadSigmaIeIe","subleadSigmaIeIe",
"leadHoE","subleadHoE",
])
from diphotons.Analysis.DiPhotonAnalysis import DiPhotonAnalysis
analysis = DiPhotonAnalysis(diphotonDumper,
"subleadPt:leadPt>>ptSubVsLead(145,100,3000:145,100,3000)",
"minR9>>minR9(110,0,1.1)",
"maxEta>>maxEta(250,0,2.5)"
]
)
minimalDumper = diphotonDumper.clone()
cfgTools.dumpOnly(minimalDumper,
["mass","pt",
"leadPt","leadEta","leadScEta","leadPhi",
"subleadPt","subleadEta","subleadScEta","subleadPhi",
"leadBlockPhoIso","subleadBlockPhoIso",
"leadBlockChIso","subleadBlockChIso",
"leadRndConePhoIso","leadRndConeChIso",
"subleadRndConePhoIso","subleadRndConeChIso",
"leadMatchType","leadGenIso",
"subleadMatchType","subleadGenIso",
"leadPhoIsoEA","subleadPhoIsoEA",
"leadPhoIso","subleadPhoIso",
"leadChIso","subleadChIso",
"leadSigmaIeIe","subleadSigmaIeIe",
"leadHoE","subleadHoE",
])
# single photon dumpoer
photonDumper.processId = "test"
photonDumper.dumpTrees = False
photonDumper.dumpWorkspace = False
photonDumper.quietRooFit = True
photonDumper.maxCandPerEvent=2
genDiphotonDumper,
[("EB", "max(abs(leadingPhoton.eta),abs(subLeadingPhoton.eta))<1.4442", 0),
("EE", "1", 0)],
variables=[
"genMass := mass",
"pt",
"genLeadPt := leadingPhoton.pt",
"genSubeadPt := subLeadingPhoton.pt",
],
histograms=[
"genMass>>genmass(1500,0,15000)",
])
# restricted set of variables to be dumped for ancillary event selections
minimalDumper = diphotonDumper.clone()
cfgTools.dumpOnly(minimalDumper, minimalVariables)
# 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
