analysis.vetoGenDiphotons = 200.
# alyaws get full info for signal
if "Grav" in customize.datasetName():
dumpKinTree=True
minimalDumper=diphotonDumper
dumpTrees=True
if customize.histosOnly:
dumpKinTree=False
dumpTrees=False
dumpNm1Trees=dumpTrees and not invertEleVeto
## kinematic selection
analysis.addKinematicSelection(process,dumpTrees=dumpKinTree,splitByIso=True
)
if not dumpKinTree: minimalDumper=diphotonDumper
## analysis selections
if customize.idversion != "":
if customize.idversion == "V2":
from diphotons.Analysis.highMassCiCDiPhotons_cfi import highMassCiCDiPhotonsV2 as highMassCiCDiPhotons
from diphotons.Analysis.highMassCiCDiPhotons_cfi import highMassCiCDiPhotonsSBV2 as highMassCiCDiPhotonsSB
else:
print "Unknown ID version %s " % customize.idversion
sys.exit(-1)
else:
from diphotons.Analysis.highMassCiCDiPhotons_cfi import highMassCiCDiPhotons, highMassCiCDiPhotonsSB
if invertEleVeto:
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,
nMinusOne=[
(0, "NoChIso", True, False,
