Fork me on GitHub

Ignore:
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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/TrackCountingTauTagging.cc

    r437e999 r7e227ae  
    130130  Int_t i, size;
    131131
     132  fBitNumber = GetInt("BitNumber", 0);
     133
    132134  fDeltaR = GetDouble("DeltaR", 0.5);
    133135  fDeltaRTrack = GetDouble("DeltaRTrack", 0.2);
     
    203205  Candidate *jet, *tau, *track, *daughter;
    204206  TLorentzVector tauMomentum;
    205   Double_t pt, eta, phi;
     207  Double_t pt, eta, phi, e;
    206208  TObjArray *tauArray;
    207209  map< Int_t, DelphesFormula * >::iterator itEfficiencyMap;
     
    228230    phi = jetMomentum.Phi();
    229231    pt = jetMomentum.Pt();
     232    e = jetMomentum.E();
     233
    230234
    231235// loop over all input tracks
     
    288292
    289293    // apply an efficency formula
    290     jet->TauTag = gRandom->Uniform() <= formula->Eval(pt, eta);
     294    jet->TauTag |= (gRandom->Uniform() <= formula->Eval(pt, eta, phi, e)) << fBitNumber;
     295   
     296   
    291297    // set tau charge
    292298    jet->Charge = charge;
Note: See TracChangeset for help on using the changeset viewer.