#include "AliPHOSGeometry.h"
#include "AliPHOSPpsdRecPoint.h"
+#include "AliPHOSCpvRecPoint.h"
#include "AliRun.h"
-#include "AliPHOSIndexToObject.h"
+#include "AliPHOSGetter.h"
ClassImp(AliPHOSPpsdRecPoint)
// ctor
fMulDigit = 0 ;
- AliPHOSGeometry * geom = AliPHOSGeometry::GetInstance() ;
- fDelta = geom->GetCrystalSize(0) ;
fLocPos.SetX(1000000.) ; //Local position should be evaluated
}
fDigitsList[fMulDigit++] = digit.GetIndexInList() ;
fAmp += Energy ;
+ EvalPHOSMod(&digit) ;
}
//____________________________________________________________________________
-Int_t AliPHOSPpsdRecPoint::Compare(TObject * obj)
+Int_t AliPHOSPpsdRecPoint::Compare(const TObject * obj) const
{
// Compares according to the position
-
+ Float_t delta = 1 ; //width of the "Sorting row"
+
Int_t rv ;
- AliPHOSPpsdRecPoint * clu = (AliPHOSPpsdRecPoint *)obj ;
-
-
- Float_t x1 , z1 ;
- Float_t x2 , z2 ;
-
- Int_t phosmod1 ;
- Int_t phosmod2 ;
-
- Int_t up1 ;
- Int_t up2 ;
+ if( (strcmp(obj->ClassName() , "AliPHOSPpsdRecPoint" )) == 0) // PPSD Rec Point
+ {
+ AliPHOSPpsdRecPoint * clu = (AliPHOSPpsdRecPoint *)obj ;
+
+ Float_t x1 , z1 ; //This rec point
+ Float_t x2 , z2 ; //
+
+ Int_t phosmod1 ;
+ Int_t phosmod2 ;
+
+ Int_t up1 ;
+ Int_t up2 ;
- if(GetUp()) // upper layer
- up1 = 0 ;
- else // lower layer
- up1 = 1 ;
-
- if(clu->GetUp()) // upper layer
- up2 = 0 ;
- else // lower layer
- up2 = 1 ;
-
- TVector3 posloc ;
- this->GetLocalPosition(posloc) ;
- x1 = posloc.X() ;
- z1 = posloc.Z() ;
- phosmod1 = this->GetPHOSMod();
- clu->GetLocalPosition(posloc) ;
- x2 = posloc.X() ;
- z2 = posloc.Z() ;
- phosmod2 = clu->GetPHOSMod();
-
- if(phosmod1 == phosmod2 ) {
-
- if(up1 == up2 ){
- Int_t rowdif = (Int_t)TMath::Ceil(x1/fDelta) - (Int_t) TMath::Ceil(x2/fDelta) ;
-
- if (rowdif> 0)
- rv = -1 ;
- else if(rowdif < 0)
- rv = 1 ;
- else if(z1>z2)
- rv = -1 ;
- else
- rv = 1 ;
+ if(GetUp()) // upper layer
+ up1 = 0 ;
+ else // lower layer
+ up1 = 1 ;
+
+ if(clu->GetUp()) // upper layer
+ up2 = 0 ;
+ else // lower layer
+ up2 = 1 ;
+
+ TVector3 posloc ;
+ GetLocalPosition(posloc) ;
+ x1 = posloc.X() ;
+ z1 = posloc.Z() ;
+ phosmod1 = GetPHOSMod();
+ clu->GetLocalPosition(posloc) ;
+ x2 = posloc.X() ;
+ z2 = posloc.Z() ;
+ phosmod2 = clu->GetPHOSMod();
+
+ if(phosmod1 == phosmod2 ) {
+
+ if(up1 == up2 ){
+ Int_t rowdif = (Int_t)TMath::Ceil(x1/delta) - (Int_t) TMath::Ceil(x2/delta) ;
+
+ if (rowdif> 0)
+ rv = 1 ;
+ else if(rowdif < 0)
+ rv = -1 ;
+ else if(z1>z2)
+ rv = -1 ;
+ else
+ rv = 1 ;
+ }
+
+ else {
+
+ if(up1 < up2 ) // Upper level first (up = True or False, True > False)
+ rv = -1 ;
+ else
+ rv = 1 ;
+ }
+
+ } // if phosmod1 == phosmod2
+
+ else {
+
+ if(phosmod1 < phosmod2 )
+ rv = -1 ;
+ else
+ rv = 1 ;
+
+ }
+
+ return rv ;
}
-
- else {
-
- if(up1 < up2 ) // Upper level first (up = True or False, True > False)
- rv = 1 ;
+ else
+ {
+ AliPHOSCpvRecPoint * clu = (AliPHOSCpvRecPoint *) obj ;
+ if(GetPHOSMod() < clu->GetPHOSMod() )
+ rv = -1 ;
else
- rv = - 1 ;
+ rv = 1 ;
+ return rv ;
}
-
- } // if phosmod1 == phosmod2
- else {
-
- if(phosmod1 < phosmod2 )
- rv = -1 ;
- else
- rv = 1 ;
}
-
- return rv ;
+
+//____________________________________________________________________________
+void AliPHOSPpsdRecPoint::EvalAll(Float_t logWeight,TClonesArray * digits )
+{
+ // calculates all the characteristics of the RecPoint
+ AliPHOSRecPoint::EvalAll(logWeight,digits) ;
+ EvalLocalPosition(logWeight,digits) ;
+ EvalUp(digits) ;
}
//____________________________________________________________________________
-void AliPHOSPpsdRecPoint::GetLocalPosition(TVector3 &LPos)
+void AliPHOSPpsdRecPoint::EvalLocalPosition(Float_t logWeight,TClonesArray * digits )
{
- // Calculates the local position in the PHOS-PPSD-module corrdinates
+ // Calculates the local position in the PHOS-PPSD-module coordinates
- AliPHOSIndexToObject * please = AliPHOSIndexToObject::GetInstance() ;
-
- if( fLocPos.X() < 1000000.) { //allready evaluated
- LPos = fLocPos ;
- return ;
- }
-
Int_t relid[4] ;
Float_t x = 0. ;
Float_t z = 0. ;
- AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ;
+ AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ;
+ AliPHOSGeometry * phosgeom = (AliPHOSGeometry*)gime->PHOSGeometry();
AliPHOSDigit * digit ;
Int_t iDigit;
for(iDigit = 0; iDigit < fMulDigit; iDigit++) {
- digit = (AliPHOSDigit *) ( please->GimeDigit(fDigitsList[iDigit]) );
+ digit = (AliPHOSDigit *) digits->At(fDigitsList[iDigit]) ;
Float_t xi ;
Float_t zi ;
fLocPos.SetY(0.) ;
fLocPos.SetZ(z) ;
- LPos = fLocPos ;
-
}
//____________________________________________________________________________
-Bool_t AliPHOSPpsdRecPoint::GetUp()
+void AliPHOSPpsdRecPoint::EvalUp(TClonesArray * digits)
{
// Are we in the uper PPSD module ?
Int_t relid[4] ;
- AliPHOSIndexToObject * please = AliPHOSIndexToObject::GetInstance() ;
-
- AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ;
+ AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ;
+ AliPHOSGeometry * phosgeom = (AliPHOSGeometry*)gime->PHOSGeometry();
- AliPHOSDigit *digit = (AliPHOSDigit *) ( please->GimeDigit(fDigitsList[0]) ) ;
+ AliPHOSDigit *digit = (AliPHOSDigit *) digits->At(fDigitsList[0]) ;
phosgeom->AbsToRelNumbering(digit->GetId(),relid);
- Bool_t up ;
if((Int_t)TMath::Ceil((Float_t)relid[1]/
(phosgeom->GetNumberOfModulesPhi()*phosgeom->GetNumberOfModulesZ())-0.0001 ) > 1)
- up = kFALSE ;
+ fUp = kFALSE ;
else
- up = kTRUE ;
+ fUp = kTRUE ;
- return up ;
}
//______________________________________________________________________________
void AliPHOSPpsdRecPoint::Paint(Option_t *)
cout << "AliPHOSPpsdRecPoint: " << endl ;
- AliPHOSDigit * digit ;
- Int_t iDigit;
- AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ;
-
- Float_t xi ;
- Float_t zi ;
- Int_t relid[4] ;
-
- AliPHOSIndexToObject * please = AliPHOSIndexToObject::GetInstance() ;
- for(iDigit=0; iDigit<fMulDigit; iDigit++) {
- digit = please->GimeDigit( fDigitsList[iDigit] ) ;
- if (digit) {
- phosgeom->AbsToRelNumbering(digit->GetId(), relid) ;
- phosgeom->RelPosInModule(relid, xi, zi);
- cout << " Id = " << digit->GetId() ;
- cout << " Phos mod = " << relid[0] ;
- cout << " PPSD mod = " << relid[1] ;
- cout << " x = " << xi ;
- cout << " z = " << zi ;
- cout << " Energy = " << digit->GetAmp() << endl ;
- }
- }
+ Int_t iDigit;
+ cout << " Digit{s} # " ;
+ for(iDigit=0; iDigit<fMulDigit; iDigit++)
+ cout << fDigitsList[iDigit] << " " ;
+ cout << endl ;
+
cout << " Multiplicity = " << fMulDigit << endl ;
cout << " Stored at position " << fIndexInList << endl ;
}