ClassImp(AliAnalysisTwoParticleResonanceFlowTask)
AliAnalysisTwoParticleResonanceFlowTask::AliAnalysisTwoParticleResonanceFlowTask() : AliAnalysisTaskSE(),
- fSpeciesA(0), fSpeciesB(0), fChargeA(0), fChargeB(0), fMassA(0), fMassB(0), fDebug(0), fIsMC(0), fEventMixing(0), fPhiMinusPsiMethod(0), fQA(0), fV0(0), fMassBins(1), fMinMass(-1.), fMaxMass(0.), fCutsRP(NULL), fNullCuts(0), fPIDResponse(0), fFlowEvent(0), fBayesianResponse(0), fCandidates(0), fCandidateEtaPtCut(0), fCandidateMinEta(0), fCandidateMaxEta(0), fCandidateMinPt(0), fCandidateMaxPt(0), fNPtBins(18), fNdPhiBins(18), fCentrality(999), fVertex(999), fAOD(0), fPoolManager(0), fOutputList(0), fEventStats(0), fCentralityPass(0), fCentralityNoPass(0), fNOPID(0), fPIDk(0), fPIDp(0), fNOPIDTOF(0), fPIDTOF(0), fPtP(0), fPtN(0), fPtKP(0), fPtKN(0), fCentralityMin(0), fCentralityMax(100), fkCentralityMethod(0), fPOICuts(0), fVertexRange(0), fPhi(0), fPt(0), fEta(0), fVZEROA(0), fVZEROC(0), fTPCM(0), fDeltaDipAngle(0), fDeltaDipPt(0), fApplyDeltaDipCut(0), fDCAAll(0), fDCAXYQA(0), fDCAZQA(0), fDCAPrim(0), fDCASecondaryWeak(0), fDCAMaterial(0), fSubEventDPhiv2(0), fAnalysisSummary(0)
+ fSpeciesA(0), fSpeciesB(0), fChargeA(0), fChargeB(0), fMassA(0), fMassB(0), fMinPtA(0), fMaxPtA(0), fMinPtB(0), fMaxPtB(0), fIsMC(0), fEventMixing(0), fPhiMinusPsiMethod(0), fQA(0), fV0(0), fMassBins(1), fMinMass(-1.), fMaxMass(0.), fCutsRP(NULL), fNullCuts(0), fPIDResponse(0), fFlowEvent(0), fBayesianResponse(0), fCandidates(0), fCandidateEtaPtCut(0), fCandidateMinEta(0), fCandidateMaxEta(0), fCandidateMinPt(0), fCandidateMaxPt(0), fNPtBins(18), fNdPhiBins(18), fCentrality(999), fVertex(999), fAOD(0), fPoolManager(0), fOutputList(0), fEventStats(0), fCentralityPass(0), fCentralityNoPass(0), fNOPID(0), fPIDk(0), fPIDp(0), fPtP(0), fPtN(0), fPtSpeciesA(0), fPtSpeciesB(0), fCentralityMin(0), fCentralityMax(100), fkCentralityMethod(0), fPOICuts(0), fVertexRange(0), fPhi(0), fEta(0), fVZEROA(0), fVZEROC(0), fTPCM(0), fDeltaDipAngle(0), fDeltaDipPt(0), fApplyDeltaDipCut(0), fDCAAll(0), fDCAXYQA(0), fDCAZQA(0), fDCAPrim(0), fDCASecondaryWeak(0), fDCAMaterial(0), fSubEventDPhiv2(0), fAnalysisSummary(0)
{
// Default constructor
for(Int_t i(0); i < 7; i++) fPIDConfig[i] = 0.;
}
}
//_____________________________________________________________________________
-AliAnalysisTwoParticleResonanceFlowTask::AliAnalysisTwoParticleResonanceFlowTask(const char *name) : AliAnalysisTaskSE(name), fSpeciesA(0), fSpeciesB(0), fChargeA(0), fChargeB(0), fMassA(0), fMassB(0), fDebug(0), fIsMC(0), fEventMixing(0), fPhiMinusPsiMethod(0), fQA(0), fV0(0), fMassBins(1), fMinMass(-1.), fMaxMass(0.), fCutsRP(NULL), fNullCuts(0), fPIDResponse(0), fFlowEvent(0), fBayesianResponse(0), fCandidates(0), fCandidateEtaPtCut(0), fCandidateMinEta(0), fCandidateMaxEta(0), fCandidateMinPt(0), fCandidateMaxPt(0), fNPtBins(18), fNdPhiBins(18), fCentrality(999), fVertex(999), fAOD(0), fPoolManager(0), fOutputList(0), fEventStats(0), fCentralityPass(0), fCentralityNoPass(0), fNOPID(0), fPIDk(0), fPIDp(0), fNOPIDTOF(0), fPIDTOF(0), fPtP(0), fPtN(0), fPtKP(0), fPtKN(0), fCentralityMin(0), fCentralityMax(100), fkCentralityMethod(0), fPOICuts(0), fVertexRange(0), fPhi(0), fPt(0), fEta(0), fVZEROA(0), fVZEROC(0), fTPCM(0), fDeltaDipAngle(0), fDeltaDipPt(0), fApplyDeltaDipCut(0), fDCAAll(0), fDCAXYQA(0), fDCAZQA(0), fDCAPrim(0), fDCASecondaryWeak(0), fDCAMaterial(0), fSubEventDPhiv2(0), fAnalysisSummary(0)
+AliAnalysisTwoParticleResonanceFlowTask::AliAnalysisTwoParticleResonanceFlowTask(const char *name) : AliAnalysisTaskSE(name), fSpeciesA(0), fSpeciesB(0), fChargeA(0), fChargeB(0), fMassA(0), fMassB(0), fMinPtA(0), fMaxPtA(0), fMinPtB(0), fMaxPtB(0), fIsMC(0), fEventMixing(0), fPhiMinusPsiMethod(0), fQA(0), fV0(0), fMassBins(1), fMinMass(-1.), fMaxMass(0.), fCutsRP(NULL), fNullCuts(0), fPIDResponse(0), fFlowEvent(0), fBayesianResponse(0), fCandidates(0), fCandidateEtaPtCut(0), fCandidateMinEta(0), fCandidateMaxEta(0), fCandidateMinPt(0), fCandidateMaxPt(0), fNPtBins(18), fNdPhiBins(18), fCentrality(999), fVertex(999), fAOD(0), fPoolManager(0), fOutputList(0), fEventStats(0), fCentralityPass(0), fCentralityNoPass(0), fNOPID(0), fPIDk(0), fPIDp(0), fPtP(0), fPtN(0), fPtSpeciesA(0), fPtSpeciesB(0), fCentralityMin(0), fCentralityMax(100), fkCentralityMethod(0), fPOICuts(0), fVertexRange(0), fPhi(0), fEta(0), fVZEROA(0), fVZEROC(0), fTPCM(0), fDeltaDipAngle(0), fDeltaDipPt(0), fApplyDeltaDipCut(0), fDCAAll(0), fDCAXYQA(0), fDCAZQA(0), fDCAPrim(0), fDCASecondaryWeak(0), fDCAMaterial(0), fSubEventDPhiv2(0), fAnalysisSummary(0)
{
// Constructor
for(Int_t i(0); i < 7; i++) fPIDConfig[i] = 0.;
DefineInput(0, TChain::Class());
DefineOutput(1, TList::Class());
DefineOutput(2, AliFlowEventSimple::Class());
- if(fDebug > 0) cout << " === Created instance of AliAnaysisTaskGenericResonanceFlowAnalysis === " << endl;
}
//_____________________________________________________________________________
AliAnalysisTwoParticleResonanceFlowTask::~AliAnalysisTwoParticleResonanceFlowTask()
if (fFlowEvent) delete fFlowEvent;
if (fBayesianResponse) delete fBayesianResponse;
if (fEventMixing) delete fPoolManager;
- if (fDebug > 0) cout << " === Deleted instance of AliAnalysisTwoParticleResonanceFlowTask === " << endl;
}
//_____________________________________________________________________________
void AliAnalysisTwoParticleResonanceFlowTask::ForceExit(Int_t type, const char* message)
TH1F* AliAnalysisTwoParticleResonanceFlowTask::BookHistogram(const char* name)
{
// Return a pointer to a TH1 with predefined binning
- if(fDebug > 0) cout << " *** BookHistogram() *** " << name << endl;
TH1F *hist = new TH1F(name, Form("M_{INV} (%s)", name), 2*fMassBins, fMinMass, fMaxMass);
hist->GetXaxis()->SetTitle("M_{INV} (GeV / c^{2})");
hist->GetYaxis()->SetTitle("No. of pairs");
TH2F* AliAnalysisTwoParticleResonanceFlowTask::BookPIDHistogram(const char* name, Bool_t TPC)
{
// Return a pointer to a TH2 with predefined binning
- if(fDebug > 0) cout << " *** BookPIDHisotgram() *** " << endl;
TH2F *hist = 0x0;
if(TPC) {
hist = new TH2F(name, Form("PID (%s)", name), 100, 0, 5, 100, 0, 1000);
TH1F* AliAnalysisTwoParticleResonanceFlowTask::InitPtSpectraHistograms(Float_t nmin, Float_t nmax)
{
// intialize p_t histograms for each p_t bin
- if(fDebug > 0) cout << " *** InitPtSpectraHistograms() *** " << endl;
TH1F* hist = new TH1F(Form("%4.2f p_{t} %4.2f", nmin, nmax), Form("%f p_{t} %f", nmin, nmax), 2*fMassBins, nmin, nmax);
hist->GetXaxis()->SetTitle("p_{T} GeV / c");
fOutputList->Add(hist);
TH1F* AliAnalysisTwoParticleResonanceFlowTask::BookPtHistogram(const char* name)
{
// Return a pointer to a p_T spectrum histogram
- if(fDebug > 0) cout << " *** BookPtHistogram() *** " << endl;
TH1F* ratio = new TH1F(name, name, 100, 0, 7);
ratio->GetXaxis()->SetTitle("p_{T} ( GeV / c^{2} )");
ratio->GetYaxis()->SetTitle("No. of events");
{
// Create user defined output objects
if(!InitializeAnalysis()) ForceExit(0, " > Couldn't initialize the analysis !!! <");
- if(fDebug > 0) cout << " *** UserCreateOutputObjects() *** " << endl;
// Create all output objects and store them to a list
fOutputList = new TList();
fOutputList->SetOwner(kTRUE);
fCentralityNoPass = new TH1F("fCentralityNoPass", "Centrality No Pass", 101, -1, 100);
fOutputList->Add(fCentralityNoPass);
fNOPID = BookPIDHistogram("TPC signal, all particles", kTRUE);
- fPIDk = BookPIDHistogram("TPC signal, kaons", kTRUE);
- fPIDp = BookPIDHistogram("TPC signal, pions", kTRUE);
- fNOPIDTOF = BookPIDHistogram("TOF signal, all particles", kFALSE);
- fPIDTOF = BookPIDHistogram("TOF signal, kaons", kFALSE);
+ fPIDk = BookPIDHistogram("TPC signal, Species A", kTRUE);
+ fPIDp = BookPIDHistogram("TPC signal, Species B", kTRUE);
}
for(Int_t ptbin(0); ptbin < fNPtBins; ptbin++) {
if(!fPhiMinusPsiMethod) {
for(Int_t ptbin(0); ptbin < fNPtBins; ptbin++) fPtSpectra[ptbin] = InitPtSpectraHistograms(fPtBins[ptbin], fPtBins[ptbin+1]);
fPtP = BookPtHistogram("i^{+}");
fPtN = BookPtHistogram("i^{-}");
- fPtKP = BookPtHistogram("K^{+}");
- fPtKN = BookPtHistogram("K^{-}");
+ fPtSpeciesA = BookPtHistogram("p_{T} spectrum Species A");
+ fPtSpeciesB = BookPtHistogram("p_{T} spectrum Species B");
fPhi = new TH1F("fPhi", "#phi distribution", 100, -.5, 7);
fOutputList->Add(fPhi);
- fPt = new TH1F("fPt", "p_{T}", 100, 0, 5.5);
- fOutputList->Add(fPt);
fEta = new TH1F("fEta", "#eta distribution", 100, -1.1, 1.1);
fOutputList->Add(fEta);
fVZEROA = new TH1F("fVZEROA", "VZERO A Multiplicity", 1000, 0, 10000);
fOutputList->Add(fV0Data[ptbin][iV0]);
}
}
- fAnalysisSummary = new TH1F("Analysis summary","Analysis summary", 26, 0, 26);
+ TString name = "fAnalysisSummary_";
+ name+=this->GetName();
+ fAnalysisSummary = new TH1F(name.Data(), name.Data(), 34, 0, 34);
fAnalysisSummary->GetXaxis()->SetBinLabel(1, "fIsMC");
fAnalysisSummary->GetXaxis()->SetBinLabel(2, "fEventMixing");
fAnalysisSummary->GetXaxis()->SetBinLabel(3, "fAODAnalysis");
- fAnalysisSummary->GetXaxis()->SetBinLabel(4, "bayesian purity");
+ if(fPIDConfig[6] > 0) fAnalysisSummary->GetXaxis()->SetBinLabel(4, "bayesian purity");
+ else fAnalysisSummary->GetXaxis()->SetBinLabel(4, "TPC TOF PID");
fAnalysisSummary->GetXaxis()->SetBinLabel(5, "dca mode");
fAnalysisSummary->GetXaxis()->SetBinLabel(6, "fVertex");
fAnalysisSummary->GetXaxis()->SetBinLabel(7, "fCentralityMin");
fAnalysisSummary->GetXaxis()->SetBinLabel(8, "fCentralityMax");
- fAnalysisSummary->GetXaxis()->SetBinLabel(9, "V0 centrality");
- fAnalysisSummary->GetXaxis()->SetBinLabel(10, "TPC centrality");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(9, "centrality");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(10, fkCentralityMethod);
fAnalysisSummary->GetXaxis()->SetBinLabel(11, "fApplyDeltaDipCut");
fAnalysisSummary->GetXaxis()->SetBinLabel(12, "POI filterbit");
fAnalysisSummary->GetXaxis()->SetBinLabel(13, "RP filterbit");
fAnalysisSummary->GetXaxis()->SetBinLabel(22, "harm");
fAnalysisSummary->GetXaxis()->SetBinLabel(23, "highPtMode");
fAnalysisSummary->GetXaxis()->SetBinLabel(24, "deltaDipAngle");
- fAnalysisSummary->GetXaxis()->SetBinLabel(25, "POImaxPt");
- fAnalysisSummary->GetXaxis()->SetBinLabel(26, "iterator");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(26, "SpeciesA");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(27, "chargeA");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(28, "fMinPtA");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(29, "fMaxPtA");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(30, "SpeciesB");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(31, "chargeB");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(32, "fMinPtB");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(33, "fMaxPtB");
+ fAnalysisSummary->GetXaxis()->SetBinLabel(34, "iterator");
fOutputList->Add(fAnalysisSummary);
fAnalysisSummary->SetBinContent(1, (Float_t)fIsMC);
fAnalysisSummary->SetBinContent(2, (Float_t)fEventMixing);
fAnalysisSummary->SetBinContent(3, (Float_t)1);
- fAnalysisSummary->SetBinContent(4, (Float_t)fPIDConfig[6]);
+ if(fPIDConfig[6] > 0) fAnalysisSummary->SetBinContent(4, (Float_t)fPIDConfig[6]);
fAnalysisSummary->SetBinContent(5, (Float_t)fDCAConfig[0]);
fAnalysisSummary->SetBinContent(6, (Float_t)fVertexRange);
fAnalysisSummary->SetBinContent(7, (Float_t)fCentralityMin);
fAnalysisSummary->SetBinContent(8, (Float_t)fCentralityMax);
- const char* v = "V0M"; const char* t = "TRK"; // dont compare string literals
- if(fkCentralityMethod==v) fAnalysisSummary->SetBinContent(9, (Float_t)1);
- if(fkCentralityMethod==t) fAnalysisSummary->SetBinContent(10, (Float_t)1);
fAnalysisSummary->SetBinContent(11, (Float_t)fApplyDeltaDipCut);
fAnalysisSummary->SetBinContent(12, (Float_t)fPOICuts->GetAODFilterBit());
fAnalysisSummary->SetBinContent(13, (Float_t)fCutsRP->GetAODFilterBit());
- for(Int_t i(0); i<12;i++) fAnalysisSummary->SetBinContent(14+i, fAddTaskMacroSummary[i]);
- fAnalysisSummary->SetBinContent(26, (Float_t)1);
+ for(Int_t i(0); i<13;i++) fAnalysisSummary->SetBinContent(14+i, fAddTaskMacroSummary[i]);
+ fAnalysisSummary->SetBinContent(26, (Float_t)fSpeciesA);
+ fAnalysisSummary->SetBinContent(27, (Float_t)fChargeA);
+ fAnalysisSummary->SetBinContent(28, (Float_t)fMinPtA);
+ fAnalysisSummary->SetBinContent(29, (Float_t)fMaxPtA);
+ fAnalysisSummary->SetBinContent(30, (Float_t)fSpeciesB);
+ fAnalysisSummary->SetBinContent(31, (Float_t)fChargeB);
+ fAnalysisSummary->SetBinContent(32, (Float_t)fMinPtB);
+ fAnalysisSummary->SetBinContent(33, (Float_t)fMaxPtB);
+ fAnalysisSummary->SetBinContent(34, (Float_t)1);
PostData(1, fOutputList);
// create candidate array
fCandidates = new TObjArray(1000);
AliEventPoolManager* AliAnalysisTwoParticleResonanceFlowTask::InitializeEventMixing()
{
// initialize event mixing
- if(fDebug > 0) cout << " *** InitializeEventMixing() *** " << endl;
Int_t _c(0), _v(0);
for(Int_t i(0); i < 19; i++) {
if (fCentralityMixingBins[i+1] < fCentralityMixingBins[i]) { _c = i; break; }
if (fVertexMixingBins[i+1] < fVertexMixingBins[i]) { _v = i; break; }
else _v = 19;
}
- if(fDebug > 0 ) cout << Form(" --> found %d centrality bins for mixing, %d vertex bins for mixing", _c, _v) << endl;
Float_t centralityBins[_c];
Float_t vertexBins[_v];
for(Int_t i(0); i < _c + 1; i++) centralityBins[i] = fCentralityMixingBins[i];
template <typename T> Float_t AliAnalysisTwoParticleResonanceFlowTask::InvariantMass(const T* track1, const T* track2) const
{
// Return the invariant mass of two tracks
- if(fDebug > 1) cout << " *** InvariantMass() *** " << endl;
if ((!track2) || (!track1)) return 0.;
Float_t pxs = TMath::Power((track1->Px() + track2->Px()), 2);
Float_t pys = TMath::Power((track1->Py() + track2->Py()), 2);
template <typename T> Float_t AliAnalysisTwoParticleResonanceFlowTask::DeltaDipAngle(const T* track1, const T* track2) const
{
// Calculate the delta dip angle between two particles
- if(fDebug > 1) cout << " *** DeltaDipAngle() *** " << endl;
if (track1->P()*track2->P() == 0) return 999;
return TMath::ACos(((track1->Pt() * track2->Pt()) + (track1->Pz() * track2->Pz())) / (track1->P() * track2->P()));
}
template <typename T> Bool_t AliAnalysisTwoParticleResonanceFlowTask::CheckDeltaDipAngle(const T* track1, const T* track2) const
{
// Check if pair passes delta dip angle cut within 0 < p_t < fDeltaDipPt
- if(fDebug > 1) cout << " *** CheckDeltaDipAngle() *** " << endl;
if ((TMath::Abs(DeltaDipAngle(track1, track2)) < fDeltaDipAngle) && (PairPt(track1, track2) < fDeltaDipPt)) return kFALSE;
return kTRUE;
}
template <typename T> Bool_t AliAnalysisTwoParticleResonanceFlowTask::CheckCandidateEtaPtCut(const T* track1, const T* track2) const
{
// Check if pair passes eta and pt cut
- if(fDebug > 1) cout << " *** CheckCandidateEtaPtCut() *** " << endl;
if (fCandidateMinPt > PairPt(track1, track2) || fCandidateMaxPt < PairPt(track1, track2)) return kFALSE;
TVector3 a(track1->Px(), track1->Py(), track1->Pz());
TVector3 b(track2->Px(), track2->Py(), track2->Pz());
template <typename T> Bool_t AliAnalysisTwoParticleResonanceFlowTask::EventCut(T* event)
{
// Impose event cuts
- if(fDebug > 0) cout << " *** EventCut() *** " << endl;
if (!event) return kFALSE;
if (!CheckVertex(event)) return kFALSE;
if (!CheckCentrality(event)) return kFALSE;
template <typename T> void AliAnalysisTwoParticleResonanceFlowTask::PlotMultiplcities(const T* event) const
{
// QA multiplicity plots
- if(fDebug > 1) cout << " *** PlotMultiplcities() *** " << endl;
fVZEROA->Fill(event->GetVZEROData()->GetMTotV0A());
fVZEROC->Fill(event->GetVZEROData()->GetMTotV0C());
fTPCM->Fill(event->GetNumberOfTracks());
template <typename T> Bool_t AliAnalysisTwoParticleResonanceFlowTask::CheckVertex(const T* event)
{
// Check if event vertex is within given range
- if(fDebug > 0) cout << " *** CheckVertex() *** " << endl;
if (!event->GetPrimaryVertex()) return 0x0;
fVertex = event->GetPrimaryVertex()->GetZ();
if (TMath::Abs(fVertex) > fVertexRange) return 0x0;
template <typename T> Bool_t AliAnalysisTwoParticleResonanceFlowTask::CheckCentrality(T* event)
{
// Check if event is within the set centrality range. Falls back to V0 centrality determination if no method is set
- if(fDebug > 0) cout << " *** CheckCentrality() *** " << endl;
if (!fkCentralityMethod) AliFatal("No centrality method set! FATAL ERROR!");
fCentrality = event->GetCentrality()->GetCentralityPercentile(fkCentralityMethod);
if ((fCentrality <= fCentralityMin) || (fCentrality > fCentralityMax)) {
void AliAnalysisTwoParticleResonanceFlowTask::InitializeBayesianPID(AliAODEvent* event)
{
// Initialize the Bayesian PID object for AOD
- if(fDebug > 0) cout << " *** InitializeBayesianPID() *** " << endl;
fBayesianResponse->SetDetResponse(event, fCentrality);
}
//_____________________________________________________________________________
template <typename T> Bool_t AliAnalysisTwoParticleResonanceFlowTask::PassesTPCbayesianCut(T* track, Int_t species) const
{
// Check if the particle passes the TPC TOF bayesian cut.
- if(fDebug > 1) cout << " *** PassesTPCbayesianCut() *** " << endl;
fBayesianResponse->ComputeProb(track);
if (!fBayesianResponse->GetCurrentMask(0)) return kFALSE; // return false if TPC has no response
Int_t contamination(0); //type of particle we don't want in the sample
if(species==3) contamination = 2; //we want to reject pions when selection kaons
if(species==2) contamination = 3; // we want to reject kaons when selecting pions
Float_t *probabilities = fBayesianResponse->GetProb();
- if (probabilities[species] > fPIDConfig[6] && (probabilities[contamination] < probabilities[species])) { // 3 is kaon, 2 is pion
- if(fQA) {fPhi->Fill(track->Phi()); fPt->Fill(track->Pt()); fEta->Fill(track->Eta());}
- return kTRUE;
- }
- return kFALSE;
+ return (probabilities[species] > fPIDConfig[6] && (probabilities[contamination] < probabilities[species]));// 3 is kaon, 2 is pion
}
//_____________________________________________________________________________
Bool_t AliAnalysisTwoParticleResonanceFlowTask::PassesDCACut(AliAODTrack* track) const
// fDCAConfig[0] < -1 no pt dependence
// fDCAConfig[0] = 0 do nothing
// fDCAConfig[0] > 1 pt dependent dca cut
- if(fDebug > 1) cout << " *** PassesDCACut() *** " << endl;
if(fIsMC) return kTRUE;
if( (fDCAConfig[0] < 0.1) && (fDCAConfig[0] > -0.1) ) {
if(fQA) {
return kTRUE;
}
//_____________________________________________________________________________
+Bool_t AliAnalysisTwoParticleResonanceFlowTask::DoOwnPID(AliAODTrack* track, Int_t species) const
+{
+ // do custom pid based on tpc and tof response
+ Float_t nSigmasTPC = (species==fSpeciesA) ? fPIDConfig[0] : fPIDConfig[3]; // number that is used when only tpc is available (5)
+ Float_t nSigmasTPCwithTOF = (species==fSpeciesA) ? fPIDConfig[1] : fPIDConfig[4]; // number that is used for tpc when tof is available (5)
+ Float_t nSigmasTOF = (species==fSpeciesA) ? fPIDConfig[2] : fPIDConfig[5]; // number that is used for tof (if available) (3)
+ if ((track->GetStatus() & AliESDtrack::kTOFout)&&(track->GetStatus() & AliESDtrack::kTIME)) { // check if the track is matched by tof signal
+ if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC(track, (AliPID::EParticleType)species)) > nSigmasTPCwithTOF) return kFALSE; // reject with tpc
+ if (track->P() < 1.5) return (TMath::Abs(fPIDResponse->NumberOfSigmasTOF(track, (AliPID::EParticleType)species)) <= nSigmasTOF);
+ return (TMath::Abs(fPIDResponse->NumberOfSigmasTOF(track, (AliPID::EParticleType)species)) <= (nSigmasTOF-.333*nSigmasTOF));
+ }
+ else { // switch to tpc pid if no tof hit is found
+ if(track->Pt() < .35) return (TMath::Abs(fPIDResponse->NumberOfSigmasTPC(track, (AliPID::EParticleType)species)) <= nSigmasTPC);
+ else if(track->Pt() < .5) return (TMath::Abs(fPIDResponse->NumberOfSigmasTPC(track, (AliPID::EParticleType)species)) < (nSigmasTPC-.4*nSigmasTPC));
+ return (TMath::Abs(fPIDResponse->NumberOfSigmasTPC(track, (AliPID::EParticleType)species)) <= (nSigmasTPC-.6*nSigmasTPC));
+ }
+ return kFALSE;
+}
+//_____________________________________________________________________________
Bool_t AliAnalysisTwoParticleResonanceFlowTask::AcceptTrack(AliAODTrack* track, Int_t species) const
{
- if(fDebug > 1) cout << " *** IsTrack() *** " << endl;
- if (!PassesDCACut(track)) return kFALSE;
- if (PassesTPCbayesianCut(track, species)) {
- if(fQA) if(species==3) fPIDk->Fill(track->P(), track->GetTPCsignal());
- if(fQA) if(species==2) fPIDp->Fill(track->P(), track->GetTPCsignal());
- return kTRUE;
+ if(fQA) TrackQA(track, species, kTRUE); // do qa before cuts
+ if(species==fSpeciesA&&((track->Pt()<fMinPtA)||(track->Pt()>fMaxPtA))) return kFALSE; // least expensive check first
+ if(species==fSpeciesB&&((track->Pt()<fMinPtB)||(track->Pt()>fMaxPtB))) return kFALSE;
+ if(!PassesDCACut(track)) return kFALSE;
+ if(fPIDConfig[6] < 0) {
+ if(DoOwnPID(track, species)) {
+ if(fQA) TrackQA(track, species, kFALSE);
+ return kTRUE;
+ }
+ return kFALSE;
+ }
+ else {
+ if (PassesTPCbayesianCut(track, species)) {
+ if(fQA) TrackQA(track, species, kFALSE);
+ return kTRUE;
+ }
}
return kFALSE;
}
return fPOICuts->IsSelected(track);
}
//_____________________________________________________________________________
+void AliAnalysisTwoParticleResonanceFlowTask::TrackQA(AliAODTrack* track, Int_t species, Bool_t allChargedParticles) const
+{
+ // fill qa histo's
+ if(allChargedParticles) { // track didn't pass identification yet
+ fNOPID->Fill(track->P(), track->GetTPCsignal());
+ if(track->Charge() > 0) {fEventStats->Fill(1); fPtP->Fill(track->Pt());}
+ if(track->Charge() < 0) {fEventStats->Fill(2); fPtN->Fill(track->Pt());}
+ }
+ else if(!allChargedParticles) { // track has been identified
+ if(species==fSpeciesA) {
+ fPtSpeciesA->Fill(track->Pt());
+ fPIDk->Fill(track->P(), track->GetTPCsignal());
+ }
+ else if(species==fSpeciesB) {
+ fPtSpeciesB->Fill(track->Pt());
+ fPIDp->Fill(track->P(), track->GetTPCsignal());
+ }
+ fPhi->Fill(track->Phi());
+ fEta->Fill(track->Eta());
+ }
+}
+//_____________________________________________________________________________
template <typename T> void AliAnalysisTwoParticleResonanceFlowTask::SetNullCuts(T* event)
{
// Set null cuts
- if (fDebug > 0) cout << " *** SetNullCuts() *** " << fCutsRP << endl;
fCutsRP->SetEvent(event, MCEvent());
fNullCuts->SetParamType(AliFlowTrackCuts::kGlobal);
fNullCuts->SetPtRange(+1, -1); // select nothing QUICK
void AliAnalysisTwoParticleResonanceFlowTask::PrepareFlowEvent(Int_t iMulti)
{
// Prepare flow events
- if (fDebug > 0 ) cout << " *** PrepareFlowEvent() *** " << endl;
fFlowEvent->ClearFast();
fFlowEvent->Fill(fCutsRP, fNullCuts);
fFlowEvent->SetReferenceMultiplicity(iMulti);
void AliAnalysisTwoParticleResonanceFlowTask::PhiMinusPsiMethod(TObjArray* MixingCandidates)
{
// extract v2 using the phi minus psi method
- if(fDebug > 0) cout << " ** PhiMinusPsiMethod() *** " << endl;
if(!AliAnalysisTaskVnV0::IsPsiComputed()) { // AliAnalysisTaskVnV0 must be added to analysis que before this task !!!
- if(fDebug > 0 ) cout << "Fatal error: unable to retrieve VZERO task output !!! Exiting ..." << endl;
return;
}
// retrieve data
fSubEventDPhiv2->Fill(0.5, TMath::Cos(2.*(abcPsi2[0]-abcPsi2[1]))); // vzeroa - tpc
fSubEventDPhiv2->Fill(1.5, TMath::Cos(2.*(abcPsi2[0]-abcPsi2[2]))); // vzeroa - vzeroc
fSubEventDPhiv2->Fill(2.5, TMath::Cos(2.*(abcPsi2[1]-abcPsi2[2]))); // tpc - vzeroc
-// if(abcPsi2[0] < 0) abcPsi2[0] += TMath::Pi();
-// if(abcPsi2[2] < 0) abcPsi2[2] += TMath::Pi();
// if event plane stuff went ok, fill the histograms necessary for flow analysis with phi - psi method
- if(fDebug > 0) cout << " *** ReconstructionWithEventMixing() *** " << endl;
AliEventPool* pool = fPoolManager->GetEventPool(fCentrality, fVertex);
if(!pool) AliFatal(Form("No pool found for centrality = %f, zVtx = %f", fCentrality, fVertex));
if(pool->IsReady() || pool->NTracksInPool() > fMixingParameters[1] / 10 || pool->GetCurrentNEvents() >= fMixingParameters[2]) {
Int_t nEvents = pool->GetCurrentNEvents();
- if(fDebug > 0) cout << " --> " << nEvents << " events in mixing buffer ... " << endl;
for (Int_t iEvent(0); iEvent < nEvents; iEvent++) {
TObjArray* mixed_candidates = pool->GetEvent(iEvent);
if(!mixed_candidates) continue; // this should NEVER happen
Int_t bufferTracks = mixed_candidates->GetEntriesFast(); // buffered candidates
Int_t candidates = MixingCandidates->GetEntriesFast(); // mixing candidates
- if(fDebug > 0) cout << Form(" - mixing %d tracks with %d buffered tracks from event %d ... ", candidates, bufferTracks, iEvent) << endl;
TObjArray* SpeciesA = new TObjArray();
SpeciesA->SetOwner(kTRUE);
TObjArray* SpeciesB = new TObjArray();
} // end of mixed events loop
} // end of checking to see whether pool is filled correctly
pool->UpdatePool(MixingCandidates); // update pool with current mixing candidates (for next event)
- if(fDebug > 0 ) pool->PrintInfo();
}
//_____________________________________________________________________________
void AliAnalysisTwoParticleResonanceFlowTask::PhiMinusPsiMethodWriteData(Bool_t signal, TObjArray* SpeciesA, TObjArray* SpeciesB, Float_t* abcPsi2)
{
// write data for phi minus psi method into correct histograms
- if(fDebug > 0 ) cout << " *** PhiMinusPsiMethodWriteData *** " << endl;
for (Int_t i = 0; i < SpeciesA->GetEntriesFast(); i++) { // create pairs
for(Int_t j = 0; j < SpeciesB->GetEntriesFast(); j++) {
AliResonanceFlowHelperTrack* trackA = (AliResonanceFlowHelperTrack*)(SpeciesA->At(i));
void AliAnalysisTwoParticleResonanceFlowTask::VZEROSubEventAnalysis()
{
// vzero event plane analysis using three subevents
- if(fDebug > 0) cout << " ** VZEROSubEventAnalysis() *** " << endl;
if(!AliAnalysisTaskVnV0::IsPsiComputed()) { // AliAnalysisTaskVnV0 must be added to analysis que before this task !!!
- if(fDebug > 0 ) cout << "Fatal error: unable to retrieve VZERO task output !!! Exiting ..." << endl;
return;
}
// retrieve data
Float_t abcPsi2[] = {AliAnalysisTaskVnV0::GetPsi2V0A(), AliAnalysisTaskVnV0::GetPsi2TPC(), AliAnalysisTaskVnV0::GetPsi2V0C()};
for(Int_t i(0); i < 3; i++) if(abcPsi2[i] == 0) {
- if(fDebug > 0) cout << " Warning: I've encountered 0 in one of the symmetry panes (TPC, VZERO) - skipping event !" << endl;
return;
}
// save info for resolution calculations
for(Int_t iCand(0); iCand < fCandidates->GetEntriesFast(); iCand++) { // loop over unlike sign tracks
AliFlowTrackSimple *track = dynamic_cast<AliFlowTrackSimple*>(fCandidates->At(iCand));
if(!track) {
- if(fDebug > 1) cout << " --> dynamic_cast returned null-pointer, skipping candidate <-- " << endl;
continue;
}
for(Int_t mb(0); mb < fMassBins; mb++) { // loop over massbands
void AliAnalysisTwoParticleResonanceFlowTask::UserExec(Option_t *)
{
// UserExec: called for each event. Commented where necessary
- if(fDebug > 0 ) cout << " *** UserExec() *** " << endl;
TObjArray* MixingCandidates = 0x0; // init null pointer for event mixing
if(fEventMixing) {
MixingCandidates = new TObjArray();
ocSpeciesB->SetOwner(kTRUE);
if (!fPIDResponse) { // kill the event if there isn't a pid response
- if(fDebug > 0 ) cout << " Could not get PID response " << endl;
return;
}
for (Int_t i = 0; i < iTracks; i++) { // select analysis candidates
AliAODTrack* track = fAOD->GetTrack(i);
if (!QualityCheck(track)) continue; // reject poor quality tracks
- if (fQA) { // fill qa histo's
- if(track->Charge() > 0) {fEventStats->Fill(1); fPtP->Fill(track->Pt());}
- if(track->Charge() < 0) {fEventStats->Fill(2); fPtN->Fill(track->Pt());}
- }
if (track->Charge() == fChargeA && AcceptTrack(track, fSpeciesA)) { // store species a
SpeciesA->Add(new AliResonanceFlowHelperTrack(track->Eta(), track->Phi(), track->P(), track->Px(), track->Py(), track->Pz(),
track->Pt(), track->Charge(), fMassA, track->GetID(), fSpeciesA));
track->GetID(), fSpeciesB));
}
}
-// cout << SpeciesA->GetEntriesFast() << endl;
-// cout << SpeciesB->GetEntriesFast() << endl;
-// cout << ocSpeciesA->GetEntriesFast() << endl;
-// cout << ocSpeciesB->GetEntriesFast() << endl;
// do the phi minus psi method if that's specified FIXME currenlty only supports event mixing
if (fPhiMinusPsiMethod) { // for the phi minus psi method, no need for flow package etc
PhiMinusPsiMethod(MixingCandidates); // call the method
ResonanceSignal(SpeciesA, SpeciesB);
//
if(fV0) VZEROSubEventAnalysis();
- if (fDebug > 0) printf("I received %d candidates\n", fCandidates->GetEntriesFast()); // insert candidates into flow events
for (int iCand = 0; iCand != fCandidates->GetEntriesFast(); ++iCand) {
AliFlowCandidateTrack *cand = dynamic_cast<AliFlowCandidateTrack*>(fCandidates->At(iCand));
if (!cand) continue;
- if (fDebug > 1) printf(" --> Checking at candidate %d with %d daughters: mass %f\n", iCand, cand->GetNDaughters(), cand->Mass());
for (int iDau = 0; iDau != cand->GetNDaughters(); ++iDau) {
- if (fDebug>1) printf(" *** Daughter %d with fID %d ***", iDau, cand->GetIDDaughter(iDau));
for (int iRPs = 0; iRPs != fFlowEvent->NumberOfTracks(); ++iRPs) {
AliFlowTrack *iRP = dynamic_cast<AliFlowTrack*>(fFlowEvent->GetTrack(iRPs));
if (!iRP) continue;
if (!iRP->InRPSelection()) continue;
if (cand->GetIDDaughter(iDau) == iRP->GetID()) {
- if (fDebug > 1) printf(" was in RP set");
iRP->SetForRPSelection(kFALSE);
fFlowEvent->SetNumberOfRPs(fFlowEvent->GetNumberOfRPs() - 1);
}
}
- if (fDebug > 1) printf("\n");
}
cand->SetForPOISelection(kTRUE);
fFlowEvent->InsertTrack(((AliFlowTrack*) cand));
}
if(!fEventMixing) { // do the combinatorial background
ResonanceBackground(ocSpeciesA, SpeciesB); // mix opposite charge species A with species B
- ResonanceBackground(SpeciesA, ocSpeciesB); // mix species A with opposite charge species B
+ if(fSpeciesA!=fSpeciesB) ResonanceBackground(SpeciesA, ocSpeciesB); // mix species A with opposite charge species B
}
delete SpeciesA;
delete SpeciesB;
void AliAnalysisTwoParticleResonanceFlowTask::ReconstructionWithEventMixing(TObjArray* MixingCandidates) const
{
// perform reconstruction with event mixing
- if(fDebug > 0) cout << " *** ReconstructionWithEventMixing() *** " << endl;
AliEventPool* pool = fPoolManager->GetEventPool(fCentrality, fVertex);
if(!pool) AliFatal(Form("No pool found for centrality = %f, zVtx = %f", fCentrality, fVertex));
if(pool->IsReady() || pool->NTracksInPool() > fMixingParameters[1] / 10 || pool->GetCurrentNEvents() >= fMixingParameters[2]) {
Int_t nEvents = pool->GetCurrentNEvents();
- if(fDebug > 0) cout << " --> " << nEvents << " events in mixing buffer ... " << endl;
for (Int_t iEvent(0); iEvent < nEvents; iEvent++) {
TObjArray* mixed_candidates = pool->GetEvent(iEvent);
if(!mixed_candidates) continue; // this should NEVER happen
Int_t bufferTracks = mixed_candidates->GetEntriesFast(); // buffered candidates
Int_t candidates = MixingCandidates->GetEntriesFast(); // mixing candidates
- if(fDebug > 0) cout << Form(" - mixing %d tracks with %d buffered tracks from event %d ... ", candidates, bufferTracks, iEvent) << endl;
TObjArray* SpeciesA = new TObjArray();
SpeciesA->SetOwner(kTRUE);
TObjArray* SpeciesB = new TObjArray();
} // end of mixed events loop
} // end of checking to see whether pool is filled correctly
pool->UpdatePool(MixingCandidates); // update pool with current mixing candidates (for next event)
- if(fDebug > 0 ) pool->PrintInfo();
}
//_____________________________________________________________________________
void AliAnalysisTwoParticleResonanceFlowTask::Terminate(Option_t *)
{
// terminate
- if(fDebug > 0) cout << " *** Terminate() *** " << endl;
}
//______________________________________________________________________________
void AliAnalysisTwoParticleResonanceFlowTask::MakeTrack(Float_t mass, Float_t pt, Float_t phi, Float_t eta, Int_t nDau, Int_t iID[]) const
{
// Construct Flow Candidate Track from two selected candidates
- if(fDebug > 1 ) cout << " *** MakeTrack() *** " << endl;
Bool_t overwrite = kTRUE;
AliFlowCandidateTrack *sTrack = static_cast<AliFlowCandidateTrack*>(fCandidates->At(fCandidates->GetLast() + 1));
if (!sTrack) {
virtual ~AliAnalysisTwoParticleResonanceFlowTask();
// technical aspects of the analysis
void ForceExit(Int_t type, const char* message);
- Int_t SetDebugLevelResonanceTask(Int_t debug) {fDebug = debug; return fDebug; }
Bool_t SetIsMC(Bool_t ismc) {fIsMC = ismc; return fIsMC; }
void IsMC();
Bool_t UseEventMixing(Bool_t mix) { fEventMixing = mix; return mix; }
void SetAddTaskMacroSummary(Float_t m[12]) {for(Int_t i(0); i < 12; i++) fAddTaskMacroSummary[i] = m[i];}
void SetPOIDCAXYZ(Float_t dca[5]) { for(Int_t i = 0; i < 5; i++) fDCAConfig[i] = dca[i]; }
void SetMixingParameters(Int_t p[3]) { for(Int_t i = 0; i < 3; i++) fMixingParameters[i] = p[i]; }
- void SetupSpeciesA(Int_t species, Int_t charge, Float_t mass) {fSpeciesA = species; fChargeA = charge; fMassA = mass;}
- void SetupSpeciesB(Int_t species, Int_t charge, Float_t mass) {fSpeciesB = species; fChargeB = charge; fMassB = mass;}
+ void SetupSpeciesA(Int_t species, Int_t charge, Float_t mass, Float_t minPtA, Float_t maxPtA) {fSpeciesA = species; fChargeA = charge; fMassA = mass; fMinPtA = minPtA; fMaxPtA = maxPtA;}
+ void SetupSpeciesB(Int_t species, Int_t charge, Float_t mass, Float_t minPtB, Float_t maxPtB) {fSpeciesB = species; fChargeB = charge; fMassB = mass; fMinPtB = minPtB; fMaxPtB = maxPtB;}
void SetCandidateEtaAndPt(Float_t mineta, Float_t maxeta, Float_t minpt, Float_t maxpt) { fCandidateMinEta = mineta;
fCandidateMaxEta = maxeta;
fCandidateMinPt = minpt;
void InitializeBayesianPID(AliAODEvent* event);
template <typename T> Bool_t PassesTPCbayesianCut(T* track, Int_t species) const;
Bool_t PassesDCACut(AliAODTrack* track) const;
+ Bool_t DoOwnPID(AliAODTrack* track, Int_t species) const;
Bool_t AcceptTrack(AliAODTrack* track, Int_t species) const;
template <typename T> Float_t PairPt(const T* track_1, const T* track_2, Bool_t phi = kFALSE) const;
template <typename T> Float_t PtSelector(Int_t _track_type, const T* track_1, const T* track_2, Float_t mass) const;
template <typename T> Bool_t QualityCheck(T* track) const;
+ void TrackQA(AliAODTrack* track, Int_t species, Bool_t allChargedParticles) const;
template <typename T> void SetNullCuts(T* esd);
void PrepareFlowEvent(Int_t iMulti);
void PhiMinusPsiMethod(TObjArray* MixingCandidates);
Int_t fChargeB; // charge for species b
Float_t fMassA; // mass for species a
Float_t fMassB; // mass for species b
- Int_t fDebug; // debug level (0 none, 1 fcn calls, 2 verbose)
+ Float_t fMinPtA; // min pt for species a
+ Float_t fMaxPtA; // max pt for species a
+ Float_t fMinPtB; // min pt for species b
+ Float_t fMaxPtB; // max pt for species b
Bool_t fIsMC; // use mc mode
Bool_t fEventMixing; // use event mixing
Bool_t fPhiMinusPsiMethod; // use phi minus psi method (default is invariant mass fit method)
TH2F *fNOPID;//! QA histogram of TPC response of all charged particles
TH2F *fPIDk;//! QA histogram of TPC response of kaons
TH2F *fPIDp; //! QA histogram of TPC response of pions
- TH2F *fNOPIDTOF; //! QA histo of TOF repsonse charged particles
- TH2F *fPIDTOF; //! QA histo of TOF response kaons
TH1F *fResonanceSignal[18]; //! signal histograms
TH1F *fResonanceBackground[18]; //! like-sign kaon pairs
TH1F *fPtSpectra[18]; //! pt spectra
TH1F *fPtP; //! QA histogram of p_t distribution of positive particles
TH1F *fPtN; //! QA histogram of p_t distribution of negative particles
- TH1F *fPtKP; //! QA histogram of p_t distribution of positive kaons
- TH1F *fPtKN; //! QA histogram of p_t distribution of negative kaons
+ TH1F *fPtSpeciesA; //! QA histogram of p_t distribution of species A
+ TH1F *fPtSpeciesB; //! QA histogram of p_t distribution of species B
Float_t fCentralityMin; // lower bound of cenrality bin
Float_t fCentralityMax; // upper bound of centrality bin
const char *fkCentralityMethod; // method used to determine centrality (V0 by default)
AliFlowTrackCuts *fPOICuts; // cuts for particles of interest (flow package)
Float_t fVertexRange; // absolute value of maximum distance of vertex along the z-axis
- TH1F *fPhi; //! QA plot of azimuthal distribution of tracks used for event plane estimation
- TH1F *fPt; //! QA plot of p_t sectrum of tracks used for event plane estimation
- TH1F *fEta; //! QA plot of eta distribution of tracks used for event plane estimation
+ TH1F *fPhi; //! QA plot of azimuthal distribution of POI daughters
+ TH1F *fEta; //! QA plot of eta distribution of POI daughters
TH1F *fVZEROA; //! QA plot vzeroa multiplicity (all tracks in event)
TH1F *fVZEROC; //! QA plot vzeroc multiplicity (all tracks in event)
TH1F *fTPCM; //! QA plot TPC multiplicity (tracks used for event plane estimation)
AliAnalysisTwoParticleResonanceFlowTask& operator=(const AliAnalysisTwoParticleResonanceFlowTask&); // Not implemented
void MakeTrack(Float_t, Float_t, Float_t, Float_t, Int_t , Int_t[]) const;
- ClassDef(AliAnalysisTwoParticleResonanceFlowTask, 1);
+ ClassDef(AliAnalysisTwoParticleResonanceFlowTask, 2);
};
#endif
Bool_t EP3sub = 0,
Bool_t VZERO_SP = 0, // use vzero sp method
Float_t centrMin = 20., // centrality selection
- Float_t centrMax = 30.,
- Float_t ITSsigma = 0., // pid mode (see task implementation)
- Float_t ITSrange = 0.,
- Float_t TPCcontrol = 1.,
- Float_t TPCsigma = 3.,
- Float_t TPCrange = 0.,
- Float_t ITScontrol = -1.,
- Float_t Bpurity = 0.3,
+ Float_t centrMax = 40.,
+ Float_t TPConlyA = 3.,
+ Float_t TPCTOFA = 5.,
+ Float_t TOFA = 3.,
+ Float_t TPConlyB = 3.,
+ Float_t TPCTOFB = 5.,
+ Float_t TOFB = 3.,
+ Float_t Bpurity = 0.9, // set to 0 to disable bayesian PID
TString suffixName = "UniqueID", // unique id for output objects
Bool_t bCentralTrigger = kTRUE, // trigger selection
Float_t EtaGap = 0., // eta gap for SPSUB
Float_t maxMaxx = 1.13, // max mass for flow analysis
Int_t speciesA = 3, // species a
Int_t speciesB = 2, // species b
- Int_t chargeA = -1, // charge a
- Int_t chargeB = 1, // charge b
+ Int_t chargeA = 1, // charge a
+ Int_t chargeB = -1, // charge b
Float_t massA = 4.93676999999999977e-01, // mass species a
Float_t massB = 1.39570e-01, // mass species b
+ Float_t minPtA = 0.15, // min pt for track a
+ Float_t maxPtA = 10., // max pt for track a
+ Float_t minPtB = 0.15, // min pt for track b
+ Float_t maxPtB = 10, // max pt for track b
Bool_t PhiMinusPsiMethod = 0, // use phi minus psi method
Bool_t event_mixing = kTRUE, // use event mixing
Bool_t highPtMode = kFALSE, // use with caution !!! disables invariant mass fit method
- Float_t deltaPhiMass = 0.0003, // dM in which to look for phi
- Float_t POIPtMax = 5., // max pt of daughterp particles
- Bool_t shrinkSP = kTRUE, // shrink output
- Bool_t debug = kTRUE) // macro debug mode, for task's debug mode see header
+ Float_t mass = 1.092, // mass for high pt analysis
+ Float_t deltaMass = 0.0003) // dM in which to look for resonance
{
+ // removed because CINT cannot cope with more than 40 function arguments
+ Bool_t shrinkSP = kTRUE;
+ Bool_t debug = kFALSE;
/* some defaults that are used frequently:
* for phi:
* - minMass = 0.99
* pion mass = 1.39570e-01
*/
// some defaults that have been removed as function arguments (august 30 2012)
- Float_t POIPtMin = 0.2; // pt of daughters
Float_t deltaDip = 0.;
Float_t deltaDipMaxPt = 0.;
- Bool_t TPCStandAloneTracks = kFALSE; // deprecated
- Bool_t useGlobalRPCuts = kTRUE; // deprecated
Float_t vertexZ = 10.;
Float_t POIEtaMin = -0.8;
Float_t POIEtaMax = 0.8;
// start of macro
- Float_t PIDconfig[] = {ITSsigma, ITSrange, TPCcontrol, TPCsigma, TPCrange, ITScontrol, Bpurity};
+ Float_t PIDconfig[] = {TPConlyA, TPCTOFA, TOFA, TPConlyB, TPCTOFB, TOFB, Bpurity};
// main function, create and add tasks
- if(debug) cout << " === Adding Task PhiFlow === " << endl;
+ if(debug) cout << " === Adding Task TwoParticleResonanceFlow === " << endl;
// set up main output container's name
TString fileName = AliAnalysisManager::GetCommonFileName();
- fileName += ":PhiReconstruction";
+ fileName += ":Reconstruction";
suffixName += Form("%.0f", centrMin);
suffixName += Form("%.0f", centrMax);
fileName+=suffixName;
AddTaskVZERO(0,0,0,0);
}
// create the main task
- AliAnalysisTwoParticleResonanceFlowTask *task = new AliAnalysisTwoParticleResonanceFlowTask("TaskPhiFlow");
- if(debug) cout << " === AliAnalysisTaskPhiFlow === " << task << endl;
+ AliAnalysisTwoParticleResonanceFlowTask *task = new AliAnalysisTwoParticleResonanceFlowTask(suffixName.Data());
+ if(debug) cout << " === AliAnalysisTwoParticleResonanceFlowTask === " << task << endl;
if(!task) {
if(debug) cout << " --> Unexpected error occurred: NO TASK WAS CREATED! (could be a library problem!) " << endl;
return 0x0;
cout << " --> Using phi - Psi method <-- ... " << endl;
Float_t dphibins[] = {0., 0.63, 1.26, 1.89, 2.52, 3.15};
task->SetdPhiBins(dphibins, (Int_t)(sizeof(dphibins)/sizeof(dphibins[1])-1));
- event_mixing = kTRUE;
+ event_mixing = kTRUE;// only implemented for event mixing FIXME
}
- task->SetupSpeciesA(speciesA, chargeA, massA);
- task->SetupSpeciesB(speciesB, chargeB, massB);
+ task->SetupSpeciesA(speciesA, chargeA, massA, minPtA, maxPtA);
+ task->SetupSpeciesB(speciesB, chargeB, massB, minPtB, maxPtB);
if(task->SetVZEROSubEvents(EP3sub)) cout << " --> Setting up VZERO subevents method ... " << endl;
if(event_mixing) {
- if(debug) cout << " --> Enabeling event mixing for phi reconstruction - try at your own risk !!! <-- " << endl;
+ if(debug) cout << " --> Enabeling event mixing for reconstruction - try at your own risk !!! <-- " << endl;
// set vertex and mixing bins - arrays MUST have length 20!
Int_t c[] = {0, 2, 4, 6, 8, 10, 20, 30, 40, 50, 60, 70, 80, 90, 101, 0, 0, 0, 0, 0};
Int_t v[] = {-10, -8, -6, -4, -2, 0, 2, 4, 6, 8, 10, 0, 0, 0, 0, 0, 0, 0, 0, 0};
if(debug) cout << " Fatal error: no RP cuts found, could be a library problem! " << endl;
return 0x0;
}
- if(useGlobalRPCuts&&(!VZERO_SP)) {
+ if(!VZERO_SP) {
AliFlowTrackCuts::trackParameterType rptype = AliFlowTrackCuts::kGlobal;
cutsRP->SetParamType(rptype);
cutsRP->SetPtRange(0.2, 5.0);
return 0x0;
}
cutsPOI = cutsPOI->GetStandardGlobalTrackCuts2010();
- cutsPOI->SetPtRange(POIPtMin, POIPtMax); // pt range of DAUGHTERS !!!
cutsPOI->SetMaxDCAToVertexXY(0.3); // FIXME not implemented in aod086 aod095 see PassesDCACut() in implementation
cutsPOI->SetMaxDCAToVertexZ(0.3);
if(poi_filter < 9999 ) {
POIDCA[0] = 1.; POIDCA[1] = 0.0105; POIDCA[2] = 0.0350; POIDCA[3] = 1.1; POIDCA[4] = 2.;
}
task->SetPOIDCAXYZ(POIDCA);
- Float_t AddTaskMacroSummary[] = {(Float_t)0, (Float_t)SP, (Float_t)SPSUB, (Float_t)QC, (Float_t)EP, (Float_t)EP3sub, (Float_t)VZERO_SP, (Float_t)EtaGap, (Float_t)harm, (Float_t)highPtMode, (Float_t)deltaDip, (Float_t)POIPtMax};
+ Float_t AddTaskMacroSummary[] = {(Float_t)0, (Float_t)SP, (Float_t)SPSUB, (Float_t)QC, (Float_t)EP, (Float_t)EP3sub, (Float_t)VZERO_SP, (Float_t)EtaGap, (Float_t)harm, (Float_t)highPtMode, (Float_t)deltaDip};
task->SetAddTaskMacroSummary(AddTaskMacroSummary);
if(highPtMode) { // high pt loop - loop will end macro
AliFlowTrackSimpleCuts* HighPtSubEventFilterA = new AliFlowTrackSimpleCuts("HighPtSubEventFilterA");
HighPtSubEventFilterA->SetEtaMin(-0.8);
HighPtSubEventFilterA->SetEtaMax(0.0);
- HighPtSubEventFilterA->SetMassMin(1.019445 - deltaPhiMass);
- HighPtSubEventFilterA->SetMassMax(1.019445 + deltaPhiMass);
+ HighPtSubEventFilterA->SetMassMin(mass - deltaMass);
+ HighPtSubEventFilterA->SetMassMax(mass + deltaMass);
AliFlowTrackSimpleCuts* HighPtSubEventFilterB = new AliFlowTrackSimpleCuts("HighPtSubEventFilterB");
HighPtSubEventFilterB->SetEtaMin(0.0);
HighPtSubEventFilterB->SetEtaMax(+0.8);
- HighPtSubEventFilterB->SetMassMin(1.019445 - deltaPhiMass);
- HighPtSubEventFilterB->SetMassMax(1.019445 + deltaPhiMass);
+ HighPtSubEventFilterB->SetMassMin(mass - deltaMass);
+ HighPtSubEventFilterB->SetMassMax(mass + deltaMass);
AliFlowTrackSimpleCuts* HighPtGenericFilter = new AliFlowTrackSimpleCuts("HighPtGenericFilter");
HighPtGenericFilter->SetEtaMin(-0.8);
HighPtGenericFilter->SetEtaMax(+0.8);
- HighPtGenericFilter->SetMassMin(1.019445 - deltaPhiMass);
- HighPtGenericFilter->SetMassMax(1.019445 + deltaPhiMass);
+ HighPtGenericFilter->SetMassMin(mass - deltaMass);
+ HighPtGenericFilter->SetMassMax(mass + deltaMass);
if(debug) cout << " --> Created poi filters " << endl;
// set pair and event cuts
if((deltaDip>0.005)&&(deltaDipMaxPt>0.005)) task->SetMaxDeltaDipAngleAndPt(deltaDip, deltaDipMaxPt);
// specify the PID procedure which will be used
task->SetPIDConfiguration(PIDconfig);
AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
- AliAnalysisDataContainer *coutHist = mgr->CreateContainer(Form("PhiV2_%s", OutputName(centrMin, centrMax,PIDconfig, suffixName, bCentralTrigger, EtaGap, POIEtaMin, POIEtaMax, 0., 15., deltaDip, deltaDipMaxPt, DCA, harm, TPCStandAloneTracks, vertexZ, debug, useGlobalRPCuts).Data()), TList::Class(), AliAnalysisManager::kOutputContainer, fileName);
+ AliAnalysisDataContainer *coutHist = mgr->CreateContainer(Form("%s", TwoParticleResonanceFlow::OutputName(centrMin, centrMax,PIDconfig, suffixName, bCentralTrigger, EtaGap, POIEtaMin, POIEtaMax, deltaDip, deltaDipMaxPt, DCA, harm, vertexZ, debug).Data()), TList::Class(), AliAnalysisManager::kOutputContainer, fileName);
if(debug) cout << " --> Created IO containers " << cinput << ", " << coutHist << endl;
mgr->AddTask(task);
if(debug) cout << " === ADDING MAIN TASK == " << endl;
if(debug) cout << " --> Created IO containers " << flowEvent << endl;
if(debug) cout << " --> suffixName " << suffixName << endl;
if (QC) { // add qc tasks
- AddQCmethod(Form("QCTPCMB_%d_%s", mb, suffixName.Data()), harm, flowEvent, HighPtGenericFilter, debug, 0x0, suffixName.Data());
+ TwoParticleResonanceFlow::AddQCmethod(Form("QCTPCMB_%d_%s", mb, suffixName.Data()), harm, flowEvent, HighPtGenericFilter, debug, 0x0, suffixName.Data());
if(debug) cout << " --> Hanging QC task ... " << mb << " succes! "<< endl;
}
if (SPSUB) { // add sp subevent tasks
- AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -.5*EtaGap, .5*EtaGap, +0.8, "Qa", harm, flowEvent, HighPtSubEventFilterB, NULL, false, shrinkSP, debug, true, suffixName.Data());
+ TwoParticleResonanceFlow::AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -.5*EtaGap, .5*EtaGap, +0.8, "Qa", harm, flowEvent, HighPtSubEventFilterB, NULL, false, shrinkSP, debug, true, suffixName.Data());
if(debug) cout << " --> Hanging SP Qa task " << mb << " succes!" << endl;
- AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -.5*EtaGap, .5*EtaGap, +0.8, "Qb", harm, flowEvent, HighPtSubEventFilterA, NULL, false, shrinkSP, debug, true, suffixName.Data());
+ TwoParticleResonanceFlow::AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -.5*EtaGap, .5*EtaGap, +0.8, "Qb", harm, flowEvent, HighPtSubEventFilterA, NULL, false, shrinkSP, debug, true, suffixName.Data());
if(debug) cout << " --> Hanging SP Qb task ..." << mb << " succes!"<< endl;
}
if (SP) { // add sp tasks
- AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -0.0, +0.0, +0.8, "QaQb", harm, flowEvent, HighPtGenericFilter, NULL, false, shrinkSP, debug, 0x0, suffixName.Data());
+ TwoParticleResonanceFlow::AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -0.0, +0.0, +0.8, "QaQb", harm, flowEvent, HighPtGenericFilter, NULL, false, shrinkSP, debug, 0x0, suffixName.Data());
if(debug) cout << " --> Hanging SP task ... " << mb << " succes!" << endl;
}
if (EP) { // add ep tasks
- AddSPmethod(Form("EPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -0.0, +0.0, +0.8, "QaQb", harm, flowEvent, HighPtGenericFilter, NULL, true, shrinkSP, debug, 0x0, suffixName.Data());
+ TwoParticleResonanceFlow::AddSPmethod(Form("EPTPCMBSP_%d_%s", mb, suffixName.Data()), -0.8, -0.0, +0.0, +0.8, "QaQb", harm, flowEvent, HighPtGenericFilter, NULL, true, shrinkSP, debug, 0x0, suffixName.Data());
if(debug) cout << " --> Hanging EP task ... " << mb << " succes!" << endl;
}
}
// print summary to screen
- cout << endl << endl << " ==== AddTaskPhiFlow launched ===== " << endl;
+ cout << endl << endl << " ==== AddTwoParticleResonanceFlowTask launched ===== " << endl;
cout << " ************ Configured PID routine ************ " << endl;
cout << " 0 < " << PIDconfig[1] << " p_t, ITS || TPC with s < " << PIDconfig[0] << endl;
if(PIDconfig[2] < 0.) cout << " --> TPC control disabled " << endl;
cout << " DCA type: " << DCA;
if (DCA == "") cout << "default";
cout << endl << " ************* Task statisti:q!cs ****************** " << endl;
- cout << " -> Launched PHI reconstruction " << endl;
+ cout << " -> Launched resonance reconstruction " << endl;
if(SP) cout << " --> Launched QaQb SP filters and corresponding SP task " << endl;
if(EP) cout << " --> Launched QaQb QC filters and corresponding EP task " << endl;
if(SPSUB) cout << " --> Launched Qa&&Qb SP filters and corresponding SP tasks " << endl;
cout << " ************************************************ " << endl;
return task;
} // end of high pt loop - high-pt task is ready to go at this point
- Float_t _pt[] = {0.0, 0.3, 0.6, 0.9, 1.2, 1.5, 1.8, 2.1, 2.4, 2.7, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0};
+ Float_t _pt[] = {0.0, 0.6, 1.2, 1.8, 2.4, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0};
task->SetPtBins(_pt, (Int_t)(sizeof(_pt)/sizeof(_pt[1]))-1);
// POI filter cuts, will filter invm mass bands and subevents
AliFlowTrackSimpleCuts* POIfilterQC[30];
cout << " " << PIDconfig[4] << " < 7 p_t, TPC / TOF Bayesian with p < " << PIDconfig[6] << endl;
cout << " ************************************************ " << endl;
}
- if (TPCStandAloneTracks)
- { // switch to tpc standalone tracks for POI selection
- if(debug) cout << " --> Switching to TPC standalone analsis " << endl;
- task->SetRequireStrictKaonCuts();
- task->SetRequireTPCStandAloneKaons();
- }
// create and configure IO containers
AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
- AliAnalysisDataContainer *coutHist = mgr->CreateContainer(Form("PhiV2_%s", OutputName(centrMin, centrMax,PIDconfig, suffixName, bCentralTrigger, EtaGap, POIEtaMin, POIEtaMax, 0., 10., deltaDip, deltaDipMaxPt, DCA, harm, TPCStandAloneTracks, vertexZ, debug, useGlobalRPCuts).Data()), TList::Class(), AliAnalysisManager::kOutputContainer, fileName);
+ AliAnalysisDataContainer *coutHist = mgr->CreateContainer(Form("%s", TwoParticleResonanceFlow::OutputName(centrMin, centrMax,PIDconfig, suffixName, bCentralTrigger, EtaGap, POIEtaMin, POIEtaMax, deltaDip, deltaDipMaxPt, DCA, harm, vertexZ, debug).Data()), TList::Class(), AliAnalysisManager::kOutputContainer, fileName);
if(debug) cout << " --> Created IO containers " << cinput << ", " << coutHist << endl;
mgr->AddTask(task);
if(debug) cout << " === ADDING MAIN TASK == " << endl;
if(debug) cout << " --> suffixName " << suffixName << endl;
for (int mb = 0; mb != 30; ++mb) {
if (QC) { // add qc tasks
- AddQCmethod(Form("QCTPCMB_%d_%s", mb, suffixName.Data()), harm, flowEvent, POIfilterQC[mb], debug, 0x0, suffixName.Data());
+ TwoParticleResonanceFlow::AddQCmethod(Form("QCTPCMB_%d_%s", mb, suffixName.Data()), harm, flowEvent, POIfilterQC[mb], debug, 0x0, suffixName.Data());
if(debug) cout << " --> Hanging QC task ... " << mb << " succes! "<< endl;
}
if (SPSUB) { // add sp subevent tasks
- AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -.5*EtaGap, .5*EtaGap, +0.8, "Qa", harm, flowEvent, POIfilterSP[mb][1], NULL, false, shrinkSP, debug, true, suffixName.Data(), VZERO_SP);
+ TwoParticleResonanceFlow::AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -.5*EtaGap, .5*EtaGap, +0.8, "Qa", harm, flowEvent, POIfilterSP[mb][1], NULL, false, shrinkSP, debug, true, suffixName.Data(), VZERO_SP);
if(debug) cout << " --> Hanging SP Qa task " << mb << " succes!" << endl;
- AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -.5*EtaGap, .5*EtaGap, +0.8, "Qb", harm, flowEvent, POIfilterSP[mb][0], NULL, false, shrinkSP, debug, true, suffixName.Data(), VZERO_SP);
+ TwoParticleResonanceFlow::AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -.5*EtaGap, .5*EtaGap, +0.8, "Qb", harm, flowEvent, POIfilterSP[mb][0], NULL, false, shrinkSP, debug, true, suffixName.Data(), VZERO_SP);
if(debug) cout << " --> Hanging SP Qb task ..." << mb << " succes!"<< endl;
}
if (SP) { // add sp tasks
- AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -0.0, +0.0, +0.8, "QaQb", harm, flowEvent, POIfilterQC[mb], NULL, false, shrinkSP, debug, 0x0, suffixName.Data());
+ TwoParticleResonanceFlow::AddSPmethod(Form("SPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -0.0, +0.0, +0.8, "QaQb", harm, flowEvent, POIfilterQC[mb], NULL, false, shrinkSP, debug, 0x0, suffixName.Data());
if(debug) cout << " --> Hanging SP task ... " << mb << " succes!" << endl;
}
if (EP) { // add ep tasks
- AddSPmethod(Form("EPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -0.0, +0.0, +0.8, "QaQb", harm, flowEvent, POIfilterQC[mb], NULL, true, shrinkSP, debug, 0x0, suffixName.Data());
+ TwoParticleResonanceFlow::AddSPmethod(Form("EPTPCMB_%d_%s", mb, suffixName.Data()), -0.8, -0.0, +0.0, +0.8, "QaQb", harm, flowEvent, POIfilterQC[mb], NULL, true, shrinkSP, debug, 0x0, suffixName.Data());
if(debug) cout << " --> Hanging EP task ... " << mb << " succes!" << endl;
}
}
}
// print summary to screen
- cout << endl << endl << " ==== AddTaskPhiFlow launched ===== " << endl;
+ cout << endl << endl << " ==== AddTwoParticleResonanceFlowTask launched ===== " << endl;
cout << " ************ Configured PID routine ************ " << endl;
cout << " 0 < " << PIDconfig[1] << " p_t, ITS || TPC with s < " << PIDconfig[0] << endl;
if(PIDconfig[2] < 0.) cout << " --> TPC control disabled " << endl;
cout << " DCA type: " << DCA;
if (DCA == "") cout << "default";
cout << endl << " ************* Task statisti:q!cs ****************** " << endl;
- cout << " -> Launched PHI reconstruction " << endl;
+ cout << " -> Launched resonance reconstruction " << endl;
if(SP) cout << " --> Launched 30 QaQb SP filters and corresponding 30 SP tasks " << endl;
if(EP) cout << " --> Launched 30 QaQb QC filters and corresponding 30 EP tasks " << endl;
if(SPSUB) cout << " --> Launched 30+30 Qa&&Qb SP filters and corresponding 30+30 SP tasks " << endl;
cout << " ************************************************ " << endl;
return task;
}
+namespace TwoParticleResonanceFlow {
//_____________________________________________________________________________
-void AddSPmethod(char *name, double minEtaA, double maxEtaA, double minEtaB, double maxEtaB, char *Qvector, int harmonic, AliAnalysisDataContainer *flowEvent, AliFlowTrackSimpleCuts *cutsPOI = NULL, AliFlowTrackSimpleCuts *cutsRFP = NULL, bool bEP, bool shrink = false, bool debug, bool etagap = false, char *suffixName, Bool_t VZERO_SP = kFALSE)
-{
- // add sp task and invm filter tasks
- if(debug) (bEP) ? cout << " ****** Reveived request for EP task ****** " << endl : cout << " ******* Switching to SP task ******* " << endl;
- TString fileName = AliAnalysisManager::GetCommonFileName();
- (bEP) ? fileName+=":EP" : fileName+=":SP";
- fileName+=suffixName;
- if(etagap) {
- fileName+="_SUBEVENTS";
- if(debug) cout << " --> Setting up subevent analysis <-- " << endl;
- }
- if(debug) cout << " --> fileName " << fileName << endl;
- TString myNameSP;
- (bEP) ? myNameSP = Form("%sEPv%d%s", name, harmonic, Qvector): myNameSP = Form("%sSPv%d%s", name, harmonic, Qvector);
- if(debug) cout << " Task and filter name: " << myNameSP << endl;
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- AliAnalysisDataContainer *flowEvent2 = mgr->CreateContainer(Form("Filter_%s", myNameSP.Data()), AliFlowEventSimple::Class(), AliAnalysisManager::kExchangeContainer);
- AliAnalysisTaskFilterFE *tskFilter = new AliAnalysisTaskFilterFE(Form("TaskFilter_%s", myNameSP.Data()), cutsRFP, cutsPOI);
- tskFilter->SetSubeventEtaRange(minEtaA, maxEtaA, minEtaB, maxEtaB);
- if(VZERO_SP) tskFilter->SetSubeventEtaRange(-10, 0, 0, 10);
- mgr->AddTask(tskFilter);
- mgr->ConnectInput(tskFilter, 0, flowEvent);
- mgr->ConnectOutput(tskFilter, 1, flowEvent2);
- AliAnalysisDataContainer *outSP = mgr->CreateContainer(myNameSP.Data(), TList::Class(), AliAnalysisManager::kOutputContainer, fileName);
- AliAnalysisTaskScalarProduct *tskSP = new AliAnalysisTaskScalarProduct(Form("TaskScalarProduct_%s", myNameSP.Data()), kFALSE);
- tskSP->SetApplyCorrectionForNUA(kTRUE);
- tskSP->SetHarmonic(harmonic);
- tskSP->SetTotalQvector(Qvector);
- if (bEP) tskSP->SetBehaveAsEP();
- if (shrink) tskSP->SetBookOnlyBasicCCH(kTRUE);
- mgr->AddTask(tskSP);
- mgr->ConnectInput(tskSP, 0, flowEvent2);
- mgr->ConnectOutput(tskSP, 1, outSP);
-}
+ void AddSPmethod(char *name, double minEtaA, double maxEtaA, double minEtaB, double maxEtaB, char *Qvector, int harmonic, AliAnalysisDataContainer *flowEvent, AliFlowTrackSimpleCuts *cutsPOI = NULL, AliFlowTrackSimpleCuts *cutsRFP = NULL, bool bEP, bool shrink = false, bool debug, bool etagap = false, char *suffixName, Bool_t VZERO_SP = kFALSE)
+ {
+ // add sp task and invm filter tasks
+ if(debug) (bEP) ? cout << " ****** Reveived request for EP task ****** " << endl : cout << " ******* Switching to SP task ******* " << endl;
+ TString fileName = AliAnalysisManager::GetCommonFileName();
+ (bEP) ? fileName+=":EP" : fileName+=":SP";
+ fileName+=suffixName;
+ if(etagap) {
+ fileName+="_SUBEVENTS";
+ if(debug) cout << " --> Setting up subevent analysis <-- " << endl;
+ }
+ if(debug) cout << " --> fileName " << fileName << endl;
+ TString myNameSP;
+ (bEP) ? myNameSP = Form("%sSPasEPv%d%s", name, harmonic, Qvector): myNameSP = Form("%sSPv%d%s", name, harmonic, Qvector);
+ if(debug) cout << " Task and filter name: " << myNameSP << endl;
+ AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+ AliAnalysisDataContainer *flowEvent2 = mgr->CreateContainer(Form("Filter_%s", myNameSP.Data()), AliFlowEventSimple::Class(), AliAnalysisManager::kExchangeContainer);
+ AliAnalysisTaskFilterFE *tskFilter = new AliAnalysisTaskFilterFE(Form("TaskFilter_%s", myNameSP.Data()), cutsRFP, cutsPOI);
+ tskFilter->SetSubeventEtaRange(minEtaA, maxEtaA, minEtaB, maxEtaB);
+ if(VZERO_SP) tskFilter->SetSubeventEtaRange(-10, 0, 0, 10);
+ mgr->AddTask(tskFilter);
+ mgr->ConnectInput(tskFilter, 0, flowEvent);
+ mgr->ConnectOutput(tskFilter, 1, flowEvent2);
+ AliAnalysisDataContainer *outSP = mgr->CreateContainer(myNameSP.Data(), TList::Class(), AliAnalysisManager::kOutputContainer, fileName);
+ AliAnalysisTaskScalarProduct *tskSP = new AliAnalysisTaskScalarProduct(Form("TaskScalarProduct_%s", myNameSP.Data()), kFALSE);
+ tskSP->SetApplyCorrectionForNUA(kTRUE);
+ tskSP->SetHarmonic(harmonic);
+ tskSP->SetTotalQvector(Qvector);
+ if (bEP) tskSP->SetBehaveAsEP();
+ if (shrink) tskSP->SetBookOnlyBasicCCH(kTRUE);
+ mgr->AddTask(tskSP);
+ mgr->ConnectInput(tskSP, 0, flowEvent2);
+ mgr->ConnectOutput(tskSP, 1, outSP);
+ }
//_____________________________________________________________________________
-void AddQCmethod(char *name, int harmonic, AliAnalysisDataContainer *flowEvent, AliFlowTrackSimpleCuts *cutsPOI = NULL, Bool_t debug, AliFlowTrackSimpleCuts *cutsRFP = NULL, char *suffixName)
-{
- // add qc task and invm filter tasks
- if(debug) cout << " ****** Received request for QC v" << harmonic << " task " << name << ", POI " << cutsPOI << ", IO ****** " << flowEvent << endl;
- TString fileName = AliAnalysisManager::GetCommonFileName();
- fileName+=":QC";
- fileName+=suffixName;
- if(debug) cout << " --> Common filename: " << fileName << endl;
- TString myName = Form("%s", name);
- if(debug) cout << " --> myName: " << myName << endl;
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- AliAnalysisDataContainer *flowEvent2 = mgr->CreateContainer(Form("Filter_%s", myName.Data()), AliFlowEventSimple::Class(),AliAnalysisManager::kExchangeContainer);
- AliAnalysisTaskFilterFE *tskFilter = new AliAnalysisTaskFilterFE(Form("TaskFilter_%s", myName.Data()), cutsRFP, cutsPOI);
- mgr->AddTask(tskFilter);
- mgr->ConnectInput(tskFilter, 0, flowEvent);
- mgr->ConnectOutput(tskFilter, 1, flowEvent2);
- AliAnalysisDataContainer *outQC = mgr->CreateContainer(myName.Data(), TList::Class(), AliAnalysisManager::kOutputContainer, fileName);
- AliAnalysisTaskQCumulants *tskQC = new AliAnalysisTaskQCumulants(Form("TaskQCumulants_%s", myName.Data()), kFALSE);
- tskQC->SetApplyCorrectionForNUA(kTRUE);
- tskQC->SetHarmonic(harmonic);
- tskQC->SetBookOnlyBasicCCH(kTRUE);
- mgr->AddTask(tskQC);
- mgr->ConnectInput(tskQC, 0, flowEvent2);
- mgr->ConnectOutput(tskQC, 1, outQC);
-}
+ void AddQCmethod(char *name, int harmonic, AliAnalysisDataContainer *flowEvent, AliFlowTrackSimpleCuts *cutsPOI = NULL, Bool_t debug, AliFlowTrackSimpleCuts *cutsRFP = NULL, char *suffixName)
+ {
+ // add qc task and invm filter tasks
+ if(debug) cout << " ****** Received request for QC v" << harmonic << " task " << name << ", POI " << cutsPOI << ", IO ****** " << flowEvent << endl;
+ TString fileName = AliAnalysisManager::GetCommonFileName();
+ fileName+=":QC";
+ fileName+=suffixName;
+ if(debug) cout << " --> Common filename: " << fileName << endl;
+ TString myName = Form("%s", name);
+ if(debug) cout << " --> myName: " << myName << endl;
+ AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+ AliAnalysisDataContainer *flowEvent2 = mgr->CreateContainer(Form("Filter_%s", myName.Data()), AliFlowEventSimple::Class(),AliAnalysisManager::kExchangeContainer);
+ AliAnalysisTaskFilterFE *tskFilter = new AliAnalysisTaskFilterFE(Form("TaskFilter_%s", myName.Data()), cutsRFP, cutsPOI);
+ mgr->AddTask(tskFilter);
+ mgr->ConnectInput(tskFilter, 0, flowEvent);
+ mgr->ConnectOutput(tskFilter, 1, flowEvent2);
+ AliAnalysisDataContainer *outQC = mgr->CreateContainer(myName.Data(), TList::Class(), AliAnalysisManager::kOutputContainer, fileName);
+ AliAnalysisTaskQCumulants *tskQC = new AliAnalysisTaskQCumulants(Form("TaskQCumulants_%s", myName.Data()), kFALSE);
+ tskQC->SetApplyCorrectionForNUA(kTRUE);
+ tskQC->SetHarmonic(harmonic);
+ tskQC->SetBookOnlyBasicCCH(kTRUE);
+ mgr->AddTask(tskQC);
+ mgr->ConnectInput(tskQC, 0, flowEvent2);
+ mgr->ConnectOutput(tskQC, 1, outQC);
+ }
//_____________________________________________________________________________
-TString OutputName( Float_t centrMin,
- Float_t centrMax,
- Float_t PIDconfig[7],
- TString suffixName,
- Bool_t bCentralTrigger,
- Float_t EtaGap,
- Float_t POIEtaMin,
- Float_t POIEtaMax,
- Float_t POIPtMin,
- Float_t POIPtMax,
- Float_t deltaDip,
- Float_t deltaDipMaxPt,
- TString DCA,
- Int_t harm,
- Bool_t TPCStandAloneTracks,
- Float_t vertexZ,
- Bool_t debug,
- Bool_t useGlobalRPCuts)
-{
- // generate output name
- TString centralityName = "";
- centralityName += suffixName;
- centralityName += "_DCA";
- centralityName += DCA;
- centralityName += Form("_vZ%.f", vertexZ);
- centralityName += "_";
- for(Int_t i = 0; i < 7; i++) centralityName += Form("%.1f_", PIDconfig[i]);
- centralityName += "_POIEta";
- centralityName += Form("%.1f", POIEtaMin);
- centralityName += Form("%.1f", POIEtaMax);
- centralityName += "_gap";
- centralityName += Form("%.1f", -0.5*EtaGap);
- centralityName += Form("%.1f", 0.5*EtaGap);
- centralityName += "_";
- centralityName += Form("dDip%.2f", deltaDip);
- centralityName += "-";
- centralityName += Form("dDipPt%.2f", deltaDipMaxPt);
- if (TPCStandAloneTracks) {
- centralityName += "-";
- centralityName += "TPCStandAloneTracks";
- }
- if (bCentralTrigger) {
- centralityName += "-";
- centralityName += "kMBkCkSC";
- }
- if (!useGlobalRPCuts) {
+ TString OutputName( Float_t centrMin,
+ Float_t centrMax,
+ Float_t PIDconfig[7],
+ TString suffixName,
+ Bool_t bCentralTrigger,
+ Float_t EtaGap,
+ Float_t POIEtaMin,
+ Float_t POIEtaMax,
+ Float_t deltaDip,
+ Float_t deltaDipMaxPt,
+ TString DCA,
+ Int_t harm,
+ Float_t vertexZ,
+ Bool_t debug)
+ {
+ // generate output name
+ TString centralityName = "";
+ centralityName += suffixName;
+ centralityName += "_DCA";
+ centralityName += DCA;
+ centralityName += Form("_vZ%.f", vertexZ);
+ centralityName += "_";
+ for(Int_t i = 0; i < 7; i++) centralityName += Form("%.1f_", PIDconfig[i]);
+ centralityName += "_POIEta";
+ centralityName += Form("%.1f", POIEtaMin);
+ centralityName += Form("%.1f", POIEtaMax);
+ centralityName += "_gap";
+ centralityName += Form("%.1f", -0.5*EtaGap);
+ centralityName += Form("%.1f", 0.5*EtaGap);
+ centralityName += "_";
+ centralityName += Form("dDip%.2f", deltaDip);
centralityName += "-";
- centralityName += "TPCRP";
- }
- if(debug) cout << " --> centralityName " << centralityName << endl;
- return centralityName;
-}
+ centralityName += Form("dDipPt%.2f", deltaDipMaxPt);
+ if (bCentralTrigger) {
+ centralityName += "-";
+ centralityName += "kMBkCkSC";
+ }
+ if(debug) cout << " --> centralityName " << centralityName << endl;
+ return centralityName;
+ }
//_____________________________________________________________________________
-
+} // end of namespace TwoParticleResonanceFlow
if (!directory) continue;
TList* listTemp = directory->GetListOfKeys();
+ if(!listTemp) continue;
for (Int_t icent=0; icent<listTemp->GetEntries(); icent++)
{
TList* list = dynamic_cast<TList*>(directory->Get(listTemp->At(icent)->GetName()));