Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • classes/DelphesClasses.cc

    r9cc5aeb r3051ea17  
    102102{
    103103  TLorentzVector vec;
    104   vec.SetPtEtaPhiM(PT, Eta, Phi, Mass);
    105   return vec;
    106 }
    107 
    108 //------------------------------------------------------------------------------
    109 
    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, 2.);
    117   Cv(1, 1)=TMath::Power(ErrorPhi, 2.);
    118   Cv(2, 2)=TMath::Power(ErrorC, 2.);
    119   Cv(3, 3)=TMath::Power(ErrorDZ, 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 
     104  vec.SetPtEtaPhiM(PT, Eta, Phi, 0.0);
     105  return vec;
     106}
    148107
    149108//------------------------------------------------------------------------------
     
    161120{
    162121  TLorentzVector vec;
    163   vec.SetPtEtaPhiM(PT, Eta, Phi, Mass);
    164   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, 2.);
    176   Cv(1, 1)=TMath::Power(ErrorPhi, 2.);
    177   Cv(2, 2)=TMath::Power(ErrorC, 2.);
    178   Cv(3, 3)=TMath::Power(ErrorDZ, 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;
     122  vec.SetPtEtaPhiM(PT, Eta, Phi, 0.0);
     123  return vec;
    205124}
    206125
     
    234153  NNeutrals(0),
    235154  NeutralEnergyFraction(0),  // charged energy fraction
    236   ChargedEnergyFraction(0),  // neutral energy fraction
     155  ChargedEnergyFraction(0),  // neutral energy fraction 
    237156  Beta(0),
    238157  BetaStar(0),
Note: See TracChangeset for help on using the changeset viewer.