/* $Id$ */
//_________________________________________________________________________
-// EMCAL digit: Id
-// energy
-// 3 identifiers for the primary particle(s) at the origine of the digit
-// The digits are made in FinishEvent() by summing all the hits in a single EMCAL crystal or PPSD gas cell
-// It would be nice to replace the 3 identifiers by an array, but, because digits are kept in a TClonesQArray,
-// it is not possible to stream such an array... (beyond my understqnding!)
+// EMCAL digit:
+//
+// A Digit is the sum of energy in a Tower (Hit sum) and stores information, about primaries
+// and entering particle contributing to a Digit
//
-//*-- Author: Laurent Aphecetche & Yves Schutz (SUBATECH)
+//*-- Author: Sahal Yacoob (LBL)
+// based on : AliPHOSDigit
+//___________________________________________________________________________
// --- ROOT system ---
// --- Standard library ---
// --- AliRoot header files ---
-
#include "AliDigitNew.h"
class AliEMCALDigit : public AliDigitNew {
public:
AliEMCALDigit() ;
- AliEMCALDigit(Int_t primary, Int_t iparent, Int_t id, Int_t DigEnergy, Int_t index = -1) ;
+ AliEMCALDigit(Int_t primary, Int_t iparent, Int_t id, Int_t DigEnergy, Float_t Time, Int_t index = -1, Float_t dE = 0) ;
AliEMCALDigit(const AliEMCALDigit & digit) ;
virtual ~AliEMCALDigit() ;
Bool_t operator==(const AliEMCALDigit &rValue) const;
- AliEMCALDigit& operator+(AliEMCALDigit const &rValue) ;
-
+ AliEMCALDigit operator+(const AliEMCALDigit &rValue) ;
+ AliEMCALDigit operator*(Float_t factor) ;
+ const AliEMCALDigit& operator = (const AliEMCALDigit &) {return *this;}
+
Int_t Compare(const TObject * obj) const ;
- Int_t GetNprimary() const {
- // returns the number of primaries
- return fNprimary ; }
+ 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 ;
- Bool_t IsSortable() const {
- // says that AliEMCALDigits are sortable (needed for Sort method
- return kTRUE ; }
- void SetAmp(Int_t Amp) {
- // sets the amplitude data member
- fAmp=Amp ; }
- void ShiftPrimary(Int_t shift); // shift to semarate different TreeK in merging
-
- private:
+ 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 ;}
+ Bool_t IsSortable() const { return kTRUE ; }
+ void SetAmp(Int_t amp) { fAmp= amp ; }
+ void SetId(Int_t id) {fId = id ;}
+ void SetTime(Float_t time) {fTime = time ;}
+ void SetTimeR(Float_t time) {fTimeR = time ;}
+ void ShiftPrimary(Int_t shift); // shift to separate different TreeK in merging
+
+ private:
Int_t fNprimary ; // Number of primaries
- Int_t fNMaxPrimary ; //! Max Number of primaries
- Int_t fPrimary[5] ; // Array of primaries
+ Int_t fNMaxPrimary ; // Max Number of primaries
+ Int_t *fPrimary ; //[fNMaxPrimary] Array of primaries
+ Float_t *fDEPrimary; //[fNMaxPrimary] Array of primary energy contributions
Int_t fNiparent ; // Number of initial parents
- Int_t fNMaxiparent ; //! Max Number of parents
- Int_t fIparent[50] ; // Array of parents
-
- ClassDef(AliEMCALDigit,1) // Digit in EMCAL
+ Int_t fNMaxiparent ; // Max Number of parents
+ Int_t *fIparent ; //[fNMaxiparent] Array of parents
+ Float_t *fDEParent; //[fNMaxiparent] Array of parent energy contributions
+ Int_t fMaxIter ; // Number to Increment Maxiparent, and MaxPrimary if default is not sufficient
+ Float_t fTime ; // Calculated time
+ Float_t fTimeR ; // Earliest time: to be used by Digits2Raw
+
+ ClassDef(AliEMCALDigit,3) // Digit in EMCAL
} ;