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)