cms.PSet(TagName = cms.InputTag('flashggTTHHadronicTag')) )
print "customize.processId:",customize.processId
# load appropriate scale and smearing bins here
# systematics customization scripts will take care of adjusting flashggDiPhotonSystematics
#process.load("flashgg.Systematics.escales.escale76X_16DecRereco_2015")
# Or use the official tool instead
useEGMTools(process)
# Only run systematics for signal events
if customize.processId.count("h_") or customize.processId.count("vbf_") or customize.processId.count("Acceptance"): # convention: ggh vbf wzh (wh zh) tth
print "Signal MC, so adding systematics and dZ"
variablesToUse = minimalVariables
if customize.doFiducial:
variablesToUse.extend(fc.getGenVariables(True))
variablesToUse.extend(fc.getRecoVariables(True))
variablesToUse.append("genLeadGenIso := ? diPhoton().leadingPhoton().hasMatchedGenPhoton() ? diPhoton().leadingPhoton().userFloat(\"genIso\") : -99")
variablesToUse.append("decorrSigmarv := diPhotonMVA().decorrSigmarv")
variablesToUse.append("leadmva := diPhotonMVA().leadmva")
variablesToUse.append("subleadmva := diPhotonMVA().subleadmva")
if customize.doSystematics:
for direction in ["Up","Down"]:
phosystlabels.append("MvaShift%s01sigma" % direction)
# phosystlabels.append("MvaLinearSyst%s01sigma" % direction)
phosystlabels.append("SigmaEOverEShift%s01sigma" % direction)
phosystlabels.append("MaterialCentral%s01sigma" % direction)
phosystlabels.append("MaterialForward%s01sigma" % direction)
phosystlabels.append("FNUFEB%s01sigma" % direction)
phosystlabels.append("FNUFEE%s01sigma" % direction)
# print'checking customize options'
# print'customize.processId.count("ggF_X250_WWgg_qqlnugg") = ',customize.processId.count("ggF_X250_WWgg_qqlnugg")
# for thing in customize.processId.count(0):
# print'thing = ',thing
is_signal = reduce(
lambda y, z: y or z,
map(lambda x: customize.processId.count(x), signal_processes))
#if customize.processId.count("h_") or customize.processId.count("vbf_") or customize.processId.count("Acceptance") or customize.processId.count("hh_"):
if is_signal:
print "Signal MC, so adding systematics and dZ"
if customize.doHTXS:
variablesToUse = minimalVariablesHTXS
else:
variablesToUse = minimalVariables
if customize.doFiducial:
variablesToUse.extend(fc.getGenVariables(True))
variablesToUse.extend(fc.getRecoVariables(True))
variablesToUse.append(
"genLeadGenIso := ? diPhoton().leadingPhoton().hasMatchedGenPhoton() ? diPhoton().leadingPhoton().userFloat(\"genIso\") : -99"
)
variablesToUse.append("decorrSigmarv := diPhotonMVA().decorrSigmarv")
variablesToUse.append("leadmva := diPhotonMVA().leadmva")
variablesToUse.append("subleadmva := diPhotonMVA().subleadmva")
if customize.doSystematics:
for direction in ["Up", "Down"]:
phosystlabels.append("MvaShift%s01sigma" % direction)
# phosystlabels.append("MvaLinearSyst%s01sigma" % direction)
phosystlabels.append("SigmaEOverEShift%s01sigma" % direction)
phosystlabels.append("MaterialCentralBarrel%s01sigma" % direction)
phosystlabels.append("MaterialOuterBarrel%s01sigma" % direction)
