PHOS/AliPHOSPpsdRecPoint.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15 /* $Id$ */
  16 //_________________________________________________________________________
  17 //  A RecPoint (cluster) in the PPSD
  18 //  A PPSD RecPoint ends up to be a single digit
  19 //  Why should I put meaningless comments
  20 //  just to satisfy
  21 //  the code checker
  22 //
  23 //*--  Yves Schutz (SUBATECH)
  24
  25 // --- ROOT system ---
  26
  27 #include "TPad.h"
  28
  29 // --- Standard library ---
  30
  31 #include <iostream.h>
  32
  33 // --- AliRoot header files ---
  34
  35 #include "AliPHOSGeometry.h"
  36 #include "AliPHOSPpsdRecPoint.h"
  37 #include "AliPHOSCpvRecPoint.h"
  38 #include "AliRun.h"
  39
  40 ClassImp(AliPHOSPpsdRecPoint)
  41
  42 //____________________________________________________________________________
  43 AliPHOSPpsdRecPoint::AliPHOSPpsdRecPoint(void)
  44 {
  45   // ctor
  46
  47   fMulDigit = 0 ;
  48   fGeom = AliPHOSGeometry::GetInstance() ;
  49   fLocPos.SetX(1000000.)  ;      //Local position should be evaluated
  50 }
  51
  52 //____________________________________________________________________________
  53 void AliPHOSPpsdRecPoint::AddDigit(AliPHOSDigit & digit, Float_t Energy)
  54 {
  55   // adds a digit to the digits list
  56   // and accumulates the total amplitude and the multiplicity
  57
  58
  59   if ( fMulDigit >= fMaxDigit ) { // increase the size of the lists
  60     fMaxDigit*=2 ;
  61     int * tempo = new ( int[fMaxDigit] ) ;
  62     Int_t index ;
  63
  64     for ( index = 0 ; index < fMulDigit ; index++ )
  65       tempo[index] = fDigitsList[index] ;
  66
  67     delete []  fDigitsList ;
  68     fDigitsList =  new ( int[fMaxDigit] ) ;
  69
  70     for ( index = 0 ; index < fMulDigit ; index++ )
  71       fDigitsList[index] = tempo[index] ;
  72
  73     delete [] tempo ;
  74   }
  75
  76   fDigitsList[fMulDigit++]  =  digit.GetIndexInList() ;
  77   fAmp += Energy ;
  78   EvalPHOSMod(&digit) ;
  79 }
  80
  81 //____________________________________________________________________________
  82 Int_t AliPHOSPpsdRecPoint::Compare(const TObject * obj) const
  83 {
  84   // Compares according to the position
  85   Float_t delta = 1 ; //width of the "Sorting row"
  86
  87   Int_t rv ;
  88
  89   if( (strcmp(obj->ClassName() , "AliPHOSPpsdRecPoint" )) == 0)  // PPSD Rec Point
  90     {
  91      AliPHOSPpsdRecPoint * clu = (AliPHOSPpsdRecPoint *)obj ;
  92
  93      Float_t x1 , z1 ;    //This rec point
  94      Float_t x2 , z2 ;    //
  95
  96      Int_t phosmod1 ;
  97      Int_t phosmod2 ;
  98
  99      Int_t up1 ;
 100      Int_t up2 ;
 101
 102      if(GetUp()) // upper layer
 103        up1 = 0 ;
 104      else        // lower layer
 105        up1 = 1 ;
 106
 107      if(clu->GetUp()) // upper layer
 108        up2 = 0 ;
 109      else            // lower layer
 110        up2 = 1 ;
 111
 112      TVector3 posloc ;
 113      GetLocalPosition(posloc) ;
 114      x1 = posloc.X() ;
 115      z1 = posloc.Z() ;
 116      phosmod1 = GetPHOSMod();
 117      clu->GetLocalPosition(posloc) ;
 118      x2 = posloc.X() ;
 119      z2 = posloc.Z() ;
 120      phosmod2 = clu->GetPHOSMod();
 121
 122      if(phosmod1 == phosmod2 ) {
 123
 124        if(up1 == up2 ){
 125          Int_t rowdif = (Int_t)TMath::Ceil(x1/delta) - (Int_t) TMath::Ceil(x2/delta) ;
 126
 127          if (rowdif> 0)
 128            rv = 1 ;
 129          else if(rowdif < 0)
 130            rv = -1 ;
 131          else if(z1>z2)
 132            rv = -1 ;
 133          else
 134            rv = 1 ;
 135        }
 136
 137        else {
 138
 139          if(up1 < up2 ) // Upper level first (up = True or False, True > False)
 140            rv = -1 ;
 141          else
 142            rv = 1 ;
 143        }
 144
 145      } // if phosmod1 == phosmod2
 146
 147      else {
 148
 149        if(phosmod1 < phosmod2 )
 150          rv = -1 ;
 151        else
 152          rv = 1 ;
 153
 154      }
 155
 156      return rv ;
 157     }
 158   else
 159     {
 160       AliPHOSCpvRecPoint * clu  = (AliPHOSCpvRecPoint *) obj ;
 161       if(GetPHOSMod()  < clu->GetPHOSMod() )
 162         rv = -1 ;
 163       else
 164         rv = 1 ;
 165       return rv ;
 166     }
 167
 168
 169 }
 170
 171 //____________________________________________________________________________
 172 void AliPHOSPpsdRecPoint::EvalAll(Float_t logWeight,TClonesArray * digits )
 173 {
 174   // calculates all the characteristics of the RecPoint
 175   AliPHOSRecPoint::EvalAll(logWeight,digits) ;
 176   EvalLocalPosition(logWeight,digits) ;
 177   EvalUp(digits) ;
 178 }
 179
 180 //____________________________________________________________________________
 181 void AliPHOSPpsdRecPoint::EvalLocalPosition(Float_t logWeight,TClonesArray * digits )
 182 {
 183   // Calculates the local position in the PHOS-PPSD-module coordinates
 184
 185   Int_t relid[4] ;
 186
 187   Float_t x = 0. ;
 188   Float_t z = 0. ;
 189
 190   AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ;
 191
 192   AliPHOSDigit * digit ;
 193   Int_t iDigit;
 194
 195   for(iDigit = 0; iDigit < fMulDigit; iDigit++) {
 196     digit = (AliPHOSDigit *) digits->At(fDigitsList[iDigit]) ;
 197
 198     Float_t xi ;
 199     Float_t zi ;
 200     phosgeom->AbsToRelNumbering(digit->GetId(), relid) ;
 201     phosgeom->RelPosInModule(relid, xi, zi);
 202     x += xi ;
 203     z += zi ;
 204   }
 205
 206   x /= fMulDigit ;
 207   z /= fMulDigit ;
 208
 209   fLocPos.SetX(x)  ;
 210   fLocPos.SetY(0.) ;
 211   fLocPos.SetZ(z)  ;
 212
 213 }
 214
 215 //____________________________________________________________________________
 216 void AliPHOSPpsdRecPoint::EvalUp(TClonesArray * digits)
 217 {
 218   // Are we in the uper PPSD module ?
 219
 220   Int_t relid[4] ;
 221
 222   AliPHOSGeometry * phosgeom = (AliPHOSGeometry *) fGeom ;
 223
 224
 225   AliPHOSDigit *digit = (AliPHOSDigit *) digits->At(fDigitsList[0]) ;
 226
 227   phosgeom->AbsToRelNumbering(digit->GetId(),relid);
 228
 229   if((Int_t)TMath::Ceil((Float_t)relid[1]/
 230                         (phosgeom->GetNumberOfModulesPhi()*phosgeom->GetNumberOfModulesZ())-0.0001 ) > 1)
 231     fUp = kFALSE ;
 232   else
 233     fUp = kTRUE ;
 234
 235 }
 236 //______________________________________________________________________________
 237 void AliPHOSPpsdRecPoint::Paint(Option_t *)
 238 {
 239   //*-*-*-*-*-*-*-*-*-*-*Paint this ALiRecPoint as a TMarker  with its current attributes*-*-*-*-*-*-*
 240   //*-*                  =============================================
 241
 242   TVector3 pos(0.,0.,0.) ;
 243   GetLocalPosition(pos) ;
 244   Coord_t x = pos.X() ;
 245   Coord_t y = pos.Z() ;
 246   Color_t markercolor = 1 ;
 247   Size_t  markersize  = 1. ;
 248   Style_t markerstyle = 2 ;
 249   if (GetUp())
 250     markerstyle = 3 ;
 251
 252   if (!gPad->IsBatch()) {
 253     gVirtualX->SetMarkerColor(markercolor);
 254     gVirtualX->SetMarkerSize (markersize);
 255     gVirtualX->SetMarkerStyle(markerstyle);
 256   }
 257   gPad->SetAttMarkerPS(markercolor,markerstyle,markersize);
 258   gPad->PaintPolyMarker(1,&x,&y,"");
 259
 260
 261 }
 262
 263 //____________________________________________________________________________
 264 void AliPHOSPpsdRecPoint::Print(Option_t * option)
 265 {
 266   // Print the digits information
 267
 268   cout << "AliPHOSPpsdRecPoint: " << endl ;
 269
 270   Int_t iDigit;
 271   cout << " Digit{s} # " ;
 272   for(iDigit=0; iDigit<fMulDigit; iDigit++)
 273     cout  << fDigitsList[iDigit] << "  " ;
 274   cout << endl   ;
 275
 276   cout << "       Multiplicity    = " << fMulDigit  << endl ;
 277   cout << "       Stored at position " << fIndexInList << endl ;
 278 }
 279
 280