ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Utils/interface/IsolationTools.h
(Generate patch)

Comparing UserCode/MitPhysics/Utils/interface/IsolationTools.h (file contents):
Revision 1.15 by ceballos, Fri May 13 08:27:15 2011 UTC vs.
Revision 1.27 by bendavid, Fri May 25 19:41:11 2012 UTC

# Line 13 | Line 13
13  
14   #include <TMath.h>
15   #include "MitAna/DataTree/interface/Track.h"
16 + #include "MitAna/DataTree/interface/Photon.h"
17   #include "MitAna/DataTree/interface/BasicCluster.h"
18   #include "MitAna/DataTree/interface/SuperCluster.h"
19   #include "MitAna/DataTree/interface/CaloTower.h"
# Line 22 | Line 23
23   #include "MitAna/DataTree/interface/PFCandidateCol.h"
24   #include "MitAna/DataTree/interface/TrackCol.h"
25   #include "MitAna/DataTree/interface/DecayParticleCol.h"
26 + #include "MitAna/DataTree/interface/PileupEnergyDensityCol.h"
27 + #include "MitPhysics/Utils/interface/ElectronTools.h"
28  
29   namespace mithep
30   {
# Line 39 | Line 42 | namespace mithep
42        static Double_t CaloTowerEmIsolation(const ThreeVector *p, Double_t extRadius,
43                                             Double_t intRadius, Double_t etLow,
44                                             const mithep::Collection<mithep::CaloTower> *caloTowers);
45 <      static Double_t PFMuonIsolation(const Muon *p, const PFCandidateCol *PFCands,
46 <                                      const Vertex *vertex, Double_t  delta_z, Double_t ptMin,
44 <                                      Double_t extRadius, Double_t intRadius);
45 >      static Double_t PFRadialMuonIsolation(const Muon *p, const PFCandidateCol *PFCands,
46 >                                            Double_t ptMin = 1.0, Double_t extRadius = 0.3);
47        static Double_t PFMuonIsolation(const Muon *p, const PFCandidateCol *PFCands, const Vertex *vertex,
48 <                                      const MuonCol *goodMuons, const ElectronCol *goodElectrons,
49 <                                      Double_t  delta_z = 0.1, Double_t ptMin = 1.0, Double_t extRadius = 0.4,
50 <                                      Double_t intRadius = 0.0, int isoType = 0, Double_t beta = 1.0);
48 >                                      Double_t  delta_z = 0.1, Double_t ptMin = 1.0,
49 >                                      Double_t extRadius = 0.4, Double_t intRadiusGamma = 0.07, Double_t intRadius = 0.0);
50 >      static Double_t PFMuonIsolation(const Muon *p, const PFCandidateCol *PFCands,
51 >                                      const MuonCol *goodMuons, const ElectronCol *goodElectrons,
52 >                                      const Vertex *vertex, Double_t  delta_z, Double_t ptMin,
53 >                                      Double_t extRadius, Double_t intRadiusGamma, Double_t intRadius);
54        static Double_t PFElectronIsolation(const Electron *p, const PFCandidateCol *PFCands,
55                                            const Vertex *vertex, Double_t  delta_z, Double_t ptMin,
56 <                                          Double_t extRadius, Double_t intRadius);
57 <      static Double_t PFElectronIsolation(const Electron *p, const PFCandidateCol *PFCands, const Vertex *vertex,
58 <                                          const MuonCol *goodMuons, const ElectronCol *goodElectrons,
59 <                                          Double_t  delta_z = 0.1 , Double_t ptMin = 1.0,Double_t extRadius = 0.4,                                        
60 <                                          Double_t intRadius = 0.0, int isoType = 0, Double_t beta = 1.0);
61 <      static Double_t BetaM(const TrackCol *tracks, const Muon *p, const Vertex *vertex,
56 >                                          Double_t extRadius, Double_t intRadius, Int_t PFCandidateType = -1);
57 >      static Double_t PFElectronIsolation(const Electron *p, const PFCandidateCol *PFCands,
58 >                                          const MuonCol *goodMuons, const ElectronCol *goodElectrons,
59 >                                          const Vertex *vertex, Double_t  delta_z, Double_t ptMin,
60 >                                          Double_t extRadius, Double_t intRadius, Int_t PFCandidateType = -1);
61 >      static Double_t PFElectronIsolation2012(const Electron *ele, const Vertex *vertex,
62 >                                              const PFCandidateCol *PFCands,
63 >                                              const PileupEnergyDensityCol *PileupEnergyDensity,
64 >                                              ElectronTools::EElectronEffectiveAreaTarget EffectiveAreaTarget,
65 >                                              const ElectronCol *goodElectrons,
66 >                                              const MuonCol *goodMuons, Double_t dRMax = 0.4, Bool_t isDebug = kFALSE);
67 >       static Double_t BetaM(const TrackCol *tracks, const Muon *p, const Vertex *vertex,
68                              Double_t ptMin, Double_t  delta_z, Double_t extRadius,
69                              Double_t intRadius);
70        static Double_t BetaE(const TrackCol *tracks, const Electron *p, const Vertex *vertex,
# Line 72 | Line 83 | namespace mithep
83                                           UInt_t maxNExpectedHitsInner = 999,
84                                           const mithep::DecayParticleCol *conversions = 0);
85  
86 <    ClassDef(IsolationTools, 0) // Isolation tools
87 <  };
86 >      // methods for Hgg BaseLien Selection. These isolation are stupid, but what can we do.... ;(
87 >      static Double_t CiCTrackIsolation(const mithep::Photon*,
88 >                                        const BaseVertex*,
89 >                                        Double_t extRadius,
90 >                                        Double_t intRadius,
91 >                                        Double_t ptLow,
92 >                                        Double_t etaStrip,
93 >                                        Double_t maxD0,
94 >                                        Double_t maxDZ,
95 >                                        const mithep::Collection<mithep::Track> *tracks,
96 >                                        unsigned int* worstVtxIdx = NULL,
97 >                                        const mithep::Collection<mithep::Vertex> *vtxs = NULL,
98 >                                        const mithep::Collection<mithep::Electron> *eles = NULL,
99 >                                        bool print=false,
100 >                                        double* ptmax=NULL,
101 >                                        double* dRmax=NULL);
102 >      
103 >      static Double_t PFChargedIsolation(const mithep::Photon *p,
104 >                                         const BaseVertex *theVtx,
105 >                                         Double_t extRadius,
106 >                                         Double_t intRadius,
107 >                                         const PFCandidateCol *PFCands,
108 >                                         unsigned int* worstVtxIndex = NULL,
109 >                                         const mithep::Collection<mithep::Vertex> *vtxs = NULL,
110 >                                         bool print = NULL);
111 >
112 >      static Double_t PFGammaIsolation(const mithep::Photon *p,
113 >                                         Double_t extRadius,
114 >                                         Double_t intRadius,
115 >                                         const PFCandidateCol *PFCands);                                        
116 >
117 >      static Double_t PFNeutralHadronIsolation(const mithep::Photon *p,
118 >                                         Double_t extRadius,
119 >                                         Double_t intRadius,
120 >                                         const PFCandidateCol *PFCands);                                                                                  
121 >                                        
122 >      static Float_t PFChargedCount(const mithep::Photon*,
123 >                                const BaseVertex*,
124 >                                Double_t extRadius,
125 >                                Double_t intRadius,
126 >                                Double_t ptLow,
127 >                                Double_t etaStrip,
128 >                                Double_t maxD0,
129 >                                Double_t maxDZ,
130 >                                const PFCandidateCol *PFCands,
131 >                                unsigned int* worstVtxIndex = NULL,
132 >                                const mithep::Collection<mithep::Vertex> *vtxs = NULL,
133 >                                const mithep::Collection<mithep::Electron> *eles = NULL,
134 >                                bool print = NULL,
135 >                                double* ptmax = NULL,
136 >                                double* dRmax = NULL);
137 >      
138 >      ClassDef(IsolationTools, 0) // Isolation tools
139 >        };
140   }
141   #endif

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines