# use custum TMVA weights
process.flashggVBFMVA.vbfMVAweightfile = cms.FileInPath("flashgg/Taggers/data/Flashgg_VBF_CHS_STD_BDTG.weights.xml")
process.flashggVBFMVA.MVAMethod = cms.string("BDTG")
# QCD Recovery
# process.flashggVBFMVA.merge3rdJet = cms.untracked.bool(False)
# process.flashggVBFMVA.thirdJetDRCut = cms.untracked.double(1.5)
# combined MVA boundary set
process.flashggVBFTag.Boundaries = cms.vdouble(-2,0,2)
process.systematicsTagSequences = cms.Sequence()
from flashgg.MetaData.JobConfig import customize
customize.parse()
print "customize.processId:",customize.processId
# Do systematics on all MC (n.b. different from ordinary workspaces) but not on data
if doSystematics:
process.load("flashgg.Systematics.flashggDiPhotonSystematics_cfi")
massSearchReplaceAnyInputTag(process.flashggTagSequence,cms.InputTag("flashggDiPhotons"),cms.InputTag("flashggDiPhotonSystematics"))
systlabels = [""]
jetsystlabels = []
if customize.processId == "Data":
systprodlist = [getattr(process,"flashggJetSystematics%i"%i) for i in range(len(UnpackedJetCollectionVInputTag))]
for systprod in systprodlist:
# systprod.SystMethods = cms.VPSet() # do no systematics
# For any MicroAOD up to 1_3_0 the JEC in Data MicroAOD were bugged and this line makes sure they are fixed
# It should be a noop in cases where they are already correct
## HLT filter
from HLTrigger.HLTfilters.hltHighLevel_cfi import hltHighLevel
process.hltHighLevel = hltHighLevel.clone(
HLTPaths=cms.vstring("HLT_DoubleEle33_CaloIdL_MW_v*",
"HLT_DoubleEle33_CaloIdL_GsfTrkIdVL_MW_v*",
"HLT_Mu50_v*", "HLT_TkMu50_v*",
"HLT_Mu33_Ele33_CaloIdL_GsfTrkIdVL_v*",
"HLT_Ele27_WPTight_Gsf_v*", "HLT_IsoMu24_v*",
"HLT_IsoMu27_v*"),
throw=False
# andOr = cms.bool(True) # True = or between triggers
)
# import flashgg customization to check if we have signal or background
from flashgg.MetaData.JobConfig import customize
customize.parse()
# # flashgg tag sequence (for dipho MVA) and jet collections
process.load("flashgg/Taggers/flashggTagSequence_cfi")
# process.load("dafne/Taggers/flashggDiLeptonDiJetTagSequence_cfi") #uso questa?
# process.flashggTagSequence.remove(process.flashggUpdatedIdMVADiPhotons) # Needs to be run before systematics
# massSearchReplaceAnyInputTag(process.flashggTagSequence,cms.InputTag("flashggUpdatedIdMVADiPhotons"),cms.InputTag("flashggDiPhotonSystematics"))
# #remove un-necessary tags ...
# process.flashggTagSequence.remove(process.flashggVBFTag)
# process.flashggTagSequence.remove(process.flashggTTHLeptonicTag)
# process.flashggTagSequence.remove(process.flashggTTHHadronicTag)
# process.flashggTagSequence.remove(process.flashggVHEtTag)
# process.flashggTagSequence.remove(process.flashggVHLooseTag)
# process.flashggTagSequence.remove(process.flashggVHTightTag)
# process.flashggTagSequence.remove(process.flashggTagSorter)
