///////////////////////////////////////////////////////////////////////////////
#include "AliESDMuonTrack.h"
+#include "AliESDMuonCluster.h"
+#include <TClonesArray.h>
#include <TLorentzVector.h>
#include <TMath.h>
//_____________________________________________________________________________
AliESDMuonTrack::AliESDMuonTrack ():
- TObject(),
- fInverseBendingMomentum(0),
+ AliVParticle(),
+ fInverseBendingMomentum(FLT_MAX),
fThetaX(0),
fThetaY(0),
fZ(0),
fBendingCoor(0),
fNonBendingCoor(0),
- fInverseBendingMomentumUncorrected(0),
+ fInverseBendingMomentumAtDCA(FLT_MAX),
+ fThetaXAtDCA(0),
+ fThetaYAtDCA(0),
+ fBendingCoorAtDCA(0),
+ fNonBendingCoorAtDCA(0),
+ fInverseBendingMomentumUncorrected(FLT_MAX),
fThetaXUncorrected(0),
fThetaYUncorrected(0),
fZUncorrected(0),
fBendingCoorUncorrected(0),
fNonBendingCoorUncorrected(0),
fChi2(0),
- fNHit(0),
- fLocalTrigger(234),
fChi2MatchTrigger(0),
- fHitsPatternInTrigCh(0)
+ fLocalTrigger(0),
+ fX1Pattern(0),
+ fY1Pattern(0),
+ fX2Pattern(0),
+ fY2Pattern(0),
+ fX3Pattern(0),
+ fY3Pattern(0),
+ fX4Pattern(0),
+ fY4Pattern(0),
+ fMuonClusterMap(0),
+ fHitsPatternInTrigCh(0),
+ fNHit(0),
+ fClusters(0x0)
{
+ //
// Default constructor
+ //
+ for (Int_t i = 0; i < 15; i++) fCovariances[i] = 0;
}
//_____________________________________________________________________________
-AliESDMuonTrack::AliESDMuonTrack (const AliESDMuonTrack& MUONTrack):
- TObject(MUONTrack),
- fInverseBendingMomentum(MUONTrack.fInverseBendingMomentum),
- fThetaX(MUONTrack.fThetaX),
- fThetaY(MUONTrack.fThetaY),
- fZ(MUONTrack.fZ),
- fBendingCoor(MUONTrack.fBendingCoor),
- fNonBendingCoor(MUONTrack.fNonBendingCoor),
- fInverseBendingMomentumUncorrected(MUONTrack.fInverseBendingMomentumUncorrected),
- fThetaXUncorrected(MUONTrack.fThetaXUncorrected),
- fThetaYUncorrected(MUONTrack.fThetaYUncorrected),
- fZUncorrected(MUONTrack.fZUncorrected),
- fBendingCoorUncorrected(MUONTrack.fBendingCoorUncorrected),
- fNonBendingCoorUncorrected(MUONTrack.fNonBendingCoorUncorrected),
- fChi2(MUONTrack.fChi2),
- fNHit(MUONTrack.fNHit),
- fLocalTrigger(MUONTrack.fLocalTrigger),
- fChi2MatchTrigger(MUONTrack.fChi2MatchTrigger),
- fHitsPatternInTrigCh(MUONTrack.fHitsPatternInTrigCh)
+AliESDMuonTrack::AliESDMuonTrack (const AliESDMuonTrack& muonTrack):
+ AliVParticle(muonTrack),
+ fInverseBendingMomentum(muonTrack.fInverseBendingMomentum),
+ fThetaX(muonTrack.fThetaX),
+ fThetaY(muonTrack.fThetaY),
+ fZ(muonTrack.fZ),
+ fBendingCoor(muonTrack.fBendingCoor),
+ fNonBendingCoor(muonTrack.fNonBendingCoor),
+ fInverseBendingMomentumAtDCA(muonTrack.fInverseBendingMomentumAtDCA),
+ fThetaXAtDCA(muonTrack.fThetaXAtDCA),
+ fThetaYAtDCA(muonTrack.fThetaYAtDCA),
+ fBendingCoorAtDCA(muonTrack.fBendingCoorAtDCA),
+ fNonBendingCoorAtDCA(muonTrack.fNonBendingCoorAtDCA),
+ fInverseBendingMomentumUncorrected(muonTrack.fInverseBendingMomentumUncorrected),
+ fThetaXUncorrected(muonTrack.fThetaXUncorrected),
+ fThetaYUncorrected(muonTrack.fThetaYUncorrected),
+ fZUncorrected(muonTrack.fZUncorrected),
+ fBendingCoorUncorrected(muonTrack.fBendingCoorUncorrected),
+ fNonBendingCoorUncorrected(muonTrack.fNonBendingCoorUncorrected),
+ fChi2(muonTrack.fChi2),
+ fChi2MatchTrigger(muonTrack.fChi2MatchTrigger),
+ fLocalTrigger(muonTrack.fLocalTrigger),
+ fX1Pattern(muonTrack.fX1Pattern),
+ fY1Pattern(muonTrack.fY1Pattern),
+ fX2Pattern(muonTrack.fX2Pattern),
+ fY2Pattern(muonTrack.fY2Pattern),
+ fX3Pattern(muonTrack.fX3Pattern),
+ fY3Pattern(muonTrack.fY3Pattern),
+ fX4Pattern(muonTrack.fX4Pattern),
+ fY4Pattern(muonTrack.fY4Pattern),
+ fMuonClusterMap(muonTrack.fMuonClusterMap),
+ fHitsPatternInTrigCh(muonTrack.fHitsPatternInTrigCh),
+ fNHit(muonTrack.fNHit),
+ fClusters(0x0)
{
//
// Copy constructor
// Deep copy implemented
//
+ for (Int_t i = 0; i < 15; i++) fCovariances[i] = muonTrack.fCovariances[i];
+
+ // necessary to make a copy of the objects and not only the pointers in TClonesArray
+ if (muonTrack.fClusters) {
+ fClusters = new TClonesArray("AliESDMuonCluster",muonTrack.fClusters->GetEntriesFast());
+ AliESDMuonCluster *cluster = (AliESDMuonCluster*) muonTrack.fClusters->First();
+ while (cluster) {
+ new ((*fClusters)[fClusters->GetEntriesFast()]) AliESDMuonCluster(*cluster);
+ cluster = (AliESDMuonCluster*) muonTrack.fClusters->After(cluster);
+ }
+ }
}
//_____________________________________________________________________________
-AliESDMuonTrack& AliESDMuonTrack::operator=(const AliESDMuonTrack& MUONTrack)
+AliESDMuonTrack& AliESDMuonTrack::operator=(const AliESDMuonTrack& muonTrack)
{
//
// Equal operator for a deep copy
//
- if (this == &MUONTrack)
+ if (this == &muonTrack)
return *this;
- fInverseBendingMomentum = MUONTrack.fInverseBendingMomentum;
- fThetaX = MUONTrack.fThetaX;
- fThetaY = MUONTrack.fThetaY;
- fZ = MUONTrack.fZ;
- fBendingCoor = MUONTrack.fBendingCoor;
- fNonBendingCoor = MUONTrack.fNonBendingCoor;
+ AliVParticle::operator=(muonTrack); // don't forget to invoke the base class' assignment operator
+
+ fInverseBendingMomentum = muonTrack.fInverseBendingMomentum;
+ fThetaX = muonTrack.fThetaX;
+ fThetaY = muonTrack.fThetaY;
+ fZ = muonTrack.fZ;
+ fBendingCoor = muonTrack.fBendingCoor;
+ fNonBendingCoor = muonTrack.fNonBendingCoor;
- fInverseBendingMomentumUncorrected = MUONTrack.fInverseBendingMomentumUncorrected;
- fThetaXUncorrected = MUONTrack.fThetaXUncorrected;
- fThetaYUncorrected = MUONTrack.fThetaYUncorrected;
- fZUncorrected = MUONTrack.fZUncorrected;
- fBendingCoorUncorrected = MUONTrack.fBendingCoorUncorrected;
- fNonBendingCoorUncorrected = MUONTrack.fNonBendingCoorUncorrected;
+ fInverseBendingMomentumAtDCA = muonTrack.fInverseBendingMomentumAtDCA;
+ fThetaXAtDCA = muonTrack.fThetaXAtDCA;
+ fThetaYAtDCA = muonTrack.fThetaYAtDCA;
+ fBendingCoorAtDCA = muonTrack.fBendingCoorAtDCA;
+ fNonBendingCoorAtDCA = muonTrack.fNonBendingCoorAtDCA;
- fChi2 = MUONTrack.fChi2;
- fNHit = MUONTrack.fNHit;
+ fInverseBendingMomentumUncorrected = muonTrack.fInverseBendingMomentumUncorrected;
+ fThetaXUncorrected = muonTrack.fThetaXUncorrected;
+ fThetaYUncorrected = muonTrack.fThetaYUncorrected;
+ fZUncorrected = muonTrack.fZUncorrected;
+ fBendingCoorUncorrected = muonTrack.fBendingCoorUncorrected;
+ fNonBendingCoorUncorrected = muonTrack.fNonBendingCoorUncorrected;
+
+ for (Int_t i = 0; i < 15; i++) fCovariances[i] = muonTrack.fCovariances[i];
+
+ fChi2 = muonTrack.fChi2;
+ fNHit = muonTrack.fNHit;
- fLocalTrigger = MUONTrack.fLocalTrigger;
- fChi2MatchTrigger = MUONTrack.fChi2MatchTrigger;
+ fLocalTrigger = muonTrack.fLocalTrigger;
+ fX1Pattern = muonTrack.fX1Pattern;
+ fY1Pattern = muonTrack.fY1Pattern;
+ fX2Pattern = muonTrack.fX2Pattern;
+ fY2Pattern = muonTrack.fY2Pattern;
+ fX3Pattern = muonTrack.fX3Pattern;
+ fY3Pattern = muonTrack.fY3Pattern;
+ fX4Pattern = muonTrack.fX4Pattern;
+ fY4Pattern = muonTrack.fY4Pattern;
+ fChi2MatchTrigger = muonTrack.fChi2MatchTrigger;
- fHitsPatternInTrigCh = MUONTrack.fHitsPatternInTrigCh;
+ fHitsPatternInTrigCh = muonTrack.fHitsPatternInTrigCh;
+ fMuonClusterMap = muonTrack.fMuonClusterMap;
+
+ // necessary to make a copy of the objects and not only the pointers in TClonesArray
+ delete fClusters;
+ if (muonTrack.fClusters) {
+ fClusters = new TClonesArray("AliESDMuonCluster",muonTrack.fClusters->GetEntriesFast());
+ AliESDMuonCluster *cluster = (AliESDMuonCluster*) muonTrack.fClusters->First();
+ while (cluster) {
+ new ((*fClusters)[fClusters->GetEntriesFast()]) AliESDMuonCluster(*cluster);
+ cluster = (AliESDMuonCluster*) muonTrack.fClusters->After(cluster);
+ }
+ } else fClusters = 0x0;
+
return *this;
}
+void AliESDMuonTrack::Copy(TObject &obj) const {
+
+ // this overwrites the virtual TOBject::Copy()
+ // to allow run time copying without casting
+ // in AliESDEvent
+
+ if(this==&obj)return;
+ AliESDMuonTrack *robj = dynamic_cast<AliESDMuonTrack*>(&obj);
+ if(!robj)return; // not an AliESDMuonTrack
+ *robj = *this;
+
+}
+
+
+//__________________________________________________________________________
+AliESDMuonTrack::~AliESDMuonTrack()
+{
+ /// Destructor
+ delete fClusters;
+}
+
+//__________________________________________________________________________
+void AliESDMuonTrack::Clear(Option_t* opt)
+{
+ /// Clear arrays
+ if (fClusters) fClusters->Clear(opt);
+}
+
+//__________________________________________________________________________
+void AliESDMuonTrack::Reset()
+{
+ /// Reset to default values
+ SetUniqueID(0);
+ fInverseBendingMomentum = FLT_MAX;
+ fThetaX = 0.;
+ fThetaY = 0.;
+ fZ = 0.;
+ fBendingCoor = 0.;
+ fNonBendingCoor = 0.;
+ fInverseBendingMomentumAtDCA = FLT_MAX;
+ fThetaXAtDCA = 0.;
+ fThetaYAtDCA = 0.;
+ fBendingCoorAtDCA = 0.;
+ fNonBendingCoorAtDCA = 0.;
+ fInverseBendingMomentumUncorrected = FLT_MAX;
+ fThetaXUncorrected = 0.;
+ fThetaYUncorrected = 0.;
+ fZUncorrected = 0.;
+ fBendingCoorUncorrected = 0.;
+ fNonBendingCoorUncorrected = 0.;
+ fChi2 = 0.;
+ fChi2MatchTrigger = 0.;
+ fLocalTrigger = 0;
+ fX1Pattern = 0;
+ fY1Pattern = 0;
+ fX2Pattern = 0;
+ fY2Pattern = 0;
+ fX3Pattern = 0;
+ fY3Pattern = 0;
+ fX4Pattern = 0;
+ fY4Pattern = 0;
+ fMuonClusterMap = 0;
+ fHitsPatternInTrigCh = 0;
+ fNHit = 0;
+ delete fClusters; fClusters = 0x0;
+ for (Int_t i = 0; i < 15; i++) fCovariances[i] = 0.;
+}
+
+//_____________________________________________________________________________
+void AliESDMuonTrack::GetCovariances(TMatrixD& cov) const
+{
+ // return covariance matrix of uncorrected parameters
+ cov.ResizeTo(5,5);
+ for (Int_t i = 0; i < 5; i++)
+ for (Int_t j = 0; j <= i; j++)
+ cov(i,j) = cov (j,i) = fCovariances[i*(i+1)/2 + j];
+}
+
+//_____________________________________________________________________________
+void AliESDMuonTrack::SetCovariances(const TMatrixD& cov)
+{
+ // set reduced covariance matrix of uncorrected parameters
+ for (Int_t i = 0; i < 5; i++)
+ for (Int_t j = 0; j <= i; j++)
+ fCovariances[i*(i+1)/2 + j] = cov(i,j);
+
+}
+
+//_____________________________________________________________________________
+void AliESDMuonTrack::GetCovarianceXYZPxPyPz(Double_t cov[21]) const
+{
+ // return reduced covariance matrix of uncorrected parameters in (X,Y,Z,Px,Py,Pz) coordinate system
+ //
+ // Cov(x,x) ... : cov[0]
+ // Cov(y,x) ... : cov[1] cov[2]
+ // Cov(z,x) ... : cov[3] cov[4] cov[5]
+ // Cov(px,x)... : cov[6] cov[7] cov[8] cov[9]
+ // Cov(py,x)... : cov[10] cov[11] cov[12] cov[13] cov[14]
+ // Cov(pz,x)... : cov[15] cov[16] cov[17] cov[18] cov[19] cov[20]
+ //
+ // Get ESD covariance matrix into a TMatrixD
+ TMatrixD covESD(5,5);
+ GetCovariances(covESD);
+
+ // compute Jacobian to change the coordinate system
+ // from (X,thetaX,Y,thetaY,c/pYZ) to (X,Y,Z,pX,pY,pZ)
+ Double_t tanThetaX = TMath::Tan(fThetaXUncorrected);
+ Double_t tanThetaY = TMath::Tan(fThetaYUncorrected);
+ Double_t cosThetaX2 = TMath::Cos(fThetaXUncorrected) * TMath::Cos(fThetaXUncorrected);
+ Double_t cosThetaY2 = TMath::Cos(fThetaYUncorrected) * TMath::Cos(fThetaYUncorrected);
+ Double_t pZ = PzUncorrected();
+ Double_t dpZdthetaY = - fInverseBendingMomentumUncorrected * fInverseBendingMomentumUncorrected *
+ pZ * pZ * pZ * tanThetaY / cosThetaY2;
+ Double_t dpZdinvpYZ = (fInverseBendingMomentumUncorrected != 0.) ? - pZ / fInverseBendingMomentumUncorrected : - FLT_MAX;
+ TMatrixD jacob(6,5);
+ jacob.Zero();
+ jacob(0,0) = 1.;
+ jacob(1,2) = 1.;
+ jacob(3,1) = pZ / cosThetaX2;
+ jacob(3,3) = dpZdthetaY * tanThetaX;
+ jacob(3,4) = dpZdinvpYZ * tanThetaX;
+ jacob(4,3) = dpZdthetaY * tanThetaY + pZ / cosThetaY2;
+ jacob(4,4) = dpZdinvpYZ * tanThetaY;
+ jacob(5,3) = dpZdthetaY;
+ jacob(5,4) = dpZdinvpYZ;
+
+ // compute covariance matrix in AOD coordinate system
+ TMatrixD tmp(covESD,TMatrixD::kMultTranspose,jacob);
+ TMatrixD covAOD(jacob,TMatrixD::kMult,tmp);
+
+ // Get AOD covariance matrix into co[21]
+ for (Int_t i = 0; i < 6; i++)
+ for (Int_t j = 0; j <= i; j++)
+ cov[i*(i+1)/2 + j] = covAOD(i,j);
+
+}
+
//_____________________________________________________________________________
Double_t AliESDMuonTrack::Px() const
{
// return p_x from track parameters
Double_t nonBendingSlope = TMath::Tan(fThetaX);
Double_t bendingSlope = TMath::Tan(fThetaY);
- Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : 0.;
+ Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : - FLT_MAX;
Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
return pZ * nonBendingSlope;
}
{
// return p_y from track parameters
Double_t bendingSlope = TMath::Tan(fThetaY);
- Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : 0.;
+ Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : - FLT_MAX;
Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
return pZ * bendingSlope;
}
{
// return p_z from track parameters
Double_t bendingSlope = TMath::Tan(fThetaY);
- Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : 0.;
+ Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : - FLT_MAX;
return -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
}
// return p from track parameters
Double_t nonBendingSlope = TMath::Tan(fThetaX);
Double_t bendingSlope = TMath::Tan(fThetaY);
- Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : 0.;
+ Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : - FLT_MAX;
Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
return -pZ * TMath::Sqrt(1.0 + bendingSlope*bendingSlope + nonBendingSlope*nonBendingSlope);
}
void AliESDMuonTrack::LorentzP(TLorentzVector& vP) const
{
// return Lorentz momentum vector from track parameters
- Double_t muonMass = 0.105658369;
+ Double_t muonMass = M();
Double_t nonBendingSlope = TMath::Tan(fThetaX);
Double_t bendingSlope = TMath::Tan(fThetaY);
- Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : 0.;
+ Double_t pYZ = (fInverseBendingMomentum != 0.) ? TMath::Abs(1. / fInverseBendingMomentum) : - FLT_MAX;
+ Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
+ Double_t pX = pZ * nonBendingSlope;
+ Double_t pY = pZ * bendingSlope;
+ Double_t e = TMath::Sqrt(muonMass*muonMass + pX*pX + pY*pY + pZ*pZ);
+ vP.SetPxPyPzE(pX, pY, pZ, e);
+}
+
+//_____________________________________________________________________________
+Double_t AliESDMuonTrack::PxAtDCA() const
+{
+ // return p_x from track parameters
+ Double_t nonBendingSlope = TMath::Tan(fThetaXAtDCA);
+ Double_t bendingSlope = TMath::Tan(fThetaYAtDCA);
+ Double_t pYZ = (fInverseBendingMomentumAtDCA != 0.) ? TMath::Abs(1. / fInverseBendingMomentumAtDCA) : - FLT_MAX;
+ Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
+ return pZ * nonBendingSlope;
+}
+
+//_____________________________________________________________________________
+Double_t AliESDMuonTrack::PyAtDCA() const
+{
+ // return p_y from track parameters
+ Double_t bendingSlope = TMath::Tan(fThetaYAtDCA);
+ Double_t pYZ = (fInverseBendingMomentumAtDCA != 0.) ? TMath::Abs(1. / fInverseBendingMomentumAtDCA) : - FLT_MAX;
+ Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
+ return pZ * bendingSlope;
+}
+
+//_____________________________________________________________________________
+Double_t AliESDMuonTrack::PzAtDCA() const
+{
+ // return p_z from track parameters
+ Double_t bendingSlope = TMath::Tan(fThetaYAtDCA);
+ Double_t pYZ = (fInverseBendingMomentumAtDCA != 0.) ? TMath::Abs(1. / fInverseBendingMomentumAtDCA) : - FLT_MAX;
+ return -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
+}
+
+//_____________________________________________________________________________
+Double_t AliESDMuonTrack::PAtDCA() const
+{
+ // return p from track parameters
+ Double_t nonBendingSlope = TMath::Tan(fThetaXAtDCA);
+ Double_t bendingSlope = TMath::Tan(fThetaYAtDCA);
+ Double_t pYZ = (fInverseBendingMomentumAtDCA != 0.) ? TMath::Abs(1. / fInverseBendingMomentumAtDCA) : - FLT_MAX;
+ Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
+ return -pZ * TMath::Sqrt(1.0 + bendingSlope*bendingSlope + nonBendingSlope*nonBendingSlope);
+}
+
+//_____________________________________________________________________________
+void AliESDMuonTrack::LorentzPAtDCA(TLorentzVector& vP) const
+{
+ // return Lorentz momentum vector from track parameters
+ Double_t muonMass = M();
+ Double_t nonBendingSlope = TMath::Tan(fThetaXAtDCA);
+ Double_t bendingSlope = TMath::Tan(fThetaYAtDCA);
+ Double_t pYZ = (fInverseBendingMomentumAtDCA != 0.) ? TMath::Abs(1. / fInverseBendingMomentumAtDCA) : - FLT_MAX;
Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
Double_t pX = pZ * nonBendingSlope;
Double_t pY = pZ * bendingSlope;
// return p_x from track parameters
Double_t nonBendingSlope = TMath::Tan(fThetaXUncorrected);
Double_t bendingSlope = TMath::Tan(fThetaYUncorrected);
- Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : 0.;
+ Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : - FLT_MAX;
Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
return pZ * nonBendingSlope;
}
{
// return p_y from track parameters
Double_t bendingSlope = TMath::Tan(fThetaYUncorrected);
- Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : 0.;
+ Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : - FLT_MAX;
Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
return pZ * bendingSlope;
}
{
// return p_z from track parameters
Double_t bendingSlope = TMath::Tan(fThetaYUncorrected);
- Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : 0.;
+ Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : - FLT_MAX;
return -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
}
// return p from track parameters
Double_t nonBendingSlope = TMath::Tan(fThetaXUncorrected);
Double_t bendingSlope = TMath::Tan(fThetaYUncorrected);
- Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : 0.;
+ Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : - FLT_MAX;
Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
return -pZ * TMath::Sqrt(1.0 + bendingSlope*bendingSlope + nonBendingSlope*nonBendingSlope);
}
void AliESDMuonTrack::LorentzPUncorrected(TLorentzVector& vP) const
{
// return Lorentz momentum vector from track parameters
- Double_t muonMass = 0.105658369;
+ Double_t muonMass = M();
Double_t nonBendingSlope = TMath::Tan(fThetaXUncorrected);
Double_t bendingSlope = TMath::Tan(fThetaYUncorrected);
- Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : 0.;
+ Double_t pYZ = (fInverseBendingMomentumUncorrected != 0.) ? TMath::Abs(1. / fInverseBendingMomentumUncorrected) : - FLT_MAX;
Double_t pZ = -pYZ / TMath::Sqrt(1.0 + bendingSlope*bendingSlope); // spectro. (z<0)
Double_t pX = pZ * nonBendingSlope;
Double_t pY = pZ * bendingSlope;
// 2 track match Low pt cut
// 3 track match High pt cut
- if (LoCircuit() == -1) {
+ if (!LoCircuit()) {
return 0;
} else if (LoLpt() == 0 && LoHpt() == 0) {
return 1;
}
+//_____________________________________________________________________________
+Int_t AliESDMuonTrack::GetNClusters() const
+{
+ // return the number of clusters associated to the track
+ if (!fClusters) return 0;
+
+ return fClusters->GetEntriesFast();
+}
+
+//_____________________________________________________________________________
+TClonesArray& AliESDMuonTrack::GetClusters() const
+{
+ // return the array of clusters associated to the track
+ if (!fClusters) fClusters = new TClonesArray("AliESDMuonCluster",10);
+
+ return *fClusters;
+}
+
+//_____________________________________________________________________________
+void AliESDMuonTrack::AddCluster(const AliESDMuonCluster &cluster)
+{
+ // add a cluster to the TClonesArray of clusters associated to the track
+ if (!fClusters) fClusters = new TClonesArray("AliESDMuonCluster",10);
+
+ new ((*fClusters)[fClusters->GetEntriesFast()]) AliESDMuonCluster(cluster);
+}
+
+//_____________________________________________________________________________
+Bool_t AliESDMuonTrack::ClustersStored() const
+{
+ // return kTRUE if the clusters associated to the track are registered
+ if (GetNClusters() == 0) return kFALSE;
+
+ return kTRUE;
+}
+