]> git.uio.no Git - u/mrichter/AliRoot.git/blame - EMCAL/AliEMCALEMCGeometry.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Create wrapper for JSON macro that calls traditional macro, and add development JSON...
[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//
171d2441 41// Adding V1 (EMCAL_FIRSTYEARV1, EMCAL_COMPLETEV1) - geometry from December 2009 ;
3d841a9f 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).
d7f5c01a 48// 3. COMPLETEV1 contains now only 10 SM for runs for year 2011
49// 4. COMPLETE12SMV1 contains 12 SM for runs from year 2012 and on
8cc543cb 50// 5. COMPLETE12SMV1_DCAL contains 12 SM and 6 DCal SM
51// 6. COMPLETE12SMV1_DCAL_DEV contains 12 SM and 10 DCal SM
52// 7. COMPLETE12SMV1_DCAL_8SM contains 12 SM and 6 DCal SM and 2 extentions
3d841a9f 53//
0c5b726e 54// EMCAL_WSUC (Wayne State test stand)
55// = no definite equivalent in old notation, was only used by
56// Aleksei, but kept for testing purposes
57//
58// etc.
59//
60//
61//
62//*-- Author: Sahal Yacoob (LBL / UCT)
63// and : Yves Schutz (SUBATECH)
64// and : Jennifer Klay (LBL)
65// and : Aleksei Pavlinov (WSU)
66// and : Magali Estienne (SUBATECH)
8cc543cb 67// and : Adapted for DCAL by M.L. Wang CCNU Wuhan & Subatech Oct-23-2009
0c5b726e 68
69// --- Root header files ---
70#include <TObjArray.h>
71#include <TObjString.h>
72#include <TRegexp.h>
73
74// -- ALICE Headers.
75#include "AliLog.h"
76
77// --- EMCAL headers
78#include "AliEMCALEMCGeometry.h"
79#include <cassert>
80
81ClassImp(AliEMCALEMCGeometry)
82
83// these initialisations are needed for a singleton
84Bool_t AliEMCALEMCGeometry::fgInit = kFALSE;
d7f5c01a 85const Char_t* AliEMCALEMCGeometry::fgkDefaultGeometryName = "EMCAL_COMPLETE12SMV1";
0c5b726e 86
87
88AliEMCALEMCGeometry::AliEMCALEMCGeometry()
89 : TNamed(),
8cc543cb 90 fGeoName(0),fEMCSMSystem(0x0),fArrayOpts(0),fNAdditionalOpts(0),fECPbRadThickness(0.),fECScintThick(0.),
0c5b726e 91 fNECLayers(0),fArm1PhiMin(0.),fArm1PhiMax(0.),fArm1EtaMin(0.),fArm1EtaMax(0.),fIPDistance(0.),
8cc543cb 92 fShellThickness(0.),fZLength(0.),fDCALInnerEdge(0.),fDCALPhiMin(0),fDCALPhiMax(0),fEMCALPhiMax(0),
93 fDCALStandardPhiMax(0),fDCALInnerExtandedEta(0),fNZ(0),fNPhi(0),fSampling(0.),fNumberOfSuperModules(0),
0c5b726e 94 fFrontSteelStrip(0.),fLateralSteelStrip(0.),fPassiveScintThick(0.),fPhiModuleSize(0.),
8cc543cb 95 fEtaModuleSize(0.),fPhiTileSize(0.),fEtaTileSize(0.),fLongModuleSize(0.),fPhiSuperModule(0),fNPhiSuperModule(0),
0c5b726e 96 fNPHIdiv(0),fNETAdiv(0), fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),
97 // Trigger staff
8cc543cb 98 fNTRUEta(0), fNTRUPhi(0), fNModulesInTRUEta(0), fNModulesInTRUPhi(0), fNEtaSubOfTRU(0), fNTotalTRU(0),
0c5b726e 99 //
8cc543cb 100 fTrd1Angle(0.),f2Trd1Dx2(0.),fPhiGapForSM(0.),fKey110DEG(0),fnSupModInDCAL(0),fPhiBoundariesOfSM(0),
101 fPhiCentersOfSM(0),fPhiCentersOfSMSec(0),fEtaMaxOfTRD1(0),fTrd1AlFrontThick(0.0), fTrd1BondPaperThick(0.),
0c5b726e 102 fCentersOfCellsEtaDir(0), fCentersOfCellsXDir(0),fCentersOfCellsPhiDir(0),
103 fEtaCentersOfCells(0),fPhiCentersOfCells(0),fShishKebabTrd1Modules(0),
3d841a9f 104 fParSM(), fILOSS(-1), fIHADR(-1),
0c5b726e 105 //obsolete member data
3d841a9f 106 fGap2Active(0.), fSteelFrontThick(0.), fTrd2AngleY(0.),
0c5b726e 107 f2Trd2Dy2(0.), fEmptySpace(0.), fTubsR(0.), fTubsTurnAngle(0.)
108{
109 // default ctor only for internal usage (singleton)
110 // must be kept public for root persistency purposes,
111 // but should never be called by the outside world
7e1d9a9b 112 fParSM[0]=0; fParSM[1]=0; fParSM[2]=0;
a5f0a2ba 113 fEnvelop[0] = 0; fEnvelop[1] = 0; fEnvelop[2] = 0;
a51e676d 114 for(Int_t i = 0; i < 6; i++) fkAdditionalOpts[i] = "";
a5f0a2ba 115
0c5b726e 116 AliDebug(2, "AliEMCALEMCGeometry : default ctor ");
117}
118//______________________________________________________________________
126215a5 119AliEMCALEMCGeometry::AliEMCALEMCGeometry(const Text_t* name, const Text_t* title,
120 const Text_t* mcname, const Text_t* mctitle ) :
0c5b726e 121 TNamed(name,title),
8cc543cb 122 fGeoName(0),fEMCSMSystem(0x0),fArrayOpts(0),fNAdditionalOpts(0),fECPbRadThickness(0.),fECScintThick(0.),
0c5b726e 123 fNECLayers(0),fArm1PhiMin(0.),fArm1PhiMax(0.),fArm1EtaMin(0.),fArm1EtaMax(0.),fIPDistance(0.),
8cc543cb 124 fShellThickness(0.),fZLength(0.),fDCALInnerEdge(0.),fDCALPhiMin(0),fDCALPhiMax(0),fEMCALPhiMax(0),
125 fDCALStandardPhiMax(0),fDCALInnerExtandedEta(0),fNZ(0),fNPhi(0),fSampling(0.),fNumberOfSuperModules(0),
0c5b726e 126 fFrontSteelStrip(0.),fLateralSteelStrip(0.),fPassiveScintThick(0.),fPhiModuleSize(0.),
8cc543cb 127 fEtaModuleSize(0.),fPhiTileSize(0.),fEtaTileSize(0.),fLongModuleSize(0.),fPhiSuperModule(0),fNPhiSuperModule(0),
0c5b726e 128 fNPHIdiv(0),fNETAdiv(0), fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),
129 // Trigger staff
8cc543cb 130 fNTRUEta(0), fNTRUPhi(0), fNModulesInTRUEta(0), fNModulesInTRUPhi(0), fNEtaSubOfTRU(0), fNTotalTRU(0),
0c5b726e 131 //
8cc543cb 132 fTrd1Angle(0.),f2Trd1Dx2(0.),fPhiGapForSM(0.),fKey110DEG(0),fnSupModInDCAL(0),fPhiBoundariesOfSM(0),
133 fPhiCentersOfSM(0),fPhiCentersOfSMSec(0),fEtaMaxOfTRD1(0),fTrd1AlFrontThick(0.0), fTrd1BondPaperThick(0.),
0c5b726e 134 fCentersOfCellsEtaDir(0),fCentersOfCellsXDir(0),fCentersOfCellsPhiDir(0),
135 fEtaCentersOfCells(0),fPhiCentersOfCells(0),fShishKebabTrd1Modules(0),
3d841a9f 136 fParSM(),fILOSS(-1), fIHADR(-1),
0c5b726e 137 //obsolete member data
3d841a9f 138 fGap2Active(0.), fSteelFrontThick(0.), fTrd2AngleY(0.),
0c5b726e 139 f2Trd2Dy2(0.), fEmptySpace(0.), fTubsR(0.), fTubsTurnAngle(0.)
140{
141 // ctor only for internal usage (singleton)
126215a5 142 AliDebug(2, Form("AliEMCALEMCGeometry(%s,%s,%s,%s) ", name,title,mcname,mctitle));
0c5b726e 143
126215a5 144 Init(mcname,mctitle);
0c5b726e 145
146 // CreateListOfTrd1Modules();
147
148 if (AliDebugLevel()>=2) {
149 PrintGeometry();
150 }
151
152}
153//______________________________________________________________________
154AliEMCALEMCGeometry::AliEMCALEMCGeometry(const AliEMCALEMCGeometry& geom)
155 : TNamed(geom),
156 fGeoName(geom.fGeoName),
8cc543cb 157 fEMCSMSystem(geom.fEMCSMSystem),
0c5b726e 158 fArrayOpts(geom.fArrayOpts),
159 fNAdditionalOpts(geom.fNAdditionalOpts),
160 fECPbRadThickness(geom.fECPbRadThickness),
161 fECScintThick(geom.fECScintThick),
162 fNECLayers(geom.fNECLayers),
163 fArm1PhiMin(geom.fArm1PhiMin),
164 fArm1PhiMax(geom.fArm1PhiMax),
165 fArm1EtaMin(geom.fArm1EtaMin),
166 fArm1EtaMax(geom.fArm1EtaMax),
167 fIPDistance(geom.fIPDistance),
168 fShellThickness(geom.fShellThickness),
169 fZLength(geom.fZLength),
8cc543cb 170 fDCALInnerEdge(geom.fDCALInnerEdge),
171 fDCALPhiMin(geom.fDCALPhiMin),
172 fDCALPhiMax(geom.fDCALPhiMax),
173 fEMCALPhiMax(geom.fEMCALPhiMax),
174 fDCALStandardPhiMax(geom.fDCALStandardPhiMax),
175 fDCALInnerExtandedEta(geom.fDCALInnerExtandedEta),
0c5b726e 176 fNZ(geom.fNZ),
177 fNPhi(geom.fNPhi),
178 fSampling(geom.fSampling),
179 fNumberOfSuperModules(geom.fNumberOfSuperModules),
180 fFrontSteelStrip(geom.fFrontSteelStrip),
181 fLateralSteelStrip(geom.fLateralSteelStrip),
182 fPassiveScintThick(geom.fPassiveScintThick),
183 fPhiModuleSize(geom.fPhiModuleSize),
184 fEtaModuleSize(geom.fEtaModuleSize),
185 fPhiTileSize(geom.fPhiTileSize),
186 fEtaTileSize(geom.fEtaTileSize),
187 fLongModuleSize(geom.fLongModuleSize),
8cc543cb 188 fPhiSuperModule(geom.fPhiSuperModule),
0c5b726e 189 fNPhiSuperModule(geom.fNPhiSuperModule),
190 fNPHIdiv(geom.fNPHIdiv),
191 fNETAdiv(geom.fNETAdiv),
192 fNCells(geom.fNCells),
193 fNCellsInSupMod(geom.fNCellsInSupMod),
194 fNCellsInModule(geom.fNCellsInModule),
195 // Trigger staff
196 fNTRUEta(geom.fNTRUEta),
197 fNTRUPhi(geom.fNTRUPhi),
198 fNModulesInTRUEta(geom.fNModulesInTRUEta),
199 fNModulesInTRUPhi(geom.fNModulesInTRUPhi),
200 fNEtaSubOfTRU(geom.fNEtaSubOfTRU),
8cc543cb 201 fNTotalTRU(geom.fNTotalTRU),
0c5b726e 202 //
203 fTrd1Angle(geom.fTrd1Angle),
204 f2Trd1Dx2(geom.f2Trd1Dx2),
205 fPhiGapForSM(geom.fPhiGapForSM),
206 fKey110DEG(geom.fKey110DEG),
8cc543cb 207 fnSupModInDCAL(geom.fnSupModInDCAL),
0c5b726e 208 fPhiBoundariesOfSM(geom.fPhiBoundariesOfSM),
209 fPhiCentersOfSM(geom.fPhiCentersOfSM),
8cc543cb 210 fPhiCentersOfSMSec(geom.fPhiCentersOfSMSec),
0c5b726e 211 fEtaMaxOfTRD1(geom.fEtaMaxOfTRD1),
3d841a9f 212 fTrd1AlFrontThick(geom.fTrd1AlFrontThick),
213 fTrd1BondPaperThick(geom.fTrd1BondPaperThick),
0c5b726e 214 fCentersOfCellsEtaDir(geom.fCentersOfCellsEtaDir),
215 fCentersOfCellsXDir(geom.fCentersOfCellsXDir),
216 fCentersOfCellsPhiDir(geom.fCentersOfCellsPhiDir),
217 fEtaCentersOfCells(geom.fEtaCentersOfCells),
218 fPhiCentersOfCells(geom.fPhiCentersOfCells),
219 fShishKebabTrd1Modules(geom.fShishKebabTrd1Modules),
220 fILOSS(geom.fILOSS), fIHADR(geom.fIHADR),
221 //obsolete member data
0c5b726e 222 fGap2Active(geom.fGap2Active),
223 fSteelFrontThick(geom.fSteelFrontThick),
224 fTrd2AngleY(geom.fTrd2AngleY),
225 f2Trd2Dy2(geom.f2Trd2Dy2),
226 fEmptySpace(geom.fEmptySpace),
227 fTubsR(geom.fTubsR),
228 fTubsTurnAngle(geom.fTubsTurnAngle)
229{
230 //copy ctor
7e1d9a9b 231 fParSM[0]=geom.fParSM[0];
232 fParSM[1]=geom.fParSM[1];
233 fParSM[2]=geom.fParSM[2];
a5f0a2ba 234 fEnvelop[0] = geom.fEnvelop[0];
235 fEnvelop[1] = geom.fEnvelop[1];
236 fEnvelop[2] = geom.fEnvelop[2];
a51e676d 237 for(Int_t i = 0; i < 6; i++) fkAdditionalOpts[i] = geom.fkAdditionalOpts[i];
a5f0a2ba 238
0c5b726e 239}
240
241//______________________________________________________________________
242AliEMCALEMCGeometry::~AliEMCALEMCGeometry(void){
243 // dtor
8cc543cb 244 delete fEMCSMSystem;
0c5b726e 245}
246
247//______________________________________________________________________
126215a5 248void AliEMCALEMCGeometry::Init(const Text_t* mcname, const Text_t* mctitle){
0c5b726e 249 //
250 // Initializes the EMCAL parameters based on the name
251 // Only Shashlyk geometry is available, but various combinations of
252 // layers and number of supermodules can be selected with additional
253 // options or geometry name
254 //
255
256 fkAdditionalOpts[0] = "nl="; // number of sampling layers (fNECLayers)
257 fkAdditionalOpts[1] = "pbTh="; // cm, Thickness of the Pb (fECPbRadThick)
258 fkAdditionalOpts[2] = "scTh="; // cm, Thickness of the Sc (fECScintThick)
259 fkAdditionalOpts[3] = "latSS="; // cm, Thickness of lateral steel strip (fLateralSteelStrip)
260 fkAdditionalOpts[4] = "allILOSS="; // = 0,1,2,3,4 (4 - energy loss without fluctuation)
261 fkAdditionalOpts[5] = "allIHADR="; // = 0,1,2 (0 - no hadronic interaction)
262
263 fNAdditionalOpts = sizeof(fkAdditionalOpts) / sizeof(char*);
264
265 // geometry
266 fgInit = kFALSE; // Assume failed until proven otherwise.
267 fGeoName = GetName();
268 fGeoName.ToUpper();
269
270 //Convert old geometry names to new ones
271 if(fGeoName.Contains("SHISH_77_TRD1_2X2_FINAL_110DEG")) {
272 if(fGeoName.Contains("PBTH=0.144") && fGeoName.Contains("SCTH=0.176")) {
273 fGeoName = "EMCAL_COMPLETE";
274 } else {
275 fGeoName = "EMCAL_PDC06";
276 }
277 }
d7f5c01a 278
0c5b726e 279 if(fGeoName.Contains("WSUC")) fGeoName = "EMCAL_WSUC";
280
281 //check that we have a valid geometry name
8cc543cb 282 if(!( fGeoName.Contains("EMCAL_PDC06")
283 || fGeoName.Contains("EMCAL_WSUC")
284 || fGeoName.Contains("EMCAL_COMPLETE")
285 || fGeoName.Contains("EMCAL_COMPLETEV1")
286 || fGeoName.Contains("EMCAL_COMPLETE12SMV1")
287 || fGeoName.Contains("EMCAL_FIRSTYEAR")
288 || fGeoName.Contains("EMCAL_FIRSTYEARV1") )) {
0c5b726e 289 Fatal("Init", "%s is an undefined geometry!", fGeoName.Data()) ;
290 }
291
292 // Option to know whether we have the "half" supermodule(s) or not
293 fKey110DEG = 0;
8cc543cb 294 if( fGeoName.Contains("COMPLETE")
295 || fGeoName.Contains("PDC06")
296 || fGeoName.Contains("12SM")) fKey110DEG = 1; // for GetAbsCellId
171d2441 297 if(fGeoName.Contains("COMPLETEV1")) fKey110DEG = 0;
0c5b726e 298 fShishKebabTrd1Modules = 0;
299
8cc543cb 300 fnSupModInDCAL = 0;
301 if(fGeoName.Contains("DCAL_DEV")){
302 fnSupModInDCAL = 10;
303 } else if(fGeoName.Contains("DCAL_8SM")){
304 fnSupModInDCAL = 8;
305 } else if(fGeoName.Contains("DCAL")){
306 fnSupModInDCAL = 6;
307 }
308
0c5b726e 309 // JLK 13-Apr-2008
310 //default parameters are those of EMCAL_COMPLETE geometry
311 //all others render variations from these at the end of
312 //geometry-name specific options
313
314 fNumberOfSuperModules = 12; // 12 = 6 * 2 (6 in phi, 2 in Z)
8cc543cb 315 fNPhi = 12; // module granularity in phi within smod (azimuth)
0c5b726e 316 fNZ = 24; // module granularity along Z within smod (eta)
317 fNPHIdiv = fNETAdiv = 2; // tower granularity within module
318 fArm1PhiMin = 80.0; // degrees, Starting EMCAL Phi position
319 fArm1PhiMax = 200.0; // degrees, Ending EMCAL Phi position
320 fArm1EtaMin = -0.7; // pseudorapidity, Starting EMCAL Eta position
321 fArm1EtaMax = +0.7; // pseudorapidity, Ending EMCAL Eta position
322 fIPDistance = 428.0; // cm, radial distance to front face from nominal vertex point
323 fPhiGapForSM = 2.; // cm, only for final TRD1 geometry
324 fFrontSteelStrip = 0.025; // 0.025cm = 0.25mm (13-may-05 from V.Petrov)
325 fPassiveScintThick = 0.8; // 0.8cm = 8mm (13-may-05 from V.Petrov)
326 fLateralSteelStrip = 0.01; // 0.01cm = 0.1mm (13-may-05 from V.Petrov) - was 0.025
327 fTrd1Angle = 1.5; // in degrees
328
329 fSampling = 1.; // should be calculated with call to DefineSamplingFraction()
330 fNECLayers = 77; // (13-may-05 from V.Petrov) - can be changed with additional options
331 fECScintThick = 0.176; // scintillator layer thickness
332 fECPbRadThickness = 0.144; // lead layer thickness
333
334 fPhiModuleSize = 12.26 - fPhiGapForSM / Float_t(fNPhi); // first assumption
335 fEtaModuleSize = fPhiModuleSize;
336
337 fZLength = 700.; // Z coverage (cm)
8cc543cb 338 fPhiSuperModule = 20. ; // phi in degree
339 fDCALInnerEdge = fIPDistance * TMath::Tan( fTrd1Angle * 8.* TMath::DegToRad());
0c5b726e 340
341 //needs to be called for each geometry and before setting geometry
342 //parameters which can depend on the outcome
343 CheckAdditionalOptions();
344
345 //modifications to the above for PDC06 geometry
346 if(fGeoName.Contains("PDC06")){ // 18-may-05 - about common structure
347 fECScintThick = fECPbRadThickness = 0.16;// (13-may-05 from V.Petrov)
348 CheckAdditionalOptions();
349 }
350
351 //modifications to the above for WSUC geometry
352 if(fGeoName.Contains("WSUC")){ // 18-may-05 - about common structure
2d7a9f00 353 fNumberOfSuperModules = 2; // 27-may-05; Nov 24,2010 for TB
0c5b726e 354 fNPhi = fNZ = 4;
2d7a9f00 355 fTrd1AlFrontThick = 1.0; // one cm
356 // Bond paper - two sheets around Sc tile
357 fTrd1BondPaperThick = 0.01; // 0.01cm = 0.1 mm
358
359 fPhiModuleSize = 12.0;
360 fEtaModuleSize = fPhiModuleSize;
361 fLateralSteelStrip = 0.015; // 0.015cm = 0.15mm
362
0c5b726e 363 CheckAdditionalOptions();
364 }
365
9f467289 366 //In 2009-2010 data taking runs only 4 SM, in the upper position.
367 if(fGeoName.Contains("FIRSTYEAR")){
3d841a9f 368 fNumberOfSuperModules = 4;
8cc543cb 369 fArm1PhiMax = 120.0;
3d841a9f 370 CheckAdditionalOptions();
8cc543cb 371 }
171d2441 372
d7f5c01a 373 if(fGeoName.Contains("FIRSTYEARV1") || fGeoName.Contains("COMPLETEV1") || fGeoName.Contains("COMPLETE12SMV1") ){
171d2441 374 // Oct 26,2010 : First module has tilt = 0.75 degree :
375 // look to AliEMCALShishKebabTrd1Module::DefineFirstModule(key)
376 // New sizes from production drawing, added Al front plate.
377 // The thickness of sampling is change due to existing two sheets of paper.
378
379 // Will replace fFrontSteelStrip
380 fTrd1AlFrontThick = 1.0; // one cm
381 // Bond paper - two sheets around Sc tile
382 fTrd1BondPaperThick = 0.01; // 0.01cm = 0.1 mm
383
384 fPhiModuleSize = 12.0;
385 fEtaModuleSize = fPhiModuleSize;
386 fLateralSteelStrip = 0.015; // 0.015cm = 0.15mm
387
d7f5c01a 388 if(fGeoName.Contains("COMPLETEV1"))
389 {
171d2441 390 fNumberOfSuperModules = 10;
8cc543cb 391 fArm1PhiMax = 180.0;
171d2441 392 }
d7f5c01a 393 else if (fGeoName.Contains("COMPLETE12SMV1"))
394 {
395 fNumberOfSuperModules = 12;
396 fArm1PhiMax = 200.0;
397 }
8cc543cb 398 if (fGeoName.Contains("DCAL"))
399 {
400 fNumberOfSuperModules = 12 + fnSupModInDCAL;
401 fArm1PhiMax = 320.0;
402 if(fGeoName.Contains("DCAL_8SM")) fArm1PhiMax = 340.0; // degrees, End of DCAL Phi position
403 else if(fGeoName.Contains("DCAL_DEV")) fArm1PhiMin = 40.0; // degrees, Starting EMCAL(shifted) Phi position
404 fDCALPhiMin = fArm1PhiMax - 10.*fnSupModInDCAL;
405 }
171d2441 406 CheckAdditionalOptions();
407 }
408
8cc543cb 409 fEMCSMSystem = new Int_t[fNumberOfSuperModules];
410 Int_t iSM = 0;
411
412 // BASIC EMCAL SM
413 if(fGeoName.Contains("WSUC") ){
414 for(int i = 0; i<2; i++){
415 fEMCSMSystem[iSM] = kEMCAL_Standard;
416 iSM++;
417 }
418 } else if(fGeoName.Contains("FIRSTYEAR") ){
419 for(int i = 0; i<4; i++){
420 fEMCSMSystem[iSM] = kEMCAL_Standard;
421 iSM++;
422 }
423 } else if( fGeoName.Contains("PDC06")
424 || fGeoName.Contains("COMPLETE") ){
425 for(int i = 0; i<10; i++){
426 fEMCSMSystem[iSM] = kEMCAL_Standard;
427 iSM++;
428 }
429 }
430 // EMCAL 110SM
431 if(fKey110DEG && fGeoName.Contains("12SM") ){
432 for(int i = 0; i<2; i++){
433 fEMCSMSystem[iSM] = kEMCAL_Half;
434 if(fGeoName.Contains("12SMV1") ){
435 fEMCSMSystem[iSM] = kEMCAL_3rd;
436 }
437 iSM++;
438 }
439 }
440 // DCAL SM
441 if(fnSupModInDCAL && fGeoName.Contains("DCAL")){
442 if(fGeoName.Contains("8SM")) {
443 for(int i = 0; i<fnSupModInDCAL-2; i++){
444 fEMCSMSystem[iSM] = kDCAL_Standard;
445 iSM++;
446 }
447 for(int i = 0; i<2; i++){
448 fEMCSMSystem[iSM] = kDCAL_Ext;
449 iSM++;
450 }
451 } else {
452 for(int i = 0; i<fnSupModInDCAL; i++){
453 fEMCSMSystem[iSM] = kDCAL_Standard;
454 iSM++;
455 }
456 }
457 }
458
0c5b726e 459 // constant for transition absid <--> indexes
9f467289 460 fNCellsInModule = fNPHIdiv*fNETAdiv;
0c5b726e 461 fNCellsInSupMod = fNCellsInModule*fNPhi*fNZ;
8cc543cb 462 fNCells = 0;
463 for( int i=0; i<fNumberOfSuperModules; i++) {
464 if( GetSMType(i) == kEMCAL_Standard) fNCells += fNCellsInSupMod ;
465 else if( GetSMType(i) == kEMCAL_Half) fNCells += fNCellsInSupMod/2 ;
466 else if( GetSMType(i) == kEMCAL_3rd) fNCells += fNCellsInSupMod/3 ;
467 else if( GetSMType(i) == kDCAL_Standard) fNCells += 2*fNCellsInSupMod/3 ;
468 else if( GetSMType(i) == kDCAL_Ext) fNCells += fNCellsInSupMod/3 ;
469 else AliError(Form("Uknown SuperModule Type !!"));
470 }
0c5b726e 471
472 fNPhiSuperModule = fNumberOfSuperModules/2;
473 if(fNPhiSuperModule < 1) fNPhiSuperModule = 1;
474
475 fPhiTileSize = fPhiModuleSize/double(fNPHIdiv) - fLateralSteelStrip; // 13-may-05
476 fEtaTileSize = fEtaModuleSize/double(fNETAdiv) - fLateralSteelStrip; // 13-may-05
477
478 fLongModuleSize = fNECLayers*(fECScintThick + fECPbRadThickness);
171d2441 479 if(fGeoName.Contains("V1")){
3d841a9f 480 Double_t ws = fECScintThick + fECPbRadThickness + 2.*fTrd1BondPaperThick; // sampling width
481 // Number of Pb tiles = Number of Sc tiles - 1
482 fLongModuleSize = fTrd1AlFrontThick + (ws*fNECLayers - fECPbRadThickness);
483 }
0c5b726e 484 f2Trd1Dx2 = fEtaModuleSize + 2.*fLongModuleSize*TMath::Tan(fTrd1Angle*TMath::DegToRad()/2.);
3d841a9f 485
0c5b726e 486 if(!fGeoName.Contains("WSUC")) fShellThickness = TMath::Sqrt(fLongModuleSize*fLongModuleSize + f2Trd1Dx2*f2Trd1Dx2);
487
488 //These parameters are used to create the mother volume to hold the supermodules
489 //2cm padding added to allow for misalignments - JLK 30-May-2008
490 fEnvelop[0] = fIPDistance - 1.; // mother volume inner radius
491 fEnvelop[1] = fIPDistance + fShellThickness + 1.; // mother volume outer r.
492 fEnvelop[2] = fZLength + 2.; //mother volume length
493
494 // Local coordinates
495 fParSM[0] = GetShellThickness()/2.;
496 fParSM[1] = GetPhiModuleSize() * GetNPhi()/2.;
497 fParSM[2] = fZLength/4.; //divide by 4 to get half-length of SM
498
8cc543cb 499 // SM phi boundaries - (0,1),(2,3) ... - has the same boundaries;
0c5b726e 500 fPhiBoundariesOfSM.Set(fNumberOfSuperModules);
501 fPhiCentersOfSM.Set(fNumberOfSuperModules/2);
8cc543cb 502 fPhiCentersOfSMSec.Set(fNumberOfSuperModules/2);
503 Double_t kfSupermodulePhiWidth = fPhiSuperModule*TMath::DegToRad();
504 fPhiCentersOfSM[0] = (fArm1PhiMin + fPhiSuperModule/2.) * TMath::DegToRad(); // Define from First SM
505 fPhiCentersOfSMSec[0] = fPhiCentersOfSM[0]; // the same in the First SM
506 fPhiBoundariesOfSM[0] = fPhiCentersOfSM[0] - TMath::ATan2(fParSM[1] , fIPDistance); // 1th and 2th modules)
507 fPhiBoundariesOfSM[1] = fPhiCentersOfSM[0] + TMath::ATan2(fParSM[1] , fIPDistance);
508 if(fNumberOfSuperModules > 2) { // 2 to Max
509 Int_t tmpSMType = GetSMType(2);
510 for(int i = 1; i<fNPhiSuperModule; i++) {
511 fPhiBoundariesOfSM[2*i] += fPhiBoundariesOfSM[2*i-2] + kfSupermodulePhiWidth;
512 if(tmpSMType == GetSMType(2*i)) {
513 fPhiBoundariesOfSM[2*i+1] += fPhiBoundariesOfSM[2*i-1] + kfSupermodulePhiWidth;
514 } else {
515 //changed SM Type, redefine the [2*i+1] Boundaries
516 tmpSMType = GetSMType(2*i);
517 if( GetSMType(2*i) == kEMCAL_Standard) {
518 fPhiBoundariesOfSM[2*i+1] = fPhiBoundariesOfSM[2*i] + kfSupermodulePhiWidth;
519 } else if( GetSMType(2*i) == kEMCAL_Half) {
520 fPhiBoundariesOfSM[2*i+1] = fPhiBoundariesOfSM[2*i] + 2.*TMath::ATan2((fParSM[1])/2, fIPDistance);
521 } else if( GetSMType(2*i) == kEMCAL_3rd ) {
522 fPhiBoundariesOfSM[2*i+1] = fPhiBoundariesOfSM[2*i] + 2.*TMath::ATan2((fParSM[1])/3, fIPDistance);
523 } else if( GetSMType(2*i) == kDCAL_Standard ) { // jump the gap
524 fPhiBoundariesOfSM[2*i] = (fDCALPhiMin - fArm1PhiMin)*TMath::DegToRad() + fPhiBoundariesOfSM[0];
525 fPhiBoundariesOfSM[2*i+1] = (fDCALPhiMin - fArm1PhiMin)*TMath::DegToRad() + fPhiBoundariesOfSM[1];
526 } else if( GetSMType(2*i) == kDCAL_Ext) {
527 fPhiBoundariesOfSM[2*i+1] = fPhiBoundariesOfSM[2*i] + 2.*TMath::ATan2((fParSM[1])/3, fIPDistance);
528 }
529 }
530 fPhiCentersOfSM[i] = (fPhiBoundariesOfSM[2*i]+fPhiBoundariesOfSM[2*i+1])/2.;
531 fPhiCentersOfSMSec[i] = fPhiBoundariesOfSM[2*i] + TMath::ATan2(fParSM[1] , fIPDistance);
0c5b726e 532 }
533 }
0c5b726e 534
8cc543cb 535 //inner extend in eta (same as outer part) for DCal (0.189917), //calculated from the smallest gap (1# cell to the 80-degree-edge),
536 Double_t fInnerExtandedPhi = 1.102840997; //calculated from the smallest gap (1# cell to the 80-degree-edge), too complicatd to explain...
537 fDCALInnerExtandedEta = -TMath::Log(TMath::Tan( (TMath::Pi()/2. - 8*fTrd1Angle*TMath::DegToRad() + (TMath::Pi()/2 - fNZ*fTrd1Angle*TMath::DegToRad() - TMath::ATan(TMath::Exp(fArm1EtaMin))*2))/2.));
538
539 fEMCALPhiMax = fArm1PhiMin;
540 fDCALPhiMax = fDCALPhiMin;// DCAl extention will not be included
541 for( Int_t i = 0; i < fNumberOfSuperModules; i+=2) {
542 if( GetSMType(i) == kEMCAL_Standard ) fEMCALPhiMax += 20.;
543 else if( GetSMType(i) == kEMCAL_Half ) fEMCALPhiMax += fPhiSuperModule/2. + fInnerExtandedPhi;
544 else if( GetSMType(i) == kEMCAL_3rd ) fEMCALPhiMax += fPhiSuperModule/3. + 4.0*fInnerExtandedPhi/3.0;
545 else if( GetSMType(i) == kDCAL_Standard ) {fDCALPhiMax += 20.; fDCALStandardPhiMax = fDCALPhiMax;}
546 else if( GetSMType(i) == kDCAL_Ext ) fDCALPhiMax += fPhiSuperModule/3. + 4.0*fInnerExtandedPhi/3.0;
547 else AliError("Unkown SM Type!!");
548 }
549// for compatible reason
550// if(fNumberOfSuperModules == 4) {fEMCALPhiMax = fArm1PhiMax ;}
551 if(fNumberOfSuperModules == 12) {fEMCALPhiMax = fArm1PhiMax ;}
552
0c5b726e 553 //called after setting of scintillator and lead layer parameters
126215a5 554 DefineSamplingFraction(mcname,mctitle);
0c5b726e 555
556
557 // TRU parameters - Apr 29,08 by PAI.
558 // These parameters values was updated at Nov 05, 2007
559 // As is on Olivier BOURRION (LPSC) ppt preasentation
560 // at ALICE trigger meeting at 13th-14th March
561 fNTRUEta = 1; // was 3
562 fNTRUPhi = 3; // was 1
563 fNModulesInTRUEta = 24; // was 8
564 fNModulesInTRUPhi = 4; // was 12
565 // Jet trigger
566 // 3*6*10 + 2*6*2 = 204 -> matrix (nphi(17), neta(12))
567 fNEtaSubOfTRU = 6;
0c5b726e 568
8cc543cb 569 fNTotalTRU = 0;
570 for(Int_t i = 0; i < GetNumberOfSuperModules(); i++){
571 if( GetSMType(i) == kEMCAL_Standard) fNTotalTRU += 3;
572 else if( GetSMType(i) == kEMCAL_Half) fNTotalTRU += 1;
573 else if( GetSMType(i) == kEMCAL_3rd) fNTotalTRU += 1;
574 else if( GetSMType(i) == kDCAL_Standard) fNTotalTRU += 3;
575 else if( GetSMType(i) == kDCAL_Ext) fNTotalTRU += 1;
576 else {
577 AliError(Form("Uknown SuperModule Type !!"));
578 }
579 }
580
0c5b726e 581 fgInit = kTRUE;
582}
583
584//___________________________________________________________________
585void AliEMCALEMCGeometry::PrintGeometry()
586{
587 // Separate routine is callable from broswer; Nov 7,2006
588 printf("\nInit: geometry of EMCAL named %s :\n", fGeoName.Data());
589 if(fArrayOpts) {
590 for(Int_t i=0; i<fArrayOpts->GetEntries(); i++){
591 TObjString *o = (TObjString*)fArrayOpts->At(i);
592 printf(" %i : %s \n", i, o->String().Data());
593 }
594 }
8cc543cb 595 if(fGeoName.Contains("DCAL")) {
596 printf("Phi min of DCAL SuperModule: %7.1f, DCAL has %d SuperModule\n", fDCALPhiMin, fnSupModInDCAL);
597 printf("The DCAL inner edge is +- %7.1f\n", fDCALInnerEdge);
598 if(fGeoName.Contains("DCAL_8SM")) printf("DCAL has its 2 EXTENTION SM\n");
599 }
0c5b726e 600 printf("Granularity: %d in eta and %d in phi\n", GetNZ(), GetNPhi()) ;
601 printf("Layout: phi = (%7.1f, %7.1f), eta = (%5.2f, %5.2f), IP = %7.2f -> for EMCAL envelope only\n",
602 GetArm1PhiMin(), GetArm1PhiMax(),GetArm1EtaMin(), GetArm1EtaMax(), GetIPDistance() );
603
604 printf( " ECAL : %d x (%f cm Pb, %f cm Sc) \n",
605 GetNECLayers(), GetECPbRadThick(), GetECScintThick() ) ;
606 printf(" fSampling %5.2f \n", fSampling );
607 printf(" fIPDistance %6.3f cm \n", fIPDistance);
608 printf(" fNPhi %i | fNZ %i \n", fNPhi, fNZ);
609 printf(" fNCellsInModule %i : fNCellsInSupMod %i : fNCells %i\n",fNCellsInModule, fNCellsInSupMod, fNCells);
610 printf(" X:Y module size %6.3f , %6.3f cm \n", fPhiModuleSize, fEtaModuleSize);
611 printf(" X:Y tile size %6.3f , %6.3f cm \n", fPhiTileSize, fEtaTileSize);
612 printf(" #of sampling layers %i(fNECLayers) \n", fNECLayers);
613 printf(" fLongModuleSize %6.3f cm \n", fLongModuleSize);
614 printf(" #supermodule in phi direction %i \n", fNPhiSuperModule );
8cc543cb 615 printf(" supermodule width in phi direction %f \n", fPhiSuperModule );
0c5b726e 616 printf(" fILOSS %i : fIHADR %i \n", fILOSS, fIHADR);
617 printf(" fTrd1Angle %7.4f\n", fTrd1Angle);
618 printf(" f2Trd1Dx2 %7.4f\n", f2Trd1Dx2);
3d841a9f 619 printf(" fTrd1AlFrontThick %7.4f \n", fTrd1AlFrontThick);
620 printf(" fTrd1BondPaperThick %5.4f \n", fTrd1BondPaperThick);
0c5b726e 621 printf("SM dimensions(TRD1) : dx %7.2f dy %7.2f dz %7.2f (SMOD, BOX)\n",
622 fParSM[0],fParSM[1],fParSM[2]);
623 printf(" fPhiGapForSM %7.4f cm (%7.4f <- phi size in degree)\n",
624 fPhiGapForSM, TMath::ATan2(fPhiGapForSM,fIPDistance)*TMath::RadToDeg());
8cc543cb 625 if( fKey110DEG && !fGeoName.Contains("12SMV1") ) printf(" Last two modules have size 10 degree in phi (180<phi<190)\n");
626 if( fKey110DEG && fGeoName.Contains("12SMV1")) printf(" Last two modules have size 6.6 degree in phi (180<phi<186.6)\n");
0c5b726e 627 printf(" phi SM boundaries \n");
628 for(int i=0; i<fPhiBoundariesOfSM.GetSize()/2.; i++) {
8cc543cb 629 printf(" %i : %7.15f(%7.12f) -> %7.15f(%7.12f) : center %7.15f(%7.12f) \n", i,
0c5b726e 630 fPhiBoundariesOfSM[2*i], fPhiBoundariesOfSM[2*i]*TMath::RadToDeg(),
631 fPhiBoundariesOfSM[2*i+1], fPhiBoundariesOfSM[2*i+1]*TMath::RadToDeg(),
632 fPhiCentersOfSM[i], fPhiCentersOfSM[i]*TMath::RadToDeg());
633 }
634
635}
636
637//______________________________________________________________________
638void AliEMCALEMCGeometry::CheckAdditionalOptions()
639{
640 // Feb 06,2006
641 // Additional options that
642 // can be used to select
643 // the specific geometry of
644 // EMCAL to run
645 // Dec 27,2006
646 // adeed allILOSS= and allIHADR= for MIP investigation
647 fArrayOpts = new TObjArray;
648 Int_t nopt = ParseString(fGeoName, *fArrayOpts);
649 if(nopt==1) { // no aditional option(s)
650 fArrayOpts->Delete();
651 delete fArrayOpts;
652 fArrayOpts = 0;
653 return;
654 }
655 for(Int_t i=1; i<nopt; i++){
656 TObjString *o = (TObjString*)fArrayOpts->At(i);
657
658 TString addOpt = o->String();
659 Int_t indj=-1;
660 for(Int_t j=0; j<fNAdditionalOpts; j++) {
661 TString opt = fkAdditionalOpts[j];
662 if(addOpt.Contains(opt,TString::kIgnoreCase)) {
663 indj = j;
664 break;
665 }
666 }
667 if(indj<0) {
668 AliDebug(2,Form("<E> option |%s| unavailable : ** look to the file AliEMCALGeometry.h **\n",
669 addOpt.Data()));
670 assert(0);
671 } else {
672 AliDebug(2,Form("<I> option |%s| is valid : number %i : |%s|\n",
673 addOpt.Data(), indj, fkAdditionalOpts[indj]));
674 if (addOpt.Contains("NL=",TString::kIgnoreCase)) {// number of sampling layers
675 sscanf(addOpt.Data(),"NL=%i", &fNECLayers);
676 AliDebug(2,Form(" fNECLayers %i (new) \n", fNECLayers));
677 } else if(addOpt.Contains("PBTH=",TString::kIgnoreCase)) {//Thickness of the Pb(fECPbRadThicknes)
678 sscanf(addOpt.Data(),"PBTH=%f", &fECPbRadThickness);
679 } else if(addOpt.Contains("SCTH=",TString::kIgnoreCase)) {//Thickness of the Sc(fECScintThick)
680 sscanf(addOpt.Data(),"SCTH=%f", &fECScintThick);
681 } else if(addOpt.Contains("LATSS=",TString::kIgnoreCase)) {// Thickness of lateral steel strip (fLateralSteelStrip)
682 sscanf(addOpt.Data(),"LATSS=%f", &fLateralSteelStrip);
683 AliDebug(2,Form(" fLateralSteelStrip %f (new) \n", fLateralSteelStrip));
684 } else if(addOpt.Contains("ILOSS=",TString::kIgnoreCase)) {// As in Geant
685 sscanf(addOpt.Data(),"ALLILOSS=%i", &fILOSS);
686 AliDebug(2,Form(" fILOSS %i \n", fILOSS));
687 } else if(addOpt.Contains("IHADR=",TString::kIgnoreCase)) {// As in Geant
688 sscanf(addOpt.Data(),"ALLIHADR=%i", &fIHADR);
689 AliDebug(2,Form(" fIHADR %i \n", fIHADR));
690 }
691 }
692 }
693}
694
695//__________________________________________________________________
126215a5 696void AliEMCALEMCGeometry::DefineSamplingFraction(const Text_t* mcname, const Text_t* mctitle)
0c5b726e 697{
698 // Jun 05,2006
699 // Look http://rhic.physics.wayne.edu/~pavlinov/ALICE/SHISHKEBAB/RES/linearityAndResolutionForTRD1.html
700 // Keep for compatibilty
701 //
126215a5 702
703 // Sampling factor for G3
71e3c6ea 704 fSampling = 10.87; // Default value - Nov 25,2010
0c5b726e 705 if(fNECLayers == 69) { // 10% layer reduction
706 fSampling = 12.55;
707 } else if(fNECLayers == 61) { // 20% layer reduction
708 fSampling = 12.80;
709 } else if(fNECLayers == 77) {
171d2441 710 if(fGeoName.Contains("V1")){
71e3c6ea 711 fSampling = 10.87; //Adding paper sheets and cover plate; Nov 25,2010
171d2441 712 } else if (fECScintThick>0.159 && fECScintThick<0.161) { // original sampling fraction, equal layers
0c5b726e 713 fSampling = 12.327; // fECScintThick = fECPbRadThickness = 0.160;
714 } else if (fECScintThick>0.175 && fECScintThick<0.177) { // 10% Pb thicknes reduction
715 fSampling = 10.5; // fECScintThick = 0.176, fECPbRadThickness=0.144;
716 } else if(fECScintThick>0.191 && fECScintThick<0.193) { // 20% Pb thicknes reduction
717 fSampling = 8.93; // fECScintThick = 0.192, fECPbRadThickness=0.128;
718 }
0c5b726e 719 }
126215a5 720
721 // Default sampling factor for G3, modify it for other transport model
722 TString mcName = mcname;
723 TString mcTitle = mctitle;
724
725 Float_t samplingFactorTranportModel = 1. ;
726 if (mcName.Contains("Geant3")) samplingFactorTranportModel = 1.;//0.988 // Do nothing
727 else if(mcName.Contains("Fluka") ) samplingFactorTranportModel = 1.; // To be set
728 else if(mcName.Contains("Geant4")){
729 if(mcTitle.Contains("EMV")) samplingFactorTranportModel = 1.096; // 0.906, 0.896 (OPT)
730 else samplingFactorTranportModel = 0.86; // 1.15 (CHIPS), 1.149 (BERT), 1.147 (BERT_CHIPS)
731 }
732
733 AliDebug(2,Form("MC modeler <%s>, Title <%s>: Sampling %f, model fraction with respect to G3 %f, final sampling %f \n",
734 mcName.Data(),mcTitle.Data(),fSampling,samplingFactorTranportModel,fSampling*samplingFactorTranportModel));
735
736
737 fSampling*=samplingFactorTranportModel;
738
0c5b726e 739}
740
8cc543cb 741//________________________________________________________________________________________________
742Double_t AliEMCALEMCGeometry::GetPhiCenterOfSMSec(Int_t nsupmod) const
743{
744 //returns center of supermodule in phi
745 int i = nsupmod/2;
746 return fPhiCentersOfSMSec[i];
747
748}
749
0c5b726e 750//________________________________________________________________________________________________
751Double_t AliEMCALEMCGeometry::GetPhiCenterOfSM(Int_t nsupmod) const
752{
753 //returns center of supermodule in phi
754 int i = nsupmod/2;
755 return fPhiCentersOfSM[i];
756
757}
758
759//________________________________________________________________________________________________
760Bool_t AliEMCALEMCGeometry::GetPhiBoundariesOfSM(Int_t nSupMod, Double_t &phiMin, Double_t &phiMax) const
761{
8cc543cb 762 // 0<= nSupMod <=17; phi in rad
763 static int i;
764 if(nSupMod<0 || nSupMod >12+fnSupModInDCAL-1) return kFALSE;
0c5b726e 765 i = nSupMod/2;
766 phiMin = (Double_t)fPhiBoundariesOfSM[2*i];
767 phiMax = (Double_t)fPhiBoundariesOfSM[2*i+1];
768 return kTRUE;
769}
770
771//________________________________________________________________________________________________
772Bool_t AliEMCALEMCGeometry::GetPhiBoundariesOfSMGap(Int_t nPhiSec, Double_t &phiMin, Double_t &phiMax) const
773{
8cc543cb 774 // 0<= nPhiSec <=max; phi in rad
0c5b726e 775 // 0; gap boundaries between 0th&2th | 1th&3th SM
776 // 1; gap boundaries between 2th&4th | 3th&5th SM
777 // 2; gap boundaries between 4th&6th | 5th&7th SM
778 // 3; gap boundaries between 6th&8th | 7th&9th SM
779 // 4; gap boundaries between 8th&10th | 9th&11th SM
8cc543cb 780 // 5; gap boundaries between 10th&12th | 11h&13th SM
781 // ...
782 if(nPhiSec<0 || nPhiSec >5+fnSupModInDCAL/2-1) return kFALSE;
0c5b726e 783 phiMin = fPhiBoundariesOfSM[2*nPhiSec+1];
784 phiMax = fPhiBoundariesOfSM[2*nPhiSec+2];
785 return kTRUE;
786}
787
788//________________________________________________________________________________________________
789int AliEMCALEMCGeometry::ParseString(const TString &topt, TObjArray &Opt)
790{
791 //Parse string, does what? GCB 08/09
792 Ssiz_t begin, index, end, end2;
793 begin = index = end = end2 = 0;
794 TRegexp separator("[^ ;,\\t\\s/]+");
795 while ( (begin < topt.Length()) && (index != kNPOS) ) {
796 // loop over given options
797 index = topt.Index(separator,&end,begin);
798 if (index >= 0 && end >= 1) {
799 TString substring(topt(index,end));
800 Opt.Add(new TObjString(substring.Data()));
801 }
802 begin += end+1;
803 }
804 return Opt.GetEntries();
805}
806