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.14 by ceballos, Thu May 12 23:16:11 2011 UTC vs.
Revision 1.26 by mingyang, Tue May 22 23:47:02 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.2, 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.2 , 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*,
104 >                                         const BaseVertex*,
105 >                                         Double_t extRadius,
106 >                                         Double_t intRadius,
107 >                                         Double_t ptLow,
108 >                                         Double_t etaStrip,
109 >                                         Double_t maxD0,
110 >                                         Double_t maxDZ,
111 >                                         const PFCandidateCol *PFCands,
112 >                                         unsigned int* worstVtxIndex = NULL,
113 >                                         const mithep::Collection<mithep::Vertex> *vtxs = NULL,
114 >                                         const mithep::Collection<mithep::Electron> *eles = NULL,
115 >                                         bool print = NULL,
116 >                                         double* ptmax = NULL,
117 >                                         double* dRmax = NULL);
118 >
119 >      static Float_t PFChargedCount(const mithep::Photon*,
120 >                                const BaseVertex*,
121 >                                Double_t extRadius,
122 >                                Double_t intRadius,
123 >                                Double_t ptLow,
124 >                                Double_t etaStrip,
125 >                                Double_t maxD0,
126 >                                Double_t maxDZ,
127 >                                const PFCandidateCol *PFCands,
128 >                                unsigned int* worstVtxIndex = NULL,
129 >                                const mithep::Collection<mithep::Vertex> *vtxs = NULL,
130 >                                const mithep::Collection<mithep::Electron> *eles = NULL,
131 >                                bool print = NULL,
132 >                                double* ptmax = NULL,
133 >                                double* dRmax = NULL);
134 >      
135 >      ClassDef(IsolationTools, 0) // Isolation tools
136 >        };
137   }
138   #endif

Diff Legend

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