]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
updated
authorssakai <ssakai@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 28 Jan 2013 13:07:14 +0000 (13:07 +0000)
committerssakai <ssakai@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 28 Jan 2013 13:07:14 +0000 (13:07 +0000)
PWGHF/hfe/AliAnalysisTaskHFECal.cxx
PWGHF/hfe/AliAnalysisTaskHFECal.h

index 1ec49318769528048b3846ff28b09e7ac02e275c..ea4fae2e29ead12dcbbb01fd6d218cb49eb58e86 100644 (file)
@@ -104,6 +104,7 @@ AliAnalysisTaskHFECal::AliAnalysisTaskHFECal(const char *name)
   ,fInvmassCut(0.01)   
   ,fNoEvents(0)
   ,fEMCAccE(0)
+  ,hEMCAccE(0)
   ,fTrkpt(0)
   ,fTrkEovPBef(0)       
   ,fTrkEovPAft(0)      
@@ -221,6 +222,7 @@ AliAnalysisTaskHFECal::AliAnalysisTaskHFECal()
   ,fInvmassCut(0.01)   
   ,fNoEvents(0)
   ,fEMCAccE(0)
+  ,hEMCAccE(0)
   ,fTrkpt(0)
   ,fTrkEovPBef(0)       
   ,fTrkEovPAft(0)       
@@ -434,7 +436,7 @@ void AliAnalysisTaskHFECal::UserExec(Option_t*)
   if(TMath::Abs(pVtxZ)>10) return;
   fNoEvents->Fill(1);
   
-  if(fNOtrks<2) return;
+  if(fNOtrks<1) return;
   fNoEvents->Fill(2);
   
   AliPIDResponse *pidResponse = fInputHandler->GetPIDResponse();
@@ -472,9 +474,10 @@ void AliAnalysisTaskHFECal::UserExec(Option_t*)
          double emceta = clustpos.Eta();
          double calInfo[5];
          calInfo[0] = emcphi; calInfo[1] = emceta; calInfo[2] = clustE; calInfo[3] = cent; calInfo[4] = clust->Chi2(); 
-         fEMCAccE->Fill(calInfo); 
-         //if(clustE>1.5)fEMCAccE->Fill(calInfo); 
+         //fEMCAccE->Fill(calInfo); 
+         //if(clustE>3.0)fEMCAccE->Fill(calInfo); 
          //if(fqahist==1 && clustE>1.5)fEMCAccE->Fill(calInfo); 
+         hEMCAccE->Fill(cent,clustE); 
         }
    }
 
@@ -925,6 +928,9 @@ void AliAnalysisTaskHFECal::UserCreateOutputObjects()
   fEMCAccE = new THnSparseD("fEMCAccE","EMC acceptance & E;#eta;#phi;Energy;Centrality;trugCondition;",5,binsE,xminE,xmaxE);
   if(fqahist==1)fOutputList->Add(fEMCAccE);
 
+  hEMCAccE = new TH2F("hEMCAccE","Cluster Energy",200,0,100,100,0,20);
+  fOutputList->Add(hEMCAccE);
+
   fTrkpt = new TH1F("fTrkpt","track pt",100,0,50);
   fOutputList->Add(fTrkpt);
   
index 5b301f4dcd8c8de687caaca0a77c55cf78a8da4e..977b1b7ccfcaab55a6ec0c333ef3d9ef18d69b10 100644 (file)
@@ -91,6 +91,7 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE {
 
   TH1F                 *fNoEvents;             //! no of events
   THnSparseD           *fEMCAccE;              //! EMC acc
+  TH2F                 *hEMCAccE;              //! EMC acc
   TH1F                 *fTrkpt;                //! track pt
   TH2F                 *fTrkEovPBef;           //! track E/p before HFE pid
   TH2F                 *fTrkEovPAft;           //! track E/p after HFE pid