Adding the array of the recosntruction parameters (Marian)
[u/mrichter/AliRoot.git] / PHOS / AliPHOSGeometry.cxx
CommitLineData
d15a28e7 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
b2a60966 16/* $Id$ */
17
d15a28e7 18//_________________________________________________________________________
b2a60966 19// Geometry class for PHOS : singleton
a3dfe79c 20// PHOS consists of the electromagnetic calorimeter (EMCA)
21// and a charged particle veto either in the Subatech's version (PPSD)
22// or in the IHEP's one (CPV).
23// The EMCA/PPSD/CPV modules are parametrized so that any configuration
24// can be easily implemented
25// The title is used to identify the version of CPV used.
b2a60966 26//
05d33a3d 27// -- Author: Yves Schutz (SUBATECH) & Dmitri Peressounko (RRC "KI" & SUBATECH)
d15a28e7 28
29// --- ROOT system ---
30
31#include "TVector3.h"
32#include "TRotation.h"
e957fea8 33#include "TParticle.h"
91daaf24 34#include <TGeoManager.h>
268f57b1 35#include <TGeoMatrix.h>
d15a28e7 36
37// --- Standard library ---
38
d15a28e7 39// --- AliRoot header files ---
351dd634 40#include "AliLog.h"
d15a28e7 41#include "AliPHOSGeometry.h"
468794ea 42#include "AliPHOSEMCAGeometry.h"
710f859a 43#include "AliPHOSRecPoint.h"
d15a28e7 44
925e6570 45ClassImp(AliPHOSGeometry)
d15a28e7 46
a4e98857 47// these initialisations are needed for a singleton
05d33a3d 48AliPHOSGeometry * AliPHOSGeometry::fgGeom = 0 ;
49Bool_t AliPHOSGeometry::fgInit = kFALSE ;
9ec91567 50
d15a28e7 51//____________________________________________________________________________
3663622c 52AliPHOSGeometry::AliPHOSGeometry() :
bcc9e89e 53 AliPHOSGeoUtils(),
3663622c 54 fAngle(0.f),
55 fPHOSAngle(0),
56 fIPtoUpperCPVsurface(0),
308fb942 57 fCrystalShift(0),
58 fCryCellShift(0),
bcc9e89e 59 fRotMatrixArray(0)
3663622c 60{
e957fea8 61 // default ctor
62 // must be kept public for root persistency purposes, but should never be called by the outside world
e957fea8 63 fgGeom = 0 ;
e957fea8 64}
65
66//____________________________________________________________________________
3663622c 67AliPHOSGeometry::AliPHOSGeometry(const AliPHOSGeometry & rhs)
bcc9e89e 68 : AliPHOSGeoUtils(rhs),
3663622c 69 fAngle(rhs.fAngle),
70 fPHOSAngle(0),
71 fIPtoUpperCPVsurface(rhs.fIPtoUpperCPVsurface),
308fb942 72 fCrystalShift(rhs.fCrystalShift),
73 fCryCellShift(rhs.fCryCellShift),
bcc9e89e 74 fRotMatrixArray(0)
3663622c 75{
76 Fatal("cpy ctor", "not implemented") ;
77}
78
91daaf24 79//____________________________________________________________________________
3663622c 80AliPHOSGeometry::AliPHOSGeometry(const Text_t* name, const Text_t* title)
bcc9e89e 81 : AliPHOSGeoUtils(name, title),
3663622c 82 fAngle(0.f),
83 fPHOSAngle(0),
84 fIPtoUpperCPVsurface(0),
308fb942 85 fCrystalShift(0),
86 fCryCellShift(0),
bcc9e89e 87 fRotMatrixArray(0)
3663622c 88{
89 // ctor only for internal usage (singleton)
90 Init() ;
9a2cdbdf 91 fgGeom = this;
3663622c 92}
93
3663622c 94//____________________________________________________________________________
d15a28e7 95AliPHOSGeometry::~AliPHOSGeometry(void)
96{
b2a60966 97 // dtor
98
52a36ffd 99 if (fRotMatrixArray) fRotMatrixArray->Delete() ;
100 if (fRotMatrixArray) delete fRotMatrixArray ;
fa0bc588 101 if (fPHOSAngle ) delete[] fPHOSAngle ;
52a36ffd 102}
52a36ffd 103
91daaf24 104//____________________________________________________________________________
52a36ffd 105void AliPHOSGeometry::Init(void)
106{
a4e98857 107 // Initializes the PHOS parameters :
108 // IHEP is the Protvino CPV (cathode pad chambers)
710f859a 109
110 fgInit = kTRUE ;
05d33a3d 111
85698486 112 fAngle = 20;
05d33a3d 113
710f859a 114
115 fPHOSAngle = new Float_t[fNModules] ;
116
117 Float_t * emcParams = fGeometryEMCA->GetEMCParams() ;
118
581e32d4 119 fPHOSParams[0] = TMath::Max((Double_t)fGeometryCPV->GetCPVBoxSize(0)/2.,
a83b6179 120 (Double_t)(emcParams[0] - (emcParams[1]-emcParams[0])*
121 fGeometryCPV->GetCPVBoxSize(1)/2/emcParams[3]));
710f859a 122 fPHOSParams[1] = emcParams[1] ;
581e32d4 123 fPHOSParams[2] = TMath::Max((Double_t)emcParams[2], (Double_t)fGeometryCPV->GetCPVBoxSize(2)/2.);
710f859a 124 fPHOSParams[3] = emcParams[3] + fGeometryCPV->GetCPVBoxSize(1)/2. ;
125
126 fIPtoUpperCPVsurface = fGeometryEMCA->GetIPtoOuterCoverDistance() - fGeometryCPV->GetCPVBoxSize(1) ;
8e20650f 127
128 //calculate offset to crystal surface
129 Float_t * inthermo = fGeometryEMCA->GetInnerThermoHalfSize() ;
130 Float_t * strip = fGeometryEMCA->GetStripHalfSize() ;
131 Float_t* splate = fGeometryEMCA->GetSupportPlateHalfSize();
132 Float_t * crystal = fGeometryEMCA->GetCrystalHalfSize() ;
133 Float_t * pin = fGeometryEMCA->GetAPDHalfSize() ;
134 Float_t * preamp = fGeometryEMCA->GetPreampHalfSize() ;
135 fCrystalShift=-inthermo[1]+strip[1]+splate[1]+crystal[1]-fGeometryEMCA->GetAirGapLed()/2.+pin[1]+preamp[1] ;
136 fCryCellShift=crystal[1]-(fGeometryEMCA->GetAirGapLed()-2*pin[1]-2*preamp[1])/2;
137
710f859a 138 Int_t index ;
139 for ( index = 0; index < fNModules; index++ )
52a36ffd 140 fPHOSAngle[index] = 0.0 ; // Module position angles are set in CreateGeometry()
710f859a 141
710f859a 142 fRotMatrixArray = new TObjArray(fNModules) ;
05d33a3d 143
85698486 144 // Geometry parameters are calculated
145
146 SetPHOSAngles();
147 Double_t const kRADDEG = 180.0 / TMath::Pi() ;
148 Float_t r = GetIPtoOuterCoverDistance() + fPHOSParams[3] - GetCPVBoxSize(1) ;
149 for (Int_t iModule=0; iModule<fNModules; iModule++) {
150 fModuleCenter[iModule][0] = r * TMath::Sin(fPHOSAngle[iModule] / kRADDEG );
151 fModuleCenter[iModule][1] =-r * TMath::Cos(fPHOSAngle[iModule] / kRADDEG );
152 fModuleCenter[iModule][2] = 0.;
153
154 fModuleAngle[iModule][0][0] = 90;
155 fModuleAngle[iModule][0][1] = fPHOSAngle[iModule];
156 fModuleAngle[iModule][1][0] = 0;
157 fModuleAngle[iModule][1][1] = 0;
158 fModuleAngle[iModule][2][0] = 90;
159 fModuleAngle[iModule][2][1] = 270 + fPHOSAngle[iModule];
05d33a3d 160 }
161
52a36ffd 162}
163
5ccb0008 164//____________________________________________________________________________
52a36ffd 165AliPHOSGeometry * AliPHOSGeometry::GetInstance()
166{
a4e98857 167 // Returns the pointer of the unique instance; singleton specific
168
809cd394 169 return static_cast<AliPHOSGeometry *>( fgGeom ) ;
52a36ffd 170}
171
172//____________________________________________________________________________
173AliPHOSGeometry * AliPHOSGeometry::GetInstance(const Text_t* name, const Text_t* title)
174{
175 // Returns the pointer of the unique instance
a4e98857 176 // Creates it with the specified options (name, title) if it does not exist yet
177
52a36ffd 178 AliPHOSGeometry * rv = 0 ;
179 if ( fgGeom == 0 ) {
180 if ( strcmp(name,"") == 0 )
181 rv = 0 ;
182 else {
183 fgGeom = new AliPHOSGeometry(name, title) ;
184 if ( fgInit )
185 rv = (AliPHOSGeometry * ) fgGeom ;
186 else {
187 rv = 0 ;
188 delete fgGeom ;
189 fgGeom = 0 ;
190 }
191 }
192 }
193 else {
21cd0c07 194 if ( strcmp(fgGeom->GetName(), name) != 0 )
351dd634 195 ::Error("GetInstance", "Current geometry is %s. You cannot call %s",
196 fgGeom->GetName(), name) ;
52a36ffd 197 else
198 rv = (AliPHOSGeometry *) fgGeom ;
199 }
200 return rv ;
201}
4697edca 202
52a36ffd 203//____________________________________________________________________________
204void AliPHOSGeometry::SetPHOSAngles()
205{
a4e98857 206 // Calculates the position of the PHOS modules in ALICE global coordinate system
91daaf24 207 // in ideal geometry
52a36ffd 208
a8c47ab6 209 Double_t const kRADDEG = 180.0 / TMath::Pi() ;
710f859a 210 Float_t pphi = 2 * TMath::ATan( GetOuterBoxSize(0) / ( 2.0 * GetIPtoUpperCPVsurface() ) ) ;
52a36ffd 211 pphi *= kRADDEG ;
710f859a 212 if (pphi > fAngle){
351dd634 213 AliError(Form("PHOS modules overlap!\n pphi = %f fAngle = %f",
214 pphi, fAngle));
710f859a 215
216 }
ed4205d8 217 pphi = fAngle;
52a36ffd 218
219 for( Int_t i = 1; i <= fNModules ; i++ ) {
ed4205d8 220 Float_t angle = pphi * ( i - fNModules / 2.0 - 0.5 ) ;
52a36ffd 221 fPHOSAngle[i-1] = - angle ;
222 }
d15a28e7 223}
9f616d61 224//____________________________________________________________________________
308fb942 225void AliPHOSGeometry::GetGlobal(const AliRecPoint* , TVector3 & ) const
d15a28e7 226{
9ee9f78d 227 AliFatal(Form("Please use GetGlobalPHOS(recPoint,gpos) instead of GetGlobal!"));
228}
229
230//____________________________________________________________________________
231void AliPHOSGeometry::GetGlobalPHOS(const AliPHOSRecPoint* recPoint, TVector3 & gpos) const
232{
a4e98857 233 // Calculates the coordinates of a RecPoint and the error matrix in the ALICE global coordinate system
b2a60966 234
9ee9f78d 235 const AliPHOSRecPoint * tmpPHOS = recPoint ;
92862013 236 TVector3 localposition ;
d15a28e7 237
238 tmpPHOS->GetLocalPosition(gpos) ;
239
91daaf24 240 if (!gGeoManager){
241 AliFatal("Geo manager not initialized\n");
d15a28e7 242 }
91daaf24 243 //construct module name
244 char path[100] ;
245 Double_t dy ;
246 if(tmpPHOS->IsEmc()){
247 sprintf(path,"/ALIC_1/PHOS_%d/PEMC_1/PCOL_1/PTIO_1/PCOR_1/PAGA_1/PTII_1",tmpPHOS->GetPHOSMod()) ;
8e20650f 248// sprintf(path,"/ALIC_1/PHOS_%d",tmpPHOS->GetPHOSMod()) ;
249 dy=fCrystalShift ;
d15a28e7 250 }
91daaf24 251 else{
252 sprintf(path,"/ALIC_1/PHOS_%d/PCPV_1",tmpPHOS->GetPHOSMod()) ;
253 dy= GetCPVBoxSize(1)/2. ; //center of CPV module
254 }
255 Double_t pos[3]={gpos.X(),gpos.Y()-dy,gpos.Z()} ;
8e20650f 256 if(tmpPHOS->IsEmc())
257 pos[2]=-pos[2] ; //Opposite z directions in EMC matrix and local frame!!!
91daaf24 258 Double_t posC[3];
259 //now apply possible shifts and rotations
260 if (!gGeoManager->cd(path)){
261 AliFatal("Geo manager can not find path \n");
262 }
263 TGeoHMatrix *m = gGeoManager->GetCurrentMatrix();
8e20650f 264 if (m){
265 m->LocalToMaster(pos,posC);
266 }
91daaf24 267 else{
268 AliFatal("Geo matrixes are not loaded \n") ;
269 }
8e20650f 270 gpos.SetXYZ(posC[0],posC[1],posC[2]) ;
7b51037f 271
aa35fc01 272}
aa35fc01 273//____________________________________________________________________________
aa0b9641 274
7b51037f 275void AliPHOSGeometry::GetModuleCenter(TVector3& center,
276 const char *det,
277 Int_t module) const
aa0b9641 278{
bfc17d18 279 // Returns a position of the center of the CPV or EMC module
91daaf24 280 // in ideal (not misaligned) geometry
1504ac47 281 Float_t rDet = 0.;
e77bb310 282 if (strcmp(det,"CPV") == 0) rDet = GetIPtoCPVDistance ();
283 else if (strcmp(det,"EMC") == 0) rDet = GetIPtoCrystalSurface();
351dd634 284 else
285 AliFatal(Form("Wrong detector name %s",det));
bfc17d18 286
395f4eea 287 Float_t angle = GetPHOSAngle(module); // (40,20,0,-20,-40) degrees
288 angle *= TMath::Pi()/180;
289 angle += 3*TMath::Pi()/2.;
7b51037f 290 center.SetXYZ(rDet*TMath::Cos(angle), rDet*TMath::Sin(angle), 0.);
aa0b9641 291}
292