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