+}
+//______________________________________________________________________
+void AliITSv5::InitAliITSgeom(){
+// Based on the geometry tree defined in Geant 3.21, this
+// routine initilizes the Class AliITSgeom from the Geant 3.21 ITS geometry
+// sturture.
+// if(gMC->IsA()!=TGeant3::Class()) {
+ if(strcmp(gMC->GetName(),"TGeant3")) {
+ Error("InitAliITSgeom",
+ "Wrong Monte Carlo. InitAliITSgeom uses TGeant3 calls");
+ return;
+ } // end if
+ cout << "Reading Geometry transformation directly from Geant 3." << endl;
+ const Int_t nlayers = 6;
+ const Int_t ndeep = 7;
+ Int_t itsGeomTreeNames[nlayers][ndeep],lnam[20],lnum[20];
+ Int_t nlad[nlayers],ndet[nlayers];
+ Double_t t[3],r[10];
+ Float_t par[20],att[20];
+ Int_t npar,natt,idshape,imat,imed;
+ AliITSGeant3Geometry *ig = new AliITSGeant3Geometry();
+ Int_t mod,lay,lad,det,i,j,k;
+ char *names[nlayers][ndeep] = {
+ {"ALIC","ITSV","ITSD","IT12","I132","I186","ITS1"}, // lay=1
+ {"ALIC","ITSV","ITSD","IT12","I132","I131","ITS2"}, // lay=2
+ {"ALIC","ITSV","ITSD","IT34","I004","I302","ITS3"}, // lay=3
+ {"ALIC","ITSV","ITSD","IT34","I005","I402","ITS4"}, // lay=4
+ {"ALIC","ITSV","ITSD","IT56","I565","I562","ITS5"}, // lay=5
+ {"ALIC","ITSV","ITSD","IT56","I569","I566","ITS6"}};// lay=6
+ Int_t itsGeomTreeCopys[nlayers][ndeep] = {{1,1,1,1,10, 2,4}, // lay=1
+ {1,1,1,1,10, 4,4}, // lay=2
+ {1,1,1,1,14, 6,1}, // lay=3
+ {1,1,1,1,22, 8,1}, // lay=4
+ {1,1,1,1,34,23,1}, // lay=5
+ {1,1,1,1,38,26,1}};// lay=6
+
+ // Sorry, but this is not very pritty code. It should be replaced
+ // at some point with a version that can search through the geometry
+ // tree its self.
+ for(i=0;i<20;i++) lnam[i] = lnum[i] = 0;
+ for(i=0;i<nlayers;i++)for(j=0;j<ndeep;j++)
+ itsGeomTreeNames[i][j] = ig->StringToInt(names[i][j]);
+ mod = 0;
+ for(i=0;i<nlayers;i++){
+ k = 1;
+ for(j=0;j<ndeep;j++) if(itsGeomTreeCopys[i][j]!=0)
+ k *= TMath::Abs(itsGeomTreeCopys[i][j]);
+ mod += k;
+ } // end for i
+
+ if(fITSgeom!=0) delete fITSgeom;
+ nlad[0]=20;nlad[1]=40;nlad[2]=14;nlad[3]=22;nlad[4]=34;nlad[5]=38;
+ ndet[0]=4;ndet[1]=4;ndet[2]=6;ndet[3]=8;ndet[4]=22;ndet[5]=25;
+ fITSgeom = new AliITSgeom(0,6,nlad,ndet,mod);
+ mod = -1;
+ for(lay=1;lay<=nlayers;lay++){
+ for(j=0;j<ndeep;j++) lnam[j] = itsGeomTreeNames[lay-1][j];
+ for(j=0;j<ndeep;j++) lnum[j] = itsGeomTreeCopys[lay-1][j];
+ switch (lay){
+ case 1: case 2: // layers 1 and 2 are a bit special
+ lad = 0;
+ for(j=1;j<=itsGeomTreeCopys[lay-1][4];j++){
+ lnum[4] = j;
+ for(k=1;k<=itsGeomTreeCopys[lay-1][5];k++){
+ lad++;
+ lnum[5] = k;
+ for(det=1;det<=itsGeomTreeCopys[lay-1][6];det++){
+ lnum[6] = det;
+ mod++;
+ ig->GetGeometry(ndeep,lnam,lnum,t,r,idshape,npar,natt,
+ par,att,imat,imed);
+ fITSgeom->CreatMatrix(mod,lay,lad,det,kSPD,t,r);
+ if(!(fITSgeom->IsShapeDefined((Int_t)kSPD)))
+ fITSgeom->ReSetShape(kSPD,
+ new AliITSgeomSPD300());
+ } // end for det
+ } // end for k
+ } // end for j
+ break;
+ case 3: case 4: case 5: case 6: // layers 3-6
+ lnum[6] = 1;
+ for(lad=1;lad<=itsGeomTreeCopys[lay-1][4];lad++){
+ lnum[4] = lad;
+ for(det=1;det<=itsGeomTreeCopys[lay-1][5];det++){
+ lnum[5] = det;
+ mod++;
+ ig->GetGeometry(7,lnam,lnum,t,r,idshape,npar,natt,
+ par,att,imat,imed);
+ switch (lay){
+ case 3: case 4:
+ fITSgeom->CreatMatrix(mod,lay,lad,det,kSDD,t,r);
+ if(!(fITSgeom->IsShapeDefined(kSDD)))
+ fITSgeom->ReSetShape(kSDD,new AliITSgeomSDD300());
+ break;
+ case 5: case 6:
+ fITSgeom->CreatMatrix(mod,lay,lad,det,kSSD,t,r);
+ if(!(fITSgeom->IsShapeDefined(kSSD)))
+ fITSgeom->ReSetShape(kSSD,new AliITSgeomSSD175());
+ break;
+ } // end switch
+ } // end for det
+ } // end for lad
+ break;
+ } // end switch
+ } // end for lay
+ return;