Fork me on GitHub

Changeset 1348 in svn


Ignore:
Timestamp:
Jan 3, 2014, 4:41:28 PM (11 years ago)
Author:
Michele Selvaggi
Message:

added PileUpJetID module from Seth Senz

Location:
trunk
Files:
2 added
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/Makefile

    r1345 r1348  
    327327        modules/JetPileUpSubtractor.h \
    328328        modules/TrackPileUpSubtractor.h \
     329        modules/PileUpJetID.h \
    329330        modules/ConstituentFilter.h \
    330331        modules/StatusPidFilter.h \
    … …  
    511512        modules/EnergyScale.$(SrcSuf) \
    512513        modules/EnergyScale.h \
     514        classes/DelphesClasses.h \
     515        classes/DelphesFactory.h \
     516        classes/DelphesFormula.h \
     517        external/ExRootAnalysis/ExRootResult.h \
     518        external/ExRootAnalysis/ExRootFilter.h \
     519        external/ExRootAnalysis/ExRootClassifier.h
     520tmp/modules/PileUpJetID.$(ObjSuf): \
     521        modules/PileUpJetID.$(SrcSuf) \
     522        modules/PileUpJetID.h \
    513523        classes/DelphesClasses.h \
    514524        classes/DelphesFactory.h \
    … …  
    948958        tmp/modules/ExampleModule.$(ObjSuf) \
    949959        tmp/modules/EnergyScale.$(ObjSuf) \
     960        tmp/modules/PileUpJetID.$(ObjSuf) \
    950961        tmp/modules/Cloner.$(ObjSuf) \
    951962        tmp/modules/TreeWriter.$(ObjSuf) \
    … …  
    12961307        @touch $@
    12971308
     1309modules/PileUpJetID.h: \
     1310        classes/DelphesModule.h
     1311        @touch $@
     1312
    12981313external/fastjet/version.hh: \
    12991314        external/fastjet/config.h

  • trunk/classes/DelphesClasses.cc

    r1318 r1348  
    112112  Area(0.0, 0.0, 0.0, 0.0),
    113113  fFactory(0),
    114   fArray(0)
     114  fArray(0),
     115  NCharged(0),
     116  NNeutrals(0),
     117  Beta(0),
     118  BetaStar(0),
     119  MeanSqDeltaR(0),
     120  PTD(0)
    115121{
    116122  Edges[0] = 0.0;
    … …  
    118124  Edges[2] = 0.0;
    119125  Edges[3] = 0.0;
     126  FracPt[0] = 0.0;
     127  FracPt[1] = 0.0;
     128  FracPt[2] = 0.0;
     129  FracPt[3] = 0.0;
     130  FracPt[4] = 0.0;
    120131}
    121132
    … …  
    206217  object.Area = Area;
    207218
     219  object.NCharged = NCharged;
     220  object.NNeutrals = NNeutrals;
     221  object.Beta = Beta;
     222  object.BetaStar = BetaStar;
     223  object.MeanSqDeltaR = MeanSqDeltaR;
     224  object.PTD = PTD;
     225  object.FracPt[0] = FracPt[0];
     226  object.FracPt[1] = FracPt[1];
     227  object.FracPt[2] = FracPt[2];
     228  object.FracPt[3] = FracPt[3];
     229  object.FracPt[4] = FracPt[4];
     230
    208231  object.fFactory = fFactory;
    209232  object.fArray = 0;
    … …  
    246269  Position.SetXYZT(0.0, 0.0, 0.0, 0.0);
    247270  Area.SetXYZT(0.0, 0.0, 0.0, 0.0);
     271
     272  NCharged = 0;
     273  NNeutrals = 0;
     274  Beta = 0.0;
     275  BetaStar = 0.0;
     276  MeanSqDeltaR = 0.0;
     277  PTD = 0.0;
     278  FracPt[0] = 0.0;
     279  FracPt[1] = 0.0;
     280  FracPt[2] = 0.0;
     281  FracPt[3] = 0.0;
     282  FracPt[4] = 0.0;
     283
    248284  fArray = 0;
    249285}

  • trunk/classes/DelphesClasses.h

    r1345 r1348  
    308308  TLorentzVector P4();
    309309
     310  // -- PileUpJetID variables ---
     311
     312  Int_t    NCharged;
     313  Int_t    NNeutrals;
     314  Float_t  Beta;
     315  Float_t  BetaStar;
     316  Float_t  MeanSqDeltaR;
     317  Float_t  PTD;
     318  Float_t  FracPt[5];
     319
     320
     321
    310322  ClassDef(Jet, 2)
    311323};
    … …  
    421433  virtual void Clear(Option_t* option = "");
    422434
     435// -- PileUpJetID variables ---
     436
     437  Int_t    NCharged;
     438  Int_t    NNeutrals;
     439  Float_t  Beta;
     440  Float_t  BetaStar;
     441  Float_t  MeanSqDeltaR;
     442  Float_t  PTD;
     443  Float_t  FracPt[5];
     444
    423445private:
    424446  DelphesFactory *fFactory; //!

  • trunk/modules/ModulesLinkDef.h

    r1345 r1348  
    3232#include "modules/JetPileUpSubtractor.h"
    3333#include "modules/TrackPileUpSubtractor.h"
     34#include "modules/PileUpJetID.h"
    3435#include "modules/ConstituentFilter.h"
    3536#include "modules/StatusPidFilter.h"
    … …  
    6465#pragma link C++ class JetPileUpSubtractor+;
    6566#pragma link C++ class TrackPileUpSubtractor+;
     67#pragma link C++ class PileUpJetID+;
    6668#pragma link C++ class ConstituentFilter+;
    6769#pragma link C++ class StatusPidFilter+;

  • trunk/modules/TrackPileUpSubtractor.cc

    r1345 r1348  
    102102  TIterator *iterator;
    103103  TObjArray *array;
    104   Double_t z, zvtx;
     104  Double_t z, zvtx=0;
    105105
    106106 

  • trunk/modules/TreeWriter.cc

    r1345 r1348  
    533533
    534534    entry->EhadOverEem = ecalEnergy > 0.0 ? hcalEnergy/ecalEnergy : 999.9;
     535 
     536    //---   Pile-Up Jet ID variables ----
     537
     538    entry->NCharged = candidate->NCharged;
     539    entry->NNeutrals = candidate->NNeutrals;
     540    entry->Beta = candidate->Beta;
     541    entry->BetaStar = candidate->BetaStar;
     542    entry->MeanSqDeltaR = candidate->MeanSqDeltaR;
     543    entry->PTD = candidate->PTD;
     544    entry->FracPt[0] = candidate->FracPt[0];
     545    entry->FracPt[1] = candidate->FracPt[1];
     546    entry->FracPt[2] = candidate->FracPt[2];
     547    entry->FracPt[3] = candidate->FracPt[3];
     548    entry->FracPt[4] = candidate->FracPt[4];
    535549
    536550    FillParticles(candidate, &entry->Particles);
Note: See TracChangeset for help on using the changeset viewer.