From d683f3d57a07c8a59d76fffd38458f25f07e4224 Mon Sep 17 00:00:00 2001 From: barbera Date: Sat, 7 Oct 2000 15:46:38 +0000 Subject: [PATCH] Version 5 of the geometry with asymmetric services --- ITS/AliITSv5asymm.cxx | 779 ++++++++++++++++++++++++++++++++++++++++++ ITS/AliITSv5asymm.h | 39 +++ 2 files changed, 818 insertions(+) create mode 100644 ITS/AliITSv5asymm.cxx create mode 100644 ITS/AliITSv5asymm.h diff --git a/ITS/AliITSv5asymm.cxx b/ITS/AliITSv5asymm.cxx new file mode 100644 index 00000000000..2c7fe339731 --- /dev/null +++ b/ITS/AliITSv5asymm.cxx @@ -0,0 +1,779 @@ +/************************************************************************** + * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * + * * + * Author: The ALICE Off-line Project. * + * Contributors are mentioned in the code where appropriate. * + * * + * Permission to use, copy, modify and distribute this software and its * + * documentation strictly for non-commercial purposes is hereby granted * + * without fee, provided that the above copyright notice appears in all * + * copies and that both the copyright notice and this permission notice * + * appear in the supporting documentation. The authors make no claims * + * about the suitability of this software for any purpose. It is * + * provided "as is" without express or implied warranty. * + **************************************************************************/ + +/* +$Log$ +Revision 1.25 2000/10/05 20:50:00 nilsen +Now using root generated streamers. + +Revision 1.14.4.12 2000/10/02 16:04:03 barbera +Forward declarations added + +Revision 1.22 2000/07/10 16:07:19 fca +Release version of ITS code + +Revision 1.14.4.4 2000/05/19 10:10:21 nilsen +fix for bug with HP and Sun unix + fix for event display in ITS-working branch + +Revision 1.14.4.3 2000/03/04 23:46:38 nilsen +Fixed up the comments/documentation. + +Revision 1.14.4.2 2000/03/02 21:53:02 nilsen +To make it compatable with the changes in AliRun/AliModule. + +Revision 1.14.4.1 2000/01/12 19:03:33 nilsen +This is the version of the files after the merging done in December 1999. +See the ReadMe110100.txt file for details + +Revision 1.14 1999/10/22 08:16:49 fca +Correct destructors, thanks to I.Hrivnacova + +Revision 1.13 1999/10/06 10:15:19 fca +Correct bug in allocation of layer name and add destructor + +Revision 1.12 1999/10/05 08:05:09 fca +Minor corrections for uninitialised variables. + +Revision 1.11 1999/09/29 09:24:20 fca +Introduction of the Copyright and cvs Log + +*/ + +/////////////////////////////////////////////////////////////////////////////// +// +// Inner Traking System version 5 with asymmetric services +// This class contains the base procedures for the Inner Tracking System +// +// Authors: R. Barbera, B. S. Nilsen. +// version 5. +// Created October 7 2000. +// +/////////////////////////////////////////////////////////////////////////////// + +// See AliITSv5asymm::StepManager(). +#define ALIITSPRINTGEOM 0 // default. don't print out gemetry information +//#define ALIITSPRINTGEOM 1 // print out geometry information + +#include +#include +#include +#include +#include +#include +#include // only required for Tracking function? +#include +#include +#include +#include +#include +#include + +#include "AliMC.h" +#include "AliRun.h" +#if ALIITSPRINTGEOM==1 +#include "../TGeant3/TGeant3.h" +#endif +#include "AliITShit.h" +#include "AliITS.h" +#include "AliITSv5asymm.h" +#include "AliITSgeom.h" + +ClassImp(AliITSv5asymm) + +//_____________________________________________________________________________ +AliITSv5asymm::AliITSv5asymm() { +//////////////////////////////////////////////////////////////////////// +// Standard default constructor for the ITS version 5. +//////////////////////////////////////////////////////////////////////// + + fIdN = 0; + fIdName = 0; + fIdSens = 0; + fMajorVersion = 5; + fMinorVersion = -1; +} +//____________________________________________________________________________ +AliITSv5asymm::AliITSv5asymm(const AliITSv5asymm &source){ +//////////////////////////////////////////////////////////////////////// +// Copy Constructor for ITS version 5. +//////////////////////////////////////////////////////////////////////// + if(&source == this) return; + printf("Not allowed to copy AliITSv5asymm\n"); + return; +} +//_____________________________________________________________________________ +AliITSv5asymm& AliITSv5asymm::operator=(const AliITSv5asymm &source){ +//////////////////////////////////////////////////////////////////////// +// Assignment operator for the ITS version 1. +//////////////////////////////////////////////////////////////////////// + + if(&source == this) return *this; + printf("Not allowed to copy AliITSv3\n"); + return *this; + +} +//_____________________________________________________________________________ +AliITSv5asymm::~AliITSv5asymm() { +//////////////////////////////////////////////////////////////////////// +// Standard destructor for the ITS version 5. +//////////////////////////////////////////////////////////////////////// +} +//_____________________________________________________________________________ +AliITSv5asymm::AliITSv5asymm(const char *name, const char *title) : AliITS(name, title){ +//////////////////////////////////////////////////////////////////////// +// Standard constructor for the ITS version 5. +//////////////////////////////////////////////////////////////////////// + fIdN = 6; +/* +// TObjArray of TObjStrings + fIdName = new TObjArray(fIdN); + fIdName->AddAt(new TObjString("ITS1"),0); + fIdName->AddAt(new TObjString("ITS2"),1); + fIdName->AddAt(new TObjString("ITS3"),2); + fIdName->AddAt(new TObjString("ITS4"),3); + fIdName->AddAt(new TObjString("ITS5"),4); + fIdName->AddAt(new TObjString("ITS6"),5); +*/ +// Array of TStrings. + fIdName = new TString[fIdN]; + fIdName[0] = "ITS1"; + fIdName[1] = "ITS2"; + fIdName[2] = "ITS3"; + fIdName[3] = "ITS4"; + fIdName[4] = "ITS5"; + fIdName[5] = "ITS6"; + fIdSens = new Int_t[fIdN]; + for (Int_t i=0;iGetGeometry()->GetNode("alice"); + AliITSgeom *gm = this->GetITSgeom(); + + Int_t lay,lad,det,i; + Text_t name[10]; + Float_t xg[3]; + Float_t rt[9]; + Double_t rtd[9]; + TBRIK *box; + TRotMatrix *rm; + //TCanvas *c1 = new TCanvas("c1","ITS"); + + for(lay=1;lay<=2;lay++) + for(lad=1;lad<=gm->GetNladders(lay);lad++) + for(det=1;det<=gm->GetNdetectors(lay);det++){ + try { + box = new TBRIK ("ActiveSPD","Active volume of SPD","SPD SI DET", + 0.64,0.0075,4.19); + } catch (...) { + cout << "EXCEPTION in box = new TBRIK" << endl; + return; + } + gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]); + gm->GetRotMatrix(lay,lad,det,rt); + //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det); + for(i=0;i<9;i++) rtd[i] = rt[i]; + try { + rm = new TRotMatrix(name,name,rtd); + } catch (...) { + cout << "EXCEPTION in new TRotMatrix" << endl; + return; + } + top->cd(); + //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); + try { + nd = new TNode("SPD"," ",box,xg[0],xg[1],xg[2],rm); + } catch (...) { + cout << "EXCEPTION in new TNode" << endl; + return; + } + nd->SetLineColor(kColorITSSSD); + fNodes->Add(nd); + } + + for(lay=3;lay<=3;lay++) + for(lad=1;lad<=gm->GetNladders(lay);lad++) + for(det=1;det<=gm->GetNdetectors(lay);det++){ + try { + box = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET", + 3.5,0.014,3.763); + } catch (...) { + cout << "EXCEPTION in box = new TBRIK" << endl; + return; + } + gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]); + gm->GetRotMatrix(lay,lad,det,rt); + //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det); + for(i=0;i<9;i++) rtd[i] = rt[i]; + try { + rm = new TRotMatrix(name,name,rtd); + } catch (...) { + cout << "EXCEPTION in new TRotMatrix" << endl; + return; + } + top->cd(); + //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); + try { + nd = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm); + } catch (...) { + cout << "EXCEPTION in new TNode" << endl; + return; + } + nd->SetLineColor(kColorITSSSD); + fNodes->Add(nd); + } + + for(lay=4;lay<=4;lay++) + for(lad=1;lad<=gm->GetNladders(lay);lad++) + for(det=1;det<=gm->GetNdetectors(lay);det++){ + try { + box = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET", + 3.5,0.014,3.763); + } catch (...) { + cout << "EXCEPTION in box = new TBRIK" << endl; + return; + } + gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]); + gm->GetRotMatrix(lay,lad,det,rt); + //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det); + for(i=0;i<9;i++) rtd[i] = rt[i]; + try { + rm = new TRotMatrix(name,name,rtd); + } catch (...) { + cout << "EXCEPTION in new TRotMatrix" << endl; + return; + } + top->cd(); + //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); + try { + nd = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm); + } catch (...) { + cout << "EXCEPTION in new TNode" << endl; + return; + } + nd->SetLineColor(kColorITSSSD); + fNodes->Add(nd); + } + for(lay=5;lay<=5;lay++) + for(lad=1;lad<=gm->GetNladders(lay);lad++) + for(det=1;det<=gm->GetNdetectors(lay);det++){ + try { + box = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET", + 3.65,0.015,2.0); + } catch (...) { + cout << "EXCEPTION in box = new TBRIK" << endl; + return; + } + gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]); + gm->GetRotMatrix(lay,lad,det,rt); + //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det); + for(i=0;i<9;i++) rtd[i] = rt[i]; + try { + rm = new TRotMatrix(name,name,rtd); + } catch (...) { + cout << "EXCEPTION in new TRotMatrix" << endl; + return; + } + top->cd(); + //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); + try { + nd = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm); + } catch (...) { + cout << "EXCEPTION in new TNode" << endl; + return; + } + nd->SetLineColor(kColorITSSSD); + fNodes->Add(nd); + } + + for(lay=6;lay<=6;lay++) + for(lad=1;lad<=gm->GetNladders(lay);lad++) + for(det=1;det<=gm->GetNdetectors(lay);det++){ + try { + box = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET", + 3.65,0.015,2.0); + } catch (...) { + cout << "EXCEPTION in box = new TBRIK" << endl; + return; + } + + gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]); + gm->GetRotMatrix(lay,lad,det,rt); + //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det); + for(i=0;i<9;i++) rtd[i] = rt[i]; + try { + rm = new TRotMatrix(name,name,rtd); + } catch (...) { + cout << "EXCEPTION in new TRotMatrix" << endl; + return; + } + top->cd(); + //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det); + try { + nd = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm); + } catch (...) { + cout << "EXCEPTION in new TNode" << endl; + return; + } + nd->SetLineColor(kColorITSSSD); + fNodes->Add(nd); + } + + +} +//_____________________________________________________________________________ +void AliITSv5asymm::CreateMaterials(){ +//////////////////////////////////////////////////////////////////////// +// Read a file containing the materials for the ITS version 5. +//////////////////////////////////////////////////////////////////////// + char *filtmp; + + filtmp = gSystem->ExpandPathName(fEuclidMaterial.Data()); + + FILE *file = fopen(filtmp,"r"); + if(file) { + fclose(file); + ReadEuclidMedia(filtmp); + } else { + Error("CreateMaterials"," THE MEDIA FILE %s DOES NOT EXIST !",filtmp); + exit(1); + } // end if(file) +} +//_____________________________________________________________________________ +void AliITSv5asymm::CreateGeometry(){ +////////////////////////////////////////////////////////////////////// +// This is the geometry used for the ITS Pre-TDR and comes from an +// Euclid to Geant conversion. The only difference +// is in the details of the ITS supports. The detectors elements, +// detector numbering, and local and global reference frames are shown in +// the following figures. +//Begin_Html +/* + + +
+ +

This shows the front view of the SPDs. + +

+
+
+
+ +

This shows the perspective view of the SPDs. + + + +
+ +

This shows the geometry Tree for the SPDs. + +

+
+
+
+

+

+
+
This shows the front view of the SDDs.
+
+
+
+

+

+
+
This shows the perspective view of the SDDs.
+
+
+
+
+ +

This shows the geometry Tree for the SDDs. + +

+
+
+
+

+

+
+
This shows the front view of the SSDs.
+
+
+
+

+

+
+
This shows the perspective view of the SSDs.
+
+
+
+

+

+
+
This shows the geometry Tree for the SSDs.
+
+
+
+
+
+

+

+
+
This shows the front view of the whole ITS..
+
+
+
+
+

+

+
+
This shows the perspective view of the whole ITS..
+
+
+
+
+

+

+
+
This shows the geometry Tree for the whole ITS.
+
+
+*/
+//End_Html
+//
+//
+//      Here are shown the details of the ITS support cones and services.
+// First is a GEANT tree showing the organization of all of the volumes
+// that make up the ITS supports and services.
+//Begin_Html
+/*
+
+ */
+//End_Html
+//     What follows are a number of figures showing what these support
+// structures look like.
+//Begin_Html
+/*
+
+
+

+

+
+
This shows the geometry of the supports for the Drift and Strip layers only.
+
+
+
+
+

+

+
+
This shows the geometry of the supports for the Drift and Strip layers in front cut out.
+
+
+
+
+

+

+
+
This shows the geometry of the supports for the Drift and Strip layers in a back cut out..
+
+
+
+
+

+

+
+
This shows the geometry for the Strip layers supports.
+
+
+
+
+

+

+
+
This shows the geometry for the Drift layers supports.
+
+
+ */
+//End_Html
+//
+//    Read a file containing the geometry for the ITS version 5.
+////////////////////////////////////////////////////////////////////////
+
+    Int_t size;
+    char topvol[5];
+    char *filtmp;
+
+  filtmp = gSystem->ExpandPathName(fEuclidGeometry.Data());
+  FILE *file = fopen(filtmp,"r");
+  delete [] filtmp;
+  if(file) {
+    fclose(file);
+    printf("Ready to read Euclid geometry file\n");
+    ReadEuclid(fEuclidGeometry.Data(),topvol);
+    printf("Read in euclid geometries\n");
+  } else {
+    Error("CreateGeometry"," THE GEOM FILE %s DOES NOT EXIST !",
+	  fEuclidGeometry.Data());
+    exit(1);
+  } // end if(file)
+  //
+  // Place the ITS ghost volume ITSV in its mother volume (ALIC) and make it
+  // invisible
+  //
+  gMC->Gspos("ITSV",1,"ALIC",0,0,0,0,"ONLY");
+  //
+  // Outputs the geometry tree in the EUCLID/CAD format if requested to do so
+  
+    if (fEuclidOut) {
+      gMC->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
+    } // end if (fEuclidOut)
+
+    // read in the file containing the transformations for the active
+    // volumes for the ITS version 5. This is expected to be in a file
+    // ending in .det. This geometry is kept in the AliITSgeom class.
+    filtmp = gSystem->ExpandPathName(fEuclidGeometry.Data());
+    size = strlen(filtmp);
+    if(size>4){
+	filtmp[size-3] = 'd'; // change from .euc to .det
+        filtmp[size-2] = 'e';
+        filtmp[size-1] = 't';
+	file = fopen(filtmp,"r");
+	if(file){ // if file exists use it to fill AliITSgeom structure.
+	    fclose(file);
+	    printf("ready to read .det file %s\n",filtmp);
+	    fITSgeom = new AliITSgeom(filtmp);
+	}else{
+	    fITSgeom = 0;
+	    // fill AliITSgeom structure from geant structure just filled above
+	}// end if(file)
+        delete [] filtmp;
+    }// end if(size>4)
+    printf("finished with euclid geometrys\n");
+}
+//_____________________________________________________________________________
+void AliITSv5asymm::Init(){
+////////////////////////////////////////////////////////////////////////
+//     Initialise the ITS after it has been created.
+////////////////////////////////////////////////////////////////////////
+
+    AliITS::Init();
+    fMajorVersion = 5;
+    fMinorVersion = 0;
+} 
+//_____________________________________________________________________________
+void AliITSv5asymm::StepManager(){
+////////////////////////////////////////////////////////////////////////
+//    Called for every step in the ITS, then calles the AliITShit class
+// creator with the information to be recoreded about that hit.
+//     The value of the macro ALIITSPRINTGEOM if set to 1 will allow the
+// printing of information to a file which can be used to create a .det
+// file read in by the routine CreateGeometry(). If set to 0 or any other
+// value except 1, the default behavior, then no such file is created nor
+// it the extra variables and the like used in the printing allocated.
+////////////////////////////////////////////////////////////////////////
+  Int_t          copy, id;
+  Int_t          copy1,copy2;
+  Float_t        hits[8];
+  Int_t          vol[4];
+  TLorentzVector position, momentum;
+  TClonesArray   &lhits = *fHits;
+#if ALIITSPRINTGEOM==1
+  FILE          *fp;
+  Int_t         i;
+  Float_t       xl[3],xt[3],angl[6];
+//  Float_t       par[20],att[20];
+  Float_t      mat[9];
+  static Bool_t first=kTRUE,printit[6][50][50];
+  if(first){ for(copy1=0;copy1<6;copy1++)for(copy2=0;copy2<50;copy2++)
+      for(id=0;id<50;id++) printit[copy1][copy2][id] = kTRUE;
+  first = kFALSE;
+  }
+  // end if first
+#endif
+  //
+  // Track status
+  vol[3] = 0;
+  if(gMC->IsTrackInside())      vol[3] +=  1;
+  if(gMC->IsTrackEntering())    vol[3] +=  2;
+  if(gMC->IsTrackExiting())     vol[3] +=  4;
+  if(gMC->IsTrackOut())         vol[3] +=  8;
+  if(gMC->IsTrackDisappeared()) vol[3] += 16;
+  if(gMC->IsTrackStop())        vol[3] += 32;
+  if(gMC->IsTrackAlive())       vol[3] += 64;
+  //
+  // Fill hit structure.
+  if(!(gMC->TrackCharge())) return;
+  //
+  // Only entering charged tracks
+  if((id = gMC->CurrentVolID(copy)) == fIdSens[0]) {
+      vol[0] = 1;
+      id = gMC->CurrentVolOffID(0,copy);
+      //detector copy in the ladder = 1<->4  (ITS1)
+      vol[1] = copy;
+      gMC->CurrentVolOffID(1,copy1);
+      //ladder copy in the module   = 1<->2  (I186)
+      gMC->CurrentVolOffID(2,copy2);
+      //module copy in the layer    = 1<->10 (I132)
+      vol[2] = copy1+(copy2-1)*2;//# of ladders in one module  = 2
+  } else if(id == fIdSens[1]){
+      vol[0] = 2;
+      id = gMC->CurrentVolOffID(0,copy);
+      //detector copy in the ladder = 1<->4  (ITS2)
+      vol[1] = copy;
+      gMC->CurrentVolOffID(1,copy1);
+      //ladder copy in the module   = 1<->4  (I131)
+      gMC->CurrentVolOffID(2,copy2);
+      //module copy in the layer    = 1<->10 (I132)
+      vol[2] = copy1+(copy2-1)*4;//# of ladders in one module  = 4
+  } else if(id == fIdSens[2]){
+      vol[0] = 3;
+      id = gMC->CurrentVolOffID(1,copy);
+      //detector copy in the ladder = 1<->5  (ITS3 is inside I314)
+      vol[1] = copy;
+      id = gMC->CurrentVolOffID(2,copy);
+      //ladder copy in the layer    = 1<->12 (I316)
+      vol[2] = copy;
+  } else if(id == fIdSens[3]){
+      vol[0] = 4;
+      id = gMC->CurrentVolOffID(1,copy);
+      //detector copy in the ladder = 1<->8  (ITS4 is inside I414)
+      vol[1] = copy;
+      id = gMC->CurrentVolOffID(2,copy);
+      //ladder copy in the layer    = 1<->22 (I417)
+      vol[2] = copy;
+  }else if(id == fIdSens[4]){
+      vol[0] = 5;
+      id = gMC->CurrentVolOffID(1,copy);
+      //detector copy in the ladder = 1<->23  (ITS5 is inside I562)
+      vol[1] = copy;
+      id = gMC->CurrentVolOffID(2,copy);
+     //ladder copy in the layer    = 1<->34 (I565)
+      vol[2] = copy;
+  }else if(id == fIdSens[5]){
+      vol[0] = 6;
+      id = gMC->CurrentVolOffID(1,copy);
+      //detector copy in the ladder = 1<->26  (ITS6 is inside I566)
+      vol[1] = copy;
+      id = gMC->CurrentVolOffID(2,copy);
+      //ladder copy in the layer = 1<->38 (I569)
+      vol[2] = copy;
+  } else {
+      return; // not an ITS volume?
+  } // end if/else if (gMC->CurentVolID(copy) == fIdSens[i])
+//
+  gMC->TrackPosition(position);
+  gMC->TrackMomentum(momentum);
+  hits[0]=position[0];
+  hits[1]=position[1];
+  hits[2]=position[2];
+  hits[3]=momentum[0];
+  hits[4]=momentum[1];
+  hits[5]=momentum[2];
+  hits[6]=gMC->Edep();
+  hits[7]=gMC->TrackTime();
+  // Fill hit structure with this new hit.
+  new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
+#if ALIITSPRINTGEOM==1
+  if(printit[vol[0]][vol[2]][vol[1]]){
+      printit[vol[0]][vol[2]][vol[1]] = kFALSE;
+      xl[0] = xl[1] = xl[2] = 0.0;
+      gMC->Gdtom(xl,xt,1);
+      for(i=0;i<9;i++) mat[i] = 0.0;
+      mat[0] = mat[4] = mat[8] = 1.0;  // default with identity matrix
+      xl[0] = 1.0;
+      xl[1] = xl[2] =0.0;
+      gMC->Gdtom(xl,&(mat[0]),2);
+      xl[1] = 1.0;
+      xl[0] = xl[2] =0.0;
+      gMC->Gdtom(xl,&(mat[3]),2);
+      xl[2] = 1.0;
+      xl[1] = xl[0] =0.0;
+      gMC->Gdtom(xl,&(mat[6]),2);
+
+      angl[0] = TMath::ACos(mat[2]);
+      if(mat[2]==1.0) angl[0] = 0.0;
+      angl[1] = TMath::ATan2(mat[1],mat[0]);
+      if(angl[1]<0.0) angl[1] += 2.0*TMath::Pi();
+
+      angl[2] = TMath::ACos(mat[5]);
+      if(mat[5]==1.0) angl[2] = 0.0;
+      angl[3] = TMath::ATan2(mat[4],mat[3]);
+      if(angl[3]<0.0) angl[3] += 2.0*TMath::Pi();
+
+      angl[4] = TMath::ACos(mat[8]);
+      if(mat[8]==1.0) angl[4] = 0.0;
+      angl[5] = TMath::ATan2(mat[7],mat[6]);
+      if(angl[5]<0.0) angl[5] += 2.0*TMath::Pi();
+
+      for(i=0;i<6;i++) angl[i] *= 180.0/TMath::Pi(); // degrees
+      fp = fopen("ITSgeometry_v5.det","a");
+      fprintf(fp,"%2d %2d %2d %9e %9e %9e %9e %9e %9e %9e %9e %9e ",
+	      vol[0],vol[2],vol[1], // layer ladder detector
+	      xt[0],xt[1],xt[2],    // Translation vector
+	      angl[0],angl[1],angl[2],angl[3],angl[4],angl[5] // Geant rotaion
+                                                           // angles (degrees)
+	      );
+      fprintf(fp,"%9e %9e %9e %9e %9e %9e %9e %9e %9e",
+	     mat[0],mat[1],mat[2],mat[3],mat[4],mat[5],mat[6],mat[7],mat[8]
+	  );  // Adding the rotation matrix.
+      fprintf(fp,"\n");
+      fclose(fp);
+  } // end if printit[layer][ladder][detector]
+#endif
+  return;
+}
+/*
+//____________________________________________________________________________
+void AliITSv5asymm::Streamer(TBuffer &R__b){
+////////////////////////////////////////////////////////////////////////
+//    A dummy Streamer function for this class AliITSv5asymm. By default it
+// only streams the AliITS class as it is required. Since this class
+// dosen't contain any "real" data to be saved, it doesn't.
+////////////////////////////////////////////////////////////////////////
+
+   if (R__b.IsReading()) {
+      Version_t R__v = R__b.ReadVersion(); 
+      if (R__v==1) {
+	  AliITS::Streamer(R__b);
+      }else{
+	  AliITS::Streamer(R__b);
+      } // end if R__v==1
+   } else {
+      R__b.WriteVersion(AliITSv5asymm::IsA());
+      AliITS::Streamer(R__b);
+   } // end if R__b.IsReading()
+}
+*/
diff --git a/ITS/AliITSv5asymm.h b/ITS/AliITSv5asymm.h
new file mode 100644
index 00000000000..0b8d4b94c40
--- /dev/null
+++ b/ITS/AliITSv5asymm.h
@@ -0,0 +1,39 @@
+#ifndef ALIITSV5ASYMM_H
+#define ALIITSV5ASYMM_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/* $Id$ */
+
+////////////////////////////////////////////////////////////////////////
+//     Manager and hits classes for  ITS version 5 with symm services
+////////////////////////////////////////////////////////////////////////
+
+
+#include "AliITS.h"
+
+class TBRIK;
+class AliITSv5asymm : public AliITS {
+
+ public:
+    AliITSv5asymm();
+    AliITSv5asymm(const char *name, const char *title);
+    AliITSv5asymm(const AliITSv5asymm &source); // copy constructor
+    AliITSv5asymm& operator=(const AliITSv5asymm &source); // assignment operator	 
+    virtual       ~AliITSv5asymm() ;
+    virtual void   BuildGeometry();
+    virtual void  CreateGeometry();
+    virtual void  CreateMaterials();
+    virtual void  Init();   
+    virtual Int_t IsVersion() const {
+      // returns the ITS version number
+      return 5;
+    }
+    virtual void  StepManager();
+
+
+    ClassDef(AliITSv5asymm,1)//Hits manager for ITS version 5 Official detailed 
+	                         //geometry with asymmetric services
+};
+ 
+#endif
-- 
2.43.0