* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
+/* $Id$ */
//_________________________________________________________________________
-// RecPoint in the PHOS PPSD: a list of AliPHOSDigit's
-//*-- Author : Yves Schutz SUBATECH
-//////////////////////////////////////////////////////////////////////////////
+// A RecPoint (cluster) in the PPSD
+// A PPSD RecPoint ends up to be a single digit
+// Why should I put meaningless comments
+// just to satisfy
+// the code checker
+//
+//*-- Yves Schutz (SUBATECH)
// --- ROOT system ---
+
#include "TPad.h"
-// --- Standard library ---
-#include <iostream>
-#include <cassert>
+// --- Standard library ---
// --- AliRoot header files ---
#include "AliPHOSGeometry.h"
#include "AliPHOSPpsdRecPoint.h"
+#include "AliPHOSCpvRecPoint.h"
#include "AliRun.h"
+#include "AliPHOSLoader.h"
ClassImp(AliPHOSPpsdRecPoint)
//____________________________________________________________________________
AliPHOSPpsdRecPoint::AliPHOSPpsdRecPoint(void)
{
+ // ctor
+
fMulDigit = 0 ;
- AliPHOSGeometry * geom = AliPHOSGeometry::GetInstance() ;
- fDelta = geom->GetCrystalSize(0) ;
fLocPos.SetX(1000000.) ; //Local position should be evaluated
}
//____________________________________________________________________________
-AliPHOSPpsdRecPoint::~AliPHOSPpsdRecPoint(void) // dtor
-{
- //dtor
-}
-
-//____________________________________________________________________________
-void AliPHOSPpsdRecPoint::AddDigit(AliDigitNew & digit, Float_t Energy)
+void AliPHOSPpsdRecPoint::AddDigit(AliPHOSDigit & digit, Float_t Energy)
{
// adds a digit to the digits list
// and accumulates the total amplitude and the multiplicity
delete [] tempo ;
}
- fDigitsList[fMulDigit++] = (int) &digit ;
+ 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::GetLocalPosition(TVector3 &LPos){
-
- if( fLocPos.X() < 1000000.) { //allready evaluated
- LPos = fLocPos ;
- return ;
- }
+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::EvalLocalPosition(Float_t logWeight,TClonesArray * digits )
+{
+ // Calculates the local position in the PHOS-PPSD-module coordinates
+
Int_t relid[4] ;
Float_t x = 0. ;
Float_t z = 0. ;
- AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ;
+ AliPHOSLoader * gime = AliPHOSLoader::GetInstance() ;
+ AliPHOSGeometry * phosgeom = (AliPHOSGeometry*)gime->PHOSGeometry();
AliPHOSDigit * digit ;
Int_t iDigit;
for(iDigit = 0; iDigit < fMulDigit; iDigit++) {
- digit = (AliPHOSDigit *) 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] ;
- AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ;
+ AliPHOSLoader * gime = AliPHOSLoader::GetInstance() ;
+ AliPHOSGeometry * phosgeom = (AliPHOSGeometry*)gime->PHOSGeometry();
+
- AliPHOSDigit *digit = (AliPHOSDigit *)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 *)
{
-//*-*-*-*-*-*-*-*-*-*-*Paint this ALiRecPoint as a TMarker with its current attributes*-*-*-*-*-*-*
-//*-* =============================================
- TVector3 pos(0.,0.,0.) ;
- GetLocalPosition(pos) ;
- Coord_t x = pos.X() ;
- Coord_t y = pos.Z() ;
- Color_t markercolor = 1 ;
- Size_t markersize = 1. ;
- Style_t markerstyle = 2 ;
- if (GetUp())
- markerstyle = 3 ;
-
- if (!gPad->IsBatch()) {
- gVirtualX->SetMarkerColor(markercolor);
- gVirtualX->SetMarkerSize (markersize);
- gVirtualX->SetMarkerStyle(markerstyle);
- }
- gPad->SetAttMarkerPS(markercolor,markerstyle,markersize);
- gPad->PaintPolyMarker(1,&x,&y,"");
-
-
+ //*-*-*-*-*-*-*-*-*-*-*Paint this ALiRecPoint as a TMarker with its current attributes*-*-*-*-*-*-*
+ //*-* =============================================
+
+ TVector3 pos(0.,0.,0.) ;
+ GetLocalPosition(pos) ;
+ Coord_t x = pos.X() ;
+ Coord_t y = pos.Z() ;
+ Color_t markercolor = 1 ;
+ Size_t markersize = 1. ;
+ Style_t markerstyle = 2 ;
+ if (GetUp())
+ markerstyle = 3 ;
+
+ if (!gPad->IsBatch()) {
+ gVirtualX->SetMarkerColor(markercolor);
+ gVirtualX->SetMarkerSize (markersize);
+ gVirtualX->SetMarkerStyle(markerstyle);
+ }
+ gPad->SetAttMarkerPS(markercolor,markerstyle,markersize);
+ gPad->PaintPolyMarker(1,&x,&y,"");
+
+
}
//____________________________________________________________________________
void AliPHOSPpsdRecPoint::Print(Option_t * option)
{
- cout << "AliPHOSPpsdRecPoint: " << endl ;
+ // Print the digits information
- AliPHOSDigit * digit ;
- Int_t iDigit;
- AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ;
-
- Float_t xi ;
- Float_t zi ;
- Int_t relid[4] ;
-
- for(iDigit=0; iDigit<fMulDigit; iDigit++) {
- digit = (AliPHOSDigit *) fDigitsList[iDigit];
- 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 ;
- }
- cout << " Multiplicity = " << fMulDigit << endl ;
-}
-
-
-// //____________________________________________________________________________
-// AliPHOSPpsdRecPoint& AliPHOSPpsdRecPoint::operator = (AliPHOSPpsdRecPoint Clu)
-// {
-// int* DL = Clu.GetDigitsList() ;
+ Info("Print","AliPHOSPpsdRecPoint: " ) ;
-// if(fDigitsList)
-// delete fDigitsList ;
-
-// AliPHOSDigit * digit ;
-
-// Int_t iDigit;
-
-// for(iDigit=0; iDigit<fMulDigit; iDigit++) {
-// digit = (AliPHOSDigit *) DL[iDigit];
-// AddDigit(*digit) ;
-// }
-
-// fDelta = Clu.GetDelta() ;
-// delete DL ;
+ Int_t iDigit;
+ Info("Print"," Digit{s} # ") ;
+ for(iDigit=0; iDigit<fMulDigit; iDigit++)
+ Info("Print", "%d", fDigitsList[iDigit]) ;
+ Info("Print"," Multiplicity = %d", fMulDigit ) ;
+ Info("Print"," Stored at position %d", fIndexInList ) ;
+}
-// fAmp = Clu.GetEnergy() ;
-// fGeom = Clu.GetGeom() ;
-// TVector3 LocPos;
-// Clu.GetLocalPosition(LocPos) ;
-// fLocPos = LocPos;
-// fMulDigit = Clu.GetMultiplicity() ;
-// fMaxDigit = Clu.GetMaximumMultiplicity() ;
-// fPHOSMod = Clu.GetPHOSMod() ;
-
-// return *this ;
-// }