fTracks(0), fV1Clus(0), fV2Clus(0),
fRunNumber(-999),fInterRunNumber(-999),
fVtxCut(15.),
- fNcellCut(2.), fECut(1.), fEtaCut(0.65), fM02Cut(0.5),fDrCut(0.025), fPi0AsyCut(0), isV1Clus(1), isPhosCali(0),
+ fNcellCut(2.), fECut(1.), fEtaCut(0.65), fM02Cut(0.5),fDrCut(0.025), fPi0AsyCut(0), isV1Clus(1), isPhosCali(0), isCentFlat(0), isFullHist(0),
fCentrality(99.),
fEPTPC(-999.),
fEPTPCreso(0.),
fEPV0(-999.), fEPV0A(-999.), fEPV0C(-999.), fEPV0Ar(-999.), fEPV0Cr(-999.), fEPV0r(-999.),
fEPV0AR4(-999.), fEPV0AR5(-999.), fEPV0AR6(-999.), fEPV0AR7(-999.), fEPV0CR0(-999.), fEPV0CR1(-999.), fEPV0CR2(-999.), fEPV0CR3(-999.),
- hEvtCount(0),
+ hEvtCount(0), hCent(0),
h2DcosV0A(0), h2DsinV0A(0), h2DcosV0C(0), h2DsinV0C(0), h2DcosTPC(0), h2DsinTPC(0),
hEPTPC(0), hresoTPC(0),
hEPV0(0), hEPV0A(0), hEPV0C(0), hEPV0Ar(0), hEPV0Cr(0), hEPV0r(0), hEPV0AR4(0), hEPV0AR7(0), hEPV0CR0(0), hEPV0CR3(0),
hdifV0Ar_V0Cr(0), hdifV0A_V0CR0(0), hdifV0A_V0CR3(0), hdifV0ACR0_V0CR3(0), hdifV0C_V0AR4(0), hdifV0C_V0AR7(0), hdifV0AR4_V0AR7(0),
hdifV0A_V0C(0), hdifV0A_TPC(0), hdifV0C_TPC(0), hdifV0C_V0A(0),
hM02vsPtA(0), hM02vsPtB(0), hClusDxDZA(0), hClusDxDZB(0),
- hdifEMC_EPV0(0), hdifEMC_EPV0A(0), hdifEMC_EPV0C(0), hdifful_EPV0(0), hdifful_EPV0A(0), hdifful_EPV0C(0),
- hdifout_EPV0(0), hdifout_EPV0A(0), hdifout_EPV0C(0),
+ hdifEMC_EPV0A(0), hdifEMC_EPV0C(0),
+ hdifful_EPV0A(0), hdifful_EPV0C(0),
+ hdifout_EPV0A(0), hdifout_EPV0C(0),
+ hCv2EMC_EPV0A(0), hCv2EMC_EPV0C(0), hCv2ful_EPV0A(0), hCv2ful_EPV0C(0), hCv2out_EPV0A(0), hCv2out_EPV0C(0),
+ hclusDif_EPV0A(0), hclusDif_EPV0C(0), hclusv2_EPV0A(0), hclusv2_EPV0C(0),
fEPcalibFileName("$ALICE_ROOT/OADB/PHOS/PHOSflat.root"), fTPCFlat(0x0), fV0AFlat(0x0), fV0CFlat(0x0),
fClusterPbV0(0), fClusterPbV0A(0), fClusterPbV0C(0), fClusterPbTPC(0),
- fHEPV0r(0x0), fHEPV0A(0x0), fHEPV0C(0x0), fHEPTPC(0x0)
+ fHEPV0A(0x0), fHEPV0C(0x0), fHEPTPC(0x0),
+ fHEPV0AM2(0x0), fHEPV0CM2(0x0), fHEPTPCM2(0x0)
{
// Dummy constructor ALWAYS needed for I/O.
DefineInput(0, TChain::Class());
fTracks(0), fV1Clus(0), fV2Clus(0),
fRunNumber(-999),fInterRunNumber(-999),
fVtxCut(15.),
- fNcellCut(2.), fECut(1.), fEtaCut(0.65), fM02Cut(0.5), fDrCut(0.025), fPi0AsyCut(0), isV1Clus(1),isPhosCali(0),
+ fNcellCut(2.), fECut(1.), fEtaCut(0.65), fM02Cut(0.5), fDrCut(0.025), fPi0AsyCut(0), isV1Clus(1),isPhosCali(0),isCentFlat(0), isFullHist(0),
fCentrality(99.),
fEPTPC(-999.),
fEPTPCreso(0.),
fEPV0(-999.), fEPV0A(-999.), fEPV0C(-999.), fEPV0Ar(-999.), fEPV0Cr(-999.), fEPV0r(-999.),
fEPV0AR4(-999.), fEPV0AR5(-999.), fEPV0AR6(-999.), fEPV0AR7(-999.), fEPV0CR0(-999.), fEPV0CR1(-999.), fEPV0CR2(-999.), fEPV0CR3(-999.),
- hEvtCount(0),
+ hEvtCount(0), hCent(0),
h2DcosV0A(0), h2DsinV0A(0), h2DcosV0C(0), h2DsinV0C(0), h2DcosTPC(0), h2DsinTPC(0),
hEPTPC(0), hresoTPC(0),
hEPV0(0), hEPV0A(0), hEPV0C(0), hEPV0Ar(0), hEPV0Cr(0), hEPV0r(0), hEPV0AR4(0), hEPV0AR7(0), hEPV0CR0(0), hEPV0CR3(0),
hdifV0Ar_V0Cr(0), hdifV0A_V0CR0(0), hdifV0A_V0CR3(0), hdifV0ACR0_V0CR3(0), hdifV0C_V0AR4(0), hdifV0C_V0AR7(0), hdifV0AR4_V0AR7(0),
hdifV0A_V0C(0), hdifV0A_TPC(0), hdifV0C_TPC(0), hdifV0C_V0A(0),
hM02vsPtA(0), hM02vsPtB(0), hClusDxDZA(0), hClusDxDZB(0),
- hdifEMC_EPV0(0), hdifEMC_EPV0A(0), hdifEMC_EPV0C(0), hdifful_EPV0(0), hdifful_EPV0A(0), hdifful_EPV0C(0),
- hdifout_EPV0(0), hdifout_EPV0A(0), hdifout_EPV0C(0),
+ hdifEMC_EPV0A(0), hdifEMC_EPV0C(0),
+ hdifful_EPV0A(0), hdifful_EPV0C(0),
+ hdifout_EPV0A(0), hdifout_EPV0C(0),
+ hCv2EMC_EPV0A(0), hCv2EMC_EPV0C(0), hCv2ful_EPV0A(0), hCv2ful_EPV0C(0), hCv2out_EPV0A(0), hCv2out_EPV0C(0),
+ hclusDif_EPV0A(0), hclusDif_EPV0C(0), hclusv2_EPV0A(0), hclusv2_EPV0C(0),
fEPcalibFileName("$ALICE_ROOT/OADB/PHOS/PHOSflat.root"), fTPCFlat(0x0), fV0AFlat(0x0), fV0CFlat(0x0),
fClusterPbV0(0), fClusterPbV0A(0), fClusterPbV0C(0), fClusterPbTPC(0),
- fHEPV0r(0x0), fHEPV0A(0x0), fHEPV0C(0x0), fHEPTPC(0x0)
+ fHEPV0A(0x0), fHEPV0C(0x0), fHEPTPC(0x0),
+ fHEPV0AM2(0x0), fHEPV0CM2(0x0), fHEPTPCM2(0x0)
{
// Constructor
// Define input and output slots here (never in the dummy constructor)
Double_t cos2phiV0C = TMath::Cos(2.*(dphiV0C));
Double_t cos2phiTPC = TMath::Cos(2.*(dphiTPC));
- dphiV0A = TVector2::Phi_0_2pi(dphiV0A); if(dphiV0A >TMath::Pi()) dphiV0A -= TMath::Pi();
- dphiV0C = TVector2::Phi_0_2pi(dphiV0C); if(dphiV0C >TMath::Pi()) dphiV0C -= TMath::Pi();
- dphiTPC = TVector2::Phi_0_2pi(dphiTPC); if(dphiTPC >TMath::Pi()) dphiTPC -= TMath::Pi();
+ while(dphiV0A<0.) dphiV0A+=TMath::Pi(); while(dphiV0A>TMath::Pi()) dphiV0A-=TMath::Pi();
+ while(dphiV0C<0.) dphiV0C+=TMath::Pi(); while(dphiV0C>TMath::Pi()) dphiV0C-=TMath::Pi();
+ while(dphiTPC<0.) dphiTPC+=TMath::Pi(); while(dphiTPC>TMath::Pi()) dphiTPC-=TMath::Pi();
- Double_t xV0A[5]; // Match ndims in fH V0A EP
+
+ Double_t xV0A[4]; // Match ndims in fH V0A EP for method 1
xV0A[0] = mass;
xV0A[1] = pt;
xV0A[2] = fCentrality;
xV0A[3] = dphiV0A;
- xV0A[4] = cos2phiV0A;
fHEPV0A->Fill(xV0A);
- Double_t xV0C[5]; // Match ndims in fH V0C EP
+
+ Double_t xV0AM2[4]; // Match ndims in fH V0A EP for method 2
+ xV0AM2[0] = mass;
+ xV0AM2[1] = pt;
+ xV0AM2[2] = fCentrality;
+ xV0AM2[3] = cos2phiV0A;
+ fHEPV0AM2->Fill(xV0AM2);
+
+
+ Double_t xV0C[4]; // Match ndims in fH V0C EP for method 1
xV0C[0] = mass;
xV0C[1] = pt;
xV0C[2] = fCentrality;
xV0C[3] = dphiV0C;
- xV0C[4] = cos2phiV0C;
fHEPV0C->Fill(xV0C);
+ Double_t xV0CM2[4]; // Match ndims in fH V0C EP for method 2
+ xV0CM2[0] = mass;
+ xV0CM2[1] = pt;
+ xV0CM2[2] = fCentrality;
+ xV0CM2[3] = cos2phiV0C;
+ fHEPV0CM2->Fill(xV0CM2);
+
+
if (fEPTPC!=-999.){
- Double_t xTPC[5]; // Match ndims in fH TPC EP
+ Double_t xTPC[4]; // Match ndims in fH TPC EP for method 1
xTPC[0] = mass;
xTPC[1] = pt;
xTPC[2] = fCentrality;
xTPC[3] = dphiTPC;
- xTPC[4] = cos2phiTPC;
fHEPTPC->Fill(xTPC);
+
+ Double_t xTPCM2[4]; // Match ndims in fH TPC EP
+ xTPCM2[0] = mass;
+ xTPCM2[1] = pt;
+ xTPCM2[2] = fCentrality;
+ xTPCM2[3] = cos2phiTPC;
+ fHEPTPCM2->Fill(xTPCM2);
+
}
}
Double_t Et = p1.Et();
Double_t Phi = p1.Phi();
Double_t M02 = c->GetM02();
- Double_t DxClus = c->GetTrackDx();
- Double_t DzClus = c->GetTrackDz();
- Double_t dr = TMath::Sqrt(DxClus*DxClus + DzClus*DzClus);
- Double_t difClusV0A = TVector2::Phi_0_2pi(Phi-EPV0A); if(difClusV0A >TMath::Pi()) difClusV0A -= TMath::Pi();
- Double_t difClusV0C = TVector2::Phi_0_2pi(Phi-EPV0C); if(difClusV0C >TMath::Pi()) difClusV0C -= TMath::Pi();
- Double_t difClusTPC = TVector2::Phi_0_2pi(Phi-EPTPC); if(difClusTPC >TMath::Pi()) difClusTPC -= TMath::Pi();
+ Double_t difClusV0A = Phi-EPV0A;
+ Double_t difClusV0C = Phi-EPV0C;
+ Double_t difClusTPC = Phi-EPTPC;
+ while(difClusV0A<0.) difClusV0A+=TMath::Pi(); while(difClusV0A>TMath::Pi()) difClusV0A-=TMath::Pi();
+ while(difClusV0C<0.) difClusV0C+=TMath::Pi(); while(difClusV0C>TMath::Pi()) difClusV0C-=TMath::Pi();
+ while(difClusTPC<0.) difClusTPC+=TMath::Pi(); while(difClusTPC>TMath::Pi()) difClusTPC-=TMath::Pi();
- Double_t DataV0A[5];
+ Double_t DataV0A[4];
DataV0A[0] = Et;
DataV0A[1] = M02;
DataV0A[2] = fCentrality;
DataV0A[3] = difClusV0A;
- DataV0A[4] = dr;
fClusterPbV0A->Fill(DataV0A);
- Double_t DataV0C[5];
+ Double_t DataV0C[4];
DataV0C[0] = Et;
DataV0C[1] = M02;
DataV0C[2] = fCentrality;
DataV0C[3] = difClusV0C;
- DataV0C[4] = dr;
fClusterPbV0C->Fill(DataV0C);
- Double_t DataTPC[5];
+ Double_t DataTPC[4];
DataTPC[0] = Et;
DataTPC[1] = M02;
DataTPC[2] = fCentrality;
DataTPC[3] = difClusTPC;
- DataTPC[4] = dr;
fClusterPbTPC->Fill(DataTPC);
}
hEvtCount->GetXaxis()->SetBinLabel(7,"ClusterTask");
hEvtCount->GetXaxis()->SetBinLabel(8,"Pass");
fOutput->Add(hEvtCount);
-
+
+ hCent = new TH1F("hCent", "centrality dist. before App. flat cut", 100, 0., 100.);
+ fOutput->Add(hCent);
+
hEPTPC = new TH2F("hEPTPC", "EPTPC vs cent", 100, 0., 100., 100, 0., TMath::Pi());
hresoTPC = new TH2F("hresoTPC", "TPc reso vs cent", 100, 0., 100., 100, 0., 1.);
- hEPV0 = new TH2F("hEPV0", "EPV0 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
hEPV0A = new TH2F("hEPV0A", "EPV0A vs cent", 100, 0., 100., 100, 0., TMath::Pi());
hEPV0C = new TH2F("hEPV0C", "EPV0C vs cent", 100, 0., 100., 100, 0., TMath::Pi());
- hEPV0Ar = new TH2F("hEPV0Ar", "EPV0Ar vs cent", 100, 0., 100., 100, 0., TMath::Pi());
- hEPV0Cr = new TH2F("hEPV0Cr", "EPV0Cr vs cent", 100, 0., 100., 100, 0., TMath::Pi());
- hEPV0r = new TH2F("hEPV0r", "EPV0r vs cent", 100, 0., 100., 100, 0., TMath::Pi());
- hEPV0AR4 = new TH2F("hEPV0AR4", "EPV0AR4 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
- hEPV0AR7 = new TH2F("hEPV0AR7", "EPV0AR7 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
- hEPV0CR0 = new TH2F("hEPV0CR0", "EPV0CR0 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
- hEPV0CR3 = new TH2F("hEPV0CR3", "EPV0CR3 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
fOutput->Add(hEPTPC);
fOutput->Add(hresoTPC);
- fOutput->Add(hEPV0);
fOutput->Add(hEPV0A);
fOutput->Add(hEPV0C);
- fOutput->Add(hEPV0Ar);
- fOutput->Add(hEPV0Cr);
- fOutput->Add(hEPV0r);
- fOutput->Add(hEPV0AR4);
- fOutput->Add(hEPV0AR7);
- fOutput->Add(hEPV0CR0);
- fOutput->Add(hEPV0CR3);
+
+ if(isFullHist){
+ hEPV0 = new TH2F("hEPV0", "EPV0 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
+ hEPV0Ar = new TH2F("hEPV0Ar", "EPV0Ar vs cent", 100, 0., 100., 100, 0., TMath::Pi());
+ hEPV0Cr = new TH2F("hEPV0Cr", "EPV0Cr vs cent", 100, 0., 100., 100, 0., TMath::Pi());
+ hEPV0r = new TH2F("hEPV0r", "EPV0r vs cent", 100, 0., 100., 100, 0., TMath::Pi());
+ hEPV0AR4 = new TH2F("hEPV0AR4", "EPV0AR4 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
+ hEPV0AR7 = new TH2F("hEPV0AR7", "EPV0AR7 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
+ hEPV0CR0 = new TH2F("hEPV0CR0", "EPV0CR0 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
+ hEPV0CR3 = new TH2F("hEPV0CR3", "EPV0CR3 vs cent", 100, 0., 100., 100, 0., TMath::Pi());
+ fOutput->Add(hEPV0);
+ fOutput->Add(hEPV0Ar);
+ fOutput->Add(hEPV0Cr);
+ fOutput->Add(hEPV0r);
+ fOutput->Add(hEPV0AR4);
+ fOutput->Add(hEPV0AR7);
+ fOutput->Add(hEPV0CR0);
+ fOutput->Add(hEPV0CR3);
+ }
hEPTPCCor = new TH2F("hEPTPCCor", "EPTPC vs cent after PHOS Correct", 100, 0., 100., 100, 0., TMath::Pi());
hEPV0ACor = new TH2F("hEPV0ACor", "EPV0A vs cent after PHOS Correct", 100, 0., 100., 100, 0., TMath::Pi());
fOutput->Add(hEPV0ACor);
fOutput->Add(hEPV0CCor);
- hdifV0Ar_V0Cr = new TH2F("hdifV0Ar_V0Cr", "EP Ar-Cr ", 100, 0., 100., 100, -1., 1.);
hdifV0A_V0CR0 = new TH2F("hdifV0A_V0CR0", "EP A-R0 ", 100, 0., 100., 100, -1., 1.);
hdifV0A_V0CR3 = new TH2F("hdifV0A_V0CR3", "EP A-R3 ", 100, 0., 100., 100, -1., 1.);
hdifV0ACR0_V0CR3 = new TH2F("hdifV0ACR0_V0CR3", "EP R0-R3 ", 100, 0., 100., 100, -1., 1.);
hdifV0C_V0AR4 = new TH2F("hdifV0C_V0AR4", "EP C-R4 ", 100, 0., 100., 100, -1., 1.);
hdifV0C_V0AR7 = new TH2F("hdifV0C_V0AR7", "EP C-R7 ", 100, 0., 100., 100, -1., 1.);
hdifV0AR4_V0AR7 = new TH2F("hdifV0AR4_V0AR7", "EP R4-R7 ", 100, 0., 100., 100, -1., 1.);
- fOutput->Add(hdifV0Ar_V0Cr);
fOutput->Add(hdifV0A_V0CR0);
fOutput->Add(hdifV0A_V0CR3);
fOutput->Add(hdifV0ACR0_V0CR3);
fOutput->Add(hdifV0C_V0AR7);
fOutput->Add(hdifV0AR4_V0AR7);
- hdifV0A_V0C = new TH2F("hdifV0A_V0C", "EP A-C ", 100, 0., 100., 100, -1., 1.);
- hdifV0A_TPC = new TH2F("hdifV0A_TPC", "EP A-TPC", 100, 0., 100., 100, -1., 1.);
- hdifV0C_TPC = new TH2F("hdifV0C_TPC", "EP C-TPC", 100, 0., 100., 100, -1., 1.);
- hdifV0C_V0A = new TH2F("hdifV0C_V0A", "EP C-A ", 100, 0., 100., 100, -1., 1.);
- fOutput->Add(hdifV0A_V0C);
- fOutput->Add(hdifV0A_TPC);
- fOutput->Add(hdifV0C_TPC);
- fOutput->Add(hdifV0C_V0A);
-
- hdifEMC_EPV0 = new TH3F("hdifEMC_EPV0", "dif phi in EMC with EP", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
- hdifEMC_EPV0A = new TH3F("hdifEMC_EPV0A", "dif phi in EMC with EP", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
- hdifEMC_EPV0C = new TH3F("hdifEMC_EPV0C", "dif phi in EMC with EP", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
- fOutput->Add(hdifEMC_EPV0);
+ if(isFullHist){
+ hdifV0Ar_V0Cr = new TH2F("hdifV0Ar_V0Cr", "EP Ar-Cr ", 100, 0., 100., 100, -1., 1.);
+ fOutput->Add(hdifV0Ar_V0Cr);
+
+ hdifV0A_V0C = new TH2F("hdifV0A_V0C", "EP A-C ", 100, 0., 100., 100, -1., 1.);
+ hdifV0A_TPC = new TH2F("hdifV0A_TPC", "EP A-TPC", 100, 0., 100., 100, -1., 1.);
+ hdifV0C_TPC = new TH2F("hdifV0C_TPC", "EP C-TPC", 100, 0., 100., 100, -1., 1.);
+ hdifV0C_V0A = new TH2F("hdifV0C_V0A", "EP C-A ", 100, 0., 100., 100, -1., 1.);
+ fOutput->Add(hdifV0A_V0C);
+ fOutput->Add(hdifV0A_TPC);
+ fOutput->Add(hdifV0C_TPC);
+ fOutput->Add(hdifV0C_V0A);
+ }
+
+ hdifEMC_EPV0A = new TH3F("hdifEMC_EPV0A", "dif phi in EMC with EPV0A", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
+ hdifEMC_EPV0C = new TH3F("hdifEMC_EPV0C", "dif phi in EMC with EPV0C", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
fOutput->Add(hdifEMC_EPV0A);
fOutput->Add(hdifEMC_EPV0C);
- hdifful_EPV0 = new TH3F("hdifful_EPV0", "dif phi in full with EP", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
- hdifful_EPV0A = new TH3F("hdifful_EPV0A", "dif phi in full with EP", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
- hdifful_EPV0C = new TH3F("hdifful_EPV0C", "dif phi in full with EP", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
- fOutput->Add(hdifful_EPV0);
+ hdifful_EPV0A = new TH3F("hdifful_EPV0A", "dif phi in full with EPV0A", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
+ hdifful_EPV0C = new TH3F("hdifful_EPV0C", "dif phi in full with EPV0C", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
fOutput->Add(hdifful_EPV0A);
fOutput->Add(hdifful_EPV0C);
- hdifout_EPV0 = new TH3F("hdifout_EPV0", "dif phi NOT in EMC with EP", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
- hdifout_EPV0A = new TH3F("hdifout_EPV0A", "dif phi NOT in EMC with EP", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
- hdifout_EPV0C = new TH3F("hdifout_EPV0C", "dif phi NOT in EMC with EP", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
- fOutput->Add(hdifout_EPV0);
+ hdifout_EPV0A = new TH3F("hdifout_EPV0A", "dif phi NOT in EMC with EPV0A", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
+ hdifout_EPV0C = new TH3F("hdifout_EPV0C", "dif phi NOT in EMC with EPV0C", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
fOutput->Add(hdifout_EPV0A);
fOutput->Add(hdifout_EPV0C);
+ hCv2EMC_EPV0A = new TH3F("hCv2EMC_EPV0A", " raw v2 of charged trc in EMC with V0A", 100, 0, 100, 50, -1., 1., 15, 0., 15.);
+ hCv2EMC_EPV0C = new TH3F("hCv2EMC_EPV0C", " raw v2 of charged trc in EMC with V0C", 100, 0, 100, 50, -1., 1., 15, 0., 15.);
+ fOutput->Add(hCv2EMC_EPV0A);
+ fOutput->Add(hCv2EMC_EPV0C);
+
+ hCv2ful_EPV0A = new TH3F("hCv2ful_EPV0A", " raw v2 of charged trc in ful with V0A", 100, 0, 100, 50, -1., 1., 15, 0., 15.);
+ hCv2ful_EPV0C = new TH3F("hCv2ful_EPV0C", " raw v2 of charged trc in ful with V0C", 100, 0, 100, 50, -1., 1., 15, 0., 15.);
+ fOutput->Add(hCv2ful_EPV0A);
+ fOutput->Add(hCv2ful_EPV0C);
+
+ hCv2out_EPV0A = new TH3F("hCv2out_EPV0A", " raw v2 of charged trc out with V0A", 100, 0, 100, 50, -1., 1., 15, 0., 15.);
+ hCv2out_EPV0C = new TH3F("hCv2out_EPV0C", " raw v2 of charged trc out with V0A", 100, 0, 100, 50, -1., 1., 15, 0., 15.);
+ fOutput->Add(hCv2out_EPV0A);
+ fOutput->Add(hCv2out_EPV0C);
+
+ hclusDif_EPV0A = new TH3F("hclusDif_EPV0A", "dif phi of clus with EP V0A", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
+ hclusDif_EPV0C = new TH3F("hclusDif_EPV0C", "dif phi of clus with EP V0C", 100, 0., 100., 100, 0., TMath::Pi(), 15, 0., 15.);
+ fOutput->Add(hclusDif_EPV0A);
+ fOutput->Add(hclusDif_EPV0C);
+
+ hclusv2_EPV0A = new TH3F("hclusv2_EPV0A", " raw v2 of clus in ful with V0A", 100, 0, 100, 50, -1., 1., 15, 0., 15.);
+ hclusv2_EPV0C = new TH3F("hclusv2_EPV0C", " raw v2 of clus in ful with V0C", 100, 0, 100, 50, -1., 1., 15, 0., 15.);
+ fOutput->Add(hclusv2_EPV0A);
+ fOutput->Add(hclusv2_EPV0C);
+
if (isV1Clus) {
- // Et M02 spdcent DeltaPhi Dr
- Int_t bins[5] = { 500, 350, 100, 100, 100}; // binning
- Double_t min[5] = { 0.0, 0.0, 0, 0.0, 0 }; // min x
- Double_t max[5] = { 50.0, 3.5, 100, TMath::Pi(), 0.1}; // max x
+ // Et M02 spdcent DeltaPhi
+ Int_t bins[4] = { 40, 350, 60, 100 }; // binning
+ Double_t min[4] = { 0.0, 0.0, 0, 0.0 }; // min x
+ Double_t max[4] = { 40.0, 3.5, 60, TMath::Pi()}; // max x
fClusterPbV0A = new THnSparseF("fClusterPbV0A","",5,bins,min,max);
fClusterPbV0A->GetAxis(0)->SetTitle("Transverse Energy [GeV]");
fClusterPbV0A->GetAxis(1)->SetTitle("M02");
fClusterPbV0A->GetAxis(2)->SetTitle("V0M Centrality");
fClusterPbV0A->GetAxis(3)->SetTitle("Delta(#phi) [rad]");
- fClusterPbV0A->GetAxis(4)->SetTitle("Dr");
fOutput->Add(fClusterPbV0A);
fClusterPbV0C = new THnSparseF("fClusterPbV0C","",5,bins,min,max);
fClusterPbV0C->GetAxis(1)->SetTitle("M02");
fClusterPbV0C->GetAxis(2)->SetTitle("V0M Centrality");
fClusterPbV0C->GetAxis(3)->SetTitle("Delta(#phi) [rad]");
- fClusterPbV0C->GetAxis(4)->SetTitle("Dr");
fOutput->Add(fClusterPbV0C);
fClusterPbTPC = new THnSparseF("fClusterPbTPC","",5,bins,min,max);
fClusterPbTPC->GetAxis(1)->SetTitle("M02");
fClusterPbTPC->GetAxis(2)->SetTitle("V0M Centrality");
fClusterPbTPC->GetAxis(3)->SetTitle("Delta(#phi) [rad]");
- fClusterPbTPC->GetAxis(4)->SetTitle("Dr");
fOutput->Add(fClusterPbTPC);
}
+ if(isFullHist){
+ h2DcosV0C = new TProfile("h2DcosV0C", "cos(Phi) V0r vs Run NUmber", 200, 0., 200.);
+ h2DsinV0C = new TProfile("h2DsinV0C", "sin(Phi) V0r vs Run NUmber", 200, 0., 200.);
+ h2DcosTPC = new TProfile("h2DcosTPC", "cos(Phi) V0r vs Run NUmber", 200, 0., 200.);
+ h2DsinTPC = new TProfile("h2DsinTPC", "sin(Phi) V0r vs Run NUmber", 200, 0., 200.);
+ fOutput->Add(h2DcosV0C);
+ fOutput->Add(h2DsinV0C);
+ fOutput->Add(h2DcosTPC);
+ fOutput->Add(h2DsinTPC);
+ }
+
h2DcosV0A = new TProfile("h2DcosV0A", "cos(Phi) V0r vs Run NUmber", 200, 0., 200.);
h2DsinV0A = new TProfile("h2DsinV0A", "sin(Phi) V0r vs Run NUmber", 200, 0., 200.);
- h2DcosV0C = new TProfile("h2DcosV0C", "cos(Phi) V0r vs Run NUmber", 200, 0., 200.);
- h2DsinV0C = new TProfile("h2DsinV0C", "sin(Phi) V0r vs Run NUmber", 200, 0., 200.);
- h2DcosTPC = new TProfile("h2DcosTPC", "cos(Phi) V0r vs Run NUmber", 200, 0., 200.);
- h2DsinTPC = new TProfile("h2DsinTPC", "sin(Phi) V0r vs Run NUmber", 200, 0., 200.);
fOutput->Add(h2DcosV0A);
fOutput->Add(h2DsinV0A);
- fOutput->Add(h2DcosV0C);
- fOutput->Add(h2DsinV0C);
- fOutput->Add(h2DcosTPC);
- fOutput->Add(h2DsinTPC);
if (isV1Clus) {
hM02vsPtA = new TH2F("hM02vsPtA", "M02 vs Et before cut", 5000, 0, 50, 400, 0, 4.);
fOutput->Add(hM02vsPtA);
fOutput->Add(hM02vsPtB);
}
- hClusDxDZA = new TH2F("hClusDxDZA", "clus Dx vs Dz", 1000, -1., 1., 1000, -1., 1);
- hClusDxDZB = new TH2F("hClusDxDZB", "clus Dx vs Dz", 1000, -1., 1., 1000, -1., 1);
- fOutput->Add(hClusDxDZA);
- fOutput->Add(hClusDxDZB);
-
+ if(isFullHist){
+ hClusDxDZA = new TH2F("hClusDxDZA", "clus Dx vs Dz", 1000, -1., 1., 1000, -1., 1);
+ hClusDxDZB = new TH2F("hClusDxDZB", "clus Dx vs Dz", 1000, -1., 1., 1000, -1., 1);
+ fOutput->Add(hClusDxDZA);
+ fOutput->Add(hClusDxDZB);
+ }
+
if (!isV1Clus) {
- const Int_t ndims = 5;
- Int_t nMgg=500, nPt=40, nCent=20, nDeltaPhi=315, ncos2phi=500;
- Int_t binsv1[ndims] = {nMgg, nPt, nCent, nDeltaPhi, ncos2phi};
- Double_t xmin[ndims] = { 0, 0., 0, 0., -1.};
- Double_t xmax[ndims] = { 0.5, 20., 100, 3.15, 1.};
+ const Int_t ndims = 4;
+ Int_t nMgg=500, nPt=40, nCent=20, nDeltaPhi=315, ncos2phi=200;
+ Int_t binsv1[ndims] = {nMgg, nPt, nCent, nDeltaPhi};
+ Double_t xmin[ndims] = { 0, 0., 0, 0. };
+ Double_t xmax[ndims] = { 0.5, 20., 100, 3.15 };
fHEPV0A = new THnSparseF("fHEPV0A", "Flow histogram EPV0A", ndims, binsv1, xmin, xmax);
fHEPV0C = new THnSparseF("fHEPV0C", "Flow histogram EPV0C", ndims, binsv1, xmin, xmax);
fHEPTPC = new THnSparseF("fHEPTPC", "Flow histogram EPTPC", ndims, binsv1, xmin, xmax);
fHEPV0A->GetAxis(1)->SetTitle("p_{T}[GeV]");
fHEPV0A->GetAxis(2)->SetTitle("centrality");
fHEPV0A->GetAxis(3)->SetTitle("#delta #phi");
- fHEPV0A->GetAxis(4)->SetTitle("cos(2*#delta #phi)");
fHEPV0C->GetAxis(0)->SetTitle("m_{#gamma#gamma} ");
fHEPV0C->GetAxis(1)->SetTitle("p_{T}[GeV]");
fHEPV0C->GetAxis(2)->SetTitle("centrality");
fHEPV0C->GetAxis(3)->SetTitle("#delta #phi");
- fHEPV0C->GetAxis(4)->SetTitle("cos(2*#delta #phi)");
fHEPTPC->GetAxis(0)->SetTitle("m_{#gamma#gamma} ");
fHEPTPC->GetAxis(1)->SetTitle("p_{T}[GeV]");
fHEPTPC->GetAxis(2)->SetTitle("centrality");
fHEPTPC->GetAxis(3)->SetTitle("#delta #phi");
- fHEPTPC->GetAxis(4)->SetTitle("cos(2*#delta #phi)");
fOutput->Add(fHEPV0A);
fOutput->Add(fHEPV0C);
fOutput->Add(fHEPTPC);
+
+ Int_t binsv2[ndims] = {nMgg, nPt, nCent, ncos2phi};
+ Double_t xmin2[ndims] = { 0, 0., 0, -1.};
+ Double_t xmax2[ndims] = { 0.5, 20., 100, 1.};
+ fHEPV0AM2 = new THnSparseF("fHEPV0AM2", "Flow histogram EPV0A M2", ndims, binsv2, xmin2, xmax2);
+ fHEPV0CM2 = new THnSparseF("fHEPV0CM2", "Flow histogram EPV0C M2", ndims, binsv2, xmin2, xmax2);
+ fHEPTPCM2 = new THnSparseF("fHEPTPCM2", "Flow histogram EPTPC M2", ndims, binsv2, xmin2, xmax2);
+ fHEPV0AM2->GetAxis(0)->SetTitle("m_{#gamma#gamma} ");
+ fHEPV0AM2->GetAxis(1)->SetTitle("p_{T}[GeV]");
+ fHEPV0AM2->GetAxis(2)->SetTitle("centrality");
+ fHEPV0AM2->GetAxis(3)->SetTitle("cos(2*#delta #phi)");
+ fHEPV0CM2->GetAxis(0)->SetTitle("m_{#gamma#gamma} ");
+ fHEPV0CM2->GetAxis(1)->SetTitle("p_{T}[GeV]");
+ fHEPV0CM2->GetAxis(2)->SetTitle("centrality");
+ fHEPV0CM2->GetAxis(3)->SetTitle("cos(2*#delta #phi)");
+ fHEPTPCM2->GetAxis(0)->SetTitle("m_{#gamma#gamma} ");
+ fHEPTPCM2->GetAxis(1)->SetTitle("p_{T}[GeV]");
+ fHEPTPCM2->GetAxis(2)->SetTitle("centrality");
+ fHEPTPCM2->GetAxis(3)->SetTitle("cos(2*#delta #phi)");
+ fOutput->Add(fHEPV0AM2);
+ fOutput->Add(fHEPV0CM2);
+ fOutput->Add(fHEPTPCM2);
+
}
PostData(1, fOutput); // Post data for ALL output slots >0 here, to get at least an empty histogram
}
hEvtCount->Fill(4);
fCentrality = event->GetCentrality()->GetCentralityPercentile("CL1"); //spd vertex
-
+ hCent->Fill(fCentrality);
+ if(isCentFlat){
+ Bool_t bIsNot = kFALSE;
+ if (fCentrality<=10){ //0-10%
+ TRandom3 *rndm = new TRandom3(0);
+ Double_t Nrndm = rndm->Uniform(0.,1.);
+ if(fCentrality<=1){
+ if(Nrndm > 0.77308) bIsNot = kTRUE;
+ } else if(1<fCentrality && fCentrality<=2) {
+ if(Nrndm > 0.75863) bIsNot = kTRUE;
+ } else if (2<fCentrality && fCentrality<=3){
+ if(Nrndm > 0.76365) bIsNot = kTRUE;
+ } else if (3<fCentrality && fCentrality<=4){
+ if(Nrndm > 0.76763) bIsNot = kTRUE;
+ } else if (4<fCentrality && fCentrality<=5){
+ if(Nrndm > 0.76251) bIsNot = kTRUE;
+ } else if (5<fCentrality && fCentrality<=6){
+ if(Nrndm > 0.79069) bIsNot = kTRUE;
+ } else if (6<fCentrality && fCentrality<=7){
+ if(Nrndm > 0.77669) bIsNot = kTRUE;
+ } else if (7<fCentrality && fCentrality<=8){
+ if(Nrndm > 0.78537) bIsNot = kTRUE;
+ } else if (8<fCentrality && fCentrality<=9){
+ if(Nrndm > 0.82727) bIsNot = kTRUE;
+ } else if (9<fCentrality && fCentrality<=10){
+ if(Nrndm > 1) bIsNot = kTRUE;
+ }
+ delete rndm; rndm = 0;
+ if(bIsNot)
+ return;
+ }
+ }
+ if (10<fCentrality && fCentrality<=50){ //10-50%
+ TString centfired;
+ if (fESD) {
+ centfired = fESD->GetFiredTriggerClasses();
+ } else {
+ centfired = fAOD->GetFiredTriggerClasses();
+ }
+ if(!centfired.Contains("CVLN_B2-B-NOPF-ALLNOTRD") && !centfired.Contains("CVLN_R1-B-NOPF-ALLNOTRD") && !centfired.Contains("CSEMI_R1-B-NOPF-ALLNOTRD"))
+ return;
+ }
hEvtCount->Fill(5);
AliEventplane *ep = event->GetEventplane();
hEvtCount->Fill(6);
- fEPV0 = TVector2::Phi_0_2pi(fEPV0);
- if (fEPV0>TMath::Pi())
- fEPV0 = fEPV0 - TMath::Pi();
- fEPV0r = TVector2::Phi_0_2pi(fEPV0r);
- if (fEPV0r>TMath::Pi())
- fEPV0r = fEPV0r - TMath::Pi();
- fEPV0A = TVector2::Phi_0_2pi(fEPV0A);
- if (fEPV0A>TMath::Pi())
- fEPV0A = fEPV0A - TMath::Pi();
- fEPV0C = TVector2::Phi_0_2pi(fEPV0C);
- if (fEPV0C>TMath::Pi())
- fEPV0C = fEPV0C - TMath::Pi();
- fEPV0Ar = TVector2::Phi_0_2pi(fEPV0Ar);
- if (fEPV0Ar>TMath::Pi())
- fEPV0Ar = fEPV0Ar - TMath::Pi();
- fEPV0Cr = TVector2::Phi_0_2pi(fEPV0Cr);
- if (fEPV0Cr>TMath::Pi())
- fEPV0Cr = fEPV0Cr - TMath::Pi();
- fEPV0AR4 = TVector2::Phi_0_2pi(fEPV0AR4);
- if (fEPV0AR4>TMath::Pi())
- fEPV0AR4 = fEPV0AR4 - TMath::Pi();
- fEPV0AR7 = TVector2::Phi_0_2pi(fEPV0AR7);
- if (fEPV0AR7>TMath::Pi())
- fEPV0AR7 = fEPV0AR7 - TMath::Pi();
- fEPV0CR0 = TVector2::Phi_0_2pi(fEPV0CR0);
- if (fEPV0CR0>TMath::Pi())
- fEPV0CR0 = fEPV0CR0 - TMath::Pi();
- fEPV0CR3 = TVector2::Phi_0_2pi(fEPV0CR3);
- if (fEPV0CR3>TMath::Pi())
- fEPV0CR3 = fEPV0CR3 - TMath::Pi();
+ while(fEPV0<0.) fEPV0+=TMath::Pi(); while(fEPV0>TMath::Pi()) fEPV0-=TMath::Pi();
+ while(fEPV0r<0.) fEPV0r+=TMath::Pi(); while(fEPV0r>TMath::Pi()) fEPV0r-=TMath::Pi();
+ while(fEPV0A<0.) fEPV0A+=TMath::Pi(); while(fEPV0A>TMath::Pi()) fEPV0A-=TMath::Pi();
+ while(fEPV0C<0.) fEPV0C+=TMath::Pi(); while(fEPV0C>TMath::Pi()) fEPV0C-=TMath::Pi();
+ while(fEPV0Ar<0.) fEPV0Ar+=TMath::Pi(); while(fEPV0Ar>TMath::Pi()) fEPV0Ar-=TMath::Pi();
+ while(fEPV0Cr<0.) fEPV0Cr+=TMath::Pi(); while(fEPV0Cr>TMath::Pi()) fEPV0Cr-=TMath::Pi();
+
+ while(fEPV0AR4<0.) fEPV0AR4+=TMath::Pi(); while(fEPV0AR4>TMath::Pi()) fEPV0AR4-=TMath::Pi();
+ while(fEPV0AR7<0.) fEPV0AR7+=TMath::Pi(); while(fEPV0AR7>TMath::Pi()) fEPV0AR7-=TMath::Pi();
+ while(fEPV0CR0<0.) fEPV0CR0+=TMath::Pi(); while(fEPV0CR0>TMath::Pi()) fEPV0CR0-=TMath::Pi();
+ while(fEPV0CR3<0.) fEPV0CR3+=TMath::Pi(); while(fEPV0CR3>TMath::Pi()) fEPV0CR3-=TMath::Pi();
+
+ while(fEPTPC<0.) fEPTPC+=TMath::Pi(); while(fEPTPC>TMath::Pi()) fEPTPC-=TMath::Pi();
if (fEPTPC != -999. )
hEPTPC->Fill(fCentrality, fEPTPC);
if (fEPTPCreso!=-1)
hresoTPC->Fill(fCentrality, fEPTPCreso);
- hEPV0->Fill(fCentrality, fEPV0);
+ if(isFullHist){
+ hEPV0->Fill(fCentrality, fEPV0);
+ hEPV0Ar->Fill(fCentrality, fEPV0Ar);
+ hEPV0Cr->Fill(fCentrality, fEPV0Cr);
+ hEPV0r->Fill(fCentrality, fEPV0r);
+ hEPV0AR4->Fill(fCentrality, fEPV0AR4);
+ hEPV0AR7->Fill(fCentrality, fEPV0AR7);
+ hEPV0CR0->Fill(fCentrality, fEPV0CR0);
+ hEPV0CR3->Fill(fCentrality, fEPV0CR3);
+ }
hEPV0A->Fill(fCentrality, fEPV0A);
hEPV0C->Fill(fCentrality, fEPV0C);
- hEPV0Ar->Fill(fCentrality, fEPV0Ar);
- hEPV0Cr->Fill(fCentrality, fEPV0Cr);
- hEPV0r->Fill(fCentrality, fEPV0r);
- hEPV0AR4->Fill(fCentrality, fEPV0AR4);
- hEPV0AR7->Fill(fCentrality, fEPV0AR7);
- hEPV0CR0->Fill(fCentrality, fEPV0CR0);
- hEPV0CR3->Fill(fCentrality, fEPV0CR3);
if (isPhosCali) {
// PHOS Flattening
fEPV0C = ApplyFlatteningV0C(fEPV0C, fCentrality); //V0C after Phos flatten
if(fEPTPC != -999.)
fEPTPC = ApplyFlattening(fEPTPC, fCentrality); //TPC after Phos flatten
+ while(fEPV0A <0.) fEPV0A+=TMath::Pi(); while(fEPV0A >TMath::Pi()) fEPV0A-=TMath::Pi();
+ while(fEPV0C <0.) fEPV0C+=TMath::Pi(); while(fEPV0C >TMath::Pi()) fEPV0C-=TMath::Pi();
+ while(fEPTPC <0.) fEPTPC+=TMath::Pi(); while(fEPTPC >TMath::Pi()) fEPTPC-=TMath::Pi();
}
if (!isPhosCali) {
+ Double_t EPV0ACor = ApplyFlattening(fEPTPC, fCentrality);
+ Double_t EPV0CCor = ApplyFlattening(fEPTPC, fCentrality);
+ Double_t EPTPCCor = ApplyFlattening(fEPTPC, fCentrality);
+ while(EPV0ACor <0.) EPV0ACor+=TMath::Pi(); while(EPV0ACor >TMath::Pi()) EPV0ACor-=TMath::Pi();
+ while(EPV0CCor <0.) EPV0CCor+=TMath::Pi(); while(EPV0CCor >TMath::Pi()) EPV0CCor-=TMath::Pi();
+ while(EPTPCCor <0.) EPTPCCor+=TMath::Pi(); while(EPTPCCor >TMath::Pi()) EPTPCCor-=TMath::Pi();
+
if(fEPTPC != -999.)
- hEPTPCCor->Fill(fCentrality, ApplyFlattening(fEPTPC, fCentrality));
- hEPV0ACor->Fill(fCentrality, ApplyFlatteningV0A(fEPV0A, fCentrality));
- hEPV0CCor->Fill(fCentrality, ApplyFlatteningV0C(fEPV0C, fCentrality));
+ hEPTPCCor->Fill(fCentrality, EPTPCCor);
+ hEPV0ACor->Fill(fCentrality, EPV0ACor);
+ hEPV0CCor->Fill(fCentrality, EPV0CCor);
} else {
if(fEPTPC != -999.)
hEPTPCCor->Fill(fCentrality, fEPTPC);
hEPV0CCor->Fill(fCentrality, fEPV0C);
}
- hdifV0Ar_V0Cr->Fill(fCentrality, TMath::Cos(2.*(fEPV0Ar - fEPV0Cr)));
hdifV0A_V0CR0->Fill(fCentrality, TMath::Cos(2.*(fEPV0A - fEPV0CR0)));
hdifV0A_V0CR3->Fill(fCentrality, TMath::Cos(2.*(fEPV0A - fEPV0CR3)));
hdifV0ACR0_V0CR3->Fill(fCentrality, TMath::Cos(2*(fEPV0CR0 - fEPV0CR3)));
hdifV0C_V0AR7->Fill(fCentrality, TMath::Cos(2*(fEPV0C - fEPV0AR7)));
hdifV0AR4_V0AR7->Fill(fCentrality, TMath::Cos(2*(fEPV0AR4 - fEPV0AR7)));
- hdifV0A_V0C->Fill(fCentrality, TMath::Cos(2*(fEPV0A - fEPV0C)));
- if (fEPTPC!=-999.){
- hdifV0A_TPC->Fill(fCentrality, TMath::Cos(2*(fEPV0A - fEPTPC)));
- hdifV0C_TPC->Fill(fCentrality, TMath::Cos(2*(fEPV0C - fEPTPC)));
+ if(isFullHist){
+ hdifV0Ar_V0Cr->Fill(fCentrality, TMath::Cos(2.*(fEPV0Ar - fEPV0Cr)));
+ hdifV0A_V0C->Fill(fCentrality, TMath::Cos(2*(fEPV0A - fEPV0C)));
+ if (fEPTPC!=-999.){
+ hdifV0A_TPC->Fill(fCentrality, TMath::Cos(2*(fEPV0A - fEPTPC)));
+ hdifV0C_TPC->Fill(fCentrality, TMath::Cos(2*(fEPV0C - fEPTPC)));
+ }
+ hdifV0C_V0A->Fill(fCentrality, TMath::Cos(2*(fEPV0C - fEPV0A)));
}
- hdifV0C_V0A->Fill(fCentrality, TMath::Cos(2*(fEPV0C - fEPV0A)));
-
// Cluster loop for reconstructed event
//================ for v2 clusterize analysis==============================================
AliVCluster *c1 = static_cast<AliVCluster*>(fV2Clus->At(i));
if (!c1)
continue;
- hClusDxDZA->Fill(c1->GetTrackDz(), c1->GetTrackDx());
+ if(isFullHist) hClusDxDZA->Fill(c1->GetTrackDz(), c1->GetTrackDx());
if (!c1->IsEMCAL())
continue;
if (!IsGoodCluster(c1))
continue;
- hClusDxDZB->Fill(c1->GetTrackDz(), c1->GetTrackDx());
+ if(isFullHist) hClusDxDZB->Fill(c1->GetTrackDz(), c1->GetTrackDx());
TLorentzVector p1;
GetMom(p1, c1, vertex);
+ Double_t cluPhi = p1.Phi();
+ Double_t cluPt = p1.Pt();
+ Double_t difclusV0A = cluPhi-fEPV0A;
+ if(difclusV0A<0.) difclusV0A+=TMath::Pi(); if(difclusV0A>TMath::Pi()) difclusV0A-=TMath::Pi();
+ Double_t difclusV0C = cluPhi-fEPV0C;
+ if(difclusV0C<0.) difclusV0C+=TMath::Pi(); if(difclusV0C>TMath::Pi()) difclusV0C-=TMath::Pi();
+ hclusDif_EPV0A->Fill(fCentrality, difclusV0A, cluPt);
+ hclusDif_EPV0C->Fill(fCentrality, difclusV0C, cluPt);
+ hclusv2_EPV0A->Fill(fCentrality, TMath::Cos(2.*difclusV0A), cluPt);
+ hclusv2_EPV0C->Fill(fCentrality, TMath::Cos(2.*difclusV0C), cluPt);
for (Int_t j=i+1; j<nCluster; ++j) {
AliVCluster *c2 = static_cast<AliVCluster*>(fV2Clus->At(j));
if (!c2)
Double_t Dxc3 = c3->GetTrackDx();
Double_t Dzc3 = c3->GetTrackDz();
- hClusDxDZA->Fill(Dzc3, Dxc3);
+ if(isFullHist) hClusDxDZA->Fill(Dzc3, Dxc3);
Float_t clsPosEt[3] = {0,0,0};
c3->GetPosition(clsPosEt);
TVector3 clsVec(clsPosEt);
if (!IsGoodClusterV1(c3))
continue;
hM02vsPtB->Fill(Et, M02c3);
- hClusDxDZB->Fill(Dzc3, Dxc3);
+ if(isFullHist) hClusDxDZB->Fill(Dzc3, Dxc3);
TLorentzVector p3;
GetMom(p3, c3, vertex);
FillCluster(p3, fEPV0A, fEPV0C, fEPTPC, c3);
Double_t tPt = track->Pt();
Double_t Eta = track->Eta();
- Double_t difTrackV0 = TVector2::Phi_0_2pi(tPhi-fEPV0);
- if (difTrackV0 >TMath::Pi())
- difTrackV0 -= TMath::Pi();
- Double_t difTrackV0A = TVector2::Phi_0_2pi(tPhi-fEPV0A);
- if (difTrackV0A >TMath::Pi())
- difTrackV0A -= TMath::Pi();
- Double_t difTrackV0C = TVector2::Phi_0_2pi(tPhi-fEPV0C);
- if (difTrackV0C >TMath::Pi())
- difTrackV0C -= TMath::Pi();
- Double_t difTrackTPC = TVector2::Phi_0_2pi(tPhi-fEPTPC);
- if (difTrackTPC >TMath::Pi())
- difTrackTPC -= TMath::Pi();
+ Double_t difTrackV0 = tPhi-fEPV0;
+ while(difTrackV0 <0.) difTrackV0+=TMath::Pi(); while(difTrackV0 >TMath::Pi()) difTrackV0-=TMath::Pi();
+ Double_t difTrackV0A = tPhi-fEPV0A;
+ while(difTrackV0A <0.) difTrackV0A+=TMath::Pi(); while(difTrackV0A >TMath::Pi()) difTrackV0A-=TMath::Pi();
+ Double_t difTrackV0C = tPhi-fEPV0C;
+ while(difTrackV0C <0.) difTrackV0C+=TMath::Pi(); while(difTrackV0C >TMath::Pi()) difTrackV0C-=TMath::Pi();
+ Double_t difTrackTPC = tPhi-fEPTPC;
+ while(difTrackTPC <0.) difTrackTPC+=TMath::Pi(); while(difTrackTPC >TMath::Pi()) difTrackTPC-=TMath::Pi();
if (tPhi*TMath::RadToDeg()>80. && tPhi*TMath::RadToDeg()<180. && Eta <0.7 && Eta >(-0.7)){
- hdifEMC_EPV0->Fill(fCentrality, difTrackV0, tPt);
hdifEMC_EPV0A->Fill(fCentrality, difTrackV0A, tPt);
hdifEMC_EPV0C->Fill(fCentrality, difTrackV0C, tPt);
+ hCv2EMC_EPV0A->Fill(fCentrality, TMath::Cos(2.*difTrackV0A), tPt);
+ hCv2EMC_EPV0C->Fill(fCentrality, TMath::Cos(2.*difTrackV0C), tPt);
} else {
- hdifout_EPV0->Fill(fCentrality, difTrackV0, tPt);
hdifout_EPV0A->Fill(fCentrality, difTrackV0A, tPt);
hdifout_EPV0C->Fill(fCentrality, difTrackV0C, tPt);
+ hCv2out_EPV0A->Fill(fCentrality, TMath::Cos(2.*difTrackV0A), tPt);
+ hCv2out_EPV0C->Fill(fCentrality, TMath::Cos(2.*difTrackV0C), tPt);
}
- hdifful_EPV0->Fill(fCentrality, difTrackV0, tPt);
hdifful_EPV0A->Fill(fCentrality, difTrackV0A, tPt);
hdifful_EPV0C->Fill(fCentrality, difTrackV0C, tPt);
+ hCv2ful_EPV0A->Fill(fCentrality, TMath::Cos(2.*difTrackV0A), tPt);
+ hCv2ful_EPV0C->Fill(fCentrality, TMath::Cos(2.*difTrackV0C), tPt);
}
hEvtCount->Fill(8);
{
h2DcosV0A->Fill(fInterRunNumber, TMath::Cos(fEPV0A));
h2DsinV0A->Fill(fInterRunNumber, TMath::Sin(fEPV0A));
- h2DcosV0C->Fill(fInterRunNumber, TMath::Cos(fEPV0C));
- h2DsinV0C->Fill(fInterRunNumber, TMath::Sin(fEPV0C));
- if (fEPTPC!=-999.){
- h2DcosTPC->Fill(fInterRunNumber, TMath::Cos(fEPTPC));
- h2DsinTPC->Fill(fInterRunNumber, TMath::Sin(fEPTPC));
+ if(isFullHist){
+ h2DcosV0C->Fill(fInterRunNumber, TMath::Cos(fEPV0C));
+ h2DsinV0C->Fill(fInterRunNumber, TMath::Sin(fEPV0C));
+ if (fEPTPC!=-999.){
+ h2DcosTPC->Fill(fInterRunNumber, TMath::Cos(fEPTPC));
+ h2DsinTPC->Fill(fInterRunNumber, TMath::Sin(fEPTPC));
+ }
}
}