Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/FastJetFinder.cc

    r01f457a r7c0fcd5  
    241241  Candidate *candidate, *constituent;
    242242  TLorentzVector momentum;
    243  
    244   TLorentzVector constmomentum;
    245  
    246243  Double_t deta, dphi, detaMax, dphiMax;
    247244  Double_t time, weightTime, avTime;
     
    307304  for(itOutputList = outputList.begin(); itOutputList != outputList.end(); ++itOutputList)
    308305  {
    309     jet = *itOutputList;
    310     if(fJetAlgorithm == 7) jet = join(jet.constituents());
    311    
    312     momentum.SetPxPyPzE(jet.px(), jet.py(), jet.pz(), jet.E());
    313    
     306    momentum.SetPxPyPzE(itOutputList->px(), itOutputList->py(), itOutputList->pz(), itOutputList->E());
    314307    area.reset(0.0, 0.0, 0.0, 0.0);
    315308    if(fAreaDefinition) area = itOutputList->area_4vector();
     
    322315    inputList.clear();
    323316    inputList = sequence->constituents(*itOutputList);
    324 
    325     constmomentum.SetPxPyPzE(0.0,0.0,0.0,0.0);;
    326317
    327318    for(itInputList = inputList.begin(); itInputList != inputList.end(); ++itInputList)
    328319    {
    329320      constituent = static_cast<Candidate*>(fInputArray->At(itInputList->user_index()));
    330 
    331       constmomentum += constituent->Momentum;
    332321
    333322      deta = TMath::Abs(momentum.Eta() - constituent->Momentum.Eta());
Note: See TracChangeset for help on using the changeset viewer.