X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOSDigit.cxx;h=193c61e21570461b5f9f158b3b97a571b91c19cd;hb=b671cde111ea27cfcc84c2768d76414018c75c38;hp=aa492a98cc83ff464319ee82eb0ecd2fa5683107;hpb=6a3f1304e9544d50138605fa6467b55773963293;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOSDigit.cxx b/PHOS/AliPHOSDigit.cxx index aa492a98cc8..193c61e2157 100644 --- a/PHOS/AliPHOSDigit.cxx +++ b/PHOS/AliPHOSDigit.cxx @@ -13,73 +13,151 @@ * provided "as is" without express or implied warranty. * **************************************************************************/ +/* $Id$ */ + +/* History of cvs commits: + * + * $Log$ + * Revision 1.39 2006/04/22 15:04:24 hristov + * Effective C++ initialization of data members in the default constructor + * + * Revision 1.38 2006/04/22 10:30:17 hristov + * Add fEnergy to AliPHOSDigit and operate with EMC amplitude in energy units (Yu.Kharlov) + * + * Revision 1.37 2005/05/28 14:19:04 schutz + * Compilation warnings fixed by T.P. + * + */ + //_________________________________________________________________________ -// Digit class for PHOS that contains an absolute ID and an energy -//*-- Author : Laurent Aphecetche SUBATECH -////////////////////////////////////////////////////////////////////////////// +// PHOS 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 PHOS crystal or PPSD gas cell +// +//*-- Author: Laurent Aphecetche & Yves Schutz (SUBATECH) & Dmitri Peressounko (RRC KI & SUBATECH) + // --- ROOT system --- -// --- Standard library --- +#include "TMath.h" -#include -#include +// --- Standard library --- // --- AliRoot header files --- #include "AliPHOSDigit.h" + + ClassImp(AliPHOSDigit) //____________________________________________________________________________ - AliPHOSDigit::AliPHOSDigit() : fPrimary(0) +AliPHOSDigit::AliPHOSDigit() : + AliDigitNew(), + fNprimary(0), + fPrimary(0x0), + fEnergy(0.), + fTime(0.), + fTimeR(0.) + { - fNprimary = 0 ; - Int_t index ; - for (index = 1 ; index <= 3 ; index++) - SetPrimary(index, -1) ; + // default ctor } //____________________________________________________________________________ -AliPHOSDigit::AliPHOSDigit(Int_t primary, Int_t id, Int_t DigEnergy) +AliPHOSDigit::AliPHOSDigit(Int_t primary, Int_t id, Int_t digEnergy, Float_t time, Int_t index) : + fNprimary(0), + fPrimary(0), + fEnergy(0.f), + fTime(0.f), + fTimeR(0.f) { - fId = id ; - fAmp = DigEnergy ; - fPrimary = new Int_t[1] ; - fPrimary[0] = primary ; - fNprimary = 1 ; - SetPrimary(1, primary) ; - + // ctor with all data + + fAmp = digEnergy ; + fEnergy = 0 ; + fTime = time ; + fTimeR = fTime ; + fId = id ; + fIndexInList = index ; + if( primary != -1){ + fNprimary = 1 ; + fPrimary = new Int_t[fNprimary] ; + fPrimary[0] = primary ; + } + else{ //If the contribution of this primary smaller than fDigitThreshold (AliPHOSv1) + fNprimary = 0 ; + fPrimary = 0 ; + } } //____________________________________________________________________________ -AliPHOSDigit::AliPHOSDigit(const AliPHOSDigit & digit) +AliPHOSDigit::AliPHOSDigit(Int_t primary, Int_t id, Float_t energy, Float_t time, Int_t index) : + fNprimary(0), + fPrimary(0), + fEnergy(0.f), + fTime(0.f), + fTimeR(0.f) { - fId = digit.fId; - fAmp = digit.fAmp ; - fNprimary = digit.GetNprimary() ; - fPrimary = new Int_t[fNprimary] ; - Int_t * primary = digit.GetPrimary() ; - Int_t index ; - for ( index = 0 ; index < fNprimary ; index++ ) { - fPrimary[index] = primary[index] ; - SetPrimary( index+1, digit.GetPrimary(index+1) ) ; + // ctor with all data + + fAmp = 0 ; + fEnergy = energy ; + fTime = time ; + fTimeR = fTime ; + fId = id ; + fIndexInList = index ; + if( primary != -1){ + fNprimary = 1 ; + fPrimary = new Int_t[fNprimary] ; + fPrimary[0] = primary ; + } + else{ //If the contribution of this primary smaller than fDigitThreshold (AliPHOSv1) + fNprimary = 0 ; + fPrimary = 0 ; } +} +//____________________________________________________________________________ +AliPHOSDigit::AliPHOSDigit(const AliPHOSDigit & digit) : + AliDigitNew(digit), + fNprimary(digit.fNprimary), + fPrimary(0), + fEnergy(digit.fEnergy), + fTime(digit.fTime), + fTimeR(digit.fTimeR) + +{ + // copy ctor + if(fNprimary){ + fPrimary = new Int_t[fNprimary] ; + for (Int_t i = 0; i < fNprimary ; i++) + fPrimary[i] = digit.fPrimary[i] ; + } + else + fPrimary = 0 ; + fAmp = digit.fAmp ; + fId = digit.fId; + fIndexInList = digit.fIndexInList ; } //____________________________________________________________________________ -AliPHOSDigit::~AliPHOSDigit() +AliPHOSDigit::~AliPHOSDigit() { - if ( fPrimary != 0 ) - delete fPrimary ; + // Delete array of primiries if any + if(fPrimary) + delete [] fPrimary ; } //____________________________________________________________________________ -Int_t AliPHOSDigit::Compare(TObject * obj) +Int_t AliPHOSDigit::Compare(const TObject * obj) const { - Int_t rv ; + // Compares two digits with respect to its Id + // to sort according increasing Id + + Int_t rv ; AliPHOSDigit * digit = (AliPHOSDigit *)obj ; @@ -99,30 +177,37 @@ Int_t AliPHOSDigit::Compare(TObject * obj) //____________________________________________________________________________ Int_t AliPHOSDigit::GetPrimary(Int_t index) const { - Int_t rv = -1 ; - if ( index > 3 ) - cout << "AliPHOSDigit ERROR > only 3 primaries allowed" << endl ; - else { - switch (index) { - case 1 : - rv = fPrimary1 ; - break ; - case 2 : - rv = fPrimary2 ; - break ; - case 3 : - rv = fPrimary3 ; - break ; - } + // retrieves the primary particle number given its index in the list + Int_t rv = -1 ; + if ( index <= fNprimary && index > 0){ + rv = fPrimary[index-1] ; } return rv ; - + +} +//____________________________________________________________________________ +void AliPHOSDigit::Print(const Option_t *) const +{ + // Print the digit together with list of primaries + printf("PHOS digit: Amp=%d/E=%.3f, Id=%d, Time=%.3e, TimeR=%.3e, NPrim=%d ",fAmp,fEnergy,fId,fTime,fTimeR,fNprimary); + for(Int_t index = 0; index 0){ + Int_t *tmp = new Int_t[fNprimary+digit.fNprimary] ; + if(fAmp < digit.fAmp || fEnergy < digit.fEnergy){//most energetic primary in second digit => first primaries in list from second digit + for (Int_t index = 0 ; index < digit.fNprimary ; index++) + tmp[index]=(digit.fPrimary)[index] ; + for (Int_t index = 0 ; index < fNprimary ; index++) + tmp[index+digit.fNprimary]=fPrimary[index] ; + } + else{ //add new primaries to the end + for (Int_t index = 0 ; index < fNprimary ; index++) + tmp[index]=fPrimary[index] ; + for (Int_t index = 0 ; index < digit.fNprimary ; index++) + tmp[index+fNprimary]=(digit.fPrimary)[index] ; + } + if(fPrimary) + delete []fPrimary ; + fPrimary = tmp ; + } + fNprimary+=digit.fNprimary ; + fAmp += digit.fAmp ; + fEnergy += digit.fEnergy ; + if(fTime > digit.fTime) + fTime = digit.fTime ; + fTimeR = fTime ; + return *this ; +} +//____________________________________________________________________________ +AliPHOSDigit& AliPHOSDigit::operator *= (Float_t factor) +{ + // Multiplies the amplitude by a factor - for ( index = 0 ; index < oldfNprimary ; index++ ) { - fPrimary[index] = tempo[index] ; - } - - Int_t jndex = 0 ; - for ( index = oldfNprimary ; index < fNprimary ; index++ ) { - fPrimary[index] = digit.fPrimary[jndex] ; - jndex++ ; - } - - // Here comes something crummy ... but I do not know how to stream pointers - // because AliPHOSDigit is in a TCLonesArray - - if ( fNprimary > 3 ) { - cout << "AliPHOSDigit + operator ERROR > too many primaries, modify AliPHOSDigit" << endl ; - fNprimary = 3 ; - } - else { - switch (fNprimary) { - case 1 : - SetPrimary(1, fPrimary[0]) ; - break ; - case 2 : - SetPrimary(2, fPrimary[1]) ; - break ; - case 3: - SetPrimary(3, fPrimary[2]) ; - break ; - } - } - - // end of crummy stuff - + Float_t tempo = static_cast(fAmp) ; + tempo *= factor ; + fAmp = static_cast(TMath::Ceil(tempo)) ; return *this ; } //____________________________________________________________________________ ostream& operator << ( ostream& out , const AliPHOSDigit & digit) { - out << "ID " << digit.fId << " Energy = " << digit.fAmp ; - + // Prints the data of the digit + +// out << "ID " << digit.fId << " Energy = " << digit.fAmp << " Time = " << digit.fTime << endl ; +// Int_t i ; +// for(i=0;i