Fork me on GitHub

Changeset b286067 in git


Ignore:
Timestamp:
Aug 19, 2013, 3:40:54 PM (11 years ago)
Author:
pavel <pavel@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
76d3973
Parents:
82e07cf
Message:

added UsePTSum and PTSumMax parameters to the Isolation module

Files:
3 edited

Legend:

Unmodified
Added
Removed
  • CHANGELOG

    r82e07cf rb286067  
    113.0.10:
     2  added UsePTSum and PTSumMax parameters to the Isolation module
    23  added posibility to read reweighting information from LHEF
    34  updated ATLAS and CMS example cards
    45  added DelphesProMC.py reader
    5   added calculation of time after propagation in ParticlePropagator
     6  added calculation of time after propagation in the ParticlePropagator module
    67  improved energy flow algorithm
    7   fixed tau-jet matching in TauTagging module
     8  fixed tau-jet matching in the TauTagging module
    89  modified all the readers to read particle mass from input file
    910  adapted LHCO output files to counter a bug in MadAnalysis
  • modules/Isolation.cc

    r82e07cf rb286067  
    9090  fPTRatioMax = GetDouble("PTRatioMax", 0.1);
    9191
     92  fPTSumMax = GetDouble("PTSumMax", 5.0);
     93
     94  fUsePTSum = GetBool("UsePTSum", false);
     95
    9296  fClassifier->fPTMin = GetDouble("PTMin", 0.5);
    9397
     
    132136  Candidate *candidate, *isolation;
    133137  TObjArray *isolationArray;
    134   Double_t sumPT, ratio;
     138  Double_t sum, ratio;
    135139  Int_t counter;
    136140  Double_t rho = 0.0;
     
    157161
    158162    // loop over all input tracks
    159     sumPT = 0.0;
     163    sum = 0.0;
    160164    counter = 0;
    161165    itIsolationArray.Reset();
     
    167171         !candidate->Overlaps(isolation))
    168172      {
    169         sumPT += isolationMomentum.Pt();
     173        sum += isolationMomentum.Pt();
    170174        ++counter;
    171175      }
    172176    }
    173177
    174     // correct sumPT for pile-up contamination
    175     sumPT = sumPT - rho*fDeltaRMax*fDeltaRMax*TMath::Pi(); 
     178    // correct sum for pile-up contamination
     179    sum = sum - rho*fDeltaRMax*fDeltaRMax*TMath::Pi(); 
    176180
    177     ratio = sumPT/candidateMomentum.Pt();
    178     if(ratio > fPTRatioMax) continue;
     181    ratio = sum/candidateMomentum.Pt();
     182    if((fUsePTSum && sum > fPTSumMax) || ratio > fPTRatioMax) continue;
    179183
    180184    fOutputArray->Add(candidate);
  • modules/Isolation.h

    r82e07cf rb286067  
    4141  Double_t fPTRatioMax;
    4242
     43  Double_t fPTSumMax;
     44
     45  Bool_t fUsePTSum;
     46
    4347  IsolationClassifier *fClassifier; //!
    4448
Note: See TracChangeset for help on using the changeset viewer.