A lot of changes here:
[u/mrichter/AliRoot.git] / 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 #include "AliPHOSGetter.h"
40
41 ClassImp(AliPHOSPpsdRecPoint)
42
43 //____________________________________________________________________________
44 AliPHOSPpsdRecPoint::AliPHOSPpsdRecPoint(void)
45
46   // ctor
47
48   fMulDigit = 0 ;  
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   AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ; 
191   AliPHOSGeometry * phosgeom =  (AliPHOSGeometry*)gime->PHOSGeometry();
192   
193   AliPHOSDigit * digit ;
194   Int_t iDigit;
195
196   for(iDigit = 0; iDigit < fMulDigit; iDigit++) {
197     digit = (AliPHOSDigit *) digits->At(fDigitsList[iDigit]) ; 
198  
199     Float_t xi ;
200     Float_t zi ;
201     phosgeom->AbsToRelNumbering(digit->GetId(), relid) ;
202     phosgeom->RelPosInModule(relid, xi, zi);
203     x += xi ;
204     z += zi ;
205   }
206
207   x /= fMulDigit ;
208   z /= fMulDigit ;
209
210   fLocPos.SetX(x)  ;
211   fLocPos.SetY(0.) ;
212   fLocPos.SetZ(z)  ;
213
214 }
215
216 //____________________________________________________________________________
217 void AliPHOSPpsdRecPoint::EvalUp(TClonesArray * digits)
218 {
219   // Are we in the uper PPSD module ?
220
221   Int_t relid[4] ;
222   
223   AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ; 
224   AliPHOSGeometry * phosgeom =  (AliPHOSGeometry*)gime->PHOSGeometry();
225   
226   
227   AliPHOSDigit *digit = (AliPHOSDigit *) digits->At(fDigitsList[0]) ; 
228   
229   phosgeom->AbsToRelNumbering(digit->GetId(),relid);
230
231   if((Int_t)TMath::Ceil((Float_t)relid[1]/
232                         (phosgeom->GetNumberOfModulesPhi()*phosgeom->GetNumberOfModulesZ())-0.0001 ) > 1) 
233     fUp = kFALSE ;
234   else  
235     fUp = kTRUE ;
236   
237 }
238 //______________________________________________________________________________
239 void AliPHOSPpsdRecPoint::Paint(Option_t *)
240 {
241   //*-*-*-*-*-*-*-*-*-*-*Paint this ALiRecPoint as a TMarker  with its current attributes*-*-*-*-*-*-*
242   //*-*                  =============================================
243
244   TVector3 pos(0.,0.,0.) ;
245   GetLocalPosition(pos) ;
246   Coord_t x = pos.X() ;
247   Coord_t y = pos.Z() ;
248   Color_t markercolor = 1 ;
249   Size_t  markersize  = 1. ;
250   Style_t markerstyle = 2 ;
251   if (GetUp()) 
252     markerstyle = 3 ;
253   
254   if (!gPad->IsBatch()) {
255     gVirtualX->SetMarkerColor(markercolor);
256     gVirtualX->SetMarkerSize (markersize);
257     gVirtualX->SetMarkerStyle(markerstyle);
258   }
259   gPad->SetAttMarkerPS(markercolor,markerstyle,markersize);
260   gPad->PaintPolyMarker(1,&x,&y,"");
261   
262   
263 }
264
265 //____________________________________________________________________________
266 void AliPHOSPpsdRecPoint::Print(Option_t * option) 
267 {
268   // Print the digits information 
269   
270   cout << "AliPHOSPpsdRecPoint: " << endl ;
271   
272   Int_t iDigit; 
273   cout << " Digit{s} # " ; 
274   for(iDigit=0; iDigit<fMulDigit; iDigit++) 
275     cout  << fDigitsList[iDigit] << "  " ;  
276   cout << endl   ;  
277
278   cout << "       Multiplicity    = " << fMulDigit  << endl ;
279   cout << "       Stored at position " << fIndexInList << endl ; 
280 }
281
282