+ // 2. TPC constrained umulation
+ AliExternalTrackParam *tdConstrained = new AliExternalTrackParam(*td1);
+ tdConstrained->Rotate(ot1->GetAlpha());
+ tdConstrained->PropagateTo(ot1->GetX(), AliTrackerBase::GetBz());
+ Double_t pointPos[2]={ot1->GetY(),ot1->GetZ()}; // local y and local z of point
+ Double_t pointCov[3]={0.0001,0,0.0001}; //
+ tdConstrained->Update(pointPos,pointCov);
+ // 3. TPC+ITS constrained umulation
+ AliExternalTrackParam *tdITS = new AliExternalTrackParam(*td0);
+ AliExternalTrackParam *tdITSOrig = new AliExternalTrackParam(*ot0);
+ Bool_t itsOK=kTRUE;
+ //
+ for (Int_t ilayer=6; ilayer<=0; ilayer--){
+ if (!AliTrackerBase::PropagateTrackTo(tdITSOrig,xITSlayer[ilayer],kMass,5.,kTRUE,kMaxSnp)) itsOK=kFALSE;
+ if (!AliTrackerBase::PropagateTrackTo(tdITS,xITSlayer[ilayer],kMass,5.,kTRUE,kMaxSnp)) itsOK=kFALSE;
+ //
+ tdITS->Rotate(tdITSOrig->GetAlpha());
+ tdITS->PropagateTo(tdITSOrig->GetX(), AliTrackerBase::GetBz());
+ Double_t itspointPos[2]={tdITS->GetY(),tdITS->GetZ()}; // local y and local z of point
+ Double_t itspointCov[3]={resITSlayer[ilayer]*resITSlayer[ilayer],0,resITSlayer[ilayer]*resITSlayer[ilayer]};
+ tdITS->Update(itspointPos,itspointCov);
+ }
+ //