ITS/AliITSv5.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 /*
  17 $Log$
  18 Revision 1.36  2002/10/14 14:57:08  hristov
  19 Merging the VirtualMC branch to the main development branch (HEAD)
  20
  21 Revision 1.35.10.2  2002/08/30 15:45:54  alibrary
  22 Adding geant4vmc support
  23
  24 Revision 1.35.10.1  2002/06/10 17:51:15  hristov
  25 Merged with v3-08-02
  26
  27 Revision 1.35  2001/05/30 16:15:47  fca
  28 Correct comparison wiht AliGeant3::Class() introduced. Thanks to I.Hrivnacova
  29
  30 Revision 1.34  2001/05/30 15:55:35  hristov
  31 Strings compared instead of pointers
  32
  33 Revision 1.33  2001/05/30 14:04:31  hristov
  34 Dynamic cast replaced (F.Carminati)
  35
  36 Revision 1.32  2001/03/23 00:12:23  nilsen
  37 Set Reading of AliITSgeom data from Geant3 common blocks as the default and
  38 not a .det file. Removed redundent calls to BuildGeometry.
  39
  40 Revision 1.31  2001/02/13 16:53:35  nilsen
  41 Fixed a but when trying to use GEANT4. Needed to replace
  42 if(!((TGeant3*)gMC)) with if(!(dynamic_casst<TGeant3*>(gMC)))
  43 because just casting gMC to be TGeant3* even when it realy is a TGeant3 pointer
  44 did not result in a zero value. For AliITSv5asymm and AliITSv5symm, needed
  45 to fix a bug in the initilizers and a bug in BuildGeometry. This is now done
  46 in the same way as in AliITSv5.cxx.
  47
  48 Revision 1.30  2001/02/09 20:06:26  nilsen
  49 Fixed bug in distructor. Can't distroy fixxed length arrays. Thanks Peter.
  50
  51 Revision 1.29  2001/02/09 00:05:31  nilsen
  52 Added fMajor/MinorVersion variables and made other changes to better make
  53 use of the new code changes in AliITSgeom related classes.
  54
  55 Revision 1.28  2001/02/02 23:57:28  nilsen
  56 Added include file that are no londer included in AliITSgeom.h
  57
  58 Revision 1.27  2001/01/30 09:23:13  hristov
  59 Streamers removed (R.Brun)
  60
  61 Revision 1.26  2000/11/30 11:13:11  barbera
  62  Added changes suggested by Federico Carminati on nov, 30, 2000
  63
  64 Revision 1.25  2000/10/05 20:50:00  nilsen
  65 Now using root generated streamers.
  66
  67 Revision 1.14.4.12  2000/10/02 16:04:03  barbera
  68 Forward declarations added
  69
  70 Revision 1.22  2000/07/10 16:07:19  fca
  71 Release version of ITS code
  72
  73 Revision 1.14.4.4  2000/05/19 10:10:21  nilsen
  74 fix for bug with HP and Sun unix + fix for event display in ITS-working branch
  75
  76 Revision 1.14.4.3  2000/03/04 23:46:38  nilsen
  77 Fixed up the comments/documentation.
  78
  79 Revision 1.14.4.2  2000/03/02 21:53:02  nilsen
  80 To make it compatable with the changes in AliRun/AliModule.
  81
  82 Revision 1.14.4.1  2000/01/12 19:03:33  nilsen
  83 This is the version of the files after the merging done in December 1999.
  84 See the ReadMe110100.txt file for details
  85
  86 Revision 1.14  1999/10/22 08:16:49  fca
  87 Correct destructors, thanks to I.Hrivnacova
  88
  89 Revision 1.13  1999/10/06 10:15:19  fca
  90 Correct bug in allocation of layer name and add destructor
  91
  92 Revision 1.12  1999/10/05 08:05:09  fca
  93 Minor corrections for uninitialised variables.
  94
  95 Revision 1.11  1999/09/29 09:24:20  fca
  96 Introduction of the Copyright and cvs Log
  97
  98 */
  99
 100 ///////////////////////////////////////////////////////////////////////////////
 101 //
 102 //  Inner Traking System version 5
 103 //  This class contains the base procedures for the Inner Tracking System
 104 //
 105 // Authors: R. Barbera, B. S. Nilsen.
 106 // version 5.
 107 // Created September 17 1999.
 108 //
 109 ///////////////////////////////////////////////////////////////////////////////
 110 #include <Riostream.h>
 111 #include <stdio.h>
 112 #include <stdlib.h>
 113 #include <TMath.h>
 114 #include <TGeometry.h>
 115 #include <TNode.h>
 116 #include <TTUBE.h>
 117 #include <TFile.h>    // only required for Tracking function?
 118 #include <TCanvas.h>
 119 #include <TObjArray.h>
 120 #include <TLorentzVector.h>
 121 #include <TObjString.h>
 122 #include <TClonesArray.h>
 123 #include <TBRIK.h>
 124 #include <TSystem.h>
 125
 126 #include "AliMC.h"
 127 #include "AliRun.h"
 128 #include "AliITShit.h"
 129 #include "AliITSGeant3Geometry.h"
 130 #include "AliITS.h"
 131 #include "AliITSv5.h"
 132 #include "AliITSgeom.h"
 133 #include "AliITSgeomSPD.h"
 134 #include "AliITSgeomSDD.h"
 135 #include "AliITSgeomSSD.h"
 136
 137 ClassImp(AliITSv5)
 138
 139 //_____________________________________________________________________________
 140 AliITSv5::AliITSv5() {
 141 ////////////////////////////////////////////////////////////////////////
 142 //    Standard default constructor for the ITS version 5.
 143 ////////////////////////////////////////////////////////////////////////
 144     Int_t i;
 145
 146     fIdN    = 0;
 147     fIdName = 0;
 148     fIdSens = 0;
 149     fEuclidOut    = kFALSE; // Don't write Euclide file
 150     fGeomDetOut   = kFALSE; // Don't write .det file
 151     fGeomDetIn    = kFALSE; // Don't Read .det file
 152     fGeomOldDetIn = kFALSE; // Don't Read old formatted .det file
 153     fMajorVersion = IsVersion();
 154     fMinorVersion = 1;
 155     for(i=0;i<60;i++) fRead[i] = '\0';
 156     for(i=0;i<60;i++) fWrite[i] = '\0';
 157     for(i=0;i<60;i++) fEuclidGeomDet[i] = '\0';
 158 }
 159 //_____________________________________________________________________________
 160 AliITSv5::AliITSv5(const char *name, const char *title) : AliITS(name, title){
 161 ////////////////////////////////////////////////////////////////////////
 162 //    Standard constructor for the ITS version 5.
 163 ////////////////////////////////////////////////////////////////////////
 164     Int_t i;
 165
 166     fIdN    = 6;
 167     fIdName    = new TString[fIdN];
 168     fIdName[0] = "ITS1";
 169     fIdName[1] = "ITS2";
 170     fIdName[2] = "ITS3";
 171     fIdName[3] = "ITS4";
 172     fIdName[4] = "ITS5";
 173     fIdName[5] = "ITS6";
 174     fIdSens    = new Int_t[fIdN];
 175     for (i=0;i<fIdN;i++) fIdSens[i] = 0;
 176     fEuclidOut    = kFALSE; // Don't write Euclide file
 177     fGeomDetOut   = kFALSE; // Don't write .det file
 178     fGeomDetIn    = kFALSE; // Don't Read .det file
 179     fGeomOldDetIn = kFALSE; // Don't Read old formatted .det file
 180     fMajorVersion = IsVersion();
 181     fMinorVersion = 1;
 182     for(i=0;i<60;i++) fRead[i] = '\0';
 183     for(i=0;i<60;i++) fWrite[i] = '\0';
 184
 185     fEuclidMaterial = "$ALICE_ROOT/Euclid/ITSgeometry_5.tme";
 186     fEuclidGeometry = "$ALICE_ROOT/Euclid/ITSgeometry_5.euc";
 187     strncpy(fEuclidGeomDet,"$ALICE_ROOT/ITS/ITSgeometry_v5.det",60);
 188     strncpy(fRead,fEuclidGeomDet,60);
 189     strncpy(fWrite,fEuclidGeomDet,60);
 190 }
 191 //____________________________________________________________________________
 192 AliITSv5::AliITSv5(const AliITSv5 &source){
 193 ////////////////////////////////////////////////////////////////////////
 194 //     Copy Constructor for ITS version 5.
 195 ////////////////////////////////////////////////////////////////////////
 196     if(&source == this) return;
 197     Warning("Copy Constructor","Not allowed to copy AliITSv5");
 198     return;
 199 }
 200 //_____________________________________________________________________________
 201 AliITSv5& AliITSv5::operator=(const AliITSv5 &source){
 202 ////////////////////////////////////////////////////////////////////////
 203 //    Assignment operator for the ITS version 5.
 204 ////////////////////////////////////////////////////////////////////////
 205
 206     if(&source == this) return *this;
 207     Warning("= operator","Not allowed to copy AliITSv5");
 208     return *this;
 209
 210 }
 211 //_____________________________________________________________________________
 212 AliITSv5::~AliITSv5() {
 213 ////////////////////////////////////////////////////////////////////////
 214 //    Standard destructor for the ITS version 5.
 215 ////////////////////////////////////////////////////////////////////////
 216 }
 217 //______________________________________________________________________
 218 void AliITSv5::BuildGeometry(){
 219 ////////////////////////////////////////////////////////////////////////
 220 //    Geometry builder for the ITS version 5.
 221 ////////////////////////////////////////////////////////////////////////
 222   //
 223   // Build ITS TNODE geometry for event display using detailed geometry.
 224   // This function builds a simple ITS geometry used by the ROOT macro
 225   // ITSdisplay.C.
 226
 227   TNode *top;
 228   TNode *nd;
 229   //const int kColorITSSPD=kRed;
 230   //const int kColorITSSDD=kGreen;
 231   const int kColorITSSSD=kBlue;
 232   //
 233   AliITSgeom  *gm = this->GetITSgeom();
 234   if(gm==0) return;
 235   top=gAlice->GetGeometry()->GetNode("alice");
 236
 237   Int_t       lay,lad,det,i;
 238   Text_t      name[10];
 239   Float_t     xg[3];
 240   Float_t     rt[9];
 241   Double_t    rtd[9];
 242   TBRIK       *box;
 243   TRotMatrix  *rm;
 244   //TCanvas     *c1 = new TCanvas("c1","ITS");
 245
 246   for(lay=1;lay<=2;lay++)
 247    for(lad=1;lad<=gm->GetNladders(lay);lad++)
 248     for(det=1;det<=gm->GetNdetectors(lay);det++){
 249           try {
 250               box  = new TBRIK ("ActiveSPD","Active volume of SPD","SPD SI DET",
 251                                     0.64,0.0075,4.19);
 252           } catch (...) {
 253               cout << "EXCEPTION in box = new TBRIK" << endl;
 254               return;
 255           }
 256           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
 257           gm->GetRotMatrix(lay,lad,det,rt);
 258           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
 259           for(i=0;i<9;i++) rtd[i] = rt[i];
 260           try {
 261                 rm  = new TRotMatrix(name,name,rtd);
 262           } catch (...) {
 263                 cout << "EXCEPTION in   new TRotMatrix" << endl;
 264                 return;
 265           }
 266          top->cd();
 267           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
 268          try {
 269               nd  = new TNode("SPD"," ",box,xg[0],xg[1],xg[2],rm);
 270          } catch (...) {
 271               cout << "EXCEPTION in new TNode" << endl;
 272               return;
 273          }
 274          nd->SetLineColor(kColorITSSSD);
 275          fNodes->Add(nd);
 276     }
 277
 278   for(lay=3;lay<=3;lay++)
 279    for(lad=1;lad<=gm->GetNladders(lay);lad++)
 280     for(det=1;det<=gm->GetNdetectors(lay);det++){
 281           try {
 282               box  = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET",
 283                                     3.5,0.014,3.763);
 284           } catch (...) {
 285               cout << "EXCEPTION in box = new TBRIK" << endl;
 286               return;
 287           }
 288           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
 289           gm->GetRotMatrix(lay,lad,det,rt);
 290           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
 291           for(i=0;i<9;i++) rtd[i] = rt[i];
 292           try {
 293                 rm  = new TRotMatrix(name,name,rtd);
 294           } catch (...) {
 295                 cout << "EXCEPTION in   new TRotMatrix" << endl;
 296                 return;
 297           }
 298          top->cd();
 299           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
 300          try {
 301               nd  = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm);
 302          } catch (...) {
 303               cout << "EXCEPTION in new TNode" << endl;
 304               return;
 305          }
 306          nd->SetLineColor(kColorITSSSD);
 307          fNodes->Add(nd);
 308     }
 309
 310   for(lay=4;lay<=4;lay++)
 311    for(lad=1;lad<=gm->GetNladders(lay);lad++)
 312     for(det=1;det<=gm->GetNdetectors(lay);det++){
 313           try {
 314               box  = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET",
 315                                     3.5,0.014,3.763);
 316           } catch (...) {
 317               cout << "EXCEPTION in box = new TBRIK" << endl;
 318               return;
 319           }
 320           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
 321           gm->GetRotMatrix(lay,lad,det,rt);
 322           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
 323           for(i=0;i<9;i++) rtd[i] = rt[i];
 324           try {
 325                 rm  = new TRotMatrix(name,name,rtd);
 326           } catch (...) {
 327                 cout << "EXCEPTION in   new TRotMatrix" << endl;
 328                 return;
 329           }
 330          top->cd();
 331           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
 332          try {
 333               nd  = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm);
 334          } catch (...) {
 335               cout << "EXCEPTION in new TNode" << endl;
 336               return;
 337          }
 338          nd->SetLineColor(kColorITSSSD);
 339          fNodes->Add(nd);
 340     }
 341  for(lay=5;lay<=5;lay++)
 342    for(lad=1;lad<=gm->GetNladders(lay);lad++)
 343     for(det=1;det<=gm->GetNdetectors(lay);det++){
 344           try {
 345               box  = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET",
 346                                     3.65,0.015,2.0);
 347           } catch (...) {
 348               cout << "EXCEPTION in box = new TBRIK" << endl;
 349               return;
 350           }
 351           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
 352           gm->GetRotMatrix(lay,lad,det,rt);
 353           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
 354           for(i=0;i<9;i++) rtd[i] = rt[i];
 355           try {
 356                 rm  = new TRotMatrix(name,name,rtd);
 357           } catch (...) {
 358                 cout << "EXCEPTION in   new TRotMatrix" << endl;
 359                 return;
 360           }
 361          top->cd();
 362           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
 363          try {
 364               nd  = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm);
 365          } catch (...) {
 366               cout << "EXCEPTION in new TNode" << endl;
 367               return;
 368          }
 369          nd->SetLineColor(kColorITSSSD);
 370          fNodes->Add(nd);
 371     }
 372
 373  for(lay=6;lay<=6;lay++)
 374    for(lad=1;lad<=gm->GetNladders(lay);lad++)
 375     for(det=1;det<=gm->GetNdetectors(lay);det++){
 376           try {
 377               box  = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET",
 378                                     3.65,0.015,2.0);
 379           } catch (...) {
 380               cout << "EXCEPTION in box = new TBRIK" << endl;
 381               return;
 382           }
 383
 384           gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
 385           gm->GetRotMatrix(lay,lad,det,rt);
 386           //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
 387           for(i=0;i<9;i++) rtd[i] = rt[i];
 388           try {
 389                 rm  = new TRotMatrix(name,name,rtd);
 390           } catch (...) {
 391                 cout << "EXCEPTION in   new TRotMatrix" << endl;
 392                 return;
 393           }
 394          top->cd();
 395           //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
 396          try {
 397               nd  = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm);
 398          } catch (...) {
 399               cout << "EXCEPTION in new TNode" << endl;
 400               return;
 401          }
 402          nd->SetLineColor(kColorITSSSD);
 403          fNodes->Add(nd);
 404     }
 405
 406
 407 }
 408 //_____________________________________________________________________________
 409 void AliITSv5::CreateMaterials(){
 410 ////////////////////////////////////////////////////////////////////////
 411 //     Read a file containing the materials for the ITS version 5.
 412 ////////////////////////////////////////////////////////////////////////
 413     char *filtmp;
 414
 415   filtmp = gSystem->ExpandPathName(fEuclidMaterial.Data());
 416
 417   FILE *file = fopen(filtmp,"r");
 418   if(file) {
 419     fclose(file);
 420     ReadEuclidMedia(filtmp);
 421   } else {
 422     Error("CreateMaterials"," THE MEDIA FILE %s DOES NOT EXIST !",filtmp);
 423     exit(1);
 424   } // end if(file)
 425 }
 426 //_____________________________________________________________________________
 427 void AliITSv5::CreateGeometry(){
 428 //////////////////////////////////////////////////////////////////////
 429 //    This is the geometry used for the ITS Pre-TDR and comes from an
 430 // Euclid to Geant conversion. The only difference
 431 // is in the details of the ITS supports. The detectors elements,
 432 // detector numbering, and local and global reference frames are shown in
 433 //  the following figures.
 434 //Begin_Html
 435 /*
 436 <img src="picts/ITS/its1+2_convention_front_5.gif">
 437 </pre>
 438 <br clear=left>
 439 <font size=+2 color=red>
 440 <p>This shows the front view of the SPDs.
 441 </font>
 442 <pre>
 443 <img src="picts/ITS/its1+2_convention_side_5.gif">
 444 </pre>
 445 <br clear=left>
 446 <font size=+2 color=red>
 447 <p>This shows the perspective view of the SPDs.
 448 </font>
 449 <img src="picts/ITS/its1+2_tree.gif">
 450 </pre>
 451 <br clear=left>
 452 <font size=+2 color=red>
 453 <p>This shows the geometry Tree for the SPDs.
 454 </font>
 455 <pre>
 456
 457 <pre>
 458 <img src="picts/ITS/its3+4_convention_front_5.gif">
 459 </pre>
 460 <br clear=left>
 461 <font size=+2 color=red>
 462 <p>This shows the front view of the SDDs.
 463 </font>
 464 <pre>
 465 <img src="picts/ITS/its3+4_convention_side_5.gif">
 466 </pre>
 467 <br clear=left>
 468 <font size=+2 color=red>
 469 <p>This shows the perspective view of the SDDs.
 470 </font>
 471 <img src="picts/ITS/its3+4_tree.gif">
 472 </pre>
 473 <br clear=left>
 474 <font size=+2 color=red>
 475 <p>This shows the geometry Tree for the SDDs.
 476 </font>
 477 <pre>
 478
 479 <pre>
 480 <img src="picts/ITS/its5+6_convention_front_5.gif">
 481 </pre>
 482 <br clear=left>
 483 <font size=+2 color=red>
 484 <p>This shows the front view of the SSDs.
 485 </font>
 486 <pre>
 487 <img src="picts/ITS/its5+6_convention_side_5.gif">
 488 </pre>
 489 <br clear=left>
 490 <font size=+2 color=red>
 491 <p>This shows the perspective view of the SSDs.
 492 </font>
 493 <pre>
 494 <img src="picts/ITS/its5+6_tree.gif">
 495 </pre>
 496 <br clear=left>
 497 <font size=+2 color=red>
 498 <p>This shows the geometry Tree for the SSDs.
 499 </font>
 500 <pre>
 501
 502
 503 <img src="picts/ITS/its_layer1-6_2.gif">
 504 </pre>
 505 <br clear=left>
 506 <font size=+2 color=red>
 507 <p>This shows the front view of the whole ITS..
 508 </font>
 509 <pre>
 510
 511 <img src="picts/ITS/its_layer1-6_1.gif">
 512 </pre>
 513 <br clear=left>
 514 <font size=+2 color=red>
 515 <p>This shows the perspective view of the whole ITS..
 516 </font>
 517 <pre>
 518
 519 <img src="picts/ITS/its1-6_tree.gif">
 520 </pre>
 521 <br clear=left>
 522 <font size=+2 color=red>
 523 <p>This shows the geometry Tree for the whole ITS.
 524 </font>
 525 <pre>
 526 */
 527 //End_Html
 528 //
 529 //
 530 //      Here are shown the details of the ITS support cones and services.
 531 // First is a GEANT tree showing the organization of all of the volumes
 532 // that make up the ITS supports and services.
 533 //Begin_Html
 534 /*
 535 <img src="picts/ITS/supports_tree.gif">
 536  */
 537 //End_Html
 538 //     What follows are a number of figures showing what these support
 539 // structures look like.
 540 //Begin_Html
 541 /*
 542
 543 <img src="picts/ITS/supports_3.gif">
 544 </pre>
 545 <br clear=left>
 546 <font size=+2 color=red>
 547 <p>This shows the geometry of the supports for the Drift and Strip layers only.
 548 </font>
 549 <pre>
 550
 551 <img src="picts/ITS/supports_2.gif">
 552 </pre>
 553 <br clear=left>
 554 <font size=+2 color=red>
 555 <p>This shows the geometry of the supports for the Drift and Strip layers in front cut out.
 556 </font>
 557 <pre>
 558
 559 <img src="picts/ITS/supports_1.gif">
 560 </pre>
 561 <br clear=left>
 562 <font size=+2 color=red>
 563 <p>This shows the geometry of the supports for the Drift and Strip layers in a back cut out..
 564 </font>
 565 <pre>
 566
 567 <img src="picts/ITS/suppssd.gif">
 568 </pre>
 569 <br clear=left>
 570 <font size=+2 color=red>
 571 <p>This shows the geometry for the Strip layers supports.
 572 </font>
 573 <pre>
 574
 575 <img src="picts/ITS/suppsdd.gif">
 576 </pre>
 577 <br clear=left>
 578 <font size=+2 color=red>
 579 <p>This shows the geometry for the Drift layers supports.
 580 </font>
 581 <pre>
 582  */
 583 //End_Html
 584 //
 585 //    Read a file containing the geometry for the ITS version 5.
 586 ////////////////////////////////////////////////////////////////////////
 587
 588
 589     char topvol[5];
 590     char *filtmp;
 591
 592   filtmp = gSystem->ExpandPathName(fEuclidGeometry.Data());
 593   FILE *file = fopen(filtmp,"r");
 594   delete [] filtmp;
 595   if(file) {
 596     fclose(file);
 597     cout << "Ready to read Euclid geometry file" << endl;
 598     ReadEuclid(fEuclidGeometry.Data(),topvol);
 599     printf("Read in euclid geometries\n");
 600   } else {
 601     Error("CreateGeometry"," THE GEOM FILE %s DOES NOT EXIST !",
 602           fEuclidGeometry.Data());
 603     exit(1);
 604   } // end if(file)
 605   //
 606   // Place the ITS ghost volume ITSV in its mother volume (ALIC) and make it
 607   // invisible
 608   //
 609   gMC->Gspos("ITSV",1,"ALIC",0,0,0,0,"ONLY");
 610   //
 611   // Outputs the geometry tree in the EUCLID/CAD format if requested to do so
 612
 613     if (fEuclidOut) {
 614       gMC->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
 615     } // end if (fEuclidOut)
 616
 617     cout << "finished with euclid geometrys" << endl;
 618 }
 619 //______________________________________________________________________
 620 void AliITSv5::ReadOldGeometry(const char *filename){
 621     // read in the file containing the transformations for the active
 622     // volumes for the ITS version 5. This is expected to be in a file
 623     // ending in .det. This geometry is kept in the AliITSgeom class.
 624     Int_t size;
 625     char *filtmp;
 626     FILE *file;
 627
 628     if(fITSgeom!=0) delete fITSgeom;
 629     filtmp = gSystem->ExpandPathName(filename);
 630     size = strlen(filtmp);
 631     if(size>4 && fGeomDetIn){
 632         filtmp[size-3] = 'd'; // change from .euc to .det
 633         filtmp[size-2] = 'e';
 634         filtmp[size-1] = 't';
 635         file = fopen(filtmp,"r");
 636         if(file){ // if file exists use it to fill AliITSgeom structure.
 637             fclose(file);
 638             fITSgeom = new AliITSgeom(filtmp);
 639             fITSgeom->DefineShapes(3); // if fShape isn't defined define it.
 640             // Now define the detector types/shapes.
 641             fITSgeom->ReSetShape(kSPD,new AliITSgeomSPD300());
 642             fITSgeom->ReSetShape(kSDD,new AliITSgeomSDD300());
 643             fITSgeom->ReSetShape(kSSD,new AliITSgeomSSD175());
 644         }else{
 645             fITSgeom = 0;
 646             // fill AliITSgeom structure from geant structure just filled above
 647         }// end if(file)
 648         delete [] filtmp;
 649     }// end if(size>4)
 650 }
 651 //______________________________________________________________________
 652 void AliITSv5::InitAliITSgeom(){
 653 //     Based on the geometry tree defined in Geant 3.21, this
 654 // routine initilizes the Class AliITSgeom from the Geant 3.21 ITS geometry
 655 // sturture.
 656 //    if(gMC->IsA()!=TGeant3::Class()) {
 657     if(strcmp(gMC->GetName(),"TGeant3")) {
 658         Error("InitAliITSgeom",
 659                 "Wrong Monte Carlo. InitAliITSgeom uses TGeant3 calls");
 660         return;
 661     } // end if
 662     cout << "Reading Geometry transformation directly from Geant 3." << endl;
 663     const Int_t nlayers = 6;
 664     const Int_t ndeep = 7;
 665     Int_t itsGeomTreeNames[nlayers][ndeep],lnam[20],lnum[20];
 666     Int_t nlad[nlayers],ndet[nlayers];
 667     Double_t t[3],r[10];
 668     Float_t  par[20],att[20];
 669     Int_t    npar,natt,idshape,imat,imed;
 670     AliITSGeant3Geometry *ig = new AliITSGeant3Geometry();
 671     Int_t mod,lay,lad,det,i,j,k;
 672     char *names[nlayers][ndeep] = {
 673         {"ALIC","ITSV","ITSD","IT12","I132","I186","ITS1"}, // lay=1
 674         {"ALIC","ITSV","ITSD","IT12","I132","I131","ITS2"}, // lay=2
 675         {"ALIC","ITSV","ITSD","IT34","I004","I302","ITS3"}, // lay=3
 676         {"ALIC","ITSV","ITSD","IT34","I005","I402","ITS4"}, // lay=4
 677         {"ALIC","ITSV","ITSD","IT56","I565","I562","ITS5"}, // lay=5
 678         {"ALIC","ITSV","ITSD","IT56","I569","I566","ITS6"}};// lay=6
 679     Int_t itsGeomTreeCopys[nlayers][ndeep] = {{1,1,1,1,10, 2,4}, // lay=1
 680                                               {1,1,1,1,10, 4,4}, // lay=2
 681                                               {1,1,1,1,14, 6,1}, // lay=3
 682                                               {1,1,1,1,22, 8,1}, // lay=4
 683                                               {1,1,1,1,34,23,1}, // lay=5
 684                                               {1,1,1,1,38,26,1}};// lay=6
 685
 686     // Sorry, but this is not very pritty code. It should be replaced
 687     // at some point with a version that can search through the geometry
 688     // tree its self.
 689     for(i=0;i<20;i++) lnam[i] = lnum[i] = 0;
 690     for(i=0;i<nlayers;i++)for(j=0;j<ndeep;j++)
 691         itsGeomTreeNames[i][j] = ig->StringToInt(names[i][j]);
 692     mod = 0;
 693     for(i=0;i<nlayers;i++){
 694         k = 1;
 695         for(j=0;j<ndeep;j++) if(itsGeomTreeCopys[i][j]!=0)
 696             k *= TMath::Abs(itsGeomTreeCopys[i][j]);
 697         mod += k;
 698     } // end for i
 699
 700     if(fITSgeom!=0) delete fITSgeom;
 701     nlad[0]=20;nlad[1]=40;nlad[2]=14;nlad[3]=22;nlad[4]=34;nlad[5]=38;
 702     ndet[0]=4;ndet[1]=4;ndet[2]=6;ndet[3]=8;ndet[4]=22;ndet[5]=25;
 703     fITSgeom = new AliITSgeom(0,6,nlad,ndet,mod);
 704     mod = -1;
 705     for(lay=1;lay<=nlayers;lay++){
 706         for(j=0;j<ndeep;j++) lnam[j] = itsGeomTreeNames[lay-1][j];
 707         for(j=0;j<ndeep;j++) lnum[j] = itsGeomTreeCopys[lay-1][j];
 708         switch (lay){
 709         case 1: case 2: // layers 1 and 2 are a bit special
 710             lad = 0;
 711             for(j=1;j<=itsGeomTreeCopys[lay-1][4];j++){
 712                 lnum[4] = j;
 713                 for(k=1;k<=itsGeomTreeCopys[lay-1][5];k++){
 714                     lad++;
 715                     lnum[5] = k;
 716                     for(det=1;det<=itsGeomTreeCopys[lay-1][6];det++){
 717                         lnum[6] = det;
 718                         mod++;
 719                         ig->GetGeometry(ndeep,lnam,lnum,t,r,idshape,npar,natt,
 720                                         par,att,imat,imed);
 721                         fITSgeom->CreatMatrix(mod,lay,lad,det,kSPD,t,r);
 722                         if(!(fITSgeom->IsShapeDefined((Int_t)kSPD)))
 723                              fITSgeom->ReSetShape(kSPD,
 724                                                   new AliITSgeomSPD300());
 725                     } // end for det
 726                 } // end for k
 727             } // end for j
 728             break;
 729         case 3: case 4: case 5: case 6: // layers 3-6
 730             lnum[6] = 1;
 731             for(lad=1;lad<=itsGeomTreeCopys[lay-1][4];lad++){
 732                 lnum[4] = lad;
 733                 for(det=1;det<=itsGeomTreeCopys[lay-1][5];det++){
 734                     lnum[5] = det;
 735                     mod++;
 736                     ig->GetGeometry(7,lnam,lnum,t,r,idshape,npar,natt,
 737                                     par,att,imat,imed);
 738                     switch (lay){
 739                     case 3: case 4:
 740                         fITSgeom->CreatMatrix(mod,lay,lad,det,kSDD,t,r);
 741                         if(!(fITSgeom->IsShapeDefined(kSDD)))
 742                             fITSgeom->ReSetShape(kSDD,new AliITSgeomSDD300());
 743                         break;
 744                     case 5: case 6:
 745                         fITSgeom->CreatMatrix(mod,lay,lad,det,kSSD,t,r);
 746                         if(!(fITSgeom->IsShapeDefined(kSSD)))
 747                             fITSgeom->ReSetShape(kSSD,new AliITSgeomSSD175());
 748                         break;
 749                         } // end switch
 750                 } // end for det
 751             } // end for lad
 752             break;
 753         } // end switch
 754     } // end for lay
 755     return;
 756 }
 757 //_____________________________________________________________________________
 758 void AliITSv5::Init(){
 759 ////////////////////////////////////////////////////////////////////////
 760 //     Initialise the ITS after it has been created.
 761 ////////////////////////////////////////////////////////////////////////
 762     Int_t i;
 763
 764     cout << endl;
 765     for(i=0;i<30;i++) cout << "*";cout << " ITSv5_Init ";
 766     for(i=0;i<30;i++) cout << "*";cout << endl;
 767 //
 768     if(fRead[0]=='\0') strncpy(fRead,fEuclidGeomDet,60);
 769     if(fWrite[0]=='\0') strncpy(fWrite,fEuclidGeomDet,60);
 770     if(fGeomDetIn && !fGeomOldDetIn){
 771         if(fITSgeom!=0) delete fITSgeom;
 772         fITSgeom = new AliITSgeom();
 773         fITSgeom->ReadNewFile(fRead);
 774     } // end if
 775     if(fGeomDetIn &&  fGeomOldDetIn) ReadOldGeometry(fEuclidGeometry.Data());
 776
 777     if(!fGeomDetIn) this->InitAliITSgeom();
 778     if(fGeomDetOut) fITSgeom->WriteNewFile(fWrite);
 779     AliITS::Init();
 780 //
 781     for(i=0;i<72;i++) cout << "*";
 782     cout << endl;
 783 }
 784 //_____________________________________________________________________________
 785 void AliITSv5::StepManager(){
 786 ////////////////////////////////////////////////////////////////////////
 787 //    Called for every step in the ITS, then calles the AliITShit class
 788 // creator with the information to be recoreded about that hit.
 789 //     The value of the macro ALIITSPRINTGEOM if set to 1 will allow the
 790 // printing of information to a file which can be used to create a .det
 791 // file read in by the routine CreateGeometry(). If set to 0 or any other
 792 // value except 1, the default behavior, then no such file is created nor
 793 // it the extra variables and the like used in the printing allocated.
 794 ////////////////////////////////////////////////////////////////////////
 795   Int_t          copy, id;
 796   Int_t          copy1,copy2;
 797   Float_t        hits[8];
 798   Int_t          vol[4];
 799   TLorentzVector position, momentum;
 800   TClonesArray   &lhits = *fHits;
 801   //
 802   // Track status
 803   vol[3] = 0;
 804   if(gMC->IsTrackInside())      vol[3] +=  1;
 805   if(gMC->IsTrackEntering())    vol[3] +=  2;
 806   if(gMC->IsTrackExiting())     vol[3] +=  4;
 807   if(gMC->IsTrackOut())         vol[3] +=  8;
 808   if(gMC->IsTrackDisappeared()) vol[3] += 16;
 809   if(gMC->IsTrackStop())        vol[3] += 32;
 810   if(gMC->IsTrackAlive())       vol[3] += 64;
 811   //
 812   // Fill hit structure.
 813   if(!(gMC->TrackCharge())) return;
 814   //
 815   // Only entering charged tracks
 816   if((id = gMC->CurrentVolID(copy)) == fIdSens[0]) {
 817       vol[0] = 1;
 818       id = gMC->CurrentVolOffID(0,copy);
 819       //detector copy in the ladder = 1<->4  (ITS1)
 820       vol[1] = copy;
 821       gMC->CurrentVolOffID(1,copy1);
 822       //ladder copy in the module   = 1<->2  (I186)
 823       gMC->CurrentVolOffID(2,copy2);
 824       //module copy in the layer    = 1<->10 (I132)
 825       vol[2] = copy1+(copy2-1)*2;//# of ladders in one module  = 2
 826   } else if(id == fIdSens[1]){
 827       vol[0] = 2;
 828       id = gMC->CurrentVolOffID(0,copy);
 829       //detector copy in the ladder = 1<->4  (ITS2)
 830       vol[1] = copy;
 831       gMC->CurrentVolOffID(1,copy1);
 832       //ladder copy in the module   = 1<->4  (I131)
 833       gMC->CurrentVolOffID(2,copy2);
 834       //module copy in the layer    = 1<->10 (I132)
 835       vol[2] = copy1+(copy2-1)*4;//# of ladders in one module  = 4
 836   } else if(id == fIdSens[2]){
 837       vol[0] = 3;
 838       id = gMC->CurrentVolOffID(1,copy);
 839       //detector copy in the ladder = 1<->5  (ITS3 is inside I314)
 840       vol[1] = copy;
 841       id = gMC->CurrentVolOffID(2,copy);
 842       //ladder copy in the layer    = 1<->12 (I316)
 843       vol[2] = copy;
 844   } else if(id == fIdSens[3]){
 845       vol[0] = 4;
 846       id = gMC->CurrentVolOffID(1,copy);
 847       //detector copy in the ladder = 1<->8  (ITS4 is inside I414)
 848       vol[1] = copy;
 849       id = gMC->CurrentVolOffID(2,copy);
 850       //ladder copy in the layer    = 1<->22 (I417)
 851       vol[2] = copy;
 852   }else if(id == fIdSens[4]){
 853       vol[0] = 5;
 854       id = gMC->CurrentVolOffID(1,copy);
 855       //detector copy in the ladder = 1<->23  (ITS5 is inside I562)
 856       vol[1] = copy;
 857       id = gMC->CurrentVolOffID(2,copy);
 858      //ladder copy in the layer    = 1<->34 (I565)
 859       vol[2] = copy;
 860   }else if(id == fIdSens[5]){
 861       vol[0] = 6;
 862       id = gMC->CurrentVolOffID(1,copy);
 863       //detector copy in the ladder = 1<->26  (ITS6 is inside I566)
 864       vol[1] = copy;
 865       id = gMC->CurrentVolOffID(2,copy);
 866       //ladder copy in the layer = 1<->38 (I569)
 867       vol[2] = copy;
 868   } else {
 869       return; // not an ITS volume?
 870   } // end if/else if (gMC->CurentVolID(copy) == fIdSens[i])
 871 //
 872   gMC->TrackPosition(position);
 873   gMC->TrackMomentum(momentum);
 874   hits[0]=position[0];
 875   hits[1]=position[1];
 876   hits[2]=position[2];
 877   hits[3]=momentum[0];
 878   hits[4]=momentum[1];
 879   hits[5]=momentum[2];
 880   hits[6]=gMC->Edep();
 881   hits[7]=gMC->TrackTime();
 882   // Fill hit structure with this new hit.
 883   new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
 884   return;
 885 }