Changes in classes/DelphesClasses.cc [2671df6:3051ea17] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
classes/DelphesClasses.cc
r2671df6 r3051ea17 108 108 //------------------------------------------------------------------------------ 109 109 110 TMatrixDSym Track::CovarianceMatrix() const111 {112 TMatrixDSym Cv;113 Cv.ResizeTo(5, 5);114 115 // convert diagonal term to original units116 Cv(0, 0)=TMath::Power(ErrorD0*1.0E-3, 2.);117 Cv(1, 1)=TMath::Power(ErrorPhi, 2.);118 Cv(2, 2)=TMath::Power(ErrorC*1.0E-3, 2.);119 Cv(3, 3)=TMath::Power(ErrorDZ*1.0E-3, 2.);120 Cv(4, 4)=TMath::Power(ErrorCtgTheta, 2.);121 122 // off diagonal terms123 Cv(0, 1)=ErrorD0Phi;124 Cv(0, 2)=ErrorD0C;125 Cv(0, 3)=ErrorD0DZ;126 Cv(0, 4)=ErrorD0CtgTheta;127 Cv(1, 2)=ErrorPhiC;128 Cv(1, 3)=ErrorPhiDZ;129 Cv(1, 4)=ErrorPhiCtgTheta;130 Cv(2, 3)=ErrorCDZ;131 Cv(2, 4)=ErrorCCtgTheta;132 Cv(3, 4)=ErrorDZCtgTheta;133 134 Cv(1, 0)=Cv(0, 1);135 Cv(2, 0)=Cv(0, 2);136 Cv(3, 0)=Cv(0, 3);137 Cv(4, 0)=Cv(0, 4);138 Cv(2, 1)=Cv(1, 2);139 Cv(3, 1)=Cv(1, 3);140 Cv(4, 1)=Cv(1, 4);141 Cv(3, 2)=Cv(2, 3);142 Cv(4, 2)=Cv(2, 4);143 Cv(4, 3)=Cv(3, 4);144 145 return Cv;146 }147 148 149 //------------------------------------------------------------------------------150 151 110 TLorentzVector Tower::P4() const 152 111 { … … 163 122 vec.SetPtEtaPhiM(PT, Eta, Phi, 0.0); 164 123 return vec; 165 }166 167 //------------------------------------------------------------------------------168 169 TMatrixDSym ParticleFlowCandidate::CovarianceMatrix() const170 {171 TMatrixDSym Cv;172 Cv.ResizeTo(5, 5);173 174 // convert diagonal term to original units175 Cv(0, 0)=TMath::Power(ErrorD0*1.0E-3, 2.);176 Cv(1, 1)=TMath::Power(ErrorPhi, 2.);177 Cv(2, 2)=TMath::Power(ErrorC*1.0E-3, 2.);178 Cv(3, 3)=TMath::Power(ErrorDZ*1.0E-3, 2.);179 Cv(4, 4)=TMath::Power(ErrorCtgTheta, 2.);180 181 // off diagonal terms182 Cv(0, 1)=ErrorD0Phi;183 Cv(0, 2)=ErrorD0C;184 Cv(0, 3)=ErrorD0DZ;185 Cv(0, 4)=ErrorD0CtgTheta;186 Cv(1, 2)=ErrorPhiC;187 Cv(1, 3)=ErrorPhiDZ;188 Cv(1, 4)=ErrorPhiCtgTheta;189 Cv(2, 3)=ErrorCDZ;190 Cv(2, 4)=ErrorCCtgTheta;191 Cv(3, 4)=ErrorDZCtgTheta;192 193 Cv(1, 0)=Cv(0, 1);194 Cv(2, 0)=Cv(0, 2);195 Cv(3, 0)=Cv(0, 3);196 Cv(4, 0)=Cv(0, 4);197 Cv(2, 1)=Cv(1, 2);198 Cv(3, 1)=Cv(1, 3);199 Cv(4, 1)=Cv(1, 4);200 Cv(3, 2)=Cv(2, 3);201 Cv(4, 2)=Cv(2, 4);202 Cv(4, 3)=Cv(3, 4);203 204 return Cv;205 124 } 206 125
Note:
See TracChangeset
for help on using the changeset viewer.