}
//_______________________________________________________________
-AliTRDtrackV1::AliTRDtrackV1(const AliTRDseedV1 * const trklts, const Double_t p[5], const Double_t cov[15]
+AliTRDtrackV1::AliTRDtrackV1(AliTRDseedV1 * const trklts, const Double_t p[5], const Double_t cov[15]
, Double_t x, Double_t alpha) : AliKalmanTrack()
,fStatus(0)
,fDE(0.)
fTrackletIndex[iplane] = 0xffff;
if(!trklts[iplane].IsOK()) fTracklet[iplane] = NULL;
else{
- fTracklet[iplane] = const_cast<AliTRDseedV1 *>(&trklts[iplane]);
+ fTracklet[iplane] = &trklts[iplane];
ncls += fTracklet[iplane]->GetN();
}
}
nPoints++;
}
}
- if(fitter->Eval()) return 1.E10;
+ fitter->Eval();
// Calculate curvature
Double_t a = fitter->GetParameter(0);
nPoints++;
}
}
- if(fitter->Eval()) return 1.E10;
+ fitter->Eval();
zfitter.Eval();
Double_t offset = fitter->GetParameter(3);
Double_t zmf = zfitter.GetFunctionValue(&xref);
fgTiltedRieman->FixParameter(3, zmf);
fgTiltedRieman->FixParameter(4, dzmf);
- if (fitter->Eval()) return 1.E10;
+ fitter->Eval();
fitter->ReleaseParameter(3);
fitter->ReleaseParameter(4);
offset = fitter->GetParameter(3);
Double_t zmf = zfitter.GetFunctionValue(&xref);
fitter->FixParameter(3, zmf);
fitter->FixParameter(4, dzmf);
- if(fitter->Eval()) return 1.E10;
+ fitter->Eval();
fitter->ReleaseParameter(3);
fitter->ReleaseParameter(4);
z0 = fitter->GetParameter(3); // = zmf ?
//____________________________________________________________________
-Double_t AliTRDtrackerV1::FitKalman(AliTRDtrackV1 *track, const AliTRDseedV1 * const tracklets, Bool_t up, Int_t np, AliTrackPoint *points)
+Double_t AliTRDtrackerV1::FitKalman(AliTRDtrackV1 *track, AliTRDseedV1 * const tracklets, Bool_t up, Int_t np, AliTrackPoint *points)
{
// Kalman filter implementation for the TRD.
// It returns the positions of the fit in the array "points"
// GET TRACKLET OR BUILT IT
Int_t iplane = up ? jplane : kNPlanes - 1 - jplane;
if(tracklets){
- if(!(ptrTracklet = const_cast<AliTRDseedV1 *>(&tracklets[iplane]))) continue;
+ if(!(ptrTracklet = &tracklets[iplane])) continue;
}else{
if(!(ptrTracklet = track->GetTracklet(iplane))){
/*AliTRDtrackerV1 *tracker = NULL;
}
//____________________________________________________________________
-Int_t AliTRDtrackerV1::MakeSeeds(AliTRDtrackingChamber **stack, const AliTRDseedV1 * const sseed, const Int_t * const ipar)
+Int_t AliTRDtrackerV1::MakeSeeds(AliTRDtrackingChamber **stack, AliTRDseedV1 * const sseed, const Int_t * const ipar)
{
//
// Seed tracklets and build candidate TRD tracks. The procedure is used during barrel tracking to account for tracks which are
AliTRDtrackingChamber *chamber = NULL;
AliTRDcluster *c[kNSeedPlanes] = {NULL, NULL, NULL, NULL}; // initilize seeding clusters
- AliTRDseedV1 *cseed = const_cast<AliTRDseedV1 *>(&sseed[0]); // initialize tracklets for first track
+ AliTRDseedV1 *cseed = &sseed[0]; // initialize tracklets for first track
Int_t ncl, mcl; // working variable for looping over clusters
Int_t index[AliTRDchamberTimeBin::kMaxClustersLayer], jndex[AliTRDchamberTimeBin::kMaxClustersLayer];
// chi2 storage
}
//_____________________________________________________________________________
-AliTRDtrackV1* AliTRDtrackerV1::MakeTrack(const AliTRDseedV1 * const seeds, Double_t *params)
+AliTRDtrackV1* AliTRDtrackerV1::MakeTrack(AliTRDseedV1 * const seeds, Double_t *params)
{
//
// Build a TRD track out of tracklet candidates
if(fkReconstructor->IsHLT()){
for (Int_t jLayer = 0; jLayer < AliTRDgeometry::kNlayer; jLayer++) {
track.UnsetTracklet(jLayer);
- ptrTracklet = const_cast<AliTRDseedV1 *>(&seeds[jLayer]);
+ ptrTracklet = &seeds[jLayer];
if(!ptrTracklet->IsOK()) continue;
//if(TMath::Abs(ptrTracklet->GetYref(1) - ptrTracklet->GetYfit(1)) >= .2) continue; // check this condition with Marian
ptrTracklet = SetTracklet(ptrTracklet);
static Double_t FitRiemanTilt(const AliTRDtrackV1 *trk, AliTRDseedV1 *tracklets = NULL, Bool_t err=0, Int_t np = 0, AliTrackPoint *points = NULL);
static Double_t FitLine(const AliTRDtrackV1 *trk, AliTRDseedV1 *tracklets = NULL, Bool_t err=0, Int_t np = 0, AliTrackPoint *points = NULL);
- static Double_t FitKalman(AliTRDtrackV1 *trk, const AliTRDseedV1 * const tracklets = NULL, Bool_t up=0, Int_t np = 0, AliTrackPoint *points = NULL);
+ static Double_t FitKalman(AliTRDtrackV1 *trk, AliTRDseedV1 * const tracklets = NULL, Bool_t up=0, Int_t np = 0, AliTrackPoint *points = NULL);
Bool_t IsClustersOwner() const { return TestBit(kOwner);}
Bool_t HasRemoveContainers() const { return TestBit(kRemoveContainers);}
AliTRDseedV1* GetTracklet(AliTRDtrackV1 *const trk, Int_t plane, Int_t &idx);
Bool_t GetTrackPoint(Int_t index, AliTrackPoint &p) const;
Float_t GetR4Layer(Int_t ly) const { return fR[ly];}
- Int_t MakeSeeds(AliTRDtrackingChamber **stack, const AliTRDseedV1 * const sseed, const Int_t * const ipar);
- AliTRDtrackV1* MakeTrack(const AliTRDseedV1 * const seeds, Double_t *params);
+ Int_t MakeSeeds(AliTRDtrackingChamber **stack, AliTRDseedV1 * const sseed, const Int_t * const ipar);
+ AliTRDtrackV1* MakeTrack(AliTRDseedV1 * const seeds, Double_t *params);
AliTRDtrackV1* SetTrack(const AliTRDtrackV1 * const track);
AliTRDseedV1* SetTracklet(const AliTRDseedV1 * const tracklet);