Double_t two = fQCorr->GetBinContent(1); //<<2>>_{n|n}
Double_t four = fQCorr->GetBinContent(11); //<<4>>_{n,n|n,n}
Double_t six = fQCorr->GetBinContent(24); //<<6>>_{n,n,n|n,n,n}
- Double_t eight = fQCorr->GetBinContent(31); //<<8>>_{n,n,n,n|n,n,n,n}
+ Double_t eight = fQCorr->GetBinContent(31); //<<8>>_{n,n,n,n|n,n,n,n}
//errors of 2-, 4- and 6-particle azimuthal correlation:
Double_t twoErr = fQCorr->GetBinError(1); //sigma_{<<2>>_{n|n}}
Double_t fourErr = fQCorr->GetBinError(11); //sigma_{<<4>>_{n,n|n,n}}
- Double_t sixErr = fQCorr->GetBinError(21); //sigma_{<<6>>_{n,n,n|n,n,n}}
+ Double_t sixErr = fQCorr->GetBinError(24); //sigma_{<<6>>_{n,n,n|n,n,n}}
//Double_t eightErr = fQCorr->GetBinError(); //sigma_{<<8>>_{n,n,n,n|n,n,n,n}}
//covariances of multi-particle correlations
cout<<endl;
cout<<"**************************************"<<endl;
cout<<"**************************************"<<endl;
- cout<<"flow estimates from Q-cumulants (RP):"<<endl;
+ cout<<"flow estimates from Q-cumulants :"<<endl;
cout<<endl;
Double_t vn2=0.,vn4=0.,vn6=0.,vn8=0.;
fIntRes->SetBinContent(1,vn2);
fIntRes->SetBinError(1,sd2);
//common histograms:
- fchr2nd->FillIntegratedFlow(vn2,sd2);//to be removed
+ fchr2nd->FillIntegratedFlow(vn2,sd2);
fchr2nd->FillChi(vn2*pow(AvM,0.5));//to be removed
- fchr2nd->FillIntegratedFlowRP(vn2,sd2);
}else{
cout<<" v_"<<n<<"{2} = Im"<<endl;
}
fIntRes->SetBinContent(2,vn4);
fIntRes->SetBinError(2,sd4);
//common histograms:
- fchr4th->FillIntegratedFlow(vn4,sd4);//to be removed
+ fchr4th->FillIntegratedFlow(vn4,sd4);
fchr4th->FillChi(vn4*pow(AvM,0.5));//to be removed
- fchr4th->FillIntegratedFlowRP(vn4,sd4);
}else{
cout<<" v_"<<n<<"{4} = Im"<<endl;
}
fIntRes->SetBinContent(3,vn6);
fIntRes->SetBinError(3,sd6);
//common histograms:
- fchr6th->FillIntegratedFlow(vn6,sd6);//to be removed
+ fchr6th->FillIntegratedFlow(vn6,sd6);
fchr6th->FillChi(vn6*pow(AvM,0.5));//to be removed
- fchr6th->FillIntegratedFlowRP(vn6,sd6);
}else{
cout<<" v_"<<n<<"{6} = Im"<<endl;
}
fIntRes->SetBinContent(4,vn8);
fIntRes->SetBinError(4,sd8);
//common histograms:
- fchr8th->FillIntegratedFlow(vn8,sd8);//to be removed
+ fchr8th->FillIntegratedFlow(vn8,sd8);
fchr8th->FillChi(vn8*pow(AvM,0.5));//to be removed
- fchr8th->FillIntegratedFlowRP(vn8,sd8);
}else{
cout<<" v_"<<n<<"{8} = Im"<<endl;
}
cout<<endl;
//---------------------------------------------------------------------------------------------------------
-//---------------------------------------------------------------------------------------------------------
+
+
+//---------------------------------------------------------------------------------------------------------
//differential flow (RP)
+Int_t nBinsPtRP = fbinPt2p1n1nRP->GetNbinsX();
+Int_t nBinsEtaRP = fbinEta2p1n1nRP->GetNbinsX();
//Pt:
Double_t secondOrderQCumulantDiffFlowPtRP = 0.;
Double_t fourthOrderQCumulantDiffFlowPtRP = 0.;
-Int_t nBinsPtRP = fbinPt2p1n1nRP->GetNbinsX();
+Double_t dVn2ndRP=0.,dSd2ndRP=0.,dDiffvn2ndRP=0.,dYield2ndRP=0.,dSum2ndRP=0.;
+Double_t dVn4thRP=0.,dSd4thRP=0.,dDiffvn4thRP=0.,dYield4thRP=0.,dSum4thRP=0.;
for(Int_t bb=1;bb<nBinsPtRP+1;bb++)
{
+ //QC{2}
if(fbinPt2p1n1nRP->GetBinEntries(bb)>0.&&vn2!=0)
{
+ //cout<<"bin = "<<bb<<" : "<<(fch2nd->GetHistPtDiff())->GetBinContent(bb)<<endl;
+ //cout<<endl;
secondOrderQCumulantDiffFlowPtRP = fbinPt2p1n1nRP->GetBinContent(bb);
fDiffRes2nd->SetBinContent(bb,secondOrderQCumulantDiffFlowPtRP/vn2);
//common histogram:
- fchr2nd->FillDifferentialFlowPtRP(bb,secondOrderQCumulantDiffFlowPtRP/vn2, 0.);//to be improved (errors)
+ fchr2nd->FillDifferentialFlowPtRP(bb,secondOrderQCumulantDiffFlowPtRP/vn2, 0.);//to be improved (errors) && bb or bb+1
+ //-------------------------------------------------------------
+ //integrated flow (RP, Pt, 2nd order):
+ dDiffvn2ndRP=(fchr2nd->GetHistDiffFlowPtRP())->GetBinContent(bb);
+ dYield2ndRP=(fch2nd->GetHistPtInt())->GetBinContent(bb);
+ dVn2ndRP+=dDiffvn2ndRP*dYield2ndRP;
+ dSum2ndRP+=dYield2ndRP;
+ //-------------------------------------------------------------
}
+ //QC{4]
if(fbinPt4p1n1n1n1nRP->GetBinEntries(bb)>0.&&vn4!=0.)
{
fourthOrderQCumulantDiffFlowPtRP = fbinPt4p1n1n1n1nRP->GetBinContent(bb)-2.*fbinPt2p1n1nRP->GetBinContent(bb)*pow(vn2,2.);
fDiffRes4th->SetBinContent(bb,-1.*fourthOrderQCumulantDiffFlowPtRP/pow(vn4,3.));
//common histogram:
fchr4th->FillDifferentialFlowPtRP(bb,-1.*fourthOrderQCumulantDiffFlowPtRP/pow(vn4,3.), 0.);//to be improved (errors)
+ //-------------------------------------------------------------
+ //integrated flow (RP, Pt, 4th order):
+ dDiffvn4thRP=(fchr4th->GetHistDiffFlowPtRP())->GetBinContent(bb);
+ dYield4thRP=(fch4th->GetHistPtInt())->GetBinContent(bb);
+ dVn4thRP+=dDiffvn4thRP*dYield4thRP;
+ dSum4thRP+=dYield4thRP;
+ //-------------------------------------------------------------
}
-}
+}
+
+//integrated flow estimates from Q-cumulants (RP):
+cout<<endl;
+cout<<"**************************************"<<endl;
+cout<<"**************************************"<<endl;
+cout<<"flow estimates from Q-cumulants (RP):"<<endl;
+cout<<endl;
+//storing the final results for integrated flow (RP):
+//QC{2}
+if(dSum2ndRP && fchr2nd)
+{
+ dVn2ndRP/=dSum2ndRP;
+ fchr2nd->FillIntegratedFlowRP(dVn2ndRP,0.);//to be improved (errors)
+ cout<<" v_"<<n<<"{2} = "<<dVn2ndRP<<" +/- "<<dSd2ndRP<<endl;
+}else
+ {
+ cout<<" v_"<<n<<"{2} = Im"<<endl;
+ }
+
+//QC{4}
+if(dSum4thRP && fchr4th)
+{
+ dVn4thRP/=dSum4thRP;
+ fchr4th->FillIntegratedFlowRP(dVn4thRP,0.);//to be improved (errors)
+ cout<<" v_"<<n<<"{4} = "<<dVn4thRP<<" +/- "<<dSd4thRP<<endl;
+}else
+ {
+ cout<<" v_"<<n<<"{4} = Im"<<endl;
+ }
+
+cout<<endl;
+cout<<" nEvts = "<<nEvts<<", AvM = "<<AvM<<endl;
+cout<<"**************************************"<<endl;
+cout<<"**************************************"<<endl;
+cout<<endl;
//Eta:
Double_t secondOrderQCumulantDiffFlowEtaRP = 0.;
Double_t fourthOrderQCumulantDiffFlowEtaRP = 0.;
-Int_t nBinsEtaRP = fbinEta2p1n1nRP->GetNbinsX();
-
for(Int_t bb=1;bb<nBinsEtaRP+1;bb++)
{
if(fbinEta2p1n1nRP->GetBinEntries(bb)>0.&&vn2!=0)
//common histogram:
fchr4th->FillDifferentialFlowEtaRP(bb,-1.*fourthOrderQCumulantDiffFlowEtaRP/pow(vn4,3.), 0.);//to be improved (errors)
}
-}
-
+}
+//---------------------------------------------------------------------------------------------------
+
+
+
+
+//---------------------------------------------------------------------------------------------------------
//differential flow (POI)
Int_t nBinsPtPOI = fbinPt2p1n1nPOI->GetNbinsX();
Int_t nBinsEtaPOI = fbinEta2p1n1nPOI->GetNbinsX();
Double_t secondOrderQCumulantDiffFlowPtPOI = 0.;
Double_t fourthOrderQCumulantDiffFlowPtPOI = 0.;
-Double_t dVn2ndPOI=0.,dSd2ndPOI=0.,dDiffvn2nd=0.,dYield2nd=0.,dSum2nd=0.;
-Double_t dVn4thPOI=0.,dSd4thPOI=0.,dDiffvn4th=0.,dYield4th=0.,dSum4th=0.;
-
-Double_t dVn2ndPOIEta=0.,dDiffvn2ndEta=0.,dYield2ndEta=0.,dSum2ndEta=0.;//to be removed
+Double_t dVn2ndPOI=0.,dSd2ndPOI=0.,dDiffvn2ndPOI=0.,dYield2ndPOI=0.,dSum2ndPOI=0.;
+Double_t dVn4thPOI=0.,dSd4thPOI=0.,dDiffvn4thPOI=0.,dYield4thPOI=0.,dSum4thPOI=0.;
for(Int_t bb=1;bb<nBinsPtPOI+1;bb++)
{
fchr2nd->FillDifferentialFlowPtPOI(bb,secondOrderQCumulantDiffFlowPtPOI/vn2, 0.);//to be improved (errors) && bb or bb+1
//-------------------------------------------------------------
//integrated flow (POI, Pt, 2nd order):
- dDiffvn2nd=(fchr2nd->GetHistDiffFlowPtPOI())->GetBinContent(bb);
- dYield2nd=(fch2nd->GetHistPtDiff())->GetBinContent(bb);
- dVn2ndPOI+=dDiffvn2nd*dYield2nd;
- dSum2nd+=dYield2nd;
+ dDiffvn2ndPOI=(fchr2nd->GetHistDiffFlowPtPOI())->GetBinContent(bb);
+ dYield2ndPOI=(fch2nd->GetHistPtDiff())->GetBinContent(bb);
+ dVn2ndPOI+=dDiffvn2ndPOI*dYield2ndPOI;
+ dSum2ndPOI+=dYield2ndPOI;
//-------------------------------------------------------------
}
//QC{4]
fchr4th->FillDifferentialFlowPtPOI(bb,-1.*fourthOrderQCumulantDiffFlowPtPOI/pow(vn4,3.), 0.);//to be improved (errors)
//-------------------------------------------------------------
//integrated flow (POI, Pt, 4th order):
- dDiffvn4th=(fchr4th->GetHistDiffFlowPtPOI())->GetBinContent(bb);
- dYield4th=(fch4th->GetHistPtDiff())->GetBinContent(bb);
- dVn4thPOI+=dDiffvn4th*dYield4th;
- dSum4th+=dYield4th;
+ dDiffvn4thPOI=(fchr4th->GetHistDiffFlowPtPOI())->GetBinContent(bb);
+ dYield4thPOI=(fch4th->GetHistPtDiff())->GetBinContent(bb);
+ dVn4thPOI+=dDiffvn4thPOI*dYield4thPOI;
+ dSum4thPOI+=dYield4thPOI;
//-------------------------------------------------------------
}
}
-//---------------------------------------------------------------------------------------------------
-//integrated flow estimates from Q-cumulants (POI):
cout<<endl;
cout<<"**************************************"<<endl;
cout<<"**************************************"<<endl;
cout<<endl;
//storing the final results for integrated flow (POI):
//QC{2}
-if(dSum2nd && fchr2nd)
+if(dSum2ndPOI && fchr2nd)
{
- dVn2ndPOI/=dSum2nd;
+ dVn2ndPOI/=dSum2ndPOI;
fchr2nd->FillIntegratedFlowPOI(dVn2ndPOI,0.);//to be improved (errors)
cout<<" v_"<<n<<"{2} = "<<dVn2ndPOI<<" +/- "<<dSd2ndPOI<<endl;
}else
}
//QC{4}
-if(dSum4th && fchr4th)
+if(dSum4thPOI && fchr4th)
{
- dVn4thPOI/=dSum4th;
+ dVn4thPOI/=dSum4thPOI;
fchr4th->FillIntegratedFlowPOI(dVn4thPOI,0.);//to be improved (errors)
cout<<" v_"<<n<<"{4} = "<<dVn4thPOI<<" +/- "<<dSd4thPOI<<endl;
}else
cout<<"**************************************"<<endl;
cout<<"**************************************"<<endl;
cout<<endl;
-//---------------------------------------------------------------------------------------------------
-
//Eta:
Double_t secondOrderQCumulantDiffFlowEtaPOI = 0.;
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-//xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-// !!!! to be removed !!!!
-for(Int_t bb=1;bb<nBinsEtaPOI+1;bb++)
-{
- //integrated flow (POI, Eta, 2nd order):
- dDiffvn2ndEta=(fchr2nd->GetHistDiffFlowEtaPOI())->GetBinContent(bb);
- dYield2ndEta=(fch2nd->GetHistEtaDiff())->GetBinContent(bb);
- // cout<<"bin = "<<bb<<" : "<<dDiffvn2ndEta<<" "<<dYield2ndEta<<endl;
- // cout<<endl;
- dVn2ndPOIEta+=dDiffvn2ndEta*dYield2ndEta;
- dSum2ndEta+=dYield2ndEta;
-
-
-//check why for small number of events there is a mismatch here (should be exactly the same)://to be removed
-//cout<<dSum2nd<<" to be improved "<<(fch2nd->GetHistPtDiff())->Integral()<<endl;//to be removed
-}
-if(dSum2ndEta)
-{
- dVn2ndPOIEta/=dSum2ndEta;
- //cout<<"POI Pt vs POI Eta:"<<endl;
- //cout<<dVn2ndPOI<<" vs "<<dVn2ndPOIEta<<endl;
- //cout<<dSum2nd<<" vs "<<dSum2ndEta<<endl;
-// cout<<"eta = "<<(fchr2nd->GetHistDiffFlowEtaPOI())->GetBinContent(44)<<endl;
-// cout<<"poi = "<<(fchr2nd->GetHistDiffFlowPtPOI())->GetBinContent(14)<<endl;
-// cout<<"+++++++++"<<endl;
-}
-
-//xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- /*
- Double_t first=16.*pow(two*twoErr,2.);
- Double_t second=pow(fourErr,2.);
- Double_t third=-8.*two*cov24;
-
- cout<<endl;
- cout<<endl;
- cout<<endl;
- cout<<" 1st = "<<(1./16.)*first/(pow(-fourthOrderQCumulant,1.5))<<endl;
- cout<<" 2nd = "<<(1./16.)*second/(pow(-fourthOrderQCumulant,1.5))<<endl;
- cout<<" 3rd = "<<(1./16.)*third/(pow(-fourthOrderQCumulant,1.5))<<endl;
- cout<<endl;
- cout<<" nEvts = "<<nEvts<<", AvM = "<<AvM<<endl;
- cout<<endl;
- cout<<endl;
- cout<<"c-c = "<<pow((1./16.)*((first+second+third)/(pow(-fourthOrderQCumulant,1.5))),0.5)<<endl;
- */
-
-
-
- /*
-
- //<Q-vector components>
- Double_t AvQx = fQVector->GetBinContent(1); //<Q_x>
- Double_t AvQy = fQVector->GetBinContent(2); //<Q_y>
- Double_t AvQ2x = fQVector->GetBinContent(3); //<(Q_x)^2>
- Double_t AvQ2y = fQVector->GetBinContent(4); //<(Q_y)^2>
- Double_t AvQ2 = fQVector->GetBinContent(5); //<|Q|^2>
- Double_t AvQ4 = fQVector->GetBinContent(6); //<|Q|^4>
- Double_t AvQ_2n2 = fQVector->GetBinContent(7); //<|Q_2n|^2>
- Double_t ReQ2nQnstarQnstar = fQVector->GetBinContent(8); //<Re[Q_n^2 * Q_2n^*]>
- Double_t ImQ2nQnstarQnstar = fQVector->GetBinContent(9); //<Im[Q_n^2 * Q_2n^*]>
- Double_t AvQ_3n2 = fQVector->GetBinContent(10); //<|Q_3n|^2>
- Double_t AvQ_4n2 = fQVector->GetBinContent(11); //<|Q_4n|^2>
- Double_t AvQ_5n2 = fQVector->GetBinContent(12); //<|Q_5n|^2>
- Double_t ReQ3nQ2nstarQnstar = fQVector->GetBinContent(13); //Re<Q_{3n} Q_{2n}^* Q_{n}^*>
- Double_t ReQ3nQnstarQnstarQnstar = fQVector->GetBinContent(14); //Re<Q_{3n} Q_{n}^* Q_{n}^* Q_{n}^*>
- Double_t HereQ2nQnQ2nstarQnstar = fQVector->GetBinContent(15); //<|Q_{2n}|^2 |Q_{n}|^2>
- Double_t HereQ2nQnQnstarQnstarQnstar = fQVector->GetBinContent(16); //Re<Q_{2n} Q_{n} Q_{n}^* Q_{n}^* Q_{n}^*>
- Double_t AvQ6 = fQVector->GetBinContent(17); //<|Q|^6>
-
-
-
-
-*/
-
- /*
- //direct particle correlations
- Double_t twoDirect_n_n = fDirect->GetBinContent(1);
- Double_t fourDirect = fDirect->GetBinContent(2);
- Double_t threeDirect_2n_n_n = fDirect->GetBinContent(3);
- Double_t threeDirect_3n_2n_n = fDirect->GetBinContent(4);
- Double_t fourDirect_3n_n_n_n = fDirect->GetBinContent(5);
- Double_t fourDirect_2n_n_2n_n = fDirect->GetBinContent(6);
- Double_t fiveDirect_2n_n_n_n_n = fDirect->GetBinContent(7);
- Double_t sixDirect_n_n_n_n_n_n = fDirect->GetBinContent(8);
- Double_t twoDirect_2n_2n = fDirect->GetBinContent(9);
- Double_t twoDirect_3n_3n = fDirect->GetBinContent(10);
- Double_t twoDirect_4n_4n = fDirect->GetBinContent(11);
- Double_t twoDirect_5n_5n = fDirect->GetBinContent(12);
-
-*/
-
-
-
-
-
-
-
-/*
-
-
- //CORRELATIONS
- //two particle correlations
- Double_t two_n_n=(AvQ2-AvM)/(AvM*(AvM-1));
- Double_t two_2n_2n=(AvQ_2n2-AvM)/(AvM*(AvM-1));
- Double_t two_3n_3n=(AvQ_3n2-AvM)/(AvM*(AvM-1));
- Double_t two_4n_4n=(AvQ_4n2-AvM)/(AvM*(AvM-1));
- Double_t two_5n_5n=(AvQ_5n2-AvM)/(AvM*(AvM-1));
-
-
-
-
-
-
-
- //four particle correlations:
- //<4>_{n,n|n,n}
- Double_t four_n_n_n_n = (AvQ4+AvQ_2n2-4.*(AvM-2.)*AvQ2-2.*ReQ2nQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))+2./((AvM-1)*(AvM-2));
-
-
- //<4>_{3n|n,n,n}
- //Double_t four_3n_n_n_n=(ReQ3nQnstarQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))-3.*(three_3n_2n_n+three_2n_n_n)/(AvM-3.)-(3.*two+3.*two_2n+two_3n)/((AvM-2.)*(AvM-3.))-1./((AvM-1.)*(AvM-2.)*(AvM-3.)); OK!!!
- //Double_t four_3n_n_n_n_=(ReQ3nQnstarQnstarQnstar-3.*ReQ3nQ2nstarQnstar-3.*ReQ2nQnstarQnstar)/(AvM*(AvM-1.)*(AvM-2.)*(AvM-3.))+ (2.*AvQ_3n2 + 3.*AvQ_2n2 + 6.*AvQ2-6.*AvM ) /(AvM*(AvM-1.)*(AvM-2.)*(AvM-3.)); OK!!! (final version)
-
-
- //<4>_{2n,n|2n,n}
- //Double_t four_2n_n_2n_n = (HereQ2nQnQ2nstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))-2.*(three_2n_n_n+three_3n_2n_n)/(AvM-3.)-((AvM+1.)*two+AvM*two_2n+two_3n)/((AvM-2.)*(AvM-3.))-(AvM)/((AvM-1)*(AvM-2)*(AvM-3)); OK!!!
- //Double_t four_2n_n_2n_n_Alternative = (HereQ2nQnQ2nstarQnstar-2.*ReQ3nQ2nstarQnstar-2.*ReQ2nQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))-((AvM-5)*AvQ2+(AvM-4)*AvQ_2n2-AvQ_3n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))+(AvM-6)/((AvM-1)*(AvM-2)*(AvM-3)); OK!!! (final version)
-
-
- //<6>_{n,n,n|n,n,n}
- //Double_t six_n_n_n_n_n_n=AvQ6/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) - 3.*(2.*five_2n_n_n_n_n)/(AvM-5) - (2.*four_3n_n_n_n+9.*(AvM-2.)*four+9.*four_2n_n_2n_n)/((AvM-4)*(AvM-5)) - 3.*(2.*(3.*AvM-5.)*three_2n_n_n+2.*three_3n_2n_n)/((AvM-3)*(AvM-4)*(AvM-5)) - ((18.*AvM*AvM-45.*AvM+33.)*two+3.*(3.*AvM-4.)*two_2n+two_3n)/((AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) - (6.*AvM*AvM-9.*AvM+4.)/((AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)); OK!!!
- //Double_t six_n_n_n_n_n_n = (AvQ6+9.*HereQ2nQnQ2nstarQnstar-6.*HereQ2nQnQnstarQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) + 4.*(ReQ3nQnstarQnstarQnstar-3.*ReQ3nQ2nstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) + 2.*(9.*(AvM-4.)*ReQ2nQnstarQnstar+2.*AvQ_3n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) - 9.*(AvQ4+AvQ_2n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-5)) + (18.*AvQ2 )/(AvM*(AvM-1)*(AvM-3)*(AvM-4)) - (6.)/((AvM-1)*(AvM-2)*(AvM-3)); OK!!! (final version)
-
-
-
-
-
-*/
-
- /*
-
-
- cout<<"<4>_{3n,n,n,n} correlation from Q-vector = "<<four_3n_n_n_n<<endl;AvQ2
- cout<<"<4>_{3n,n,n,n} correlation directly = "<<fourDirect_3n_n_n_n<<endl;
- cout<<endl;
- cout<<"<4>_{2n,n,2n,n} correlation from Q-vector = "<<four_2n_n_2n_n<<endl;
- cout<<"<4>_{2n,n,2n,n} correlation directly = "<<fourDirect_2n_n_2n_n<<endl;
- cout<<endl;
- cout<<"<5>_{2n,n,n,n,n} correlation from Q-vector = "<<five_2n_n_n_n_n<<endl;
- cout<<"<5>_{2n,n,n,n,n} correlation directly = "<<fiveDirect_2n_n_n_n_n<<endl;
- cout<<endl;
- cout<<"<6>_{n,n,n,n,n,n} correlation from Q-vector = "<<six_n_n_n_n_n_n<<endl;
- cout<<"<6>_{n,n,n,n,n,n} correlation directly = "<<sixDirect_n_n_n_n_n_n<<endl;
- cout<<endl;
- */
-
-
-
-
- /*
-
-//Q-CUMULANTS
-Double_t fourthOrderQCumulant = (AvQ4+AvQ_2n2-2.*ReQ2nQnstarQnstar)/(AvM*(AvM-1.)*(AvM-2.)*(AvM-3.)) - (2.*AvQ2*AvQ2)/(AvM*AvM*(AvM-1.)*(AvM-1.)) - (8.*AvQ2)/(AvM*(AvM-1.)*(AvM-1.)*(AvM-3.)) + 2./((AvM-1.)*(AvM-1.)*(AvM-2.));
-
-Double_t sixthOrderQCumulant = (AvQ6+9.*HereQ2nQnQ2nstarQnstar-6.*HereQ2nQnQnstarQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) + 4.*(ReQ3nQnstarQnstarQnstar - 3.*ReQ3nQ2nstarQnstar+AvQ_3n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) + 18.*(4.*AvM+(AvM-5.)*AvQ2)*ReQ2nQnstarQnstar/(AvM*AvM*(AvM-1)*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-5)) - 9.*(4.*AvM+(AvM-5.)*AvQ2)*(AvQ4+AvQ_2n2)/(AvM*AvM*(AvM-1)*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-5)) + 36.*(5.*AvM-11.)*AvQ2/(AvM*(AvM-1)*(AvM-1)*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)) + 12.*AvQ2*AvQ2*(6.*AvM-3.*AvQ2+AvM*AvQ2)/(AvM*AvM*AvM*(AvM-1)*(AvM-1)*(AvM-1)*(AvM-3)) - 24./((AvM-1)*(AvM-1)*(AvM-1)*(AvM-2)*(AvM-3));
-
-
-
-
-
-
-
-
- //cout<<endl;
- //cout<<"should be the same? "<<threeDirect_3n_2n_n<<" "<<three_2n_n_nTemp<<endl;
- //cout<<endl;
-
-
-
-
-
-
-
- //cout<<endl;
- cout<<"***************************"<<endl;
- cout<<endl;
- cout<<"multiplicity = "<<AvM<<endl;
- cout<<endl;
-
- Double_t two=(AvQ2-AvM)/(AvM*(AvM-1));
- Double_t two_2n=(AvQ_2n2-AvM)/(AvM*(AvM-1));
- Double_t two_3n=(AvQ_3n2-AvM)/(AvM*(AvM-1));
- Double_t two_4n=(AvQ_4n2-AvM)/(AvM*(AvM-1));
- Double_t two_5n=(AvQ_5n2-AvM)/(AvM*(AvM-1));
-
-
- //cout<<"two's = "<<two<<" "<<two_2n<<" "<<two_3n<<" "<<two_4n<<" "<<two_5n<<endl;
- cout<<endl;
- Double_t four=(2.*AvM*(AvM-3.)+AvQ4-4.*(AvM-2.)*AvQ2-2.*ReQ2nQnstarQnstar+AvQ_2n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3));
- Double_t three_2n_n_n=(ReQ2nQnstarQnstar-AvM-2*AvM*(AvM-1)*two-AvM*(AvM-1)*two_2n)/(AvM*(AvM-1)*(AvM-2));
-
- //<3>_{3n,2n,n}
- Double_t three_3n_2n_n=(ReQ3nQ2nstarQnstar-AvM*(AvM-1)*(two+two_2n+two_3n)-AvM)/(AvM*(AvM-1)*(AvM-2));
-
- //<4>_{3n,n,n,n}
- //Double_t four_3n_n_n_n=(ReQ3nQnstarQnstarQnstar-3.*AvM*(AvM-1)*(AvM-2)*(three_2n_n_n+three_3n_2n_n)-AvM*(AvM-1)*(3.*two+3.*two_2n+two_3n)-AvM)/(AvM*(AvM-1)*(AvM-2)*(AvM-3));
-
- Double_t four_3n_n_n_n=(ReQ3nQnstarQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))-3.*(three_3n_2n_n+three_2n_n_n)/(AvM-3.)-(3.*two+3.*two_2n+two_3n)/((AvM-2.)*(AvM-3.))-1./((AvM-1.)*(AvM-2.)*(AvM-3.));
-
-
-
- //<4>_{2n,n,2n,n}
- Double_t four_2n_n_2n_n=(HereQ2nQnQ2nstarQnstar-2.*AvM*(AvM-1)*(AvM-2)*(three_2n_n_n+three_3n_2n_n)-AvM*(AvM-1)*((AvM+1)*two+AvM*two_2n+two_3n)-AvM*AvM)/(AvM*(AvM-1)*(AvM-2)*(AvM-3));
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- cout<<"<2> correlation from Q-vector = "<<two<<endl;
- cout<<"<2> correlation directly = "<<twoDirect_n_n<<endl;
- cout<<endl;
- cout<<"<4> correlation from Q-vector = "<<four<<endl;
- cout<<"<4> correlation directly = "<<fourDirect<<endl;
- cout<<endl;
- cout<<"<3>_{2n,n,n} correlation from Q-vector = "<<three_2n_n_n<<endl;
- cout<<"<3>_{2n,n,n} correlation directly = "<<threeDirect_2n_n_n<<endl;
- cout<<endl;
- cout<<"<3>_{3n,2n,n} correlation from Q-vector = "<<three_3n_2n_n<<endl;
- cout<<"<3>_{3n,2n,n} correlation directly = "<<threeDirect_3n_2n_n<<endl;
- cout<<endl;
- cout<<"<4>_{3n,n,n,n} correlation from Q-vector = "<<four_3n_n_n_n<<endl;
- cout<<"<4>_{3n,n,n,n} correlation directly = "<<fourDirect_3n_n_n_n<<endl;
- cout<<endl;
- cout<<"<4>_{2n,n,2n,n} correlation from Q-vector = "<<four_2n_n_2n_n<<endl;
- cout<<"<4>_{2n,n,2n,n} correlation directly = "<<fourDirect_2n_n_2n_n<<endl;
- cout<<endl;
- cout<<"<5>_{2n,n,n,n,n} correlation from Q-vector = "<<five_2n_n_n_n_n<<endl;
- cout<<"<5>_{2n,n,n,n,n} correlation directly = "<<fiveDirect_2n_n_n_n_n<<endl;
- cout<<endl;
- cout<<"<6>_{n,n,n,n,n,n} correlation from Q-vector = "<<six_n_n_n_n_n_n<<endl;
- cout<<"<6>_{n,n,n,n,n,n} correlation directly = "<<sixDirect_n_n_n_n_n_n<<endl;
- cout<<endl;
-
-
-
-
- //have so far: <2>, <2>_{2n}, <3>_{2n,n,n}, <4>
-
-
- //cout<<"temp re = "<<ReQ2nQnstarQnstar<<endl;
- //cout<<"temp im = "<<ImQ2nQnstarQnstar<<endl;
- //cout<<AvQ2<<" "<<sqrt(AvQ4)<<endl;
- cout<<endl;
-
-
-
-
- if(AvQ2>AvM){
- cout<<"v_2{2} = "<<100*sqrt(two)<<"%"<<endl;
- }else{
- cout<<"v_2{2} = Im"<<endl;
- }
-
- if(four>0){
- cout<<"v_2{4} = "<<100*pow(four,1./4.)<<"%"<<endl;
- }else{
- cout<<"v_2{4} = Im"<<endl;
- }
-
-
-
-
- cout<<"AvM = "<<AvM<<endl;
- cout<<endl;
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<"flow estimates from Q-cumulants:"<<endl;
- Double_t cumulant2Q=two;
- Double_t cumulant4Q=four-2.*two*two;
- Double_t cumulant6Q=12.*two*two*two-9.*two*four+six_n_n_n_n_n_n;
-
- //cout<<endl;
- if(cumulant2Q>0.){
- cout<<"v_2{2} = "<<100*pow(cumulant2Q,1./2.)<<"%"<<endl;
- }else{
- cout<<"v_2{2} = Im"<<endl;
- }
- if(cumulant4Q<0.){
- cout<<"v_2{4} = "<<100*pow(-cumulant4Q,1./4.)<<"%"<<endl;
- }else{
- cout<<"v_2{4} = Im"<<endl;
- }
- if(cumulant6Q>0.){
- cout<<"v_2{6} = "<<100*pow((1./4.)*cumulant6Q,1./6.)<<"%"<<endl;
- }else{
- cout<<"v_2{6} = Im"<<endl;
- }
-
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<endl;
-
-
- cout<<endl;
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<"flow estimates from Q-cumulants"<<endl;
- cout<<"without multiplicities fluctuations:"<<endl;
- Double_t cumulant2QnoFluct=fQCorr->GetBinContent(1);
- Double_t cumulant4QnoFluct=fQCorr->GetBinContent(11)-2.*fQCorr->GetBinContent(1)*fQCorr->GetBinContent(1);
- Double_t cumulant6QnoFluct=12.*fQCorr->GetBinContent(1)*fQCorr->GetBinContent(1)*fQCorr->GetBinContent(1) - 9.*fQCorr->GetBinContent(1)*fQCorr->GetBinContent(11)+fQCorr->GetBinContent(21);
-
- //cout<<endl;
- if(cumulant2QnoFluct>0.){
- cout<<"v_2{2} = "<<100*pow(cumulant2QnoFluct,1./2.)<<"%"<<endl;
- fIntRes->SetBinContent(1,100*pow(cumulant2QnoFluct,1./2.));
- }else{
- cout<<"v_2{2} = Im"<<endl;
- }
- if(cumulant4QnoFluct<0.){
- cout<<"v_2{4} = "<<100*pow(-cumulant4QnoFluct,1./4.)<<"%"<<endl;
- fIntRes->SetBinContent(2,100*pow(-cumulant4QnoFluct,1./4.));
- }else{
- cout<<"v_2{4} = Im"<<endl;
- }
- if(cumulant6QnoFluct>0.){
- cout<<"v_2{6} = "<<100*pow((1./4.)*cumulant6QnoFluct,1./6.)<<"%"<<endl;
- fIntRes->SetBinContent(3,100*pow((1./4.)*cumulant6QnoFluct,1./6.));
- }else{
- cout<<"v_2{6} = Im"<<endl;
- }
-
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<endl;
-
-
-
-
-
-
- cout<<"*****************************************************************"<<endl;
- cout<<"*****************************************************************"<<endl;
- cout<<"direct (nested loops) correlations vs correlations from Q-vectors"<<endl;
- cout<<"for ARBITRARY multiplicity:"<<endl;
- cout<<endl;
- cout<<"avarage multiplicity = "<<AvM<<endl;
- cout<<endl;
- cout<<"<2>_{n|n} correlation from Q-vector = "<<fQCorr->GetBinContent(1)<<endl;
- cout<<"<2>_{n|n} correlation nested loops = "<<fDirect->GetBinContent(1)<<endl;
- cout<<endl;
- cout<<"<4>_{n,n|n,n} correlation from Q-vector = "<<fQCorr->GetBinContent(11)<<endl;
- cout<<"<4>_{n,n|n,n} correlation nested loops = "<<fDirect->GetBinContent(2)<<endl;
- cout<<endl;
- cout<<"<6>_{n,n,n|n,n,n} correlation from Q-vector = "<<fQCorr->GetBinContent(21)<<endl;
- cout<<"<6>_{n,n,n|n,n,n} correlation nested loops = "<<fDirect->GetBinContent(8)<<endl;
- cout<<"*****************************************************************"<<endl;
- cout<<"*****************************************************************"<<endl;
-
-
-
-
-
-
-
- cout<<"****** BDO ******"<<endl;
- if(BDO4<0.){
- cout<<"v_2{4} = "<<100*pow(-BDO4,1./4.)<<"%"<<endl;
- }else{
- cout<<"v_2{4} = Im"<<endl;
- }
- if(BDO6>0.){
- cout<<"v_2{6} = "<<100*pow((1./4.)*BDO6,1./6.)<<"%"<<endl;
- }else{
- cout<<"v_2{6} = Im"<<endl;
- }
-
-
-
- cout<<endl;
- cout<<"fourth order Q-cumulant should be the same and BDO4? "<<fourthOrderQCumulant<<" "<<cumulant4Q<<" "<<BDO4<<endl;
- cout<<endl;
- cout<<endl;
- cout<<"sixth order Q-cumulant should be the same and BDO6? "<<sixthOrderQCumulant<<" "<<cumulant6Q<<" "<<BDO6<<endl;
- cout<<endl;
-
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<endl;
-
-
-
- cout<<"*********************"<<endl;
- cout<<"multiplicity = "<<AvM<<endl;
- cout<<"correlations: "<<endl;
- cout<<"direct = "<<twoDirect_n_n<<endl;
- cout<<"approx = "<<two_n_n<<endl;
- cout<<"weight = "<<fQCorr->GetBinContent(1)<<endl;
- cout<<"*********************"<<endl;
- cout<<"test = "<<fQCorr->GetBinContent(11)<<" "<<four<<endl;
-
-
-
-
-
-
-
-
-
-
-
-
-*/
-
-
-
-
-
-
-
-
-
-
//type of analysis can be: ESD, AOD, MC, ESDMC0, ESDMC1
//const TString type = "ESD";
-void CompareFlowResults(TString type="ESD")
+void compareFlowResults(TString type="ESD")
{
//load needed libraries:
gSystem->AddIncludePath("-I$ROOTSYS/include");
// set here which plots will be shown by default
//==================================================================================
//RP = particles used to determine the reaction plane
- Bool_t plotIntFlowRP = kFALSE; //integrated flow RP
- Bool_t plotDiffFlowPtRP = kFALSE; //differential flow (Pt,RP)
- Bool_t plotDiffFlowEtaRP = kFALSE; //differential flow (Eta,RP)
+ Bool_t plotIntFlowRP = kTRUE; //integrated flow RP
+ Bool_t plotDiffFlowPtRP = kTRUE; //differential flow (Pt,RP)
+ Bool_t plotDiffFlowEtaRP = kTRUE; //differential flow (Eta,RP)
//POI = particle of interest
Bool_t plotIntFlowPOI = kTRUE; //integrated flow POI
Bool_t plotDiffFlowPtPOI = kTRUE; //differential flow (Pt,POI)
//==================================================================================
// accessing output files
//==================================================================================
-
//open the output files:
TString inputFileNameMCEP = "outputMCEPanalysis";
TFile* fileMCEP = NULL;
Int_t meshStyle = 1001;
Int_t meshColor = kRed-10;
- //marker style and color (int. flow)
+ //marker style and color (plots for int. flow)
Int_t markerStyle = 21;
Int_t markerColor = kRed-3;
//==================================================================================
//booking the graph to store flow values and errors from all methods:
Double_t x[nMethods] = {0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5};
Double_t xError[nMethods] = {0.};
- Double_t flowValue[nMethods] = {0.};//to be removed
- Double_t flowError[nMethods] = {0.};//to be removed
+ Double_t flowValue[nMethods] = {0.};
+ Double_t flowError[nMethods] = {0.};
Double_t flowValueRP[nMethods] = {0.};
Double_t flowErrorRP[nMethods] = {0.};
Double_t flowValuePOI[nMethods] = {0.};
mcepCommonHist = dynamic_cast<AliFlowCommonHist*> (pListMCEP->FindObject("AliFlowCommonHistMCEP"));
mcepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListMCEP->FindObject("AliFlowCommonHistResultsMCEP"));
if(mcepCommonHistRes) {
- flowValue[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
- flowError[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinContent(1);
+ flowError[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinError(1);
flowValueRP[0] = (mcepCommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
flowErrorRP[0] = (mcepCommonHistRes->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[0] = (mcepCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
lyz1CommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZ1->FindObject("AliFlowCommonHistLYZ1"));
lyz1CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZ1->FindObject("AliFlowCommonHistResultsLYZ1"));
if(lyz1CommonHistRes) {
- flowValue[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
- flowError[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinContent(1);
+ flowError[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinError(1);
flowValueRP[10] = (lyz1CommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
flowErrorRP[10] = (lyz1CommonHistRes->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[10] = (lyz1CommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
lyzepCommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZEP->FindObject("AliFlowCommonHistLYZEP"));
lyzepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZEP->FindObject("AliFlowCommonHistResultsLYZEP"));
if(lyzepCommonHistRes) {
- flowValue[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
- //flowError[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinContent(1);
+ //flowError[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinError(1);
flowValueRP[11] = (lyzepCommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
//flowErrorRP[11] = (lyzepCommonHistRes->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[11] = (lyzepCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
fqdCommonHist = dynamic_cast<AliFlowCommonHist*> (pListFQD->FindObject("AliFlowCommonHistFQD"));
fqdCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListFQD->FindObject("AliFlowCommonHistResultsFQD"));
if(fqdCommonHistRes) {
- flowValue[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
- flowError[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinContent(1);
+ flowError[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinError(1);
flowValueRP[9] = (fqdCommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
flowErrorRP[9] = (fqdCommonHistRes->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[9] = (fqdCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
gfcCommonHist = dynamic_cast<AliFlowCommonHist*> (pListGFC->FindObject("AliFlowCommonHistGFC"));
gfcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults2ndOrderGFC"));
if(gfcCommonHistRes2) {
- flowValue[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);//to be removed
- flowError[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);
+ flowError[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinError(1);
flowValueRP[1] = (gfcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1);
flowErrorRP[1] = (gfcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[1] = (gfcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1);
}
gfcCommonHistRes4 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults4thOrderGFC"));
if(gfcCommonHistRes4) {
- flowValue[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);//to be removed
- flowError[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);
+ flowError[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinError(1);
flowValueRP[3] = (gfcCommonHistRes4->GetHistIntFlowRP())->GetBinContent(1);
flowErrorRP[3] = (gfcCommonHistRes4->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[3] = (gfcCommonHistRes4->GetHistIntFlowPOI())->GetBinContent(1);
}
gfcCommonHistRes6 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults6thOrderGFC"));
if(gfcCommonHistRes6) {
- flowValue[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);//to be removed
- flowError[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);
+ flowError[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinError(1);
flowValueRP[5] = (gfcCommonHistRes6->GetHistIntFlowRP())->GetBinContent(1);
flowErrorRP[5] = (gfcCommonHistRes6->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[5] = (gfcCommonHistRes6->GetHistIntFlowPOI())->GetBinContent(1);
}
gfcCommonHistRes8 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults8thOrderGFC"));
if(gfcCommonHistRes8) {
- flowValue[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);//to be removed
- flowError[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);
+ flowError[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinError(1);
flowValueRP[7] = (gfcCommonHistRes8->GetHistIntFlowRP())->GetBinContent(1);
flowErrorRP[7] = (gfcCommonHistRes8->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[7] = (gfcCommonHistRes8->GetHistIntFlowPOI())->GetBinContent(1);
qcCommonHist2 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist2ndOrderQC"));
qcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults2ndOrderQC"));
if(qcCommonHistRes2) {
- flowValue[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);//to be removed
- //flowError[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);
+ //flowError[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinError(1);
flowValueRP[2] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1);
//flowErrorRP[2] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[2] = (qcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1);
qcCommonHist4 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist4thOrderQC"));
qcCommonHistRes4 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults4thOrderQC"));
if(qcCommonHistRes4) {
- flowValue[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);//to be removed
- //flowError[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);
+ //flowError[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinError(1);
flowValueRP[4] = (qcCommonHistRes4->GetHistIntFlowRP())->GetBinContent(1);
//flowErrorRP[4] = (qcCommonHistRes4->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[4] = (qcCommonHistRes4->GetHistIntFlowPOI())->GetBinContent(1);
qcCommonHist6 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist6thOrderQC"));
qcCommonHistRes6 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults6thOrderQC"));
if(qcCommonHistRes6) {
- flowValue[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);//to be removed
- //flowError[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);
+ //flowError[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinError(1);
flowValueRP[6] = (qcCommonHistRes6->GetHistIntFlowRP())->GetBinContent(1);
//flowErrorRP[6] = (qcCommonHistRes6->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[6] = (qcCommonHistRes6->GetHistIntFlowPOI())->GetBinContent(1);
qcCommonHist8 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist8thOrderQC"));
qcCommonHistRes8 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults8thOrderQC"));
if(qcCommonHistRes8) {
- flowValue[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);//to be removed
- //flowError[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinError(1);//to be removed
+ flowValue[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);
+ //flowError[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinError(1);
flowValueRP[8] = (qcCommonHistRes8->GetHistIntFlowRP())->GetBinContent(1);
//flowErrorRP[8] = (qcCommonHistRes8->GetHistIntFlowRP())->GetBinError(1);
flowValuePOI[8] = (qcCommonHistRes8->GetHistIntFlowPOI())->GetBinContent(1);
}
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+ //no-name:
Double_t dMax=flowValue[0]+flowError[0];
Double_t dMin=flowValue[0]-flowError[0];
if(dMin>flowValue[i]-flowError[i]) dMin=flowValue[i]-flowError[i];
}
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
//RP:
Double_t dMaxRP=flowValueRP[0]+flowErrorRP[0];
}
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+ //no-name:
TGraph* flowResults = new TGraphErrors(nMethods, x, flowValue, xError, flowError);
-
flowResults->SetMarkerStyle(markerStyle);
flowResults->SetMarkerColor(markerColor);
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
//RP:
- TGraph* flowResultsRP = new TGraphErrors(nMethods, x, flowValueRP, xError, flowErrorRP);
-
+ TGraph* flowResultsRP = new TGraphErrors(nMethods, x, flowValueRP, xError, flowErrorRP);
flowResultsRP->SetMarkerStyle(markerStyle);
flowResultsRP->SetMarkerColor(markerColor);
//POI:
TGraph* flowResultsPOI = new TGraphErrors(nMethods, x, flowValuePOI, xError, flowErrorPOI);
-
flowResultsPOI->SetMarkerStyle(markerStyle);
flowResultsPOI->SetMarkerColor(markerColor);
//----------------------------------------------------------------------------------
//cosmetics: mesh for MC error bands (integrated flow)
- TGraph* pMesh = NULL;//to be removed
+ TGraph* pMesh = NULL;
TGraph* pMeshRP = NULL;
TGraph* pMeshPOI = NULL;
if(intFlowAll && mcepCommonHistRes) {
- //Double_t valueMC = intFlowAll->GetBinContent(1);
- //Double_t errorMC = intFlowAll->GetBinError(1);
- //Int_t nPts = intFlowAll->GetNbinsX();
-
Int_t nPts = nMethods;
- Double_t valueMC = flowValue[0];//to be removed
- Double_t errorMC = flowError[0];//to be removed
+ Double_t valueMC = flowValue[0];
+ Double_t errorMC = flowError[0];
Double_t valueMCRP = flowValueRP[0];
Double_t errorMCRP = flowErrorRP[0];
Double_t valueMCPOI = flowValuePOI[0];
Double_t errorMCPOI = flowErrorPOI[0];
- pMesh = new TGraph(nPts);//to be removed
+ pMesh = new TGraph(nPts);
pMeshRP = new TGraph(nPts);
pMeshPOI = new TGraph(nPts);
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+ //no-name:
pMesh->SetPoint(1,0,valueMC+errorMC);
pMesh->SetPoint(2,nPts+1,valueMC+errorMC);
pMesh->SetPoint(3,nPts+1,valueMC-errorMC);
pMesh->SetPoint(4,0,valueMC-errorMC);
- pMesh->SetPoint(5,0,valueMC+errorMC);
-
+ pMesh->SetPoint(5,0,valueMC+errorMC);
pMesh->SetFillStyle(meshStyle);
pMesh->SetFillColor(meshColor);
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-
+
//RP:
pMeshRP->SetPoint(1,0,valueMCRP+errorMCRP);
pMeshRP->SetPoint(2,nPts+1,valueMCRP+errorMCRP);
pMeshRP->SetPoint(3,nPts+1,valueMCRP-errorMCRP);
pMeshRP->SetPoint(4,0,valueMCRP-errorMCRP);
- pMeshRP->SetPoint(5,0,valueMCRP+errorMCRP);
-
+ pMeshRP->SetPoint(5,0,valueMCRP+errorMCRP);
pMeshRP->SetFillStyle(meshStyle);
pMeshRP->SetFillColor(meshColor);
pMeshPOI->SetPoint(2,nPts+1,valueMCPOI+errorMCPOI);
pMeshPOI->SetPoint(3,nPts+1,valueMCPOI-errorMCPOI);
pMeshPOI->SetPoint(4,0,valueMCPOI-errorMCPOI);
- pMeshPOI->SetPoint(5,0,valueMCPOI+errorMCPOI);
-
+ pMeshPOI->SetPoint(5,0,valueMCPOI+errorMCPOI);
pMeshPOI->SetFillStyle(meshStyle);
pMeshPOI->SetFillColor(meshColor);
}
TPaveText *textDefault = new TPaveText(0.05,0.77,0.95,0.90,"NDC");
textDefault->SetTextFont(72);
textDefault->SetTextSize(0.08);
- //textDefault->SetLineColor(kFALSE);
- //textDefault->SetShadowColor(kFALSE);
TString *entryDefaultAvM = new TString("Average Multiplicity");
TString *entryDefaultAnd = new TString("and");
textDefault->AddText(entryDefaultAnd->Data());
textDefault->AddText(entryDefaultNumOfEvts->Data());
- //results:
- TPaveText *textResults = new TPaveText(0.05,0.12,0.95,0.70,"NDC");//to be removed
- textResults->SetTextFont(72);//to be removed
- textResults->SetTextSize(0.06);//to be removed
- //textResults->SetLineColor(kFALSE);
- //textResults->SetShadowColor(kFALSE);
+ //results (no-name):
+ TPaveText *textResults = new TPaveText(0.05,0.12,0.95,0.70,"NDC");
+ textResults->SetTextFont(72);
+ textResults->SetTextSize(0.06);
//results (RP):
TPaveText *textResultsRP = new TPaveText(0.05,0.12,0.95,0.70,"NDC");
textResultsPOI->SetTextFont(72);
textResultsPOI->SetTextSize(0.06);
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!!!
+ //no-name:
TString *entryMC = new TString("MC ...... ");
TString *entryGFC = new TString("GFC ..... ");
TString *entryQC2 = new TString("QC{2} ... ");
TString *entryFQD = new TString("FQD ..... ");
TString *entryLYZ1 = new TString("LYZ ..... ");
TString *entryLYZEP = new TString("LYZEP ... ");
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
//RP:
TString *entryMCRP = new TString("MC ...... ");
TString *entryLYZ1POI = new TString("LYZ ..... ");
TString *entryLYZEPPOI = new TString("LYZEP ... ");
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+ //no-name:
Double_t avMultMC=0.;
Long_t nEvtsMC=0;
-
Double_t avMultGFC=0.;
Long_t nEvtsGFC=0;
-
Double_t avMultQC2=0., avMultQC4=0., avMultQC6=0., avMultQC8=0.;
Long_t nEvtsQC2=0, nEvtsQC4=0, nEvtsQC6=0, nEvtsQC8=0;
-
Double_t avMultFQD=0.;
Long_t nEvtsFQD=0;
-
Double_t avMultLYZ1=0.;
Long_t nEvtsLYZ1=0;
-
Double_t avMultLYZEP=0.;
Long_t nEvtsLYZEP=0;
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
//RP:
Double_t avMultMCRP=0.;
Long_t nEvtsMCRP=0;
-
Double_t avMultGFCRP=0.;
Long_t nEvtsGFCRP=0;
-
Double_t avMultQC2RP=0., avMultQC4RP=0., avMultQC6RP=0., avMultQC8RP=0.;
Long_t nEvtsQC2RP=0, nEvtsQC4RP=0, nEvtsQC6RP=0, nEvtsQC8RP=0;
-
Double_t avMultFQDRP=0.;
Long_t nEvtsFQDRP=0;
-
Double_t avMultLYZ1RP=0.;
Long_t nEvtsLYZ1RP=0;
-
Double_t avMultLYZEPRP=0.;
Long_t nEvtsLYZEPRP=0;
//POI:
Double_t avMultMCPOI=0.;
Long_t nEvtsMCPOI=0;
-
Double_t avMultGFCPOI=0.;
Long_t nEvtsGFCPOI=0;
-
Double_t avMultQC2POI=0., avMultQC4POI=0., avMultQC6POI=0., avMultQC8POI=0.;
Long_t nEvtsQC2POI=0, nEvtsQC4POI=0, nEvtsQC6POI=0, nEvtsQC8POI=0;
-
Double_t avMultFQDPOI=0.;
Long_t nEvtsFQDPOI=0;
-
Double_t avMultLYZ1POI=0.;
Long_t nEvtsLYZ1POI=0;
-
Double_t avMultLYZEPPOI=0.;
Long_t nEvtsLYZEPPOI=0;
//MC:
if(mcepCommonHist) {
- avMultMC = (mcepCommonHist->GetHistMultInt())->GetMean();//to be removed
- nEvtsMC = (mcepCommonHist->GetHistMultInt())->GetEntries();//to be removed
+ avMultMC = (mcepCommonHist->GetHistMultInt())->GetMean();
+ nEvtsMC = (mcepCommonHist->GetHistMultInt())->GetEntries();
avMultMCRP = (mcepCommonHist->GetHistMultInt())->GetMean();
nEvtsMCRP = (mcepCommonHist->GetHistMultInt())->GetEntries();
avMultMCPOI = (mcepCommonHist->GetHistMultDiff())->GetMean();
nEvtsMCPOI = (mcepCommonHist->GetHistMultDiff())->GetEntries();
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
if(entryMC) {
entryMC->Append("M = ");
(*entryMC)+=(Long_t)avMultMC;
entryMC->Append(", N = ");
(*entryMC)+=(Long_t)nEvtsMC;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
if(entryMCRP) {
entryMCRP->Append("M = ");
(*entryMCRP)+=(Long_t)avMultMCRP;
entryMCRP->Append(", N = ");
(*entryMCRP)+=(Long_t)nEvtsMCRP;
}
- if(entryMCPOI) {
+
+ if(entryMCPOI) {
entryMCPOI->Append("M = ");
(*entryMCPOI)+=(Long_t)avMultMCPOI;
entryMCPOI->Append(", N = ");
(*entryMCPOI)+=(Long_t)nEvtsMCPOI;
- }
+ }
//GFC:
if(gfcCommonHist) {
- avMultGFC = (gfcCommonHist->GetHistMultInt())->GetMean();//to be removed
- nEvtsGFC = (gfcCommonHist->GetHistMultInt())->GetEntries();//to be removed
+ avMultGFC = (gfcCommonHist->GetHistMultInt())->GetMean();
+ nEvtsGFC = (gfcCommonHist->GetHistMultInt())->GetEntries();
avMultGFCRP = (gfcCommonHist->GetHistMultInt())->GetMean();
nEvtsGFCRP = (gfcCommonHist->GetHistMultInt())->GetEntries();
avMultGFCPOI = (gfcCommonHist->GetHistMultDiff())->GetMean();
nEvtsGFCPOI = (gfcCommonHist->GetHistMultDiff())->GetEntries();
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
if(entryGFC) {
entryGFC->Append("M = ");
(*entryGFC)+=(Long_t)avMultGFC;
entryGFC->Append(", N = ");
(*entryGFC)+=(Long_t)nEvtsGFC;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
if(entryGFCRP) {
entryGFCRP->Append("M = ");
(*entryGFCRP)+=(Long_t)avMultGFCRP;
//QC:
if(qcCommonHist2) {
- avMultQC2 = (qcCommonHist2->GetHistMultInt())->GetMean();//to be removed
- nEvtsQC2 = (qcCommonHist2->GetHistMultInt())->GetEntries();//to be removed
+ avMultQC2 = (qcCommonHist2->GetHistMultInt())->GetMean();
+ nEvtsQC2 = (qcCommonHist2->GetHistMultInt())->GetEntries();
avMultQC2RP = (qcCommonHist2->GetHistMultInt())->GetMean();
nEvtsQC2RP = (qcCommonHist2->GetHistMultInt())->GetEntries();
avMultQC2POI = (qcCommonHist2->GetHistMultDiff())->GetMean();
nEvtsQC2POI = (qcCommonHist2->GetHistMultDiff())->GetEntries();
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
if(entryQC2)
{
entryQC2->Append("M = ");
entryQC2->Append(", N = ");
(*entryQC2)+=(Long_t)nEvtsQC2;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
if(entryQC2RP)
{
entryQC2RP->Append("M = ");
entryQC2RP->Append(", N = ");
(*entryQC2RP)+=(Long_t)nEvtsQC2RP;
}
+
if(entryQC2POI)
{
entryQC2POI->Append("M = ");
if(qcCommonHist4)
{
- avMultQC4 = (qcCommonHist4->GetHistMultInt())->GetMean();//to be removed
- nEvtsQC4 = (qcCommonHist4->GetHistMultInt())->GetEntries();//to be removed
+ avMultQC4 = (qcCommonHist4->GetHistMultInt())->GetMean();
+ nEvtsQC4 = (qcCommonHist4->GetHistMultInt())->GetEntries();
avMultQC4RP = (qcCommonHist4->GetHistMultInt())->GetMean();
nEvtsQC4RP = (qcCommonHist4->GetHistMultInt())->GetEntries();
avMultQC4POI = (qcCommonHist4->GetHistMultDiff())->GetMean();
nEvtsQC4POI = (qcCommonHist4->GetHistMultDiff())->GetEntries();
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
if(entryQC4)
{
entryQC4->Append("M = ");
entryQC4->Append(", N = ");
(*entryQC4)+=(Long_t)nEvtsQC4;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
if(entryQC4RP)
{
entryQC4RP->Append("M = ");
entryQC4RP->Append(", N = ");
(*entryQC4RP)+=(Long_t)nEvtsQC4RP;
}
+
if(entryQC4POI)
{
entryQC4POI->Append("M = ");
if(qcCommonHist6)
{
- avMultQC6 = (qcCommonHist6->GetHistMultInt())->GetMean();//to be removed
- nEvtsQC6 = (qcCommonHist6->GetHistMultInt())->GetEntries();//to be removed
+ avMultQC6 = (qcCommonHist6->GetHistMultInt())->GetMean();
+ nEvtsQC6 = (qcCommonHist6->GetHistMultInt())->GetEntries();
avMultQC6RP = (qcCommonHist6->GetHistMultInt())->GetMean();
nEvtsQC6RP = (qcCommonHist6->GetHistMultInt())->GetEntries();
avMultQC6POI = (qcCommonHist6->GetHistMultDiff())->GetMean();
nEvtsQC6POI = (qcCommonHist6->GetHistMultDiff())->GetEntries();
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
if(entryQC6)
{
entryQC6->Append("M = ");
entryQC6->Append(", N = ");
(*entryQC6)+=(Long_t)nEvtsQC6;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
if(entryQC6RP)
{
entryQC6RP->Append("M = ");
entryQC6RP->Append(", N = ");
(*entryQC6RP)+=(Long_t)nEvtsQC6RP;
}
+
if(entryQC6POI)
{
entryQC6POI->Append("M = ");
if(qcCommonHist8)
{
- avMultQC8 = (qcCommonHist8->GetHistMultInt())->GetMean();//to be removed
- nEvtsQC8 = (qcCommonHist8->GetHistMultInt())->GetEntries();//to be removed
+ avMultQC8 = (qcCommonHist8->GetHistMultInt())->GetMean();
+ nEvtsQC8 = (qcCommonHist8->GetHistMultInt())->GetEntries();
avMultQC8RP = (qcCommonHist8->GetHistMultInt())->GetMean();
nEvtsQC8RP = (qcCommonHist8->GetHistMultInt())->GetEntries();
avMultQC8POI = (qcCommonHist8->GetHistMultDiff())->GetMean();
nEvtsQC8POI = (qcCommonHist8->GetHistMultDiff())->GetEntries();
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
if(entryQC8)
{
entryQC8->Append("M = ");
entryQC8->Append(", N = ");
(*entryQC8)+=(Long_t)nEvtsQC8;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
if(entryQC8RP)
{
entryQC8RP->Append("M = ");
entryQC8RP->Append(", N = ");
(*entryQC8RP)+=(Long_t)nEvtsQC8RP;
}
+
if(entryQC8POI)
{
entryQC8POI->Append("M = ");
//FQD:
if(fqdCommonHist)
{
- avMultFQD = (fqdCommonHist->GetHistMultInt())->GetMean();//to be removed
- nEvtsFQD = (fqdCommonHist->GetHistMultInt())->GetEntries();//to be removed
+ avMultFQD = (fqdCommonHist->GetHistMultInt())->GetMean();
+ nEvtsFQD = (fqdCommonHist->GetHistMultInt())->GetEntries();
avMultFQDRP = (fqdCommonHist->GetHistMultInt())->GetMean();
nEvtsFQDRP = (fqdCommonHist->GetHistMultInt())->GetEntries();
avMultFQDPOI = (fqdCommonHist->GetHistMultDiff())->GetMean();
nEvtsFQDPOI = (fqdCommonHist->GetHistMultDiff())->GetEntries();
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
if(entryFQD)
{
entryFQD->Append("M = ");
entryFQD->Append(", N = ");
(*entryFQD)+=(Long_t)nEvtsFQD;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
if(entryFQDRP)
{
entryFQDRP->Append("M = ");
entryFQDRP->Append(", N = ");
(*entryFQDRP)+=(Long_t)nEvtsFQDRP;
}
+
if(entryFQDPOI)
{
entryFQDPOI->Append("M = ");
//LYZ1:
if(lyz1CommonHist)
{
- avMultLYZ1 = (lyz1CommonHist->GetHistMultInt())->GetMean();//to be removed
- nEvtsLYZ1 = (lyz1CommonHist->GetHistMultInt())->GetEntries();//to be removed
+ avMultLYZ1 = (lyz1CommonHist->GetHistMultInt())->GetMean();
+ nEvtsLYZ1 = (lyz1CommonHist->GetHistMultInt())->GetEntries();
avMultLYZ1RP = (lyz1CommonHist->GetHistMultInt())->GetMean();
nEvtsLYZ1RP = (lyz1CommonHist->GetHistMultInt())->GetEntries();
avMultLYZ1POI = (lyz1CommonHist->GetHistMultDiff())->GetMean();
nEvtsLYZ1POI = (lyz1CommonHist->GetHistMultDiff())->GetEntries();
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
if(entryLYZ1)
{
entryLYZ1->Append("M = ");
entryLYZ1->Append(", N = ");
(*entryLYZ1)+=(Long_t)nEvtsLYZ1;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
if(entryLYZ1RP)
{
entryLYZ1RP->Append("M = ");
entryLYZ1RP->Append(", N = ");
(*entryLYZ1RP)+=(Long_t)nEvtsLYZ1RP;
}
+
if(entryLYZ1POI)
{
entryLYZ1POI->Append("M = ");
//LYZEP:
if(lyzepCommonHist)
{
- avMultLYZEP = (lyzepCommonHist->GetHistMultInt())->GetMean();//to be removed
- nEvtsLYZEP = (lyzepCommonHist->GetHistMultInt())->GetEntries();//to be removed
+ avMultLYZEP = (lyzepCommonHist->GetHistMultInt())->GetMean();
+ nEvtsLYZEP = (lyzepCommonHist->GetHistMultInt())->GetEntries();
avMultLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetMean();
nEvtsLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetEntries();
avMultLYZEPPOI = (lyzepCommonHist->GetHistMultDiff())->GetMean();
nEvtsLYZEPPOI = (lyzepCommonHist->GetHistMultDiff())->GetEntries();
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
if(entryLYZEP)
{
entryLYZEP->Append("M = ");
entryLYZEP->Append(", N = ");
(*entryLYZEP)+=(Long_t)nEvtsLYZEP;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
if(entryLYZEPRP)
{
entryLYZEPRP->Append("M = ");
entryLYZEPRP->Append(", N = ");
(*entryLYZEPRP)+=(Long_t)nEvtsLYZEPRP;
}
+
if(entryLYZEPPOI)
{
entryLYZEPPOI->Append("M = ");
(*entryLYZEPPOI)+=(Long_t)nEvtsLYZEPPOI;
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+ //no-name:
if(textResults)
{
textResults->AddText(entryMC->Data());
textResults->AddText(entryLYZ1->Data());
textResults->AddText(entryLYZEP->Data());
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
//RP:
if(textResultsRP)
{
}
//----------------------------------------------------------------------------------
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
+
+
+
//----------------------------------------------------------------------------------
- //final drawing for integrated flow:
+ //final drawing for integrated flow (no-name):
TCanvas* intFlowAllCanvas = new TCanvas("Integrated Flow","Integrated Flow",1000,600);
intFlowAllCanvas->Divide(2,1);
textDefault->Draw();
textResults->Draw();
}
-
//----------------------------------------------------------------------------------
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
}//end of if(plotIntFlowRP}
//----------------------------------------------------------------------------------
+
+
+
//----------------------------------------------------------------------------------
//final drawing for integrated flow of POI (i.e. of particles of interest):
if(plotIntFlowPOI)
//----------------------------------------------------------------------------------
-
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
- //----------------------------------------------------------------------------------
- //cosmetics: Monte Carlo error bands for differential flow
- TGraph* pMeshDiffFlow = NULL;
- if(mcepCommonHistRes)
- {
- Int_t nPtsDiffFlow = (mcepCommonHistRes->GetHistDiffFlow())->GetNbinsX();
- Double_t binWidth = (mcepCommonHistRes->GetHistDiffFlow())->GetBinWidth(1);//assuming that all bins have the same width
-
- pMeshDiffFlow = new TGraph(2*nPtsDiffFlow+1);
-
- Double_t valueMC=0., errorMC=0.;
- for(Int_t i=1;i<nPtsDiffFlow+1;i++)
- {
- valueMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinContent(i);
- errorMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinError(i);
- pMeshDiffFlow->SetPoint(i,(i-0.5)*binWidth,valueMC+errorMC);
- }
- for(Int_t i=nPtsDiffFlow+1;i<2*nPtsDiffFlow+1;i++)
- {
- valueMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinContent(2*nPtsDiffFlow+1-i);
- errorMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinError(2*nPtsDiffFlow+1-i);
- pMeshDiffFlow->SetPoint(i,(2*nPtsDiffFlow-i+0.5)*binWidth,valueMC-errorMC);
- }
- pMeshDiffFlow->SetPoint(2*nPtsDiffFlow+1,0.5*binWidth,valueMC+errorMC);
- pMeshDiffFlow->SetFillStyle(meshStyle);
- pMeshDiffFlow->SetFillColor(meshColor);
- }
- //----------------------------------------------------------------------------------
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
//----------------------------------------------------------------------------------
//RP:
//----------------------------------------------------------------------------------
//MCEP = Monte Carlo Event Plane
- Double_t avMultDiffFlowMC=0.;//to be removed
- Double_t nEvtsDiffFlowMC=0;//to be removed
Double_t avMultDiffFlowMCRP=0.;
Double_t nEvtsDiffFlowMCRP=0;
Double_t avMultDiffFlowMCPOI=0.;
{
if(mcepCommonHistRes)
{
- (mcepCommonHistRes->GetHistDiffFlow())->SetMarkerColor(markerColorMC);//to be removed
- (mcepCommonHistRes->GetHistDiffFlow())->SetMarkerStyle(markerStyleMC);//to be removed
(mcepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorMC);
(mcepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleMC);
(mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorMC);
}
if(mcepCommonHist)
{
- avMultDiffFlowMC = (mcepCommonHist->GetHistMultDiff())->GetMean();//to be removed
- nEvtsDiffFlowMC = (mcepCommonHist->GetHistMultDiff())->GetEntries();//to be removed
avMultDiffFlowMCRP = (mcepCommonHist->GetHistMultInt())->GetMean();
nEvtsDiffFlowMCRP = (mcepCommonHist->GetHistMultInt())->GetEntries();
avMultDiffFlowMCPOI = (mcepCommonHist->GetHistMultDiff())->GetMean();
{
if(gfcCommonHistRes2)
{
- (gfcCommonHistRes2->GetHistDiffFlow())->SetMarkerColor(markerColorGFC2);//to be removed
- (gfcCommonHistRes2->GetHistDiffFlow())->SetMarkerStyle(markerStyleGFC2);//to be removed
(gfcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC2);
(gfcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC2);
(gfcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC2);
}
if(gfcCommonHistRes4)
{
- (gfcCommonHistRes4->GetHistDiffFlow())->SetMarkerColor(markerColorGFC4);//to be removed
- (gfcCommonHistRes4->GetHistDiffFlow())->SetMarkerStyle(markerStyleGFC4);//to be removed
(gfcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC4);
(gfcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC4);
(gfcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC4);
}
if(gfcCommonHistRes6)
{
- (gfcCommonHistRes6->GetHistDiffFlow())->SetMarkerColor(markerColorGFC6);//to be removed
- (gfcCommonHistRes6->GetHistDiffFlow())->SetMarkerStyle(markerStyleGFC6);//to be removed
(gfcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC6);
(gfcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC6);
(gfcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC6);
}
if(gfcCommonHistRes8)
{
- (gfcCommonHistRes8->GetHistDiffFlow())->SetMarkerColor(markerColorGFC8);//to be removed
- (gfcCommonHistRes8->GetHistDiffFlow())->SetMarkerStyle(markerStyleGFC8);//to be removed
(gfcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC8);
(gfcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC8);
(gfcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC8);
}
if(gfcCommonHist)
{
- avMultDiffFlowGFC = (gfcCommonHist->GetHistMultDiff())->GetMean();//to be removed
- nEvtsDiffFlowGFC = (gfcCommonHist->GetHistMultDiff())->GetEntries();//to be removed
avMultDiffFlowGFCRP = (gfcCommonHist->GetHistMultInt())->GetMean();
nEvtsDiffFlowGFCRP = (gfcCommonHist->GetHistMultInt())->GetEntries();
avMultDiffFlowGFCPOI = (gfcCommonHist->GetHistMultDiff())->GetMean();
}
//QC = Q-cumulants
- Double_t avMultDiffFlowQC2=0., avMultDiffFlowQC4=0.;//to be removed
- Double_t nEvtsDiffFlowQC2=0., nEvtsDiffFlowQC4=0.;//to be removed
- Double_t avMultDiffFlowQC6=0., avMultDiffFlowQC8=0.;//to be removed
- Double_t nEvtsDiffFlowQC6=0., nEvtsDiffFlowQC8=0.;//to be removed
Double_t avMultDiffFlowQC2RP=0.;
Double_t nEvtsDiffFlowQC2RP=0.;
Double_t avMultDiffFlowQC2POI=0.;
//QC{2}
if(qcCommonHistRes2)
{
- (qcCommonHistRes2->GetHistDiffFlow())->SetMarkerColor(markerColorQC2);//to be removed
- (qcCommonHistRes2->GetHistDiffFlow())->SetMarkerStyle(markerStyleQC2);//to be removed
(qcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC2);
(qcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC2);
(qcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC2);
}
if(qcCommonHist2)
{
- avMultDiffFlowQC2 = (qcCommonHist2->GetHistMultDiff())->GetMean();//to be removed
- nEvtsDiffFlowQC2 = (qcCommonHist2->GetHistMultDiff())->GetEntries();//to be removed
avMultDiffFlowQC2RP = (qcCommonHist2->GetHistMultInt())->GetMean();
nEvtsDiffFlowQC2RP = (qcCommonHist2->GetHistMultInt())->GetEntries();
avMultDiffFlowQC2POI = (qcCommonHist2->GetHistMultDiff())->GetMean();
//QC{4}
if(qcCommonHistRes4)
{
- (qcCommonHistRes4->GetHistDiffFlow())->SetMarkerColor(markerColorQC4);//to be removed
- (qcCommonHistRes4->GetHistDiffFlow())->SetMarkerStyle(markerStyleQC4);//to be removed
(qcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC4);
(qcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC4);
(qcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC4);
}
if(qcCommonHist4)
{
- avMultDiffFlowQC4 = (qcCommonHist4->GetHistMultDiff())->GetMean();//to be removed
- nEvtsDiffFlowQC4 = (qcCommonHist4->GetHistMultDiff())->GetEntries();//to be removed
avMultDiffFlowQC4RP = (qcCommonHist4->GetHistMultInt())->GetMean();
nEvtsDiffFlowQC4RP = (qcCommonHist4->GetHistMultInt())->GetEntries();
avMultDiffFlowQC4POI = (qcCommonHist4->GetHistMultDiff())->GetMean();
//QC{6}
if(qcCommonHistRes6)
{
- (qcCommonHistRes6->GetHistDiffFlow())->SetMarkerColor(markerColorQC6);//to be removed
- (qcCommonHistRes6->GetHistDiffFlow())->SetMarkerStyle(markerStyleQC6);//to be removed
(qcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC6);
(qcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC6);
(qcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC6);
}
if(qcCommonHist6)
{
- avMultDiffFlowQC6 = (qcCommonHist6->GetHistMultDiff())->GetMean();//to be removed
- nEvtsDiffFlowQC6 = (qcCommonHist6->GetHistMultDiff())->GetEntries();//to be removed
avMultDiffFlowQC6RP = (qcCommonHist6->GetHistMultInt())->GetMean();
nEvtsDiffFlowQC6RP = (qcCommonHist6->GetHistMultInt())->GetEntries();
avMultDiffFlowQC6POI = (qcCommonHist6->GetHistMultDiff())->GetMean();
//QC{8}
if(qcCommonHistRes8)
{
- (qcCommonHistRes8->GetHistDiffFlow())->SetMarkerColor(markerColorQC8);//to be removed
- (qcCommonHistRes8->GetHistDiffFlow())->SetMarkerStyle(markerStyleQC8);//to be removed
(qcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC8);
(qcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC8);
(qcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC8);
}
if(qcCommonHist8)
{
- avMultDiffFlowQC8 = (qcCommonHist8->GetHistMultDiff())->GetMean();//to be removed
- nEvtsDiffFlowQC8 = (qcCommonHist8->GetHistMultDiff())->GetEntries();//to be removed
avMultDiffFlowQC8RP = (qcCommonHist8->GetHistMultInt())->GetMean();
nEvtsDiffFlowQC8RP = (qcCommonHist8->GetHistMultInt())->GetEntries();
avMultDiffFlowQC8POI = (qcCommonHist8->GetHistMultDiff())->GetMean();
}
//LYZ2 = Lee-Yang Zeros (2nd run)
- Double_t avMultDiffFlowLYZ2=0.;//to be removed
- Double_t nEvtsDiffFlowLYZ2=0;//to be removed
Double_t avMultDiffFlowLYZ2RP=0.;
Double_t nEvtsDiffFlowLYZ2RP=0;
Double_t avMultDiffFlowLYZ2POI=0.;
{
if(lyz2CommonHistRes)
{
- (lyz2CommonHistRes->GetHistDiffFlow())->SetMarkerColor(markerColorLYZ2);//to be removed
- (lyz2CommonHistRes->GetHistDiffFlow())->SetMarkerStyle(markerStyleLYZ2);//to be removed
(lyz2CommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorLYZ2);
(lyz2CommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleLYZ2);
(lyz2CommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorLYZ2);
}
if(lyz2CommonHist)
{
- avMultDiffFlowLYZ2 = (lyz2CommonHist->GetHistMultDiff())->GetMean();//to be removed
- nEvtsDiffFlowLYZ2 = (lyz2CommonHist->GetHistMultDiff())->GetEntries();//to be removed
avMultDiffFlowLYZ2RP = (lyz2CommonHist->GetHistMultInt())->GetMean();
nEvtsDiffFlowLYZ2RP = (lyz2CommonHist->GetHistMultInt())->GetEntries();
avMultDiffFlowLYZ2POI = (lyz2CommonHist->GetHistMultDiff())->GetMean();
}
//LYZEP = Lee-Yang Zeros Event Plane
- Double_t avMultDiffFlowLYZEP=0.;//to be removed
- Double_t nEvtsDiffFlowLYZEP=0;//to be removed
Double_t avMultDiffFlowLYZEPRP=0.;
Double_t nEvtsDiffFlowLYZEPRP=0;
Double_t avMultDiffFlowLYZEPPOI=0.;
{
if(lyzepCommonHistRes)
{
- (lyzepCommonHistRes->GetHistDiffFlow())->Scale(0.01);//to be improved
- (lyzepCommonHistRes->GetHistDiffFlow())->SetMarkerColor(markerColorLYZEP);//to be removed
- (lyzepCommonHistRes->GetHistDiffFlow())->SetMarkerStyle(markerStyleLYZEP);//to be removed
(lyzepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorLYZEP);
(lyzepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleLYZEP);
(lyzepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorLYZEP);
}
if(lyzepCommonHist)
{
- avMultDiffFlowLYZEP = (lyzepCommonHist->GetHistMultDiff())->GetMean();//to be removed
- nEvtsDiffFlowLYZEP = (lyzepCommonHist->GetHistMultDiff())->GetEntries();//to be removed
avMultDiffFlowLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetMean();
nEvtsDiffFlowLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetEntries();
avMultDiffFlowLYZEPPOI = (lyzepCommonHist->GetHistMultDiff())->GetMean();
}
}
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
- //----------------------------------------------------------------------------------
- //final drawing for differential flow:
- TCanvas* diffFlowAllCanvas = new TCanvas("Differential Flow","Differential Flow",1000,600);
-
- diffFlowAllCanvas->Divide(2,1);
-
- //1st pad is for plot:
- (diffFlowAllCanvas->cd(1))->SetPad(0.0,0.0,0.75,1.0);
-
- if(styleHistPt)
- {
- styleHistPt->Draw();
- }
-
- if(pMeshDiffFlow)
- {
- pMeshDiffFlow->Draw("LFSAME");
- }
- //MC
- if(mcepCommonHistRes)
- {
- (mcepCommonHistRes->GetHistDiffFlow())->Draw("E1PSAME");
- }
- //GFC
- if(gfcCommonHistRes2)
- {
- (gfcCommonHistRes2->GetHistDiffFlow())->Draw("E1PSAME");
- }
- if(gfcCommonHistRes4)
- {
- (gfcCommonHistRes4->GetHistDiffFlow())->Draw("E1PSAME");
- }
- if(gfcCommonHistRes6)
- {
- (gfcCommonHistRes6->GetHistDiffFlow())->Draw("E1PSAME");
- }
- if(gfcCommonHistRes8)
- {
- (gfcCommonHistRes8->GetHistDiffFlow())->Draw("E1PSAME");
- }
- //QC
- if(qcCommonHistRes2)
- {
- (qcCommonHistRes2->GetHistDiffFlow())->Draw("E1PSAME");
- }
- if(qcCommonHistRes4)
- {
- (qcCommonHistRes4->GetHistDiffFlow())->Draw("E1PSAME");
- }
- //LYZ2
- if(lyz2CommonHistRes)
- {
- (lyz2CommonHistRes->GetHistDiffFlow())->Draw("E1PSAME");
- }
- //LYZEP
- if(lyzepCommonHistRes)
- {
- (lyzepCommonHistRes->GetHistDiffFlow())->Draw("E1PSAME");
- }
-
- //2nd pad is for legend:
- (diffFlowAllCanvas->cd(2))->SetPad(0.75,0.0,1.0,1.0);
-
- TLegend* legendDiffFlow = new TLegend(0.02,0.25,0.97,0.75);
- legendDiffFlow->SetTextFont(72);
- legendDiffFlow->SetTextSize(0.06);
-
- //legend's entries:
- TString *entryDiffMC = new TString("MC ....... ");
- TString *entryDiffGFC2 = new TString("GFC{2} ... ");
- TString *entryDiffGFC4 = new TString("GFC{4} ... ");
- TString *entryDiffGFC6 = new TString("GFC{6} ... ");
- TString *entryDiffGFC8 = new TString("GFC{8} ... ");
- TString *entryDiffQC2 = new TString("QC{2} .... ");
- TString *entryDiffQC4 = new TString("QC{4} .... ");
- TString *entryDiffLYZ2 = new TString("LYZ ...... ");
- TString *entryDiffLYZEP = new TString("LYZEP ... ");
-
- //MC
- if(mcepCommonHistRes)
- {
- (mcepCommonHistRes->GetHistDiffFlow())->SetFillStyle(meshStyle);
- (mcepCommonHistRes->GetHistDiffFlow())->SetFillColor(meshColor);
- entryDiffMC->Append("M = ");
- (*entryDiffMC)+=(Long_t)avMultDiffFlowMC;
- entryDiffMC->Append(", N = ");
- (*entryDiffMC)+=(Long_t)nEvtsDiffFlowMC;
- legendDiffFlow->AddEntry(mcepCommonHistRes->GetHistDiffFlow(),entryDiffMC->Data(),"f");
- }
-
- //GFC
- if(gfcCommonHistRes2)
- {
- entryDiffGFC2->Append("M = ");
- (*entryDiffGFC2)+=(Long_t)avMultDiffFlowGFC;
- entryDiffGFC2->Append(", N = ");
- (*entryDiffGFC2)+=(Long_t)nEvtsDiffFlowGFC;
- legendDiffFlow->AddEntry(gfcCommonHistRes2->GetHistDiffFlow(),entryDiffGFC2->Data(),"p");
- }
- if(gfcCommonHistRes4)
- {
- entryDiffGFC4->Append("M = ");
- (*entryDiffGFC4)+=(Long_t)avMultDiffFlowGFC;
- entryDiffGFC4->Append(", N = ");
- (*entryDiffGFC4)+=(Long_t)nEvtsDiffFlowGFC;
- legendDiffFlow->AddEntry(gfcCommonHistRes4->GetHistDiffFlow(),entryDiffGFC4->Data(),"p");
- }
- if(gfcCommonHistRes6)
- {
- entryDiffGFC6->Append("M = ");
- (*entryDiffGFC6)+=(Long_t)avMultDiffFlowGFC;
- entryDiffGFC6->Append(", N = ");
- (*entryDiffGFC6)+=(Long_t)nEvtsDiffFlowGFC;
- legendDiffFlow->AddEntry(gfcCommonHistRes6->GetHistDiffFlow(),entryDiffGFC6->Data(),"p");
- }
- if(gfcCommonHistRes8)
- {
- entryDiffGFC8->Append("M = ");
- (*entryDiffGFC8)+=(Long_t)avMultDiffFlowGFC;
- entryDiffGFC8->Append(", N = ");
- (*entryDiffGFC8)+=(Long_t)nEvtsDiffFlowGFC;
- legendDiffFlow->AddEntry(gfcCommonHistRes8->GetHistDiffFlow(),entryDiffGFC8->Data(),"p");
- }
- //QC
- if(qcCommonHistRes2)
- {
- entryDiffQC2->Append("M = ");
- (*entryDiffQC2)+=(Long_t)avMultDiffFlowQC2;
- entryDiffQC2->Append(", N = ");
- (*entryDiffQC2)+=(Long_t)nEvtsDiffFlowQC2;
- legendDiffFlow->AddEntry(qcCommonHistRes2->GetHistDiffFlow(),entryDiffQC2->Data(),"p");
- }
- if(qcCommonHistRes4)
- {
- entryDiffQC4->Append("M = ");
- (*entryDiffQC4)+=(Long_t)avMultDiffFlowQC4;
- entryDiffQC4->Append(", N = ");
- (*entryDiffQC4)+=(Long_t)nEvtsDiffFlowQC4;
- legendDiffFlow->AddEntry(qcCommonHistRes4->GetHistDiffFlow(),entryDiffQC4->Data(),"p");
- }
-
- //LYZ
- if(lyz2CommonHistRes)
- {
- entryDiffLYZ2->Append("M = ");
- (*entryDiffLYZ2)+=(Long_t)avMultDiffFlowLYZ2;
- entryDiffLYZ2->Append(", N = ");
- (*entryDiffLYZ2)+=(Long_t)nEvtsDiffFlowLYZ2;
- legendDiffFlow->AddEntry(lyz2CommonHistRes->GetHistDiffFlow(),entryDiffLYZ2->Data(),"p");
- }
-
- //LYZEP
- if(lyzepCommonHistRes)
- {
- entryDiffLYZEP->Append("M = ");
- (*entryDiffLYZEP)+=(Long_t)avMultDiffFlowLYZEP;
- entryDiffLYZEP->Append(", N = ");
- (*entryDiffLYZEP)+=(Long_t)nEvtsDiffFlowLYZEP;
- legendDiffFlow->AddEntry(lyzepCommonHistRes->GetHistDiffFlow(),entryDiffLYZEP->Data(),"p");
- }
-
-
- //drawing finally the legend in the 2nd pad:
- if(legendDiffFlow)
- {
- legendDiffFlow->SetMargin(0.15);
- legendDiffFlow->Draw();
- }
- //----------------------------------------------------------------------------------
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-
//----------------------------------------------------------------------------------
//final drawing for differential flow (Pt, RP):
//set here the results of which methods will be plotted by default:
//MC
if(plotMCPtPOI && mcepCommonHistRes)
- {
+ { kTRUE;
(mcepCommonHistRes->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
}
//GFC
(qcCommonHistRes4->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
}
if(plotQC6PtPOI && qcCommonHistRes6)
- {
+ { kTRUE;
//(qcCommonHistRes6->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
}
if(plotQC8PtPOI && qcCommonHistRes8)