Merge pull request #184 from susy2015/finishSFAndSystImplementation

Finish sf and syst implementation
susy2015 · Aug 12, 2019 · 5792182 · 5792182
2 parents 6ce3962 + 26cf185
commit 5792182
Show file tree

Hide file tree

Showing 3 changed files with 13 additions and 24 deletions.
diff --git a/TopTagger/python/TopTagger.py b/TopTagger/python/TopTagger.py
@@ -132,7 +132,7 @@ def run(self, saveCandidates = False, saveSFAndSyst = False, *args, **kwargs):
                 return TopTaggerResult(results, sfAndSyst)
         return TopTaggerResult(results)
 
-    def runFromNanoAOD(self, event, isFirstEvent = False):
+    def runFromNanoAOD(self, event, isFirstEvent = False, saveSFAndSyst = False):
         #This is a hack for the nanoAOD postprocessor to force it to read all necessary variables before passing them to C because each new branch accessed causes all branches to be reallocated 
         nHackLoop = 1
         if isFirstEvent:
@@ -188,7 +188,7 @@ def runFromNanoAOD(self, event, isFirstEvent = False):
 
 
         if isMC:
-            results = self.run(ak4Inputs = ak4Inputs, resolvedTopInputs=resTopInputs, ak8Inputs=ak8Inputs, genInputs=genInputs, saveSFAndSyst = True)
+            results = self.run(ak4Inputs = ak4Inputs, resolvedTopInputs=resTopInputs, ak8Inputs=ak8Inputs, genInputs=genInputs, saveSFAndSyst = saveSFAndSyst)
         else:
             results = self.run(ak4Inputs = ak4Inputs, resolvedTopInputs=resTopInputs, ak8Inputs=ak8Inputs)
 

diff --git a/TopTagger/python/TopTaggerProducer.py b/TopTagger/python/TopTaggerProducer.py
@@ -127,40 +127,29 @@ def runFromNanoAOD(self, event):
             nJets = event.nJet
             jetInputs = (getattr(event, self.AK4JetInputs[0]), getattr(event, self.AK4JetInputs[1]), getattr(event, self.AK4JetInputs[2]), getattr(event, self.AK4JetInputs[3]), nJets)
 
+            inputDict = {}
+
             if self.doLepCleaning:
                 nElec = event.nElectron
                 nMuon = event.nMuon
-                ak4Inputs = (nJets, jetInputs, event.Jet_btagCSVV2, supplementaryFloatVariables, supplementaryIntVariables, event.Jet_electronIdx1, event.Jet_muonIdx1, nElec, (event.Electron_pt, event.Electron_eta, event.Electron_phi, event.Electron_mass, nElec), event.Electron_vidNestedWPBitmap, event.Electron_miniPFRelIso_all, nMuon, (event.Muon_pt, event.Muon_eta, event.Muon_phi, event.Muon_mass, nMuon), None, event.Muon_miniPFRelIso_all)
+                inputDict["ak4Inputs"] = (nJets, jetInputs, event.Jet_btagCSVV2, supplementaryFloatVariables, supplementaryIntVariables, event.Jet_electronIdx1, event.Jet_muonIdx1, nElec, (event.Electron_pt, event.Electron_eta, event.Electron_phi, event.Electron_mass, nElec), event.Electron_vidNestedWPBitmap, event.Electron_miniPFRelIso_all, nMuon, (event.Muon_pt, event.Muon_eta, event.Muon_phi, event.Muon_mass, nMuon), None, event.Muon_miniPFRelIso_all)
             else:
-                ak4Inputs = (nJets, jetInputs, event.Jet_btagCSVV2, supplementaryFloatVariables, supplementaryIntVariables)
+                inputDict["ak4Inputs"] = (nJets, jetInputs, event.Jet_btagCSVV2, supplementaryFloatVariables, supplementaryIntVariables)
 
             if self.useAK8:
                 nFatJet = event.nFatJet
                 nSubJet = event.nSubJet
-                ak8Inputs = (nFatJet, (event.FatJet_pt, event.FatJet_eta, event.FatJet_phi, event.FatJet_mass, nFatJet), event.FatJet_msoftdrop, event.FatJet_deepTag_TvsQCD, event.FatJet_deepTag_WvsQCD, nSubJet, (event.SubJet_pt, event.SubJet_eta, event.SubJet_phi, event.SubJet_mass, nSubJet), event.FatJet_subJetIdx1, event.FatJet_subJetIdx2)
-            else:
-                ak8Inputs = None
+                inputDict["ak8Inputs"] = (nFatJet, (event.FatJet_pt, event.FatJet_eta, event.FatJet_phi, event.FatJet_mass, nFatJet), event.FatJet_msoftdrop, event.FatJet_deepTag_TvsQCD, event.FatJet_deepTag_WvsQCD, nSubJet, (event.SubJet_pt, event.SubJet_eta, event.SubJet_phi, event.SubJet_mass, nSubJet), event.FatJet_subJetIdx1, event.FatJet_subJetIdx2)
 
             if not self.recalculateFromRawInputs:
                 nResCand = event.nResolvedTopCandidate
-                resTopInputs = (nResCand, (event.ResolvedTopCandidate_pt, event.ResolvedTopCandidate_eta, event.ResolvedTopCandidate_phi, event.ResolvedTopCandidate_mass, nResCand), event.ResolvedTopCandidate_discriminator, event.ResolvedTopCandidate_j1Idx, event.ResolvedTopCandidate_j2Idx, event.ResolvedTopCandidate_j3Idx)
-            else:
-                resTopInputs = None
+                inputDict["resTopInputs"] = (nResCand, (event.ResolvedTopCandidate_pt, event.ResolvedTopCandidate_eta, event.ResolvedTopCandidate_phi, event.ResolvedTopCandidate_mass, nResCand), event.ResolvedTopCandidate_discriminator, event.ResolvedTopCandidate_j1Idx, event.ResolvedTopCandidate_j2Idx, event.ResolvedTopCandidate_j3Idx)
 
             if self.isMC:
                 nGenPart = event.nGenPart
-                genInputs = (nGenPart, (event.GenPart_pt, event.GenPart_eta, event.GenPart_phi, event.GenPart_mass, nGenPart), event.GenPart_pdgId, event.GenPart_statusFlags, event.GenPart_genPartIdxMother)
-            else:
-                genInputs = None
-
-        if ak8Inputs and resTopInputs:
-            return self.tt.run(ak4Inputs = ak4Inputs, resolvedTopInputs=resTopInputs, ak8Inputs=ak8Inputs, saveCandidates=self.saveCandidates, genInputs=genInputs, saveSFAndSyst=self.saveSFAndSyst)
-        elif ak8Inputs and not resTopInputs:
-            return self.tt.run(ak4Inputs = ak4Inputs, ak8Inputs=ak8Inputs, saveCandidates=self.saveCandidates, genInputs=genInputs, saveSFAndSyst=self.saveSFAndSyst)
-        elif not ak8Inputs and resTopInputs:
-            return self.tt.run(ak4Inputs = ak4Inputs, resolvedTopInputs=resTopInputs, saveCandidates=self.saveCandidates, genInputs=genInputs, saveSFAndSyst=self.saveSFAndSyst)
-        else:
-            return self.tt.run(ak4Inputs = ak4Inputs, saveCandidates=self.saveCandidates, genInputs=genInputs, saveSFAndSyst=self.saveSFAndSyst)
+                inputDict["genInputs"] = (nGenPart, (event.GenPart_pt, event.GenPart_eta, event.GenPart_phi, event.GenPart_mass, nGenPart), event.GenPart_pdgId, event.GenPart_statusFlags, event.GenPart_genPartIdxMother)
+
+        return self.tt.run(saveSFAndSyst=self.saveSFAndSyst, **inputDict)
 
     def analyze(self, event):
         """process event, return True (go to next module) or False (fail, go to next event)"""
@@ -199,7 +188,7 @@ def analyze(self, event):
                         reorgSyst[syst].append(topSysts[syst])
                     except KeyError:
                         reorgSyst[syst].append(0.0)
-                        
+
             self.out.fillBranch("ResolvedTop%s_sf"%self.suffixResolved, ttr.sfCol()[resolvedFilter])
             for syst in self.systToSave:
                 self.out.fillBranch("ResolvedTop%s_syst_%s"%(self.suffixResolved, syst), reorgSyst[syst])

diff --git a/TopTagger/scripts/getTaggerCfg.sh b/TopTagger/scripts/getTaggerCfg.sh
@@ -187,7 +187,7 @@ MVAFILES=
 
 if [ -f TopTagger.cfg ]
 then
-    MVAFILES=$(grep "modelFile" TopTagger.cfg | sed 's/[^"]*"\([^"]*\)"/\1/')
+    MVAFILES=$(grep -e "modelFile" -e "inputFile"  TopTagger.cfg | sed 's/[^"]*"\([^"]*\)"/\1/')
     MISSING=
     if [[ ! -z ${MVAFILES// } ]]
     then