Fork me on GitHub

Changeset 7aea88d in git


Ignore:
Timestamp:
Jul 28, 2015, 3:12:14 PM (9 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
d1942df
Parents:
fb4337d
git-author:
Michele Selvaggi <michele.selvaggi@…> (07/28/15 15:11:03)
git-committer:
Michele Selvaggi <michele.selvaggi@…> (07/28/15 15:12:14)
Message:

now accounts for rhos calculated in different eta regions in Isolation module

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/Isolation.cc

    rfb4337d r7aea88d  
    121121  fItCandidateInputArray = fCandidateInputArray->MakeIterator();
    122122
    123   rhoInputArrayName = GetString("RhoInputArray", "");
     123  rhoInputArrayName = GetString("RhoInputArray", "Rho/rho");
    124124  if(rhoInputArrayName[0] != '\0')
    125125  {
     
    158158  Double_t rho = 0.0;
    159159
    160   if(fRhoInputArray && fRhoInputArray->GetEntriesFast() > 0)
    161   {
    162     candidate = static_cast<Candidate*>(fRhoInputArray->At(0));
    163     rho = candidate->Momentum.Pt();
    164   }
    165 
    166160  // select isolation objects
    167161  fFilter->Reset();
     
    178172    const TLorentzVector &candidateMomentum = candidate->Momentum;
    179173    eta = TMath::Abs(candidateMomentum.Eta());
     174
     175    // find rho
     176    rho = 0.0;
     177    if(fRhoInputArray)
     178    {
     179      fItRhoInputArray->Reset();
     180      while((object = static_cast<Candidate*>(fItRhoInputArray->Next())))
     181      {
     182        if(eta >= object->Edges[0] && eta < object->Edges[1])
     183        {
     184          rho = object->Momentum.Pt();
     185        }
     186      }
     187    }
    180188
    181189    // loop over all input tracks
Note: See TracChangeset for help on using the changeset viewer.