New control histo.

Alberica

diff --git a/ANALYSIS/AliCentralitySelectionTask.cxx b/ANALYSIS/AliCentralitySelectionTask.cxx
index 7a2708b1dc8e7adf12f02d3002cfbdce2bd00e40..9c5b8c8385270f2fe523acc25aafe40427e23cc4 100644 (file)
--- a/ANALYSIS/AliCentralitySelectionTask.cxx
+++ b/ANALYSIS/AliCentralitySelectionTask.cxx
@@ -265,6 +265,7 @@ AliAnalysisTaskSE(),
   fHOutCentZNAvsCentTRK(0),
   fHOutCentZNAvsCentCND(0),
   fHOutCentZNAvsCentCL1(0),
+  fHOutCentZNAvsCentZPA(0),
   fHOutMultV0AC(0),
   fHOutMultV0M(0),
   fHOutMultV0A(0),
@@ -300,6 +301,7 @@ AliAnalysisTaskSE(),
   fHOutMultV0OvsCL1(0),
   fHOutMultV0OvsTRK(0),
   fHOutMultCL1vsTKL(0),
+  fHOutMultZNAvsZPA(0),
   fHOutCentV0Mqual1(0),
   fHOutCentTRKqual1(0),
   fHOutCentCL1qual1(0),
@@ -513,6 +515,7 @@ AliCentralitySelectionTask::AliCentralitySelectionTask(const char *name):
   fHOutCentZNAvsCentTRK(0),
   fHOutCentZNAvsCentCND(0),
   fHOutCentZNAvsCentCL1(0),
+  fHOutCentZNAvsCentZPA(0),
   fHOutMultV0AC(0),
   fHOutMultV0M(0),
   fHOutMultV0A(0),
@@ -548,6 +551,7 @@ AliCentralitySelectionTask::AliCentralitySelectionTask(const char *name):
   fHOutMultV0OvsCL1(0),
   fHOutMultV0OvsTRK(0),
   fHOutMultCL1vsTKL(0),
+  fHOutMultZNAvsZPA(0),
   fHOutCentV0Mqual1(0),
   fHOutCentTRKqual1(0),
   fHOutCentCL1qual1(0),
@@ -771,6 +775,7 @@ AliCentralitySelectionTask::AliCentralitySelectionTask(const AliCentralitySelect
   fHOutCentZNAvsCentTRK(ana.fHOutCentZNAvsCentTRK),
   fHOutCentZNAvsCentCND(ana.fHOutCentZNAvsCentCND),
   fHOutCentZNAvsCentCL1(ana.fHOutCentZNAvsCentCL1),
+  fHOutCentZNAvsCentZPA(ana.fHOutCentZNAvsCentZPA),
   fHOutMultV0AC(ana.fHOutMultV0AC),
   fHOutMultV0M(ana.fHOutMultV0M),
   fHOutMultV0A(ana.fHOutMultV0A),
@@ -806,6 +811,7 @@ AliCentralitySelectionTask::AliCentralitySelectionTask(const AliCentralitySelect
   fHOutMultV0OvsCL1(ana.fHOutMultV0OvsCL1),
   fHOutMultV0OvsTRK(ana.fHOutMultV0OvsTRK),
   fHOutMultCL1vsTKL(ana.fHOutMultCL1vsTKL),
+  fHOutMultZNAvsZPA(ana.fHOutMultZNAvsZPA),
   fHOutCentV0Mqual1(ana.fHOutCentV0Mqual1),
   fHOutCentTRKqual1(ana.fHOutCentTRKqual1),
   fHOutCentCL1qual1(ana.fHOutCentCL1qual1),
@@ -910,6 +916,7 @@ void AliCentralitySelectionTask::UserCreateOutputObjects()
     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);
@@ -950,6 +957,7 @@ void AliCentralitySelectionTask::UserCreateOutputObjects()
     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);
@@ -1032,6 +1040,7 @@ void AliCentralitySelectionTask::UserCreateOutputObjects()
     fOutputList->Add(fHOutCentZNAvsCentTRK);
     fOutputList->Add(fHOutCentZNAvsCentCND);
     fOutputList->Add(fHOutCentZNAvsCentCL1);
+    fOutputList->Add(fHOutCentZNAvsCentZPA);
 
     fOutputList->Add(fHOutMultV0AC); 
     fOutputList->Add(fHOutMultV0M); 
@@ -1068,6 +1077,7 @@ void AliCentralitySelectionTask::UserCreateOutputObjects()
     fOutputList->Add(fHOutMultV0OvsCL1);
     fOutputList->Add(fHOutMultV0OvsTRK);
     fOutputList->Add(fHOutMultCL1vsTKL);
+    fOutputList->Add(fHOutMultZNAvsZPA);
     fOutputList->Add(fHOutCentV0Mqual1);
     fOutputList->Add(fHOutCentTRKqual1);
     fOutputList->Add(fHOutCentCL1qual1);                   
@@ -1422,15 +1432,6 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/)
     zem1Energy = (Float_t) (esdZDC->GetZDCEMEnergy(0))/8.;
     zem2Energy = (Float_t) (esdZDC->GetZDCEMEnergy(1))/8.;
 
-    const Double_t *ZNAtower = esdZDC->GetZN2TowerEnergy(); 
-    const Double_t *ZNCtower = esdZDC->GetZN1TowerEnergy();
-    const Double_t *ZPAtower = esdZDC->GetZP2TowerEnergy(); 
-    const Double_t *ZPCtower = esdZDC->GetZP1TowerEnergy();
-    znaTower = ZNAtower[0];
-    zncTower = ZNCtower[0];
-    zpaTower = ZPAtower[0];
-    zpcTower = ZPCtower[0];
-
     for (Int_t j = 0; j < 4; ++j) 
       if (esdZDC->GetZDCTDCData(12,j) != 0) 
        znaFired = kTRUE;
@@ -1447,6 +1448,15 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/)
       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());
@@ -1727,6 +1737,7 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/)
        fHOutCentZNAvsCentTRK->Fill(fCentZNA,fCentTRK);
        fHOutCentZNAvsCentCND->Fill(fCentZNA,fCentCND);
        fHOutCentZNAvsCentCL1->Fill(fCentZNA,fCentCL1);
+       fHOutCentZNAvsCentZPA->Fill(fCentZNA,fCentZPA);
 
        fHOutMultV0AC->Fill(multV0A,multV0C);
        fHOutMultV0M->Fill(multV0ACorr+multV0CCorr);
@@ -1751,7 +1762,7 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/)
        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,znaTower);

diff --git a/ANALYSIS/AliCentralitySelectionTask.h b/ANALYSIS/AliCentralitySelectionTask.h
index 43c1112f4a98a87c9a41e4ace8708281220c2d7f..0d56668225a857790499074c2ed0ee310606d218 100644 (file)
--- a/ANALYSIS/AliCentralitySelectionTask.h
+++ b/ANALYSIS/AliCentralitySelectionTask.h
@@ -248,6 +248,7 @@ class AliCentralitySelectionTask : public AliAnalysisTaskSE {
   TH2F *fHOutCentZNAvsCentTRK;  //control histogram for centrality
   TH2F *fHOutCentZNAvsCentCND;  //control histogram for centrality
   TH2F *fHOutCentZNAvsCentCL1;  //control histogram for centrality
+  TH2F *fHOutCentZNAvsCentZPA;  //control histogram for centrality
 
   TH2F *fHOutMultV0AC;          //control histogram for multiplicity
   TH1F *fHOutMultV0M ;          //control histogram for multiplicity
@@ -285,6 +286,7 @@ class AliCentralitySelectionTask : public AliAnalysisTaskSE {
   TH2F *fHOutMultV0OvsCL1;      //control histogram for multiplicity
   TH2F *fHOutMultV0OvsTRK;      //control histogram for multiplicity
   TH2F *fHOutMultCL1vsTKL;      //control histogram for multiplicity
+  TH2F *fHOutMultZNAvsZPA;      //control histogram for multiplicity
 
   TH1F *fHOutCentV0Mqual1     ; //control histogram for centrality quality 1
   TH1F *fHOutCentTRKqual1     ; //control histogram for centrality quality 1
@@ -304,7 +306,7 @@ class AliCentralitySelectionTask : public AliAnalysisTaskSE {
   TH1F *fHOutVertex ;           //control histogram for vertex SPD
   TH1F *fHOutVertexT0 ;         //control histogram for vertex T0
 
-  ClassDef(AliCentralitySelectionTask, 27); 
+  ClassDef(AliCentralitySelectionTask, 28); 
 };
 
 #endif