Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • external/TrackCovariance/ObsTrk.cc

    r00b14d5 rebf40fd  
    7474        fCovILC.ResizeTo(5, 5);
    7575        fGenPar = XPtoPar(fGenX, fGenP, Q);
    76         fGenParMm = ParToMm(fGenPar);
    7776        fGenParACTS = ParToACTS(fGenPar);
    7877        fGenParILC = ParToILC(fGenPar);
    7978        //
    8079        fObsPar = GenToObsPar(fGenPar);
    81         fObsParMm = ParToMm(fObsPar);
    8280        fObsParACTS = ParToACTS(fObsPar);
    8381        fObsParILC = ParToILC(fObsPar);
     
    8583        fObsP = ParToP(fObsPar);
    8684        fObsQ = ParToQ(fObsPar);
    87         fCovMm = CovToMm(fCov);
    8885        fCovACTS = CovToACTS(fObsPar, fCov);
    8986        fCovILC = CovToILC(fCov);
     
    134131        Double_t ZinNeg = fG->GetZminNeg();
    135132        Bool_t inside = TrkUtil::IsInside(fGenX, Rin, ZinNeg, ZinPos); // Check if in inner box
    136         SolTrack* trk = new SolTrack(fGenX, fGenP, fG);
    137         Double_t Xfirst, Yfirst, Zfirst;
    138         Int_t iLay = trk->FirstHit(Xfirst, Yfirst, Zfirst);
    139         TVector3 fXfirst(Xfirst, Yfirst, Zfirst);
    140133        if (inside)
    141134        {
     
    151144                //std::cout<<"ObsTrk:: outside: x= "<<fGenX(0)<<", y= "<<fGenX(1)
    152145                //                         <<", z= "<<fGenX(2)<<std::endl;
     146                SolTrack* trk = new SolTrack(fGenX, fGenP, fG);
    153147                Bool_t Res = kTRUE; Bool_t MS = kTRUE;
    154148                trk->CovCalc(Res, MS);                                  // Calculate covariance matrix
    155                 Cov = trk->Cov();
    156         }                                       // Track covariance
    157         delete trk;
     149                Cov = trk->Cov();                                       // Track covariance
     150                delete trk;
     151        }
    158152        //
    159153        fCov = Cov;
Note: See TracChangeset for help on using the changeset viewer.