New PWG library
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALEMCGeometry.cxx
CommitLineData
0c5b726e 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: AliEMCALEMCGeometry.cxx 29514 2008-10-26 10:24:38Z hristov $*/
17
18//_________________________________________________________________________
19// Geometry class for EMCAL : singleton
20// EMCAL consists of layers of scintillator and lead
21// with scintillator fiber arranged as "shish-kebab" skewers
22// Places the the Barrel Geometry of The EMCAL at Midrapidity
23// between 80 and 180(or 190) degrees of Phi and
24// -0.7 to 0.7 in eta
25//
26// EMCAL geometry tree:
27// EMCAL -> superModule -> module -> tower(cell)
28// Indexes
29// absId -> nSupMod -> nModule -> (nIphi,nIeta)
30//
31// Name choices:
32// EMCAL_PDC06 (geometry used for PDC06 simulations, kept for backward compatibility)
33// = equivalent to SHISH_77_TRD1_2X2_FINAL_110DEG in old notation
34// EMCAL_COMPLETE (geometry for expected complete detector)
35// = equivalent to SHISH_77_TRD1_2X2_FINAL_110DEG scTh=0.176 pbTh=0.144
36// in old notation
3d841a9f 37//
38// EMCAL_FIRSTYEAR - geometry for December 2009 to December 2010 run period
39// with four Super Modules
40//
41// EMCAL_FIRSTYEARV1 - geometry for December 2009 to December 2010 run period;
42// 1. Fixed bug for positions of modules inside SM
43// (first module has tilt 0.75 degree);
44// 2. Added Al front plate (width 1 cm) and 2 paper sheets per sampling
45// layer (additional 0.2 mm)
46// The sizes have updated with last information from production
47// drawing (end of October 2010).
48//
0c5b726e 49// EMCAL_WSUC (Wayne State test stand)
50// = no definite equivalent in old notation, was only used by
51// Aleksei, but kept for testing purposes
52//
53// etc.
54//
55//
56//
57//*-- Author: Sahal Yacoob (LBL / UCT)
58// and : Yves Schutz (SUBATECH)
59// and : Jennifer Klay (LBL)
60// and : Aleksei Pavlinov (WSU)
61// and : Magali Estienne (SUBATECH)
62
63// --- Root header files ---
64#include <TObjArray.h>
65#include <TObjString.h>
66#include <TRegexp.h>
67
68// -- ALICE Headers.
69#include "AliLog.h"
70
71// --- EMCAL headers
72#include "AliEMCALEMCGeometry.h"
73#include <cassert>
74
75ClassImp(AliEMCALEMCGeometry)
76
77// these initialisations are needed for a singleton
78Bool_t AliEMCALEMCGeometry::fgInit = kFALSE;
79const Char_t* AliEMCALEMCGeometry::fgkDefaultGeometryName = "EMCAL_COMPLETE";
80
81
82AliEMCALEMCGeometry::AliEMCALEMCGeometry()
83 : TNamed(),
84 fGeoName(0),fArrayOpts(0),fNAdditionalOpts(0),fECPbRadThickness(0.),fECScintThick(0.),
85 fNECLayers(0),fArm1PhiMin(0.),fArm1PhiMax(0.),fArm1EtaMin(0.),fArm1EtaMax(0.),fIPDistance(0.),
86 fShellThickness(0.),fZLength(0.),fNZ(0),fNPhi(0),fSampling(0.),fNumberOfSuperModules(0),
87 fFrontSteelStrip(0.),fLateralSteelStrip(0.),fPassiveScintThick(0.),fPhiModuleSize(0.),
88 fEtaModuleSize(0.),fPhiTileSize(0.),fEtaTileSize(0.),fLongModuleSize(0.),fNPhiSuperModule(0),
89 fNPHIdiv(0),fNETAdiv(0), fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),
90 // Trigger staff
91 fNTRUEta(0), fNTRUPhi(0), fNModulesInTRUEta(0), fNModulesInTRUPhi(0), fNEtaSubOfTRU(0),
92 //
93 fTrd1Angle(0.),f2Trd1Dx2(0.),
94 fPhiGapForSM(0.),fKey110DEG(0),fPhiBoundariesOfSM(0), fPhiCentersOfSM(0),fEtaMaxOfTRD1(0),
3d841a9f 95 fTrd1AlFrontThick(0.0), fTrd1BondPaperThick(0.),
0c5b726e 96 fCentersOfCellsEtaDir(0), fCentersOfCellsXDir(0),fCentersOfCellsPhiDir(0),
97 fEtaCentersOfCells(0),fPhiCentersOfCells(0),fShishKebabTrd1Modules(0),
3d841a9f 98 fParSM(), fILOSS(-1), fIHADR(-1),
0c5b726e 99 //obsolete member data
3d841a9f 100 fGap2Active(0.), fSteelFrontThick(0.), fTrd2AngleY(0.),
0c5b726e 101 f2Trd2Dy2(0.), fEmptySpace(0.), fTubsR(0.), fTubsTurnAngle(0.)
102{
103 // default ctor only for internal usage (singleton)
104 // must be kept public for root persistency purposes,
105 // but should never be called by the outside world
7e1d9a9b 106 fParSM[0]=0; fParSM[1]=0; fParSM[2]=0;
a5f0a2ba 107 fEnvelop[0] = 0; fEnvelop[1] = 0; fEnvelop[2] = 0;
a51e676d 108 for(Int_t i = 0; i < 6; i++) fkAdditionalOpts[i] = "";
a5f0a2ba 109
0c5b726e 110 AliDebug(2, "AliEMCALEMCGeometry : default ctor ");
111}
112//______________________________________________________________________
113AliEMCALEMCGeometry::AliEMCALEMCGeometry(const Text_t* name, const Text_t* title) :
114 TNamed(name,title),
115 fGeoName(0),fArrayOpts(0),fNAdditionalOpts(0),fECPbRadThickness(0.),fECScintThick(0.),
116 fNECLayers(0),fArm1PhiMin(0.),fArm1PhiMax(0.),fArm1EtaMin(0.),fArm1EtaMax(0.),fIPDistance(0.),
117 fShellThickness(0.),fZLength(0.),fNZ(0),fNPhi(0),fSampling(0.),fNumberOfSuperModules(0),
118 fFrontSteelStrip(0.),fLateralSteelStrip(0.),fPassiveScintThick(0.),fPhiModuleSize(0.),
119 fEtaModuleSize(0.),fPhiTileSize(0.),fEtaTileSize(0.),fLongModuleSize(0.),fNPhiSuperModule(0),
120 fNPHIdiv(0),fNETAdiv(0), fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),
121 // Trigger staff
122 fNTRUEta(0), fNTRUPhi(0), fNModulesInTRUEta(0), fNModulesInTRUPhi(0), fNEtaSubOfTRU(0),
123 //
124 fTrd1Angle(0.),f2Trd1Dx2(0.),
125 fPhiGapForSM(0.),fKey110DEG(0),fPhiBoundariesOfSM(0), fPhiCentersOfSM(0), fEtaMaxOfTRD1(0),
3d841a9f 126 fTrd1AlFrontThick(0.0), fTrd1BondPaperThick(0.),
0c5b726e 127 fCentersOfCellsEtaDir(0),fCentersOfCellsXDir(0),fCentersOfCellsPhiDir(0),
128 fEtaCentersOfCells(0),fPhiCentersOfCells(0),fShishKebabTrd1Modules(0),
3d841a9f 129 fParSM(),fILOSS(-1), fIHADR(-1),
0c5b726e 130 //obsolete member data
3d841a9f 131 fGap2Active(0.), fSteelFrontThick(0.), fTrd2AngleY(0.),
0c5b726e 132 f2Trd2Dy2(0.), fEmptySpace(0.), fTubsR(0.), fTubsTurnAngle(0.)
133{
134 // ctor only for internal usage (singleton)
135 AliDebug(2, Form("AliEMCALEMCGeometry(%s,%s) ", name,title));
136
137 Init();
138
139 // CreateListOfTrd1Modules();
140
141 if (AliDebugLevel()>=2) {
142 PrintGeometry();
143 }
144
145}
146//______________________________________________________________________
147AliEMCALEMCGeometry::AliEMCALEMCGeometry(const AliEMCALEMCGeometry& geom)
148 : TNamed(geom),
149 fGeoName(geom.fGeoName),
150 fArrayOpts(geom.fArrayOpts),
151 fNAdditionalOpts(geom.fNAdditionalOpts),
152 fECPbRadThickness(geom.fECPbRadThickness),
153 fECScintThick(geom.fECScintThick),
154 fNECLayers(geom.fNECLayers),
155 fArm1PhiMin(geom.fArm1PhiMin),
156 fArm1PhiMax(geom.fArm1PhiMax),
157 fArm1EtaMin(geom.fArm1EtaMin),
158 fArm1EtaMax(geom.fArm1EtaMax),
159 fIPDistance(geom.fIPDistance),
160 fShellThickness(geom.fShellThickness),
161 fZLength(geom.fZLength),
162 fNZ(geom.fNZ),
163 fNPhi(geom.fNPhi),
164 fSampling(geom.fSampling),
165 fNumberOfSuperModules(geom.fNumberOfSuperModules),
166 fFrontSteelStrip(geom.fFrontSteelStrip),
167 fLateralSteelStrip(geom.fLateralSteelStrip),
168 fPassiveScintThick(geom.fPassiveScintThick),
169 fPhiModuleSize(geom.fPhiModuleSize),
170 fEtaModuleSize(geom.fEtaModuleSize),
171 fPhiTileSize(geom.fPhiTileSize),
172 fEtaTileSize(geom.fEtaTileSize),
173 fLongModuleSize(geom.fLongModuleSize),
174 fNPhiSuperModule(geom.fNPhiSuperModule),
175 fNPHIdiv(geom.fNPHIdiv),
176 fNETAdiv(geom.fNETAdiv),
177 fNCells(geom.fNCells),
178 fNCellsInSupMod(geom.fNCellsInSupMod),
179 fNCellsInModule(geom.fNCellsInModule),
180 // Trigger staff
181 fNTRUEta(geom.fNTRUEta),
182 fNTRUPhi(geom.fNTRUPhi),
183 fNModulesInTRUEta(geom.fNModulesInTRUEta),
184 fNModulesInTRUPhi(geom.fNModulesInTRUPhi),
185 fNEtaSubOfTRU(geom.fNEtaSubOfTRU),
186 //
187 fTrd1Angle(geom.fTrd1Angle),
188 f2Trd1Dx2(geom.f2Trd1Dx2),
189 fPhiGapForSM(geom.fPhiGapForSM),
190 fKey110DEG(geom.fKey110DEG),
191 fPhiBoundariesOfSM(geom.fPhiBoundariesOfSM),
192 fPhiCentersOfSM(geom.fPhiCentersOfSM),
193 fEtaMaxOfTRD1(geom.fEtaMaxOfTRD1),
3d841a9f 194 fTrd1AlFrontThick(geom.fTrd1AlFrontThick),
195 fTrd1BondPaperThick(geom.fTrd1BondPaperThick),
0c5b726e 196 fCentersOfCellsEtaDir(geom.fCentersOfCellsEtaDir),
197 fCentersOfCellsXDir(geom.fCentersOfCellsXDir),
198 fCentersOfCellsPhiDir(geom.fCentersOfCellsPhiDir),
199 fEtaCentersOfCells(geom.fEtaCentersOfCells),
200 fPhiCentersOfCells(geom.fPhiCentersOfCells),
201 fShishKebabTrd1Modules(geom.fShishKebabTrd1Modules),
202 fILOSS(geom.fILOSS), fIHADR(geom.fIHADR),
203 //obsolete member data
0c5b726e 204 fGap2Active(geom.fGap2Active),
205 fSteelFrontThick(geom.fSteelFrontThick),
206 fTrd2AngleY(geom.fTrd2AngleY),
207 f2Trd2Dy2(geom.f2Trd2Dy2),
208 fEmptySpace(geom.fEmptySpace),
209 fTubsR(geom.fTubsR),
210 fTubsTurnAngle(geom.fTubsTurnAngle)
211{
212 //copy ctor
7e1d9a9b 213 fParSM[0]=geom.fParSM[0];
214 fParSM[1]=geom.fParSM[1];
215 fParSM[2]=geom.fParSM[2];
a5f0a2ba 216 fEnvelop[0] = geom.fEnvelop[0];
217 fEnvelop[1] = geom.fEnvelop[1];
218 fEnvelop[2] = geom.fEnvelop[2];
a51e676d 219 for(Int_t i = 0; i < 6; i++) fkAdditionalOpts[i] = geom.fkAdditionalOpts[i];
a5f0a2ba 220
0c5b726e 221}
222
223//______________________________________________________________________
224AliEMCALEMCGeometry::~AliEMCALEMCGeometry(void){
225 // dtor
226}
227
228//______________________________________________________________________
229void AliEMCALEMCGeometry::Init(void){
230 //
231 // Initializes the EMCAL parameters based on the name
232 // Only Shashlyk geometry is available, but various combinations of
233 // layers and number of supermodules can be selected with additional
234 // options or geometry name
235 //
236
237 fkAdditionalOpts[0] = "nl="; // number of sampling layers (fNECLayers)
238 fkAdditionalOpts[1] = "pbTh="; // cm, Thickness of the Pb (fECPbRadThick)
239 fkAdditionalOpts[2] = "scTh="; // cm, Thickness of the Sc (fECScintThick)
240 fkAdditionalOpts[3] = "latSS="; // cm, Thickness of lateral steel strip (fLateralSteelStrip)
241 fkAdditionalOpts[4] = "allILOSS="; // = 0,1,2,3,4 (4 - energy loss without fluctuation)
242 fkAdditionalOpts[5] = "allIHADR="; // = 0,1,2 (0 - no hadronic interaction)
243
244 fNAdditionalOpts = sizeof(fkAdditionalOpts) / sizeof(char*);
245
246 // geometry
247 fgInit = kFALSE; // Assume failed until proven otherwise.
248 fGeoName = GetName();
249 fGeoName.ToUpper();
250
251 //Convert old geometry names to new ones
252 if(fGeoName.Contains("SHISH_77_TRD1_2X2_FINAL_110DEG")) {
253 if(fGeoName.Contains("PBTH=0.144") && fGeoName.Contains("SCTH=0.176")) {
254 fGeoName = "EMCAL_COMPLETE";
255 } else {
256 fGeoName = "EMCAL_PDC06";
257 }
258 }
259 if(fGeoName.Contains("WSUC")) fGeoName = "EMCAL_WSUC";
260
261 //check that we have a valid geometry name
1d806785 262 if(!(fGeoName.Contains("EMCAL_PDC06") || fGeoName.Contains("EMCAL_COMPLETE") || fGeoName.Contains("EMCAL_WSUC")
263 || fGeoName.Contains("EMCAL_FIRSTYEAR") || fGeoName.Contains("EMCAL_FIRSTYEARV1"))) {
0c5b726e 264 Fatal("Init", "%s is an undefined geometry!", fGeoName.Data()) ;
265 }
266
267 // Option to know whether we have the "half" supermodule(s) or not
268 fKey110DEG = 0;
269 if(fGeoName.Contains("COMPLETE") || fGeoName.Contains("PDC06")) fKey110DEG = 1; // for GetAbsCellId
270 fShishKebabTrd1Modules = 0;
271
272 // JLK 13-Apr-2008
273 //default parameters are those of EMCAL_COMPLETE geometry
274 //all others render variations from these at the end of
275 //geometry-name specific options
276
277 fNumberOfSuperModules = 12; // 12 = 6 * 2 (6 in phi, 2 in Z)
278 fNPhi = 12; // module granularity in phi within smod (azimuth)
279 fNZ = 24; // module granularity along Z within smod (eta)
280 fNPHIdiv = fNETAdiv = 2; // tower granularity within module
281 fArm1PhiMin = 80.0; // degrees, Starting EMCAL Phi position
282 fArm1PhiMax = 200.0; // degrees, Ending EMCAL Phi position
283 fArm1EtaMin = -0.7; // pseudorapidity, Starting EMCAL Eta position
284 fArm1EtaMax = +0.7; // pseudorapidity, Ending EMCAL Eta position
285 fIPDistance = 428.0; // cm, radial distance to front face from nominal vertex point
286 fPhiGapForSM = 2.; // cm, only for final TRD1 geometry
287 fFrontSteelStrip = 0.025; // 0.025cm = 0.25mm (13-may-05 from V.Petrov)
288 fPassiveScintThick = 0.8; // 0.8cm = 8mm (13-may-05 from V.Petrov)
289 fLateralSteelStrip = 0.01; // 0.01cm = 0.1mm (13-may-05 from V.Petrov) - was 0.025
290 fTrd1Angle = 1.5; // in degrees
291
292 fSampling = 1.; // should be calculated with call to DefineSamplingFraction()
293 fNECLayers = 77; // (13-may-05 from V.Petrov) - can be changed with additional options
294 fECScintThick = 0.176; // scintillator layer thickness
295 fECPbRadThickness = 0.144; // lead layer thickness
296
297 fPhiModuleSize = 12.26 - fPhiGapForSM / Float_t(fNPhi); // first assumption
298 fEtaModuleSize = fPhiModuleSize;
299
300 fZLength = 700.; // Z coverage (cm)
301
302
303 //needs to be called for each geometry and before setting geometry
304 //parameters which can depend on the outcome
305 CheckAdditionalOptions();
306
307 //modifications to the above for PDC06 geometry
308 if(fGeoName.Contains("PDC06")){ // 18-may-05 - about common structure
309 fECScintThick = fECPbRadThickness = 0.16;// (13-may-05 from V.Petrov)
310 CheckAdditionalOptions();
311 }
312
313 //modifications to the above for WSUC geometry
314 if(fGeoName.Contains("WSUC")){ // 18-may-05 - about common structure
315 fPhiModuleSize = 12.5; // 20-may-05 - rectangular shape
316 fEtaModuleSize = 11.9;
317 fECScintThick = fECPbRadThickness = 0.16;// (13-may-05 from V.Petrov)
318 fNumberOfSuperModules = 1; // 27-may-05
319 fShellThickness = 30.; // should be change
320 fNPhi = fNZ = 4;
321 CheckAdditionalOptions();
322 }
323
9f467289 324 //In 2009-2010 data taking runs only 4 SM, in the upper position.
325 if(fGeoName.Contains("FIRSTYEAR")){
3d841a9f 326 fNumberOfSuperModules = 4;
327 fArm1PhiMax = 120.0;
328 if(fGeoName.Contains("FIRSTYEARV1")){
329 // Oct 26,2010 : First module has tilt = 0.75 degree :
330 // look to AliEMCALShishKebabTrd1Module::DefineFirstModule(key)
331 // New sizes from production drawing, added Al front plate.
332 // The thickness of sampling is change due to existing two sheets of paper.
333
334 // Will replace fFrontSteelStrip
335 fTrd1AlFrontThick = 1.0; // one cm
336 // Bond paper - two sheets around Sc tile
337 fTrd1BondPaperThick = 0.01; // 0.01cm = 0.1 mm
338
339 fPhiModuleSize = 12.0;
340 fEtaModuleSize = fPhiModuleSize;
341 fLateralSteelStrip = 0.015; // 0.015cm = 0.15mm
342 }
343 CheckAdditionalOptions();
9f467289 344 }
345
0c5b726e 346 // constant for transition absid <--> indexes
9f467289 347 fNCellsInModule = fNPHIdiv*fNETAdiv;
0c5b726e 348 fNCellsInSupMod = fNCellsInModule*fNPhi*fNZ;
349 fNCells = fNCellsInSupMod*fNumberOfSuperModules;
350 if(GetKey110DEG()) fNCells -= fNCellsInSupMod;
351
352 fNPhiSuperModule = fNumberOfSuperModules/2;
353 if(fNPhiSuperModule < 1) fNPhiSuperModule = 1;
354
355 fPhiTileSize = fPhiModuleSize/double(fNPHIdiv) - fLateralSteelStrip; // 13-may-05
356 fEtaTileSize = fEtaModuleSize/double(fNETAdiv) - fLateralSteelStrip; // 13-may-05
357
358 fLongModuleSize = fNECLayers*(fECScintThick + fECPbRadThickness);
3d841a9f 359 if(fGeoName.Contains("FIRSTYEARV1")){
360 Double_t ws = fECScintThick + fECPbRadThickness + 2.*fTrd1BondPaperThick; // sampling width
361 // Number of Pb tiles = Number of Sc tiles - 1
362 fLongModuleSize = fTrd1AlFrontThick + (ws*fNECLayers - fECPbRadThickness);
363 }
0c5b726e 364 f2Trd1Dx2 = fEtaModuleSize + 2.*fLongModuleSize*TMath::Tan(fTrd1Angle*TMath::DegToRad()/2.);
3d841a9f 365
0c5b726e 366 if(!fGeoName.Contains("WSUC")) fShellThickness = TMath::Sqrt(fLongModuleSize*fLongModuleSize + f2Trd1Dx2*f2Trd1Dx2);
367
368 //These parameters are used to create the mother volume to hold the supermodules
369 //2cm padding added to allow for misalignments - JLK 30-May-2008
370 fEnvelop[0] = fIPDistance - 1.; // mother volume inner radius
371 fEnvelop[1] = fIPDistance + fShellThickness + 1.; // mother volume outer r.
372 fEnvelop[2] = fZLength + 2.; //mother volume length
373
374 // Local coordinates
375 fParSM[0] = GetShellThickness()/2.;
376 fParSM[1] = GetPhiModuleSize() * GetNPhi()/2.;
377 fParSM[2] = fZLength/4.; //divide by 4 to get half-length of SM
378
379 // SM phi boundaries - (0,1),(2,3) .. (10,11) - has the same boundaries; Nov 7, 2006
380 fPhiBoundariesOfSM.Set(fNumberOfSuperModules);
381 fPhiCentersOfSM.Set(fNumberOfSuperModules/2);
382 fPhiBoundariesOfSM[0] = TMath::PiOver2() - TMath::ATan2(fParSM[1] , fIPDistance); // 1th and 2th modules)
383 fPhiCentersOfSM[0] = TMath::PiOver2();
384 if(fNumberOfSuperModules > 1)
385 fPhiBoundariesOfSM[1] = TMath::PiOver2() + TMath::ATan2(fParSM[1] , fIPDistance);
386 if(fNumberOfSuperModules > 2) {
9f467289 387 Int_t maxPhiBlock =fNumberOfSuperModules/2-1;
388 if(fNumberOfSuperModules > 10) maxPhiBlock = 4;
389 for(int i=1; i<=maxPhiBlock; i++) { // from 2th ro 9th
0c5b726e 390 fPhiBoundariesOfSM[2*i] = fPhiBoundariesOfSM[0] + 20.*TMath::DegToRad()*i;
391 fPhiBoundariesOfSM[2*i+1] = fPhiBoundariesOfSM[1] + 20.*TMath::DegToRad()*i;
9f467289 392 fPhiCentersOfSM[i] = fPhiCentersOfSM[0] + 20.*TMath::DegToRad()*i;
0c5b726e 393 }
394 }
395 if(fNumberOfSuperModules > 10) {
396 fPhiBoundariesOfSM[11] = 190.*TMath::DegToRad();
397 fPhiBoundariesOfSM[10] = fPhiBoundariesOfSM[11] - TMath::ATan2((fParSM[1]) , fIPDistance);
9f467289 398 fPhiCentersOfSM[5] = (fPhiBoundariesOfSM[10]+fPhiBoundariesOfSM[11])/2.;
0c5b726e 399 }
400
401 //called after setting of scintillator and lead layer parameters
402 DefineSamplingFraction();
403
404
405 // TRU parameters - Apr 29,08 by PAI.
406 // These parameters values was updated at Nov 05, 2007
407 // As is on Olivier BOURRION (LPSC) ppt preasentation
408 // at ALICE trigger meeting at 13th-14th March
409 fNTRUEta = 1; // was 3
410 fNTRUPhi = 3; // was 1
411 fNModulesInTRUEta = 24; // was 8
412 fNModulesInTRUPhi = 4; // was 12
413 // Jet trigger
414 // 3*6*10 + 2*6*2 = 204 -> matrix (nphi(17), neta(12))
415 fNEtaSubOfTRU = 6;
0c5b726e 416
417 fgInit = kTRUE;
418}
419
420//___________________________________________________________________
421void AliEMCALEMCGeometry::PrintGeometry()
422{
423 // Separate routine is callable from broswer; Nov 7,2006
424 printf("\nInit: geometry of EMCAL named %s :\n", fGeoName.Data());
425 if(fArrayOpts) {
426 for(Int_t i=0; i<fArrayOpts->GetEntries(); i++){
427 TObjString *o = (TObjString*)fArrayOpts->At(i);
428 printf(" %i : %s \n", i, o->String().Data());
429 }
430 }
431 printf("Granularity: %d in eta and %d in phi\n", GetNZ(), GetNPhi()) ;
432 printf("Layout: phi = (%7.1f, %7.1f), eta = (%5.2f, %5.2f), IP = %7.2f -> for EMCAL envelope only\n",
433 GetArm1PhiMin(), GetArm1PhiMax(),GetArm1EtaMin(), GetArm1EtaMax(), GetIPDistance() );
434
435 printf( " ECAL : %d x (%f cm Pb, %f cm Sc) \n",
436 GetNECLayers(), GetECPbRadThick(), GetECScintThick() ) ;
437 printf(" fSampling %5.2f \n", fSampling );
438 printf(" fIPDistance %6.3f cm \n", fIPDistance);
439 printf(" fNPhi %i | fNZ %i \n", fNPhi, fNZ);
440 printf(" fNCellsInModule %i : fNCellsInSupMod %i : fNCells %i\n",fNCellsInModule, fNCellsInSupMod, fNCells);
441 printf(" X:Y module size %6.3f , %6.3f cm \n", fPhiModuleSize, fEtaModuleSize);
442 printf(" X:Y tile size %6.3f , %6.3f cm \n", fPhiTileSize, fEtaTileSize);
443 printf(" #of sampling layers %i(fNECLayers) \n", fNECLayers);
444 printf(" fLongModuleSize %6.3f cm \n", fLongModuleSize);
445 printf(" #supermodule in phi direction %i \n", fNPhiSuperModule );
446 printf(" fILOSS %i : fIHADR %i \n", fILOSS, fIHADR);
447 printf(" fTrd1Angle %7.4f\n", fTrd1Angle);
448 printf(" f2Trd1Dx2 %7.4f\n", f2Trd1Dx2);
3d841a9f 449 printf(" fTrd1AlFrontThick %7.4f \n", fTrd1AlFrontThick);
450 printf(" fTrd1BondPaperThick %5.4f \n", fTrd1BondPaperThick);
0c5b726e 451 printf("SM dimensions(TRD1) : dx %7.2f dy %7.2f dz %7.2f (SMOD, BOX)\n",
452 fParSM[0],fParSM[1],fParSM[2]);
453 printf(" fPhiGapForSM %7.4f cm (%7.4f <- phi size in degree)\n",
454 fPhiGapForSM, TMath::ATan2(fPhiGapForSM,fIPDistance)*TMath::RadToDeg());
455 if(GetKey110DEG()) printf(" Last two modules have size 10 degree in phi (180<phi<190)\n");
456 printf(" phi SM boundaries \n");
457 for(int i=0; i<fPhiBoundariesOfSM.GetSize()/2.; i++) {
458 printf(" %i : %7.5f(%7.2f) -> %7.5f(%7.2f) : center %7.5f(%7.2f) \n", i,
459 fPhiBoundariesOfSM[2*i], fPhiBoundariesOfSM[2*i]*TMath::RadToDeg(),
460 fPhiBoundariesOfSM[2*i+1], fPhiBoundariesOfSM[2*i+1]*TMath::RadToDeg(),
461 fPhiCentersOfSM[i], fPhiCentersOfSM[i]*TMath::RadToDeg());
462 }
463
464}
465
466//______________________________________________________________________
467void AliEMCALEMCGeometry::CheckAdditionalOptions()
468{
469 // Feb 06,2006
470 // Additional options that
471 // can be used to select
472 // the specific geometry of
473 // EMCAL to run
474 // Dec 27,2006
475 // adeed allILOSS= and allIHADR= for MIP investigation
476 fArrayOpts = new TObjArray;
477 Int_t nopt = ParseString(fGeoName, *fArrayOpts);
478 if(nopt==1) { // no aditional option(s)
479 fArrayOpts->Delete();
480 delete fArrayOpts;
481 fArrayOpts = 0;
482 return;
483 }
484 for(Int_t i=1; i<nopt; i++){
485 TObjString *o = (TObjString*)fArrayOpts->At(i);
486
487 TString addOpt = o->String();
488 Int_t indj=-1;
489 for(Int_t j=0; j<fNAdditionalOpts; j++) {
490 TString opt = fkAdditionalOpts[j];
491 if(addOpt.Contains(opt,TString::kIgnoreCase)) {
492 indj = j;
493 break;
494 }
495 }
496 if(indj<0) {
497 AliDebug(2,Form("<E> option |%s| unavailable : ** look to the file AliEMCALGeometry.h **\n",
498 addOpt.Data()));
499 assert(0);
500 } else {
501 AliDebug(2,Form("<I> option |%s| is valid : number %i : |%s|\n",
502 addOpt.Data(), indj, fkAdditionalOpts[indj]));
503 if (addOpt.Contains("NL=",TString::kIgnoreCase)) {// number of sampling layers
504 sscanf(addOpt.Data(),"NL=%i", &fNECLayers);
505 AliDebug(2,Form(" fNECLayers %i (new) \n", fNECLayers));
506 } else if(addOpt.Contains("PBTH=",TString::kIgnoreCase)) {//Thickness of the Pb(fECPbRadThicknes)
507 sscanf(addOpt.Data(),"PBTH=%f", &fECPbRadThickness);
508 } else if(addOpt.Contains("SCTH=",TString::kIgnoreCase)) {//Thickness of the Sc(fECScintThick)
509 sscanf(addOpt.Data(),"SCTH=%f", &fECScintThick);
510 } else if(addOpt.Contains("LATSS=",TString::kIgnoreCase)) {// Thickness of lateral steel strip (fLateralSteelStrip)
511 sscanf(addOpt.Data(),"LATSS=%f", &fLateralSteelStrip);
512 AliDebug(2,Form(" fLateralSteelStrip %f (new) \n", fLateralSteelStrip));
513 } else if(addOpt.Contains("ILOSS=",TString::kIgnoreCase)) {// As in Geant
514 sscanf(addOpt.Data(),"ALLILOSS=%i", &fILOSS);
515 AliDebug(2,Form(" fILOSS %i \n", fILOSS));
516 } else if(addOpt.Contains("IHADR=",TString::kIgnoreCase)) {// As in Geant
517 sscanf(addOpt.Data(),"ALLIHADR=%i", &fIHADR);
518 AliDebug(2,Form(" fIHADR %i \n", fIHADR));
519 }
520 }
521 }
522}
523
524//__________________________________________________________________
525void AliEMCALEMCGeometry::DefineSamplingFraction()
526{
527 // Jun 05,2006
528 // Look http://rhic.physics.wayne.edu/~pavlinov/ALICE/SHISHKEBAB/RES/linearityAndResolutionForTRD1.html
529 // Keep for compatibilty
530 //
531 if(fNECLayers == 69) { // 10% layer reduction
532 fSampling = 12.55;
533 } else if(fNECLayers == 61) { // 20% layer reduction
534 fSampling = 12.80;
535 } else if(fNECLayers == 77) {
536 if (fECScintThick>0.159 && fECScintThick<0.161) { // original sampling fraction, equal layers
537 fSampling = 12.327; // fECScintThick = fECPbRadThickness = 0.160;
538 } else if (fECScintThick>0.175 && fECScintThick<0.177) { // 10% Pb thicknes reduction
539 fSampling = 10.5; // fECScintThick = 0.176, fECPbRadThickness=0.144;
540 } else if(fECScintThick>0.191 && fECScintThick<0.193) { // 20% Pb thicknes reduction
541 fSampling = 8.93; // fECScintThick = 0.192, fECPbRadThickness=0.128;
542 }
543
544 }
545}
546
547//________________________________________________________________________________________________
548Double_t AliEMCALEMCGeometry::GetPhiCenterOfSM(Int_t nsupmod) const
549{
550 //returns center of supermodule in phi
551 int i = nsupmod/2;
552 return fPhiCentersOfSM[i];
553
554}
555
556//________________________________________________________________________________________________
557Bool_t AliEMCALEMCGeometry::GetPhiBoundariesOfSM(Int_t nSupMod, Double_t &phiMin, Double_t &phiMax) const
558{
559 // 0<= nSupMod <=11; phi in rad
560 static int i;
561 if(nSupMod<0 || nSupMod >11) return kFALSE;
562 i = nSupMod/2;
563 phiMin = (Double_t)fPhiBoundariesOfSM[2*i];
564 phiMax = (Double_t)fPhiBoundariesOfSM[2*i+1];
565 return kTRUE;
566}
567
568//________________________________________________________________________________________________
569Bool_t AliEMCALEMCGeometry::GetPhiBoundariesOfSMGap(Int_t nPhiSec, Double_t &phiMin, Double_t &phiMax) const
570{
571 // 0<= nPhiSec <=4; phi in rad
572 // 0; gap boundaries between 0th&2th | 1th&3th SM
573 // 1; gap boundaries between 2th&4th | 3th&5th SM
574 // 2; gap boundaries between 4th&6th | 5th&7th SM
575 // 3; gap boundaries between 6th&8th | 7th&9th SM
576 // 4; gap boundaries between 8th&10th | 9th&11th SM
577 if(nPhiSec<0 || nPhiSec >4) return kFALSE;
578 phiMin = fPhiBoundariesOfSM[2*nPhiSec+1];
579 phiMax = fPhiBoundariesOfSM[2*nPhiSec+2];
580 return kTRUE;
581}
582
583//________________________________________________________________________________________________
584int AliEMCALEMCGeometry::ParseString(const TString &topt, TObjArray &Opt)
585{
586 //Parse string, does what? GCB 08/09
587 Ssiz_t begin, index, end, end2;
588 begin = index = end = end2 = 0;
589 TRegexp separator("[^ ;,\\t\\s/]+");
590 while ( (begin < topt.Length()) && (index != kNPOS) ) {
591 // loop over given options
592 index = topt.Index(separator,&end,begin);
593 if (index >= 0 && end >= 1) {
594 TString substring(topt(index,end));
595 Opt.Add(new TObjString(substring.Data()));
596 }
597 begin += end+1;
598 }
599 return Opt.GetEntries();
600}
601