//-----------------------------------------------------------------
#include <TMath.h>
+#include <TParticle.h>
#include "AliESDVertex.h"
#include "AliESDtrack.h"
fOp(0),
fITSchi2(0),
fITSncls(0),
+ fITSClusterMap(0),
fITSsignal(0),
fITSLabel(0),
fTPCchi2(0),
fTOFCalChannel(-1),
fTOFsignal(-1),
fTOFsignalToT(0),
- fRICHchi2(1e10),
- fRICHqn(-1),
- fRICHcluIdx(-1),
- fRICHsignal(-1),
- fRICHtrkTheta(-1),
- fRICHtrkPhi(-1),
- fRICHtrkX(-1),
- fRICHtrkY(-1),
- fRICHmipX(-1),
- fRICHmipY(-1),
+ fHMPIDchi2(1e10),
+ fHMPIDqn(-1),
+ fHMPIDcluIdx(-1),
+ fHMPIDsignal(-1),
+ fHMPIDtrkTheta(-1),
+ fHMPIDtrkPhi(-1),
+ fHMPIDtrkX(-1),
+ fHMPIDtrkY(-1),
+ fHMPIDmipX(-1),
+ fHMPIDmipY(-1),
fEMCALindex(kEMCALNoMatch),
- fPoints(0),
fFriendTrack(new AliESDfriendTrack())
{
//
fTPCr[i]=1.;
fTRDr[i]=1.;
fTOFr[i]=1.;
- fRICHr[i]=1.;
+ fHMPIDr[i]=1.;
}
for (i=0; i<3; i++) { fKinkIndexes[i]=0;}
fOp(0),
fITSchi2(track.fITSchi2),
fITSncls(track.fITSncls),
+ fITSClusterMap(track.fITSClusterMap),
fITSsignal(track.fITSsignal),
fITSLabel(track.fITSLabel),
fTPCchi2(track.fTPCchi2),
fTOFCalChannel(track.fTOFCalChannel),
fTOFsignal(track.fTOFsignal),
fTOFsignalToT(track.fTOFsignalToT),
- fRICHchi2(track.fRICHchi2),
- fRICHqn(track.fRICHqn),
- fRICHcluIdx(track.fRICHcluIdx),
- fRICHsignal(track.fRICHsignal),
- fRICHtrkTheta(track.fRICHtrkTheta),
- fRICHtrkPhi(track.fRICHtrkPhi),
- fRICHtrkX(track.fRICHtrkX),
- fRICHtrkY(track.fRICHtrkY),
- fRICHmipX(track.fRICHmipX),
- fRICHmipY(track.fRICHmipY),
+ fHMPIDchi2(track.fHMPIDchi2),
+ fHMPIDqn(track.fHMPIDqn),
+ fHMPIDcluIdx(track.fHMPIDcluIdx),
+ fHMPIDsignal(track.fHMPIDsignal),
+ fHMPIDtrkTheta(track.fHMPIDtrkTheta),
+ fHMPIDtrkPhi(track.fHMPIDtrkPhi),
+ fHMPIDtrkX(track.fHMPIDtrkX),
+ fHMPIDtrkY(track.fHMPIDtrkY),
+ fHMPIDmipX(track.fHMPIDmipX),
+ fHMPIDmipY(track.fHMPIDmipY),
fEMCALindex(track.fEMCALindex),
- fPoints(0),
fFriendTrack(0)
{
//
for (Int_t i=0;i<AliPID::kSPECIES;i++) fTOFr[i]=track.fTOFr[i];
for (Int_t i=0;i<3;i++) fTOFLabel[i]=track.fTOFLabel[i];
for (Int_t i=0;i<10;i++) fTOFInfo[i]=track.fTOFInfo[i];
- for (Int_t i=0;i<AliPID::kSPECIES;i++) fRICHr[i]=track.fRICHr[i];
+ for (Int_t i=0;i<AliPID::kSPECIES;i++) fHMPIDr[i]=track.fHMPIDr[i];
if (track.fCp) fCp=new AliExternalTrackParam(*track.fCp);
if (track.fIp) fIp=new AliExternalTrackParam(*track.fIp);
if (track.fOp) fOp=new AliExternalTrackParam(*track.fOp);
- if (track.fPoints) fPoints=new AliTrackPointArray(*(track.fPoints));
if (track.fFriendTrack) fFriendTrack=new AliESDfriendTrack(*(track.fFriendTrack));
}
+//_______________________________________________________________________
+AliESDtrack::AliESDtrack(TParticle * part) :
+ AliExternalTrackParam(),
+ fFlags(0),
+ fLabel(0),
+ fID(0),
+ fTrackLength(0),
+ fD(0),fZ(0),
+ fCdd(0),fCdz(0),fCzz(0),
+ fStopVertex(0),
+ fCp(0),
+ fCchi2(1e10),
+ fIp(0),
+ fOp(0),
+ fITSchi2(0),
+ fITSncls(0),
+ fITSClusterMap(0),
+ fITSsignal(0),
+ fITSLabel(0),
+ fTPCchi2(0),
+ fTPCncls(0),
+ fTPCnclsF(0),
+ fTPCClusterMap(159),//number of padrows
+ fTPCsignal(0),
+ fTPCsignalN(0),
+ fTPCsignalS(0),
+ fTPCLabel(0),
+ fTRDchi2(0),
+ fTRDncls(0),
+ fTRDncls0(0),
+ fTRDsignal(0),
+ fTRDLabel(0),
+ fTRDQuality(0),
+ fTRDBudget(0),
+ fTOFchi2(0),
+ fTOFindex(0),
+ fTOFCalChannel(-1),
+ fTOFsignal(-1),
+ fTOFsignalToT(0),
+ fHMPIDchi2(1e10),
+ fHMPIDqn(-1),
+ fHMPIDcluIdx(-1),
+ fHMPIDsignal(-1),
+ fHMPIDtrkTheta(-1),
+ fHMPIDtrkPhi(-1),
+ fHMPIDtrkX(-1),
+ fHMPIDtrkY(-1),
+ fHMPIDmipX(-1),
+ fHMPIDmipY(-1),
+ fEMCALindex(kEMCALNoMatch),
+ fFriendTrack(0)
+{
+ //
+ // ESD track from TParticle
+ //
+
+ // Reset all the arrays
+ Int_t i, j;
+ for (i=0; i<AliPID::kSPECIES; i++) {
+ fTrackTime[i]=0.;
+ fR[i]=0.;
+ fITSr[i]=0.;
+ fTPCr[i]=0.;
+ fTRDr[i]=0.;
+ fTOFr[i]=0.;
+ fHMPIDr[i]=0.;
+ }
+
+ for (i=0; i<3; i++) { fKinkIndexes[i]=0;}
+ for (i=0; i<3; i++) { fV0Indexes[i]=-1;}
+ for (i=0;i<kNPlane;i++) {
+ for (j=0;j<kNSlice;j++) {
+ fTRDsignals[i][j]=0.;
+ }
+ fTRDTimBin[i]=-1;
+ }
+ for (i=0;i<4;i++) {fTPCPoints[i]=-1;}
+ for (i=0;i<3;i++) {fTOFLabel[i]=-1;}
+ for (i=0;i<10;i++) {fTOFInfo[i]=-1;}
+
+ // Calculate the AliExternalTrackParam content
+
+ Double_t xref;
+ Double_t alpha;
+ Double_t param[5];
+ Double_t covar[15];
+
+ // Calculate alpha: the rotation angle of the corresponding local system (TPC sector)
+ alpha = part->Phi()*180./TMath::Pi();
+ if (alpha<0) alpha+= 360.;
+ if (alpha>360) alpha -= 360.;
+
+ Int_t sector = (Int_t)(alpha/20.);
+ alpha = 10. + 20.*sector;
+ alpha /= 180;
+ alpha *= TMath::Pi();
+
+ // Covariance matrix: no errors, the parameters are exact
+ for (i=0; i<15; i++) covar[i]=0.;
+
+ // Get the vertex of origin and the momentum
+ TVector3 ver(part->Vx(),part->Vy(),part->Vz());
+ TVector3 mom(part->Px(),part->Py(),part->Pz());
+
+ // Rotate to the local coordinate system (TPC sector)
+ ver.RotateZ(-alpha);
+ mom.RotateZ(-alpha);
+
+ // X of the referense plane
+ xref = ver.X();
+
+ Int_t pdgCode = part->GetPdgCode();
+
+ Double_t charge =
+ TDatabasePDG::Instance()->GetParticle(pdgCode)->Charge();
+
+ param[0] = ver.Y();
+ param[1] = ver.Z();
+ param[2] = TMath::Sin(mom.Phi());
+ param[3] = mom.Pz()/mom.Pt();
+ param[4] = TMath::Sign(1/mom.Pt(),charge);
+
+ // Set AliExternalTrackParam
+ Set(xref, alpha, param, covar);
+
+ // Set the PID
+ Int_t indexPID = 99;
+
+ switch (TMath::Abs(pdgCode)) {
+
+ case 11: // electron
+ indexPID = 0;
+ break;
+
+ case 13: // muon
+ indexPID = 1;
+ break;
+
+ case 211: // pion
+ indexPID = 2;
+ break;
+
+ case 321: // kaon
+ indexPID = 3;
+ break;
+
+ case 2212: // proton
+ indexPID = 4;
+ break;
+
+ default:
+ break;
+ }
+
+ // If the particle is not e,mu,pi,K or p the PID probabilities are set to 0
+ if (indexPID < AliPID::kSPECIES) {
+ fR[indexPID]=1.;
+ fITSr[indexPID]=1.;
+ fTPCr[indexPID]=1.;
+ fTRDr[indexPID]=1.;
+ fTOFr[indexPID]=1.;
+ fHMPIDr[indexPID]=1.;
+
+ }
+ // AliESD track label
+ SetLabel(part->GetUniqueID());
+
+}
+
//_______________________________________________________________________
AliESDtrack::~AliESDtrack(){
//
delete fOp;
delete fCp;
delete fFriendTrack;
- delete fPoints;
}
void AliESDtrack::AddCalibObject(TObject * object){
// Reset ITS track related information
fITSchi2 = 0;
fITSncls = 0;
+ fITSClusterMap=0;
fITSsignal = 0;
for (Int_t i=0;i<AliPID::kSPECIES;i++) fITSr[i]=0;
fITSLabel = 0;
for (Int_t i=0;i<3;i++) fTOFLabel[i] = 0;
for (Int_t i=0;i<10;i++) fTOFInfo[i] = 0;
- // Reset RICH related track information
- fRICHchi2 = 0;
- fRICHqn = -1;
- fRICHcluIdx = -1;
- fRICHsignal = 0;
- for (Int_t i=0;i<AliPID::kSPECIES;i++) fRICHr[i] = 0;
- fRICHtrkTheta = -1;
- fRICHtrkPhi = -1;
- fRICHtrkX = -1;
- fRICHtrkY = -1;
- fRICHmipX = -1;
- fRICHmipY = -1;
+ // Reset HMPID related track information
+ fHMPIDchi2 = 0;
+ fHMPIDqn = -1;
+ fHMPIDcluIdx = -1;
+ fHMPIDsignal = 0;
+ for (Int_t i=0;i<AliPID::kSPECIES;i++) fHMPIDr[i] = 0;
+ fHMPIDtrkTheta = -1;
+ fHMPIDtrkPhi = -1;
+ fHMPIDtrkX = -1;
+ fHMPIDtrkY = -1;
+ fHMPIDmipX = -1;
+ fHMPIDmipY = -1;
fEMCALindex = kEMCALNoMatch;
delete fFriendTrack; fFriendTrack = 0;
- delete fPoints; fPoints = 0;
}
//_______________________________________________________________________
Double_t AliESDtrack::GetMass() const {
switch (flags) {
case kITSin: case kITSout: case kITSrefit:
- index=fFriendTrack->GetITSindices();
- for (Int_t i=0;i<AliESDfriendTrack::kMaxITScluster;i++)
- index[i]=t->GetClusterIndex(i);
fITSncls=t->GetNumberOfClusters();
+ index=fFriendTrack->GetITSindices();
+ for (Int_t i=0;i<AliESDfriendTrack::kMaxITScluster;i++) {
+ index[i]=t->GetClusterIndex(i);
+ if (i<fITSncls) {
+ Int_t l=(index[i] & 0xf0000000) >> 28;
+ SETBIT(fITSClusterMap,l);
+ }
+ }
fITSchi2=t->GetChi2();
fITSsignal=t->GetPIDsignal();
fITSLabel = t->GetLabel();
//_______________________________________________________________________
-void AliESDtrack::SetRICHpid(const Double_t *p) {
- // Sets the probability of each particle type (in RICH)
- SetPIDValues(fRICHr,p,AliPID::kSPECIES);
- SetStatus(AliESDtrack::kRICHpid);
+void AliESDtrack::SetHMPIDpid(const Double_t *p) {
+ // Sets the probability of each particle type (in HMPID)
+ SetPIDValues(fHMPIDr,p,AliPID::kSPECIES);
+ SetStatus(AliESDtrack::kHMPIDpid);
}
//_______________________________________________________________________
-void AliESDtrack::GetRICHpid(Double_t *p) const {
- // Gets probabilities of each particle type (in RICH)
- for (Int_t i=0; i<AliPID::kSPECIES; i++) p[i]=fRICHr[i];
+void AliESDtrack::GetHMPIDpid(Double_t *p) const {
+ // Gets probabilities of each particle type (in HMPID)
+ for (Int_t i=0; i<AliPID::kSPECIES; i++) p[i]=fHMPIDr[i];
}
return kTRUE;
}
-void AliESDtrack::SetTrackPointArray(AliTrackPointArray *points) {
- fPoints=points;
- //fFriendTrack->SetTrackPointArray(points);
-}
-const AliTrackPointArray *AliESDtrack::GetTrackPointArray() const {
- return fPoints;
- //return fFriendTrack->GetTrackPointArray();
-}
-
//_______________________________________________________________________
void AliESDtrack::Print(Option_t *) const {
// Prints info on the track
printf("%f, ", p[index]) ;
printf("\n signal = %f\n", GetTOFsignal()) ;
}
- if( IsOn(kRICHpid) ){
- printf("From RICH: ") ;
- GetRICHpid(p) ;
+ if( IsOn(kHMPIDpid) ){
+ printf("From HMPID: ") ;
+ GetHMPIDpid(p) ;
for(index = 0 ; index < AliPID::kSPECIES; index++)
printf("%f, ", p[index]) ;
- printf("\n signal = %f\n", GetRICHsignal()) ;
+ printf("\n signal = %f\n", GetHMPIDsignal()) ;
}
}