fCSEMIbit(0),
fCentV0M(0),
fCentV0A(0),
+ fCentV0A123(0),
fCentV0C(0),
fCentV0MEq(0),
fCentV0AEq(0),
fCentCND(0),
fCentZNA(0),
fCentZNC(0),
+ fCentZPA(0),
+ fCentZPC(0),
fCentNPA(0),
fCentV0MvsFMD(0),
fCentTKLvsV0M(0),
fCentCNDtrue(0),
fCentZNAtrue(0),
fCentZNCtrue(0),
+ fCentZPAtrue(0),
+ fCentZPCtrue(0),
fHtempV0M(0),
fHtempV0A(0),
+ fHtempV0A123(0),
fHtempV0C(0),
fHtempV0MEq(0),
fHtempV0AEq(0),
fHtempCND(0),
fHtempZNA(0),
fHtempZNC(0),
+ fHtempZPA(0),
+ fHtempZPC(0),
fHtempV0MvsFMD(0),
fHtempTKLvsV0M(0),
fHtempZEMvsZDC(0),
fHtempCNDtrue(0),
fHtempZNAtrue(0),
fHtempZNCtrue(0),
+ fHtempZPAtrue(0),
+ fHtempZPCtrue(0),
fOutputList(0),
fHOutCentV0M(0),
fHOutCentV0A(0),
+ fHOutCentV0A123(0),
fHOutCentV0C(0),
fHOutCentV0MEq(0),
fHOutCentV0AEq(0),
fHOutCentNPA(0),
fHOutCentZNA(0),
fHOutCentZNC(0),
+ fHOutCentZPA(0),
+ fHOutCentZPC(0),
fHOutCentV0MvsFMD(0),
fHOutCentTKLvsV0M(0),
fHOutCentZEMvsZDC(0),
fHOutCentZNAvsCentTRK(0),
fHOutCentZNAvsCentCND(0),
fHOutCentZNAvsCentCL1(0),
+ fHOutCentZNAvsCentZPA(0),
fHOutMultV0AC(0),
fHOutMultV0M(0),
fHOutMultV0A(0),
+ fHOutMultV0A123(0),
fHOutMultV0C(0),
fHOutMultV0MEq(0),
fHOutMultV0AEq(0),
fHOutMultNPA(0),
fHOutMultZNA(0),
fHOutMultZNC(0),
+ fHOutMultZPA(0),
+ fHOutMultZPC(0),
fHOutMultV0MvsZDN(0),
fHOutMultZEMvsZDN(0),
fHOutMultV0MvsZDC(0),
fHOutMultV0OvsCL1(0),
fHOutMultV0OvsTRK(0),
fHOutMultCL1vsTKL(0),
+ fHOutMultZNAvsZPA(0),
fHOutCentV0Mqual1(0),
fHOutCentTRKqual1(0),
fHOutCentCL1qual1(0),
fCSEMIbit(0),
fCentV0M(0),
fCentV0A(0),
+ fCentV0A123(0),
fCentV0C(0),
fCentV0MEq(0),
fCentV0AEq(0),
fCentCND(0),
fCentZNA(0),
fCentZNC(0),
+ fCentZPA(0),
+ fCentZPC(0),
fCentNPA(0),
fCentV0MvsFMD(0),
fCentTKLvsV0M(0),
fCentCNDtrue(0),
fCentZNAtrue(0),
fCentZNCtrue(0),
+ fCentZPAtrue(0),
+ fCentZPCtrue(0),
fHtempV0M(0),
fHtempV0A(0),
+ fHtempV0A123(0),
fHtempV0C(0),
fHtempV0MEq(0),
fHtempV0AEq(0),
fHtempCND(0),
fHtempZNA(0),
fHtempZNC(0),
+ fHtempZPA(0),
+ fHtempZPC(0),
fHtempV0MvsFMD(0),
fHtempTKLvsV0M(0),
fHtempZEMvsZDC(0),
fHtempCNDtrue(0),
fHtempZNAtrue(0),
fHtempZNCtrue(0),
+ fHtempZPAtrue(0),
+ fHtempZPCtrue(0),
fOutputList(0),
fHOutCentV0M(0),
fHOutCentV0A(0),
+ fHOutCentV0A123(0),
fHOutCentV0C(0),
fHOutCentV0MEq(0),
fHOutCentV0AEq(0),
fHOutCentNPA(0),
fHOutCentZNA(0),
fHOutCentZNC(0),
+ fHOutCentZPA(0),
+ fHOutCentZPC(0),
fHOutCentV0MvsFMD(0),
fHOutCentTKLvsV0M(0),
fHOutCentZEMvsZDC(0),
fHOutCentZNAvsCentTRK(0),
fHOutCentZNAvsCentCND(0),
fHOutCentZNAvsCentCL1(0),
+ fHOutCentZNAvsCentZPA(0),
fHOutMultV0AC(0),
fHOutMultV0M(0),
fHOutMultV0A(0),
+ fHOutMultV0A123(0),
fHOutMultV0C(0),
fHOutMultV0MEq(0),
fHOutMultV0AEq(0),
fHOutMultNPA(0),
fHOutMultZNA(0),
fHOutMultZNC(0),
+ fHOutMultZPA(0),
+ fHOutMultZPC(0),
fHOutMultV0MvsZDN(0),
fHOutMultZEMvsZDN(0),
fHOutMultV0MvsZDC(0),
fHOutMultV0OvsCL1(0),
fHOutMultV0OvsTRK(0),
fHOutMultCL1vsTKL(0),
+ fHOutMultZNAvsZPA(0),
fHOutCentV0Mqual1(0),
fHOutCentTRKqual1(0),
fHOutCentCL1qual1(0),
fCSEMIbit(ana.fCSEMIbit),
fCentV0M(ana.fCentV0M),
fCentV0A(ana.fCentV0A),
+ fCentV0A123(ana.fCentV0A123),
fCentV0C(ana.fCentV0C),
fCentV0MEq(ana.fCentV0MEq),
fCentV0AEq(ana.fCentV0AEq),
fCentCND(ana.fCentCND),
fCentZNA(ana.fCentZNA),
fCentZNC(ana.fCentZNC),
+ fCentZPA(ana.fCentZPA),
+ fCentZPC(ana.fCentZPC),
fCentNPA(ana.fCentNPA),
fCentV0MvsFMD(ana.fCentV0MvsFMD),
fCentTKLvsV0M(ana.fCentTKLvsV0M),
fCentCNDtrue(ana.fCentCNDtrue),
fCentZNAtrue(ana.fCentZNAtrue),
fCentZNCtrue(ana.fCentZNCtrue),
+ fCentZPAtrue(ana.fCentZPAtrue),
+ fCentZPCtrue(ana.fCentZPCtrue),
fHtempV0M(ana.fHtempV0M),
fHtempV0A(ana.fHtempV0A),
+ fHtempV0A123(ana.fHtempV0A123),
fHtempV0C(ana.fHtempV0C),
fHtempV0MEq(ana.fHtempV0MEq),
fHtempV0AEq(ana.fHtempV0AEq),
fHtempCND(ana.fHtempCND),
fHtempZNA(ana.fHtempZNA),
fHtempZNC(ana.fHtempZNC),
+ fHtempZPA(ana.fHtempZPA),
+ fHtempZPC(ana.fHtempZPC),
fHtempV0MvsFMD(ana.fHtempV0MvsFMD),
fHtempTKLvsV0M(ana.fHtempTKLvsV0M),
fHtempZEMvsZDC(ana.fHtempZEMvsZDC),
fHtempCNDtrue(ana.fHtempCNDtrue),
fHtempZNAtrue(ana.fHtempZNAtrue),
fHtempZNCtrue(ana.fHtempZNCtrue),
+ fHtempZPAtrue(ana.fHtempZPAtrue),
+ fHtempZPCtrue(ana.fHtempZPCtrue),
fOutputList(ana.fOutputList),
fHOutCentV0M(ana.fHOutCentV0M),
fHOutCentV0A(ana.fHOutCentV0A),
+ fHOutCentV0A123(ana.fHOutCentV0A123),
fHOutCentV0C(ana.fHOutCentV0C),
fHOutCentV0MEq(ana.fHOutCentV0MEq),
fHOutCentV0AEq(ana.fHOutCentV0AEq),
fHOutCentNPA(ana.fHOutCentNPA),
fHOutCentZNA(ana.fHOutCentZNA),
fHOutCentZNC(ana.fHOutCentZNC),
+ fHOutCentZPA(ana.fHOutCentZPA),
+ fHOutCentZPC(ana.fHOutCentZPC),
fHOutCentV0MvsFMD(ana.fHOutCentV0MvsFMD),
fHOutCentTKLvsV0M(ana.fHOutCentTKLvsV0M),
fHOutCentZEMvsZDC(ana.fHOutCentZEMvsZDC),
fHOutCentZNAvsCentTRK(ana.fHOutCentZNAvsCentTRK),
fHOutCentZNAvsCentCND(ana.fHOutCentZNAvsCentCND),
fHOutCentZNAvsCentCL1(ana.fHOutCentZNAvsCentCL1),
+ fHOutCentZNAvsCentZPA(ana.fHOutCentZNAvsCentZPA),
fHOutMultV0AC(ana.fHOutMultV0AC),
fHOutMultV0M(ana.fHOutMultV0M),
fHOutMultV0A(ana.fHOutMultV0A),
+ fHOutMultV0A123(ana.fHOutMultV0A123),
fHOutMultV0C(ana.fHOutMultV0C),
fHOutMultV0MEq(ana.fHOutMultV0MEq),
fHOutMultV0AEq(ana.fHOutMultV0AEq),
fHOutMultNPA(ana.fHOutMultNPA),
fHOutMultZNA(ana.fHOutMultZNA),
fHOutMultZNC(ana.fHOutMultZNC),
+ fHOutMultZPA(ana.fHOutMultZPA),
+ fHOutMultZPC(ana.fHOutMultZPC),
fHOutMultV0MvsZDN(ana.fHOutMultV0MvsZDN),
fHOutMultZEMvsZDN(ana.fHOutMultZEMvsZDN),
fHOutMultV0MvsZDC(ana.fHOutMultV0MvsZDC),
fHOutMultV0OvsCL1(ana.fHOutMultV0OvsCL1),
fHOutMultV0OvsTRK(ana.fHOutMultV0OvsTRK),
fHOutMultCL1vsTKL(ana.fHOutMultCL1vsTKL),
+ fHOutMultZNAvsZPA(ana.fHOutMultZNAvsZPA),
fHOutCentV0Mqual1(ana.fHOutCentV0Mqual1),
fHOutCentTRKqual1(ana.fHOutCentTRKqual1),
fHOutCentCL1qual1(ana.fHOutCentCL1qual1),
fOutputList->SetOwner();
fHOutCentV0M = new TH1F("fHOutCentV0M","fHOutCentV0M; Centrality V0",505,0,101);
fHOutCentV0A = new TH1F("fHOutCentV0A","fHOutCentV0A; Centrality V0A",505,0,101);
+ fHOutCentV0A123 = new TH1F("fHOutCentV0A123","fHOutCentV0A123; Centrality V0A-123",505,0,101);
fHOutCentV0C = new TH1F("fHOutCentV0C","fHOutCentV0C; Centrality V0C",505,0,101);
fHOutCentV0MEq = new TH1F("fHOutCentV0MEq","fHOutCentV0MEq; Centrality V0 equalized",505,0,101);
fHOutCentV0AEq = new TH1F("fHOutCentV0AEq","fHOutCentV0AEq; Centrality V0A equalized",505,0,101);
fHOutCentNPA = new TH1F("fHOutCentNPA","fHOutCentNPA; Centrality Npart",505,0,101);
fHOutCentZNA = new TH1F("fHOutCentZNA","fHOutCentZNA; Centrality ZNA",505,0,101);
fHOutCentZNC = new TH1F("fHOutCentZNC","fHOutCentZNC; Centrality ZNC",505,0,101);
+ fHOutCentZPA = new TH1F("fHOutCentZPA","fHOutCentZPA; Centrality ZPA",505,0,101);
+ fHOutCentZPC = new TH1F("fHOutCentZPC","fHOutCentZPC; Centrality ZPC",505,0,101);
fHOutCentV0MvsFMD= new TH1F("fHOutCentV0MvsFMD","fHOutCentV0MvsFMD; Centrality V0 vs FMD",505,0,101);
fHOutCentTKLvsV0M= new TH1F("fHOutCentTKLvsV0M","fHOutCentTKLvsV0M; Centrality tracklets vs V0",505,0,101);
fHOutCentZEMvsZDC= new TH1F("fHOutCentZEMvsZDC","fHOutCentZEMvsZDC; Centrality ZEM vs ZDC",505,0,101);
fHOutCentZNAvsCentTRK= new TH2F("fHOutCentZNAvsCentTRK","fHOutCentZNAvsCentTRK; Cent ZNA; Cent TRK;", 505,0,101,505,0,101);
fHOutCentZNAvsCentCND= new TH2F("fHOutCentZNAvsCentCND","fHOutCentZNAvsCentCND; Cent ZNA; Cent CND;", 505,0,101,505,0,101);
fHOutCentZNAvsCentCL1= new TH2F("fHOutCentZNAvsCentCL1","fHOutCentZNAvsCentCL1; Cent ZNA; Cent CL1;", 505,0,101,505,0,101);
+ fHOutCentZNAvsCentZPA= new TH2F("fHOutCentZNAvsCentZPA","fHOutCentZNAvsCentZPA; Cent ZNA; Cent ZPA;", 505,0,101,505,0,101);
fHOutMultV0AC = new TH2F("fHOutMultV0AC","fHOutMultV0AC; Multiplicity V0A; Multiplicity V0C",1000,0,1000,1000,0,1000);
fHOutMultV0M = new TH1F("fHOutMultV0M","fHOutMultV0M; Multiplicity V0",25000,0,25000);
fHOutMultV0A = new TH1F("fHOutMultV0A","fHOutMultV0A; Multiplicity V0",25000,0,25000);
+ fHOutMultV0A123 = new TH1F("fHOutMultV0A123","fHOutMultV0A123; Multiplicity V0A-123",25000,0,25000);
fHOutMultV0C = new TH1F("fHOutMultV0C","fHOutMultV0C; Multiplicity V0",25000,0,25000);
fHOutMultV0MEq = new TH1F("fHOutMultV0MEq","fHOutMultV0MEq; Multiplicity V0",25000,0,25000);
fHOutMultV0AEq = new TH1F("fHOutMultV0AEq","fHOutMultV0AEq; Multiplicity V0",25000,0,25000);
fHOutMultCL1 = new TH1F("fHOutMultCL1","fHOutMultCL1; Multiplicity SPD outer",7000,0,7000);
fHOutMultCND = new TH1F("fHOutMultCND","fHOutMultCND; Multiplicity candle",4000,0,4000);
fHOutMultNPA = new TH1F("fHOutMultNPA","fHOutMultNPA; Nparticipants",450,0,450);
- fHOutMultZNA = new TH1F("fHOutMultZNA","fHOutMultZNA; ZNA Energy",500,0,2000);
- fHOutMultZNC = new TH1F("fHOutMultZNC","fHOutMultZNC; ZNC Energy",500,0,2000);
+ fHOutMultZNA = new TH1F("fHOutMultZNA","fHOutMultZNA; ZNA Energy",2000,0,2000);
+ fHOutMultZNC = new TH1F("fHOutMultZNC","fHOutMultZNC; ZNC Energy",2000,0,2000);
+ fHOutMultZPA = new TH1F("fHOutMultZPA","fHOutMultZPA; ZPA Energy",2000,0,2000);
+ fHOutMultZPC = new TH1F("fHOutMultZPC","fHOutMultZPC; ZPC Energy",2000,0,2000);
- fHOutMultV0MvsZDN = new TH2F("fHOutMultV0MvsZDN","fHOutMultV0MvsZDN; Multiplicity V0; Energy ZDC-N",500,0,30000,500,0,180000);
+ fHOutMultV0MvsZDN = new TH2F("fHOutMultV0MvsZDN","fHOutMultV0MvsZDN; Multiplicity V0; Energy ZDC-N",500,0,1200,500,0,2000);
fHOutMultZEMvsZDN = new TH2F("fHOutMultZEMvsZDN","fHOutMultZEMvsZDN; Energy ZEM; Energy ZDC-N",500,0,2500,500,0,180000);
fHOutMultV0MvsZDC = new TH2F("fHOutMultV0MvsZDC","fHOutMultV0MvsZDC; Multiplicity V0; Energy ZDC",500,0,30000,500,0,200000);
fHOutMultZEMvsZDC = new TH2F("fHOutMultZEMvsZDC","fHOutMultZEMvsZDC; Energy ZEM; Energy ZDC",500,0,2500,500,0,200000);
fHOutMultV0OvsCL1 = new TH2F("fHOutMultV0OvsCL1","fHOutMultV0OvsCL1; Multiplicity V0; Multiplicity SPD outer",2500,0,30000,700,0,7000);
fHOutMultV0OvsTRK = new TH2F("fHOutMultV0OvsTRK","fHOutMultV0OvsTRK; Multiplicity V0; Multiplicity TPC",2500,0,30000,400,0,4000);
fHOutMultCL1vsTKL = new TH2F ("fHOutMultCL1vsTKL","fHOutMultCL1vsTKL; Multiplicity SPD outer; Multiplicity tracklets",700,0,7000,700,0,7000);
+ fHOutMultZNAvsZPA = new TH2F ("fHOutMultZNAvsZPA","fHOutMultZNAvsZPA; Energy ZNA; Energy ZPA",500,0,2000,500,0,2000);
fHOutCentV0Mqual1 = new TH1F("fHOutCentV0M_qual1","fHOutCentV0M_qual1; Centrality V0",505,0,101);
fHOutCentTRKqual1 = new TH1F("fHOutCentTRK_qual1","fHOutCentTRK_qual1; Centrality TPC",505,0,101);
fOutputList->Add(fHOutCentV0M);
fOutputList->Add(fHOutCentV0A);
+ fOutputList->Add(fHOutCentV0A123);
fOutputList->Add(fHOutCentV0C);
fOutputList->Add(fHOutCentV0MEq);
fOutputList->Add(fHOutCentV0AEq);
fOutputList->Add(fHOutCentNPA);
fOutputList->Add(fHOutCentZNA);
fOutputList->Add(fHOutCentZNC);
+ fOutputList->Add(fHOutCentZPA);
+ fOutputList->Add(fHOutCentZPC);
fOutputList->Add(fHOutCentV0MvsFMD);
fOutputList->Add(fHOutCentTKLvsV0M);
fOutputList->Add(fHOutCentZEMvsZDC);
fOutputList->Add(fHOutCentZNAvsCentTRK);
fOutputList->Add(fHOutCentZNAvsCentCND);
fOutputList->Add(fHOutCentZNAvsCentCL1);
+ fOutputList->Add(fHOutCentZNAvsCentZPA);
fOutputList->Add(fHOutMultV0AC);
fOutputList->Add(fHOutMultV0M);
fOutputList->Add(fHOutMultV0A);
+ fOutputList->Add(fHOutMultV0A123);
fOutputList->Add(fHOutMultV0C);
fOutputList->Add(fHOutMultV0MEq);
fOutputList->Add(fHOutMultV0AEq);
fOutputList->Add(fHOutMultNPA);
fOutputList->Add(fHOutMultZNA);
fOutputList->Add(fHOutMultZNC);
+ fOutputList->Add(fHOutMultZPA);
+ fOutputList->Add(fHOutMultZPC);
fOutputList->Add(fHOutMultV0MvsZDN);
fOutputList->Add(fHOutMultZEMvsZDN);
fOutputList->Add(fHOutMultV0MvsZDC);
fOutputList->Add(fHOutMultV0OvsCL1);
fOutputList->Add(fHOutMultV0OvsTRK);
fOutputList->Add(fHOutMultCL1vsTKL);
+ fOutputList->Add(fHOutMultZNAvsZPA);
fOutputList->Add(fHOutCentV0Mqual1);
fOutputList->Add(fHOutCentTRKqual1);
fOutputList->Add(fHOutCentCL1qual1);
Float_t zem2Energy = 0.; // ZEM2 Energy
Bool_t zdcEnergyCal = kFALSE; // if zdc is calibrated (in pass2)
Double_t znaTower = 0.; // common PMT of ZNA
- Double_t zncTower = 0.;
+ Double_t zncTower = 0.; // common PMT of ZNC
+ Double_t zpaTower = 0.; // common PMT of ZPA
+ Double_t zpcTower = 0.; // common PMT of ZPC
Bool_t znaFired = kFALSE;
Bool_t zncFired = kFALSE;
+ Bool_t zpaFired = kFALSE;
+ Bool_t zpcFired = kFALSE;
Int_t nTracks = 0; // no. tracks
Int_t nTracklets = 0; // no. tracklets
Float_t spdCorr =0; // corrected spd2 multiplicity
Int_t multCND = 0; // no. tracks (candle condition)
+ Float_t multV0A123 = 0; // multiplicity from V0 reco side A (ring 1-2-3)
Float_t multV0A = 0; // multiplicity from V0 reco side A
Float_t multV0C = 0; // multiplicity from V0 reco side C
Float_t multV0AEq = 0; // multiplicity from V0 reco side A
multV0A=esdV0->GetMTotV0A();
multV0C=esdV0->GetMTotV0C();
+ // excluding innermost V0A ring
+ for (int i=40; i<64; i++)
+ multV0A123 += esdV0->GetMultiplicity(i);
+
multV0ACorr = AliESDUtils::GetCorrV0A(multV0A,zvtx);
multV0CCorr = AliESDUtils::GetCorrV0C(multV0C,zvtx);
Double_t mult = esd->GetVZEROEqMultiplicity(iCh);
multV0CEq += mult;
}
-
+
Bool_t kT0BB = kFALSE;
if (esd) {
// ***** T0 info
zem1Energy = (Float_t) (esdZDC->GetZDCEMEnergy(0))/8.;
zem2Energy = (Float_t) (esdZDC->GetZDCEMEnergy(1))/8.;
- const Double_t *ZNAtower = esdZDC->GetZN2TowerEnergy();
- znaTower = ZNAtower[0];
- const Double_t *ZNCtower = esdZDC->GetZN1TowerEnergy();
- zncTower = ZNCtower[0];
-
for (Int_t j = 0; j < 4; ++j)
if (esdZDC->GetZDCTDCData(12,j) != 0)
znaFired = kTRUE;
if (esdZDC->GetZDCTDCData(10,j) != 0)
zncFired = kTRUE;
+ for (Int_t j = 0; j < 4; ++j)
+ if (esdZDC->GetZDCTDCData(13,j) != 0)
+ zpaFired = kTRUE;
+
+ for (Int_t j = 0; j < 4; ++j)
+ if (esdZDC->GetZDCTDCData(11,j) != 0)
+ zpcFired = kTRUE;
+
+ const Double_t *ZNAtower = esdZDC->GetZN2TowerEnergy();
+ const Double_t *ZNCtower = esdZDC->GetZN1TowerEnergy();
+ const Double_t *ZPAtower = esdZDC->GetZP2TowerEnergy();
+ const Double_t *ZPCtower = esdZDC->GetZP1TowerEnergy();
+ if (znaFired) znaTower = ZNAtower[0];
+ if (zncFired) zncTower = ZNCtower[0];
+ if (zpaFired) zpaTower = ZPAtower[0];
+ if (zpcFired) zpcTower = ZPCtower[0];
+
} else {
AliAODHeader *h = aod->GetHeader();
- zncEnergy = (Float_t) (h->GetZDCN1Energy());
- zpcEnergy = (Float_t) (h->GetZDCP1Energy());
- znaEnergy = (Float_t) (h->GetZDCN2Energy());
- zpaEnergy = (Float_t) (h->GetZDCP2Energy());
+ zncEnergy = (Float_t) (h->GetZDCN1Energy());
+ zpcEnergy = (Float_t) (h->GetZDCP1Energy());
+ znaEnergy = (Float_t) (h->GetZDCN2Energy());
+ zpaEnergy = (Float_t) (h->GetZDCP2Energy());
zem1Energy = (Float_t) (h->GetZDCEMEnergy(0))/8.; //Todo: C.L. Should we devide here by 8? It is done in the ESD case!
zem2Energy = (Float_t) (h->GetZDCEMEnergy(1))/8.;
AliAODZDC *aodZDC = aod->GetZDCData();
const Double_t *ZNAtower = aodZDC->GetZNATowerEnergy();
+ const Double_t *ZNCtower = aodZDC->GetZNCTowerEnergy();
+ const Double_t *ZPAtower = aodZDC->GetZPATowerEnergy();
+ const Double_t *ZPCtower = aodZDC->GetZPCTowerEnergy();
znaTower = ZNAtower[0];
+ zncTower = ZNCtower[0];
+ zpaTower = ZPAtower[0];
+ zpcTower = ZPCtower[0];
znaFired = kFALSE; // trick because info is not stored in AOD
- if (esdCent->GetCentralityPercentile("ZNA") != 101)
- znaFired = kTRUE;
+ zncFired = kFALSE; // trick because info is not stored in AOD
+ zpaFired = kFALSE; // trick because info is not stored in AOD
+ zpcFired = kFALSE; // trick because info is not stored in AOD
+ if (esdCent->GetCentralityPercentile("ZNA") != 101) znaFired = kTRUE;
+ if (esdCent->GetCentralityPercentile("ZNC") != 101) zncFired = kTRUE;
+ if (esdCent->GetCentralityPercentile("ZPA") != 101) zpaFired = kTRUE;
+ if (esdCent->GetCentralityPercentile("ZPC") != 101) zpcFired = kTRUE;
}
if (esd) {
hHijing = (AliGenHijingEventHeader*)mcGenH;
else if (mcGenH->InheritsFrom(AliGenCocktailEventHeader::Class())) {
TList* headers = ((AliGenCocktailEventHeader*)mcGenH)->GetHeaders();
+ //headers->Print();
hHijing = dynamic_cast<AliGenHijingEventHeader*>(headers->FindObject("Hijing"));
if (!hHijing) hHijing = dynamic_cast<AliGenHijingEventHeader*>(headers->FindObject("Hijing pPb_0"));
+ if (!hHijing) hHijing = dynamic_cast<AliGenHijingEventHeader*>(headers->FindObject("Hijing_0"));
}
else if (mcGenH->InheritsFrom(AliGenDPMjetEventHeader::Class())) {
dpmHeader = (AliGenDPMjetEventHeader*)mcGenH;
// ***** Centrality Selection
if(fHtempV0M) fCentV0M = fHtempV0M->GetBinContent(fHtempV0M->FindBin((v0Corr)));
if(fHtempV0A) fCentV0A = fHtempV0A->GetBinContent(fHtempV0A->FindBin((multV0ACorr)));
+ if(fHtempV0A123) fCentV0A123 = fHtempV0A123->GetBinContent(fHtempV0A123->FindBin((multV0A123)));
if(fHtempV0C) fCentV0C = fHtempV0C->GetBinContent(fHtempV0C->FindBin((multV0CCorr)));
if(fHtempV0MEq) fCentV0MEq = fHtempV0MEq->GetBinContent(fHtempV0MEq->FindBin((multV0AEq+multV0CEq)));
if(fHtempV0AEq) fCentV0AEq = fHtempV0AEq->GetBinContent(fHtempV0AEq->FindBin((multV0AEq)));
if(zncFired) fCentZNC = fHtempZNC->GetBinContent(fHtempZNC->FindBin(zncTower));
else fCentZNC = 101;
}
+ if(fHtempZPA) {
+ if(znaFired) fCentZPA = fHtempZPA->GetBinContent(fHtempZPA->FindBin(zpaTower));
+ else fCentZPA = 101;
+ }
+ if(fHtempZPC) {
+ if(zpcFired) fCentZPC = fHtempZPC->GetBinContent(fHtempZPC->FindBin(zpcTower));
+ else fCentZPC = 101;
+ }
+
+
if(fHtempV0MvsFMD) fCentV0MvsFMD = fHtempV0MvsFMD->GetBinContent(fHtempV0MvsFMD->FindBin((multV0A+multV0C)));
if(fHtempTKLvsV0M) fCentTKLvsV0M = fHtempTKLvsV0M->GetBinContent(fHtempTKLvsV0M->FindBin(nTracklets));
if(fHtempZEMvsZDC) fCentZEMvsZDC = fHtempZEMvsZDC->GetBinContent(fHtempZEMvsZDC->FindBin(zem1Energy+zem2Energy,zncEnergy+znaEnergy+zpcEnergy+zpaEnergy));
esdCent->SetQuality(fQuality);
esdCent->SetCentralityV0M(fCentV0M);
esdCent->SetCentralityV0A(fCentV0A);
+ esdCent->SetCentralityV0A123(fCentV0A123);
esdCent->SetCentralityV0C(fCentV0C);
esdCent->SetCentralityV0MEq(fCentV0MEq);
esdCent->SetCentralityV0AEq(fCentV0AEq);
esdCent->SetCentralityNPA(fCentNPA);
esdCent->SetCentralityZNA(fCentZNA);
esdCent->SetCentralityZNC(fCentZNC);
+ esdCent->SetCentralityZPA(fCentZPA);
+ esdCent->SetCentralityZPC(fCentZPC);
esdCent->SetCentralityV0MvsFMD(fCentV0MvsFMD);
esdCent->SetCentralityTKLvsV0M(fCentTKLvsV0M);
esdCent->SetCentralityZEMvsZDC(fCentZEMvsZDC);
+
+ esdCent->SetCentralityV0Mtrue(fCentV0Mtrue);
+ esdCent->SetCentralityV0Atrue(fCentV0Atrue);
+ esdCent->SetCentralityV0Ctrue(fCentV0Ctrue);
+ // no setters in AliCentrality for
+ // fCentV0MEqtrue
+ // fCentV0AEqtrue
+ // fCentV0CEqtrue
+ esdCent->SetCentralityFMDtrue(fCentFMDtrue);
+ esdCent->SetCentralityTRKtrue(fCentTRKtrue);
+ esdCent->SetCentralityTKLtrue(fCentTKLtrue);
+ esdCent->SetCentralityCL0true(fCentCL0true);
+ esdCent->SetCentralityCL1true(fCentCL1true);
+ esdCent->SetCentralityCNDtrue(fCentCNDtrue);
+ esdCent->SetCentralityZNAtrue(fCentZNAtrue);
+ esdCent->SetCentralityZNCtrue(fCentZNCtrue);
+ esdCent->SetCentralityZPAtrue(fCentZPAtrue);
+ esdCent->SetCentralityZPCtrue(fCentZPCtrue);
}
// filling QA histograms
fCentCND = fCentCNDtrue;
fCentZNA = fCentZNAtrue;
fCentZNC = fCentZNCtrue;
+ fCentZPA = fCentZPAtrue;
+ fCentZPC = fCentZPCtrue;
}
if (fQuality==0) {
fHOutCentV0M->Fill(fCentV0M);
fHOutCentV0A->Fill(fCentV0A);
+ fHOutCentV0A123->Fill(fCentV0A123);
fHOutCentV0C->Fill(fCentV0C);
fHOutCentV0MEq->Fill(fCentV0MEq);
fHOutCentV0AEq->Fill(fCentV0AEq);
fHOutCentNPA->Fill(fCentNPA);
fHOutCentZNA->Fill(fCentZNA);
fHOutCentZNC->Fill(fCentZNC);
+ fHOutCentZPA->Fill(fCentZPA);
+ fHOutCentZPC->Fill(fCentZPC);
fHOutCentV0MvsFMD->Fill(fCentV0MvsFMD);
fHOutCentTKLvsV0M->Fill(fCentTKLvsV0M);
fHOutCentZEMvsZDC->Fill(fCentZEMvsZDC);
fHOutCentZNAvsCentTRK->Fill(fCentZNA,fCentTRK);
fHOutCentZNAvsCentCND->Fill(fCentZNA,fCentCND);
fHOutCentZNAvsCentCL1->Fill(fCentZNA,fCentCL1);
+ fHOutCentZNAvsCentZPA->Fill(fCentZNA,fCentZPA);
fHOutMultV0AC->Fill(multV0A,multV0C);
fHOutMultV0M->Fill(multV0ACorr+multV0CCorr);
fHOutMultV0A->Fill(multV0ACorr);
+ fHOutMultV0A123->Fill(multV0A123);
fHOutMultV0C->Fill(multV0CCorr);
fHOutMultV0MEq->Fill(multV0AEq+multV0CEq);
fHOutMultV0AEq->Fill(multV0AEq);
fHOutMultNPA->Fill(Npart);
if(znaFired)fHOutMultZNA->Fill(znaTower);
if(zncFired)fHOutMultZNC->Fill(zncTower);
+ if(zpaFired)fHOutMultZPA->Fill(zpaTower);
+ if(zpcFired)fHOutMultZPC->Fill(zpcTower);
+ fHOutMultZNAvsZPA->Fill(znaTower,zpaTower);
- fHOutMultV0MvsZDN->Fill(v0Corr,(zncEnergy+znaEnergy));
+ //fHOutMultV0MvsZDN->Fill(v0Corr,(zncEnergy+znaEnergy));
+ fHOutMultV0MvsZDN->Fill(v0Corr,znaTower);
fHOutMultZEMvsZDN->Fill((zem1Energy+zem2Energy),(zncEnergy+znaEnergy));
fHOutMultV0MvsZDC->Fill(v0Corr,(zncEnergy+znaEnergy+zpcEnergy+zpaEnergy));
fHOutMultZEMvsZDC->Fill((zem1Energy+zem2Energy),(zncEnergy+znaEnergy+zpcEnergy+zpaEnergy));
// centrality histos
fHtempV0M = centOADB->V0hist();
fHtempV0A = centOADB->V0Ahist();
+ fHtempV0A123 = centOADB->V0A123hist();
fHtempV0C = centOADB->V0Chist();
fHtempV0MEq = centOADB->V0Eqhist();
fHtempV0AEq = centOADB->V0AEqhist();
fHtempFMD = centOADB->FMDhist();
fHtempZNA = centOADB->ZNAhist();
fHtempZNC = centOADB->ZNChist();
+ fHtempZPA = centOADB->ZPAhist();
+ fHtempZPC = centOADB->ZPChist();
fHtempZEMvsZDC = centOADB->ZEMvsZDChist();
if (isHijing) {
fHtempFMDtrue = centOADB->FMDhisttrue();
fHtempZNAtrue = centOADB->ZNAhisttrue();
fHtempZNCtrue = centOADB->ZNChisttrue();
+ fHtempZPAtrue = centOADB->ZPAhisttrue();
+ fHtempZPCtrue = centOADB->ZPChisttrue();
} else if (isDpmjet) {
fHtempNPA = centOADB->NPAhistDPM();
fHtempV0Mtrue = centOADB->V0histtrueDPM();
fHtempFMDtrue = centOADB->FMDhisttrueDPM();
fHtempZNAtrue = centOADB->ZNAhisttrueDPM();
fHtempZNCtrue = centOADB->ZNChisttrueDPM();
+ fHtempZPAtrue = centOADB->ZPAhisttrueDPM();
+ fHtempZPCtrue = centOADB->ZPChisttrueDPM();
}
if (!fHtempCND) AliWarning(Form("Calibration for CND does not exist in %s", path.Data()));
if (!fHtempZNA) AliWarning(Form("Calibration for ZNA does not exist in %s", path.Data()));
if (!fHtempZNC) AliWarning(Form("Calibration for ZNC does not exist in %s", path.Data()));
+ if (!fHtempZPA) AliWarning(Form("Calibration for ZPA does not exist in %s", path.Data()));
+ if (!fHtempZPC) AliWarning(Form("Calibration for ZPC does not exist in %s", path.Data()));
if (!fHtempFMD) AliWarning(Form("Calibration for FMD does not exist in %s", path.Data()));
if (!fHtempZEMvsZDC) AliWarning(Form("Calibration for ZEMvsZDC does not exist in %s", path.Data()));
if (!fHtempNPA) AliWarning(Form("Calibration for NPA does not exist in %s", path.Data()));
if (!fHtempCNDtrue) AliWarning(Form("Calibration for CNDtrue does not exist in %s", path.Data()));
if (!fHtempZNAtrue) AliWarning(Form("Calibration for ZNAtrue does not exist in %s", path.Data()));
if (!fHtempZNCtrue) AliWarning(Form("Calibration for ZNCtrue does not exist in %s", path.Data()));
+ if (!fHtempZPAtrue) AliWarning(Form("Calibration for ZPAtrue does not exist in %s", path.Data()));
+ if (!fHtempZPCtrue) AliWarning(Form("Calibration for ZPCtrue does not exist in %s", path.Data()));
if (!fHtempFMDtrue) AliWarning(Form("Calibration for FMDtrue does not exist in %s", path.Data()));