Changeset 1180bc1 in git

Timestamp:

Jun 30, 2015, 2:38:29 PM (9 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

c627b07

Parents:

6153fb0

Message:

fixed BTagging without LHEF

Location:

modules

Files:

• BTagging.cc (modified) (7 diffs)
• BTagging.h (modified) (2 diffs)
• JetFlavorAssociation.cc (modified) (8 diffs)
• JetFlavorAssociation.h (modified) (1 diff)

modules/BTagging.cc

@@ -34 +34 @@
 #include "classes/DelphesFormula.h"
 
-#include "ExRootAnalysis/ExRootResult.h"
-#include "ExRootAnalysis/ExRootFilter.h"
-#include "ExRootAnalysis/ExRootClassifier.h"
-
 #include "TMath.h"
 #include "TString.h"
@@ -55 +51 @@
 //------------------------------------------------------------------------------
 
-class BTaggingPartonClassifier : public ExRootClassifier
+BTagging::BTagging() :
+  fItJetInputArray(0)
 {
-public:
-
-  BTaggingPartonClassifier() {}
-
-  Int_t GetCategory(TObject *object);
-
-  Double_t fEtaMax, fPTMin;
-};
-
-//------------------------------------------------------------------------------
-
-Int_t BTaggingPartonClassifier::GetCategory(TObject *object)
-{
-  Candidate *parton = static_cast<Candidate*>(object);
-  const TLorentzVector &momentum = parton->Momentum;
-  Int_t pdgCode;
-
-  if(momentum.Pt() <= fPTMin || TMath::Abs(momentum.Eta()) > fEtaMax) return -1;
-
-  pdgCode = TMath::Abs(parton->PID);
-  if(pdgCode != 21 && pdgCode > 5) return -1;
-
-  return 0;
-}
-
-//------------------------------------------------------------------------------
-
-BTagging::BTagging() :
-  fClassifier(0), fFilter(0),
-  fItPartonInputArray(0), fItJetInputArray(0)
-{
-  fClassifier = new BTaggingPartonClassifier;
 }
 
@@ -95 +60 @@
 BTagging::~BTagging()
 {
-  if(fClassifier) delete fClassifier;
 }
 
@@ -108 +72 @@
 
   fBitNumber = GetInt("BitNumber", 0);
-
-  fDeltaR = GetDouble("DeltaR", 0.5);
-
-  fClassifier->fPTMin = GetDouble("PartonPTMin", 1.0);
-  fClassifier->fEtaMax = GetDouble("PartonEtaMax", 2.5);
 
   // read efficiency formulas
@@ -139 +98 @@
   // import input array(s)
 
-  fPartonInputArray = ImportArray(GetString("PartonInputArray", "Delphes/partons"));
-  fItPartonInputArray = fPartonInputArray->MakeIterator();
-
-  fFilter = new ExRootFilter(fPartonInputArray);
-
   fJetInputArray = ImportArray(GetString("JetInputArray", "FastJetFinder/jets"));
   fItJetInputArray = fJetInputArray->MakeIterator();
@@ -155 +109 @@
   DelphesFormula *formula;
 
-  if(fFilter) delete fFilter;
   if(fItJetInputArray) delete fItJetInputArray;
-  if(fItPartonInputArray) delete fItPartonInputArray;
 
   for(itEfficiencyMap = fEfficiencyMap.begin(); itEfficiencyMap != fEfficiencyMap.end(); ++itEfficiencyMap)
@@ -172 +124 @@
   Candidate *jet;
   Double_t pt, eta, phi, e;
-  TObjArray *partonArray;
   map< Int_t, DelphesFormula * >::iterator itEfficiencyMap;
   DelphesFormula *formula;
-
-  // select quark and gluons
-  fFilter->Reset();
-  partonArray = fFilter->GetSubArray(fClassifier, 0);
-
-  if(partonArray == 0) return;
-
-  TIter itPartonArray(partonArray);
 
   // loop over all input jets

modules/BTagging.h

@@ -37 +37 @@
 class DelphesFormula;
 
-class ExRootFilter;
-class BTaggingPartonClassifier;
-
 class BTagging: public DelphesModule
 {
@@ -55 +52 @@
   Int_t fBitNumber;
 
-  Double_t fDeltaR;
-
 #if !defined(__CINT__) && !defined(__CLING__)
   std::map< Int_t, DelphesFormula * > fEfficiencyMap; //!
 #endif
 
-  BTaggingPartonClassifier *fClassifier; //!
-  
-  ExRootFilter *fFilter;
-
-  TIterator *fItPartonInputArray; //!
-  
   TIterator *fItJetInputArray; //!
-
-  const TObjArray *fPartonInputArray; //!
   
   const TObjArray *fJetInputArray; //!

modules/JetFlavorAssociation.cc

@@ -151 +151 @@
   fPartonInputArray = ImportArray(GetString("PartonInputArray", "Delphes/partons"));
   fItPartonInputArray = fPartonInputArray->MakeIterator();
+  fPartonFilter = new ExRootFilter(fPartonInputArray);
 
   fParticleInputArray = ImportArray(GetString("ParticleInputArray", "Delphes/allParticles"));
   fItParticleInputArray = fParticleInputArray->MakeIterator();
 
-  fParticleLHEFInputArray = ImportArray(GetString("ParticleLHEFInputArray", "Delphes/allParticlesLHEF"));
-  fItParticleLHEFInputArray = fParticleLHEFInputArray->MakeIterator();
-
-  fPartonFilter = new ExRootFilter(fPartonInputArray);
-  fParticleLHEFFilter = new ExRootFilter(fParticleLHEFInputArray);
+  try
+  {
+    fParticleLHEFInputArray = ImportArray(GetString("ParticleLHEFInputArray", "Delphes/allParticlesLHEF"));
+  }
+  catch(runtime_error &e)
+  {
+    fParticleLHEFInputArray = 0;
+  }
+
+  if(fParticleLHEFInputArray)
+  {
+    fItParticleLHEFInputArray = fParticleLHEFInputArray->MakeIterator();
+    fParticleLHEFFilter = new ExRootFilter(fParticleLHEFInputArray);
+  }
 
   fJetInputArray = ImportArray(GetString("JetInputArray", "FastJetFinder/jets"));
@@ -183 +193 @@
 
   Candidate *jet;
-  TObjArray *partonArray;
-  TObjArray *partonLHEFArray;
-
-  // select quark && gluons
+  TObjArray *partonArray = 0;
+  TObjArray *partonLHEFArray = 0;
+
+  // select quark and gluons
   fPartonFilter->Reset();
   partonArray = fPartonFilter->GetSubArray(fPartonClassifier, 0); // get the filtered parton array
-
   if(partonArray == 0) return;
-  TIter itPartonArray(partonArray);
-
-  fParticleLHEFFilter->Reset();
-  partonLHEFArray = fParticleLHEFFilter->GetSubArray(fParticleLHEFClassifier, 0); // get the filtered parton array
-
-  if(partonLHEFArray == 0) return;
-  TIter itPartonLHEFArray(partonLHEFArray);
-
+
+  if(fParticleLHEFInputArray) 
+  {
+    fParticleLHEFFilter->Reset();
+    partonLHEFArray = fParticleLHEFFilter->GetSubArray(fParticleLHEFClassifier, 0); // get the filtered parton array
+  }
   // loop over all input jets
   fItJetInputArray->Reset();
@@ -204 +211 @@
   {
     // get standard flavor
-    GetAlgoFlavor(jet, itPartonArray, itPartonLHEFArray);
-    GetPhysicsFlavor(jet, itPartonArray, itPartonLHEFArray);
+    GetAlgoFlavor(jet, partonArray, partonLHEFArray);
+    if(fParticleLHEFInputArray) GetPhysicsFlavor(jet, partonArray, partonLHEFArray);
   }
 }
@@ -213 +220 @@
 // https://cmssdt.cern.ch/SDT/lxr/source/PhysicsTools/JetMCAlgos/plugins/JetPartonMatcher.cc?v=CMSSW_7_3_0_pre1
 
-void JetFlavorAssociation::GetAlgoFlavor(Candidate *jet, TIter &itPartonArray, TIter &itPartonLHEFArray)
+void JetFlavorAssociation::GetAlgoFlavor(Candidate *jet, TObjArray *partonArray, TObjArray *partonLHEFArray)
 {
   float maxPt = 0;
-  bool isGoodParton = true;
   int daughterCounter = 0;
   Candidate *parton, *partonLHEF;
   Candidate *tempParton = 0, *tempPartonHighestPt = 0;
   int pdgCode, pdgCodeMax = -1;
+  
+  TIter itPartonArray(partonArray);
+  TIter itPartonLHEFArray(partonLHEFArray);
 
   itPartonArray.Reset();
@@ -232 +241 @@
       if(pdgCodeMax < pdgCode) pdgCodeMax = pdgCode;
     }
-
-    isGoodParton = true;
+ 
+    if(!fParticleLHEFInputArray) continue;
+ 
     itPartonLHEFArray.Reset();
     while((partonLHEF = static_cast<Candidate *>(itPartonLHEFArray.Next())))
@@ -240 +250 @@
          parton->PID == partonLHEF->PID &&
          partonLHEF->Charge == parton->Charge)
-      {
-         isGoodParton = false;
+      {      
          break;
       }
-
-      if(!isGoodParton) continue;
 
       // check the daugheter
@@ -285 +292 @@
 //------------------------------------------------------------------------------
 
-void JetFlavorAssociation::GetPhysicsFlavor(Candidate *jet, TIter &itPartonArray, TIter &itPartonLHEFArray)
+void JetFlavorAssociation::GetPhysicsFlavor(Candidate *jet, TObjArray *partonArray, TObjArray *partonLHEFArray)
 {
   int partonCounter = 0;
@@ -298 +305 @@
   vector<Candidate *>::iterator itContaminations;
 
+  TIter itPartonArray(partonArray);
+  TIter itPartonLHEFArray(partonLHEFArray);
+
   contaminations.clear();
 

modules/JetFlavorAssociation.h

@@ -50 +50 @@
   void Finish();
 
-  void GetAlgoFlavor(Candidate *jet, TIter &itPartonArray, TIter &itPartonLHEFArray);
-  void GetPhysicsFlavor(Candidate *jet, TIter &itPartonArray, TIter &itPartonLHEFArray);
+  void GetAlgoFlavor(Candidate *jet, TObjArray *partonArray, TObjArray *partonLHEFArray);
+  void GetPhysicsFlavor(Candidate *jet, TObjArray *partonArray, TObjArray *partonLHEFArray);
 
 private: