kNegative,
kNCharges
};
-
+Bool_t gIsTrkletsTrklets=kFALSE;
+Bool_t gIsMuTrklets=kFALSE;
+
void SetRanges(TAxis* axis)
{
if (strcmp(axis->GetTitle(), "leading p_{T} (GeV/c)") == 0)
*hist = h->GetUEHist(2)->GetSumOfRatios2(hMixed->GetUEHist(2), step, 0, ptBegin, ptEnd, centralityBeginBin, centralityEndBin, normalizePerTrigger);
TString str;
- str.Form("%.1f < p_{T,trig} < %.1f", ptBegin - 0.01, ptEnd + 0.01);
+ str.Form("%.1f < %s < %.1f", ptBegin - 0.001,(gIsTrkletsTrklets)?"Dphi_{trig}":"p_{T,trig}", ptEnd + 0.01);
TString str2;
- str2.Form("%.2f < p_{T,assoc} < %.2f", gpTMin - 0.01, gpTMax + 0.01);
+ str2.Form("%.1f < %s < %.1f", gpTMin - 0.001, (gIsTrkletsTrklets || gIsMuTrklets)?"Dphi_{ass}":"p_{T,ass}",gpTMax + 0.01);
TString newTitle;
newTitle.Form("%s - %s - %d-%d", str.Data(), str2.Data(), centralityBegin, centralityEnd);
newTitle += "%";
if ((*hist))
(*hist)->SetTitle(newTitle);
+ //Printf("\n\n\n\n %s\n\n\n",newTitle.Data());
}
void PlotDeltaPhiDistributions(const char* fileName1, const char* fileName2, Float_t yMax = 0.1, Int_t twoD = 0, Int_t centrBegin = 1, Int_t centrEnd = 1)
Bool_t symmetrizePt = kFALSE;
Int_t maxLeadingPt = 4;
Int_t maxAssocPt = 6;
- if (1)
+ if (0)
{
//PbPb, NS peak shapes
Float_t leadingPtArr[] = { 1.0, 2.0, 3.0, 4.0, 8.0, 15.0, 20.0 };
Float_t leadingPtArr[] = { 3.0, 5.0, 8.0, 16.0 };
Float_t assocPtArr[] = {0.15, 0.3, 50.0, 0.5, 50.0, 1.0, 50.0 };
}
+ else if (gIsMuTrklets)
+ {
+ //mu-trklets
+ maxLeadingPt = 8;
+ maxAssocPt = 6;
+ Float_t leadingPtArr[] = {., 0.5, 1.0, 1.5, 2.0, 3.0, 4.0, 6.0, 10.0};
+ Float_t assocPtArr[] = {0.0,0.0,0.0, 1.0, 2.0,5.0,10. };
+ gIsMuTrklets=kTRUE;
+ }
+ else if (gIsTrkletsTrklets)
+ {
+ //trklets-trklets
+ maxLeadingPt = 4;
+ maxAssocPt = 6;
+ Float_t leadingPtArr[] = {0.0, 1.0,2.0,5.0,10.0};
+ Float_t assocPtArr[] = {0.0,0.0,0.0, 1.0, 2.0,5.0,10.0};
+ gIsTrkletsTrklets=kTRUE;
+ }
else if (0)
{
//pA, trigger from all pT
Float_t leadingPtArr[] = { 0.5, 1.0, 1.5, 2.0, 2.5, 4.0, 8.0, 15.0, 20.0 };
Float_t assocPtArr[] = { 0.15, 0.5, 1.0, 1.5, 2.0, 2.5, 4.0, 8.0, 10.0, 12.0 };
}
+ else if (1)
+ {
+ //pp, for Ilya
+ maxLeadingPt = 5;
+ maxAssocPt = 6;
+ Float_t leadingPtArr[] = { 0.5, 1.0, 2.0, 3.0, 4.0, 8.0, 20. };
+ Float_t assocPtArr[] = { 0.15, 0.5, 1.0, 2.0, 3.0, 4.0, 8.0, 20. };
+ }
+ else if (1)
+ {
+ //pA, for Tim
+ maxLeadingPt = 3;
+ maxAssocPt = 4;
+ Float_t leadingPtArr[] = { 0.5, 1.0, 2.0, 4.0 };
+ Float_t assocPtArr[] = { 0.15, 0.5, 1.0, 2.0, 4.0 };
+ }
else if (0)
{
//pPb and PbPb with PID, with low pt points TPC only
Float_t leadingPtArr[] = { 0.2, 0.3, 0.5 , 0.75, 1.0, 1.25, 1.5, 2.0, 2.5, 3., 4.0};
Float_t assocPtArr[] = { 0.15, 0.2, 0.3, 0.5 , 0.75, 1.0, 1.25, 1.5, 2.0, 2.5, 3., 4.0 };
}
- else if (1)
+ else if (0)
{
//PbPb for comaprison with You
maxLeadingPt = 8;
h->GetUEHist(2)->SymmetrizepTBins();
hMixed->GetUEHist(2)->SymmetrizepTBins();
}
+ if(0){
+ h->GetUEHist(2)->SetSkipScaleMixedEvent(1);
+ hMixed->GetUEHist(2)->SetSkipScaleMixedEvent(1);
+ }
TList* list2 = 0;
AliUEHistograms* h2 = 0;
if (refMultRaw)
{
// new TCanvas; refMultRaw->Draw("COLZ");
- Int_t nCentrBins = 4;
+ //Int_t nCentrBins = 4;
+ Int_t nCentrBins = 5;
+ Double_t centrBins[] = { 0., 20., 40., 60., 80., 100. };
Double_t centrBins[] = { 0., 20., 40., 60., 100. };
-// Double_t centrBins[] = { 0., 3., 10., 50., 100. };
+
+ // Double_t centrBins[] = { 0., 3., 10., 50., 100. };
TH1* refMult = new TH1F("refMult", ";centrality;<Nch>", nCentrBins, centrBins);
for (Int_t i=0; i<nCentrBins; i++)
{
for (Int_t i=0; i<maxLeadingPt; i++)
{
- for (Int_t j=2; j<maxAssocPt; j++)
+ for (Int_t j=1; j<maxAssocPt; j++)
{
if(0){
if(j!=(i+1))continue;
gpTMin = assocPtArr[j] + 0.01;
gpTMax = assocPtArr[j+1] - 0.01;
+
+ if(gIsTrkletsTrklets || gIsMuTrklets)=gpTMin=0.01;
if(gpTMin >= gpTMax)continue;
SetupRanges(hMixed3);
// SetupRanges(hMixed3);
- if(1) if (assocPtArr[j] >= leadingPtArr[i+leadingPtOffset])
+ if(0) if (assocPtArr[j] >= leadingPtArr[i+leadingPtOffset])
continue;
+ TH1* hist0 = 0;
TH1* hist1 = 0;
TH1* hist2 = 0;
TH1* hist3 = 0;
TH1* hist6 = 0;
TH1* hist7 = 0;
TH1* hist8 = 0;
+ TH1* hist9 = 0;
+ TH1* hist10 = 0;
Bool_t equivMixedBin = 1; //kFALSE; // TODO ?
Bool_t scaleToPairs = kTRUE;
Int_t histType = 1;
- if (1)
+ if (0)
{
// PbPb
Int_t step = 8;
if (h2)
GetSumOfRatios(h2, hMixed2, &hist3, step, 0, -1, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
}
- else if (0)
+ else if (1)
+ {
+ // pp vs Mult for Ilya
+ Int_t step = 8;
+// Int_t step = 0;
+
+ GetSumOfRatios(h, hMixed, &hist0, step, 0, 1, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist1, step, 0, 10, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist2, step, 10, 20, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist3, step, 20, 30, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist4, step, 30, 40, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist5, step, 40, 50, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist6, step, 50, 60, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist7, step, 60, 70, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist8, step, 70, 80, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist9, step, 80, 90, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist10,step, 90,100, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ }
+ else if (1)
{
// pA, fine binning
Int_t step = 8;
if (h3)
GetSumOfRatios(h3, hMixed3, &hist6, step, 0, -1, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
}
+ else if (1)
+ {
+ // trkletstrklets
+ Int_t step = 6;
+ //Int_t step = 0;
+ GetSumOfRatios(h, hMixed, &hist1, step, 0, 20, (gIsTrkletsTrklets)?0.01:leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist2, step, 20, 40, (gIsTrkletsTrklets)?0.01:leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist4, step, 40, 60, (gIsTrkletsTrklets)?0.01:leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ GetSumOfRatios(h, hMixed, &hist5, step, 60, 100, (gIsTrkletsTrklets)?0.01:leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ // GetSumOfRatios(h, hMixed, &hist2, step, 20, 40, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ // GetSumOfRatios(h, hMixed, &hist4, step, 40, 60, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ // GetSumOfRatios(h, hMixed, &hist5, step, 60, 100, leadingPtArr[i] + 0.01, leadingPtArr[i+leadingPtOffset] - 0.01, kTRUE);
+ }
else if (0)
{
// pA, V0 as trigger particle detector
file = TFile::Open(outputFile, "UPDATE");
+ if (hist0)
+ {
+ hist0->SetName(Form("dphi_%d_%d_%d", i, j, 0));
+ hist0->Write();
+ }
if (hist1)
- {
- hist1->SetName(Form("dphi_%d_%d_%d", i, j, 0));
- hist1->Write();
- }
-
+ {
+ hist1->SetName(Form("dphi_%d_%d_%d", i, j, 1));
+ hist1->Write();
+ }
if (hist2)
- {
- hist2->SetName(Form("dphi_%d_%d_%d", i, j, 1));
- hist2->Write();
- }
-
+ {
+ hist2->SetName(Form("dphi_%d_%d_%d", i, j, 2));
+ hist2->Write();
+ }
+ if (hist3)
+ {
+ hist3->SetName(Form("dphi_%d_%d_%d", i, j, 3));
+ hist3->Write();
+ }
if (hist4)
- {
- hist4->SetName(Form("dphi_%d_%d_%d", i, j, 3));
- hist4->Write();
- }
-
+ {
+ hist4->SetName(Form("dphi_%d_%d_%d", i, j, 4));
+ hist4->Write();
+ }
if (hist5)
- {
- hist5->SetName(Form("dphi_%d_%d_%d", i, j, 4));
- hist5->Write();
- }
-
+ {
+ hist5->SetName(Form("dphi_%d_%d_%d", i, j, 5));
+ hist5->Write();
+ }
+ if (hist6)
+ {
+ hist6->SetName(Form("dphi_%d_%d_%d", i, j, 6));
+ hist6->Write();
+ }
if (hist7)
- {
- hist7->SetName(Form("dphi_%d_%d_%d", i, j, 6));
- hist7->Write();
- }
-
+ {
+ hist7->SetName(Form("dphi_%d_%d_%d", i, j, 7));
+ hist7->Write();
+ }
if (hist8)
- {
- hist8->SetName(Form("dphi_%d_%d_%d", i, j, 7));
- hist8->Write();
- }
-
- if (hist6)
- {
- hist6->SetName(Form("dphi_%d_%d_%d", i, j, 5));
- TString title(hist6->GetTitle());
- title.ReplaceAll("0--1%", "pp");
- hist6->SetTitle(title);
- hist6->Write();
- }
+ {
+ hist8->SetName(Form("dphi_%d_%d_%d", i, j, 8));
+ hist8->Write();
+ }
+ if (hist9)
+ {
+ hist9->SetName(Form("dphi_%d_%d_%d", i, j, 9));
+ hist9->Write();
+ }
+ if (hist10)
+ {
+ hist10->SetName(Form("dphi_%d_%d_%d", i, j, 10));
+ hist10->Write();
+ }
- if (hist3)
- {
- hist3->SetName(Form("dphi_%d_%d_%d", i, j, 2));
- TString title(hist3->GetTitle());
- title.ReplaceAll("0--1%", "pp");
- hist3->SetTitle(title);
- hist3->Write();
- }
+ // if (hist1)
+ // {
+ // hist1->SetName(Form("dphi_%d_%d_%d", i, j, 0));
+ // hist1->Write();
+ // }
+
+ // if (hist2)
+ // {
+ // hist2->SetName(Form("dphi_%d_%d_%d", i, j, 1));
+ // hist2->Write();
+ // }
+
+ // if (hist4)
+ // {
+ // hist4->SetName(Form("dphi_%d_%d_%d", i, j, 3));
+ // hist4->Write();
+ // }
+
+ // if (hist5)
+ // {
+ // hist5->SetName(Form("dphi_%d_%d_%d", i, j, 4));
+ // hist5->Write();
+ // }
+
+ // if (hist7)
+ // {
+ // hist7->SetName(Form("dphi_%d_%d_%d", i, j, 6));
+ // hist7->Write();
+ // }
+
+ // if (hist8)
+ // {
+ // hist8->SetName(Form("dphi_%d_%d_%d", i, j, 7));
+ // hist8->Write();
+ // }
+
+ // if (hist6)
+ // {
+ // hist6->SetName(Form("dphi_%d_%d_%d", i, j, 5));
+ // TString title(hist6->GetTitle());
+ // title.ReplaceAll("0--1%", "pp");
+ // hist6->SetTitle(title);
+ // hist6->Write();
+ // }
+
+ // if (hist3)
+ // {
+ // hist3->SetName(Form("dphi_%d_%d_%d", i, j, 2));
+ // TString title(hist3->GetTitle());
+ // title.ReplaceAll("0--1%", "pp");
+ // hist3->SetTitle(title);
+ // hist3->Write();
+ // }
file->Close();
+ delete hist0;
delete hist1;
delete hist2;
delete hist3;
delete hist6;
delete hist7;
delete hist8;
+ delete hist9;
+ delete hist10;
// return;
}
else
newFileName += "_zvtx.root";
- file3 = TFile::Open(newFileName, "RECREATE");
+ TFile *file3 = TFile::Open(newFileName, "RECREATE");
file3->mkdir("PWG4_PhiCorrelations");
file3->cd("PWG4_PhiCorrelations");
list->Write("histosPhiCorrelations", TObject::kSingleKey);
file3->Close();
+ delete file3;
}
void CompareZVertex(const char* fileName)
fMinTPCNcls(70),
fFillTHn(kTRUE),
fResSP(0),
+ f2dResSP(0),
+ fResSP_vs_Cent(0),
+ f2partCumQA_vs_Cent(0),
+ f2partCumQB_vs_Cent(0),
fEta_vs_Phi_bef(0),
fEta_vs_PhiA(0),
fEta_vs_PhiB(0),
fResSP_lq(0),
- fResSP_sq(0)
+ f2dResSP_lq(0),
+ fResSP_vs_Cent_lq(0),
+ f2partCumQA_vs_Cent_lq(0),
+ f2partCumQB_vs_Cent_lq(0),
+ fResSP_sq(0),
+ f2dResSP_sq(0),
+ fResSP_vs_Cent_sq(0),
+ f2partCumQA_vs_Cent_sq(0),
+ f2partCumQB_vs_Cent_sq(0)
{
for (Int_t i = 0; i< 9; i++){
if( fFillTHn ){
//dimensions of THnSparse for Q vector checks
- const Int_t nvarev=4;
- // cent Q vec Qa Qb
- Int_t binsHistRealEv[nvarev] = { 10, fnQvecBins, 100, 100};
- Double_t xminHistRealEv[nvarev] = { 0., 0., 0., 0.};
- Double_t xmaxHistRealEv[nvarev] = { 10., fQvecUpperLim, 10., 10.};
+ const Int_t nvarev=5;
+ // cent Q vec Qa Qb res
+ Int_t binsHistRealEv[nvarev] = { fnCentBins, fnQvecBins, 100, 100, 100};
+ Double_t xminHistRealEv[nvarev] = { 0., 0., 0., 0., 0.};
+ Double_t xmaxHistRealEv[nvarev] = { 100., fQvecUpperLim, 10., 10., 1.};
THnSparseF* NSparseHistEv = new THnSparseF("NSparseHistEv","NSparseHistEv",nvarev,binsHistRealEv,xminHistRealEv,xmaxHistRealEv);
NSparseHistEv->GetAxis(0)->SetTitle(Form("%s cent",fEventCuts->GetCentralityMethod().Data()));
NSparseHistEv->GetAxis(0)->SetName(Form("%s_cent",fEventCuts->GetCentralityMethod().Data()));
NSparseHistEv->GetAxis(2)->SetName("Qvec_A");
NSparseHistEv->GetAxis(3)->SetTitle("Q_vec (B)");
NSparseHistEv->GetAxis(3)->SetName("Qvec_B");
+ NSparseHistEv->GetAxis(4)->SetTitle("resolution");
+ NSparseHistEv->GetAxis(4)->SetName("res");
fOutput->Add(NSparseHistEv);
}
fResSP = new TProfile("fResSP", "Resolution; centrality; Resolution", 9, -0.5, 8.5);
fOutput->Add(fResSP);
+ f2dResSP = new TH2D("f2dResSP", "Resolution; centrality; Resolution", 9, -0.5, 8.5,100.,0.,1.);
+ fOutput->Add(f2dResSP);
+
+ fResSP_vs_Cent = new TH2D("fResSP_vs_Cent", "Resolution; centrality; Resolution", 100., 0., 100.,100.,0.,1.);
+ fOutput->Add(fResSP_vs_Cent);
+
+ f2partCumQA_vs_Cent = new TH2D("f2partCumQA_vs_Cent", "Resolution; centrality; Resolution", 100., 0., 100.,200.,-1.,1.);
+ fOutput->Add(f2partCumQA_vs_Cent);
+
+ f2partCumQB_vs_Cent = new TH2D("f2partCumQB_vs_Cent", "Resolution; centrality; Resolution", 100., 0., 100.,200.,-1.,1.);
+ fOutput->Add(f2partCumQB_vs_Cent);
+
fEta_vs_Phi_bef = new TH2D("fEta_vs_Phi_bef","eta vs phi distribution before eta gap;#eta;#phi",200.,-1.,1.,350.,0.,7.);
fOutput->Add(fEta_vs_Phi_bef);
fResSP_lq = new TProfile("fResSP_lq", "Resolution; centrality; Resolution", 9, -0.5, 8.5);
fOutput_lq->Add(fResSP_lq);
+ f2dResSP_lq = new TH2D("f2dResSP_lq", "Resolution; centrality; Resolution", 9, -0.5, 8.5,100.,0.,1.);
+ fOutput_lq->Add(f2dResSP_lq);
+
+ fResSP_vs_Cent_lq = new TH2D("fResSP_vs_Cent_lq", "Resolution; centrality; Resolution", 100., 0., 100.,100.,0.,1.);
+ fOutput_lq->Add(fResSP_vs_Cent_lq);
+
+ f2partCumQA_vs_Cent_lq = new TH2D("f2partCumQA_vs_Cent_lq", "Resolution; centrality; Resolution", 100., 0., 100.,200.,-1.,1.);
+ fOutput_lq->Add(f2partCumQA_vs_Cent_lq);
+
+ f2partCumQB_vs_Cent_lq = new TH2D("f2partCumQB_vs_Cent_lq", "Resolution; centrality; Resolution", 100., 0., 100.,200.,-1.,1.);
+ fOutput_lq->Add(f2partCumQB_vs_Cent_lq);
+
//small q resolution
fResSP_sq = new TProfile("fResSP_sq", "Resolution; centrality; Resolution", 9, -0.5, 8.5);
fOutput_sq->Add(fResSP_sq);
+
+ f2dResSP_sq = new TH2D("f2dResSP_sq", "Resolution; centrality; Resolution", 9, -0.5, 8.5,100.,0.,1.);
+ fOutput_sq->Add(f2dResSP_sq);
+
+ fResSP_vs_Cent_sq = new TH2D("fResSP_vs_Cent_sq", "Resolution; centrality; Resolution", 100., 0., 100.,100.,0.,1.);
+ fOutput_sq->Add(fResSP_vs_Cent_sq);
+
+ f2partCumQA_vs_Cent_sq = new TH2D("f2partCumQA_vs_Cent_sq", "Resolution; centrality; Resolution", 100., 0., 100.,200.,-1.,1.);
+ fOutput_sq->Add(f2partCumQA_vs_Cent_sq);
+
+ f2partCumQB_vs_Cent_sq = new TH2D("f2partCumQB_vs_Cent_sq", "Resolution; centrality; Resolution", 100., 0., 100.,200.,-1.,1.);
+ fOutput_sq->Add(f2partCumQB_vs_Cent_sq);
for (Int_t iC = 0; iC < 9; iC++){
if (multGap1A > 0 && multGap1B > 0){
Double_t res = (QxGap1A*QxGap1B + QyGap1A*QyGap1B)/(Double_t)multGap1A/(Double_t)multGap1B;
fResSP->Fill((Double_t)centV0, res);
+ f2dResSP->Fill((Double_t)centV0, res);
+ fResSP_vs_Cent->Fill((Double_t)Cent, res);
+
+ Double_t f2partCumQA = -999.;
+ if(multGap1A>1)
+ f2partCumQA = ( ( (QxGap1A*QxGap1A + QyGap1A*QyGap1A) - (Double_t)multGap1A ) / ((Double_t)multGap1A*((Double_t)multGap1A-1)) );
+
+ f2partCumQA_vs_Cent->Fill((Double_t)Cent,f2partCumQA);
+
+ Double_t f2partCumQB = -999.;
+ if(multGap1B>1)
+ f2partCumQB = ( ( (QxGap1B*QxGap1B + QyGap1B*QyGap1B) - (Double_t)multGap1B ) / ((Double_t)multGap1B*((Double_t)multGap1B-1)) );
+
+ f2partCumQB_vs_Cent->Fill((Double_t)Cent,f2partCumQB);
if (Qvec > fCutLargeQperc && Qvec < 100.){
fResSP_lq->Fill((Double_t)centV0, res);
+ f2dResSP_lq->Fill((Double_t)centV0, res);
+ fResSP_vs_Cent_lq->Fill((Double_t)Cent, res);
+ f2partCumQA_vs_Cent_lq->Fill((Double_t)Cent,f2partCumQA);
+ f2partCumQB_vs_Cent_lq->Fill((Double_t)Cent,f2partCumQB);
}
if (Qvec > 0. && Qvec < fCutSmallQperc){
fResSP_sq->Fill((Double_t)centV0, res);
+ f2dResSP_sq->Fill((Double_t)centV0, res);
+ fResSP_vs_Cent_sq->Fill((Double_t)Cent, res);
+ f2partCumQA_vs_Cent_sq->Fill((Double_t)Cent,f2partCumQA);
+ f2partCumQB_vs_Cent_sq->Fill((Double_t)Cent,f2partCumQB);
}
if( fFillTHn ){
Double_t QA = TMath::Sqrt( (QxGap1A*QxGap1A + QyGap1A*QyGap1A)/multGap1A );
Double_t QB = TMath::Sqrt( (QxGap1B*QxGap1B + QyGap1B*QyGap1B)/multGap1B );
- Double_t varEv[4];
- varEv[0]=centV0;
+ Double_t varEv[5];
+ varEv[0]=Cent;
varEv[1]=Qvec;
varEv[2]=(Double_t)QA;
- varEv[3]=(Double_t)QB;
+ varEv[2]=(Double_t)QB;
+ varEv[4]=(Double_t)res;
((THnSparseF*)fOutput->FindObject("NSparseHistEv"))->Fill(varEv);//event loop
}
}