Changes in modules/FastJetFinder.cc [c614dd7:341014c] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/FastJetFinder.cc
rc614dd7 r341014c 316 316 Double_t deta, dphi, detaMax, dphiMax; 317 317 Double_t time, timeWeight; 318 Double_t neutralEnergyFraction, chargedEnergyFraction;319 320 318 Int_t number, ncharged, nneutrals; 321 319 Int_t charge; … … 420 418 nneutrals = 0; 421 419 422 neutralEnergyFraction =0.;423 chargedEnergyFraction =0.;424 425 420 inputList.clear(); 426 421 inputList = sequence->constituents(*itOutputList); … … 437 432 438 433 if(constituent->Charge == 0) 439 {440 434 nneutrals++; 441 neutralEnergyFraction += constituent->Momentum.E();442 }443 435 else 444 {445 436 ncharged++; 446 chargedEnergyFraction += constituent->Momentum.E(); 447 } 448 437 449 438 time += TMath::Sqrt(constituent->Momentum.E()) * (constituent->Position.T()); 450 439 timeWeight += TMath::Sqrt(constituent->Momentum.E()); … … 465 454 candidate->NNeutrals = nneutrals; 466 455 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;470 456 471 457 //for exclusive clustering, access y_n,n+1 as exclusive_ymerge (fNJets);
Note:
See TracChangeset
for help on using the changeset viewer.