Fork me on GitHub

Changeset df04eb1 in git for modules


Ignore:
Timestamp:
Dec 19, 2014, 9:34:12 PM (10 years ago)
Author:
Pavel Demin <pavel.demin@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
c59be54
Parents:
df35033
Message:

remove unused variables from FastJetFinder

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/FastJetFinder.cc

    rdf35033 rdf04eb1  
    237237  Candidate *candidate, *constituent;
    238238  TLorentzVector momentum;
    239  
    240   TLorentzVector constmomentum;
    241  
     239
    242240  Double_t deta, dphi, detaMax, dphiMax;
    243   Double_t time, weightTime, avTime;
     241  Double_t time, timeWeight;
    244242  Int_t number;
    245   Double_t rho = 0;
     243  Double_t rho = 0.0;
    246244  PseudoJet jet, area;
    247245  vector<PseudoJet> inputList, outputList;
     
    305303    jet = *itOutputList;
    306304    if(fJetAlgorithm == 7) jet = join(jet.constituents());
    307    
     305
    308306    momentum.SetPxPyPzE(jet.px(), jet.py(), jet.pz(), jet.E());
    309    
     307
    310308    area.reset(0.0, 0.0, 0.0, 0.0);
    311309    if(fAreaDefinition) area = itOutputList->area_4vector();
     
    313311    candidate = factory->NewCandidate();
    314312
    315     time=0;
    316     weightTime=0;
     313    time = 0.0;
     314    timeWeight = 0.0;
    317315
    318316    inputList.clear();
    319317    inputList = sequence->constituents(*itOutputList);
    320 
    321     constmomentum.SetPxPyPzE(0.0,0.0,0.0,0.0);;
    322318
    323319    for(itInputList = inputList.begin(); itInputList != inputList.end(); ++itInputList)
    324320    {
    325321      constituent = static_cast<Candidate*>(fInputArray->At(itInputList->user_index()));
    326 
    327       constmomentum += constituent->Momentum;
    328322
    329323      deta = TMath::Abs(momentum.Eta() - constituent->Momentum.Eta());
     
    333327
    334328      time += TMath::Sqrt(constituent->Momentum.E())*(constituent->Position.T());
    335       weightTime += TMath::Sqrt(constituent->Momentum.E());
     329      timeWeight += TMath::Sqrt(constituent->Momentum.E());
    336330
    337331      candidate->AddCandidate(constituent);
    338332    }
    339333
    340     avTime = time/weightTime;
    341 
    342334    candidate->Momentum = momentum;
    343     candidate->Position.SetT(avTime);
     335    candidate->Position.SetT(time/timeWeight);
    344336    candidate->Area.SetPxPyPzE(area.px(), area.py(), area.pz(), area.E());
    345337
Note: See TracChangeset for help on using the changeset viewer.