1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
19 ////////////////////////////////////////////////////////////////
20 // This class initializes the class AliITSgeom
21 // The initialization is done starting from
22 // a geometry coded by means of the ROOT geometrical modeler
23 // This initialization can be used both for simulation and reconstruction
24 ///////////////////////////////////////////////////////////////
28 #include <TStopwatch.h>
29 #include <TGeoManager.h>
30 #include <TGeoVolume.h>
31 #include <TGeoShape.h>
38 #include <TGeoSphere.h>
47 #include "AliITSgeomSPD.h"
48 #include "AliITSgeomSDD.h"
49 #include "AliITSgeomSSD.h"
50 #include "AliITSsegmentationSPD.h"
51 #include "AliITSsegmentationSDD.h"
52 #include "AliITSsegmentationSSD.h"
53 #include "AliITSgeom.h"
54 #include "AliITSInitGeometry.h"
57 ClassImp(AliITSInitGeometry)
59 const Bool_t AliITSInitGeometry::fgkOldSPDbarrel = kTRUE;
60 const Bool_t AliITSInitGeometry::fgkOldSDDbarrel = kFALSE;
61 const Bool_t AliITSInitGeometry::fgkOldSSDbarrel = kFALSE;
62 const Bool_t AliITSInitGeometry::fgkOldSDDcone = kTRUE;
63 const Bool_t AliITSInitGeometry::fgkOldSSDcone = kTRUE;
64 const Bool_t AliITSInitGeometry::fgkOldSPDshield = kTRUE;
65 const Bool_t AliITSInitGeometry::fgkOldSDDshield = kTRUE;
66 const Bool_t AliITSInitGeometry::fgkOldSSDshield = kTRUE;
67 const Bool_t AliITSInitGeometry::fgkOldServices = kTRUE;
68 const Bool_t AliITSInitGeometry::fgkOldSupports = kTRUE;
69 //______________________________________________________________________
70 AliITSInitGeometry::AliITSInitGeometry():
71 TObject(), // Base Class
72 fName(0), // Geometry name
73 fMinorVersion(-1), // Minor version number/type
74 fMajorVersion(kvDefault), // Major versin number
75 fTiming(kFALSE), // Flag to start inilization timing
76 fSegGeom(kFALSE), // Flag to switch between the old use of
77 // AliITSgeomS?D class, or AliITSsegmentation
78 // class in fShape of AliITSgeom class.
79 fDecode(kFALSE), // Flag for new/old decoding
80 fDebug(0){ // Debug flag
87 // A default inilized AliITSInitGeometry object
91 //______________________________________________________________________
92 AliITSInitGeometry::AliITSInitGeometry(AliITSVersion_t version,
94 TObject(), // Base Class
95 fName(0), // Geometry name
96 fMinorVersion(minorversion), // Minor version number/type
97 fMajorVersion(version), // Major versin number
98 fTiming(kFALSE), // Flag to start inilization timing
99 fSegGeom(kFALSE), // Flag to switch between the old use of
100 // AliITSgeomS?D class, or AliITSsegmentation
101 // class in fShape of AliITSgeom class.
102 fDecode(kFALSE), // Flag for new/old decoding
103 fDebug(0){ // Debug flag
110 // A default inilized AliITSInitGeometry object
112 if(version == kvPPRasymmFMD && (fMinorVersion==1|| fMinorVersion==2)){
113 fName="AliITSvPPRasymmFMD";
114 }else if(version == kv11Hybrid){
115 fName="AliITSv11Hybrid";
117 AliFatal(Form("Undefined geometry: fMajorVersion=%d, "
118 "fMinorVersion= %d",(Int_t)fMajorVersion,fMinorVersion));
123 //______________________________________________________________________
124 AliITSgeom* AliITSInitGeometry::CreateAliITSgeom(){
125 // Creates and Initilizes the geometry transformation class AliITSgeom
126 // to values appropreate to this specific geometry. Now that
127 // the segmentation is part of AliITSgeom, the detector
128 // segmentations are also defined here.
134 // A pointer to a new properly inilized AliITSgeom class. If
135 // pointer = 0 then failed to init.
138 AliITSVersion_t version = kvDefault;
141 TGeoVolume *itsV = gGeoManager->GetVolume("ITSV");
143 Error("CreateAliITSgeom","Can't find ITS volume ITSV, aborting");
146 const Char_t *title = itsV->GetTitle();
147 if(!ReadVersionString(title,(Int_t)strlen(title),version,minor,
149 Warning("UpdateInternalGeometry","Can't read title=%s\n",title);
153 AliITSgeom *geom = CreateAliITSgeom(version,minor);
154 AliDebug(1,"AliITSgeom object has been initialized from TGeo\n");
157 //______________________________________________________________________
158 AliITSgeom* AliITSInitGeometry::CreateAliITSgeom(Int_t major,Int_t minor){
159 // Creates and Initilizes the geometry transformation class AliITSgeom
160 // to values appropreate to this specific geometry. Now that
161 // the segmentation is part of AliITSgeom, the detector
162 // segmentations are also defined here.
164 // Int_t major major version, see AliITSVersion_t
165 // Int_t minor minor version
169 // A pointer to a new properly inilized AliITSgeom class. If
170 // pointer = 0 then failed to init.
174 SetGeometryName("AliITSvtest");
175 SetVersion(kvtest,minor);
178 SetGeometryName("AliITSvSPD02");
179 SetVersion(kvSPD02,minor);
182 SetGeometryName("AliITSvSDD03");
183 SetVersion(kvSDD03,minor);
186 SetGeometryName("AliITSvSSD03");
187 SetVersion(kvSSD03,minor);
190 SetGeometryName("AliITSvBeamTest03");
191 SetVersion(kvITS04,minor);
193 case kvPPRcourseasymm:
194 SetGeometryName("AliITSvPPRcourseasymm");
195 SetVersion(kvPPRcourseasymm,minor);
198 SetGeometryName("AliITSvPPRasymmFMD");
199 SetVersion(kvPPRasymmFMD,minor);
202 SetGeometryName("AliITSv11");
203 SetVersion(kv11,minor);
206 SetGeometryName("AliITSv11Hybrid");
207 SetVersion(kv11Hybrid,minor);
211 SetGeometryName("Undefined");
212 SetVersion(kvDefault,minor);
215 AliITSgeom *geom = new AliITSgeom();
216 if(!InitAliITSgeom(geom)){ // Error initilization failed
222 //______________________________________________________________________
223 Bool_t AliITSInitGeometry::InitAliITSgeom(AliITSgeom *geom){
224 // Initilizes the geometry transformation class AliITSgeom
225 // to values appropreate to this specific geometry. Now that
226 // the segmentation is part of AliITSgeom, the detector
227 // segmentations are also defined here.
229 // AliITSgeom *geom A pointer to the AliITSgeom class
231 // AliITSgeom *geom This pointer recreated and properly inilized.
236 AliFatal("The geometry manager has not been initialized (e.g. "
237 "TGeoManager::Import(\"geometry.root\")should be "
238 "called in advance) - exit forced");
241 switch(fMajorVersion) {
243 if(GetMinorVersion()==1) return InitAliITSgeomPPRasymmFMD(geom);
244 else if(GetMinorVersion()==2) return InitAliITSgeomtest2(geom);
247 return InitAliITSgeomSPD02(geom);
250 return InitAliITSgeomSDD03(geom);
253 return InitAliITSgeomSSD03(geom);
256 return InitAliITSgeomITS04(geom);
258 case kvPPRasymmFMD: {
259 return InitAliITSgeomPPRasymmFMD(geom);
261 case kvPPRcourseasymm: {
262 return kTRUE; // No sensitive detectors in course geometry
265 return InitAliITSgeomV11Hybrid(geom);
268 return InitAliITSgeomV11(geom);
270 case kvDefault: default: {
271 AliFatal("Undefined geometry");
277 //______________________________________________________________________
278 void AliITSInitGeometry::TransposeTGeoHMatrix(TGeoHMatrix *m)const{
279 // Transpose the rotation matrix part of a TGeoHMatrix. This
280 // is needed because TGeo stores the transpose of the rotation
281 // matrix as compared to what AliITSgeomMatrix uses (and Geant3).
283 // TGeoHMatrix *m The matrix to be transposed
285 // TGEoHMatrix *m The transposed matrix
291 if(m==0) return; // no matrix to transpose.
292 for(i=0;i<9;i += 4) r[i] = m->GetRotationMatrix()[i]; // diagonals
293 r[1] = m->GetRotationMatrix()[3];
294 r[2] = m->GetRotationMatrix()[6];
295 r[3] = m->GetRotationMatrix()[1];
296 r[5] = m->GetRotationMatrix()[7];
297 r[6] = m->GetRotationMatrix()[2];
298 r[7] = m->GetRotationMatrix()[5];
302 //______________________________________________________________________
303 Bool_t AliITSInitGeometry::InitAliITSgeomtest2(AliITSgeom *geom){
304 // Initilizes the geometry transformation class AliITSgeom
305 // to values appropreate to this specific geometry. Now that
306 // the segmentation is part of AliITSgeom, the detector
307 // segmentations are also defined here.
309 // AliITSgeom *geom A pointer to the AliITSgeom class
311 // AliITSgeom *geom This pointer recreated and properly inilized.
314 // const Double_t kcm2micron = 1.0E4;
315 const Int_t kItype=0; // Type of transormation defined 0=> Geant
316 const Int_t klayers = 6; // number of layers in the ITS
317 const Int_t kladders[klayers] = {1,1,1,1,1,1}; // Number of ladders
318 const Int_t kdetectors[klayers] = {1,1,1,1,1,1};// number of detector/lad
319 const AliITSDetector kIdet[6] = {kSPD,kSPD,kSDD,kSDD,kSSD,kSSD};
320 const TString kNames[klayers] = {
321 "/ALIC_1/ITSV_1/ITSspd1_1/ITS1_1", // lay=1
322 "/ALIC_1/ITSV_1/ITSspd2_1/ITS2_1", // lay=2
323 "/ALIC_1/ITSV_1/ITSsdd1_1/ITS3_1", // lay=3
324 "/ALIC_1/ITSV_1/ITSsdd2_1/ITS4_1", // lay=4
325 "/ALIC_1/ITSV_1/ITSssd1_1/ITS5_1", // lay=5
326 "/ALIC_1/ITSV_1/ITSssd2_1/ITS6_1"};// Lay=6
327 Int_t mod,nmods=0,lay,lad,det,cpn0,cpn1,cpn2;
328 Double_t tran[3]={0.0,0.0,0.0},rot[10]={9*0.0,1.0};
332 Bool_t initSeg[3]={kFALSE,kFALSE,kFALSE};
333 TStopwatch *time = 0x0;if(fTiming) time=new TStopwatch();
335 if(fTiming) time->Start();
336 for(mod=0;mod<klayers;mod++) nmods += kladders[mod]*kdetectors[mod];
337 geom->Init(kItype,klayers,kladders,kdetectors,nmods);
338 for(mod=0;mod<nmods;mod++){
339 DecodeDetectorLayers(mod,lay,lad,det); // Write
340 geom->CreateMatrix(mod,lay,lad,det,kIdet[lay-1],tran,rot);
341 RecodeDetector(mod,cpn0,cpn1,cpn2); // Write reusing lay,lad,det.
342 geom->GetGeomMatrix(mod)->SetPath(kNames[lay-1]);
343 GetTransformation(kNames[lay-1].Data(),matrix);
344 geom->SetTrans(mod,matrix.GetTranslation());
345 TransposeTGeoHMatrix(&matrix); // Transpose TGeo's rotation matrixes
346 geom->SetRotMatrix(mod,matrix.GetRotationMatrix());
347 if(initSeg[kIdet[lay-1]]) continue;
348 GetShape(kNames[lay-1],shapeName,shapePar);
349 if(shapeName.CompareTo("BOX")){
350 Error("InitITSgeom2","Geometry changed without proper code update"
351 "or error in reading geometry. Shape is not BOX shape is %s",
355 InitGeomShapePPRasymmFMD(kIdet[lay-1],initSeg,shapePar,geom);
364 //______________________________________________________________________
365 Bool_t AliITSInitGeometry::InitAliITSgeomSPD02(AliITSgeom *geom){
366 // Initilizes the geometry transformation class AliITSgeom
367 // to values appropreate to this specific geometry. Now that
368 // the segmentation is part of AliITSgeom, the detector
369 // segmentations are also defined here.
371 // AliITSgeom *geom A pointer to the AliITSgeom class
373 // AliITSgeom *geom This pointer recreated and properly inilized.
376 const Int_t kltypess=2;
377 const Int_t knlayers=5;
378 const TString knames[kltypess]=
379 {"ALIC_1/ITSV_1/ITEL_%d/IMB0_1/IMBS_1",//lay=1,2,4,5
380 "ALIC_1/ITSV_1/IDET_%d/ITS0_1/ITST_1"};// lay=3
381 const Int_t kitsGeomTreeCopys[2]={4,1};
382 const Int_t knlad[knlayers]={knlayers*1},kndet[knlayers]={knlayers*1};
383 TString path,shapeName;
387 Double_t trans[3]={3*0.0},rot[10]={10*0.0};
388 Int_t npar=3,mod,i,j,lay,lad,det,cpy;
390 TStopwatch *time = 0x0;if(fTiming) time=new TStopwatch();
392 par[0]=0.64;par[1]=0.5*300.0E-4;par[2]=3.48;
394 geom->Init(0,knlayers,knlad,kndet,mod);
396 if(fTiming) time->Start();
397 for(i=0;i<kltypess;i++)for(cpy=1;cpy<=kitsGeomTreeCopys[i];cpy++){
398 path.Form(knames[i].Data(),cpy);
399 GetTransformation(path.Data(),matrix);
400 GetShape(path.Data(),shapeName,shapePar);
401 shapeParF.Set(shapePar.GetSize());
402 for(j=0;j<shapePar.GetSize();j++) shapeParF[j]=shapePar[j];
404 if(i==0&&cpy>2) lay=cpy+1;
406 DecodeDetector(mod,kitsGeomTreeCopys[i],1,cpy,0);
407 DecodeDetectorLayers(mod,lay,lad,det);
408 geom->CreateMatrix(mod,lay,lad,det,kSPD,trans,rot);
409 geom->SetTrans(mod,matrix.GetTranslation());
410 geom->SetRotMatrix(mod,matrix.GetRotationMatrix());
411 geom->GetGeomMatrix(mod)->SetPath(path.Data());
412 if(!(geom->IsShapeDefined((Int_t)kSPD)))
413 geom->ReSetShape(kSPD,new AliITSgeomSPD425Short(npar,par));
422 //______________________________________________________________________
423 Bool_t AliITSInitGeometry::InitAliITSgeomSDD03(AliITSgeom *geom){
424 // Initilizes the geometry transformation class AliITSgeom
425 // to values appropreate to this specific geometry. Now that
426 // the segmentation is part of AliITSgeom, the detector
427 // segmentations are also defined here.
429 // AliITSgeom *geom A pointer to the AliITSgeom class
431 // AliITSgeom *geom This pointer recreated and properly inilized.
434 const Int_t knlayers=12;
435 // const Int_t kndeep=6;
436 const Int_t kltypess=2;
437 const AliITSDetector kidet[knlayers]={kSSD,kSDD};
438 const TString knames[kltypess]={
439 "/ALIC_1/ITSV_1/ITEL_%d/ITAI_1/IMB0_1/IMBS_1",
440 "/ALIC_1/ITSV_1/IDET_%d/IDAI_1/ITS0_1/ITST_1"};
441 const Int_t kitsGeomTreeCopys[kltypess]={10,2};
443 const Float_t kpitch=50.E-4;/*cm*/
444 Float_t box[3]={0.5*kpitch*(Float_t)knp,150.E-4,1.0},p[knp+1],n[knp+1];
445 Int_t nlad[knlayers]={knlayers*1};
446 Int_t ndet[knlayers]={knlayers*1};
447 Int_t mod=knlayers,lay=0,lad=0,det=0,i,j,cp0;
448 TString path,shapeName;
450 Double_t trans[3]={3*0.0},rot[10]={10*0.0};
453 Bool_t isShapeDefined[kltypess]={kltypess*kFALSE};
455 geom->Init(0,knlayers,nlad,ndet,mod);
458 // Fill in anode and cathode strip locations (lower edge)
465 for(i=0;i<kltypess;i++)for(cp0=1;cp0<=kitsGeomTreeCopys[i];cp0++){
466 DecodeDetector(mod,kitsGeomTreeCopys[i],cp0,1,2);
467 DecodeDetectorLayers(mod,lay,lad,det);
468 path.Form(knames[i].Data(),cp0);
469 GetTransformation(path.Data(),matrix);
470 GetShape(path.Data(),shapeName,shapePar);
471 shapeParF.Set(shapePar.GetSize());
472 for(j=0;j<shapePar.GetSize();j++)shapeParF[j]=shapePar[j];
473 geom->CreateMatrix(mod,lay,lad,det,kidet[i],trans,rot);
474 geom->SetTrans(mod,matrix.GetTranslation());
475 geom->SetRotMatrix(mod,matrix.GetRotationMatrix());
476 geom->GetGeomMatrix(mod)->SetPath(path.Data());
478 case kSDD: if(!(geom->IsShapeDefined((Int_t)kSDD))){
479 geom->ReSetShape(kSDD,new AliITSgeomSDD256(shapeParF.GetSize(),
480 shapeParF.GetArray()));
481 isShapeDefined[i]=kTRUE;
483 case kSSD:if(!(geom->IsShapeDefined((Int_t)kSSD))){
484 geom->ReSetShape(kSSD,new AliITSgeomSSD(box,0.0,0.0,
486 isShapeDefined[i]=kTRUE;
494 //______________________________________________________________________
495 Bool_t AliITSInitGeometry::InitAliITSgeomSSD03(AliITSgeom *geom){
496 // Initilizes the geometry transformation class AliITSgeom
497 // to values appropreate to this specific geometry. Now that
498 // the segmentation is part of AliITSgeom, the detector
499 // segmentations are also defined here.
501 // AliITSgeom *geom A pointer to the AliITSgeom class
503 // AliITSgeom *geom This pointer recreated and properly inilized.
506 const Int_t knlayers=5;
507 // const Int_t kndeep=6;
508 const Int_t kltypess=3;
509 const AliITSDetector kIdet[knlayers]={kND,kSSD,kND};
510 const TString knames[kltypess]={
511 "/ALIC_1/ITSV_1/ITSA_%d/ITSS_1",
512 "/ALIC_1/ITSV_1/IGAR_%d/IAIR_1/ITST_1",
513 "/ALIC_1/ITSV_1/IFRA_%d/IFRS_1"};
514 const Int_t kitsGeomTreeCopys[kltypess]={3,1,1};
515 const Int_t kitsGeomDetTypes[kltypess]={1,2,3};
517 const Float_t kpitch=50.E-4;//cm
518 Bool_t initSeg[3]={kFALSE, kFALSE, kFALSE};
519 Float_t box[3]={0.5*kpitch*(Float_t)knp,150.E-4,1.0},p[knp+1],n[knp+1];
520 Int_t nlad[knlayers]={knlayers*1};
521 Int_t ndet[knlayers]={knlayers*1};
522 Int_t mod=knlayers,lay=0,lad=0,det=0,i,j,cp0;
523 TString path,shapeName;
525 Double_t trans[3]={3*0.0},rot[10]={10*0.0};
528 Bool_t isShapeDefined[kltypess]={kltypess*kFALSE};
530 geom->Init(0,knlayers,nlad,ndet,mod);
533 // Fill in anode and cathode strip locations (lower edge)
540 for(i=0;i<kltypess;i++)for(cp0=1;cp0<=kitsGeomTreeCopys[i];cp0++){
541 DecodeDetector(mod,kitsGeomDetTypes[i],cp0,1,0);
542 DecodeDetectorLayers(mod,lay,lad,det);
543 path.Form(knames[i].Data(),cp0);
544 GetTransformation(path.Data(),matrix);
545 GetShape(path.Data(),shapeName,shapePar);
546 shapeParF.Set(shapePar.GetSize());
547 for(j=0;j<shapePar.GetSize();j++)shapeParF[j]=shapePar[j];
548 geom->CreateMatrix(mod,lay,lad,det,kIdet[i],trans,rot);
549 geom->SetTrans(mod,matrix.GetTranslation());
550 geom->SetRotMatrix(mod,matrix.GetRotationMatrix());
551 geom->GetGeomMatrix(mod)->SetPath(path.Data());
553 case kSSD:if(!(geom->IsShapeDefined((Int_t)kSSD))){
554 InitGeomShapePPRasymmFMD(kIdet[lay-1],initSeg,shapePar,geom);
555 isShapeDefined[i]=kTRUE;
563 //______________________________________________________________________
564 Bool_t AliITSInitGeometry::InitAliITSgeomITS04(AliITSgeom *geom){
565 // Initilizes the geometry transformation class AliITSgeom
566 // to values appropreate to this specific geometry. Now that
567 // the segmentation is part of AliITSgeom, the detector
568 // segmentations are also defined here.
570 // AliITSgeom *geom A pointer to the AliITSgeom class
572 // AliITSgeom *geom This pointer recreated and properly inilized.
576 // We can not use AliITSvBeamTestITS04::fgk... data members because
577 // AliITSInitGeometry is part of the base library while AliITSvBeamTestITS04
578 // is part of the simulation library. This would introduce a dependance
579 // between the 2 libraries
582 const Int_t knlayers = 6;
583 Int_t nlad[knlayers], ndet[knlayers];
585 nlad[0] = 1; ndet[0] = 2;
586 nlad[1] = 1; ndet[1] = 2;
587 nlad[2] = 1; ndet[2] = 1;
588 nlad[3] = 1; ndet[3] = 1;
589 nlad[4] = 1; ndet[4] = 2;
590 nlad[5] = 1; ndet[5] = 2;
593 geom->Init(0,knlayers,nlad,ndet,nModTot);
596 //=== Set default shapes
597 const Float_t kDxyzSPD[] = {AliITSvBeamTestITS04::fgkSPDwidthSens/2,
598 AliITSvBeamTestITS04::fgkSPDthickSens/2,
599 AliITSvBeamTestITS04::fgkSPDlengthSens/2};
600 if(!(geom->IsShapeDefined(kSPD)))
601 geom->ReSetShape(kSPD,new AliITSgeomSPD425Short(3,(Float_t *)kDxyzSPD));
603 const Float_t kDxyzSDD[] = {AliITSvBeamTestITS04::fgkSDDwidthSens/2.,
604 AliITSvBeamTestITS04::fgkSDDthickSens/2.,
605 AliITSvBeamTestITS04::fgkSDDlengthSens/2.};
606 if(!(geom->IsShapeDefined(kSDD)))
607 geom->ReSetShape(kSDD, new AliITSgeomSDD256(3,(Float_t *)kDxyzSDD));
609 const Float_t kDxyzSSD[] = {AliITSvBeamTestITS04::fgkSSDlengthSens/2,
610 AliITSvBeamTestITS04::fgkSSDthickSens/2,
611 AliITSvBeamTestITS04::fgkSSDwidthSens/2};
612 if(!(geom->IsShapeDefined(kSSD)))
613 geom->ReSetShape(kSSD,new AliITSgeomSSD75and275(3,(Float_t *)kDxyzSSD));
615 // Creating the matrices in AliITSgeom for each sensitive volume
616 // (like in AliITSv11GeometrySDD) mln
617 // Here, each layer is one detector
620 Int_t startMod = 0,mod;
621 TGeoVolume *itsmotherVolume = gGeoManager->GetVolume("ITSV");
623 for (Int_t i=0; i<4;i++) {
624 sprintf(layerName, "ITSspdWafer_%i",i+1);
625 TGeoNode *layNode = itsmotherVolume->GetNode(layerName);
627 TGeoHMatrix layMatrix(*layNode->GetMatrix());
628 Double_t *trans = layMatrix.GetTranslation();
629 Double_t *r = layMatrix.GetRotationMatrix();
630 Double_t rot[10] = {r[0],r[1],r[2],
632 r[6],r[7],r[8], 1.0};
636 DecodeDetector(mod,layNode->GetNumber(),i+1,0,0);
637 DecodeDetectorLayers(mod,iLay,iLad,iDet);
638 geom->CreateMatrix(startMod,iLay,iLad,iDet,kSPD,trans,rot);
644 for (Int_t i=0; i<2;i++) {
645 sprintf(layerName, "ITSsddWafer_%i",i+4+1);
646 TGeoNode *layNode = itsmotherVolume->GetNode(layerName);
648 TGeoHMatrix layMatrix(*layNode->GetMatrix());
649 Double_t *trans = layMatrix.GetTranslation();
650 Double_t *r = layMatrix.GetRotationMatrix();
651 Double_t rot[10] = {r[0],r[1],r[2],
653 r[6],r[7],r[8], 1.0};
657 DecodeDetector(mod,layNode->GetNumber(),i+1,0,0);
658 DecodeDetectorLayers(mod,iLay,iLad,iDet);
659 geom->CreateMatrix(startMod,iLay,iLad,iDet,kSDD,trans,rot);
665 for (Int_t i=0; i<4;i++) {
666 sprintf(layerName, "ITSssdWafer_%i",i+4+2+1);
667 TGeoNode *layNode = itsmotherVolume->GetNode(layerName);
669 TGeoHMatrix layMatrix(*layNode->GetMatrix());
670 Double_t *trans = layMatrix.GetTranslation();
671 Double_t *r = layMatrix.GetRotationMatrix();
672 Double_t rot[10] = {r[0],r[1],r[2],
674 r[6],r[7],r[8], 1.0};
678 DecodeDetector(mod,layNode->GetNumber(),i+1,0,0);
679 DecodeDetectorLayers(mod,iLay,iLad,iDet);
680 geom->CreateMatrix(startMod,iLay,iLad,iDet,kSSD,trans,rot);
689 //______________________________________________________________________
690 Bool_t AliITSInitGeometry::InitAliITSgeomPPRasymmFMD(AliITSgeom *geom){
691 // Initilizes the geometry transformation class AliITSgeom
692 // to values appropreate to this specific geometry. Now that
693 // the segmentation is part of AliITSgeom, the detector
694 // segmentations are also defined here.
696 // AliITSgeom *geom A pointer to the AliITSgeom class
698 // AliITSgeom *geom This pointer recreated and properly inilized.
701 // const Double_t kcm2micron = 1.0E4;
702 const Int_t kItype=0; // Type of transormation defined 0=> Geant
703 const Int_t klayers = 6; // number of layers in the ITS
704 const Int_t kladders[klayers] = {20,40,14,22,34,38}; // Number of ladders
705 const Int_t kdetectors[klayers] = {4,4,6,8,22,25};// number of detector/lad
706 const AliITSDetector kIdet[6] = {kSPD,kSPD,kSDD,kSDD,kSSD,kSSD};
707 const TString kPathbase = "/ALIC_1/ITSV_1/ITSD_1/";
708 const TString kNames[2][klayers] = {
709 {"%sIT12_1/I12A_%d/I10A_%d/I103_%d/I101_1/ITS1_1", // lay=1
710 "%sIT12_1/I12A_%d/I20A_%d/I1D3_%d/I1D1_1/ITS2_1", // lay=2
711 "%sIT34_1/I004_%d/I302_%d/ITS3_%d/", // lay=3
712 "%sIT34_1/I005_%d/I402_%d/ITS4_%d/", // lay=4
713 "%sIT56_1/I565_%d/I562_%d/ITS5_%d/", // lay=5
714 "%sIT56_1/I569_%d/I566_%d/ITS6_%d/"},// lay=6
715 {"%sIT12_1/I12B_%d/I10B_%d/I107_%d/I101_1/ITS1_1", // lay=1
716 "%sIT12_1/I12B_%d/I20B_%d/I1D7_%d/I1D1_1/ITS2_1", // lay=2
717 "%sIT34_1/I004_%d/I302_%d/ITS3_%d", // lay=3
718 "%sIT34_1/I005_%d/I402_%d/ITS4_%d", // lay=4
719 "%sIT56_1/I565_%d/I562_%d/ITS5_%d", // lay=5
720 "%sIT56_1/I569_%d/I566_%d/ITS6_%d"}};// Lay=6
722 Int_t itsGeomTreeCopys[knlayers][3]= {{10, 2, 4},// lay=1
729 Int_t mod,nmods=0,lay,lad,det,cpn0,cpn1,cpn2;
730 Double_t tran[3]={0.0,0.0,0.0},rot[10]={9*0.0,1.0};
732 TString path,shapeName;
734 Bool_t initSeg[3]={kFALSE,kFALSE,kFALSE};
735 TStopwatch *time = 0x0;if(fTiming) time=new TStopwatch();
737 if(fTiming) time->Start();
738 for(mod=0;mod<klayers;mod++) nmods += kladders[mod]*kdetectors[mod];
739 geom->Init(kItype,klayers,kladders,kdetectors,nmods);
740 for(mod=0;mod<nmods;mod++){
741 DecodeDetectorLayers(mod,lay,lad,det); // Write
742 geom->CreateMatrix(mod,lay,lad,det,kIdet[lay-1],tran,rot);
743 RecodeDetector(mod,cpn0,cpn1,cpn2); // Write reusing lay,lad,det.
744 path.Form(kNames[fMinorVersion-1][lay-1].Data(),
745 kPathbase.Data(),cpn0,cpn1,cpn2);
746 geom->GetGeomMatrix(mod)->SetPath(path);
747 GetTransformation(path.Data(),matrix);
748 geom->SetTrans(mod,matrix.GetTranslation());
749 TransposeTGeoHMatrix(&matrix); //Transpose TGeo's rotation matrixes
750 geom->SetRotMatrix(mod,matrix.GetRotationMatrix());
751 if(initSeg[kIdet[lay-1]]) continue;
752 GetShape(path,shapeName,shapePar);
753 if(shapeName.CompareTo("BOX")){
754 Error("InitITSgeomPPRasymmFMD",
755 "Geometry changed without proper code update or error "
756 "in reading geometry. Shape is not BOX. Shape is %s",
760 InitGeomShapePPRasymmFMD(kIdet[lay-1],initSeg,shapePar,geom);
769 //______________________________________________________________________
770 Bool_t AliITSInitGeometry::InitAliITSgeomV11Hybrid(AliITSgeom *geom){
771 // Initilizes the geometry transformation class AliITSgeom
772 // to values appropreate to this specific geometry. Now that
773 // the segmentation is part of AliITSgeom, the detector
774 // segmentations are also defined here.
776 // AliITSgeom *geom A pointer to the AliITSgeom class
778 // AliITSgeom *geom This pointer recreated and properly inilized.
782 const Int_t kItype = 0; // Type of transformation defined 0=> Geant
783 const Int_t klayers = 6; // number of layers in the ITS
784 const Int_t kladders[klayers] = {20,40,14,22,34,38}; // Number of ladders
785 const Int_t kdetectors[klayers] = {4,4,6,8,22,25};// number of detector/lad
786 const AliITSDetector kIdet[6] = {kSPD,kSPD,kSDD,kSDD,kSSD,kSSD};
787 const TString kPathbase = "/ALIC_1/ITSV_1/";
789 char *pathSDDsens1, *pathSDDsens2;
790 if (SDDIsTGeoNative()) {
791 pathSDDsens1 = "%sITSsddLayer3_1/ITSsddLadd_%d/ITSsddSensor3_%d/ITSsddWafer3_%d/ITSsddSensitivL3_1";
792 pathSDDsens2 = "%sITSsddLayer4_1/ITSsddLadd_%d/ITSsddSensor4_%d/ITSsddWafer4_%d/ITSsddSensitivL4_1";
794 pathSDDsens1 = "%sITSD_1/IT34_1/I004_%d/I302_%d/ITS3_%d";
795 pathSDDsens2 = "%sITSD_1/IT34_1/I005_%d/I402_%d/ITS4_%d";
798 char *pathSSDsens1, *pathSSDsens2;
799 if (SSDIsTGeoNative()) {
800 pathSSDsens1 = "%sITSssdLayer5_1/ITSssdLay5Ladd_%d/ITSsddSensor5_%d/ITSsddSensitivL5_1";
801 pathSSDsens2 = "%sITSssdLayer6_1/ITSssdLay6Ladd_%d/ITSsddSensor6_%d/ITSsddSensitivL6_1";
803 pathSSDsens1 = "%sITSD_1/IT56_1/I565_%d/I562_%d/ITS5_%d";
804 pathSSDsens2 = "%sITSD_1/IT56_1/I569_%d/I566_%d/ITS6_%d";
807 const TString kNames[klayers] = {
808 "%sITSD_1/IT12_1/I12B_%d/I10B_%d/I107_%d/I101_1/ITS1_1", // lay=1
809 "%sITSD_1/IT12_1/I12B_%d/I20B_%d/I1D7_%d/I1D1_1/ITS2_1", // lay=2
810 pathSDDsens1, // lay=3
811 pathSDDsens2, // lay=4
812 pathSSDsens1, // lay=5
813 pathSSDsens2};// Lay=6
815 Int_t mod,nmods=0, lay, lad, det, cpn0, cpn1, cpn2;
816 Double_t tran[3]={0.,0.,0.}, rot[10]={9*0.0,1.0};
818 TString path, shapeName;
820 Bool_t initSeg[3]={kFALSE, kFALSE, kFALSE};
821 TStopwatch *time = 0x0;
822 if(fTiming) time = new TStopwatch();
824 if(fTiming) time->Start();
825 for(mod=0;mod<klayers;mod++) nmods += kladders[mod]*kdetectors[mod];
826 geom->Init(kItype,klayers,kladders,kdetectors,nmods);
828 for(mod=0;mod<nmods;mod++) {
830 // DecodeDetectorLayers(mod,lay,lad,det); // Write
831 // geom->CreateMatrix(mod,lay,lad,det,kIdet[lay-1],tran,rot);
832 // RecodeDetector(mod,cpn0,cpn1,cpn2); // Write reusing lay,lad,det.
833 DecodeDetectorLayers(mod,lay,lad,det); // Write
834 geom->CreateMatrix(mod,lay,lad,det,kIdet[lay-1],tran,rot);
835 RecodeDetector(mod,cpn0,cpn1,cpn2); // Write reusing lay,lad,det.
837 if (SDDIsTGeoNative())
838 if (kIdet[lay-1]==kSDD) {
843 if (SSDIsTGeoNative())
844 if (kIdet[lay-1]==kSSD) {
850 path.Form(kNames[lay-1].Data(),
851 kPathbase.Data(),cpn0,cpn1,cpn2);
852 geom->GetGeomMatrix(mod)->SetPath(path);
853 GetTransformation(path.Data(),matrix);
854 geom->SetTrans(mod,matrix.GetTranslation());
855 TransposeTGeoHMatrix(&matrix); //Transpose TGeo's rotation matrixes
856 geom->SetRotMatrix(mod,matrix.GetRotationMatrix());
857 if(initSeg[kIdet[lay-1]]) continue;
858 GetShape(path,shapeName,shapePar);
859 if(shapeName.CompareTo("BOX")){
860 Error("InitITSgeom","Geometry changed without proper code update"
861 "or error in reading geometry. Shape is not BOX.");
864 InitGeomShapePPRasymmFMD(kIdet[lay-1],initSeg,shapePar,geom);
874 //______________________________________________________________________
875 Bool_t AliITSInitGeometry::InitAliITSgeomV11(AliITSgeom *geom){
876 // Initilizes the geometry transformation class AliITSgeom
877 // Now that the segmentation is part of AliITSgeom, the detector
878 // segmentations are also defined here.
881 // AliITSgeom *geom A pointer to the AliITSgeom class
883 // AliITSgeom *geom This pointer recreated and properly inilized.
887 const Int_t kItype=0; // Type of transormation defined 0=> Geant
888 const Int_t klayers = 6; // number of layers in the ITS
889 const Int_t kladders[klayers] = {20,40,14,22,34,38}; // Number of ladders
890 const Int_t kdetectors[klayers] = {4,4,6,8,22,25};// number of detector/lad
891 const AliITSDetector kIdet[6] = {kSPD,kSPD,kSDD,kSDD,kSSD,kSSD};
893 const TString kPathbase = "/ALIC_1/ITSV_1/";
894 const TString kNames[klayers] =
895 {"AliITSInitGeometry:spd missing", // lay=1
896 "AliITSInitGeometry:spd missing", // lay=2
897 "%sITSsddLayer3_1/ITSsddLadd_%d/ITSsddSensor_%d/ITSsddWafer_1/ITSsddSensitiv_1", // lay=3
898 "%sITSsddLayer4_1/ITSsddLadd_%d/ITSsddSensor_%d/ITSsddWafer_1/ITSsddSensitiv_1", // lay=4
899 "AliITSInitGeometry:ssd missing", // lay=5
900 "AliITSInitGeometry:ssd missing"};// lay=6
902 Int_t mod,nmods=0,lay,lad,det,cpn0,cpn1,cpn2;
903 Double_t tran[3]={0.0,0.0,0.0},rot[10]={9*0.0,1.0};
905 TString path,shapeName;
907 Bool_t initSeg[3]={kFALSE,kFALSE,kFALSE};
908 TStopwatch *time = 0x0;if(fTiming) time=new TStopwatch();
910 if(fTiming) time->Start();
911 for(mod=0;mod<klayers;mod++) nmods += kladders[mod]*kdetectors[mod];
913 geom->Init(kItype,klayers,kladders,kdetectors,nmods);
914 for(mod=0;mod<nmods;mod++) {
916 DecodeDetectorLayers(mod,lay,lad,det); // Write
917 geom->CreateMatrix(mod,lay,lad,det,kIdet[lay-1],tran,rot);
918 RecodeDetector(mod,cpn0,cpn1,cpn2); // Write reusing lay,lad,det.
919 path.Form(kNames[lay-1].Data(),
920 kPathbase.Data(),cpn0,cpn1,cpn2);
921 geom->GetGeomMatrix(mod)->SetPath(path);
922 if (GetTransformation(path.Data(),matrix)) {
923 geom->SetTrans(mod,matrix.GetTranslation());
924 TransposeTGeoHMatrix(&matrix); //Transpose TGeo's rotation matrixes
925 geom->SetRotMatrix(mod,matrix.GetRotationMatrix());
928 if(initSeg[kIdet[lay-1]]) continue;
929 GetShape(path,shapeName,shapePar);
930 if(shapeName.CompareTo("BOX")){
931 Error("InitAliITSgeomV11","Geometry changed without proper code update"
932 "or error in reading geometry. Shape is not BOX.");
935 InitGeomShapePPRasymmFMD(kIdet[lay-1],initSeg,shapePar,geom);
947 //______________________________________________________________________
948 Bool_t AliITSInitGeometry::InitGeomShapePPRasymmFMD(AliITSDetector idet,
952 // Initilizes the geometry segmentation class AliITSgeomS?D, or
953 // AliITSsegmentationS?D depending on the vaule of fSegGeom,
954 // to values appropreate to this specific geometry. Now that
955 // the segmentation is part of AliITSgeom, the detector
956 // segmentations are also defined here.
958 // Int_t lay The layer number/name.
959 // AliITSgeom *geom A pointer to the AliITSgeom class
961 // AliITSgeom *geom This pointer recreated and properly inilized.
964 // const Double_t kcm2micron = 1.0E4;
965 const Double_t kmicron2cm = 1.0E-4;
969 shapeParF.Set(shapePar.GetSize());
970 for(i=0;i<shapePar.GetSize();i++) shapeParF[i]=shapePar[i];
973 initSeg[idet] = kTRUE;
974 AliITSgeomSPD *geomSPD = new AliITSgeomSPD425Short();
975 Float_t bx[256],bz[280];
976 for(i=000;i<256;i++) bx[i] = 50.0*kmicron2cm; // in x all are 50 microns.
977 for(i=000;i<160;i++) bz[i] = 425.0*kmicron2cm; // most are 425 microns
979 for(i=160;i<280;i++) bz[i] = 0.0*kmicron2cm; // Outside of detector.
980 bz[ 31] = bz[ 32] = 625.0*kmicron2cm; // first chip boundry
981 bz[ 63] = bz[ 64] = 625.0*kmicron2cm; // first chip boundry
982 bz[ 95] = bz[ 96] = 625.0*kmicron2cm; // first chip boundry
983 bz[127] = bz[128] = 625.0*kmicron2cm; // first chip boundry
984 bz[160] = 425.0*kmicron2cm;// Set so that there is no zero pixel size for fNz.
985 geomSPD->ReSetBins(shapeParF[1],256,bx,160,bz);
986 geom->ReSetShape(idet,geomSPD);
989 initSeg[idet] = kTRUE;
990 AliITSgeomSDD *geomSDD = new AliITSgeomSDD256(shapeParF.GetSize(),
991 shapeParF.GetArray());
992 geom->ReSetShape(idet,geomSDD);
995 initSeg[idet] = kTRUE;
996 AliITSgeomSSD *geomSSD = new AliITSgeomSSD275and75(
997 shapeParF.GetSize(),shapeParF.GetArray());
998 geom->ReSetShape(idet,geomSSD);
1000 default:{// Others, Note no kSDDp or kSSDp in this geometry.
1001 geom->ReSetShape(idet,0);
1002 Info("InitGeomShapePPRasymmFMD",
1003 "default Dx=%f Dy=%f Dz=%f default=%d",
1004 shapePar[0],shapePar[1],shapePar[2],idet);
1009 //______________________________________________________________________
1010 Bool_t AliITSInitGeometry::InitSegmentationPPRasymmFMD(AliITSDetector idet,
1014 // Initilizes the geometry segmentation class AliITSgeomS?D, or
1015 // AliITSsegmentationS?D depending on the vaule of fSegGeom,
1016 // to values appropreate to this specific geometry. Now that
1017 // the segmentation is part of AliITSgeom, the detector
1018 // segmentations are also defined here.
1020 // Int_t lay The layer number/name.
1021 // AliITSgeom *geom A pointer to the AliITSgeom class
1023 // AliITSgeom *geom This pointer recreated and properly inilized.
1026 const Double_t kcm2micron = 1.0E4;
1031 initSeg[idet] = kTRUE;
1032 AliITSsegmentationSPD *segSPD = new AliITSsegmentationSPD();
1033 segSPD->SetDetSize(2.*shapePar[0]*kcm2micron, // X
1034 2.*shapePar[2]*kcm2micron, // Z
1035 2.*shapePar[1]*kcm2micron);// Y Microns
1036 segSPD->SetNPads(256,160);// Number of Bins in x and z
1037 Float_t bx[256],bz[280];
1038 for(i=000;i<256;i++) bx[i] = 50.0; // in x all are 50 microns.
1039 for(i=000;i<160;i++) bz[i] = 425.0; // most are 425 microns
1041 for(i=160;i<280;i++) bz[i] = 0.0; // Outside of detector.
1042 bz[ 31] = bz[ 32] = 625.0; // first chip boundry
1043 bz[ 63] = bz[ 64] = 625.0; // first chip boundry
1044 bz[ 95] = bz[ 96] = 625.0; // first chip boundry
1045 bz[127] = bz[128] = 625.0; // first chip boundry
1046 bz[160] = 425.0;// Set so that there is no zero pixel size for fNz.
1047 segSPD->SetBinSize(bx,bz); // Based on AliITSgeomSPD for now.
1048 geom->ReSetShape(idet,segSPD);
1051 initSeg[idet] = kTRUE;
1052 AliITSsegmentationSDD *segSDD = new AliITSsegmentationSDD();
1053 segSDD->SetDetSize(shapePar[0]*kcm2micron, // X
1054 2.*shapePar[2]*kcm2micron, // Z
1055 2.*shapePar[1]*kcm2micron);// Y Microns
1056 segSDD->SetNPads(256,256);// Anodes, Samples
1057 geom->ReSetShape(idet,segSDD);
1060 initSeg[idet] = kTRUE;
1061 AliITSsegmentationSSD *segSSD = new AliITSsegmentationSSD();
1062 segSSD->SetDetSize(2.*shapePar[0]*kcm2micron, // X
1063 2.*shapePar[2]*kcm2micron, // Z
1064 2.*shapePar[1]*kcm2micron);// Y Microns.
1065 segSSD->SetPadSize(95.,0.); // strip x pitch in microns
1066 segSSD->SetNPads(768,2); // number of strips on each side, sides.
1067 segSSD->SetAngles(0.0075,0.0275); // strip angels rad P and N side.
1068 segSSD->SetAnglesLay5(0.0075,0.0275);//strip angels rad P and N
1069 segSSD->SetAnglesLay6(0.0275,0.0075);//strip angels rad P and N
1070 geom->ReSetShape(idet,segSSD);
1072 default:{// Others, Note no kSDDp or kSSDp in this geometry.
1073 geom->ReSetShape(idet,0);
1074 Info("InitSegmentationPPRasymmFMD",
1075 "default segmentation Dx=%f Dy=%f Dz=%f default=%d",
1076 shapePar[0],shapePar[1],shapePar[2],idet);
1081 //______________________________________________________________________
1082 Bool_t AliITSInitGeometry::GetTransformation(const TString &volumePath,
1084 // Returns the Transformation matrix between the volume specified
1085 // by the path volumePath and the Top or mater volume. The format
1086 // of the path volumePath is as follows (assuming ALIC is the Top volume)
1087 // "/ALIC_1/DDIP_1/S05I_2/S05H_1/S05G_3". Here ALIC is the top most
1088 // or master volume which has only 1 instance of. Of all of the daughter
1089 // volumes of ALICE, DDIP volume copy #1 is indicated. Similarly for
1090 // the daughter volume of DDIP is S05I copy #2 and so on.
1092 // TString& volumePath The volume path to the specific volume
1093 // for which you want the matrix. Volume name
1094 // hierarchy is separated by "/" while the
1095 // copy number is appended using a "_".
1097 // TGeoHMatrix &mat A matrix with its values set to those
1098 // appropriate to the Local to Master transformation
1100 // A logical value if kFALSE then an error occurred and no change to
1103 // We have to preserve the modeler state
1105 // Preserve the modeler state.
1106 gGeoManager->PushPath();
1107 if (!gGeoManager->cd(volumePath.Data())) {
1108 gGeoManager->PopPath();
1109 Error("GetTransformation","Error in cd-ing to ",volumePath.Data());
1111 } // end if !gGeoManager
1112 mat = *gGeoManager->GetCurrentMatrix();
1113 // Retstore the modeler state.
1114 gGeoManager->PopPath();
1117 //______________________________________________________________________
1118 Bool_t AliITSInitGeometry::GetShape(const TString &volumePath,
1119 TString &shapeType,TArrayD &par){
1120 // Returns the shape and its parameters for the volume specified
1123 // TString& volumeName The volume name
1125 // TString &shapeType Shape type
1126 // TArrayD &par A TArrayD of parameters with all of the
1127 // parameters of the specified shape.
1129 // A logical indicating whether there was an error in getting this
1132 gGeoManager->PushPath();
1133 if (!gGeoManager->cd(volumePath.Data())) {
1134 gGeoManager->PopPath();
1137 TGeoVolume * vol = gGeoManager->GetCurrentVolume();
1138 gGeoManager->PopPath();
1139 if (!vol) return kFALSE;
1140 TGeoShape *shape = vol->GetShape();
1141 TClass *classType = shape->IsA();
1142 if (classType==TGeoBBox::Class()) {
1146 TGeoBBox *box = (TGeoBBox*)shape;
1147 par.AddAt(box->GetDX(),0);
1148 par.AddAt(box->GetDY(),1);
1149 par.AddAt(box->GetDZ(),2);
1152 if (classType==TGeoTrd1::Class()) {
1156 TGeoTrd1 *trd1 = (TGeoTrd1*)shape;
1157 par.AddAt(trd1->GetDx1(),0);
1158 par.AddAt(trd1->GetDx2(),1);
1159 par.AddAt(trd1->GetDy(), 2);
1160 par.AddAt(trd1->GetDz(), 3);
1163 if (classType==TGeoTrd2::Class()) {
1167 TGeoTrd2 *trd2 = (TGeoTrd2*)shape;
1168 par.AddAt(trd2->GetDx1(),0);
1169 par.AddAt(trd2->GetDx2(),1);
1170 par.AddAt(trd2->GetDy1(),2);
1171 par.AddAt(trd2->GetDy2(),3);
1172 par.AddAt(trd2->GetDz(), 4);
1175 if (classType==TGeoTrap::Class()) {
1179 TGeoTrap *trap = (TGeoTrap*)shape;
1180 Double_t tth = TMath::Tan(trap->GetTheta()*TMath::DegToRad());
1181 par.AddAt(trap->GetDz(),0);
1182 par.AddAt(tth*TMath::Cos(trap->GetPhi()*TMath::DegToRad()),1);
1183 par.AddAt(tth*TMath::Sin(trap->GetPhi()*TMath::DegToRad()),2);
1184 par.AddAt(trap->GetH1(),3);
1185 par.AddAt(trap->GetBl1(),4);
1186 par.AddAt(trap->GetTl1(),5);
1187 par.AddAt(TMath::Tan(trap->GetAlpha1()*TMath::DegToRad()),6);
1188 par.AddAt(trap->GetH2(),7);
1189 par.AddAt(trap->GetBl2(),8);
1190 par.AddAt(trap->GetTl2(),9);
1191 par.AddAt(TMath::Tan(trap->GetAlpha2()*TMath::DegToRad()),10);
1194 if (classType==TGeoTube::Class()) {
1198 TGeoTube *tube = (TGeoTube*)shape;
1199 par.AddAt(tube->GetRmin(),0);
1200 par.AddAt(tube->GetRmax(),1);
1201 par.AddAt(tube->GetDz(),2);
1204 if (classType==TGeoTubeSeg::Class()) {
1208 TGeoTubeSeg *tubs = (TGeoTubeSeg*)shape;
1209 par.AddAt(tubs->GetRmin(),0);
1210 par.AddAt(tubs->GetRmax(),1);
1211 par.AddAt(tubs->GetDz(),2);
1212 par.AddAt(tubs->GetPhi1(),3);
1213 par.AddAt(tubs->GetPhi2(),4);
1216 if (classType==TGeoCone::Class()) {
1220 TGeoCone *cone = (TGeoCone*)shape;
1221 par.AddAt(cone->GetDz(),0);
1222 par.AddAt(cone->GetRmin1(),1);
1223 par.AddAt(cone->GetRmax1(),2);
1224 par.AddAt(cone->GetRmin2(),3);
1225 par.AddAt(cone->GetRmax2(),4);
1228 if (classType==TGeoConeSeg::Class()) {
1232 TGeoConeSeg *cons = (TGeoConeSeg*)shape;
1233 par.AddAt(cons->GetDz(),0);
1234 par.AddAt(cons->GetRmin1(),1);
1235 par.AddAt(cons->GetRmax1(),2);
1236 par.AddAt(cons->GetRmin2(),3);
1237 par.AddAt(cons->GetRmax2(),4);
1238 par.AddAt(cons->GetPhi1(),5);
1239 par.AddAt(cons->GetPhi2(),6);
1242 if (classType==TGeoSphere::Class()) {
1247 TGeoSphere *sphe = (TGeoSphere*)shape;
1248 par.AddAt(sphe->GetRmin(),0);
1249 par.AddAt(sphe->GetRmax(),1);
1250 par.AddAt(sphe->GetTheta1(),2);
1251 par.AddAt(sphe->GetTheta2(),3);
1252 par.AddAt(sphe->GetPhi1(),4);
1253 par.AddAt(sphe->GetPhi2(),5);
1256 if (classType==TGeoPara::Class()) {
1260 TGeoPara *para = (TGeoPara*)shape;
1261 par.AddAt(para->GetX(),0);
1262 par.AddAt(para->GetY(),1);
1263 par.AddAt(para->GetZ(),2);
1264 par.AddAt(para->GetTxy(),3);
1265 par.AddAt(para->GetTxz(),4);
1266 par.AddAt(para->GetTyz(),5);
1269 if (classType==TGeoPgon::Class()) {
1271 TGeoPgon *pgon = (TGeoPgon*)shape;
1272 Int_t nz = pgon->GetNz();
1273 const Double_t *rmin = pgon->GetRmin();
1274 const Double_t *rmax = pgon->GetRmax();
1275 const Double_t *z = pgon->GetZ();
1278 par.AddAt(pgon->GetPhi1(),0);
1279 par.AddAt(pgon->GetDphi(),1);
1280 par.AddAt(pgon->GetNedges(),2);
1281 par.AddAt(pgon->GetNz(),3);
1282 for (Int_t i=0; i<nz; i++) {
1283 par.AddAt(z[i], 4+3*i);
1284 par.AddAt(rmin[i], 4+3*i+1);
1285 par.AddAt(rmax[i], 4+3*i+2);
1289 if (classType==TGeoPcon::Class()) {
1291 TGeoPcon *pcon = (TGeoPcon*)shape;
1292 Int_t nz = pcon->GetNz();
1293 const Double_t *rmin = pcon->GetRmin();
1294 const Double_t *rmax = pcon->GetRmax();
1295 const Double_t *z = pcon->GetZ();
1298 par.AddAt(pcon->GetPhi1(),0);
1299 par.AddAt(pcon->GetDphi(),1);
1300 par.AddAt(pcon->GetNz(),2);
1301 for (Int_t i=0; i<nz; i++) {
1302 par.AddAt(z[i], 3+3*i);
1304 par.AddAt(rmin[i], 3+3*i+1);
1305 par.AddAt(rmax[i], 3+3*i+2);
1309 if (classType==TGeoEltu::Class()) {
1313 TGeoEltu *eltu = (TGeoEltu*)shape;
1314 par.AddAt(eltu->GetA(),0);
1315 par.AddAt(eltu->GetB(),1);
1316 par.AddAt(eltu->GetDz(),2);
1319 if (classType==TGeoHype::Class()) {
1323 TGeoHype *hype = (TGeoHype*)shape;
1324 par.AddAt(TMath::Sqrt(hype->RadiusHypeSq(0.,kTRUE)),0);
1325 par.AddAt(TMath::Sqrt(hype->RadiusHypeSq(0.,kFALSE)),1);
1326 par.AddAt(hype->GetDZ(),2);
1327 par.AddAt(hype->GetStIn(),3);
1328 par.AddAt(hype->GetStOut(),4);
1331 if (classType==TGeoGtra::Class()) {
1335 TGeoGtra *trap = (TGeoGtra*)shape;
1336 Double_t tth = TMath::Tan(trap->GetTheta()*TMath::DegToRad());
1337 par.AddAt(trap->GetDz(),0);
1338 par.AddAt(tth*TMath::Cos(trap->GetPhi()*TMath::DegToRad()),1);
1339 par.AddAt(tth*TMath::Sin(trap->GetPhi()*TMath::DegToRad()),2);
1340 par.AddAt(trap->GetH1(),3);
1341 par.AddAt(trap->GetBl1(),4);
1342 par.AddAt(trap->GetTl1(),5);
1343 par.AddAt(TMath::Tan(trap->GetAlpha1()*TMath::DegToRad()),6);
1344 par.AddAt(trap->GetH2(),7);
1345 par.AddAt(trap->GetBl2(),8);
1346 par.AddAt(trap->GetTl2(),9);
1347 par.AddAt(TMath::Tan(trap->GetAlpha2()*TMath::DegToRad()),10);
1348 par.AddAt(trap->GetTwistAngle(),11);
1351 if (classType==TGeoCtub::Class()) {
1355 TGeoCtub *ctub = (TGeoCtub*)shape;
1356 const Double_t *lx = ctub->GetNlow();
1357 const Double_t *tx = ctub->GetNhigh();
1358 par.AddAt(ctub->GetRmin(),0);
1359 par.AddAt(ctub->GetRmax(),1);
1360 par.AddAt(ctub->GetDz(),2);
1361 par.AddAt(ctub->GetPhi1(),3);
1362 par.AddAt(ctub->GetPhi2(),4);
1368 par.AddAt(tx[2],10);
1371 Error("GetShape","Getting shape parameters for shape %s not implemented",
1372 shape->ClassName());
1373 shapeType = "Unknown";
1376 //______________________________________________________________________
1377 void AliITSInitGeometry::DecodeDetector(
1378 Int_t &mod,Int_t layer,Int_t cpn0,Int_t cpn1,Int_t cpn2) const {
1379 // decode geometry into detector module number. There are two decoding
1380 // Scheams. Old which does not follow the ALICE coordinate system
1381 // requirements, and New which dose.
1383 // Int_t layer The ITS layer
1384 // Int_t cpn0 The lowest copy number
1385 // Int_t cpn1 The middle copy number
1386 // Int_t cpn2 the highest copy number
1388 // Int_t &mod The module number assoicated with this set
1393 // This is a FIXED switch yard function. I (Bjorn Nilsen) Don't
1394 // like them but I see not better way for the moment.
1395 switch (fMajorVersion){
1397 if(GetMinorVersion()==1)
1398 return DecodeDetectorvPPRasymmFMD(mod,layer,cpn0,cpn1,cpn2);
1399 else if(GetMinorVersion()==2)
1400 return DecodeDetectorvtest2(mod,layer,cpn0,cpn1,cpn2);
1401 Warning("DecodeDetector",
1402 "Geometry is kvtest minor version=%d is not defined",
1406 Error("DecodeDetector","Major version = kvDefault, not supported");
1409 return DecodeDetectorvSPD02(mod,layer,cpn0,cpn1,cpn2);
1412 return DecodeDetectorvSDD03(mod,layer,cpn0,cpn1,cpn2);
1415 return DecodeDetectorvSSD03(mod,layer,cpn0,cpn1,cpn2);
1418 return DecodeDetectorvITS04(mod,layer,cpn0,cpn1,cpn2);
1420 case kvPPRcourseasymm:{
1421 return DecodeDetectorvPPRcourseasymm(mod,layer,cpn0,cpn1,cpn2);
1423 case kvPPRasymmFMD:{
1424 return DecodeDetectorvPPRasymmFMD(mod,layer,cpn0,cpn1,cpn2);
1427 return DecodeDetectorv11(mod,layer,cpn0,cpn1,cpn2);
1430 return DecodeDetectorv11Hybrid(mod,layer,cpn0,cpn1,cpn2);
1433 Error("DecodeDetector","Major version = %d, not supported",
1434 (Int_t)fMajorVersion);
1440 //______________________________________________________________________
1441 void AliITSInitGeometry::RecodeDetector(Int_t mod,Int_t &cpn0,
1442 Int_t &cpn1,Int_t &cpn2){
1443 // decode geometry into detector module number. There are two decoding
1444 // Scheams. Old which does not follow the ALICE coordinate system
1445 // requirements, and New which dose.
1447 // Int_t mod The module number assoicated with this set
1450 // Int_t cpn0 The lowest copy number
1451 // Int_t cpn1 The middle copy number
1452 // Int_t cpn2 the highest copy number
1456 // This is a FIXED switch yard function. I (Bjorn Nilsen) Don't
1457 // like them but I see not better way for the moment.
1458 switch (fMajorVersion){
1460 if(GetMinorVersion()==1)
1461 return RecodeDetectorvPPRasymmFMD(mod,cpn0,cpn1,cpn2);
1462 else if(GetMinorVersion()==2)
1463 return RecodeDetectorvtest2(mod,cpn0,cpn1,cpn2);
1464 Warning("RecodeDetector",
1465 "Geometry is kvtest minor version=%d is not defined",
1470 Error("RecodeDetector","Major version = kvDefault, not supported");
1474 return RecodeDetectorvSPD02(mod,cpn0,cpn1,cpn2);
1477 return RecodeDetectorvSDD03(mod,cpn0,cpn1,cpn2);
1480 return RecodeDetectorvSSD03(mod,cpn0,cpn1,cpn2);
1483 return RecodeDetectorvITS04(mod,cpn0,cpn1,cpn2);
1485 case kvPPRcourseasymm:{
1486 return RecodeDetectorvPPRcourseasymm(mod,cpn0,cpn1,cpn2);
1488 case kvPPRasymmFMD:{
1489 return RecodeDetectorvPPRasymmFMD(mod,cpn0,cpn1,cpn2);
1492 return RecodeDetectorv11(mod,cpn0,cpn1,cpn2);
1495 return RecodeDetectorv11Hybrid(mod,cpn0,cpn1,cpn2);
1498 Error("RecodeDetector","Major version = %d, not supported",
1499 (Int_t)fMajorVersion);
1505 //______________________________________________________________________
1506 void AliITSInitGeometry::DecodeDetectorLayers(Int_t mod,Int_t &layer,
1507 Int_t &lad,Int_t &det){
1508 // decode geometry into detector module number. There are two decoding
1509 // Scheams. Old which does not follow the ALICE coordinate system
1510 // requirements, and New which dose. Note, this use of layer ladder
1511 // and detector numbers are strictly for internal use of this
1512 // specific code. They do not represent the "standard" layer ladder
1513 // or detector numbering except in a very old and obsoleate sence.
1515 // Int_t mod The module number assoicated with this set
1518 // Int_t lay The layer number
1519 // Int_t lad The ladder number
1520 // Int_t det the dettector number
1524 // This is a FIXED switch yard function. I (Bjorn Nilsen) Don't
1525 // like them but I see not better way for the moment.
1526 switch (fMajorVersion){
1528 if(GetMinorVersion()==1)
1529 return DecodeDetectorLayersvPPRasymmFMD(mod,layer,lad,det);
1530 else if(GetMinorVersion()==2)
1531 return DecodeDetectorLayersvtest2(mod,layer,lad,det);
1532 Warning("DecodeDetectorLayers",
1533 "Geometry is kvtest minor version=%d is not defined",
1538 Error("DecodeDetectorLayers",
1539 "Major version = kvDefault, not supported");
1543 return DecodeDetectorLayersvSPD02(mod,layer,lad,det);
1546 return DecodeDetectorLayersvSDD03(mod,layer,lad,det);
1549 return DecodeDetectorLayersvSSD03(mod,layer,lad,det);
1552 return DecodeDetectorLayersvITS04(mod,layer,lad,det);
1554 case kvPPRcourseasymm:{
1555 return DecodeDetectorLayersvPPRcourseasymm(mod,layer,lad,det);
1557 case kvPPRasymmFMD:{
1558 return DecodeDetectorLayersvPPRasymmFMD(mod,layer,lad,det);
1561 return DecodeDetectorLayersv11(mod,layer,lad,det);
1564 return DecodeDetectorLayersv11Hybrid(mod,layer,lad,det);
1567 Error("DecodeDetectorLayers","Major version = %d, not supported",
1568 (Int_t)fMajorVersion);
1574 //______________________________________________________________________
1575 void AliITSInitGeometry::DecodeDetectorvSPD02(
1576 Int_t &mod,Int_t ncpn,Int_t cpy0,Int_t cpy1,Int_t cpy2) const {
1577 // decode geometry into detector module number
1579 // Int_t ncpn The Number of copies of this volume
1580 // Int_t cpy0 The lowest copy number
1581 // Int_t cpy1 The middle copy number
1582 // Int_t cpy2 the highest copy number
1584 // Int_t &mod The module number assoicated with this set
1589 // detector = ladder = 1
1590 if(ncpn==4 && cpy1>2) mod = cpy1; // layer = 1,2
1591 else mod = cpy1-1; // layer = 4,5
1592 if(ncpn==1) mod = 2; // layer=3
1596 //______________________________________________________________________
1597 void AliITSInitGeometry::RecodeDetectorvSPD02(Int_t mod,Int_t &cpn0,
1598 Int_t &cpn1,Int_t &cpn2){
1599 // decode geometry into detector module number. There are two decoding
1600 // Scheams. Old which does not follow the ALICE coordinate system
1601 // requirements, and New which dose.
1603 // Int_t mod The module number assoicated with this set
1606 // Int_t cpn0 The lowest copy number
1607 // Int_t cpn1 The middle copy number
1608 // Int_t cpn2 the highest copy number
1626 //______________________________________________________________________
1627 void AliITSInitGeometry::DecodeDetectorLayersvSPD02(Int_t mod,Int_t &lay,
1628 Int_t &lad,Int_t &det){
1629 // decode geometry into detector module number. There are two decoding
1630 // Scheams. Old which does not follow the ALICE coordinate system
1631 // requirements, and New which dose. Note, this use of layer ladder
1632 // and detector numbers are strictly for internal use of this
1633 // specific code. They do not represent the "standard" layer ladder
1634 // or detector numbering except in a very old and obsoleate sence.
1636 // Int_t mod The module number assoicated with this set
1639 // Int_t lay The layer number
1640 // Int_t lad The ladder number
1641 // Int_t det the dettector number
1649 //______________________________________________________________________
1650 void AliITSInitGeometry::DecodeDetectorvSDD03(
1651 Int_t &mod,Int_t ncpys,Int_t cpy0,Int_t cpy1,Int_t cpy2) const {
1652 // decode geometry into detector module number. There are two decoding
1653 // Scheams. Old which does not follow the ALICE coordinate system
1654 // requirements, and New which dose.
1656 // Int_t ncpys The number of posible copies cpn1
1657 // Int_t cpy0 The lowest copy number
1658 // Int_t cpy1 The middle copy number
1659 // Int_t cpy2 the highest copy number
1661 // Int_t &mod The module number assoicated with this set
1666 if(ncpys==10){ // ITEL detectors
1667 if(cpy1>4) mod = cpy1+1;
1669 }else{ // IDET detectors
1670 if(cpy1==1) mod = 4;
1676 //______________________________________________________________________
1677 void AliITSInitGeometry::RecodeDetectorvSDD03(Int_t mod,Int_t &cpn0,
1678 Int_t &cpn1,Int_t &cpn2){
1679 // decode geometry into detector module number. There are two decoding
1680 // Scheams. Old which does not follow the ALICE coordinate system
1681 // requirements, and New which dose.
1683 // Int_t mod The module number assoicated with this set
1686 // Int_t cpn0 The lowest copy number
1687 // Int_t cpn1 The middle copy number
1688 // Int_t cpn2 the highest copy number
1694 if(mod<4) cpn1 = mod+1;
1695 else if(mod==4||mod==5) cpn1 = mod-3;
1699 //______________________________________________________________________
1700 void AliITSInitGeometry::DecodeDetectorLayersvSDD03(Int_t mod,Int_t &lay,
1701 Int_t &lad,Int_t &det){
1702 // decode geometry into detector module number. There are two decoding
1703 // Scheams. Old which does not follow the ALICE coordinate system
1704 // requirements, and New which dose. Note, this use of layer ladder
1705 // and detector numbers are strictly for internal use of this
1706 // specific code. They do not represent the "standard" layer ladder
1707 // or detector numbering except in a very old and obsoleate sence.
1709 // Int_t mod The module number assoicated with this set
1712 // Int_t lay The layer number
1713 // Int_t lad The ladder number
1714 // Int_t det the dettector number
1722 //______________________________________________________________________
1723 void AliITSInitGeometry::DecodeDetectorvSSD03(
1724 Int_t &mod,Int_t dtype,Int_t cpn0,Int_t cpn1,Int_t cpn2) const {
1725 // decode geometry into detector module number. There are two decoding
1726 // Scheams. Old which does not follow the ALICE coordinate system
1727 // requirements, and New which dose.
1729 // Int_t dtype The detector type 1=ITSA 2=IGAR 3=IFRA
1730 // Int_t cpn0 The lowest copy number
1731 // Int_t cpn1 The middle copy number
1732 // Int_t cpn2 the highest copy number
1734 // Int_t &mod The module number assoicated with this set
1739 if(dtype==2){mod=2; return;}
1740 if(dtype==3){mod=3; return;}
1742 if(cpn0==3) mod = 4;
1746 //______________________________________________________________________
1747 void AliITSInitGeometry::RecodeDetectorvSSD03(Int_t mod,Int_t &cpn0,
1748 Int_t &cpn1,Int_t &cpn2){
1749 // decode geometry into detector module number. There are two decoding
1750 // Scheams. Old which does not follow the ALICE coordinate system
1751 // requirements, and New which dose.
1753 // Int_t mod The module number assoicated with this set
1756 // Int_t cpn0 The lowest copy number
1757 // Int_t cpn1 The middle copy number
1758 // Int_t cpn2 the highest copy number
1764 if(mod<2) cpn0=mod+1;
1765 else if (mod==2||mod==3) cpn0=1;
1769 //______________________________________________________________________
1770 void AliITSInitGeometry::DecodeDetectorLayersvSSD03(Int_t mod,Int_t &lay,
1771 Int_t &lad,Int_t &det){
1772 // decode geometry into detector module number. There are two decoding
1773 // Scheams. Old which does not follow the ALICE coordinate system
1774 // requirements, and New which dose. Note, this use of layer ladder
1775 // and detector numbers are strictly for internal use of this
1776 // specific code. They do not represent the "standard" layer ladder
1777 // or detector numbering except in a very old and obsoleate sence.
1779 // Int_t mod The module number assoicated with this set
1782 // Int_t lay The layer number
1783 // Int_t lad The ladder number
1784 // Int_t det the dettector number
1792 //______________________________________________________________________
1793 void AliITSInitGeometry::DecodeDetectorvITS04(
1794 Int_t &mod,Int_t dtype,Int_t cpn0,Int_t cpn1,Int_t cpn2) const {
1795 // decode geometry into detector module number. There are two decoding
1796 // Scheams. Old which does not follow the ALICE coordinate system
1797 // requirements, and New which dose.
1799 // Int_t dtype The detector type 1=ITSA 2=IGAR 3=IFRA
1800 // Int_t cpn0 The lowest copy number
1801 // Int_t cpn1 The middle copy number
1802 // Int_t cpn2 the highest copy number
1804 // Int_t &mod The module number assoicated with this set
1813 //______________________________________________________________________
1814 void AliITSInitGeometry::RecodeDetectorvITS04(Int_t mod,Int_t &cpn0,
1815 Int_t &cpn1,Int_t &cpn2){
1816 // decode geometry into detector module number. There are two decoding
1817 // Scheams. Old which does not follow the ALICE coordinate system
1818 // requirements, and New which dose.
1820 // Int_t mod The module number assoicated with this set
1823 // Int_t cpn0 The lowest copy number
1824 // Int_t cpn1 The middle copy number
1825 // Int_t cpn2 the highest copy number
1831 case 0:case 1:case 2:case 3:{
1837 case 6:case 7:case 8:case 9:{
1846 //______________________________________________________________________
1847 void AliITSInitGeometry::DecodeDetectorLayersvITS04(Int_t mod,Int_t &lay,
1848 Int_t &lad,Int_t &det){
1849 // decode geometry into detector module number. There are two decoding
1850 // Scheams. Old which does not follow the ALICE coordinate system
1851 // requirements, and New which dose. Note, this use of layer ladder
1852 // and detector numbers are strictly for internal use of this
1853 // specific code. They do not represent the "standard" layer ladder
1854 // or detector numbering except in a very old and obsoleate sence.
1856 // Int_t mod The module number assoicated with this set
1859 // Int_t lay The layer number
1860 // Int_t lad The ladder number
1861 // Int_t det the dettector number
1867 case 0:case 1:case 2:case 3:{
1874 case 6:case 7:case 8:case 9:{
1885 //______________________________________________________________________
1886 void AliITSInitGeometry::DecodeDetectorvPPRasymmFMD(Int_t &mod,Int_t layer,Int_t cpn0,
1887 Int_t cpn1,Int_t cpn2) const {
1888 // decode geometry into detector module number. There are two decoding
1889 // Scheams. Old which does not follow the ALICE coordinate system
1890 // requirements, and New which dose.
1892 // Int_t layer The ITS layer
1893 // Int_t cpn0 The lowest copy number
1894 // Int_t cpn1 The middle copy number
1895 // Int_t cpn2 the highest copy number
1897 // Int_t &mod The module number assoicated with this set
1901 const Int_t kDetPerLadderSPD[2]={2,4};
1902 const Int_t kDetPerLadder[6]={4,4,6,8,22,25};
1903 const Int_t kLadPerLayer[6]={20,40,14,22,34,38};
1904 Int_t lay=-1,lad=-1,det=-1,i;
1906 if(fDecode){ // New decoding scheam
1911 if(cpn0==4&&cpn1==1) lad=1;
1912 else if(cpn0==4&&cpn1==2) lad=20;
1914 lad = 8-cpn1-kDetPerLadderSPD[layer-1]*(cpn0-1);
1916 lad = 28-cpn1-kDetPerLadderSPD[layer-1]*(cpn0-1);
1922 if(cpn0==4&&cpn1==1) lad=1;
1924 lad = 14-cpn1-kDetPerLadderSPD[layer-1]*(cpn0-1);
1926 lad = 54-cpn1-kDetPerLadderSPD[layer-1]*(cpn0-1);
1931 if(cpn0<5) lad = 5-cpn0;
1937 if(cpn0<7) lad = 7-cpn0;
1943 if(cpn0<10) lad = 10-cpn0;
1949 if(cpn0<9) lad = 9-cpn0;
1955 for(i=0;i<layer-1;i++) mod += kLadPerLayer[i]*kDetPerLadder[i];
1956 mod += kDetPerLadder[layer-1]*(lad-1)+det-1;// module start at zero.
1959 // Old decoding scheam
1963 lad = cpn1+kDetPerLadderSPD[layer-1]*(cpn0-1);
1980 for(i=0;i<layer-1;i++) mod += kLadPerLayer[i]*kDetPerLadder[i];
1981 mod += kDetPerLadder[layer-1]*(lad-1)+det-1;// module start at zero.
1984 //______________________________________________________________________
1985 void AliITSInitGeometry::RecodeDetectorvPPRasymmFMD(Int_t mod,Int_t &cpn0,
1986 Int_t &cpn1,Int_t &cpn2){
1987 // decode geometry into detector module number. There are two decoding
1988 // Scheams. Old which does not follow the ALICE coordinate system
1989 // requirements, and New which dose.
1991 // Int_t mod The module number assoicated with this set
1994 // Int_t cpn0 The lowest copy number
1995 // Int_t cpn1 The middle copy number
1996 // Int_t cpn2 the highest copy number
1999 const Int_t kITSgeoTreeCopys[6][3]= {{10, 2, 4},// lay=1
2005 const Int_t kDetPerLadderSPD[2]={2,4};
2006 // const Int_t kDetPerLadder[6]={4,4,6,8,22,25};
2007 // const Int_t kLadPerLayer[6]={20,40,14,22,34,38};
2010 cpn0 = cpn1 = cpn2 = 0;
2011 DecodeDetectorLayers(mod,lay,lad,det);
2012 if(fDecode){ // New decoding scheam
2015 cpn2 = 5-det; // Detector 1-4
2016 cpn1 = 1+(lad-1)%kDetPerLadderSPD[lay-1];
2017 cpn0 = 5-(lad+kDetPerLadderSPD[lay-1])/kDetPerLadderSPD[lay-1];
2018 if(mod>27) cpn0 = 15-(lad+kDetPerLadderSPD[lay-1])/
2019 kDetPerLadderSPD[lay-1];
2022 cpn2 = 5-det; // Detector 1-4
2023 cpn1 = 4-(lad+2)%kDetPerLadderSPD[lay-1];
2024 cpn0 = 1+(14-cpn1-lad)/kDetPerLadderSPD[lay-1];
2025 if(mod>131) cpn0 = 1+(54-lad-cpn1)/kDetPerLadderSPD[lay-1];
2029 if(lad<5) cpn0 = 5-lad;
2035 if(lad<7) cpn0 = 7-lad;
2041 if(lad<10) cpn0 = 10-lad;
2047 if(lad<9) cpn0 = 9-lad;
2052 Error("RecodeDetector","New: mod=%d lay=%d not 1-6.");
2056 if(cpn0<1||cpn1<1||cpn2<1||
2057 cpn0>kITSgeoTreeCopys[lay-1][0]||
2058 cpn1>kITSgeoTreeCopys[lay-1][1]||
2059 cpn2>kITSgeoTreeCopys[lay-1][2])
2060 Error("RecodeDetector",
2061 "cpn0=%d cpn1=%d cpn2=%d mod=%d lay=%d lad=%d det=%d",
2062 cpn0,cpn1,cpn2,mod,lay,lad,det);
2068 cpn2 = det; // Detector 1-4
2069 cpn0 = (lad+kDetPerLadderSPD[lay-1]-1)/kDetPerLadderSPD[lay-1];
2070 cpn1 = (lad+kDetPerLadderSPD[lay-1]-1)%kDetPerLadderSPD[lay-1] + 1;
2072 case 3: case 4: case 5 : case 6:{
2078 Error("RecodeDetector","Old: mod=%d lay=%d not 1-6.");
2082 if(cpn0<1||cpn1<1||cpn2<1||
2083 cpn0>kITSgeoTreeCopys[lay-1][0]||
2084 cpn1>kITSgeoTreeCopys[lay-1][1]||
2085 cpn2>kITSgeoTreeCopys[lay-1][2])
2086 Error("RecodeDetector",
2087 "cpn0=%d cpn1=%d cpn2=%d mod=%d lay=%d lad=%d det=%d",
2088 cpn0,cpn1,cpn2,mod,lay,lad,det);
2091 //______________________________________________________________________
2092 void AliITSInitGeometry::DecodeDetectorLayersvPPRasymmFMD(Int_t mod,Int_t &lay,
2093 Int_t &lad,Int_t &det){
2094 // decode geometry into detector module number. There are two decoding
2095 // Scheams. Old which does not follow the ALICE coordinate system
2096 // requirements, and New which dose. Note, this use of layer ladder
2097 // and detector numbers are strictly for internal use of this
2098 // specific code. They do not represent the "standard" layer ladder
2099 // or detector numbering except in a very old and obsoleate sence.
2101 // Int_t mod The module number assoicated with this set
2104 // Int_t lay The layer number
2105 // Int_t lad The ladder number
2106 // Int_t det the dettector number
2109 // const Int_t kDetPerLadderSPD[2]={2,4};
2110 const Int_t kDetPerLadder[6]={4,4,6,8,22,25};
2111 const Int_t kLadPerLayer[6]={20,40,14,22,34,38};
2119 mod2 += kLadPerLayer[lay]*kDetPerLadder[lay];
2121 }while(mod2<=mod); // end while
2122 if(lay>6||lay<1) Error("DecodeDetectorLayers","0<lay=%d>6",lay);
2123 mod2 -= kLadPerLayer[lay-1]*kDetPerLadder[lay-1];
2126 mod2 += kDetPerLadder[lay-1];
2127 }while(mod2<=mod); // end while
2128 if(lad>kLadPerLayer[lay-1]||lad<1) Error("DecodeDetectorLayers",
2129 "lad=%d>kLadPerLayer[lay-1=%d]=%d mod=%d mod2=%d",lad,lay-1,
2130 kLadPerLayer[lay-1],mod,mod2);
2131 mod2 -= kDetPerLadder[lay-1];
2133 if(det>kDetPerLadder[lay-1]||det<1) Error("DecodeDetectorLayers",
2134 "det=%d>detPerLayer[lay-1=%d]=%d mod=%d mod2=%d lad=%d",det,
2135 lay-1,kDetPerLadder[lay-1],mod,mod2,lad);
2138 //______________________________________________________________________
2139 void AliITSInitGeometry::DecodeDetectorv11Hybrid(Int_t &mod,Int_t layer,Int_t cpn0,
2140 Int_t cpn1,Int_t cpn2) const {
2141 // decode geometry into detector module number
2143 // Int_t layer The ITS layer
2144 // Int_t cpn0 The lowest copy number
2145 // Int_t cpn1 The middle copy number
2146 // Int_t cpn2 the highest copy number
2148 // Int_t &mod The module number assoicated with this set
2152 const Int_t kDetPerLadderSPD[2]={2,4};
2153 const Int_t kDetPerLadder[6]={4,4,6,8,22,25};
2154 const Int_t kLadPerLayer[6]={20,40,14,22,34,38};
2155 Int_t lad=-1,det=-1,i;
2159 lad = cpn1+kDetPerLadderSPD[layer-1]*(cpn0-1);
2163 if (SDDIsTGeoNative()) {
2172 if (SSDIsTGeoNative()) {
2184 for(i=0;i<layer-1;i++) mod += kLadPerLayer[i]*kDetPerLadder[i];
2185 mod += kDetPerLadder[layer-1]*(lad-1)+det-1;// module start at zero.
2188 //______________________________________________________________________
2189 void AliITSInitGeometry::RecodeDetectorv11Hybrid(Int_t mod,Int_t &cpn0,
2190 Int_t &cpn1,Int_t &cpn2) {
2191 // decode geometry into detector module number. There are two decoding
2192 // Scheams. Old which does not follow the ALICE coordinate system
2193 // requirements, and New which dose.
2195 // Int_t mod The module number assoicated with this set
2198 // Int_t cpn0 The lowest copy number
2199 // Int_t cpn1 The middle copy number
2200 // Int_t cpn2 the highest copy number
2203 const Int_t kITSgeoTreeCopys[6][3]= {{10, 2, 4},// lay=1
2209 const Int_t kDetPerLadderSPD[2]={2,4};
2210 // const Int_t kDetPerLadder[6]={4,4,6,8,22,25};
2211 // const Int_t kLadPerLayer[6]={20,40,14,22,34,38};
2214 cpn0 = cpn1 = cpn2 = 0;
2215 DecodeDetectorLayers(mod,lay,lad,det);
2219 cpn2 = det; // Detector 1-4
2220 cpn0 = (lad+kDetPerLadderSPD[lay-1]-1)/kDetPerLadderSPD[lay-1];
2221 cpn1 = (lad+kDetPerLadderSPD[lay-1]-1)%kDetPerLadderSPD[lay-1] + 1;
2223 case 3: case 4: case 5 : case 6:{
2229 Error("RecodeDetector","Old: mod=%d lay=%d not 1-6.");
2233 if(cpn0<1||cpn1<1||cpn2<1||
2234 cpn0>kITSgeoTreeCopys[lay-1][0]||
2235 cpn1>kITSgeoTreeCopys[lay-1][1]||
2236 cpn2>kITSgeoTreeCopys[lay-1][2])
2237 Error("RecodeDetector",
2238 "cpn0=%d cpn1=%d cpn2=%d mod=%d lay=%d lad=%d det=%d",
2239 cpn0,cpn1,cpn2,mod,lay,lad,det);
2242 //______________________________________________________________________
2243 void AliITSInitGeometry::DecodeDetectorLayersv11Hybrid(Int_t mod,Int_t &lay,
2244 Int_t &lad,Int_t &det) {
2246 // decode geometry into detector module number for v11Hybrid
2248 // Int_t mod The module number assoicated with this set
2251 // Int_t lay The layer number
2252 // Int_t lad The ladder number
2253 // Int_t det the dettector number
2257 const Int_t kDetPerLadder[6]={4,4,6,8,22,25};
2258 const Int_t kLadPerLayer[6]={20,40,14,22,34,38};
2266 mod2 += kLadPerLayer[lay]*kDetPerLadder[lay];
2268 } while(mod2<=mod); // end while
2269 if(lay>6||lay<1) Error("DecodeDetectorLayers","0<lay=%d>6",lay);
2270 mod2 -= kLadPerLayer[lay-1]*kDetPerLadder[lay-1];
2273 mod2 += kDetPerLadder[lay-1];
2274 } while(mod2<=mod); // end while
2275 if(lad>kLadPerLayer[lay-1]||lad<1) Error("DecodeDetectorLayers",
2276 "lad=%d>kLadPerLayer[lay-1=%d]=%d mod=%d mod2=%d",lad,lay-1,
2277 kLadPerLayer[lay-1],mod,mod2);
2278 mod2 -= kDetPerLadder[lay-1];
2280 if(det>kDetPerLadder[lay-1]||det<1) Error("DecodeDetectorLayers",
2281 "det=%d>detPerLayer[lay-1=%d]=%d mod=%d mod2=%d lad=%d",det,
2282 lay-1,kDetPerLadder[lay-1],mod,mod2,lad);
2286 //______________________________________________________________________
2287 Bool_t AliITSInitGeometry::WriteVersionString(Char_t *str,Int_t length,
2288 AliITSVersion_t maj,Int_t min,
2289 const Char_t *cvsDate,const Char_t *cvsRevision)const{
2290 // fills the string str with the major and minor version number
2292 // Char_t *str The character string to hold the major
2293 // and minor version numbers in
2294 // Int_t length The maximum number of characters which
2295 // can be accomidated by this string.
2296 // str[length-1] must exist and will be set to zero
2297 // AliITSVersion_t maj The major number
2298 // Int_t min The minor number
2299 // Char_t *cvsDate The date string from cvs
2300 // Char_t *cvsRevision The Revision string from cvs
2302 // Char_t *str The character string holding the major and minor
2303 // version numbers. str[length-1] must exist
2304 // and will be set to zero
2306 // kTRUE if no errors
2307 Int_t i,n,cvsDateLength,cvsRevisionLength;
2309 cvsDateLength = (Int_t)strlen(cvsDate);
2310 cvsRevisionLength = (Int_t)strlen(cvsRevision);
2312 n = 50+(Int_t)(TMath::Log10(TMath::Abs((Double_t)i)))+1+
2313 (Int_t)(TMath::Log10(TMath::Abs((Double_t)min)))+1
2314 +cvsDateLength-6+cvsRevisionLength-10;
2315 if(GetDebug()>1) printf("AliITSInitGeometry::WriteVersionString:"
2316 "length=%d major=%d minor=%d cvsDate=%s[%d] "
2317 "cvsRevision=%s[%d] n=%d\n",length,i,min,cvsDate,
2318 cvsDateLength,cvsRevision,cvsRevisionLength,n);
2321 if(length<n){// not enough space to write in output string.
2322 Warning("WriteVersionString","Output string not long enough "
2323 "lenght=%d must be at least %d long\n",length,n);
2325 } // end if length<n
2326 char *cvsrevision = new char[cvsRevisionLength-10];
2327 char *cvsdate = new char[cvsDateLength-6];
2328 for(i=0;i<cvsRevisionLength-10;i++)
2329 if(10+i<cvsRevisionLength-1)
2330 cvsrevision[i] = cvsRevision[10+i]; else cvsrevision[i] = 0;
2331 for(i=0;i<cvsDateLength-6;i++) if(6+i<cvsDateLength-1)
2332 cvsdate[i] = cvsDate[6+i]; else cvsdate[i] = 0;
2333 for(i=0;i<length;i++) str[i] = 0; // zero it out for now.
2335 sprintf(str,"Major Version= %d Minor Version= %d Revision: %s Date: %s",
2336 i,min,cvsrevision,cvsdate);
2337 if(GetDebug()>1)printf("AliITSInitGeometry::WriteVersionString: "
2338 "n=%d str=%s revision[%zu] date[%zu]\n",
2339 n,str,strlen(cvsrevision),strlen(cvsdate));
2340 delete[] cvsrevision;
2344 //______________________________________________________________________
2345 Bool_t AliITSInitGeometry::ReadVersionString(const Char_t *str,Int_t length,
2346 AliITSVersion_t &maj,Int_t &min,
2348 // fills the string str with the major and minor version number
2350 // Char_t *str The character string to holding the major and minor
2351 // version numbers in
2352 // Int_t length The maximum number of characters which can be
2353 // accomidated by this string. str[length-1] must exist
2355 // Char_t *str The character string holding the major and minor
2356 // version numbers unchanged. str[length-1] must exist.
2357 // AliITSVersion_t maj The major number
2358 // Int_t min The minor number
2359 // TDatime dt The date and time of the cvs commit
2361 // kTRUE if no errors
2363 Char_t cvsRevision[10],cvsDate[11],cvsTime[9];
2364 Int_t i,m,n=strlen(str),year,month,day,hours,minuits,seconds;
2366 if(GetDebug()>1)printf("AliITSInitGeometry::ReadVersionString:"
2367 "str=%s length=%d\n",
2369 if(n<35) return kFALSE; // not enough space for numbers
2370 m = sscanf(str,"Major Version= %d Minor Version= %d Revision: %s "
2371 "Date: %s %s",&i,&min,cvsRevision,cvsDate,cvsTime);
2374 m = sscanf(cvsDate,"%d/%d/%d",&year,&month,&day);
2377 m = sscanf(cvsTime,"%d:%d:%d",&hours,&minuits,&seconds);
2380 dt.Set(year,month,day,hours,minuits,seconds);
2381 if(GetDebug()>1)printf("AliITSInitGeometry::ReadVersionString: i=%d min=%d "
2382 "cvsRevision=%s cvsDate=%s cvsTime=%s m=%d\n",
2383 i,min,cvsRevision,cvsDate,cvsTime,m);
2384 if(GetDebug()>1)printf("AliITSInitGeometry::ReadVersionString: year=%d"
2385 " month=%d day=%d hours=%d minuits=%d seconds=%d\n",
2386 year,month,day,hours,minuits,seconds);
2400 case kvPPRasymmFMD:{
2401 maj = kvPPRasymmFMD;