ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitAna/DataTree/interface/Electron.h
(Generate patch)

Comparing UserCode/MitAna/DataTree/interface/Electron.h (file contents):
Revision 1.32 by phedex, Tue Apr 7 15:55:36 2009 UTC vs.
Revision 1.43 by bendavid, Mon Oct 18 01:33:37 2010 UTC

# Line 20 | Line 20 | namespace mithep
20    class Electron : public ChargedParticle
21    {
22      public:
23 <      Electron() : fESuperClusterOverP(0), fESeedClusterOverPout(0), fDeltaEtaSuperClTrkAtVtx(0),
24 <                   fDeltaEtaSeedClTrkAtCalo(0), fDeltaPhiSuperClTrkAtVtx(0),
25 <                   fDeltaPhiSeedClTrkAtCalo(0), fHadronicOverEm(0), fNumberOfClusters(0), fE33(0),
26 <                   fE55(0), fCovEtaEta(0), fCoviEtaiEta(0), fCovEtaPhi(0), fCovPhiPhi(0),
27 <                   fCaloIsolation(0), fCaloTowerIsolation(0), fTrackIsolation(0),
28 <                   fEcalJurassicIsolation(0), fHcalJurassicIsolation(0), fPassLooseID(0),
29 <                   fPassTightID(0), fIDLikelihood(0), fPIn(0), fPOut(0),
30 <                   fIsEnergyScaleCorrected(0), fIsMomentumCorrected(0), fClassification(0) {}
23 >      Electron() :
24 >        fCharge(-99), fScPixCharge(0),
25 >        fESuperClusterOverP(0), fESeedClusterOverPout(0), fDeltaEtaSuperClTrkAtVtx(0),
26 >        fDeltaEtaSeedClTrkAtCalo(0), fDeltaPhiSuperClTrkAtVtx(0),
27 >        fDeltaPhiSeedClTrkAtCalo(0), fFBrem(0), fHadronicOverEm(0), fHcalDepth1OverEcal(0),
28 >        fHcalDepth2OverEcal(0), fNumberOfClusters(0), fE15(0), fE25Max(0),
29 >        fE55(0), fCovEtaEta(0), fCoviEtaiEta(0),
30 >        fCaloIsolation(0), fHcalJurassicIsolation(0),
31 >        fHcalDepth1TowerSumEtDr04(0), fHcalDepth2TowerSumEtDr04(0),
32 >        fEcalJurassicIsolation(0), fTrackIsolationDr04(0), fCaloTowerIsolation(0),
33 >        fHcalDepth1TowerSumEtDr03(0), fHcalDepth2TowerSumEtDr03(0),
34 >        fEcalRecHitSumEtDr03(0), fTrackIsolation(0), fPassLooseID(0),
35 >        fPassTightID(0), fIDLikelihood(0), fPIn(0), fPOut(0), fFracSharedHits(0),
36 >        fMva(0), fD0PV(0), fD0PVErr(0), fIp3dPV(0), fIp3dPVErr(0),
37 >        fD0PVBS(0), fD0PVBSErr(0), fIp3dPVBS(0), fIp3dPVBSErr(0),
38 >        fGsfPVCompatibility(0), fGsfPVBSCompatibility(0),
39 >        fGsfPVCompatibilityMatched(0), fGsfPVBSCompatibilityMatched(0),
40 >        fConvPartnerDCotTheta(0), fConvPartnerDist(0), fConvPartnerRadius(0),
41 >        fIsEnergyScaleCorrected(0), fIsMomentumCorrected(0),
42 >        fClassification(0), fIsEB(), fIsEE(0), fIsEBEEGap(0), fIsEBEtaGap(0),
43 >        fIsEBPhiGap(0), fIsEEDeeGap(0), fIsEERingGap(0),
44 >        fIsEcalDriven(0), fIsTrackerDriven(0), fMatchesVertexConversion(0) {}
45  
46        const Track         *BestTrk()                const;
47 <      Double_t             CaloIsolation()          const { return fCaloIsolation;                 }
48 <      Double_t             CaloTowerIsolation()     const { return fCaloTowerIsolation;            }
49 <      Int_t                Classification()         const { return fClassification;                }
50 <      Double_t             CovEtaEta()              const { return fCovEtaEta;                     }
51 <      Double_t             CovEtaPhi()              const { return fCovEtaPhi;                     }
52 <      Double_t             CovPhiPhi()              const { return fCovPhiPhi;                     }
53 <      Double_t             CoviEtaiEta()            const { return fCoviEtaiEta;                   }
47 >      Double_t             D0PV()                   const { return fD0PV;                     }
48 >      Double_t             D0PVErr()                const { return fD0PVErr;                  }
49 >      Double_t             D0PVSignificance()       const { return fD0PV/fD0PVErr;            }
50 >      Double_t             Ip3dPV()                 const { return fIp3dPV;                   }
51 >      Double_t             Ip3dPVErr()              const { return fIp3dPVErr;                }
52 >      Double_t             Ip3dPVSignificance()     const { return fIp3dPV/fIp3dPVErr;        }
53 >      Double_t             D0PVBS()                 const { return fD0PVBS;                   }
54 >      Double_t             D0PVBSErr()              const { return fD0PVBSErr;                }
55 >      Double_t             D0PVBSSignificance()     const { return fD0PVBS/fD0PVBSErr;        }
56 >      Double_t             Ip3dPVBS()               const { return fIp3dPVBS;                 }
57 >      Double_t             Ip3dPVBSErr()            const { return fIp3dPVBSErr;              }
58 >      Double_t             Ip3dPVBSSignificance()   const { return fIp3dPVBS/fIp3dPVBSErr;    }
59 >      Double_t             GsfPVCompatibility()     const { return fGsfPVCompatibility;       }
60 >      Double_t             GsfPVBSCompatibility()   const { return fGsfPVBSCompatibility;     }
61 >      Double_t             GsfPVCompatibilityMatched()     const { return fGsfPVCompatibilityMatched;       }
62 >      Double_t             GsfPVBSCompatibilityMatched()   const { return fGsfPVBSCompatibilityMatched;     }
63 >      Double_t             ConvPartnerDCotTheta()   const { return fConvPartnerDCotTheta;     }
64 >      Double_t             ConvPartnerDist()        const { return fConvPartnerDist;          }
65 >      Double_t             ConvPartnerRadius()      const { return fConvPartnerRadius;        }
66 >      Double_t             CaloIsolation()          const { return fCaloIsolation;       } // *DEPRECATED*
67 >      Int_t                Classification()         const { return fClassification;      }
68 >      Double_t             CovEtaEta()              const { return fCovEtaEta;           }
69 >      Double_t             CoviEtaiEta()            const { return fCoviEtaiEta;         }
70        Double_t             DeltaEtaSuperClusterTrackAtVtx() const
71                               { return fDeltaEtaSuperClTrkAtVtx; }
72        Double_t             DeltaEtaSeedClusterTrackAtCalo() const
# Line 45 | Line 75 | namespace mithep
75                               { return fDeltaPhiSuperClTrkAtVtx; }
76        Double_t             DeltaPhiSeedClusterTrackAtCalo() const
77                               { return fDeltaPhiSeedClTrkAtCalo; }
78 <      Double_t             E33()                    const { return fE33;                           }
79 <      Double_t             E55()                    const { return fE55;                           }
80 <      Double_t             EcalJurassicIsolation()  const { return fEcalJurassicIsolation;         }
78 >      Double_t             E15()                    const { return fE15;                }
79 >      Double_t             E25Max()                 const { return fE25Max;             }
80 >      Double_t             E55()                    const { return fE55;                }
81        Double_t             ESuperClusterOverP()     const { return fESuperClusterOverP;            }
82        Double_t             ESeedClusterOverPout()   const { return fESeedClusterOverPout;          }
83        Double_t             ESeedClusterOverPIn()    const;
84 <      Double_t             FBrem()                  const { return (PIn() - POut())/PIn();   }
84 >      Double_t             FBrem()                  const { return fFBrem;                         }
85 >      Double_t             FBremOld()               const { return (PIn() - POut())/PIn();         }
86 >      Double_t             FracSharedHits()         const { return fFracSharedHits;                }
87        const Track         *GsfTrk()                 const { return fGsfTrackRef.Obj();             }
88        Double_t             HadronicOverEm()         const { return fHadronicOverEm;                }
89 +      Double_t             HcalDepth1OverEcal()     const { return fHcalDepth1OverEcal;            }
90 +      Double_t             HcalDepth2OverEcal()     const { return fHcalDepth2OverEcal;            }
91        Bool_t               HasGsfTrk()              const { return fGsfTrackRef.IsValid();         }
92        Bool_t               HasTrackerTrk()          const { return fTrackerTrackRef.IsValid();     }
93        Bool_t               HasSuperCluster()        const { return fSuperClusterRef.IsValid();     }
94 <      Double_t             HcalIsolation()          const { return fHcalJurassicIsolation;         }
94 >      Double_t             HcalIsolation()          const { return fHcalJurassicIsolation;         } // *DEPRECATED*
95        Double_t             IDLikelihood()           const { return fIDLikelihood;                  }
96        Bool_t               IsEnergyScaleCorrected() const { return fIsEnergyScaleCorrected;        }
97        Bool_t               IsMomentumCorrected()    const { return fIsMomentumCorrected;           }
98 +      Bool_t               IsEB()                   const { return fIsEB;                          }
99 +      Bool_t               IsEE()                   const { return fIsEE;                          }
100 +      Bool_t               IsEBEEGap()              const { return fIsEBEEGap;                     }
101 +      Bool_t               IsEBEtaGap()             const { return fIsEBEtaGap;                    }
102 +      Bool_t               IsEBPhiGap()             const { return fIsEBPhiGap;                    }
103 +      Bool_t               IsEEDeeGap()             const { return fIsEEDeeGap;                    }
104 +      Bool_t               IsEERingGap()            const { return fIsEERingGap;                   }
105 +      Bool_t               IsEcalDriven()           const { return fIsEcalDriven;                  }
106 +      Bool_t               IsTrackerDriven()        const { return fIsTrackerDriven;               }
107 +      Double_t             Mva()                    const { return fMva;                           }
108        Double_t             NumberOfClusters()       const { return fNumberOfClusters;              }
109        EObjType             ObjType()                const { return kElectron;                      }
110        Double_t             PassLooseID()            const { return fPassLooseID;                   }
# Line 68 | Line 112 | namespace mithep
112        Double_t             PIn()                    const { return fPIn;                           }
113        Double_t             POut()                   const { return fPOut;                          }
114        const SuperCluster  *SCluster()               const { return fSuperClusterRef.Obj();         }
115 <      void                 SetCaloIsolation(Double_t caloiso)     { fCaloIsolation = caloiso;      }
116 <      void                 SetCaloTowerIsolation(Double_t tiso)   { fCaloTowerIsolation = tiso;    }
115 >      Double_t             ScPixCharge()            const { return fScPixCharge;                   }
116 >      
117 >      Double_t             EcalRecHitIsoDr04()      const { return fEcalJurassicIsolation;         }
118 >      Double_t             HcalTowerSumEtDr04()     const { return HcalDepth1TowerSumEtDr04() +
119 >                                                                   HcalDepth2TowerSumEtDr04();     }
120 >      Double_t             HcalDepth1TowerSumEtDr04() const { return fHcalDepth1TowerSumEtDr04;    }
121 >      Double_t             HcalDepth2TowerSumEtDr04() const { return fHcalDepth2TowerSumEtDr04;    }
122 >      Double_t             TrackIsolationDr04()     const { return fTrackIsolationDr04;            }
123 >      Double_t             EcalRecHitIsoDr03()      const { return fEcalRecHitSumEtDr03;           }
124 >      Double_t             HcalTowerSumEtDr03()     const { return fCaloTowerIsolation;            }
125 >      Double_t             HcalDepth1TowerSumEtDr03() const { return fHcalDepth1TowerSumEtDr03;    }
126 >      Double_t             HcalDepth2TowerSumEtDr03() const { return fHcalDepth2TowerSumEtDr03;    }
127 >      Double_t             TrackIsolationDr03()     const { return fTrackIsolation;                }
128 >      Bool_t               MatchesVertexConversion() const { return fMatchesVertexConversion;      }
129 >      UInt_t               NAmbiguousGsfTracks()     const { return fAmbiguousGsfTracks.Entries();  }
130 >      Bool_t               HasAmbiguousGsfTrack(const Track *t) const { return fAmbiguousGsfTracks.HasObject(t); }
131 >      const Track         *AmbiguousGsfTrack(UInt_t i) const { return fAmbiguousGsfTracks.At(i);   }
132 >      
133 >      void                 AddAmbiguousGsfTrack(const Track *t)   { fAmbiguousGsfTracks.Add(t);    }
134 >      void                 SetCharge(Char_t x)                    { fCharge = x; ClearCharge();    }
135 >      void                 SetScPixCharge(Char_t x)               { fScPixCharge = x;              }
136 >      void                 SetD0PV(Double_t x)                    { fD0PV = x;                     }
137 >      void                 SetD0PVErr(Double_t x)                 { fD0PVErr = x;                  }
138 >      void                 SetIp3dPV(Double_t x)                  { fIp3dPV = x;                   }
139 >      void                 SetIp3dPVErr(Double_t x)               { fIp3dPVErr = x;                }
140 >      void                 SetD0PVBS(Double_t x)                  { fD0PVBS = x;                   }
141 >      void                 SetD0PVBSErr(Double_t x)               { fD0PVBSErr = x;                }
142 >      void                 SetIp3dPVBS(Double_t x)                { fIp3dPVBS = x;                 }
143 >      void                 SetIp3dPVBSErr(Double_t x)             { fIp3dPVBSErr = x;              }
144 >      void                 SetGsfPVCompatibility(Double_t x)      { fGsfPVCompatibility = x;       }
145 >      void                 SetGsfPVBSCompatibility(Double_t x)    { fGsfPVBSCompatibility = x;     }
146 >      void                 SetGsfPVCompatibilityMatched(Double_t x)      { fGsfPVCompatibilityMatched = x;   }
147 >      void                 SetGsfPVBSCompatibilityMatched(Double_t x)    { fGsfPVBSCompatibilityMatched = x; }
148 >      void                 SetConvPartnerDCotTheta(Double_t x)    { fConvPartnerDCotTheta = x;     }
149 >      void                 SetConvPartnerDist(Double_t x)         { fConvPartnerDist = x;          }
150 >      void                 SetConvPartnerRadius(Double_t x)       { fConvPartnerRadius = x;        }
151        void                 SetClassification(Int_t x)             { fClassification = x;           }
152        void                 SetCovEtaEta(Double_t CovEtaEta)       { fCovEtaEta = CovEtaEta;        }
75      void                 SetCovEtaPhi(Double_t CovEtaPhi)       { fCovEtaPhi = CovEtaPhi;        }
76      void                 SetCovPhiPhi(Double_t CovPhiPhi)       { fCovPhiPhi = CovPhiPhi;        }
153        void                 SetCoviEtaiEta(Double_t CoviEtaiEta)   { fCoviEtaiEta = CoviEtaiEta;    }
154        void                 SetDeltaEtaSuperClusterTrackAtVtx(Double_t x)  
155                               { fDeltaEtaSuperClTrkAtVtx = x;   }
# Line 83 | Line 159 | namespace mithep
159                               { fDeltaPhiSuperClTrkAtVtx = x;   }
160        void                 SetDeltaPhiSeedClusterTrackAtCalo(Double_t x)  
161                               { fDeltaPhiSeedClTrkAtCalo = x;   }
162 <      void                 SetE33(Double_t E33)                   { fE33 = E33;                    }
163 <      void                 SetE55(Double_t E55)                   { fE55 = E55;                    }
162 >      void                 SetE15(Double_t x)                     { fE15 = x;                      }
163 >      void                 SetE25Max(Double_t x)                  { fE25Max = x;                   }
164 >      void                 SetE55(Double_t x)                     { fE55 = x;                      }
165        void                 SetESeedClusterOverPout(Double_t x)    { fESeedClusterOverPout = x;     }
166        void                 SetESuperClusterOverP(Double_t x)      { fESuperClusterOverP = x;       }
167 <      void                 SetEcalJurassicIso(Double_t iso )      { fEcalJurassicIsolation = iso;  }
167 >      void                 SetFBrem(Double_t x)                   { fFBrem = x;                    }
168 >      void                 SetFracSharedHits(Double_t x)          { fFracSharedHits = x;           }
169        void                 SetGsfTrk(const Track* t)                    
170                               { fGsfTrackRef = t; ClearCharge(); }
171        void                 SetHadronicOverEm(Double_t x)          { fHadronicOverEm = x;           }
172 <      void                 SetHcalIsolation(Double_t iso )        { fHcalJurassicIsolation = iso;  }
172 >      void                 SetHcalDepth1OverEcal(Double_t x)      { fHcalDepth1OverEcal = x;       }
173 >      void                 SetHcalDepth2OverEcal(Double_t x)      { fHcalDepth2OverEcal = x;       }
174        void                 SetIDLikelihood(Double_t likelihood)   { fIDLikelihood = likelihood;    }
175        void                 SetIsEnergyScaleCorrected(Bool_t x)    { fIsEnergyScaleCorrected = x;   }
176        void                 SetIsMomentumCorrected(Bool_t x)       { fIsMomentumCorrected = x;      }
# Line 105 | Line 184 | namespace mithep
184                               { fSuperClusterRef = sc; }
185        void                 SetTrackerTrk(const Track* t)                
186                               { fTrackerTrackRef = t; ClearCharge(); }
187 <      void                 SetTrackIsolation(Double_t trkiso)     { fTrackIsolation = trkiso;      }
188 <      const Track         *TrackerTrk()            const { return fTrackerTrackRef.Obj();          }
189 <      Double_t             TrackIsolation()        const { return fTrackIsolation;                 }
190 <      const Track         *Trk()                   const { return BestTrk();                       }
187 >      void                 SetConvPartnerTrk(const Track *t)
188 >                             { fConvPartnerTrackRef = t; }
189 >      void                 SetEcalRecHitIsoDr04(Double_t x)        { fEcalJurassicIsolation = x;    }
190 >      void                 SetHcalDepth1TowerSumEtDr04(Double_t x) { fHcalDepth1TowerSumEtDr04 = x; }
191 >      void                 SetHcalDepth2TowerSumEtDr04(Double_t x) { fHcalDepth2TowerSumEtDr04 = x; }
192 >      void                 SetTrackIsolationDr04(Double_t x)       { fTrackIsolationDr04 = x;       }
193 >      void                 SetEcalRecHitIsoDr03(Double_t x)        { fEcalRecHitSumEtDr03 = x;      }
194 >      void                 SetHcalTowerSumEtDr03(Double_t x)       { fCaloTowerIsolation = x;       }
195 >      void                 SetHcalDepth1TowerSumEtDr03(Double_t x) { fHcalDepth1TowerSumEtDr03 = x; }
196 >      void                 SetHcalDepth2TowerSumEtDr03(Double_t x) { fHcalDepth2TowerSumEtDr03 = x; }
197 >      void                 SetTrackIsolationDr03(Double_t x)       { fTrackIsolation = x;   }
198 >      void                 SetMva(Double_t x)                      { fMva = x;              }
199 >      void                 SetIsEB(Bool_t b)                       { fIsEB = b;             }
200 >      void                 SetIsEE(Bool_t b)                       { fIsEE = b;             }
201 >      void                 SetIsEBEEGap(Bool_t b)                  { fIsEBEEGap = b;        }
202 >      void                 SetIsEBEtaGap(Bool_t b)                 { fIsEBEtaGap = b;       }
203 >      void                 SetIsEBPhiGap(Bool_t b)                 { fIsEBPhiGap = b;       }
204 >      void                 SetIsEEDeeGap(Bool_t b)                 { fIsEEDeeGap = b;       }
205 >      void                 SetIsEERingGap(Bool_t b)                { fIsEERingGap = b;      }
206 >      void                 SetIsEcalDriven(Bool_t b)               { fIsEcalDriven = b;     }
207 >      void                 SetIsTrackerDriven(Bool_t b)            { fIsTrackerDriven = b;  }
208 >      void                 SetMatchesVertexConversion(Bool_t b)    { fMatchesVertexConversion = b; }
209 >      void                 SetConversionXYZ(Double_t x, Double_t y, Double_t z)
210 >                                  { fConvPosition.SetXYZ(x,y,z); }
211 >      
212 >    
213 >      const Track         *TrackerTrk()            const { return fTrackerTrackRef.Obj();   }
214 >      const Track         *Trk()                   const { return BestTrk();                }
215 >      const Track         *ConvPartnerTrk()        const { return fConvPartnerTrackRef.Obj(); }
216  
217      protected:
218 <      Double_t             GetMass()               const          { return 0.51099892e-3;          }
218 >      Double_t             GetCharge()             const;
219 >      Double_t             GetMass()               const { return 0.51099892e-3;            }
220        void                 GetMom()                const;
221  
222        Vect3C               fMom;                       //stored three-momentum
223 +      Char_t               fCharge;                    //stored charge - filled with -99 when reading old files
224 +      Char_t               fScPixCharge;               //charge from supercluster-pixel matching
225        Ref<Track>           fGsfTrackRef;               //gsf track reference
226        Ref<Track>           fTrackerTrackRef;           //tracker track reference
227 +      Ref<Track>           fConvPartnerTrackRef;       //conversion partner track reference
228        Ref<SuperCluster>    fSuperClusterRef;           //reference to SuperCluster
229        Double32_t           fESuperClusterOverP;        //[0,0,14]super cluster e over p ratio
230        Double32_t           fESeedClusterOverPout;      //[0,0,14]seed cluster e over p mom
# Line 124 | Line 232 | namespace mithep
232        Double32_t           fDeltaEtaSeedClTrkAtCalo;   //[0,0,14]delta eta of seeed cluster with trk
233        Double32_t           fDeltaPhiSuperClTrkAtVtx;   //[0,0,14]delta phi of super cluster with trk
234        Double32_t           fDeltaPhiSeedClTrkAtCalo;   //[0,0,14]delta phi of seeed cluster with trk
235 <      Double32_t           fHadronicOverEm;            //[0,0,14]hadronic over em fraction
235 >      Double32_t           fFBrem;                     //[0,0,14]brem fraction
236 >      Double32_t           fHadronicOverEm;            //[0,0,14]hadronic over em fraction *DEPRECATED*
237 >      Double32_t           fHcalDepth1OverEcal;        //[0,0,14]hadronic over em fraction depth1
238 >      Double32_t           fHcalDepth2OverEcal;        //[0,0,14]hadronic over em fraction depth2
239        Double32_t           fNumberOfClusters;          //[0,0,14]number of associated clusters
240 <      Double32_t           fE33;                       //[0,0,14]3x3 crystal energy
240 >      Double32_t           fE15;                       //[0,0,14]1x5 crystal energy
241 >      Double32_t           fE25Max;                    //[0,0,14]2x5 crystal energy (max of two possible sums)
242        Double32_t           fE55;                       //[0,0,14]5x5 crystal energy
243        Double32_t           fCovEtaEta;                 //[0,0,14]variance eta-eta
244        Double32_t           fCoviEtaiEta;               //[0,0,14]covariance eta-eta (in crystals)
245 <      Double32_t           fCovEtaPhi;                 //[0,0,14]covariance eta-phi
246 <      Double32_t           fCovPhiPhi;                 //[0,0,14]covariance phi-phi
247 <      Double32_t           fCaloIsolation;             //[0,0,14]isolation based on rechits
248 <      Double32_t           fCaloTowerIsolation;        //[0,0,14]isolation based on calo towers
249 <      Double32_t           fTrackIsolation;            //[0,0,14]isolation based on tracks
250 <      Double32_t           fEcalJurassicIsolation;     //[0,0,14]ecal jura iso
251 <      Double32_t           fHcalJurassicIsolation;     //[0,0,14]hcal jura iso
245 >      Double32_t           fCaloIsolation;             //[0,0,14](non-jura) ecal isolation based on rechits dR 0.3 *DEPRECATED*
246 >      Double32_t           fHcalJurassicIsolation;     //[0,0,14]hcal jura iso dR 0.4 *DEPRECATED*
247 >      Double32_t           fHcalDepth1TowerSumEtDr04;  //[0,0,14]hcal depth1 tower based isolation dR 0.4
248 >      Double32_t           fHcalDepth2TowerSumEtDr04;  //[0,0,14]hcal depth2 tower based isolation dR 0.4
249 >      Double32_t           fEcalJurassicIsolation;     //[0,0,14]ecal jura iso dR 0.4 *RENAMING*
250 >      Double32_t           fTrackIsolationDr04;        //[0,0,14]isolation based on tracks dR 0.4
251 >      Double32_t           fCaloTowerIsolation;        //[0,0,14]hcal tower based isolation dR 0.3 *DEPRECATED*
252 >      Double32_t           fHcalDepth1TowerSumEtDr03;  //[0,0,14]hcal depth1 tower based isolation dR 0.3
253 >      Double32_t           fHcalDepth2TowerSumEtDr03;  //[0,0,14]hcal depth2 tower based isolation dR 0.3
254 >      Double32_t           fEcalRecHitSumEtDr03;       //[0,0,14]ecal jura iso dR 0.3
255 >      Double32_t           fTrackIsolation;            //[0,0,14]isolation based on tracks dR 0.3 *RENAMING*
256        Double32_t           fPassLooseID;               //[0,0,14]pass loose id
257        Double32_t           fPassTightID;               //[0,0,14]pass tight id
258        Double32_t           fIDLikelihood;              //[0,0,14]likelihood value
259        Double32_t           fPIn;                       //[0,0,14]momentum at vtx
260        Double32_t           fPOut;                      //[0,0,14]momentum at ecal surface
261 +      Double32_t           fFracSharedHits;            //[0,0,14]fraction of shared hits btw gsf and std. track
262 +      Double32_t           fMva;                       //[0,0,14] pflow mva output
263 +      Double32_t           fD0PV;                      //[0,0,14]transverse impact parameter to signal PV
264 +      Double32_t           fD0PVErr;                   //[0,0,14]transverse impact parameter uncertainty to signal PV
265 +      Double32_t           fIp3dPV;                    //[0,0,14]3d impact parameter to signal PV
266 +      Double32_t           fIp3dPVErr;                 //[0,0,14]3d impact parameter uncertainty to signal PV
267 +      Double32_t           fD0PVBS;                    //[0,0,14]transverse impact parameter to signal PV w/ bs constraint
268 +      Double32_t           fD0PVBSErr;                 //[0,0,14]transverse impact parameter uncertainty to signal PV w/ bs constraint
269 +      Double32_t           fIp3dPVBS;                  //[0,0,14]3d impact parameter to signal PV w/ bs constraint
270 +      Double32_t           fIp3dPVBSErr;               //[0,0,14]3d impact parameter uncertainty to signal PV w/ bs constraint
271 +      Double32_t           fGsfPVCompatibility;        //[0,0,14]gsf compatibility with signal PV
272 +      Double32_t           fGsfPVBSCompatibility;      //[0,0,14]gsf compatibility with signal PV w/ bs constraint
273 +      Double32_t           fGsfPVCompatibilityMatched; //[0,0,14]gsf compatibility with signal PV (matching ckf track excluded from vertex)
274 +      Double32_t           fGsfPVBSCompatibilityMatched; //[0,0,14]gsf compatibility with signal PV w/ bs constraint (matching ckf track excluded from vertex)
275 +      Double32_t           fConvPartnerDCotTheta;      //[0,0,14]delta cot theta to nearest conversion partner track
276 +      Double32_t           fConvPartnerDist;           //[0,0,14]distance in x-y plane to nearest conversion partner track
277 +      Double32_t           fConvPartnerRadius;         //[0,0,14]radius of helix intersection with conversion partner track
278 +      Vect3C               fConvPosition;
279        Bool_t               fIsEnergyScaleCorrected;    //class dependent escale correction
280        Bool_t               fIsMomentumCorrected;       //class dependent E-p combination
281        Int_t                fClassification;            //classification (see GsfElectron.h)
282 +      Bool_t               fIsEB;                      //is ECAL barrel
283 +      Bool_t               fIsEE;                      //is ECAL Endcap
284 +      Bool_t               fIsEBEEGap;                 //is in barrel-endcap gap
285 +      Bool_t               fIsEBEtaGap;                //is in EB eta module gap
286 +      Bool_t               fIsEBPhiGap;                //is in EB phi module gap
287 +      Bool_t               fIsEEDeeGap;                //is in EE dee gap
288 +      Bool_t               fIsEERingGap;               //is in EE ring gap
289 +      Bool_t               fIsEcalDriven;              //is std. egamma electron
290 +      Bool_t               fIsTrackerDriven;           //is pflow track-seeded electron
291 +      Bool_t               fMatchesVertexConversion;
292 +      RefArray<Track>      fAmbiguousGsfTracks;        //ambiguous gsf tracks for this electron
293  
294 <    ClassDef(Electron, 1) // Electron class
294 >    ClassDef(Electron, 8) // Electron class
295    };
296   }
297  
# Line 164 | Line 309 | inline const mithep::Track *mithep::Elec
309   }
310  
311   //--------------------------------------------------------------------------------------------------
312 + inline Double_t mithep::Electron::GetCharge() const
313 + {
314 +  // Return stored charge, unless it is set to invalid (-99),
315 +  // in that case get charge from track as before
316 +
317 +  if (fCharge==-99)
318 +    return mithep::ChargedParticle::GetCharge();
319 +  else
320 +    return fCharge;
321 +
322 + }
323 +
324 + //--------------------------------------------------------------------------------------------------
325   inline void mithep::Electron::GetMom() const
326   {
327    // Get momentum of the electron. We use an explicitly stored three vector, with the pdg mass,

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines