]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSv11.cxx
Reduced QA output (Yves)
[u/mrichter/AliRoot.git] / ITS / AliITSv11.cxx
index f1a264426545c8b3088a5a50e4bb71a783cf727c..cb4a1a6bf0bb7cd477f20f730f50f60195fe4f93 100644 (file)
@@ -1,5 +1,5 @@
 /**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * Copyright(c) 2007-2008, ALICE Experiment at CERN, All rights reserved. *
  *                                                                        *
  * Author: The ALICE Off-line Project.                                    *
  * Contributors are mentioned in the code where appropriate.              *
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
-/*
-$Log$
-Revision 1.6  2003/02/10 17:03:52  nilsen
-New version and structure of ITS V11 geometry. Work still in progress.
-
-Revision 1.5  2003/02/01 14:02:20  nilsen
-Work continues.
-
-Revision 1.4  2003/01/29 16:01:14  nilsen
-Update today's work.
-
-Revision 1.3  2003/01/28 17:59:54  nilsen
-Work continuing.
-
-Revision 1.2  2003/01/26 14:35:15  nilsen
-Some more geometry interface functions added and a start at the SSD support
-cone geometry. Committed to allow easy updates of partical work between 
-authors.
 
-Revision 1.1  2003/01/20 23:32:49  nilsen
-New ITS geometry. Only a Skeleton for now.
-
-$Id$
-*/
+//************************************************************************
+//
+//                 Inner Traking System geometry v11
+//
+//  Based on ROOT geometrical modeler
+//
+// B. Nilsen, L. Gaudichet
+//************************************************************************
 
-//////////////////////////////////////////////////////////////////////////////
-//                                                                          //
-//  Inner Traking System version 11                                         //
-//  This class contains the base procedures for the Inner Tracking System   //
-//                                                                          //
-// Authors: R. Barbera                                                      //
-// version 6.                                                               //
-// Created  2000.                                                           //
-//                                                                          //
-//  NOTE: THIS IS THE  SYMMETRIC PPR geometry of the ITS.                   //
-// THIS WILL NOT WORK                                                       //
-// with the geometry or module classes or any analysis classes. You are     //
-// strongly encouraged to uses AliITSv5.                                    //
-//                                                                          //
-//////////////////////////////////////////////////////////////////////////////
-// See AliITSv11::StepManager().
-#include <Riostream.h>
-#include <stdio.h>
-#include <stdlib.h>
-#include <TMath.h>
-#include <TFile.h>    // only required for Tracking function?
-#include <TObjArray.h>
-#include <TLorentzVector.h>
-#include <TObjString.h>
 #include <TClonesArray.h>
-#include <TBRIK.h>
-#include <TSystem.h>
-
+#include <TGeoManager.h>
+#include <TGeoPcon.h>
+#include <TGeoVolume.h>
+#include <TLorentzVector.h>
+#include <TVirtualMC.h>
 
-#include "AliRun.h"
-#include "AliMagF.h"
-#include "AliConst.h"
-#include "AliITSGeant3Geometry.h"
-#include "AliITShit.h"
-#include "AliITSv11.h"
+#include "AliITS.h"
+#include "AliITSCalibrationSDD.h"
+#include "AliITSDetTypeSim.h"
 #include "AliITSgeom.h"
-#include "AliITSgeomSPD.h"
-#include "AliITSgeomSDD.h"
-#include "AliITSgeomSSD.h"
-#include "AliITSDetType.h"
-#include "AliITSresponseSPD.h"
-#include "AliITSresponseSDD.h"
-#include "AliITSresponseSSD.h"
-#include "AliITSsegmentationSPD.h"
+#include "AliITShit.h"
 #include "AliITSsegmentationSDD.h"
+#include "AliITSsegmentationSPD.h"
 #include "AliITSsegmentationSSD.h"
-#include "AliITSsimulationSPD.h"
-#include "AliITSsimulationSDD.h"
-#include "AliITSsimulationSSD.h"
-#include "AliITSClusterFinderSPD.h"
-#include "AliITSClusterFinderSDD.h"
-#include "AliITSClusterFinderSSD.h"
-//
-#include "AliITSGeometryITSV.h"
-#include "AliITSGeometrySSDCone.h"
-#include "AliITSGeometrySDDCone.h"
+#include "AliITSv11.h"
+#include "AliITSv11GeometrySDD.h"
+#include "AliITSv11GeometrySPD.h"
+#include "AliITSv11GeometrySSD.h"
+#include "AliITSv11GeometrySupport.h"
+#include "AliMC.h"
+#include "AliMagF.h"
+#include "AliRun.h"
+#include "AliTrackReference.h"
 
 ClassImp(AliITSv11)
-
 //______________________________________________________________________
-AliITSv11::AliITSv11() : AliITS() {
-    // Standard default constructor for the ITS version 11.
-    // Inputs:
-    //   none.
-    // Outputs:
-    //   none.
-    // Return
-    //   A default constructed AliITSv11 class.
-
-    fITSV = 0;
-    fcS = 0;
-//   fcD = 0;
+AliITSv11::AliITSv11() : 
+fByThick(kTRUE),
+fMajorVersion(IsVersion()),
+fMinorVersion(0),
+fSPDgeom(),
+fSDDgeom(0),
+fSSDgeom(),
+fSupgeom(),
+fIgm(kv11)
+{
+  //    Standard default constructor for the ITS version 11.
+
+    fIdN          = 0;
+    fIdName       = 0;
+    fIdSens       = 0;
 }
-//______________________________________________________________________
-AliITSv11::AliITSv11(const char *title) : AliITS("ITS", title){
-    // Standard constructor for the ITS version 11.
-    // Inputs:
-    //   const char *title  The title of for this geometry.
-    // Outputs:
-    //   none.
-    // Return
-    //   A Standard constructed AliITSv11 class.
 
-    fITSV = 0;
-    fcS = 0;
-//    fcD = 0;
-}
-//______________________________________________________________________
-AliITSv11::~AliITSv11() {
-    // Standard destructor for the ITS version 11.
-    // Inputs:
-    //   none.
-    // Outputs:
-    //   none.
-    // Return
-    //   none.
 
-    if(fITSV!=0) delete fITSV;
-    if(fcS!=0) delete fcS;
-//    if(fcD!=0) delete fcD;
-}
 //______________________________________________________________________
-AliITSv11::AliITSv11(const AliITSv11 &source){
-    //     Copy Constructor for ITS version 11.
-    // Inputs:
-    //   AliITSv11 &source  class to be copied from.
-    // Outputs:
-    //   none.
-    // Return
-    //   none.
-
-    if(&source == this) return;
-    Error("Copy Constructor","Not allowed to copy AliITSv11");
-    return;
+AliITSv11::AliITSv11(const char *name, const char *title): 
+AliITS("ITS", title),
+fByThick(kTRUE),
+fMajorVersion(IsVersion()),
+fMinorVersion(0),
+fSPDgeom(),
+fSDDgeom(0),
+fSSDgeom(),
+fSupgeom(),
+fIgm(kv11)
+{
+  //    Standard constructor for the ITS version 11.
+
+  fSDDgeom = new AliITSv11GeometrySDD(0);
+
+  Int_t i;
+  fIdN = 6;
+  fIdName = new TString[fIdN];
+  fIdName[0] = name; // removes warning message
+  fIdName[0] = "ITS1";
+  fIdName[1] = "ITS2";
+  fIdName[2] = fSDDgeom->GetSenstiveVolumeName3();
+  fIdName[3] = fSDDgeom->GetSenstiveVolumeName4();
+  fIdName[4] = "ITS5";
+  fIdName[5] = "ITS6";
+  fIdSens    = new Int_t[fIdN];
+  for(i=0;i<fIdN;i++) fIdSens[i] = 0;
+  // not needed, fByThick set to kTRUE in in the member initialization lis
+  
 }
 //______________________________________________________________________
-AliITSv11& AliITSv11::operator=(const AliITSv11 &source){
-    //    Assignment operator for the ITS version 11.
-    // Inputs:
-    //   AliITSv11 &source  class to be copied from.
-    // Outputs:
-    //   none.
-    // Return
-    //   none.
-
-    if(&source == this) return *this;
-    Error("= operator","Not allowed to copy AliITSv11");
-    return *this;
+AliITSv11::AliITSv11(Int_t debugITS,Int_t debugSPD,Int_t debugSDD,
+                  Int_t debugSSD,Int_t debugSUP) :
+AliITS("ITS","ITS geometry v11"),
+fByThick(kTRUE),
+fMajorVersion(IsVersion()),
+fMinorVersion(0),
+fSPDgeom(),
+fSDDgeom(0),
+fSSDgeom(),
+fSupgeom(),
+fIgm(kv11)
+{
+  // Standard default constructor for the ITS version 11.
+
+
+  //   fSPDgeom = new AliITSv11GeometrySPD(debugSPD);
+  fSDDgeom = new AliITSv11GeometrySDD(debugSDD);
+  fSDDgeom->SetDebug(debugSDD);
+  //   fSupgeom = new AliITSv11GeometrySupport(debugSUP);
+
+  Int_t i;
+  fIdN = 6;
+  fIdName = new TString[fIdN];
+  fIdName[0] = fSPDgeom->GetSenstiveVolumeName1();
+  fIdName[1] = fSPDgeom->GetSenstiveVolumeName2();
+  fIdName[2] = fSDDgeom->GetSenstiveVolumeName3();
+  fIdName[3] = fSDDgeom->GetSenstiveVolumeName4();
+  fIdName[4] = fSSDgeom->GetSenstiveVolumeName5();
+  fIdName[5] = fSSDgeom->GetSenstiveVolumeName6();
+  fIdSens    = new Int_t[fIdN];
+  for(i=0;i<fIdN;i++) fIdSens[i] = 0;
+  debugITS = (debugSPD && debugSSD && debugSUP && debugSDD); //remove temp. warnings
 }
 //______________________________________________________________________
-void AliITSv11::BuildGeometry(){
-    // This routine defines and Creates the geometry for version 11 of 
-    // the ITS for use in the simulation display routines. This is a 
-    // very simplified geometry for speed of viewing.
-    // Inputs:
-    //   none.
-    // Outputs:
-    //   none.
-    // Return
-    //   none.
-    TVector3 t(0.0,0.0,0.0);
-
-    if(fITSV==0) fITSV = new AliITSGeometryITSV(this,"ALIC");
-    if(fcS==0) fcS = new AliITSGeometrySSDCone(this,t,"TSV",1);
-
-    fcS->BuildDisplayGeometry();
+AliITSv11::~AliITSv11() {
+  delete fSDDgeom;
 }
 //______________________________________________________________________
 void AliITSv11::CreateGeometry(){
-    // This routine defines and Creates the geometry for version 11 of 
-    // the ITS. The geometry is used by the particle trasport routines,
-    // and therefore, is very detailed.
-    // Inputs:
-    //   none.
-    // Outputs:
-    //   none.
-    // Return
-    //   none.
-    TVector3 t(0.0,0.0,0.0);
-
-    if(fITSV==0) fITSV = new AliITSGeometryITSV(this,"ALIC");
-    if(fcS==0) fcS = new AliITSGeometrySSDCone(this,t,"TSV",1);
     //
-    fITSV->CreateG3Geometry();
-    fcS->CreateG3Geometry("TSV",t);
+    // Create ROOT geometry
     //
-    fITSV->PositionGeometry("ALIC",1,t,0);
-    fcS->PositionG3Geometry(fITSV->GetParams(),1,t,0);
+    // These constant character strings are set by cvs during commit
+    // do not change them unless you know what you are doing!
+    const Char_t *cvsDate="$Date$";
+    const Char_t *cvsRevision="$Revision$";
+
+    TGeoManager *geoManager = gGeoManager;
+    TGeoVolume *vALIC = geoManager->GetTopVolume();
+
+    TGeoPcon *sITS = new TGeoPcon("ITS Top Volume",0.0,360.0,2);
+
+    // DefineSection(section number, Z, Rmin, Rmax).
+    const Double_t kcm = 1.0;
+    sITS->DefineSection(0,-300.0*kcm,0.01*kcm,50.0*kcm);
+    sITS->DefineSection(1,+300.0*kcm,0.01*kcm,50.0*kcm);
+
+    TGeoMedium *air = gGeoManager->GetMedium("ITS_AIR$");
+    TGeoVolume *vITS = new TGeoVolume("ITSV",sITS,air);
+    vITS->SetVisibility(kFALSE);
+    const Int_t length=100;
+    Char_t vstrng[length];
+    if(fIgm.WriteVersionString(vstrng,length,(AliITSVersion_t)IsVersion(),
+                               fMinorVersion,cvsDate,cvsRevision))
+        vITS->SetTitle(vstrng);
+    //printf("Title set to %s\n",vstrng);
+    vALIC->AddNode(vITS,1,0);
+
+//   fSPDgeom->CenteralSPD(vITS);
+
+  fSDDgeom->Layer3(vITS);
+  fSDDgeom->Layer4(vITS);
+
+//     fSupgeom->SPDCone(vITS);
+//     fSupgeom->SPDThermalSheald(vITS);
+//     fSupgeom->SDDCone(vITS);
+//     fSupgeom->SSDCone(vITS);
+//     fSupgeom->ServicesCableSupport(vITS);
+
 }
 //______________________________________________________________________
 void AliITSv11::CreateMaterials(){
-    // Create ITS materials
-    //     This function defines the default materials used in the Geant
-    // Monte Carlo simulations for the geometries AliITSv11.
-    // In general it is automatically replaced by
+    // Create Standard ITS Materials
     // Inputs:
-    //   none.
+    //  none.
     // Outputs:
-    //   none.
-    // Return
-    //   none.
-/*
-    TVector3 t(0.0,0.0,0.0);
-
-    if(fITSV==0) fITSV = new AliITSGeometryITSV(this,"ALIC");
-    if(fcS==0) fcS = new AliITSGeometrySSDCone(this,t,"TSV",1);
+    //  none.
+    // Return:
+    // none.
 
-    fITSV->CreateG3Materials();
-    fcS->CreateG3Materials();
+    
+    //
+    fSPDgeom->AliITSv11Geometry::CreateDefaultMaterials();
+    // Detector specific material definistions
+    fSPDgeom->CreateMaterials();
+    fSDDgeom->CreateMaterials();
+    fSSDgeom->CreateMaterials();
+    fSupgeom->CreateMaterials();
+}
+/*
+//______________________________________________________________________
+void AliITSv11::InitAliITSgeom(){
+  //
+  // Fill fITSgeom with the 3 sub-detector geometries
+  //
+
+  if (gGeoManager) gGeoManager->Export("geometry.root");
+
+    const Int_t knlayers = 6;
+    const Int_t kndeep = 3;
+    const AliITSDetector kidet[knlayers]={kSPD,kSPD,kSDD,kSDD,kSSD,kSSD};
+    const TString knames[knlayers] = {
+      "AliITSv11:spd missing",  // lay=1
+      "AliITSv11:spd missing",  // lay=2
+      "/ALIC_1/ITSV_1/ITSsddLayer3_1/ITSsddLadd_%d/ITSsddSensor_%d/ITSsddWafer_%d", // lay=3
+      "/ALIC_1/ITSV_1/ITSsddLayer4_1/ITSsddLadd_%d/ITSsddSensor_%d/ITSsddWafer_%d", // lay=4
+      "AliITSv11:ssd missing",  // lay=5
+      "AliITSv11:ssd missing"   // lay=6
+    };
+
+    const Int_t kitsGeomTreeCopys[knlayers][kndeep]= {{10, 2, 4},// lay=1
+                                                     {10, 4, 4}, // lay=2
+                                                     {14, 6, 1}, // lay=3
+                                                     {22, 8, 1}, // lay=4
+                                                     {34,22, 1}, // lay=5
+                                                     {38,25, 1}};// lay=6
+    Int_t       nlad[knlayers],ndet[knlayers];
+    Int_t       mod,lay,lad=0,det=0,i,j,k,cp0,cp1,cp2;
+    TString path,shapeName;
+    TGeoHMatrix materix;
+    Double_t trans[3]={3*0.0},rot[10]={9*0.0,1.0};
+    TArrayD shapePar;
+    TArrayF shapeParF;
+    Bool_t shapeDefined[3]={kFALSE,kFALSE,kFALSE};
+
+    AliDebug(1,"Reading Geometry transformation directly from Modler.");
+    mod = 0;
+    for(i=0;i<knlayers;i++){
+        k = 1;
+        for(j=0;j<kndeep;j++) if(kitsGeomTreeCopys[i][j]!=0)
+            k *= TMath::Abs(kitsGeomTreeCopys[i][j]);
+        mod += k;
+    } // end for i
+
+    SetITSgeom(0);
+    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;
+    AliITSgeom* geom = new AliITSgeom(0,6,nlad,ndet,mod);
+    SetITSgeom(geom);
+    mod = 0;
+    for(lay=1;lay<=knlayers;lay++){
+
+        for(cp0=0; cp0<kitsGeomTreeCopys[lay-1][0]; cp0++){
+            for(cp1=0; cp1<kitsGeomTreeCopys[lay-1][1]; cp1++){
+                for(cp2=1; cp2<=kitsGeomTreeCopys[lay-1][2]; cp2++){
+
+                    path.Form(knames[lay-1].Data(),
+                              cp0,cp1,cp2);
+                    switch (lay){
+                    case 1:{
+                        det = cp2;
+                        lad = cp1+2*(cp0-1);
+                    }break;
+                    case 2:{
+                        det = cp2;
+                        lad = cp1+4*(cp0-1);
+                    } break;
+                    case 3: case 4: case 5: case 6:{
+                        det = cp1;
+                        lad = cp0;
+                    } break;
+                    } // end switch
+                         //AliInfo(Form("path=%s lay=%d lad=%d det=%d",
+                         //             path.Data(),lay,lad,det));
+                    gMC->GetTransformation(path.Data(),materix);
+                    gMC->GetShape(path.Data(),shapeName,shapePar);
+                    shapeParF.Set(shapePar.GetSize());
+                    for(i=0;i<shapePar.GetSize();i++) shapeParF[i]=shapePar[i];
+                    geom->CreateMatrix(mod,lay,lad,det,kidet[lay-1],trans,rot);
+                    geom->SetTrans(mod,materix.GetTranslation());
+                    geom->SetRotMatrix(mod,materix.GetRotationMatrix());
+                   geom->GetGeomMatrix(mod)->SetPath(path.Data());
+                    switch (lay){
+                    case 1: case 2:
+                       if(!shapeDefined[kSPD]){
+                        geom->ReSetShape(kSPD,new AliITSgeomSPD425Short(
+                                shapeParF.GetSize(),shapeParF.GetArray()));
+                       shapeDefined[kSPD] = kTRUE;
+                    }break;
+                    case 3: case 4:
+                       if(!shapeDefined[kSDD]){
+                        geom->ReSetShape(kSDD,new AliITSgeomSDD256(
+                                shapeParF.GetSize(),shapeParF.GetArray()));
+                       shapeDefined[kSDD] = kTRUE;
+                    }break;
+                    case 5: case 6:
+                       if(!shapeDefined[kSSD]){
+                        geom->ReSetShape(kSSD,new AliITSgeomSSD75and275(
+                                shapeParF.GetSize(),shapeParF.GetArray()));
+                       shapeDefined[kSSD] = kTRUE;
+                    }break;
+                    default:{
+                    }break;
+                    } // end switch
+                    mod++;
+                } /// end for cp2
+            } // end for cp1
+        } // end for cp0
+    } // end for lay
+
+//   fSDDgeom->ExportSensorGeometry(GetITSgeom(), +3, 0);  //SDD
+}
 */
+//______________________________________________________________________
+void AliITSv11::Init(){
+  //
+  //     Initialise the ITS after it has been created.
+  //
 
+  //AliInfo(Form("Minor version %d",fMinorVersion));
+    //
+    UpdateInternalGeometry();
+    AliITS::Init();
 
-  Int_t   ifield = gAlice->Field()->Integ();
-  Float_t fieldm = gAlice->Field()->Max();
-
-  Float_t tmaxfd = 0.1; // 1.0; // Degree
-  Float_t stemax = 1.0; // cm
-  Float_t deemax = 0.1; // 30.0; // Fraction of particle's energy 0<deemax<=1
-  Float_t epsil  = 1.0E-4; // 1.0; // cm
-  Float_t stmin  = 0.0; // cm "Default value used"
-
-  Float_t tmaxfdSi = 0.1; // .10000E+01; // Degree
-  Float_t stemaxSi = 0.0075; //  .10000E+01; // cm
-  Float_t deemaxSi = 0.1; // 0.30000E-02; // Fraction of particle's energy 0<deemax<=1
-  Float_t epsilSi  = 1.0E-4;// .10000E+01;
-  Float_t stminSi  = 0.0; // cm "Default value used"
-
-  Float_t tmaxfdAir = 0.1; // .10000E+01; // Degree
-  Float_t stemaxAir = .10000E+01; // cm
-  Float_t deemaxAir = 0.1; // 0.30000E-02; // Fraction of particle's energy 0<deemax<=1
-  Float_t epsilAir  = 1.0E-4;// .10000E+01;
-  Float_t stminAir  = 0.0; // cm "Default value used"
-
-  Float_t tmaxfdServ = 1.0; // 10.0; // Degree
-  Float_t stemaxServ = 1.0; // 0.01; // cm
-  Float_t deemaxServ = 0.5; // 0.1; // Fraction of particle's energy 0<deemax<=1
-  Float_t epsilServ  = 1.0E-3; // 0.003; // cm
-  Float_t stminServ  = 0.0; //0.003; // cm "Default value used"
-
-  // Freon
-  Float_t afre[2]  = { 12.011,18.9984032 };
-  Float_t zfre[2]  = { 6., 9. };
-  Float_t wfre[2]  = { 5.,12. };
-  Float_t densfre  = 1.5;
-
-  // --- Define the various materials and media for GEANT --- 
-  // AliMaterial(Int_t imat, const char* name, Float_t a, Float_t z,
-  //              Float_t dens, Float_t radl, Float_t absl,
-  //              Float_t *buf=0, Int_t nwbuf=0)
-  //AliMedium(Int_t numed, const char *name, Int_t nmat,
-  //          Int_t isvol, Int_t ifield, Float_t fieldm,
-  //          Float_t tmaxfd, Float_t stemax, Float_t deemax,
-  //          Float_t epsil, Float_t stmin, Float_t *ubuf=0, Int_t nbuf=0)
-  AliMaterial(1,"SI$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
-  AliMedium(1,"SI$",1,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
-
-  AliMaterial(2,"SPD SI CHIP$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
-  AliMedium(2,"SPD SI CHIP$",2,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
-
-  AliMaterial(3,"SPD SI BUS$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
-  AliMedium(3,"SPD SI BUS$",3,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
-
-  AliMaterial(4,"C (M55J)$",0.12011E+02,0.60000E+01,0.1930E+01,0.22100E+02,0.99900E+03);
-  AliMedium(4,"C (M55J)$",4,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(5,"AIR$",0.14610E+02,0.73000E+01,0.12050E-02,0.30423E+05,0.99900E+03);
-  AliMedium(5,"AIR$",5,0,ifield,fieldm,tmaxfdAir,stemaxAir,deemaxAir,epsilAir,stminAir);
-
-  AliMaterial(6,"GEN AIR$",0.14610E+02,0.73000E+01,0.12050E-02,0.30423E+05,0.99900E+03);
-  AliMedium(6,"GEN AIR$",6,0,ifield,fieldm,tmaxfdAir,stemaxAir,deemaxAir,epsilAir,stminAir);
-
-  AliMaterial(7,"SDD SI CHIP$",0.374952E+02,0.178184E+02,0.24485E+01,0.76931E+01,0.99900E+03);
-  AliMedium(7,"SDD SI CHIP$",7,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
-
-  AliMaterial(9,"SDD C (M55J)$",0.123565E+02,0.64561E+01,0.18097E+01,0.229570E+02,0.99900E+03);
-  AliMedium(9,"SDD C (M55J)$",9,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(10,"SDD AIR$",0.14610E+02,0.73000E+01,0.12050E-02,0.30423E+05,0.99900E+03);
-  AliMedium(10,"SDD AIR$",10,0,ifield,fieldm,tmaxfdAir,stemaxAir,deemaxAir,epsilAir,stminAir);
-
-  AliMaterial(11,"AL$",0.26982E+02,0.13000E+02,0.26989E+01,0.89000E+01,0.99900E+03);
-  AliMedium(11,"AL$",11,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(12,"WATER$",0.14322E+02,0.72167E+01,0.10000E+01,0.35759E+02,0.94951E+02);
-  AliMedium(12,"WATER$",12,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMixture(13,"Freon$",afre,zfre,densfre,-2,wfre);
-  AliMedium(13,"Freon$",13,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(14,"COPPER$",0.63546E+02,0.29000E+02,0.89600E+01,0.14300E+01,0.99900E+03);
-  AliMedium(14,"COPPER$",14,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(15,"CERAMICS$",0.22314E+02,0.10856E+02,0.36000E+01,0.76200E+01,0.31901E+02);
-  AliMedium(15,"CERAMICS$",15,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(20,"SSD C (M55J)$",0.12011E+02,0.60000E+01,0.1930E+01,0.22100E+02,0.99900E+03);
-  AliMedium(20,"SSD C (M55J)$",20,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(21,"SSD AIR$",0.14610E+02,0.73000E+01,0.12050E-02,0.30423E+05,0.99900E+03);
-  AliMedium(21,"SSD AIR$",21,0,ifield,fieldm,tmaxfdAir,stemaxAir,deemaxAir,epsilAir,stminAir);
-
-  AliMaterial(25,"G10FR4$",0.17749E+02,0.88750E+01,0.18000E+01,0.21822E+02,0.99900E+03);
-  AliMedium(25,"G10FR4$",25,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(26,"GEN C (M55J)$",0.12011E+02,0.60000E+01,0.1930E+01,0.22100E+02,0.99900E+03);
-  AliMedium(26,"GEN C (M55J)$",26,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(27,"GEN Air$",0.14610E+02,0.73000E+01,0.12050E-02,0.30423E+05,0.99900E+03);
-  AliMedium(27,"GEN Air$",27,0,ifield,fieldm,tmaxfdAir,stemaxAir,deemaxAir,epsilAir,stminAir);
-
-  AliMaterial(51,"SPD SI$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
-  AliMedium(51,"SPD SI$",51,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
-
-  AliMaterial(52,"SPD SI CHIP$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
-  AliMedium(52,"SPD SI CHIP$",52,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
-
-  AliMaterial(53,"SPD SI BUS$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
-  AliMedium(53,"SPD SI BUS$",53,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
-
-  AliMaterial(54,"SPD C (M55J)$",0.12011E+02,0.60000E+01,0.1930E+01,0.22100E+02,0.99900E+03);
-  AliMedium(54,"SPD C (M55J)$",54,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(55,"SPD AIR$",0.14610E+02,0.73000E+01,0.12050E-02,0.30423E+05,0.99900E+03);
-  AliMedium(55,"SPD AIR$",55,0,ifield,fieldm,tmaxfdAir,stemaxAir,deemaxAir,epsilAir,stminAir);
-
-  AliMaterial(56,"SPD KAPTON(POLYCH2)$",0.14000E+02,0.71770E+01,0.13000E+01,0.31270E+02,0.99900E+03);
-  AliMedium(56,"SPD KAPTON(POLYCH2)$",56,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(61,"EPOXY$",0.17749E+02,0.88750E+01,0.18000E+01,0.21822E+02,0.99900E+03);
-  AliMedium(61,"EPOXY$",61,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(62,"SILICON$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
-  AliMedium(62,"SILICON$",62,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
-
-  AliMaterial(63,"KAPTONH(POLYCH2)$",0.14000E+02,0.71770E+01,0.13000E+01,0.31270E+02,0.99900E+03);
-  AliMedium(63,"KAPTONH(POLYCH2)$",63,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(64,"ALUMINUM$",0.26982E+02,0.13000E+02,0.26989E+01,0.89000E+01,0.99900E+03);
-  AliMedium(64,"ALUMINUM$",64,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(65,"INOX$",0.55098E+02,0.2572E+02,0.7900E+01,0.17800E+01,0.99900E+03);
-  AliMedium(65,"INOX$",65,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(68,"ROHACELL$",0.123974E+02,0.62363E+01,0.500E-01,0.80986E+03,0.99900E+03);
-  AliMedium(68,"ROHACELL$",68,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(69,"SDD C AL (M55J)$",0.138802E+02,0.71315E+01,0.19837E+01,0.176542E+02,0.99900E+03);
-  AliMedium(69,"SDD C AL (M55J)$",69,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(70,"SDDKAPTON (POLYCH2)$",0.14000E+02,0.71770E+01,0.13000E+01,0.31270E+02,0.99900E+03);
-  AliMedium(70,"SDDKAPTON (POLYCH2)$",70,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(71,"ITS SANDW A$",0.12011E+02,0.60000E+01,0.2115E+00,0.17479E+03,0.99900E+03);
-  AliMedium(71,"ITS SANDW A$",71,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(72,"ITS SANDW B$",0.12011E+02,0.60000E+01,0.27000E+00,0.18956E+03,0.99900E+03);
-  AliMedium(72,"ITS SANDW B$",72,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(73,"ITS SANDW C$",0.12011E+02,0.60000E+01,0.41000E+00,0.90868E+02,0.99900E+03);
-  AliMedium(73,"ITS SANDW C$",73,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-
-  AliMaterial(74,"HEAT COND GLUE$",0.12011E+02,0.60000E+01,0.1930E+01,0.22100E+02,0.99900E+03);
-  AliMedium(74,"HEAT COND GLUE$",74,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
+    //
+/*
+    if(GetITSgeom()!=0) SetITSgeom(0x0);
+    AliITSgeom* geom = new AliITSgeom();
+    SetITSgeom(geom);
+    AliITS::Init();
+*/    //
+}
 
-  AliMaterial(75,"ELASTO SIL$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
-  AliMedium(75,"ELASTO SIL$",75,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
 
-  AliMaterial(76,"SPDBUS(AL+KPT+EPOX)$",0.19509E+02,0.96502E+01,0.19060E+01,0.15413E+02,0.99900E+03);
-  AliMedium(76,"SPDBUS(AL+KPT+EPOX)$",76,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
-               
-  AliMaterial(77,"SDD X7R capacitors$",0.1157516E+03,0.477056E+02,0.67200E+01,0.14236E+01,0.99900E+03);
-  AliMedium(77,"SDD X7R capacitors$",77,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
+//______________________________________________________________________
+void AliITSv11::SetDefaults(){
+  //
+  // Set response and segmentation models for SPD, SDD and SSD
+  //
+
+//     if(!fDetTypeSim) fDetTypeSim = new AliITSDetTypeSim();
+//     fDetTypeSim->SetITSgeom(GetITSgeom());
+    if(!fDetTypeSim) {
+      Warning("SetDefaults","Error fDetTypeSim not defined");
+      return;
+    }
+  
+    fDetTypeSim->ResetCalibrationArray();
+    fDetTypeSim->ResetSegmentation();
+    fDetTypeSim->SetDefaults();
+    
+    if(fgkNTYPES>3){
+       Warning("SetDefaults",
+               "Only the four basic detector types are initialised!");
+    }// end if
+
+    
+    return;
+}
 
-  AliMaterial(78,"SDD ruby sph. Al2O3$",0.218101E+02,0.106467E+02,0.39700E+01,0.48539E+01,0.99900E+03);
-  AliMedium(78,"SDD ruby sph. Al2O3$",78,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
 
-  AliMaterial(79,"SDD SI insensitive$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
-  AliMedium(79,"SDD SI insensitive$",79,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
 
-  AliMaterial(80,"SDD HV microcable$",0.159379E+02,0.78598E+01,0.16087E+01,0.217906E+02,0.99900E+03);
-  AliMedium(80,"SDD HV microcable$",80,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
 
-  AliMaterial(81,"SDD LV+signal cable$",0.223689E+02,0.108531+02,0.21035E+01,0.13440E+02,0.99900E+03);
-  AliMedium(81,"SDD LV+signal cable$",81,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
 
-  AliMaterial(82,"SDD hybrid microcab$",0.218254E+02,0.106001E+02,0.20502E+01,0.137308E+02,0.99900E+03);
-  AliMedium(82,"SDD hybrid microcab$",82,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
+//______________________________________________________________________
+void AliITSv11::DrawModule() const{
 
-  AliMaterial(83,"SDD anode microcab$",0.186438E+02,0.91193E+01,0.17854E+01,0.176451E+02,0.99900E+03);
-  AliMedium(83,"SDD anode microcab$",83,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
+}
 
-  AliMaterial(84,"SDD/SSD rings$",0.123565E+02,0.64561E+01,0.18097E+01,0.229570E+02,0.99900E+03);
-  AliMedium(84,"SDD/SSD rings$",84,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
+// //______________________________________________________________________
+// void AliITSv11::StepManager(){
+//   //
+//   //    Called for every step in the ITS, then calles the AliITShit class
+//   // creator with the information to be recoreded about that hit.
+//   //
+//     Int_t         copy, id;
+//     TLorentzVector position, momentum;
+//     static TLorentzVector position0;
+//     static Int_t stat0=0;
+
+//     if(!(this->IsActive())){
+//     return;
+//     } // end if !Active volume.
+
+//     if(!(gMC->TrackCharge())) return;
+
+//     id=gMC->CurrentVolID(copy);
+
+//     Bool_t sensvol = kFALSE;
+//     for(Int_t kk=0;kk<6;kk++)if(id == fIdSens[kk])sensvol=kTRUE;
+//     if(sensvol && (gMC->IsTrackExiting())){
+//     copy = fTrackReferences->GetEntriesFast();
+//     TClonesArray &lTR = *fTrackReferences;
+//     // Fill TrackReference structure with this new TrackReference.
+//     new(lTR[copy]) AliTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber());
+//     } // if Outer ITS mother Volume
+
+
+//     Int_t   copy1,copy2;  
+//     Int_t   vol[5];
+//     TClonesArray &lhits = *fHits;
+//     //
+//     // Track status
+//     vol[3] = 0;
+//     vol[4] = 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(2,copy);
+//     //detector copy in the ladder = 1<->4  (ITS1 < I101 < I103 < I10A)
+//     vol[1] = copy;
+//     gMC->CurrentVolOffID(3,copy1);
+//     //ladder copy in the module   = 1<->2  (I10A < I12A)
+//     gMC->CurrentVolOffID(4,copy2);
+//     //module copy in the layer    = 1<->10 (I12A < IT12)
+//     vol[2] = copy1+(copy2-1)*2;//# of ladders in one module  = 2
+//     } else if(id == fIdSens[1]){
+//     vol[0] = 2;
+//     id = gMC->CurrentVolOffID(2,copy);
+//     //detector copy in the ladder = 1<->4  (ITS2 < I1D1 < I1D3 < I20A)
+//     vol[1] = copy;
+//     gMC->CurrentVolOffID(3,copy1);
+//     //ladder copy in the module   = 1<->4  (I20A < I12A)
+//     gMC->CurrentVolOffID(4,copy2);
+//     //module copy in the layer    = 1<->10 (I12A < IT12)
+//     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<->6  (ITS3 < I302 < I004)
+//     vol[1] = copy;
+//     id = gMC->CurrentVolOffID(2,copy);
+//     //ladder copy in the layer    = 1<->14 (I004 < IT34)
+//     vol[2] = copy;
+//     } else if(id == fIdSens[3]){
+//     vol[0] = 4;
+//     id = gMC->CurrentVolOffID(1,copy);
+//     //detector copy in the ladder = 1<->8  (ITS4 < I402 < I005)
+//     vol[1] = copy;
+//     id = gMC->CurrentVolOffID(2,copy);
+//     //ladder copy in the layer    = 1<->22 (I005 < IT34))
+//     vol[2] = copy;
+//     }else if(id == fIdSens[4]){
+//     vol[0] = 5;
+//     id = gMC->CurrentVolOffID(1,copy);
+//     //detector copy in the ladder = 1<->22  (ITS5 < I562 < I565)
+//     vol[1] = copy;
+//     id = gMC->CurrentVolOffID(2,copy);
+//     //ladder copy in the layer    = 1<->34 (I565 < IT56)
+//     vol[2] = copy;
+//     }else if(id == fIdSens[5]){
+//     vol[0] = 6;
+//     id = gMC->CurrentVolOffID(1,copy);
+//     //detector copy in the ladder = 1<->25  (ITS6 < I566 < I569)
+//     vol[1] = copy;
+//     id = gMC->CurrentVolOffID(2,copy);
+//     //ladder copy in the layer = 1<->38 (I569 < IT56)
+//     vol[2] = copy;
+//     } else {
+//     return; // not an ITS volume?
+//     } // end if/else if (gMC->CurentVolID(copy) == fIdSens[i])
+//     //
+//     gMC->TrackPosition(position);
+//     gMC->TrackMomentum(momentum);
+//     vol[4] = stat0;
+//     if(gMC->IsTrackEntering()){
+//     position0 = position;
+//     stat0 = vol[3];
+//     return;
+//     } // end if IsEntering
+//     // Fill hit structure with this new hit.
+    
+//     new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->GetMCApp()->GetCurrentTrackNumber(),vol,
+//                                gMC->Edep(),gMC->TrackTime(),position,
+//                                position0,momentum);
+
+//     position0 = position;
+//     stat0 = vol[3];
+
+//     return;
+// }
 
-  AliMaterial(85,"inox/alum$",0.321502E+02,0.153383E+02,0.30705E+01,0.69197E+01,0.99900E+03);
-  AliMedium(85,"inox/alum$",85,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);
 
+//______________________________________________________________________
+void AliITSv11::StepManager(){
+  //
+  //    Called for every step in the ITS, then calles the AliITShit class
+  // creator with the information to be recoreded about that hit.
+  //
+    Int_t         copy, id;
+    TLorentzVector position, momentum;
+    static TLorentzVector position0;
+    static Int_t stat0=0;
 
-  // Special media
+    if(!(this->IsActive())){
+       return;
+    } // end if !Active volume.
 
-  AliMaterial(90,"SPD shield$", 12.011, 6., 1.93/10. , 22.1*10., 999);
-  AliMedium(90,"SPD shield$",90,0,ifield,fieldm,tmaxfdServ,stemaxServ,deemaxServ,epsilServ,stminServ);
+    if(!(gMC->TrackCharge())) return;
 
-  AliMaterial(91, "SPD End ladder$", 47.0447, 21.7963, 3.6374, 4.4711, 999); 
-  AliMedium(91,"SPD End ladder$",91,0,ifield,fieldm,tmaxfdServ,stemaxServ,deemaxServ,epsilServ,stminServ);
+    id=gMC->CurrentVolID(copy);
 
-  AliMaterial(92, "SPD cone$",28.0855, 14., 2.33, 9.36, 999);    
-  AliMedium(92,"SPD cone$",92,0,ifield,fieldm,tmaxfdServ,stemaxServ,deemaxServ,epsilServ,stminServ);
+    Bool_t sensvol = kFALSE;
+    for(Int_t kk=0;kk<6;kk++)if(id == fIdSens[kk])sensvol=kTRUE;
+    if(sensvol && (gMC->IsTrackExiting())){
+       AddTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber(), AliTrackReference::kITS);
+    } // if Outer ITS mother Volume
 
-  AliMaterial(93, "SDD End ladder$", 69.9298, 29.8246, 0.3824, 36.5103, 999); 
-  AliMedium(93,"SDD End ladder$",93,0,ifield,fieldm,tmaxfdServ,stemaxServ,deemaxServ,epsilServ,stminServ);
 
-  AliMaterial(94, "SDD cone$",63.546, 29., 1.15, 1.265, 999);
-  AliMedium(94,"SDD cone$",94,0,ifield,fieldm,tmaxfdServ,stemaxServ,deemaxServ,epsilServ,stminServ);
+    Int_t   copy1,copy2;  
+    Int_t   vol[5];
+    TClonesArray &lhits = *fHits;
+    //
+    // Track status
+    vol[3] = 0;
+    vol[4] = 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(2,copy);
+       //detector copy in the ladder = 1<->4  (ITS1 < I101 < I103 < I10A)
+       vol[1] = copy;
+       gMC->CurrentVolOffID(3,copy1);
+       //ladder copy in the module   = 1<->2  (I10A < I12A)
+       gMC->CurrentVolOffID(4,copy2);
+       //module copy in the layer    = 1<->10 (I12A < IT12)
+       vol[2] = copy1+(copy2-1)*2;//# of ladders in one module  = 2
+
+    } else if(id == fIdSens[1]){
+       vol[0] = 2;
+       id = gMC->CurrentVolOffID(2,copy);
+       //detector copy in the ladder = 1<->4  (ITS2 < I1D1 < I1D3 < I20A)
+       vol[1] = copy;
+       gMC->CurrentVolOffID(3,copy1);
+       //ladder copy in the module   = 1<->4  (I20A < I12A)
+       gMC->CurrentVolOffID(4,copy2);
+       //module copy in the layer    = 1<->10 (I12A < IT12)
+       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<->6  (ITS3 < I302 < I004)
+       vol[1] = copy;
+       id = gMC->CurrentVolOffID(2,copy);
+       //ladder copy in the layer    = 1<->14 (I004 < IT34)
+       vol[2] = copy;
+
+    } else if(id == fIdSens[3]){
+       vol[0] = 4;
+       id = gMC->CurrentVolOffID(1,copy);
+       //detector copy in the ladder = 1<->8  (ITS4 < I402 < I005)
+       vol[1] = copy;
+       id = gMC->CurrentVolOffID(2,copy);
+       //ladder copy in the layer    = 1<->22 (I005 < IT34))
+       vol[2] = copy;
+
+    }else if(id == fIdSens[4]){
+       vol[0] = 5;
+       id = gMC->CurrentVolOffID(1,copy);
+       //detector copy in the ladder = 1<->22  (ITS5 < I562 < I565)
+       vol[1] = copy;
+       id = gMC->CurrentVolOffID(2,copy);
+       //ladder copy in the layer    = 1<->34 (I565 < IT56)
+       vol[2] = copy;
+
+    }else if(id == fIdSens[5]){
+       vol[0] = 6;
+       id = gMC->CurrentVolOffID(1,copy);
+       //detector copy in the ladder = 1<->25  (ITS6 < I566 < I569)
+       vol[1] = copy;
+       id = gMC->CurrentVolOffID(2,copy);
+       //ladder copy in the layer = 1<->38 (I569 < IT56)
+       vol[2] = copy;
+    } else {
+       return; // not an ITS volume?
+    } // end if/else if (gMC->CurentVolID(copy) == fIdSens[i])
+    //
+    gMC->TrackPosition(position);
+    gMC->TrackMomentum(momentum);
+    vol[4] = stat0;
+    if(gMC->IsTrackEntering()){
+       position0 = position;
+       stat0 = vol[3];
+       return;
+    } // end if IsEntering
+    // Fill hit structure with this new hit.
+    
+    new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->GetMCApp()->GetCurrentTrackNumber(),vol,
+                                  gMC->Edep(),gMC->TrackTime(),position,
+                                  position0,momentum);
+
+    position0 = position;
+    stat0 = vol[3];
 
-  AliMaterial(95, "SSD End ladder$", 32.0988, 15.4021, 0.68, 35.3238, 999); 
-  AliMedium(95,"SSD End ladder$",95,0,ifield,fieldm,tmaxfdServ,stemaxServ,deemaxServ,epsilServ,stminServ);
-  
-  AliMaterial(96, "SSD cone$",63.546, 29., 1.15, 1.265, 999);
-  AliMedium(96,"SSD cone$",96,0,ifield,fieldm,tmaxfdServ,stemaxServ,deemaxServ,epsilServ,stminServ);
-}
-//______________________________________________________________________
-void AliITSv11::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.
-    // Inputs:
-    //   none.
-    // Outputs:
-    //   none.
-    // Return
-    //   none.
-}
-//______________________________________________________________________
-void AliITSv11::Init(){
-    // Initialise the ITS after it has been created.
-    // Inputs:
-    //   none.
-    // Outputs:
-    //   none.
-    // Return
-    //   none.
-}
-//______________________________________________________________________
-void AliITSv11::SetDefaults(){
-    // Sets the default segmentation, response, digit and raw cluster 
-    // classes to be used. These defaults can be overwritten in the
-    // macros that do these later steps. Defaults are give hear for the
-    // general user.
-    // Inputs:
-    //   none.
-    // Outputs:
-    //   none.
-    // Return
-    //   none.
-}
-//______________________________________________________________________
-void AliITSv11::DrawModule(){
-    // Draw a standard set of shaded view of the ITS version 11.
-    // Inputs:
-    //   none.
-    // Outputs:
-    //   none.
-    // Return
-    //   none.
-}
-//______________________________________________________________________
-void AliITSv11::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
-    // is the extra variables and the like used in the printing allocated.
+    return;
 }
+