Changeset 7e227ae in git for modules

Timestamp:

Dec 16, 2015, 3:15:15 PM (9 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

59d6164

Parents:

437e999

Message:

added BitNumber to TauTagging

Location:

modules

Files:

• TauTagging.cc (modified) (4 diffs)
• TauTagging.h (modified) (1 diff)
• TrackCountingTauTagging.cc (modified) (4 diffs)
• TrackCountingTauTagging.h (modified) (1 diff)

modules/TauTagging.cc

@@ -126 +126 @@
   Int_t i, size;
 
+  fBitNumber = GetInt("BitNumber", 0);
+
   fDeltaR = GetDouble("DeltaR", 0.5);
 
@@ -193 +195 @@
   Candidate *jet, *tau, *daughter;
   TLorentzVector tauMomentum;
-  Double_t pt, eta, phi;
+  Double_t pt, eta, phi, e;
   TObjArray *tauArray;
   map< Int_t, DelphesFormula * >::iterator itEfficiencyMap;
@@ -213 +215 @@
     phi = jetMomentum.Phi();
     pt = jetMomentum.Pt();
+    e = jetMomentum.E();
 
     // loop over all input taus
@@ -252 +255 @@
 
     // apply an efficency formula
-    jet->TauTag = gRandom->Uniform() <= formula->Eval(pt, eta);
+    jet->TauTag |= (gRandom->Uniform() <= formula->Eval(pt, eta, phi, e)) << fBitNumber;
+    
     // set tau charge
     jet->Charge = charge;

modules/TauTagging.h

@@ -55 +55 @@
 
 private:
-
+  
+  Int_t fBitNumber;
+  
   Double_t fDeltaR;
 

modules/TrackCountingTauTagging.cc

@@ -130 +130 @@
   Int_t i, size;
 
+  fBitNumber = GetInt("BitNumber", 0);
+
   fDeltaR = GetDouble("DeltaR", 0.5);
   fDeltaRTrack = GetDouble("DeltaRTrack", 0.2);
@@ -203 +205 @@
   Candidate *jet, *tau, *track, *daughter;
   TLorentzVector tauMomentum;
-  Double_t pt, eta, phi;
+  Double_t pt, eta, phi, e;
   TObjArray *tauArray;
   map< Int_t, DelphesFormula * >::iterator itEfficiencyMap;
@@ -228 +230 @@
     phi = jetMomentum.Phi();
     pt = jetMomentum.Pt();
+    e = jetMomentum.E();
+
 
 // loop over all input tracks
@@ -288 +292 @@
 
     // apply an efficency formula
-    jet->TauTag = gRandom->Uniform() <= formula->Eval(pt, eta);
+    jet->TauTag |= (gRandom->Uniform() <= formula->Eval(pt, eta, phi, e)) << fBitNumber;
+   
+   
     // set tau charge
     jet->Charge = charge;

modules/TrackCountingTauTagging.h

@@ -39 +39 @@
 private:
 
+  Int_t fBitNumber;
+
   Double_t fDeltaR;
   Double_t fDeltaRTrack;