a quick update to correct a few coding conventions but RS2 and GC2 have been ignored
[u/mrichter/AliRoot.git] / PHOS / AliPHOSPIDv1.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
16 /* $Id$ */
17
18 //_________________________________________________________________________
19 // Implementation version v1 of the PHOS particle identifier 
20 // Identification is based on information from PPSD and EMC
21 //  Why should I put meaningless comments
22 //  just to satisfy
23 //  the code checker                
24 //                  
25 //*-- Author: Yves Schutz (SUBATECH)  & Gines Martinez (SUBATECH)
26
27
28 // --- ROOT system ---
29
30 // --- Standard library ---
31
32 #include <iostream.h>
33
34 // --- AliRoot header files ---
35
36 #include "AliPHOSPIDv1.h"
37 #include "AliPHOSTrackSegment.h"
38 #include "AliPHOSRecParticle.h"
39 #include "AliPHOSIndexToObject.h"
40
41 ClassImp( AliPHOSPIDv1) 
42
43 //____________________________________________________________________________
44 AliPHOSPIDv1::AliPHOSPIDv1():AliPHOSPID()
45
46   fCutOnDispersion = 2.0; 
47   fCutOnRelativeDistance = 3.0 ;
48 }
49
50 //____________________________________________________________________________
51 Float_t  AliPHOSPIDv1::GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * emcclu,AliPHOSPpsdRecPoint * PpsdClu, Bool_t &toofar, Option_t *  Axis)
52 {
53   // Calculates the distance between the EMC RecPoint and the PPSD RecPoint
54  
55   Float_t r = 1.e+10;
56   TVector3 vecEmc ;
57   TVector3 vecPpsd ;
58   
59   emcclu->GetLocalPosition(vecEmc) ;
60   PpsdClu->GetLocalPosition(vecPpsd)  ; 
61   if(emcclu->GetPHOSMod() == PpsdClu->GetPHOSMod())
62     { 
63
64       // Correct to difference in CPV and EMC position due to different distance to center.
65       // we assume, that particle moves from center
66       Float_t dCPV = fGeom->GetIPtoOuterCoverDistance();
67       Float_t dEMC = fGeom->GetIPtoCrystalSurface() ;
68       dEMC         = dEMC / dCPV ;
69       vecPpsd = dEMC * vecPpsd  - vecEmc ; 
70       r = vecPpsd.Mag() ;
71       if (Axis == "X") r = vecPpsd.X();
72       if (Axis == "Y") r = vecPpsd.Y();
73       if (Axis == "Z") r = vecPpsd.Z();
74       if (Axis == "R") r = vecPpsd.Mag();
75       
76     } 
77   else 
78     {
79       toofar = kTRUE ;
80     }
81   return r ;
82 }
83
84 //____________________________________________________________________________
85 void  AliPHOSPIDv1::MakeParticles(AliPHOSTrackSegment::TrackSegmentsList * trsl, 
86                                   AliPHOSRecParticle::RecParticlesList * rpl)
87 {
88   // Makes a RecParticle out of a TrackSegment
89
90   TIter next(trsl) ; 
91   AliPHOSTrackSegment * tracksegment ; 
92   Int_t index = 0 ; 
93   AliPHOSRecParticle * rp ; 
94   Bool_t tDistance;
95   Int_t type ; 
96   Int_t showerprofile;  // 0 narrow and 1 wide
97   Int_t cpvdetector ;   // 1 hit and 0 no hit
98   Int_t pcdetector ;    // 1 hit and 0 no hit
99
100   while ( (tracksegment = (AliPHOSTrackSegment *)next()) ) {
101     new( (*rpl)[index] ) AliPHOSRecParticle(tracksegment) ;
102     rp = (AliPHOSRecParticle *)rpl->At(index) ; 
103     AliPHOSEmcRecPoint * recp = tracksegment->GetEmcRecPoint() ;
104     AliPHOSPpsdRecPoint * rpcpv = tracksegment->GetPpsdUpRecPoint() ;
105     AliPHOSPpsdRecPoint * rppc  = tracksegment->GetPpsdLowRecPoint() ;
106     
107     //     Float_t * lambda = new Float_t[2]; 
108     //     recp->GetElipsAxis(lambda) ; 
109     
110     //     // Looking at the lateral development of the shower
111     //     if ( ( lambda[0] > fLambda1m && lambda[0] < fLambda1M ) && // shower profile cut
112     //   ( lambda[1] > fLambda2m && lambda[1] < fLambda2M ) )         
113     //       //    Float_t R ;
114     //       //R=(lambda[0]-1.386)*(lambda[0]-1.386)+1.707*1.707*(lambda[1]-1.008)*(lambda[1]-1.008) ;
115     //       //if(R<0.35*0.35)
116
117     Float_t dispersion;
118     dispersion = recp->GetDispersion();
119     if (dispersion < fCutOnDispersion)
120       showerprofile = 0 ;   // NARROW PROFILE   
121     else      
122       showerprofile = 1 ;// WIDE PROFILE
123     
124     
125     // Looking at the photon conversion detector
126     if( tracksegment->GetPpsdLowRecPoint() == 0 )   
127       pcdetector = 0 ;  // No hit
128     else{      
129       if (GetDistanceInPHOSPlane(recp, rppc, tDistance, "R")  < fCutOnRelativeDistance) 
130         pcdetector = 1 ;  // hit
131       else
132         pcdetector = 0 ;
133     }
134     
135     // Looking at the photon conversion detector
136     if( tracksegment->GetPpsdUpRecPoint() == 0 )
137       cpvdetector = 0 ;  // No hit
138     else{  
139       if (GetDistanceInPHOSPlane(recp, rpcpv, tDistance, "R")< fCutOnRelativeDistance) 
140         cpvdetector = 1 ;  // Hit
141       else
142         cpvdetector = 0 ;
143     }
144      
145     type = showerprofile + 2 * pcdetector + 4 * cpvdetector ;
146     rp->SetType(type) ; 
147     index++ ; 
148   }
149     
150 }
151
152 //____________________________________________________________________________
153 void  AliPHOSPIDv1:: Print(const char * opt) 
154 {
155   // Print the parameters used for the particle type identification
156   
157   cout << "AliPHOSPIDv1 : cuts for the particle idendification based on the shower profile " << endl 
158        << fLambda1m << " < value1 < " << fLambda1M << endl 
159        << fLambda2m << " < value2 < " << fLambda2M << endl ;  
160
161 }
162
163 //____________________________________________________________________________
164 void  AliPHOSPIDv1::SetShowerProfileCuts(Float_t l1m, Float_t l1M, Float_t l2m, Float_t l2M)
165 {
166   // Modifies the parameters used for the particle type identification
167
168   fLambda1m = l1m ; 
169   fLambda1M = l1M ; 
170   fLambda2m = l2m ; 
171   fLambda2M = l2M ; 
172 }
173
174 //____________________________________________________________________________
175 void  AliPHOSPIDv1::SetRelativeDistanceCut(Float_t CutOnRelativeDistance)
176 {
177   // Modifies the parameters used for the particle type identification
178
179   fCutOnRelativeDistance = CutOnRelativeDistance ; 
180  
181 }
182