added time information

diff --git a/PWGHF/hfe/AliAnalysisTaskHFEemcQA.cxx b/PWGHF/hfe/AliAnalysisTaskHFEemcQA.cxx
index 1a06ef65b0feaefb1fd81920b0cab68ad2378b86..58eb77de97dd9f7c438211f8d6b29faabd0eb949 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFEemcQA.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskHFEemcQA.cxx
@@ -78,6 +78,8 @@ ClassImp(AliAnalysisTaskHFEemcQA)
   fHistoNClsE3(0),
   fHistoNCells(0),
   fHistoCalCell(0),
+  fHistoCalTime(0),
+  fHistoCalCellTime(0),
   fNegTrkIDPt(0),
   fTrkPt(0),
   fTrketa(0),
@@ -147,6 +149,8 @@ AliAnalysisTaskHFEemcQA::AliAnalysisTaskHFEemcQA()
   fHistoNClsE3(0),
   fHistoNCells(0),
   fHistoCalCell(0),
+  fHistoCalTime(0),
+  fHistoCalCellTime(0),
   fNegTrkIDPt(0),
   fTrkPt(0),
   fTrketa(0),
@@ -271,6 +275,12 @@ void AliAnalysisTaskHFEemcQA::UserCreateOutputObjects()
   fHistoCalCell = new TH2F("fHistoCalCell","EMCAL cells in a cluster;cell ID;E (GeV)",15000,-0.5,14999.5,600,0,30);
   fOutputList->Add(fHistoCalCell);
 
+  fHistoCalTime = new TH2F("fHistoCalTime","Cluster Time;  p_{T} (GeV/c); Time (s)",300,0,30,1000,1e-8,1e-5);
+  fOutputList->Add(fHistoCalTime);
+
+  fHistoCalCellTime = new TH2F("fHistoCalCellTime","Cell Time;  ID; Time (s)",15000,-0.5,14999.5,1000,1e-8,1e-5);
+  fOutputList->Add(fHistoCalCellTime);
+
   fNegTrkIDPt = new TH1F("fNegTrkIDPt", "p_{T} distribution of tracks with negative track id;p_{T} (GeV/c);counts", 500, 0.0, 50.0); 
   fOutputList->Add(fNegTrkIDPt);
 
@@ -526,6 +536,9 @@ void AliAnalysisTaskHFEemcQA::UserExec(Option_t *)
       fEMCClsEtaPhi->Fill(emceta,emcphi);
       fHistoNCells->Fill(clustE,clust->GetNCells());
       //fHistoNCells->Fill(clust->GetNCells());
+    
+      double emctof = clust->GetTOF();
+      fHistoCalTime->Fill(clustE,emctof);
 
       NclustAll++;  
       if(clustE>0.1)NclustE1++;
@@ -546,14 +559,17 @@ void AliAnalysisTaskHFEemcQA::UserExec(Option_t *)
   Short_t cellAddr, nSACell;
   Int_t  mclabel;
   Short_t iSACell;
-  Double_t cellAmp=-1., cellTimeT=-1., clusterTime=-1., efrac=-1.;
+  Double_t cellAmp=-1., cellTimeT=-1., efrac=-1.;
 
   nSACell = fCaloCells->GetNumberOfCells();
   for(iSACell = 0; iSACell < nSACell; iSACell++ ){ 
     Bool_t haveCell = fCaloCells->GetCell(iSACell, cellAddr, cellAmp, cellTimeT , mclabel, efrac);
     //virtual Bool_t   GetCell(Short_t pos, Short_t &cellNumber, Double_t &amplitude, Double_t &time, Int_t &mclabel,    Double_t  &efrac)      
-    if(haveCell)fHistoCalCell->Fill(cellAddr,cellAmp);
-
+    if(haveCell)
+       {
+        fHistoCalCell->Fill(cellAddr,cellAmp);
+        fHistoCalCellTime->Fill(cellAddr,cellTimeT);
+        }
   }
 
   /////////////////////////////////

diff --git a/PWGHF/hfe/AliAnalysisTaskHFEemcQA.h b/PWGHF/hfe/AliAnalysisTaskHFEemcQA.h
index 4eb10d4308faf744b4e7d3392211b3473ffeb4e7..fedf97c8b8f49ab82a8ca18b5bb5fd968e5ce3d6 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFEemcQA.h
+++ b/PWGHF/hfe/AliAnalysisTaskHFEemcQA.h
@@ -66,6 +66,8 @@ class AliAnalysisTaskHFEemcQA : public AliAnalysisTaskSE {
     //TH1F        *fHistoNCells;//! No of cells per cluster
     TH2F        *fHistoNCells;//! No of cells per cluster
     TH2F        *fHistoCalCell;//! No of cells per cluster
+    TH2F        *fHistoCalTime;
+    TH2F        *fHistoCalCellTime;
     TH1F        *fNegTrkIDPt;//!neg track ID
     TH1F        *fTrkPt;//!track pt
     TH1F        *fTrketa;//!track eta