Changeset 9040259 in git for modules/JetFlavorAssociation.cc

Timestamp:

Jun 17, 2015, 4:36:20 PM (9 years ago)

Author:

Pavel Demin <pavel.demin@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

a446115

Parents:

edf10ba

Message:

add BTaggingCMS and rename JetFlavourAssociation to JetFlavorAssociation

File:

• modules/JetFlavorAssociation.cc (moved) (moved from modules/JetFlavourAssociation.cc ) (14 diffs)

modules/JetFlavorAssociation.cc

@@ -18 +18 @@
 
 
-/** \class JetFlavourAssociation
+/** \class JetFlavorAssociation
  *
  *  Find origin of jet && evaluate jet flavor
@@ -26 +26 @@
  */
 
-#include "modules/JetFlavourAssociation.h"
+#include "modules/JetFlavorAssociation.h"
 
 #include "classes/DelphesClasses.h"
@@ -117 +117 @@
 //------------------------------------------------------------------------------
 
-JetFlavourAssociation::JetFlavourAssociation() :
+JetFlavorAssociation::JetFlavorAssociation() :
   fClassifier(0), fFilter(0),
   fItPartonInputArray(0), fItPartonInputArrayLHEF(0),
   fItJetInputArray(0), fItParticleInputArray(0)
 {
-  fClassifier    = new PartonClassifier;
+  fClassifier = new PartonClassifier;
   fClassifierLHEF = new PartonClassifierLHEF;
 }
@@ -128 +128 @@
 //------------------------------------------------------------------------------
 
-JetFlavourAssociation::~JetFlavourAssociation()
+JetFlavorAssociation::~JetFlavorAssociation()
 {
   if(fClassifier) delete fClassifier;
@@ -136 +136 @@
 //------------------------------------------------------------------------------
 
-void JetFlavourAssociation::Init()
+void JetFlavorAssociation::Init()
 {
   ExRootConfParam param;
@@ -168 +168 @@
 //------------------------------------------------------------------------------
 
-void JetFlavourAssociation::Finish()
+void JetFlavorAssociation::Finish()
 {
   if(fFilter) delete fFilter;
@@ -181 +181 @@
 //------------------------------------------------------------------------------
 
-void JetFlavourAssociation::Process(){
+void JetFlavorAssociation::Process(){
 
   Candidate *jet;
@@ -205 +205 @@
   {
     // get standard flavor
-    GetAlgoFlavour(jet, itPartonArray, itPartonLHEFArray);
-    GetPhysicsFlavour(jet, itPartonArray, itPartonLHEFArray);
+    GetAlgoFlavor(jet, itPartonArray, itPartonLHEFArray);
+    GetPhysicsFlavor(jet, itPartonArray, itPartonLHEFArray);
   }
 }
@@ -214 +214 @@
 // https://cmssdt.cern.ch/SDT/lxr/source/PhysicsTools/JetMCAlgos/plugins/JetPartonMatcher.cc?v=CMSSW_7_3_0_pre1
 
-void JetFlavourAssociation::GetAlgoFlavour(Candidate *jet, TIter &itPartonArray, TIter &itPartonLHEFArray)
+void JetFlavorAssociation::GetAlgoFlavor(Candidate *jet, TIter &itPartonArray, TIter &itPartonLHEFArray)
 {
   float maxPt = 0;
@@ -254 +254 @@
       {
         // partons are only quarks || gluons
-        int daughterFlavour1 = -1;
-        int daughterFlavour2 = -1;
-        if(parton->D1 != -1) daughterFlavour1 = TMath::Abs(static_cast<Candidate *>(fParticleInputArray->At(parton->D1))->PID);
-        if(parton->D2 != -1) daughterFlavour2 = TMath::Abs(static_cast<Candidate *>(fParticleInputArray->At(parton->D2))->PID);
-        if((daughterFlavour1 == 1 || daughterFlavour1 == 2 || daughterFlavour1 == 3 || daughterFlavour1 == 4 || daughterFlavour1 == 5 || daughterFlavour1 == 21)) daughterCounter++;
-        if((daughterFlavour2 == 1 || daughterFlavour2 == 2 || daughterFlavour2 == 3 || daughterFlavour2 == 4 || daughterFlavour1 == 5 || daughterFlavour2 == 21)) daughterCounter++;
+        int daughterFlavor1 = -1;
+        int daughterFlavor2 = -1;
+        if(parton->D1 != -1) daughterFlavor1 = TMath::Abs(static_cast<Candidate *>(fParticleInputArray->At(parton->D1))->PID);
+        if(parton->D2 != -1) daughterFlavor2 = TMath::Abs(static_cast<Candidate *>(fParticleInputArray->At(parton->D2))->PID);
+        if((daughterFlavor1 == 1 || daughterFlavor1 == 2 || daughterFlavor1 == 3 || daughterFlavor1 == 4 || daughterFlavor1 == 5 || daughterFlavor1 == 21)) daughterCounter++;
+        if((daughterFlavor2 == 1 || daughterFlavor2 == 2 || daughterFlavor2 == 3 || daughterFlavor2 == 4 || daughterFlavor1 == 5 || daughterFlavor2 == 21)) daughterCounter++;
       }
       if(daughterCounter > 0) continue;
@@ -296 +296 @@
 //------------------------------------------------------------------------------
 
-void JetFlavourAssociation::GetPhysicsFlavour(Candidate *jet, TIter &itPartonArray, TIter &itPartonLHEFArray)
+void JetFlavorAssociation::GetPhysicsFlavor(Candidate *jet, TIter &itPartonArray, TIter &itPartonLHEFArray)
 {
   float minDr = 1000;
@@ -303 +303 @@
   float dist;
   bool isGoodCandidate;
-  int contaminatingFlavour = 0;
+  int contaminatingFlavor = 0;
   int motherCounter = 0;
   Candidate *parton, *partonLHEF, *mother1, *mother2;
@@ -372 +372 @@
     {
       parton = *itContaminations;
-      contaminatingFlavour = TMath::Abs(parton->PID);
+      contaminatingFlavor = TMath::Abs(parton->PID);
       motherCounter = 0;
       if(parton->M1 != -1) motherCounter++;
@@ -391 +391 @@
       {
         // keep association --> the initialParton is a c --> the contaminated parton is a c
-        if(contaminatingFlavour == 4) continue;
+        if(contaminatingFlavor == 4) continue;
         jet->FlavorPhysics = 0; // all the other cases reject!
         break;