Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/TrackCovariance.cc

    r3051ea17 r46c8df8  
    9696  Candidate *candidate, *mother;
    9797  Double_t mass, p, pt, q, ct;
    98   Double_t dd0, ddz, dphi, dct, dp, dpt, dC;
     98  Double_t dd0, ddz, dphi, dct, dp, dpt;
    9999 
    100100
     
    114114    candidate->Momentum.SetVectM(track.GetObsP(), mass);
    115115    candidate->InitialPosition.SetXYZT(track.GetObsX().X(),track.GetObsX().Y(),track.GetObsX().Z(),candidatePosition.T());
    116    
    117     // save full covariance 5x5 matrix internally (D0, phi, Curvature, dz, ctg(theta))
    118     candidate->TrackCovariance = track.GetCov();
    119116
    120117    pt = candidate->Momentum.Pt();
     
    128125   
    129126    candidate->D0 = track.GetObsPar()[0];
     127    candidate->DZ = track.GetObsPar()[3];
     128    candidate->P  = track.GetObsP().Mag();
     129    candidate->CtgTheta = track.GetObsPar()[4];
    130130    candidate->Phi = track.GetObsPar()[1];
    131     candidate->C = track.GetObsPar()[2];
    132     candidate->DZ = track.GetObsPar()[3];
    133     candidate->CtgTheta = track.GetObsPar()[4];
    134     candidate->P  = track.GetObsP().Mag();
     131
    135132    candidate->PT = pt;
    136133    candidate->Charge = q;
     
    142139    dpt       = 2 * TMath::Sqrt( track.GetCov()(2, 2))*pt*pt / (0.2998*fBz);
    143140    dp        = TMath::Sqrt((1.+ct*ct)*dpt*dpt + 4*pt*pt*ct*ct*dct*dct/(1.+ct*ct)/(1.+ct*ct));
    144     dC        = TMath::Sqrt(track.GetCov()(2, 2));
    145141
    146142    candidate->ErrorD0 = dd0;
    147143    candidate->ErrorDZ = ddz;
    148144    candidate->ErrorP = dp;
    149     candidate->ErrorC = dC;
    150145    candidate->ErrorCtgTheta = dct;
    151146    candidate->ErrorPhi = dphi;
     
    153148    //candidate->TrackResolution = dpt / pt;
    154149    candidate->TrackResolution = dp / p;
    155 
    156150
    157151    candidate->AddCandidate(mother);
Note: See TracChangeset for help on using the changeset viewer.