Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • classes/DelphesClasses.cc

    r2671df6 r3051ea17  
    108108//------------------------------------------------------------------------------
    109109
    110 TMatrixDSym Track::CovarianceMatrix() const
    111 {
    112   TMatrixDSym Cv;
    113   Cv.ResizeTo(5, 5);
    114  
    115   // convert diagonal term to original units
    116   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 terms
    123   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 
    151110TLorentzVector Tower::P4() const
    152111{
     
    163122  vec.SetPtEtaPhiM(PT, Eta, Phi, 0.0);
    164123  return vec;
    165 }
    166 
    167 //------------------------------------------------------------------------------
    168 
    169 TMatrixDSym ParticleFlowCandidate::CovarianceMatrix() const
    170 {
    171   TMatrixDSym Cv;
    172   Cv.ResizeTo(5, 5);
    173  
    174   // convert diagonal term to original units
    175   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 terms
    182   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;
    205124}
    206125
Note: See TracChangeset for help on using the changeset viewer.