Fork me on GitHub

Changeset c614dd7 in git


Ignore:
Timestamp:
09/23/19 17:29:37 (12 months ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, master
Children:
3051358a, 84a097e
Parents:
ededa33
Message:

added neutral and charged energy fraction to jets

Files:
4 edited

Legend:

Unmodified
Added
Removed
  • classes/DelphesClasses.cc

    rededa33 rc614dd7  
    151151  NCharged(0),
    152152  NNeutrals(0),
     153  NeutralEnergyFraction(0),  // charged energy fraction
     154  ChargedEnergyFraction(0),  // neutral energy fraction
    153155  Beta(0),
    154156  BetaStar(0),
     
    322324  object.NCharged = NCharged;
    323325  object.NNeutrals = NNeutrals;
     326  object.NeutralEnergyFraction = NeutralEnergyFraction;
     327  object.ChargedEnergyFraction = ChargedEnergyFraction;
    324328  object.Beta = Beta;
    325329  object.BetaStar = BetaStar;
  • classes/DelphesClasses.h

    rededa33 rc614dd7  
    387387  Int_t NCharged; // number of charged constituents
    388388  Int_t NNeutrals; // number of neutral constituents
     389
     390  Float_t NeutralEnergyFraction;  // charged energy fraction
     391  Float_t ChargedEnergyFraction;  // neutral energy fraction
     392
    389393  Float_t Beta; // (sum pt of charged pile-up constituents)/(sum pt of charged constituents)
    390394  Float_t BetaStar; // (sum pt of charged constituents coming from hard interaction)/(sum pt of charged constituents)
     
    677681  Float_t PTD;
    678682  Float_t FracPt[5];
     683  Float_t NeutralEnergyFraction;  // charged energy fraction
     684  Float_t ChargedEnergyFraction;  // neutral energy fraction
     685
    679686
    680687  // Timing information
  • modules/FastJetFinder.cc

    rededa33 rc614dd7  
    316316  Double_t deta, dphi, detaMax, dphiMax;
    317317  Double_t time, timeWeight;
     318  Double_t neutralEnergyFraction, chargedEnergyFraction;
     319
    318320  Int_t number, ncharged, nneutrals;
    319321  Int_t charge;
     
    418420    nneutrals = 0;
    419421
     422    neutralEnergyFraction =0.;
     423    chargedEnergyFraction =0.;
     424
    420425    inputList.clear();
    421426    inputList = sequence->constituents(*itOutputList);
     
    432437
    433438      if(constituent->Charge == 0)
     439      {
    434440        nneutrals++;
     441        neutralEnergyFraction += constituent->Momentum.E();
     442      }
    435443      else
     444      {
    436445        ncharged++;
    437 
     446        chargedEnergyFraction += constituent->Momentum.E();
     447      }
     448     
    438449      time += TMath::Sqrt(constituent->Momentum.E()) * (constituent->Position.T());
    439450      timeWeight += TMath::Sqrt(constituent->Momentum.E());
     
    454465    candidate->NNeutrals = nneutrals;
    455466    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;
    456470
    457471    //for exclusive clustering, access y_n,n+1 as exclusive_ymerge (fNJets);
  • modules/TreeWriter.cc

    rededa33 rc614dd7  
    775775    entry->NCharged = candidate->NCharged;
    776776    entry->NNeutrals = candidate->NNeutrals;
     777
     778    entry->NeutralEnergyFraction = candidate->NeutralEnergyFraction;
     779    entry->ChargedEnergyFraction = candidate->ChargedEnergyFraction;
    777780    entry->Beta = candidate->Beta;
    778781    entry->BetaStar = candidate->BetaStar;
Note: See TracChangeset for help on using the changeset viewer.