Adding histos to model fragment prod. from data

diff --git a/PWGLF/FORWARD/GEO/AliAnalysisTaskZDCPbPb.cxx b/PWGLF/FORWARD/GEO/AliAnalysisTaskZDCPbPb.cxx
index 1b05acdd21d285e7a6c54feb5b71a43c5ea9e6dd..2f412df56faf56c1ecb07b2d7fa3eaa9f8c6fe35 100644 (file)
--- a/PWGLF/FORWARD/GEO/AliAnalysisTaskZDCPbPb.cxx
+++ b/PWGLF/FORWARD/GEO/AliAnalysisTaskZDCPbPb.cxx
@@ -91,7 +91,11 @@ AliAnalysisTaskZDCPbPb::AliAnalysisTaskZDCPbPb():
     fhZNCvscentrality(0x0),
     fhZNAvscentrality(0x0),
     fhZPCvscentrality(0x0),
-    fhZPAvscentrality(0x0)
+    fhZPAvscentrality(0x0),
+    fhZNCpmcvscentr(0x0),
+    fhZNApmcvscentr(0x0),
+    fhZPCpmcvscentr(0x0),
+    fhZPApmcvscentr(0x0)
     
 {   
    // Default constructor
@@ -152,7 +156,11 @@ AliAnalysisTaskZDCPbPb::AliAnalysisTaskZDCPbPb(const char *name):
     fhZNCvscentrality(0x0),
     fhZNAvscentrality(0x0),
     fhZPCvscentrality(0x0),
-    fhZPAvscentrality(0x0)
+    fhZPAvscentrality(0x0),
+    fhZNCpmcvscentr(0x0),
+    fhZNApmcvscentr(0x0),
+    fhZPCpmcvscentr(0x0),
+    fhZPApmcvscentr(0x0)
 
 {
 
@@ -229,7 +237,11 @@ AliAnalysisTaskZDCPbPb::AliAnalysisTaskZDCPbPb(const AliAnalysisTaskZDCPbPb& ana
     fhZNCvscentrality(ana.fhZNCvscentrality),
     fhZNAvscentrality(ana.fhZNAvscentrality),
     fhZPCvscentrality(ana.fhZPCvscentrality),
-    fhZPAvscentrality(ana.fhZPAvscentrality)
+    fhZPAvscentrality(ana.fhZPAvscentrality),
+    fhZNCpmcvscentr(ana.fhZNCpmcvscentr),
+    fhZNApmcvscentr(ana.fhZNApmcvscentr),
+    fhZPCpmcvscentr(ana.fhZPCpmcvscentr),
+    fhZPApmcvscentr(ana.fhZPApmcvscentr)
 
 {
   //
@@ -435,6 +447,15 @@ void AliAnalysisTaskZDCPbPb::UserCreateOutputObjects()
   fhZPAvscentrality = new TH2F("hZPAvscentrality","ZPA vs. centrality",100,0.,100.,100,0.,100.);
   fOutput->Add(fhZPAvscentrality);
   
+  fhZNCpmcvscentr = new TH2F("hZNCpmcvscentr","ZNC PMC vs. centrality",100,0.,100.,100,0.,100.);
+  fOutput->Add(fhZNCpmcvscentr);
+  fhZNApmcvscentr = new TH2F("hZNApmcvscentr","ZNA PMC vs. centrality",100,0.,100.,100,0.,100.);
+  fOutput->Add(fhZNApmcvscentr);
+  fhZPCpmcvscentr = new TH2F("hZPCpmcvscentr","ZPC PMC vs. centrality",100,0.,100.,100,0.,100.);
+  fOutput->Add(fhZPCpmcvscentr);
+  fhZPApmcvscentr = new TH2F("hZPApmcvscentr","ZPA PMC vs. centrality",100,0.,100.,100,0.,100.);
+  fOutput->Add(fhZPApmcvscentr);
+  
   PostData(1, fOutput);
 }
 
@@ -618,6 +639,11 @@ void AliAnalysisTaskZDCPbPb::UserExec(Option_t */*option*/)
        fhZPCvscentrality->Fill(centrperc, energyZPC/1000.);
        fhZPAvscentrality->Fill(centrperc, energyZPA/1000.);
        
+       fhZNCpmcvscentr->Fill(centrperc, towZNC[0]/1000.);
+       fhZNApmcvscentr->Fill(centrperc, towZNA[0]/1000.);
+       fhZPCpmcvscentr->Fill(centrperc, towZPC[0]/1000.);
+       fhZPApmcvscentr->Fill(centrperc, towZPA[0]/1000.);
+       
       }
 //    } // PHYSICS SELECTION
       
@@ -738,6 +764,11 @@ void AliAnalysisTaskZDCPbPb::UserExec(Option_t */*option*/)
       fhZNAvscentrality->Fill(centrperc, energyZNA/1000.);
       fhZPCvscentrality->Fill(centrperc, energyZPC/1000.);
       fhZPAvscentrality->Fill(centrperc, energyZPA/1000.);
+       
+      fhZNCpmcvscentr->Fill(centrperc, towZNC[0]/1000.);
+      fhZNApmcvscentr->Fill(centrperc, towZNA[0]/1000.);
+      fhZPCpmcvscentr->Fill(centrperc, towZPC[0]/1000.);
+      fhZPApmcvscentr->Fill(centrperc, towZPA[0]/1000.);
     //}
   }
   

diff --git a/PWGLF/FORWARD/GEO/AliAnalysisTaskZDCPbPb.h b/PWGLF/FORWARD/GEO/AliAnalysisTaskZDCPbPb.h
index cfe52f8f95f963101f502582c1a58836e8f9d42a..a37abc1d16fe134b0e642bb7f53310d53f782655 100644 (file)
--- a/PWGLF/FORWARD/GEO/AliAnalysisTaskZDCPbPb.h
+++ b/PWGLF/FORWARD/GEO/AliAnalysisTaskZDCPbPb.h
@@ -99,6 +99,11 @@ class AliAnalysisTaskZDCPbPb : public AliAnalysisTaskSE {
   TH2F *fhZNAvscentrality;     //! ZNA vs. centrality
   TH2F *fhZPCvscentrality;     //! ZPC vs. centrality
   TH2F *fhZPAvscentrality;     //! ZPA vs. centrality
+  //
+  TH2F *fhZNCpmcvscentr;       //! ZNC vs. centrality
+  TH2F *fhZNApmcvscentr;   //! ZNA vs. centrality
+  TH2F *fhZPCpmcvscentr;   //! ZPC vs. centrality
+  TH2F *fhZPApmcvscentr;   //! ZPA vs. centrality
  
   ClassDef(AliAnalysisTaskZDCPbPb,2);