update

diff --git a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.cxx b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.cxx
index 3d638cf10e5edd8bc019006a14eed02488aa3685..ec43bffeade5b63c28b2710d75b5f3b51651d1e7 100644 (file)
--- a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.cxx
+++ b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.cxx
@@ -300,25 +300,25 @@ void AliAnalysisTaskSOH::UserCreateOutputObjects()
     fOutputList->Add(fHClsEoverMcE_Pion);
   }
   
-  fHParGenPion_p = new TH3F("fHParGenPion_p","Particle level truth Phi-Eta-p_{T} distribution of #pi+", 500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+  fHParGenPion_p = new TH3F("fHParGenPion_p","Particle level truth Phi-Eta-p_{T} distribution of #pi+",  500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
   fHParGenPion_p->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
   fHParGenPion_p->GetYaxis()->SetTitle("#eta");
   fHParGenPion_p->GetZaxis()->SetTitle("#phi");
   fOutputList->Add(fHParGenPion_p);
   
-  fHParGenPion_m = new TH3F("fHParGenPion_m", "Particle level truth Phi-Eta-p_{T} distribution of all #pi-", 500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+  fHParGenPion_m = new TH3F("fHParGenPion_m", "Particle level truth Phi-Eta-p_{T} distribution of all #pi-", 500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
   fHParGenPion_m->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
   fHParGenPion_m->GetYaxis()->SetTitle("#eta");
   fHParGenPion_m->GetZaxis()->SetTitle("#phi");
   fOutputList->Add(fHParGenPion_m);
   
-  fHParGenPion_rmInj_p = new TH3F("fHParGenPion_rmInj_p","Particle level truth Phi-Eta-p_{T} distribution of all #pi+ without injected signal", 500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+  fHParGenPion_rmInj_p = new TH3F("fHParGenPion_rmInj_p","Particle level truth Phi-Eta-p_{T} distribution of all #pi+ without injected signal",  500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
   fHParGenPion_rmInj_p->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
   fHParGenPion_rmInj_p->GetYaxis()->SetTitle("#eta");
   fHParGenPion_rmInj_p->GetZaxis()->SetTitle("#phi");
   fOutputList->Add(fHParGenPion_rmInj_p);
   
-  fHParGenPion_rmInj_m = new TH3F("fHParGenPion_rmInj_m","Particle level truth Phi-Eta-p_{T} distribution of #pi- without injected signal", 500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+  fHParGenPion_rmInj_m = new TH3F("fHParGenPion_rmInj_m","Particle level truth Phi-Eta-p_{T} distribution of #pi- without injected signal", 500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
   fHParGenPion_rmInj_m->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
   fHParGenPion_rmInj_m->GetYaxis()->SetTitle("#eta");
   fHParGenPion_rmInj_m->GetZaxis()->SetTitle("#phi");
@@ -330,26 +330,26 @@ void AliAnalysisTaskSOH::UserCreateOutputObjects()
   if(fMcProcess && fTrackProcess)
   {    
     //Fake
-    fHDetGenFakePion = new TH3F("fHDetGenFakePion", "fake charged pion track Phi-Eta-p_{T} distribution", 500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+    fHDetGenFakePion = new TH3F("fHDetGenFakePion", "fake charged pion track Phi-Eta-p_{T} distribution",500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
     fHDetGenFakePion->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
     fHDetGenFakePion->GetYaxis()->SetTitle("#eta");
     fHDetGenFakePion->GetZaxis()->SetTitle("#phi");
     fOutputList->Add(fHDetGenFakePion);
 
-    fHDetRecFakePion = new TH3F("fHDetRecFakePion", "fake charged pion track Phi-Eta-p_{T} distribution", 500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+    fHDetRecFakePion = new TH3F("fHDetRecFakePion", "fake charged pion track Phi-Eta-p_{T} distribution", 500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
     fHDetRecFakePion->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
     fHDetRecFakePion->GetYaxis()->SetTitle("#eta");
     fHDetRecFakePion->GetZaxis()->SetTitle("#phi");
     fOutputList->Add(fHDetRecFakePion);
     
     //Secondary
-    fHDetGenSecPion = new TH3F("fHDetGenSecPion", "secondary charged pion charged pion track Phi-Eta-p_{T} distribution", 500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+    fHDetGenSecPion = new TH3F("fHDetGenSecPion", "secondary charged pion charged pion track Phi-Eta-p_{T} distribution", 500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
     fHDetGenSecPion->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
     fHDetGenSecPion->GetYaxis()->SetTitle("#eta");
     fHDetGenSecPion->GetZaxis()->SetTitle("#phi");
     fOutputList->Add(fHDetGenSecPion);
 
-    fHDetRecSecPion = new TH3F("fHDetRecSecPion", "secondary charged pion charged pion track Phi-Eta-p_{T} distribution", 500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+    fHDetRecSecPion = new TH3F("fHDetRecSecPion", "secondary charged pion charged pion track Phi-Eta-p_{T} distribution", 500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
     fHDetRecSecPion->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
     fHDetRecSecPion->GetYaxis()->SetTitle("#eta");
     fHDetRecSecPion->GetZaxis()->SetTitle("#phi");
@@ -364,46 +364,46 @@ void AliAnalysisTaskSOH::UserCreateOutputObjects()
       fHDetGenPion_p[i]->GetZaxis()->SetTitle("#phi");
       fOutputList->Add(fHDetGenPion_p[i]);
       
-      fHDetRecPion_p[i] = new TH3F(Form("fHDetRecPion_p_%s", trackCut[i]), Form("%s: Reconstructed track Phi-Eta-p_{T} distribution of all #pi+", trackCut[i]),  500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+      fHDetRecPion_p[i] = new TH3F(Form("fHDetRecPion_p_%s", trackCut[i]), Form("%s: Reconstructed track Phi-Eta-p_{T} distribution of all #pi+", trackCut[i]),  500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
       fHDetRecPion_p[i]->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
       fHDetRecPion_p[i]->GetYaxis()->SetTitle("#eta");
       fHDetRecPion_p[i]->GetZaxis()->SetTitle("#phi");
       fOutputList->Add(fHDetRecPion_p[i]);
 
       // pi-
-      fHDetGenPion_m[i] = new TH3F(Form("fHDetGenPion_m_%s", trackCut[i]), Form("%s: Detector level truth Phi-Eta-p_{T} distribution of #pi-", trackCut[i]),  500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+      fHDetGenPion_m[i] = new TH3F(Form("fHDetGenPion_m_%s", trackCut[i]), Form("%s: Detector level truth Phi-Eta-p_{T} distribution of #pi-", trackCut[i]),  500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
       fHDetGenPion_m[i]->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
       fHDetGenPion_m[i]->GetYaxis()->SetTitle("#eta");
       fHDetGenPion_m[i]->GetZaxis()->SetTitle("#phi");
       fOutputList->Add(fHDetGenPion_m[i]);
       
-      fHDetRecPion_m[i] = new TH3F(Form("fHDetRecPion_m_%s", trackCut[i]), Form("%s: Reconstructed track Phi-Eta-p_{T} distribution of #pi-", trackCut[i]),  500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+      fHDetRecPion_m[i] = new TH3F(Form("fHDetRecPion_m_%s", trackCut[i]), Form("%s: Reconstructed track Phi-Eta-p_{T} distribution of #pi-", trackCut[i]),  500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
       fHDetRecPion_m[i]->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
       fHDetRecPion_m[i]->GetYaxis()->SetTitle("#eta");
       fHDetRecPion_m[i]->GetZaxis()->SetTitle("#phi");
       fOutputList->Add(fHDetRecPion_m[i]);
 
         //pi+ without injected signal
-      fHDetGenPion_rmInj_p[i] = new TH3F(Form("fHDetGenPion_rmInj_p_%s", trackCut[i]), Form("%s: Detector level truth Phi-Eta-p_{T} distribution of #pi+ without injected signal", trackCut[i]),  500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+      fHDetGenPion_rmInj_p[i] = new TH3F(Form("fHDetGenPion_rmInj_p_%s", trackCut[i]), Form("%s: Detector level truth Phi-Eta-p_{T} distribution of #pi+ without injected signal", trackCut[i]), 500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
       fHDetGenPion_rmInj_p[i]->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
       fHDetGenPion_rmInj_p[i]->GetYaxis()->SetTitle("#eta");
       fHDetGenPion_rmInj_p[i]->GetZaxis()->SetTitle("#phi");
       fOutputList->Add(fHDetGenPion_rmInj_p[i]);
       
-      fHDetRecPion_rmInj_p[i] = new TH3F(Form("fHDetRecPion_rmInj_p_%s", trackCut[i]), Form("%s: Reconstructed track Phi-Eta-p_{T} distribution of #pi+ without injected signal", trackCut[i]),  500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+      fHDetRecPion_rmInj_p[i] = new TH3F(Form("fHDetRecPion_rmInj_p_%s", trackCut[i]), Form("%s: Reconstructed track Phi-Eta-p_{T} distribution of #pi+ without injected signal", trackCut[i]), 500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
       fHDetRecPion_rmInj_p[i]->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
       fHDetRecPion_rmInj_p[i]->GetYaxis()->SetTitle("#eta");
       fHDetRecPion_rmInj_p[i]->GetZaxis()->SetTitle("#phi");
       fOutputList->Add(fHDetRecPion_rmInj_p[i]);
 
       //pi- charged particle without injected signal
-      fHDetGenPion_rmInj_m[i] = new TH3F(Form("fHDetGenPion_rmInj_m_%s", trackCut[i]), Form("%s: Detector level truth Phi-Eta-p_{T} distribution of #pi- without injected signal", trackCut[i]),  500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+      fHDetGenPion_rmInj_m[i] = new TH3F(Form("fHDetGenPion_rmInj_m_%s", trackCut[i]), Form("%s: Detector level truth Phi-Eta-p_{T} distribution of #pi- without injected signal", trackCut[i]), 500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
       fHDetGenPion_rmInj_m[i]->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
       fHDetGenPion_rmInj_m[i]->GetYaxis()->SetTitle("#eta");
       fHDetGenPion_rmInj_m[i]->GetZaxis()->SetTitle("#phi");
       fOutputList->Add(fHDetGenPion_rmInj_m[i]);
       
-      fHDetRecPion_rmInj_m[i] = new TH3F(Form("fHDetRecPion_rmInj_m_%s", trackCut[i]), Form("%s: Reconstructed track Phi-Eta-p_{T} distribution of #pi- without injected signal", trackCut[i]),  500, 0, 100, 60, -1.2, 1.2, 128, 0, 6.4);
+      fHDetRecPion_rmInj_m[i] = new TH3F(Form("fHDetRecPion_rmInj_m_%s", trackCut[i]), Form("%s: Reconstructed track Phi-Eta-p_{T} distribution of #pi- without injected signal", trackCut[i]), 500, 0, 100, 100, -1.0, 1.0, 120, 0.0,2.*TMath::Pi());
       fHDetRecPion_rmInj_m[i]->GetXaxis()->SetTitle("p_{T}^{gen} (GeV/c)");
       fHDetRecPion_rmInj_m[i]->GetYaxis()->SetTitle("#eta");
       fHDetRecPion_rmInj_m[i]->GetZaxis()->SetTitle("#phi");

diff --git a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.h b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.h
index 3da8053ecef6f3bff48e7860e9244ddd9afe2907..5dddd52a59ba6652388e36e74d4d632d2e72f2ba 100644 (file)
--- a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.h
+++ b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSOH.h
@@ -112,6 +112,6 @@ class AliAnalysisTaskSOH : public AliAnalysisTaskSE {
   AliAnalysisTaskSOH(const AliAnalysisTaskSOH&); // not implemented
   AliAnalysisTaskSOH& operator=(const AliAnalysisTaskSOH&); // not implemented
   
-  ClassDef(AliAnalysisTaskSOH, 12); // Analysis task Saehanseul Oh
+  ClassDef(AliAnalysisTaskSOH, 13); // Analysis task Saehanseul Oh
 };
 #endif

diff --git a/PWGJE/EMCALJetTasks/macros/AddTaskSOH.C b/PWGJE/EMCALJetTasks/macros/AddTaskSOH.C
index 1af24dffdda7a8b937c80eb1c21bf9883a4b5166..8eb6accb7196b506ca6cc31c365499c8c4693fd7 100644 (file)
--- a/PWGJE/EMCALJetTasks/macros/AddTaskSOH.C
+++ b/PWGJE/EMCALJetTasks/macros/AddTaskSOH.C
@@ -1,6 +1,6 @@
 // $Id$
 
-AliAnalysisTaskSOH* AddTaskSOH(const char *name, Bool_t ScaleFactorP=kFALSE, Bool_t ClusterP=kFALSE, Double_t Zvtx=10)
+AliAnalysisTaskSOH* AddTaskSOH()
 {
   // Get the pointer to the existing analysis manager via the static access method.
   //==============================================================================
@@ -23,7 +23,7 @@ AliAnalysisTaskSOH* AddTaskSOH(const char *name, Bool_t ScaleFactorP=kFALSE, Boo
   // Init the task and do settings
   //-------------------------------------------------------
 
-  AliAnalysisTaskSOH *taskSOH = new AliAnalysisTaskSOH(Form("QA_soh_%s",name));
+  AliAnalysisTaskSOH *taskSOH = new AliAnalysisTaskSOH("AliAnalysisTaskSOH");
 
   AliESDtrackCuts *esdTrackCuts = 0x0;
   AliESDtrackCuts *hybridTrackCuts1 = 0x0;
@@ -40,9 +40,9 @@ AliAnalysisTaskSOH* AddTaskSOH(const char *name, Bool_t ScaleFactorP=kFALSE, Boo
 
   taskSOH->SetMcProcess(kTRUE);
   taskSOH->SetTrackProcess(kTRUE);
-  taskSOH->SetSFProcess(ScaleFactorP);
-  taskSOH->SetClusterProcess(ClusterP);
-  taskSOH->SetZvtx(Zvtx);
+  taskSOH->SetSFProcess(kFALSE);
+  taskSOH->SetClusterProcess(kFALSE);
+  taskSOH->SetZvtx(10);
 
   // Add task(s)
   mgr->AddTask(taskSOH); 
@@ -58,7 +58,10 @@ AliAnalysisTaskSOH* AddTaskSOH(const char *name, Bool_t ScaleFactorP=kFALSE, Boo
 
   // Create containers for input/output
   AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
-  AliAnalysisDataContainer *coutputpt = mgr->CreateContainer(name, TList::Class(), AliAnalysisManager::kOutputContainer, "AnalysisResults_soh.root");
+  AliAnalysisDataContainer *coutputpt = mgr->CreateContainer("soh", 
+                                                             TList::Class(), 
+                                                             AliAnalysisManager::kOutputContainer, 
+                                                             "AnalysisResults.root");
 
   // Connect input/output
   mgr->ConnectInput(taskSOH, 0, cinput);