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