Fork me on GitHub

Changeset 7e227ae in git


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

Files:
7 edited

Legend:

Unmodified
Added
Removed

  • Makefile

    r437e999 r7e227ae  
    325325        modules/BTagging.h \
    326326        modules/TauTagging.h \
     327        modules/TrackCountingTauTagging.h \
    327328        modules/TreeWriter.h \
    328329        modules/Merger.h \
    … …  
    816817        classes/DelphesFactory.h \
    817818        classes/DelphesFormula.h
     819tmp/modules/TrackCountingTauTagging.$(ObjSuf): \
     820        modules/TrackCountingTauTagging.$(SrcSuf) \
     821        modules/TrackCountingTauTagging.h \
     822        classes/DelphesClasses.h \
     823        classes/DelphesFactory.h \
     824        classes/DelphesFormula.h \
     825        external/ExRootAnalysis/ExRootResult.h \
     826        external/ExRootAnalysis/ExRootFilter.h \
     827        external/ExRootAnalysis/ExRootClassifier.h
    818828tmp/modules/TrackPileUpSubtractor.$(ObjSuf): \
    819829        modules/TrackPileUpSubtractor.$(SrcSuf) \
    … …  
    933943        tmp/modules/TimeSmearing.$(ObjSuf) \
    934944        tmp/modules/TrackCountingBTagging.$(ObjSuf) \
     945        tmp/modules/TrackCountingTauTagging.$(ObjSuf) \
    935946        tmp/modules/TrackPileUpSubtractor.$(ObjSuf) \
    936947        tmp/modules/TreeWriter.$(ObjSuf) \
    … …  
    12921303        external/fastjet/contribs/Nsubjettiness/Nsubjettiness.hh \
    12931304        external/fastjet/contribs/Nsubjettiness/Njettiness.hh \
    1294         external/fastjet/contribs/Nsubjettiness/NjettinessPlugin.hh
     1305        external/fastjet/contribs/Nsubjettiness/NjettinessPlugin.hh \
     1306        external/fastjet/contribs/Nsubjettiness/ExtraRecombiners.hh
    12951307tmp/modules/RunPUPPI.$(ObjSuf): \
    12961308        modules/RunPUPPI.$(SrcSuf) \
    … …  
    16961708        @touch $@
    16971709
     1710modules/TrackCountingTauTagging.h: \
     1711        classes/DelphesModule.h
     1712        @touch $@
     1713
    16981714external/fastjet/RectangularGrid.hh: \
    16991715        external/fastjet/PseudoJet.hh \

  • cards/delphes_card_ATLAS.tcl

    r437e999 r7e227ae  
    553553  # -2 - two or more mistag rate
    554554 
     555  set BitNumber 0
     556 
    555557  # taken from ATL-PHYS-PUB-2015-045 (medium working point)
    556      
    557558  add EfficiencyFormula {1} {0.70}
    558559  add EfficiencyFormula {2} {0.60}

  • cards/delphes_card_ATLAS_PileUp.tcl

    r437e999 r7e227ae  
    674674  # -2 - two or more mistag rate
    675675 
     676  set BitNumber 0
     677 
    676678  # taken from ATL-PHYS-PUB-2015-045 (medium working point)
    677      
    678679  add EfficiencyFormula {1} {0.70}
    679680  add EfficiencyFormula {2} {0.60}

  • modules/TauTagging.cc

    r437e999 r7e227ae  
    126126  Int_t i, size;
    127127
     128  fBitNumber = GetInt("BitNumber", 0);
     129
    128130  fDeltaR = GetDouble("DeltaR", 0.5);
    129131
    … …  
    193195  Candidate *jet, *tau, *daughter;
    194196  TLorentzVector tauMomentum;
    195   Double_t pt, eta, phi;
     197  Double_t pt, eta, phi, e;
    196198  TObjArray *tauArray;
    197199  map< Int_t, DelphesFormula * >::iterator itEfficiencyMap;
    … …  
    213215    phi = jetMomentum.Phi();
    214216    pt = jetMomentum.Pt();
     217    e = jetMomentum.E();
    215218
    216219    // loop over all input taus
    … …  
    252255
    253256    // apply an efficency formula
    254     jet->TauTag = gRandom->Uniform() <= formula->Eval(pt, eta);
     257    jet->TauTag |= (gRandom->Uniform() <= formula->Eval(pt, eta, phi, e)) << fBitNumber;
     258   
    255259    // set tau charge
    256260    jet->Charge = charge;

  • modules/TauTagging.h

    r437e999 r7e227ae  
    5555
    5656private:
    57 
     57 
     58  Int_t fBitNumber;
     59 
    5860  Double_t fDeltaR;
    5961

  • 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;

  • modules/TrackCountingTauTagging.h

    r437e999 r7e227ae  
    3939private:
    4040
     41  Int_t fBitNumber;
     42
    4143  Double_t fDeltaR;
    4244  Double_t fDeltaRTrack;
Note: See TracChangeset for help on using the changeset viewer.