+ if(fTRDnSlices)
+ delete[] fTRDslices;
+}
+
+AliESDtrack &AliESDtrack::operator=(const AliESDtrack &source){
+
+
+ if(&source == this) return *this;
+ AliExternalTrackParam::operator=(source);
+
+
+ if(source.fCp){
+ // we have the trackparam: assign or copy construct
+ if(fCp)*fCp = *source.fCp;
+ else fCp = new AliExternalTrackParam(*source.fCp);
+ }
+ else{
+ // no track param delete the old one
+ if(fCp)delete fCp;
+ fCp = 0;
+ }
+
+ if(source.fIp){
+ // we have the trackparam: assign or copy construct
+ if(fIp)*fIp = *source.fIp;
+ else fIp = new AliExternalTrackParam(*source.fIp);
+ }
+ else{
+ // no track param delete the old one
+ if(fIp)delete fIp;
+ fIp = 0;
+ }
+
+
+ if(source.fTPCInner){
+ // we have the trackparam: assign or copy construct
+ if(fTPCInner) *fTPCInner = *source.fTPCInner;
+ else fTPCInner = new AliExternalTrackParam(*source.fTPCInner);
+ }
+ else{
+ // no track param delete the old one
+ if(fTPCInner)delete fTPCInner;
+ fTPCInner = 0;
+ }
+
+
+ if(source.fOp){
+ // we have the trackparam: assign or copy construct
+ if(fOp) *fOp = *source.fOp;
+ else fOp = new AliExternalTrackParam(*source.fOp);
+ }
+ else{
+ // no track param delete the old one
+ if(fOp)delete fOp;
+ fOp = 0;
+ }
+
+
+ if(source.fHMPIDp){
+ // we have the trackparam: assign or copy construct
+ if(fHMPIDp) *fHMPIDp = *source.fHMPIDp;
+ else fHMPIDp = new AliExternalTrackParam(*source.fHMPIDp);
+ }
+ else{
+ // no track param delete the old one
+ if(fHMPIDp)delete fHMPIDp;
+ fHMPIDp = 0;
+ }
+
+
+ // copy also the friend track
+ // use copy constructor
+ if(source.fFriendTrack){
+ // we have the trackparam: assign or copy construct
+ delete fFriendTrack; fFriendTrack=new AliESDfriendTrack(*source.fFriendTrack);
+ }
+ else{
+ // no track param delete the old one
+ delete fFriendTrack; fFriendTrack= 0;
+ }
+
+ fTPCClusterMap = source.fTPCClusterMap;
+ fTPCSharedMap = source.fTPCSharedMap;
+ // the simple stuff
+ fFlags = source.fFlags;
+ fID = source.fID;
+ fLabel = source.fLabel;
+ fITSLabel = source.fITSLabel;
+ for(int i = 0; i< 12;++i){
+ fITSModule[i] = source.fITSModule[i];
+ }
+ fTPCLabel = source.fTPCLabel;
+ fTRDLabel = source.fTRDLabel;
+ for(int i = 0; i< 3;++i){
+ fTOFLabel[i] = source.fTOFLabel[i];
+ }
+ fTOFCalChannel = source.fTOFCalChannel;
+ fTOFindex = source.fTOFindex;
+ fHMPIDqn = source.fHMPIDqn;
+ fHMPIDcluIdx = source.fHMPIDcluIdx;
+ fCaloIndex = source.fCaloIndex;
+
+ for(int i = 0; i< 3;++i){
+ fKinkIndexes[i] = source.fKinkIndexes[i];
+ fV0Indexes[i] = source.fV0Indexes[i];
+ }
+
+ for(int i = 0; i< AliPID::kSPECIES;++i){
+ fR[i] = source.fR[i];
+ fITSr[i] = source.fITSr[i];
+ fTPCr[i] = source.fTPCr[i];
+ fTRDr[i] = source.fTRDr[i];
+ fTOFr[i] = source.fTOFr[i];
+ fHMPIDr[i] = source.fHMPIDr[i];
+ fTrackTime[i] = source.fTrackTime[i];
+ }
+
+ fHMPIDtrkTheta = source.fHMPIDtrkTheta;
+ fHMPIDtrkPhi = source.fHMPIDtrkPhi;
+ fHMPIDsignal = source.fHMPIDsignal;
+
+
+ fTrackLength = source. fTrackLength;
+ fdTPC = source.fdTPC;
+ fzTPC = source.fzTPC;
+ fCddTPC = source.fCddTPC;
+ fCdzTPC = source.fCdzTPC;
+ fCzzTPC = source.fCzzTPC;
+ fCchi2TPC = source.fCchi2TPC;
+
+ fD = source.fD;
+ fZ = source.fZ;
+ fCdd = source.fCdd;
+ fCdz = source.fCdz;
+ fCzz = source.fCzz;
+ fCchi2 = source.fCchi2;
+
+ fITSchi2 = source.fITSchi2;
+ fTPCchi2 = source.fTPCchi2;
+ fTPCchi2Iter1 = source.fTPCchi2Iter1;
+ fTRDchi2 = source.fTRDchi2;
+ fTOFchi2 = source.fTOFchi2;
+ fHMPIDchi2 = source.fHMPIDchi2;
+
+ fGlobalChi2 = source.fGlobalChi2;
+
+ fITSsignal = source.fITSsignal;
+ for (Int_t i=0;i<4;i++) {fITSdEdxSamples[i]=source.fITSdEdxSamples[i];}
+ fTPCsignal = source.fTPCsignal;
+ fTPCsignalS = source.fTPCsignalS;
+ for(int i = 0; i< 4;++i){
+ fTPCPoints[i] = source.fTPCPoints[i];
+ }
+ fTRDsignal = source.fTRDsignal;
+
+ for(int i = 0;i < kTRDnPlanes;++i){
+ fTRDTimBin[i] = source.fTRDTimBin[i];
+ }
+
+ if(fTRDnSlices)
+ delete[] fTRDslices;
+ fTRDslices=0;
+ fTRDnSlices=source.fTRDnSlices;
+ if (fTRDnSlices) {
+ fTRDslices=new Double32_t[fTRDnSlices];
+ for(int j = 0;j < fTRDnSlices;++j) fTRDslices[j] = source.fTRDslices[j];
+ }
+
+ fTRDQuality = source.fTRDQuality;
+ fTRDBudget = source.fTRDBudget;
+ fTOFsignal = source.fTOFsignal;
+ fTOFsignalToT = source.fTOFsignalToT;
+ fTOFsignalRaw = source.fTOFsignalRaw;
+ fTOFsignalDz = source.fTOFsignalDz;
+ fTOFsignalDx = source.fTOFsignalDx;
+ fTOFdeltaBC = source.fTOFdeltaBC;
+ fTOFl0l1 = source.fTOFl0l1;
+
+ for(int i = 0;i<10;++i){
+ fTOFInfo[i] = source.fTOFInfo[i];
+ }
+
+ fHMPIDtrkX = source.fHMPIDtrkX;
+ fHMPIDtrkY = source.fHMPIDtrkY;
+ fHMPIDmipX = source.fHMPIDmipX;
+ fHMPIDmipY = source.fHMPIDmipY;
+
+ fTPCncls = source.fTPCncls;
+ fTPCnclsF = source.fTPCnclsF;
+ fTPCsignalN = source.fTPCsignalN;
+ fTPCnclsIter1 = source.fTPCnclsIter1;
+ fTPCnclsFIter1 = source.fTPCnclsFIter1;
+
+ fITSncls = source.fITSncls;
+ fITSClusterMap = source.fITSClusterMap;
+ fITSSharedMap = source.fITSSharedMap;
+ fTRDncls = source.fTRDncls;
+ fTRDncls0 = source.fTRDncls0;
+ fTRDntracklets = source.fTRDntracklets;
+ fVertexID = source.fVertexID;
+ return *this;
+}
+
+
+
+void AliESDtrack::Copy(TObject &obj) const {
+
+ // this overwrites the virtual TOBject::Copy()
+ // to allow run time copying without casting
+ // in AliESDEvent
+
+ if(this==&obj)return;
+ AliESDtrack *robj = dynamic_cast<AliESDtrack*>(&obj);
+ if(!robj)return; // not an AliESDtrack
+ *robj = *this;
+
+}
+
+
+
+void AliESDtrack::AddCalibObject(TObject * object){
+ //
+ // add calib object to the list
+ //
+ if (!fFriendTrack) fFriendTrack = new AliESDfriendTrack;
+ fFriendTrack->AddCalibObject(object);
+}
+
+TObject * AliESDtrack::GetCalibObject(Int_t index){
+ //
+ // return calib objct at given position
+ //
+ if (!fFriendTrack) return 0;
+ return fFriendTrack->GetCalibObject(index);
+}
+
+
+Bool_t AliESDtrack::FillTPCOnlyTrack(AliESDtrack &track){
+
+ // Fills the information of the TPC-only first reconstruction pass
+ // into the passed ESDtrack object. For consistency fTPCInner is also filled
+ // again
+
+
+
+ // For data produced before r26675
+ // RelateToVertexTPC was not properly called during reco
+ // so you'll have to call it again, before FillTPCOnlyTrack
+ // Float_t p[2],cov[3];
+ // track->GetImpactParametersTPC(p,cov);
+ // if(p[0]==0&&p[1]==0) // <- Default values
+ // track->RelateToVertexTPC(esd->GetPrimaryVertexTPC(),esd->GetMagneticField(),kVeryBig);
+
+
+ if(!fTPCInner)return kFALSE;
+
+ // fill the TPC track params to the global track parameters
+ track.Set(fTPCInner->GetX(),fTPCInner->GetAlpha(),fTPCInner->GetParameter(),fTPCInner->GetCovariance());
+ track.fD = fdTPC;
+ track.fZ = fzTPC;
+ track.fCdd = fCddTPC;
+ track.fCdz = fCdzTPC;
+ track.fCzz = fCzzTPC;
+
+ // copy the inner params
+ if(track.fIp) *track.fIp = *fIp;
+ else track.fIp = new AliExternalTrackParam(*fIp);
+
+ // copy the TPCinner parameters
+ if(track.fTPCInner) *track.fTPCInner = *fTPCInner;
+ else track.fTPCInner = new AliExternalTrackParam(*fTPCInner);
+ track.fdTPC = fdTPC;
+ track.fzTPC = fzTPC;
+ track.fCddTPC = fCddTPC;
+ track.fCdzTPC = fCdzTPC;
+ track.fCzzTPC = fCzzTPC;
+ track.fCchi2TPC = fCchi2TPC;
+
+ // copy all other TPC specific parameters
+
+ // replace label by TPC label
+ track.fLabel = fTPCLabel;
+ track.fTPCLabel = fTPCLabel;
+
+ track.fTPCchi2 = fTPCchi2;
+ track.fTPCchi2Iter1 = fTPCchi2Iter1;
+ track.fTPCsignal = fTPCsignal;
+ track.fTPCsignalS = fTPCsignalS;
+ for(int i = 0;i<4;++i)track.fTPCPoints[i] = fTPCPoints[i];
+
+ track.fTPCncls = fTPCncls;
+ track.fTPCnclsF = fTPCnclsF;
+ track.fTPCsignalN = fTPCsignalN;
+ track.fTPCnclsIter1 = fTPCnclsIter1;
+ track.fTPCnclsFIter1 = fTPCnclsFIter1;
+
+ // PID
+ for(int i=0;i<AliPID::kSPECIES;++i){
+ track.fTPCr[i] = fTPCr[i];
+ // combined PID is TPC only!
+ track.fR[i] = fTPCr[i];
+ }
+ track.fTPCClusterMap = fTPCClusterMap;
+ track.fTPCSharedMap = fTPCSharedMap;
+
+
+ // reset the flags
+ track.fFlags = kTPCin;
+ track.fID = fID;
+
+ track.fFlags |= fFlags & kTPCpid; //copy the TPCpid status flag
+
+ for (Int_t i=0;i<3;i++) track.fKinkIndexes[i] = fKinkIndexes[i];
+
+ return kTRUE;
+