PWGLF/FORWARD/analysis2/AliCentraldNdetaTask.cxx

   1 //====================================================================
   2 #include "AliCentraldNdetaTask.h"
   3 #include <TMath.h>
   4 #include <TH2D.h>
   5 #include <TH1D.h>
   6 #include <THStack.h>
   7 #include <TList.h>
   8 #include <AliAnalysisManager.h>
   9 #include <AliAODEvent.h>
  10 #include <AliAODHandler.h>
  11 #include <AliAODInputHandler.h>
  12 #include "AliForwardUtil.h"
  13 #include "AliAODForwardMult.h"
  14 #include "AliAODCentralMult.h"
  15
  16 ClassImp(AliCentraldNdetaTask)
  17 #ifdef DOXY_INPUT
  18 ;
  19 #endif
  20
  21 //____________________________________________________________________
  22 AliCentraldNdetaTask::AliCentraldNdetaTask()
  23   : AliBasedNdetaTask()
  24 {
  25   DGUARD(fDebug,3,"Default CTOR of AliCentraldNdetaTask");
  26 }
  27
  28 //____________________________________________________________________
  29 AliCentraldNdetaTask::AliCentraldNdetaTask(const char*)
  30   : AliBasedNdetaTask("Central")
  31 {
  32   DGUARD(fDebug,3,"Named CTOR of AliCentraldNdetaTask");
  33   fSymmetrice = false;
  34   fCorrEmpty  = false;
  35   SetTitle("Central");
  36 }
  37
  38 //____________________________________________________________________
  39 TH2D*
  40 AliCentraldNdetaTask::GetHistogram(const AliAODEvent* aod, Bool_t mc)
  41 {
  42   // Get objects from the event structure
  43   DGUARD(fDebug,2,"Get our histogram for AliCentraldNdetaTask");
  44   TObject* obj = 0;
  45   if (mc) obj = aod->FindListObject("CentralClustersMC");
  46   else    obj = aod->FindListObject("CentralClusters");
  47
  48   // We should have a central object at least
  49   if (!obj) {
  50     if (!mc) AliWarning("No Central object found AOD");
  51     return 0;
  52   }
  53
  54   // Cast to good types
  55   AliAODCentralMult* central   = static_cast<AliAODCentralMult*>(obj);
  56
  57   return &(central->GetHistogram());
  58 }
  59
  60 //
  61 // EOF
  62 //