Fork me on GitHub

Changeset 7429c6a in git


Ignore:
Timestamp:
Mar 15, 2018, 3:33:18 PM (6 years ago)
Author:
Martin <martin.flechl@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
f25215c
Parents:
0203656
Message:

Adding TauWeight variable

Files:
4 edited

Legend:

Unmodified
Added
Removed
  • classes/DelphesClasses.cc

    r0203656 r7429c6a  
    125125  Flavor(0), FlavorAlgo(0), FlavorPhys(0),
    126126  BTag(0), BTagAlgo(0), BTagPhys(0),
    127   TauTag(0), Eem(0.0), Ehad(0.0),
     127  TauTag(0), TauWeight(0.0), Eem(0.0), Ehad(0.0),
    128128  DeltaEta(0.0), DeltaPhi(0.0),
    129129  Momentum(0.0, 0.0, 0.0, 0.0),
     
    276276  object.BTagPhys = BTagPhys;
    277277  object.TauTag = TauTag;
     278  object.TauWeight = TauWeight;
    278279  object.Eem = Eem;
    279280  object.Ehad = Ehad;
     
    400401  BTagPhys = 0;
    401402  TauTag = 0;
     403  TauWeight = 0.0;
    402404  Eem = 0.0;
    403405  Ehad = 0.0;
  • classes/DelphesClasses.h

    r0203656 r7429c6a  
    373373
    374374  UInt_t TauTag; // 0 or 1 for a jet that has been tagged as a tau
     375  Float_t TauWeight; // probability for jet to be identified as tau
    375376
    376377  Int_t Charge; // tau charge
     
    556557
    557558  UInt_t TauTag;
     559  Float_t TauWeight;
    558560
    559561  Float_t Eem;
  • modules/TauTagging.cc

    r0203656 r7429c6a  
    255255
    256256    // apply an efficency formula
    257     jet->TauTag |= (gRandom->Uniform() <= formula->Eval(pt, eta, phi, e)) << fBitNumber;
    258    
     257    Double_t eff=formula->Eval(pt, eta, phi, e);
     258    jet->TauTag |= (gRandom->Uniform() <= eff) << fBitNumber;
     259    jet->TauWeight = eff;     
     260
    259261    // set tau charge
    260262    jet->Charge = charge;
  • modules/TreeWriter.cc

    r0203656 r7429c6a  
    650650
    651651    entry->TauTag = candidate->TauTag;
     652    entry->TauWeight = candidate->TauWeight;
    652653
    653654    entry->Charge = candidate->Charge;
Note: See TracChangeset for help on using the changeset viewer.