Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/FastJetFinder.cc

    rc614dd7 r341014c  
    316316  Double_t deta, dphi, detaMax, dphiMax;
    317317  Double_t time, timeWeight;
    318   Double_t neutralEnergyFraction, chargedEnergyFraction;
    319 
    320318  Int_t number, ncharged, nneutrals;
    321319  Int_t charge;
     
    420418    nneutrals = 0;
    421419
    422     neutralEnergyFraction =0.;
    423     chargedEnergyFraction =0.;
    424 
    425420    inputList.clear();
    426421    inputList = sequence->constituents(*itOutputList);
     
    437432
    438433      if(constituent->Charge == 0)
    439       {
    440434        nneutrals++;
    441         neutralEnergyFraction += constituent->Momentum.E();
    442       }
    443435      else
    444       {
    445436        ncharged++;
    446         chargedEnergyFraction += constituent->Momentum.E();
    447       }
    448      
     437
    449438      time += TMath::Sqrt(constituent->Momentum.E()) * (constituent->Position.T());
    450439      timeWeight += TMath::Sqrt(constituent->Momentum.E());
     
    465454    candidate->NNeutrals = nneutrals;
    466455    candidate->NCharged = ncharged;
    467 
    468     candidate->NeutralEnergyFraction = (momentum.E() > 0 ) ? neutralEnergyFraction/momentum.E() : 0.0;
    469     candidate->ChargedEnergyFraction = (momentum.E() > 0 ) ? chargedEnergyFraction/momentum.E() : 0.0;
    470456
    471457    //for exclusive clustering, access y_n,n+1 as exclusive_ymerge (fNJets);
Note: See TracChangeset for help on using the changeset viewer.