AliEMCALDigit.h

Go to the documentation of this file.

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


// --- ROOT system ---
#include "TObject.h" 

// --- AliRoot header files ---
#include "AliDigitNew.h"

using std::ostream;

class AliEMCALDigit : public AliDigitNew 
{

  friend ostream& operator << ( ostream& , const AliEMCALDigit&) ;

 public:
  
  AliEMCALDigit() ;
  AliEMCALDigit(Int_t primary, Int_t iparent, Int_t id, Float_t digEnergy, 
                Float_t time, Int_t type,Int_t index = -1, 
                Float_t chi2=0, Int_t ndf=0, Float_t dE = 0) ;
  
  AliEMCALDigit(const AliEMCALDigit & digit) ;
  
  virtual ~AliEMCALDigit() ;

  Bool_t operator==(const AliEMCALDigit &rValue) const;
  AliEMCALDigit operator+(const AliEMCALDigit &rValue) ;
  AliEMCALDigit operator*(Float_t factor) ; 
  AliEMCALDigit& operator = (const AliEMCALDigit & digit) ;
  
  enum  digitType{kUnknown=-1, kHG=0, kLG=1, kLGnoHG=2, kTrigger=3, kEmbedded = 4};

  void     Clear(Option_t*) ; 
  Int_t    Compare(const TObject * obj) const ;
  Float_t  GetAmplitude()   const { if(!fAmp)return fAmpFloat ; else return fAmp ;}//Keep backward compatibility.
  Float_t  GetEta()         const ; 
  Int_t    GetNprimary()    const { return fNprimary ;}
  Int_t    GetPrimary(Int_t index)   const ; 
  Float_t  GetDEPrimary(Int_t index) const ; 
  Int_t    GetNiparent()    const { return fNiparent ;}
  Int_t    GetIparent(Int_t index)   const ;
  Float_t  GetDEParent(Int_t index)  const ; 
  Float_t  GetPhi()         const ;
  Float_t  GetTime(void)    const { return fTime  ;}
  Float_t  GetTimeR(void)   const { return fTimeR ;}
  Float_t  GetChi2(void)    const { return fChi2  ;}
  Int_t    GetNDF(void)     const { return fNDF   ;}
  Bool_t   IsSortable()     const { return kTRUE  ;}
  Int_t    GetType()        const { return fDigitType ;}
  
  void     SetAmp(Int_t amp)         { fAmp       = amp  ; } //old
  void     SetAmplitude(Float_t amp) { fAmpFloat  = amp  ; }
  void     SetId(Int_t idt)          { fId        = idt  ; }
  void     SetTime(Float_t time)     { fTime      = time ; }
  void     SetTimeR(Float_t time)    { fTimeR     = time ; }
  void     SetChi2(Float_t chi)      { fChi2      = chi  ; }
  void     SetNDF(Int_t ndf)         { fNDF       = ndf  ; }
  void     SetType(Int_t t)          { fDigitType = t    ; }
  void     ShiftPrimary(Int_t shift); // shift to separate different TreeK in merging

  //
  // Raw time sample
  //
  
  // ALTRO
  Int_t    GetNALTROSamplesLG() const {if(fDigitType==kLG)      return fNSamples;   else return 0 ; }
  Bool_t   GetALTROSampleLG(const Int_t iSample, Int_t& timeBin, Int_t& amp) const;
  Int_t    GetNALTROSamplesHG() const {if(fDigitType==kHG)      return fNSamplesHG; else return 0 ; }
  Bool_t   GetALTROSampleHG(const Int_t iSample, Int_t& timeBin, Int_t& amp) const;
  
  // FALTRO, trigger. Same data members as Low Gain 
  Int_t    GetNFALTROSamples()  const {if(fDigitType==kTrigger) return fNSamples;   else return 0 ; }
  Bool_t   GetFALTROSample(const Int_t iSample, Int_t& timeBin, Int_t& amp)  const ;
  
  void     SetALTROSamplesHG (const Int_t nSamplesHG, Int_t *samplesHG);
  void     SetALTROSamplesLG (const Int_t nSamplesLG, Int_t *samplesLG);
  void     SetFALTROSamples  (const Int_t nSamples,   Int_t *samples) 
  { if(fDigitType==kTrigger) SetALTROSamplesLG(nSamples, samples) ; } 

  //
  // Primary/Parents array creation
  // Used at analysis level while reclusterizing
  //
  void SetListOfPrimaries(Int_t npri, Int_t * prilist, Float_t * edepList) ;  
  void SetListOfParents  (Int_t npar, Int_t * parlist, Float_t * edepList) ;  
  
  //
  // Other
  //
  void     SetCalibAmp(Float_t amp) { fAmpCalib = amp  ; }
  Double_t GetCalibAmp()   const    { return fAmpCalib ; }

  void     Print(const Option_t* /*opt*/) const;
  
 private: 
  
  Float_t  fAmpFloat;     
  
  Int_t    fNSamples;     
  
  Int_t   *fSamples;      //[fNSamples]  
  
  Int_t    fNSamplesHG;   
  
  Int_t   *fSamplesHG;    //[fNSamples] 
  
  Int_t    fNprimary ;    
  
  Int_t    fNMaxPrimary ; 
  
  Int_t   *fPrimary ;     //[fNMaxPrimary]       
  
  Float_t *fDEPrimary ;   //[fNMaxPrimary]  
    
  Int_t    fNiparent ;    
  
  Int_t    fNMaxiparent ; 
  
  Int_t   *fIparent ;     //[fNMaxiparent]       
  
  Float_t *fDEParent;     //[fNMaxiparent] 
  
  Int_t    fMaxIter  ;    
  
  Float_t  fTime ;        
  
  Float_t  fTimeR ;       
  
  Float_t  fChi2;         
  
  Int_t    fNDF;          
  
  Int_t    fDigitType;    
  
  Float_t  fAmpCalib;     

  ClassDef(AliEMCALDigit,7)  ;

} ;

// inline definitions

//____________________________________________________________________________
inline Int_t AliEMCALDigit::Compare(const TObject * obj) const
{
  AliEMCALDigit * digit = (AliEMCALDigit *)obj ;

  Int_t iddiff = fId - digit->GetId() ;

  if ( iddiff > 0 )
    return 1 ;
  else if ( iddiff < 0 )
    return -1 ;
  else
    return 0 ;
}

#endif //  ALIEMCALDIGIT_H