X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOSRecPoint.cxx;h=cef75692729028b7c9fdc37d43bd7f6030d8191f;hb=4696e514ad20b5b70890da786a9c32c92ff71742;hp=9fd62da06abccb9339f2ac3f3b0892840f830818;hpb=26d4b1416bff7602d6300ccf6c1e0fcfa2bd4f8e;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOSRecPoint.cxx b/PHOS/AliPHOSRecPoint.cxx index 9fd62da06ab..cef75692729 100644 --- a/PHOS/AliPHOSRecPoint.cxx +++ b/PHOS/AliPHOSRecPoint.cxx @@ -12,26 +12,28 @@ * about the suitability of this software for any purpose. It is * * provided "as is" without express or implied warranty. * **************************************************************************/ - +/* $Id$ */ //_________________________________________________________________________ -// PHOSRecPoint base class deriving from AliRecPoint -//*-- Author : Gines MARTINEZ SUBATECH -////////////////////////////////////////////////////////////////////////////// +// Base Class for PHOS Reconstructed Points +// Why should I put meaningless comments +// just to satisfy +// the code checker +//*-- Author: Gines Martinez (SUBATECH) // --- ROOT system --- #include "TPad.h" +#include "TClonesArray.h" // --- Standard library --- -#include -#include +#include +#include // --- AliRoot header files --- #include "AliPHOSGeometry.h" #include "AliPHOSDigit.h" #include "AliPHOSRecPoint.h" - - +#include "AliPHOSGetter.h" ClassImp(AliPHOSRecPoint) @@ -40,34 +42,41 @@ ClassImp(AliPHOSRecPoint) AliPHOSRecPoint::AliPHOSRecPoint() : AliRecPoint() { - fGeom = AliPHOSGeometry::GetInstance() ; + // ctor + + fMaxTrack = 200 ; fPHOSMod = 0; + } //____________________________________________________________________________ -AliPHOSRecPoint::~AliPHOSRecPoint() +AliPHOSRecPoint::AliPHOSRecPoint(const char * opt) : AliRecPoint(opt) { - // dtor + // ctor + + fMaxTrack = 200 ; + fPHOSMod = 0; + } + //____________________________________________________________________________ Int_t AliPHOSRecPoint::DistancetoPrimitive(Int_t px, Int_t py) { - //Compute distance from point px,py to a AliPHOSRecPoint considered as a Tmarker - // Compute the closest distance of approach from point px,py to this marker. - // The distance is computed in pixels units. - // - - TVector3 pos(0.,0.,0.) ; - GetLocalPosition( pos) ; - Float_t x = pos.X() ; - Float_t y = pos.Z() ; - const Int_t kMaxDiff = 10; - Int_t pxm = gPad->XtoAbsPixel(x); - Int_t pym = gPad->YtoAbsPixel(y); - Int_t dist = (px-pxm)*(px-pxm) + (py-pym)*(py-pym); - - if (dist > kMaxDiff) return 9999; - return dist; + // Compute distance from point px,py to a AliPHOSRecPoint considered as a Tmarker + // Compute the closest distance of approach from point px,py to this marker. + // The distance is computed in pixels units. + + TVector3 pos(0.,0.,0.) ; + GetLocalPosition( pos) ; + Float_t x = pos.X() ; + Float_t y = pos.Z() ; + const Int_t kMaxDiff = 10; + Int_t pxm = gPad->XtoAbsPixel(x); + Int_t pym = gPad->YtoAbsPixel(y); + Int_t dist = (px-pxm)*(px-pxm) + (py-pym)*(py-pym); + + if (dist > kMaxDiff) return 9999; + return dist; } //___________________________________________________________________________ @@ -86,155 +95,179 @@ void AliPHOSRecPoint::ExecuteEvent(Int_t event, Int_t px, Int_t py) // // If Left button is clicked on AliPHOSRecPoint, the digits are switched on // and switched off when the mouse button is released. - // // static Int_t pxold, pyold; - static TGraph * digitgraph = 0 ; - static TPaveText* clustertext = 0 ; - - if (!gPad->IsEditable()) return; - - switch (event) { - - - case kButton1Down:{ - AliPHOSDigit * digit ; - AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ; - Int_t iDigit; - Int_t relid[4] ; - Float_t xi[fMulDigit] ; - Float_t zi[fMulDigit] ; - - for(iDigit=0; iDigitAbsToRelNumbering(digit->GetId(), relid) ; - phosgeom->RelPosInModule(relid, xi[iDigit], zi[iDigit]) ; - } - - if (!digitgraph) { - digitgraph = new TGraph(fMulDigit,xi,zi); - digitgraph-> SetMarkerStyle(5) ; - digitgraph-> SetMarkerSize(1.) ; - digitgraph-> SetMarkerColor(1) ; - digitgraph-> Draw("P") ; - } - if (!clustertext) { + static TGraph * digitgraph = 0 ; + static TPaveText* clustertext = 0 ; - TVector3 pos(0.,0.,0.) ; - GetLocalPosition(pos) ; - clustertext = new TPaveText(pos.X()-10,pos.Z()+10,pos.X()+50,pos.Z()+35,"") ; - Text_t line1[40] ; - Text_t line2[40] ; - sprintf(line1,"Energy=%1.2f GeV",GetEnergy()) ; - sprintf(line2,"%d Digits",GetDigitsMultiplicity()) ; - clustertext ->AddText(line1) ; - clustertext ->AddText(line2) ; - clustertext ->Draw(""); - } - gPad->Update() ; - Print() ; - } - - break; - - case kButton1Up: - if (digitgraph) { - delete digitgraph ; - digitgraph = 0 ; - } - if (clustertext) { - delete clustertext ; - clustertext = 0 ; - } - - break; - + if (!gPad->IsEditable()) return; + + switch (event) { + + + case kButton1Down:{ + AliPHOSDigit * digit ; + AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ; + AliPHOSGeometry * phosgeom = const_cast(gime->PHOSGeometry()); + + Int_t iDigit; + Int_t relid[4] ; + + const Int_t kMulDigit=AliPHOSRecPoint::GetDigitsMultiplicity() ; + Float_t * xi = new Float_t [kMulDigit] ; + Float_t * zi = new Float_t [kMulDigit] ; + + for(iDigit = 0; iDigit < kMulDigit; iDigit++) { + Fatal("AliPHOSRecPoint::ExecuteEvent", "-> Something wrong with the code"); + digit = 0 ; //dynamic_cast((fDigitsList)[iDigit]); + phosgeom->AbsToRelNumbering(digit->GetId(), relid) ; + phosgeom->RelPosInModule(relid, xi[iDigit], zi[iDigit]) ; + } + + if (!digitgraph) { + digitgraph = new TGraph(fMulDigit,xi,zi); + digitgraph-> SetMarkerStyle(5) ; + digitgraph-> SetMarkerSize(1.) ; + digitgraph-> SetMarkerColor(1) ; + digitgraph-> Draw("P") ; + } + if (!clustertext) { + + TVector3 pos(0.,0.,0.) ; + GetLocalPosition(pos) ; + clustertext = new TPaveText(pos.X()-10,pos.Z()+10,pos.X()+50,pos.Z()+35,"") ; + Text_t line1[40] ; + Text_t line2[40] ; + sprintf(line1,"Energy=%1.2f GeV",GetEnergy()) ; + sprintf(line2,"%d Digits",GetDigitsMultiplicity()) ; + clustertext ->AddText(line1) ; + clustertext ->AddText(line2) ; + clustertext ->Draw(""); + } + gPad->Update() ; + Print() ; + delete[] xi ; + delete[] zi ; } + +break; + + case kButton1Up: + if (digitgraph) { + delete digitgraph ; + digitgraph = 0 ; + } + if (clustertext) { + delete clustertext ; + clustertext = 0 ; + } + + break; + + } } +//____________________________________________________________________________ +void AliPHOSRecPoint::EvalAll(Float_t logWeight,TClonesArray * digits) { + //evaluates (if necessary) all RecPoint data members + EvalPrimaries(digits) ; +} //____________________________________________________________________________ -Int_t AliPHOSRecPoint::GetPHOSMod() -{ - if(fPHOSMod > 0) - return fPHOSMod ; +void AliPHOSRecPoint::EvalPHOSMod(AliPHOSDigit * digit) +{ + // Returns the PHOS module in which the RecPoint is found - Int_t relid[4] ; + if( fPHOSMod == 0){ + Int_t relid[4] ; - AliPHOSDigit * digit ; - digit = (AliPHOSDigit *) fDigitsList[0] ; - AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ; + AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ; + AliPHOSGeometry * phosgeom = const_cast(gime->PHOSGeometry()); phosgeom->AbsToRelNumbering(digit->GetId(), relid) ; fPHOSMod = relid[0]; - return fPHOSMod ; + } } //______________________________________________________________________________ -void AliPHOSRecPoint::GetPrimaries(Int_t & number, Int_t * list) +void AliPHOSRecPoint::EvalPrimaries(TClonesArray * digits) { - AliPHOSDigit * digit ; - Int_t index ; - Int_t maxcounter = 3 ; - Int_t counter = 0 ; - Int_t * tempo = new Int_t[maxcounter] ; + // Constructs the list of primary particles (tracks) which have contributed to this RecPoint + AliPHOSDigit * digit ; + Int_t * tempo = new Int_t[fMaxTrack] ; + + Int_t index ; for ( index = 0 ; index < GetDigitsMultiplicity() ; index++ ) { // all digits - digit = (AliPHOSDigit *) fDigitsList[index] ; + digit = dynamic_cast(digits->At( fDigitsList[index] )) ; Int_t nprimaries = digit->GetNprimary() ; Int_t * newprimaryarray = new Int_t[nprimaries] ; Int_t ii ; for ( ii = 0 ; ii < nprimaries ; ii++) newprimaryarray[ii] = digit->GetPrimary(ii+1) ; + Int_t jndex ; for ( jndex = 0 ; jndex < nprimaries ; jndex++ ) { // all primaries in digit - if ( counter > maxcounter ) { - number = - 1 ; - cout << "AliPHOSRecPoint::GetNprimaries ERROR > increase maxcounter " << endl ; + if ( fMulTrack > fMaxTrack ) { + fMulTrack = - 1 ; + cout << "AliPHOSRecPoint::GetNprimaries ERROR > increase fMaxTrack " << endl ; break ; } Int_t newprimary = newprimaryarray[jndex] ; Int_t kndex ; Bool_t already = kFALSE ; - for ( kndex = 0 ; kndex < counter ; kndex++ ) { //check if not already stored + for ( kndex = 0 ; kndex < fMulTrack ; kndex++ ) { //check if not already stored if ( newprimary == tempo[kndex] ){ already = kTRUE ; break ; } } // end of check - if ( !already) { // store it - tempo[counter] = newprimary ; - counter++ ; + if ( !already) { // store it + tempo[fMulTrack] = newprimary ; + fMulTrack++ ; } // store it } // all primaries in digit delete newprimaryarray ; } // all digits - number = counter ; - for ( index = 0 ; index < number ; index ++ ) - list[index] = tempo[index] ; - delete tempo ; + fTracksList = new Int_t[fMulTrack] ; + for(index = 0; index < fMulTrack; index++) + fTracksList[index] = tempo[index] ; + + delete tempo ; + +} +//____________________________________________________________________________ +void AliPHOSRecPoint::GetGlobalPosition(TVector3 & gpos, TMatrix & gmat) const +{ + // returns the position of the cluster in the global reference system of ALICE + // and the uncertainty on this position + + + AliPHOSGetter::GetInstance()->PHOSGeometry()->GetGlobal(this, gpos, gmat) ; + } + //______________________________________________________________________________ void AliPHOSRecPoint::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 = 5 ; - - 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 = 5 ; + + if (!gPad->IsBatch()) { + gVirtualX->SetMarkerColor(markercolor) ; + gVirtualX->SetMarkerSize (markersize) ; + gVirtualX->SetMarkerStyle(markerstyle) ; + } + gPad->SetAttMarkerPS(markercolor,markerstyle,markersize) ; + gPad->PaintPolyMarker(1,&x,&y,"") ; }