AliEMCALClusterizerv3.h

Go to the documentation of this file.

#ifndef ALIEMCALCLUSTERIZERV3_H
#define ALIEMCALCLUSTERIZERV3_H

/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */


#include "AliEMCALClusterizerv1.h"
class AliEMCALRecPoint; 
class AliEMCALDigit;

namespace EMCALClusterFinder
{
  // Define numbers rows/columns for topological representation of cells
  const UInt_t kNrows    = (24+1)*(6+4);  // 10x supermodule rows (6 for EMCAL, 4 for DCAL). +1 accounts for topological gap between two supermodules
  const UInt_t kNcolumns = 48*2+1;        // 2x  supermodule columns + 1 empty space in between for DCAL (not used for EMCAL)
}

//_________________________________________________________________________
//_________________________________________________________________________
class AliEMCALClusterFinder: public TObject
{
  struct cellWithE {
    cellWithE() : energy(0.), row(0), column(0) {}
    cellWithE(Float_t e, Int_t r, Int_t c) : energy(e), row(r), column(c) {}
    // std::sort will require operator< to compile.
    bool operator<( cellWithE const& rhs ) const
       { return energy < rhs.energy; }
    Float_t energy;
    Int_t row;
    Int_t column;
  };

  public:
    AliEMCALClusterFinder(TObjArray* outputArray, AliEMCALGeometry* geometry, Double_t timeCut, Double_t timeMin, Double_t timeMax, Double_t gradientCut, Bool_t doEnergyGradientCut, Double_t thresholdSeedE, Double_t thresholdCellE);
    ~AliEMCALClusterFinder();

    Int_t               FindClusters(TClonesArray* digits);
    TObjArray*          GetFoundClusters() {return fFoundClusters;}
    
  private:
    AliEMCALRecPoint*   GetClusterFromNeighbours(AliEMCALRecPoint* recPoint, Int_t row, Int_t column);
    void                GetTopologicalRowColumn(AliEMCALDigit* digit, Int_t& row, Int_t& column);

    AliEMCALGeometry*   fEMCALGeometry;                   
    cellWithE           fSeedList[EMCALClusterFinder::kNrows*EMCALClusterFinder::kNcolumns];      
    AliEMCALDigit*      fDigitMap[EMCALClusterFinder::kNrows][EMCALClusterFinder::kNcolumns];     
    Bool_t              fCellMask[EMCALClusterFinder::kNrows][EMCALClusterFinder::kNcolumns];     

    TObjArray*          fFoundClusters;                   
    Int_t               fNumFoundClusters;                

    Double_t            fTimeCut;                         
    Double_t            fTimeMin;                         
    Double_t            fTimeMax;                         
    Double_t            fGradientCut;                     
    Bool_t              fDoEnergyGradientCut;             
    Double_t            fThresholdSeedEnergy;             
    Double_t            fThresholdCellEnergy;             

    AliEMCALClusterFinder(const AliEMCALClusterFinder&);            // not implemented
    AliEMCALClusterFinder &operator=(const AliEMCALClusterFinder&); // not implemented

  ClassDef(AliEMCALClusterFinder,1);
};


//_________________________________________________________________________
//_________________________________________________________________________
class AliEMCALClusterizerv3 : public AliEMCALClusterizerv1
{
  
public:
  AliEMCALClusterizerv3() ;         
  AliEMCALClusterizerv3(AliEMCALGeometry* geometry);
  AliEMCALClusterizerv3(AliEMCALGeometry* geometry, AliEMCALCalibData* calib, 
                        AliEMCALCalibTime * calibt, AliCaloCalibPedestal *pedestal);

  virtual ~AliEMCALClusterizerv3();

 
  virtual const char *Version() const { return "clu-v3";}  

  void                SetDoEnGradCut(Bool_t b) { fDoEnGradCut = b; }

protected:
  virtual void        MakeClusters();            

  Bool_t                  fDoEnGradCut;             
  AliEMCALClusterFinder*  fClusterFinder;           

private:

  AliEMCALClusterizerv3              (const AliEMCALClusterizerv3 &); //copy ctor
  AliEMCALClusterizerv3 & operator = (const AliEMCALClusterizerv3 &);

  ClassDef(AliEMCALClusterizerv3,1);

};

#endif // AliEMCALCLUSTERIZERV3_H