emcal cells pos vs. energy (to check bad towers)
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Sep 2012 10:59:10 +0000 (10:59 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Sep 2012 10:59:10 +0000 (10:59 +0000)
PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSAQA.cxx
PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSAQA.h

index fe46d53..ac3fc51 100644 (file)
@@ -48,7 +48,7 @@ AliAnalysisTaskSAQA::AliAnalysisTaskSAQA() :
   fHistDeltaPhiPt(0),
   fHistNCellsEnergy(0),
   fHistClusTimeEnergy(0),
-  fHistCellsEnergy(0),
+  fHistCellsAbsIdEnergy(0),
   fHistChVSneCells(0),
   fHistChVSneClus(0),
   fHistChVSneCorrCells(0)
@@ -90,7 +90,7 @@ AliAnalysisTaskSAQA::AliAnalysisTaskSAQA(const char *name) :
   fHistDeltaPhiPt(0),
   fHistNCellsEnergy(0),
   fHistClusTimeEnergy(0),
-  fHistCellsEnergy(0),
+  fHistCellsAbsIdEnergy(0),
   fHistChVSneCells(0),
   fHistChVSneClus(0),
   fHistChVSneCorrCells(0)
@@ -232,10 +232,11 @@ void AliAnalysisTaskSAQA::UserCreateOutputObjects()
     fHistNCellsEnergy->GetYaxis()->SetTitle("N_{cells}");
     fOutput->Add(fHistNCellsEnergy); 
      
-    fHistCellsEnergy = new TH1F("fHistCellsEnergy","Energy spectrum of cells", fNbins, fMinBinPt, fMaxBinPt);
-    fHistCellsEnergy->GetXaxis()->SetTitle("Energy (GeV)");
-    fHistCellsEnergy->GetYaxis()->SetTitle("counts");
-    fOutput->Add(fHistCellsEnergy);
+    fHistCellsAbsIdEnergy = new TH2F("fHistCellsAbsIdEnergy","fHistCellsAbsIdEnergy", 11600,0,11599,(Int_t)(fNbins / 2), fMinBinPt, fMaxBinPt / 2);
+    fHistCellsAbsIdEnergy->GetXaxis()->SetTitle("cell abs. Id");
+    fHistCellsAbsIdEnergy->GetYaxis()->SetTitle("Energy (GeV)");
+    fHistCellsAbsIdEnergy->GetZaxis()->SetTitle("counts");    
+    fOutput->Add(fHistCellsAbsIdEnergy);
     
     fHistChVSneCells = new TH2F("fHistChVSneCells","Charged energy vs. neutral (cells) energy", 
                                (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5, (Int_t)(fNbins * 2.5), fMinBinPt, fMaxBinPt * 2.5);
@@ -363,8 +364,9 @@ Int_t AliAnalysisTaskSAQA::DoCellLoop(Float_t &sum, Float_t &sum_cut)
   const Int_t ncells = cells->GetNumberOfCells();
 
   for (Int_t pos = 0; pos < ncells; pos++) {
-    Float_t amp = cells->GetAmplitude(pos);
-    fHistCellsEnergy->Fill(amp);
+    Float_t amp   = cells->GetAmplitude(pos);
+    Int_t   absId = cells->GetCellNumber(pos);
+    fHistCellsAbsIdEnergy->Fill(absId,amp);
     sum += amp;
     if (amp < fCellEnergyCut)
       continue;
index a40766f..79500b8 100644 (file)
@@ -62,7 +62,7 @@ class AliAnalysisTaskSAQA : public AliAnalysisTaskEmcalJet {
   TH2F                       *fHistJetsPtArea[4];        //!Pt vs. area of jets
 
   // EMCAL Cells
-  TH1F                       *fHistCellsEnergy;          //!Energy spectrum of cells
+  TH2F                       *fHistCellsAbsIdEnergy;    //!Energy spectrum of cells
 
   // Had corr QA
   TH2F                       *fHistChVSneCells;          //!Charged vs. neutral (cells) energy