* provided "as is" without express or implied warranty. *
**************************************************************************/
-/* $Id$ */
+/* $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 the energy lost in an EMCAL Tower
+// It also stores information on Primary, and enterring particle
+// tracknumbers Digits are created using AliEMCALSDigitizer, followed
+// by AliEMCALDigitizer
//
-//*-- Author: Laurent Aphecetche & Yves Schutz (SUBATECH)
-
+//*-- Author: Sahal Yacoob (LBL)
+// based on : AliPHOSDigit
+//__________________________________________________________________________
// --- ROOT system ---
// --- Standard library ---
-#include <iostream.h>
+#include <Riostream.h>
// --- AliRoot header files ---
#include "AliEMCALDigit.h"
+#include "AliEMCALGeometry.h"
+#include "AliEMCALGetter.h"
ClassImp(AliEMCALDigit)
//____________________________________________________________________________
- AliEMCALDigit::AliEMCALDigit()
+ AliEMCALDigit::AliEMCALDigit()
{
// default ctor
fIndexInList = -1 ;
fNprimary = 0 ;
fNMaxPrimary = 5 ;
- fNiparent = 0 ;
- fNMaxiparent = fNMaxPrimary*10;
+ fNiparent = 0 ;
+ fNMaxiparent = 5;
+ fPrimary = 0 ;
+ fIparent = 0 ;
+ fMaxIter = 0;
}
//____________________________________________________________________________
-AliEMCALDigit::AliEMCALDigit(Int_t primary, Int_t iparent, Int_t id, Int_t DigEnergy, Int_t index)
+AliEMCALDigit::AliEMCALDigit(Int_t primary, Int_t iparent, Int_t id, Int_t DigEnergy, Float_t time, Int_t index)
{
// ctor with all data
- fNMaxPrimary = 5 ;
- fNMaxiparent = fNMaxPrimary*10;
+ fNMaxPrimary = 25 ;
+ fNMaxiparent = 150 ;
+ fPrimary = new Int_t[fNMaxPrimary] ;
+ fIparent = new Int_t[fNMaxiparent] ;
fAmp = DigEnergy ;
+ fTime = time ;
fId = id ;
fIndexInList = index ;
+ fMaxIter = 5;
if( primary != -1){
fNprimary = 1 ;
- fPrimary[0] = primary ;
+ fPrimary[0] = primary ;
fNiparent = 1 ;
fIparent[0] = iparent ;
-
- }
+}
else{ //If the contribution of this primary smaller than fDigitThreshold (AliEMCALv1)
fNprimary = 0 ;
fPrimary[0] = -1 ;
for ( i = 1; i < fNMaxPrimary ; i++)
fPrimary[i] = -1 ;
- for ( Int_t j =1; j< fNMaxiparent ; j++)
- fIparent[j] = -1 ;
+ for ( i =1; i< fNMaxiparent ; i++)
+ fIparent[i] = -1 ;
}
//____________________________________________________________________________
-AliEMCALDigit::AliEMCALDigit(const AliEMCALDigit & digit)
+AliEMCALDigit::AliEMCALDigit(const AliEMCALDigit & digit) : AliDigitNew(digit)
{
// copy ctor
fNMaxPrimary = digit.fNMaxPrimary ;
fNMaxiparent = digit.fNMaxiparent ;
+ fPrimary = new Int_t[fNMaxPrimary] ;
+ fIparent = new Int_t[fNMaxiparent] ;
Int_t i ;
for ( i = 0; i < fNMaxPrimary ; i++)
fPrimary[i] = digit.fPrimary[i] ;
for (j = 0; j< fNMaxiparent ; j++)
fIparent[j] = digit.fIparent[j] ;
fAmp = digit.fAmp ;
+ fTime = digit.fTime ;
fId = digit.fId;
+ fMaxIter = digit.fMaxIter;
fIndexInList = digit.fIndexInList ;
fNprimary = digit.fNprimary ;
fNiparent = digit.fNiparent ;
AliEMCALDigit::~AliEMCALDigit()
{
// Delete array of primiries if any
-
+ delete [] fPrimary ;
+ delete [] fIparent ;
}
//____________________________________________________________________________
}
+//____________________________________________________________________________
+Float_t AliEMCALDigit::GetEta() const
+{
+ Float_t eta=-10., phi=-10.;
+ Int_t id = GetId();
+ const AliEMCALGeometry *g = AliEMCALGetter::Instance()->EMCALGeometry();
+ g->EtaPhiFromIndex(id,eta,phi);
+ return eta ;
+}
+
+//____________________________________________________________________________
+Float_t AliEMCALDigit::GetPhi() const
+{
+ Float_t eta=-10., phi=-10.;
+ Int_t id = GetId();
+ const AliEMCALGeometry *g = AliEMCALGetter::Instance()->EMCALGeometry();
+ g->EtaPhiFromIndex(id,eta,phi);
+ return phi ;
+}
+
//____________________________________________________________________________
Int_t AliEMCALDigit::GetPrimary(Int_t index) const
{
// retrieves the primary particle number given its index in the list
Int_t rv = -1 ;
- if ( index <= fNprimary ){
+ if ( index <= fNprimary && index > 0){
rv = fPrimary[index-1] ;
}
}
-
//____________________________________________________________________________
void AliEMCALDigit::ShiftPrimary(Int_t shift){
//shifts primary nimber to BIG offset, to separate primary in different TreeK
// if amplitude is larger than
fAmp += digit.fAmp ;
-
+ if(fTime > digit.fTime)
+ fTime = digit.fTime ;
+
Int_t max1 = fNprimary ;
Int_t max2 = fNiparent ;
Int_t index ;
- for (index = 0 ; index < digit.fNprimary ; index++){
+ for (index = 0 ; index < digit.fNprimary ; index++){
Bool_t deja = kTRUE ;
Int_t old ;
for ( old = 0 ; (old < max1) && deja; old++) { //already have this primary?
- if(fPrimary[old] == (digit.fPrimary)[index])
+ if(fPrimary[old] == digit.fPrimary[index])
deja = kFALSE;
}
if(deja){
- fPrimary[fNprimary] = (digit.fPrimary)[index] ;
+ if(max1<fNMaxPrimary){ fPrimary[max1] = digit.fPrimary[index] ;
fNprimary++ ;
- if(fNprimary>fNMaxPrimary) {
- cout << "AliEMCALDigit >> Increase NMaxPrimary "<< endl ;
- return *this ;
+ max1++;}
+ if(fNprimary==fNMaxPrimary) {
+
+ TString mess = " NMaxPrimary = " ;
+ mess += fNMaxPrimary ;
+ mess += " is too small" ;
+ Fatal("AliEMCALDigit::Operator+ -->" , mess.Data()) ;
+
}
}
}
Bool_t dejavu = kTRUE ;
Int_t old ;
for ( old = 0 ; (old < max2) && dejavu; old++) { //already have this primary?
- if(fIparent[old] == (digit.fIparent)[index])
+ if(fIparent[old] == digit.fIparent[index])
dejavu = kFALSE;
}
if(dejavu){
- fIparent[fNiparent] = (digit.fIparent)[index] ;
+ if(max2<fNMaxiparent){ fIparent[max2] = digit.fIparent[index] ;
fNiparent++ ;
- if(fNiparent>fNMaxiparent) {
- cout << "AliEMCALDigit >> Increase NMaxiparent "<< endl ;
- return *this ;
+ max2++;}
+ if(fNiparent==fNMaxiparent) {
+
+ TString mess = " NMaxiparent = " ;
+ mess += fNMaxiparent ;
+ mess += " is too small" ;
+ Fatal("AliEMCALDigit::Operator+ -->", mess.Data()) ;
+
}
}
}
return *this ;
}
+//____________________________________________________________________________
+AliEMCALDigit& AliEMCALDigit::operator*(Float_t factor)
+{
+ // Multiplies the amplitude by a factor
+
+ Float_t tempo = static_cast<Float_t>(fAmp) ;
+ tempo *= factor ;
+ fAmp = static_cast<Int_t>(TMath::Ceil(tempo)) ;
+ return *this ;
+}
+
//____________________________________________________________________________
ostream& operator << ( ostream& out , const AliEMCALDigit & digit)
{
// Prints the data of the digit
- out << "ID " << digit.fId << " Energy = " << digit.fAmp << endl ;
+ out << "ID " << digit.fId << " Energy = " << digit.fAmp << " Time = " << digit.fTime << endl ;
Int_t i,j ;
for(i=0;i<digit.fNprimary;i++)
out << "Primary " << i+1 << " = " << digit.fPrimary[i] << endl ;