From c39c35d132ae2f7966bedd9959dcc5cbd8a34b94 Mon Sep 17 00:00:00 2001 From: morsch Date: Fri, 2 Nov 2012 14:49:42 +0000 Subject: [PATCH] New faeture splitting one event "by track charge". So far we had only "by eta" and "randomly" Robert Grajcarek --- ANALYSIS/AliEPSelectionTask.cxx | 27 +++++++++++++++++++++++++++ ANALYSIS/AliEPSelectionTask.h | 2 +- 2 files changed, 28 insertions(+), 1 deletion(-) diff --git a/ANALYSIS/AliEPSelectionTask.cxx b/ANALYSIS/AliEPSelectionTask.cxx index 9ba565fcffb..99b7d61e68a 100644 --- a/ANALYSIS/AliEPSelectionTask.cxx +++ b/ANALYSIS/AliEPSelectionTask.cxx @@ -531,6 +531,33 @@ void AliEPSelectionTask::GetQsub(TVector2 &Q1, TVector2 &Q2, TObjArray* tracklis } } } + } else if (fSplitMethod == AliEPSelectionTask::kCharge) { + + for (Int_t i = 0; i < nt; i++) { + weight = 1; + track = dynamic_cast (tracklist->At(i)); + if (!track) continue; + weight = GetWeight(track); + Short_t cha = track->Charge(); + idtemp = track->GetID(); + if ((fAnalysisInput.CompareTo("AOD")==0) && (fAODfilterbit == 128)) idtemp = idtemp*(-1) - 1; + + if (cha > 0) { + mQx1 += (weight*cos(2*track->Phi())); + mQy1 += (weight*sin(2*track->Phi())); + if (fSaveTrackContribution){ + EP->GetQContributionXArraysub1()->AddAt(weight*cos(2*track->Phi()),idtemp); + EP->GetQContributionYArraysub1()->AddAt(weight*sin(2*track->Phi()),idtemp); + } + } else if (cha < 0) { + mQx2 += (weight*cos(2*track->Phi())); + mQy2 += (weight*sin(2*track->Phi())); + if (fSaveTrackContribution){ + EP->GetQContributionXArraysub2()->AddAt(weight*cos(2*track->Phi()),idtemp); + EP->GetQContributionYArraysub2()->AddAt(weight*sin(2*track->Phi()),idtemp); + } + } + } } else { printf("plane resolution determination method not available!\n\n "); return; diff --git a/ANALYSIS/AliEPSelectionTask.h b/ANALYSIS/AliEPSelectionTask.h index fc9731bfab7..646443277c1 100644 --- a/ANALYSIS/AliEPSelectionTask.h +++ b/ANALYSIS/AliEPSelectionTask.h @@ -30,7 +30,7 @@ class AliEPSelectionTask : public AliAnalysisTaskSE { public: - enum ResoMethod{kRandom,kEta}; + enum ResoMethod{kRandom,kEta,kCharge}; AliEPSelectionTask(); AliEPSelectionTask(const char *name); -- 2.43.0