X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOSRecPoint.cxx;h=d259034baa412e19acacd5c4731373bdd8f7fc99;hb=e8a76a2617f19ec07f7a5846e49d3f470164d638;hp=98229f480a12939db71e2f852bb2c7d75fe6bbf6;hpb=15605d3c0ec40003bf81418bd44b1a6eb4eae20e;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOSRecPoint.cxx b/PHOS/AliPHOSRecPoint.cxx index 98229f480a1..d259034baa4 100644 --- a/PHOS/AliPHOSRecPoint.cxx +++ b/PHOS/AliPHOSRecPoint.cxx @@ -12,63 +12,70 @@ * 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 "TGraph.h" +#include "TPaveText.h" +#include "TClonesArray.h" // --- Standard library --- -#include -#include // --- AliRoot header files --- #include "AliPHOSGeometry.h" #include "AliPHOSDigit.h" #include "AliPHOSRecPoint.h" - - +#include "AliPHOSGetter.h" ClassImp(AliPHOSRecPoint) //____________________________________________________________________________ AliPHOSRecPoint::AliPHOSRecPoint() - : AliRecPoint() + : AliRecPoint(), + fPHOSMod(0) { - fGeom = AliPHOSGeometry::GetInstance() ; - fPHOSMod = 0; + // ctor + + fMaxTrack = 0 ; } //____________________________________________________________________________ -AliPHOSRecPoint::~AliPHOSRecPoint() +AliPHOSRecPoint::AliPHOSRecPoint(const char * opt) + : AliRecPoint(opt), + fPHOSMod(0) { - // dtor + // ctor + + fMaxTrack = 200 ; } //____________________________________________________________________________ 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; } //___________________________________________________________________________ @@ -80,162 +87,193 @@ Int_t AliPHOSRecPoint::DistancetoPrimitive(Int_t px, Int_t py) } //______________________________________________________________________________ -void AliPHOSRecPoint::ExecuteEvent(Int_t event, Int_t px, Int_t py) +void AliPHOSRecPoint::ExecuteEvent(Int_t event, Int_t, Int_t) { // Execute action corresponding to one event // This member function is called when a AliPHOSRecPoint is clicked with the locator // // 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) { + static TGraph * digitgraph = 0 ; + static TPaveText* clustertext = 0 ; + + if (!gPad->IsEditable()) return; + + switch (event) { + + + case kButton1Down:{ + AliPHOSDigit * digit ; + +// Accessing geometry this way is equivalent to getting from gAlice +// to have Detector in Folder one have to load gAlice anyway +// AliPHOSLoader * gime = AliPHOSLoader::GetInstance(); +// AliPHOSGeometry * phosgeom = const_cast(gime->PHOSGeometry()); + AliPHOSGeometry * phosgeom = AliPHOSLoader::GetPHOSGeometry(); - 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) { + Int_t iDigit; + Int_t relid[4] ; - 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() ; + 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("dummy") ; + delete[] xi ; + delete[] zi ; + } + +break; + + case kButton1Up: + if (digitgraph) { + delete digitgraph ; + digitgraph = 0 ; + } + if (clustertext) { + delete clustertext ; + clustertext = 0 ; + } + + break; + } +} +//____________________________________________________________________________ +void AliPHOSRecPoint::EvalAll(TClonesArray * digits) +{ + //evaluates (if necessary) all RecPoint data members - break; - - case kButton1Up: - if (digitgraph) { - delete digitgraph ; - digitgraph = 0 ; - } - if (clustertext) { - delete clustertext ; - clustertext = 0 ; - } - - break; - - } + 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] ; - - AliPHOSDigit * digit ; - digit = (AliPHOSDigit *) fDigitsList[0] ; - AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ; + if( fPHOSMod == 0){ + Int_t relid[4] ; + + AliPHOSGeometry * phosgeom = (AliPHOSGetter::Instance())->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 ; - break ; - } - Int_t newprimary = newprimaryarray[jndex] ; - Int_t kndex ; - Bool_t already = kFALSE ; - for ( kndex = 0 ; kndex < counter ; kndex++ ) { //check if not already stored - if ( newprimary == tempo[kndex] ){ - already = kTRUE ; + if(nprimaries){ + 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 ( fMulTrack > fMaxTrack ) { + fMulTrack = - 1 ; + Error("EvalPrimaries", "GetNprimaries ERROR > increase fMaxTrack" ) ; break ; } - } // end of check - if ( !already) { // store it - tempo[counter] = newprimary ; - counter++ ; - } // store it - } // all primaries in digit - delete newprimaryarray ; + Int_t newprimary = newprimaryarray[jndex] ; + Int_t kndex ; + Bool_t already = kFALSE ; + 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[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] ; + if(fMulTrack) + fTracksList = new Int_t[fMulTrack] ; + for(index = 0; index < fMulTrack; index++) + fTracksList[index] = tempo[index] ; + + delete [] tempo ; - delete tempo ; +} +//____________________________________________________________________________ +void AliPHOSRecPoint::GetGlobalPosition(TVector3 & gpos, TMatrixF & gmat) const +{ + // returns the position of the cluster in the global reference system of ALICE + // and the uncertainty on this position + (AliPHOSGetter::Instance())->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,"") ; } +//______________________________________________________________________________ +