#include "Riostream.h" //needed as include
#include "TFile.h" //needed as include
#include "TProfile.h" //needed as include
+#include "TProfile2D.h"
#include "TComplex.h" //needed as include
#include "TList.h"
fCommonHistsRes(NULL),
fHistRP(NULL),
fHistProIntFlow(NULL),
+ fHistProDiffFlowPtEtaRP(NULL),
fHistProDiffFlowPtRP(NULL),
fHistProDiffFlowEtaRP(NULL),
+ fHistProDiffFlowPtEtaPOI(NULL),
fHistProDiffFlowPtPOI(NULL),
fHistProDiffFlowEtaPOI(NULL)
{
(fHistProIntFlow->GetXaxis())->SetBinLabel(1,"v_{n}{2}");
fHistProIntFlow->SetYTitle("");
fHistList->Add(fHistProIntFlow);
+
+ fHistProDiffFlowPtEtaRP = new TProfile2D("FlowPro_VPtEtaRP_MCEP","FlowPro_VPtEtaRP_MCEP",iNbinsPt,dPtMin,dPtMax,iNbinsEta,dEtaMin,dEtaMax);
+ fHistProDiffFlowPtEtaRP->SetXTitle("P_{t}");
+ fHistProDiffFlowPtEtaRP->SetYTitle("#eta");
+ fHistList->Add(fHistProDiffFlowPtEtaRP);
fHistProDiffFlowPtRP = new TProfile("FlowPro_VPtRP_MCEP","FlowPro_VPtRP_MCEP",iNbinsPt,dPtMin,dPtMax);
fHistProDiffFlowPtRP->SetXTitle("P_{t}");
fHistProDiffFlowEtaRP->SetYTitle("");
fHistList->Add(fHistProDiffFlowEtaRP);
+ fHistProDiffFlowPtEtaPOI = new TProfile2D("FlowPro_VPtEtaPOI_MCEP","FlowPro_VPtEtaPOI_MCEP",iNbinsPt,dPtMin,dPtMax,iNbinsEta,dEtaMin,dEtaMax);
+ fHistProDiffFlowPtEtaPOI->SetXTitle("P_{t}");
+ fHistProDiffFlowPtEtaPOI->SetYTitle("#eta");
+ fHistList->Add(fHistProDiffFlowPtEtaPOI);
+
fHistProDiffFlowPtPOI = new TProfile("FlowPro_VPtPOI_MCEP","FlowPro_VPtPOI_MCEP",iNbinsPt,dPtMin,dPtMax);
fHistProDiffFlowPtPOI->SetXTitle("P_{t}");
fHistProDiffFlowPtPOI->SetYTitle("");
dEta = pTrack->Eta();
//no-name int. flow (to be improved):
fHistProIntFlow->Fill(0.,dv2);
+ //differential flow (Pt, Eta, RP):
+ fHistProDiffFlowPtEtaRP->Fill(dPt,dEta,dv2,1.);
//differential flow (Pt, RP):
fHistProDiffFlowPtRP->Fill(dPt,dv2,1.);
//differential flow (Eta, RP):
dv2 = TMath::Cos(2*(dPhi-aRP));
dPt = pTrack->Pt();
dEta = pTrack->Eta();
+ //differential flow (Pt, Eta, POI):
+ fHistProDiffFlowPtEtaPOI->Fill(dPt,dEta,dv2,1.);
//differential flow (Pt, POI):
fHistProDiffFlowPtPOI->Fill(dPt,dv2,1.);
//differential flow (Eta, POI):
TProfile* GetHistProIntFlow() {return this->fHistProIntFlow; }
void SetHistProIntFlow(TProfile* aHistProIntFlow)
{this->fHistProIntFlow = aHistProIntFlow; }
+
+ TProfile2D* GetHistProDiffFlowPtEtaRP() {return this->fHistProDiffFlowPtEtaRP; }
+ void SetHistProDiffFlowPtEtaRP(TProfile2D* aHistProDiffFlowPtEtaRP)
+ {this->fHistProDiffFlowPtEtaRP = aHistProDiffFlowPtEtaRP; }
TProfile* GetHistProDiffFlowPtRP() {return this->fHistProDiffFlowPtRP; }
void SetHistProDiffFlowPtRP(TProfile* aHistProDiffFlowPtRP)
TProfile* GetHistProDiffFlowEtaRP() {return this->fHistProDiffFlowEtaRP; }
void SetHistProDiffFlowEtaRP(TProfile* aHistProDiffFlowEtaRP)
{this->fHistProDiffFlowEtaRP = aHistProDiffFlowEtaRP; }
+
+ TProfile2D* GetHistProDiffFlowPtEtaPOI() {return this->fHistProDiffFlowPtEtaPOI; }
+ void SetHistProDiffFlowPtEtaPOI(TProfile2D* aHistProDiffFlowPtEtaPOI)
+ {this->fHistProDiffFlowPtEtaPOI = aHistProDiffFlowPtEtaPOI; }
TProfile* GetHistProDiffFlowPtPOI() {return this->fHistProDiffFlowPtPOI; }
void SetHistProDiffFlowPtPOI(TProfile* aHistProDiffFlowPtPOI)
AliFlowCommonHist* fCommonHists; //
AliFlowCommonHistResults* fCommonHistsRes; //
- TH1F* fHistRP; //reaction plane
- TProfile* fHistProIntFlow; //profile used to calculate the integrated flow of RP particles
- TProfile* fHistProDiffFlowPtRP; //profile used to calculate the differential flow (Pt) of RP particles
- TProfile* fHistProDiffFlowEtaRP; //profile used to calculate the differential flow (Eta) of RP particles
- TProfile* fHistProDiffFlowPtPOI; //profile used to calculate the differential flow (Pt) of POI particles
- TProfile* fHistProDiffFlowEtaPOI; //profile used to calculate the differential flow (Eta) of POI particles
+ TH1F* fHistRP; // reaction plane
+ TProfile* fHistProIntFlow; // profile used to calculate the integrated flow of RP particles
+ TProfile2D* fHistProDiffFlowPtEtaRP; // profile used to calculate the differential flow (Pt,Eta) of RP particles
+ TProfile* fHistProDiffFlowPtRP; // profile used to calculate the differential flow (Pt) of RP particles
+ TProfile* fHistProDiffFlowEtaRP; // profile used to calculate the differential flow (Eta) of RP particles
+ TProfile2D* fHistProDiffFlowPtEtaPOI; // profile used to calculate the differential flow (Pt,Eta) of POI particles
+ TProfile* fHistProDiffFlowPtPOI; // profile used to calculate the differential flow (Pt) of POI particles
+ TProfile* fHistProDiffFlowEtaPOI ; // profile used to calculate the differential flow (Eta) of POI particles
ClassDef(AliFlowAnalysisWithMCEventPlane,1) // Analyse particle distribution versus MC reaction plane
};
#include "Riostream.h" //needed as include
#include "TProfile.h"
+#include "TProfile2D.h"
#include "TList.h"
TProfile *pHistProIntFlow = dynamic_cast<TProfile*>
(fListHistos->FindObject("FlowPro_V_MCEP"));
+
+ TProfile2D *pHistProDiffFlowPtEtaRP = dynamic_cast<TProfile2D*>
+ (fListHistos->FindObject("FlowPro_VPtEtaRP_MCEP"));
TProfile *pHistProDiffFlowPtRP = dynamic_cast<TProfile*>
(fListHistos->FindObject("FlowPro_VPtRP_MCEP"));
TProfile *pHistProDiffFlowEtaRP = dynamic_cast<TProfile*>
(fListHistos->FindObject("FlowPro_VetaRP_MCEP"));
-
+
+ TProfile2D *pHistProDiffFlowPtEtaPOI = dynamic_cast<TProfile2D*>
+ (fListHistos->FindObject("FlowPro_VPtEtaPOI_MCEP"));
+
TProfile *pHistProDiffFlowPtPOI = dynamic_cast<TProfile*>
(fListHistos->FindObject("FlowPro_VPtPOI_MCEP"));
fMcTerm->SetCommonHists(pCommonHists);
fMcTerm->SetCommonHistsRes(pCommonHistResults);
fMcTerm->SetHistProIntFlow(pHistProIntFlow);
+ fMcTerm->SetHistProDiffFlowPtEtaRP(pHistProDiffFlowPtEtaRP);
fMcTerm->SetHistProDiffFlowPtRP(pHistProDiffFlowPtRP);
fMcTerm->SetHistProDiffFlowEtaRP(pHistProDiffFlowEtaRP);
+ fMcTerm->SetHistProDiffFlowPtEtaPOI(pHistProDiffFlowPtEtaPOI);
fMcTerm->SetHistProDiffFlowPtPOI(pHistProDiffFlowPtPOI);
fMcTerm->SetHistProDiffFlowEtaPOI(pHistProDiffFlowEtaPOI);
fMcTerm->Finish();