/************************************************************************** * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * * * * Author: The ALICE Off-line Project. * * Contributors are mentioned in the code where appropriate. * * * * Permission to use, copy, modify and distribute this software and its * * documentation strictly for non-commercial purposes is hereby granted * * without fee, provided that the above copyright notice appears in all * * copies and that both the copyright notice and this permission notice * * appear in the supporting documentation. The authors make no claims * * about the suitability of this software for any purpose. It is * * provided "as is" without express or implied warranty. * **************************************************************************/ /* $Id$ */ //_________________________________________________________________________ // Base Class for EMCAL 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 // --- AliRoot header files --- #include "AliEMCALGeometry.h" #include "AliEMCALDigit.h" #include "AliEMCALRecPoint.h" #include "AliEMCALGetter.h" ClassImp(AliEMCALRecPoint) //____________________________________________________________________________ AliEMCALRecPoint::AliEMCALRecPoint() : AliRecPoint() { // ctor fMaxTrack = 0 ; fTheta = fPhi = 0. ; fEMCALArm = 0; } //____________________________________________________________________________ AliEMCALRecPoint::AliEMCALRecPoint(const char * opt) : AliRecPoint(opt) { // ctor fMaxTrack = 200 ; fTheta = fPhi = 0. ; fEMCALArm = 1; } //____________________________________________________________________________ Int_t AliEMCALRecPoint::DistancetoPrimitive(Int_t px, Int_t py) { // Compute distance from point px,py to a AliEMCALRecPoint 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; } //___________________________________________________________________________ void AliEMCALRecPoint::Draw(Option_t *option) { // Draw this AliEMCALRecPoint with its current attributes AppendPad(option); } //______________________________________________________________________________ void AliEMCALRecPoint::ExecuteEvent(Int_t event, Int_t px, Int_t py) { // Execute action corresponding to one event // This member function is called when a AliEMCALRecPoint is clicked with the locator // // If Left button is clicked on AliEMCALRecPoint, 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:{ AliEMCALDigit * digit ; AliEMCALGetter * gime = AliEMCALGetter::GetInstance() ; AliEMCALGeometry * emcalgeom = const_cast(gime->EMCALGeometry()); Int_t iDigit; Int_t relid[4] ; const Int_t kMulDigit=AliEMCALRecPoint::GetDigitsMultiplicity() ; Float_t * xi = new Float_t [kMulDigit] ; Float_t * zi = new Float_t [kMulDigit] ; for(iDigit = 0; iDigit < kMulDigit; iDigit++) { Fatal("AliEMCALRecPoint::ExecuteEvent", " -> Something wrong with the code"); digit = 0 ; //dynamic_cast((fDigitsList)[iDigit]); emcalgeom->AbsToRelNumbering(digit->GetId(), relid) ; emcalgeom->PosInAlice(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 AliEMCALRecPoint::EvalAll(Float_t logWeight,TClonesArray * digits) { //evaluates (if necessary) all RecPoint data members EvalPrimaries(digits) ; } //____________________________________________________________________________ void AliEMCALRecPoint::EvalEMCALArm(AliEMCALDigit * digit) { // Returns the EMCAL module in which the RecPoint is found if( fEMCALArm == 0){ Int_t relid[4] ; AliEMCALGetter * gime = AliEMCALGetter::GetInstance() ; AliEMCALGeometry * emcalgeom = const_cast(gime->EMCALGeometry()); emcalgeom->AbsToRelNumbering(digit->GetId(), relid) ; fEMCALArm = relid[0]; } } //______________________________________________________________________________ void AliEMCALRecPoint::EvalPrimaries(TClonesArray * digits) { // Constructs the list of primary particles (tracks) which have contributed to this RecPoint AliEMCALDigit * digit ; Int_t * tempo = new Int_t[fMaxTrack] ; Int_t index ; for ( index = 0 ; index < GetDigitsMultiplicity() ; index++ ) { // all digits 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 ( fMulTrack > fMaxTrack ) { fMulTrack = - 1 ; cout << "AliEMCALRecPoint::GetNprimaries ERROR > increase fMaxTrack " << endl ; break ; } 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 fTracksList = new Int_t[fMulTrack] ; for(index = 0; index < fMulTrack; index++) fTracksList[index] = tempo[index] ; delete tempo ; } //____________________________________________________________________________ void AliEMCALRecPoint::GetGlobalPosition(TVector3 & gpos) const { // returns the position of the cluster in the global reference system of ALICE // and the uncertainty on this position AliEMCALGeometry * emcalgeom = AliEMCALGetter::GetInstance()->EMCALGeometry(); gpos.SetX(fPhi) ; gpos.SetY(emcalgeom->GetIPDistance() + emcalgeom->GetAirGap()) ; gpos.SetZ(fTheta) ; } //______________________________________________________________________________ void AliEMCALRecPoint::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,"") ; }