HLT/global/physics/AliHLTCaloHistoClusterEnergy.cxx

   1 /**************************************************************************
   2  * This file is property of and copyright by the ALICE HLT Project        *
   3  * All rights reserved.                                                   *
   4  *                                                                        *
   5  * Primary Authors: Svein Lindal <slindal@fys.uio.no>                     *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 /**
  17  * @file   AliHLTCaloHistoClusterEnergy
  18  * @author Svein Lindal
  19  * @date
  20  * @brief  Produces histograms of cluster energy distributions
  21  */
  22
  23 // see header file for class documentation
  24 // or
  25 // refer to README to build package
  26 // or
  27 // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
  28
  29 #include "AliHLTCaloHistoClusterEnergy.h"
  30 #include "AliESDCaloCluster.h"
  31 #include "AliHLTCaloClusterDataStruct.h"
  32 #include "TH1F.h"
  33 #include "TH2F.h"
  34 #include "TRefArray.h"
  35
  36 ClassImp(AliHLTCaloHistoClusterEnergy);
  37
  38 AliHLTCaloHistoClusterEnergy::AliHLTCaloHistoClusterEnergy(TString det) :
  39   AliHLTCaloHistoProducer(),
  40   fHistClusterEnergy(NULL),
  41   fHistClusterEnergyVsNCells(NULL)
  42 {
  43   // See header file for documentation
  44   fHistClusterEnergy = new TH1F(Form("%s fHistClusterEnergy", det.Data()), Form("%s Distribution of total energy in clusters", det.Data()), 5000, 0, 100);
  45   fHistClusterEnergy->GetXaxis()->SetTitle("E GeV");
  46   fHistClusterEnergy->GetYaxis()->SetTitle("Number of counts");
  47   fHistArray->AddLast(fHistClusterEnergy);
  48
  49   fHistClusterEnergyVsNCells = new TH2F(Form("%s fHistClusterEnergyVsNCells", det.Data()), Form("%s Distribution of Energy vs Number of Cells in cluster", det.Data()), 1000, 0, 100, 50, 0 , 50);
  50   fHistClusterEnergyVsNCells->GetXaxis()->SetTitle("Energy in cluster (GeV)");
  51   fHistClusterEnergyVsNCells->GetYaxis()->SetTitle("Number of Cells in cluster");
  52   fHistArray->AddLast(fHistClusterEnergyVsNCells);
  53
  54 }
  55
  56 AliHLTCaloHistoClusterEnergy::~AliHLTCaloHistoClusterEnergy()
  57 {
  58   //destructor
  59   if(fHistClusterEnergy)
  60     delete fHistClusterEnergy;
  61   fHistClusterEnergy = NULL;
  62
  63   if(fHistClusterEnergyVsNCells)
  64     delete fHistClusterEnergyVsNCells;
  65   fHistClusterEnergyVsNCells = NULL;
  66 }
  67
  68 Int_t AliHLTCaloHistoClusterEnergy::FillHistograms(Int_t nc, TRefArray * clusterArray) {
  69   //See header file for documentation
  70
  71   for(int ic = 0; ic < nc; ic++) {
  72     AliESDCaloCluster * cluster = static_cast<AliESDCaloCluster*>(clusterArray->At(ic));
  73     FillClusterEnergyHistos(cluster);
  74   }
  75   return 0;
  76 }
  77
  78 Int_t AliHLTCaloHistoClusterEnergy::FillHistograms(Int_t nc, vector<AliHLTCaloClusterDataStruct*> &cVec) {
  79   //See header file for documentation
  80   // HLTInfo("histo");
  81   for(int ic = 0; ic < nc; ic++) {
  82     AliHLTCaloClusterDataStruct * cluster = cVec.at(ic);
  83     FillClusterEnergyHistos(cluster);
  84   }
  85   return 0;
  86 }
  87
  88 template <class T>
  89 Int_t AliHLTCaloHistoClusterEnergy::FillClusterEnergyHistos(T* cluster) {
  90   fHistClusterEnergy->Fill(cluster->E());
  91   fHistClusterEnergyVsNCells->Fill(cluster->E(), cluster->GetNCells());
  92   return 0;
  93 }