- //histograms to store the final results
- TH1D *intFlowResults = dynamic_cast<TH1D*>(fListHistos->FindObject("fIntFlowResultsGFC"));
- TH1D *diffFlowResults2 = dynamic_cast<TH1D*>(fListHistos->FindObject("fDiffFlowResults2ndOrderGFC"));
- TH1D *diffFlowResults4 = dynamic_cast<TH1D*>(fListHistos->FindObject("fDiffFlowResults4thOrderGFC"));
- TH1D *diffFlowResults6 = dynamic_cast<TH1D*>(fListHistos->FindObject("fDiffFlowResults6thOrderGFC"));
- TH1D *diffFlowResults8 = dynamic_cast<TH1D*>(fListHistos->FindObject("fDiffFlowResults8thOrderGFC"));
-
- //common histograms to store the final results the integrated and differential flow
- AliFlowCommonHistResults *commonHistRes2nd = dynamic_cast<AliFlowCommonHistResults*>(fListHistos->FindObject("AliFlowCommonHistResults2ndOrderGFC"));
- AliFlowCommonHistResults *commonHistRes4th = dynamic_cast<AliFlowCommonHistResults*>(fListHistos->FindObject("AliFlowCommonHistResults4thOrderGFC"));
- AliFlowCommonHistResults *commonHistRes6th = dynamic_cast<AliFlowCommonHistResults*>(fListHistos->FindObject("AliFlowCommonHistResults6thOrderGFC"));
- AliFlowCommonHistResults *commonHistRes8th = dynamic_cast<AliFlowCommonHistResults*>(fListHistos->FindObject("AliFlowCommonHistResults8thOrderGFC"));
-
- //common control histogram
- AliFlowCommonHist *commonHists = dynamic_cast<AliFlowCommonHist*>(fListHistos->FindObject("AliFlowCommonHistGFC"));
-
- //profiles with average values of generating functions for int. and diff. flow
- TProfile2D *intFlowGenFun = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fIntFlowGenFun"));
-
- TProfile2D *intFlowGenFun4 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fIntFlowGenFun4")); //only for other system of Eq.
- TProfile2D *intFlowGenFun6 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fIntFlowGenFun6")); //only for other system of Eq.
- TProfile2D *intFlowGenFun8 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fIntFlowGenFun8")); //only for other system of Eq.
- TProfile2D *intFlowGenFun16 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fIntFlowGenFun16")); //only for other system of Eq.
-
- //RP, Pt:
- TProfile3D *diffFlowPtRPGenFunRe = dynamic_cast<TProfile3D*>(fListHistos->FindObject("fDiffFlowPtRPGenFunRe"));
- TProfile3D *diffFlowPtRPGenFunIm = dynamic_cast<TProfile3D*>(fListHistos->FindObject("fDiffFlowPtRPGenFunIm"));
- TProfile *ptBinRPNoOfParticles = dynamic_cast<TProfile*>(fListHistos->FindObject("fPtBinRPNoOfParticles"));
-
- //RP, Eta:
- TProfile3D *diffFlowEtaRPGenFunRe = dynamic_cast<TProfile3D*>(fListHistos->FindObject("fDiffFlowEtaRPGenFunRe"));
- TProfile3D *diffFlowEtaRPGenFunIm = dynamic_cast<TProfile3D*>(fListHistos->FindObject("fDiffFlowEtaRPGenFunIm"));
- TProfile *etaBinRPNoOfParticles = dynamic_cast<TProfile*>(fListHistos->FindObject("fEtaBinRPNoOfParticles"));
-
- //POI, Pt:
- TProfile3D *diffFlowPtPOIGenFunRe = dynamic_cast<TProfile3D*>(fListHistos->FindObject("fDiffFlowPtPOIGenFunRe"));
- TProfile3D *diffFlowPtPOIGenFunIm = dynamic_cast<TProfile3D*>(fListHistos->FindObject("fDiffFlowPtPOIGenFunIm"));
- TProfile *ptBinPOINoOfParticles = dynamic_cast<TProfile*>(fListHistos->FindObject("fPtBinPOINoOfParticles"));
-
- //POI, Eta:
- TProfile3D *diffFlowEtaPOIGenFunRe = dynamic_cast<TProfile3D*>(fListHistos->FindObject("fDiffFlowEtaPOIGenFunRe"));
- TProfile3D *diffFlowEtaPOIGenFunIm = dynamic_cast<TProfile3D*>(fListHistos->FindObject("fDiffFlowEtaPOIGenFunIm"));
- TProfile *etaBinPOINoOfParticles = dynamic_cast<TProfile*>(fListHistos->FindObject("fEtaBinPOINoOfParticles"));
-
- //average selected multiplicity (for int. flow)
- TProfile *avMult = dynamic_cast<TProfile*>(fListHistos->FindObject("fAvMultIntFlowGFC"));
-
- TProfile *avMult4 = dynamic_cast<TProfile*>(fListHistos->FindObject("fAvMultIntFlow4GFC")); //only for other system of Eq.
- TProfile *avMult6 = dynamic_cast<TProfile*>(fListHistos->FindObject("fAvMultIntFlow6GFC")); //only for other system of Eq.
- TProfile *avMult8 = dynamic_cast<TProfile*>(fListHistos->FindObject("fAvMultIntFlow8GFC")); //only for other system of Eq.
- TProfile *avMult16 = dynamic_cast<TProfile*>(fListHistos->FindObject("fAvMultIntFlow16GFC")); //only for other system of Eq.
-
- //average values of Q-vector components (1st bin: <Q_x>, 2nd bin: <Q_y>, 3rd bin: <(Q_x)^2>, 4th bin: <(Q_y)^2>)
- TProfile *qVectorComponents = dynamic_cast<TProfile*>(fListHistos->FindObject("fQVectorComponentsGFC"));
-
- /*
- TProfile2D *diffFlowPtGenFunRe0 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunRe0"));
- TProfile2D *diffFlowPtGenFunRe1 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunRe1"));
- TProfile2D *diffFlowPtGenFunRe2 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunRe2"));
- TProfile2D *diffFlowPtGenFunRe3 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunRe3"));
- TProfile2D *diffFlowPtGenFunRe4 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunRe4"));
- TProfile2D *diffFlowPtGenFunRe5 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunRe5"));
- TProfile2D *diffFlowPtGenFunRe6 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunRe6"));
- TProfile2D *diffFlowPtGenFunRe7 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunRe7"));
- TProfile2D *diffFlowPtGenFunIm0 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunIm0"));
- TProfile2D *diffFlowPtGenFunIm1 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunIm1"));
- TProfile2D *diffFlowPtGenFunIm2 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunIm2"));
- TProfile2D *diffFlowPtGenFunIm3 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunIm3"));
- TProfile2D *diffFlowPtGenFunIm4 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunIm4"));
- TProfile2D *diffFlowPtGenFunIm5 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunIm5"));
- TProfile2D *diffFlowPtGenFunIm6 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunIm6"));
- TProfile2D *diffFlowPtGenFunIm7 = dynamic_cast<TProfile2D*>(fListHistos->FindObject("fdiffFlowPtGenFunIm7"));
- */
-
- //profile with avarage selected multiplicity for int. flow
- //TProfile *avMult = dynamic_cast<TProfile*>(fListHistos->FindObject("fAvMultIntFlow"));
-
- //profile with avarage values of Q-vector components (1st bin: <Q_x>, 2nd bin: <Q_y>, 3rd bin: <(Q_x)^2>, 4th bin: <(Q_y)^2>)
- //TProfile *QVectorComponents = dynamic_cast<TProfile*>(fListHistos->FindObject("fQVectorComponents"));
-
- //q-distribution
- //TH1D *qDist = dynamic_cast<TH1D*>(fListHistos->FindObject("fQDist"));
-
- //AliCumulantsFunctions finalResults(intFlowGenFun,NULL,NULL, intFlowResults,diffFlowResults2,diffFlowResults4,diffFlowResults6,diffFlowResults8,avMult,QVectorComponents,qDist,diffFlowPtGenFunRe0,diffFlowPtGenFunRe1,diffFlowPtGenFunRe2, diffFlowPtGenFunRe3,diffFlowPtGenFunRe4,diffFlowPtGenFunRe5,diffFlowPtGenFunRe6,diffFlowPtGenFunRe7,diffFlowPtGenFunIm0,diffFlowPtGenFunIm1, diffFlowPtGenFunIm2,diffFlowPtGenFunIm3,diffFlowPtGenFunIm4,diffFlowPtGenFunIm5,diffFlowPtGenFunIm6,diffFlowPtGenFunIm7);
-
- //AliCumulantsFunctions finalResults(intFlowGenFun,diffFlowPtGenFunRe,diffFlowPtGenFunIm, intFlowResults,diffFlowResults2,diffFlowResults4,diffFlowResults6,diffFlowResults8,avMult,QVectorComponents,qDist);
-
- //finalResults.Calculate();
-
-
-
- //----------------------------------------------------
-
- fGFC = new AliFlowAnalysisWithCumulants();
-
- fGFC->SetIntFlowResults(intFlowResults);
- fGFC->SetDiffFlowResults2nd(diffFlowResults2);
- fGFC->SetDiffFlowResults4th(diffFlowResults4);
- fGFC->SetDiffFlowResults6th(diffFlowResults6);
- fGFC->SetDiffFlowResults8th(diffFlowResults8);
-
- fGFC->SetCommonHistsResults2nd(commonHistRes2nd);
- fGFC->SetCommonHistsResults4th(commonHistRes4th);
- fGFC->SetCommonHistsResults6th(commonHistRes6th);
- fGFC->SetCommonHistsResults8th(commonHistRes8th);
-
- fGFC->SetCommonHists(commonHists);
-
- fGFC->SetIntFlowGenFun(intFlowGenFun);
-
- fGFC->SetIntFlowGenFun4(intFlowGenFun4); //only for other system of Eq.
- fGFC->SetIntFlowGenFun6(intFlowGenFun6); //only for other system of Eq.
- fGFC->SetIntFlowGenFun8(intFlowGenFun8); //only for other system of Eq.
- fGFC->SetIntFlowGenFun16(intFlowGenFun16); //only for other system of Eq.
-
- fGFC->SetDiffFlowPtRPGenFunRe(diffFlowPtRPGenFunRe);
- fGFC->SetDiffFlowPtRPGenFunIm(diffFlowPtRPGenFunIm);
- fGFC->SetNumberOfParticlesPerPtBinRP(ptBinRPNoOfParticles);
-
- fGFC->SetDiffFlowEtaRPGenFunRe(diffFlowEtaRPGenFunRe);
- fGFC->SetDiffFlowEtaRPGenFunIm(diffFlowEtaRPGenFunIm);
- fGFC->SetNumberOfParticlesPerEtaBinRP(etaBinRPNoOfParticles);
-
- fGFC->SetDiffFlowPtPOIGenFunRe(diffFlowPtPOIGenFunRe);
- fGFC->SetDiffFlowPtPOIGenFunIm(diffFlowPtPOIGenFunIm);
- fGFC->SetNumberOfParticlesPerPtBinPOI(ptBinPOINoOfParticles);
-
- fGFC->SetDiffFlowEtaPOIGenFunRe(diffFlowEtaPOIGenFunRe);
- fGFC->SetDiffFlowEtaPOIGenFunIm(diffFlowEtaPOIGenFunIm);
- fGFC->SetNumberOfParticlesPerEtaBinPOI(etaBinPOINoOfParticles);
-
- fGFC->SetAverageMultiplicity(avMult);
-
- fGFC->SetAverageMultiplicity4(avMult4); //only for other system of Eq.
- fGFC->SetAverageMultiplicity6(avMult6); //only for other system of Eq.
- fGFC->SetAverageMultiplicity8(avMult8); //only for other system of Eq.
- fGFC->SetAverageMultiplicity16(avMult16); //only for other system of Eq.
-
- fGFC->SetQVectorComponents(qVectorComponents);
-