From: morsch Date: Thu, 24 Feb 2011 15:33:09 +0000 (+0000) Subject: Number of findable cluster needed in AODTrack X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=3c01c166b54962cd5cd653a6bbb143766542c8d8 Number of findable cluster needed in AODTrack Chiara Bianchin --- diff --git a/ANALYSIS/AliAnalysisTaskESDfilter.cxx b/ANALYSIS/AliAnalysisTaskESDfilter.cxx index 0601697abc0..8881aaac7be 100644 --- a/ANALYSIS/AliAnalysisTaskESDfilter.cxx +++ b/ANALYSIS/AliAnalysisTaskESDfilter.cxx @@ -449,6 +449,7 @@ void AliAnalysisTaskESDfilter::ConvertCascades(const AliESDEvent& esd) aodTrack->SetTPCClusterMap(esdCascadeBach->GetTPCClusterMap()); aodTrack->SetTPCSharedMap (esdCascadeBach->GetTPCSharedMap()); aodTrack->SetChi2perNDF(Chi2perNDF(esdCascadeBach)); + aodTrack->SetTPCPointsF(esdCascadeBach->GetTPCNclsF()); fAODTrackRefs->AddAt(aodTrack,idxBachFromCascade); if (esdCascadeBach->GetSign() > 0) ++fNumberOfPositiveTracks; @@ -532,6 +533,7 @@ void AliAnalysisTaskESDfilter::ConvertCascades(const AliESDEvent& esd) aodTrack->SetTPCClusterMap(esdCascadePos->GetTPCClusterMap()); aodTrack->SetTPCSharedMap (esdCascadePos->GetTPCSharedMap()); aodTrack->SetChi2perNDF(Chi2perNDF(esdCascadePos)); + aodTrack->SetTPCPointsF(esdCascadePos->GetTPCNclsF()); fAODTrackRefs->AddAt(aodTrack,idxPosFromV0Dghter); if (esdCascadePos->GetSign() > 0) ++fNumberOfPositiveTracks; @@ -577,6 +579,7 @@ void AliAnalysisTaskESDfilter::ConvertCascades(const AliESDEvent& esd) aodTrack->SetTPCClusterMap(esdCascadeNeg->GetTPCClusterMap()); aodTrack->SetTPCSharedMap (esdCascadeNeg->GetTPCSharedMap()); aodTrack->SetChi2perNDF(Chi2perNDF(esdCascadeNeg)); + aodTrack->SetTPCPointsF(esdCascadeNeg->GetTPCNclsF()); fAODTrackRefs->AddAt(aodTrack,idxNegFromV0Dghter); if (esdCascadeNeg->GetSign() > 0) ++fNumberOfPositiveTracks; @@ -818,6 +821,7 @@ void AliAnalysisTaskESDfilter::ConvertV0s(const AliESDEvent& esd) aodTrack->SetTPCClusterMap(esdV0Pos->GetTPCClusterMap()); aodTrack->SetTPCSharedMap (esdV0Pos->GetTPCSharedMap()); aodTrack->SetChi2perNDF(Chi2perNDF(esdV0Pos)); + aodTrack->SetTPCPointsF(esdV0Pos->GetTPCNclsF()); fAODTrackRefs->AddAt(aodTrack,posFromV0); // if (fDebug > 0) printf("-------------------Bo: pos track from original pt %.3f \n",aodTrack->Pt()); if (esdV0Pos->GetSign() > 0) ++fNumberOfPositiveTracks; @@ -861,6 +865,7 @@ void AliAnalysisTaskESDfilter::ConvertV0s(const AliESDEvent& esd) aodTrack->SetTPCClusterMap(esdV0Neg->GetTPCClusterMap()); aodTrack->SetTPCSharedMap (esdV0Neg->GetTPCSharedMap()); aodTrack->SetChi2perNDF(Chi2perNDF(esdV0Neg)); + aodTrack->SetTPCPointsF(esdV0Neg->GetTPCNclsF()); fAODTrackRefs->AddAt(aodTrack,negFromV0); // if (fDebug > 0) printf("-------------------Bo: neg track from original pt %.3f \n",aodTrack->Pt()); @@ -993,7 +998,8 @@ void AliAnalysisTaskESDfilter::ConvertTPCOnlyTracks(const AliESDEvent& esd) aodTrack->SetTPCSharedMap (track->GetTPCSharedMap()); aodTrack->SetChi2perNDF(Chi2perNDF(track)); aodTrack->SetFlags(track->GetStatus()); - + aodTrack->SetTPCPointsF(track->GetTPCNclsF()); + delete track; } // end of loop on tracks @@ -1055,7 +1061,8 @@ void AliAnalysisTaskESDfilter::ConvertTracks(const AliESDEvent& esd) aodTrack->SetTPCClusterMap(esdTrack->GetTPCClusterMap()); aodTrack->SetTPCSharedMap (esdTrack->GetTPCSharedMap()); aodTrack->SetChi2perNDF(Chi2perNDF(esdTrack)); - + aodTrack->SetTPCPointsF(esdTrack->GetTPCNclsF()); + fAODTrackRefs->AddAt(aodTrack, nTrack); @@ -1307,7 +1314,8 @@ void AliAnalysisTaskESDfilter::ConvertKinks(const AliESDEvent& esd) mother->SetTPCClusterMap(esdTrackM->GetTPCClusterMap()); mother->SetTPCSharedMap (esdTrackM->GetTPCSharedMap()); mother->SetChi2perNDF(Chi2perNDF(esdTrackM)); - + mother->SetTPCPointsF(esdTrackM->GetTPCNclsF()); + fAODTrackRefs->AddAt(mother, imother); if (esdTrackM->GetSign() > 0) ++fNumberOfPositiveTracks; @@ -1369,6 +1377,7 @@ void AliAnalysisTaskESDfilter::ConvertKinks(const AliESDEvent& esd) selectInfo); daughter->SetTPCClusterMap(esdTrackD->GetTPCClusterMap()); daughter->SetTPCSharedMap (esdTrackD->GetTPCSharedMap()); + daughter->SetTPCPointsF(esdTrackD->GetTPCNclsF()); fAODTrackRefs->AddAt(daughter, idaughter); if (esdTrackD->GetSign() > 0) ++fNumberOfPositiveTracks; diff --git a/STEER/AliAODTrack.cxx b/STEER/AliAODTrack.cxx index 176e0527f75..50ab8a463db 100644 --- a/STEER/AliAODTrack.cxx +++ b/STEER/AliAODTrack.cxx @@ -40,6 +40,7 @@ AliAODTrack::AliAODTrack() : fFilterMap(0), fTPCClusterMap(), fTPCSharedMap(), + fTPCnclsF(0), fID(-999), fCharge(-99), fType(kUndef), @@ -83,6 +84,7 @@ AliAODTrack::AliAODTrack(Short_t id, fFilterMap(selectInfo), fTPCClusterMap(), fTPCSharedMap(), + fTPCnclsF(0), fID(id), fCharge(charge), fType(ttype), @@ -130,6 +132,7 @@ AliAODTrack::AliAODTrack(Short_t id, fFilterMap(selectInfo), fTPCClusterMap(), fTPCSharedMap(), + fTPCnclsF(0), fID(id), fCharge(charge), fType(ttype), @@ -171,6 +174,7 @@ AliAODTrack::AliAODTrack(const AliAODTrack& trk) : fFilterMap(trk.fFilterMap), fTPCClusterMap(trk.fTPCClusterMap), fTPCSharedMap(trk.fTPCSharedMap), + fTPCnclsF(trk.fTPCnclsF), fID(trk.fID), fCharge(trk.fCharge), fType(trk.fType), diff --git a/STEER/AliAODTrack.h b/STEER/AliAODTrack.h index b9f483e8e34..0363fa40a01 100644 --- a/STEER/AliAODTrack.h +++ b/STEER/AliAODTrack.h @@ -213,7 +213,10 @@ class AliAODTrack : public AliVTrack { const TBits& GetTPCSharedMap() const {return fTPCSharedMap;} void SetTPCClusterMap(const TBits amap) {fTPCClusterMap = amap;} void SetTPCSharedMap(const TBits amap) {fTPCSharedMap = amap;} - + void SetTPCPointsF(UShort_t findable){fTPCnclsF = findable;} + + UShort_t GetTPCNclsF() const { return fTPCnclsF;} + AliAODPid *GetDetPid() const { return fDetPid; } AliAODVertex *GetProdVertex() const { return (AliAODVertex*)fProdVertex.GetObject(); } @@ -303,6 +306,7 @@ class AliAODTrack : public AliVTrack { TBits fTPCClusterMap; // Map of clusters, one bit per padrow; 1 if has a cluster on given padrow TBits fTPCSharedMap; // Map of clusters, one bit per padrow; 1 if has a shared cluster on given padrow + UShort_t fTPCnclsF; // findable clusters Short_t fID; // unique track ID, points back to the ESD track @@ -313,7 +317,7 @@ class AliAODTrack : public AliVTrack { AliAODPid *fDetPid; // more detailed or detector specific pid information TRef fProdVertex; // vertex of origin - ClassDef(AliAODTrack, 11); + ClassDef(AliAODTrack, 12); }; inline Bool_t AliAODTrack::IsPrimaryCandidate() const