STEER/AliAODCaloCells.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   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 /* $Id$ */
  17
  18 //-------------------------------------------------------------------------
  19 //     AOD class to store calorimeter cell data
  20 //     Author: Markus Oldenburg, CERN
  21 //-------------------------------------------------------------------------
  22
  23 #include "AliAODCaloCells.h"
  24
  25 ClassImp(AliAODCaloCells)
  26
  27 AliAODCaloCells::AliAODCaloCells() : AliVCaloCells(), fNCells(0), fCellNumber(0), fAmplitude(0), fIsSorted(kTRUE), fType(kUndef)
  28 {
  29   // default constructor
  30 }
  31
  32 AliAODCaloCells::AliAODCaloCells(const char* name, const char* title, VCells_t ttype) :
  33     AliVCaloCells(name, title), fNCells(0), fCellNumber(0), fAmplitude(0), fIsSorted(kTRUE), fType(ttype)
  34 {
  35   //constructor
  36 }
  37
  38 AliAODCaloCells::AliAODCaloCells(const AliAODCaloCells& cells) :
  39     AliVCaloCells(cells),
  40     fNCells(cells.fNCells),
  41     fCellNumber(0),
  42     fAmplitude(0),
  43     fIsSorted(cells.fIsSorted),
  44     fType(cells.fType)
  45 {
  46 // Copy constructor
  47   fCellNumber = new Short_t[fNCells];
  48   fAmplitude  = new Double32_t[fNCells];
  49
  50   for (Int_t i = 0; i < fNCells; i++) {
  51     fCellNumber[i]    = cells.fCellNumber[i];
  52     fAmplitude[i]     = cells.fAmplitude[i];
  53   }
  54 }
  55
  56 AliAODCaloCells& AliAODCaloCells::operator=(const AliAODCaloCells& cells)
  57 {
  58     // Assignment operator
  59   if(&cells == this) return *this;
  60   fNCells = cells.fNCells;
  61   for (Int_t i = 0; i < fNCells; i++) {
  62     fCellNumber[i]    = cells.fCellNumber[i];
  63     fAmplitude[i]     = cells.fAmplitude[i];
  64   }
  65   SetName(cells.GetName()) ;
  66   SetTitle(cells.GetTitle()) ;
  67   return *this;
  68 }
  69
  70 AliAODCaloCells::~AliAODCaloCells()
  71 {
  72   // destructor
  73
  74   DeleteContainer();
  75 }
  76
  77 void AliAODCaloCells::Clear(const Option_t*)
  78 {
  79   // clear
  80
  81   DeleteContainer();
  82 }
  83
  84 void AliAODCaloCells::Copy(TObject &obj) const {
  85
  86   // this overwrites the virtual TOBject::Copy()
  87   // to allow run time copying without casting
  88   // in AliESDEvent
  89
  90   if(this==&obj)return;
  91   AliAODCaloCells *robj = dynamic_cast<AliAODCaloCells*>(&obj);
  92   if(!robj)return; // not an AliAODCaloCells
  93   *robj = *this;
  94
  95 }
  96
  97 AliVCaloCells *AliAODCaloCells::CopyCaloCells(Bool_t all = kTRUE) const {
  98
  99   // copy the calo cells into a new object. If option all=FALSE, just the object type,
 100   // for mixing
 101
 102   AliVCaloCells *obj =  new AliAODCaloCells();
 103
 104   if(all){
 105     obj->SetName (GetName()) ;
 106     obj->SetTitle(GetTitle()) ;
 107     obj->SetType (GetType()) ;
 108
 109     obj->SetNumberOfCells(fNCells);
 110     for (Short_t i = 0; i < fNCells; i++)
 111       obj->SetCell(i,fCellNumber[i],fAmplitude[i],-1);
 112   }
 113
 114   return obj;
 115
 116 }
 117
 118 void AliAODCaloCells::CreateContainer(Short_t nCells)
 119 {
 120   // function that creates container to store calorimeter cell data
 121
 122   DeleteContainer();
 123
 124   if (nCells <= 0) {
 125     fNCells = 0;
 126     return;
 127   }
 128
 129   fNCells = nCells;
 130
 131   fCellNumber = new Short_t[fNCells];
 132   fAmplitude  = new Double32_t[fNCells];
 133
 134   // set to zero
 135   for(int i = 0;i<fNCells;++i){
 136     fAmplitude[i] = fCellNumber[i] = 0 ;
 137   }
 138 }
 139
 140 void AliAODCaloCells::DeleteContainer()
 141 {
 142   // deletes allocated memory
 143
 144   if (fCellNumber)
 145   {
 146     delete[] fCellNumber;
 147     fCellNumber = NULL;
 148   }
 149
 150   if (fAmplitude)
 151   {
 152     delete[] fAmplitude;
 153     fAmplitude = NULL;
 154   }
 155
 156   fNCells = 0;
 157   fIsSorted = kFALSE;
 158 }
 159
 160 void AliAODCaloCells::Sort()
 161 {
 162   // sort the cell array by cell number
 163
 164   Int_t *idxArray = new Int_t[fNCells];
 165   TMath::Sort(fNCells,fCellNumber,idxArray,kFALSE);
 166
 167   Short_t    *newIndex     = new Short_t[fNCells];
 168   Double32_t *newAmplitude = new Double32_t[fNCells];
 169   for (Int_t i=0; i < fNCells; i++) {
 170     newIndex[i]     = fCellNumber[idxArray[i]];
 171     newAmplitude[i] = fAmplitude[idxArray[i]];
 172   }
 173   delete [] fCellNumber;
 174   delete [] fAmplitude;
 175   fCellNumber = newIndex;
 176   fAmplitude  = newAmplitude;
 177
 178   delete [] idxArray;
 179
 180   fIsSorted = kTRUE;
 181 }
 182
 183 Bool_t AliAODCaloCells::SetCell(Short_t pos, Short_t cellNumber, Double32_t amplitude, Double32_t /*time*/)
 184 {
 185   // Sets a cell at the given position
 186
 187   if (pos>=0 && pos < fNCells) {
 188     fCellNumber[pos] = cellNumber;
 189     fAmplitude[pos]  = amplitude;
 190     fIsSorted = kFALSE;
 191     return kTRUE;
 192   } else {
 193     return kFALSE;
 194   }
 195 }