ROOT logo
/**************************************************************************
 * Copyright(c) 2007-2009, 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.                  *
 **************************************************************************/


/* $Id: AliITSUv0.cxx */


//========================================================================
//
//        Geometry for the Upgrade of the Inner Tracking System
//
// Mario Sitta (sitta@to.infn.it)
// Chinorat Kobdaj (kobdaj@g.sut.ac.th)
//
//========================================================================



// $Log: AliITSUv0.cxx,v $

#include <TClonesArray.h>
#include <TGeoGlobalMagField.h>
#include <TGeoManager.h>
#include <TGeoMatrix.h>
#include <TGeoPhysicalNode.h>
#include <TGeoVolume.h>
#include <TGeoTube.h>
#include <TGeoXtru.h>
#include <TLorentzVector.h>
#include <TString.h>
#include <TVirtualMC.h>

#include "AliITSU.h"
#include "AliITSUHit.h"
#include "AliLog.h"
#include "AliMC.h"
#include "AliMagF.h"
#include "AliRun.h"
#include "AliTrackReference.h"
#include "AliITSv11Geometry.h"
#include "AliITSUv0Layer.h"
#include "AliITSUv0.h"
#include "AliITSUGeomTGeo.h"
#include "AliGeomManager.h"
using namespace TMath;


ClassImp(AliITSUv0)

//______________________________________________________________________
AliITSUv0::AliITSUv0()
:  fNWrapVol(0)
  ,fWrapRMin(0)
  ,fWrapRMax(0)
  ,fWrapZSpan(0)
  ,fLay2WrapV(0)
  ,fLayTurbo(0)
  ,fLayPhi0(0)
  ,fLayRadii(0)
  ,fLayZLength(0)
  ,fStavPerLay(0)
  ,fModPerStav(0)
  ,fStaThick(0)
  ,fStaWidth(0)
  ,fStaTilt(0)
  ,fDetThick(0)
  ,fChipTypeID(0)
  ,fBuildLevel(0)
  ,fUpGeom(0)
  ,fStaveModel(kModel0)
{
  //    Standard default constructor
  // Inputs:
  //   none.
  // Outputs:
  //   none.
  // Return:
  //   none.
}

//______________________________________________________________________
AliITSUv0::AliITSUv0(const char *title, Int_t nlay)
  :AliITSU(title,nlay)
  ,fNWrapVol(0)
  ,fWrapRMin(0)
  ,fWrapRMax(0)
  ,fWrapZSpan(0)
  ,fLay2WrapV(0)
  ,fLayTurbo(0)
  ,fLayPhi0(0)
  ,fLayRadii(0)
  ,fLayZLength(0)
  ,fStavPerLay(0)
  ,fModPerStav(0)
  ,fStaThick(0)
  ,fStaWidth(0)
  ,fStaTilt(0)
  ,fDetThick(0)
  ,fChipTypeID(0)
  ,fBuildLevel(0)
  ,fUpGeom(0)
  ,fStaveModel(kModel0)
{
  //    Standard constructor for the Upgrade geometry.
  // Inputs:
  //   const char * name   Ignored, set to "ITS"
  //   const char * title  Arbitrary title
  //   const Int_t nlay    Number of layers
  //
  fLayerName = new TString[fNLayers];
  //
  for (Int_t j=0; j<fNLayers; j++)
    fLayerName[j].Form("%s%d",AliITSUGeomTGeo::GetITSSensorPattern(),j); // See AliITSUv0Layer
  //
  fLayTurbo   = new Bool_t[fNLayers];
  fLayPhi0    = new Double_t[fNLayers];
  fLayRadii   = new Double_t[fNLayers];
  fLayZLength = new Double_t[fNLayers];
  fStavPerLay = new Int_t[fNLayers];
  fModPerStav = new Int_t[fNLayers];
  fStaThick   = new Double_t[fNLayers];
  fStaWidth   = new Double_t[fNLayers];
  fStaTilt    = new Double_t[fNLayers];
  fDetThick   = new Double_t[fNLayers];
  fChipTypeID  = new UInt_t[fNLayers];
  fBuildLevel = new Int_t[fNLayers];


  fUpGeom = new AliITSUv0Layer*[fNLayers];
  
  if (fNLayers > 0) { // if not, we'll Fatal-ize in CreateGeometry
    for (Int_t j=0; j<fNLayers; j++) {
      fLayPhi0[j] = 0;
      fLayRadii[j] = 0.;
      fLayZLength[j] = 0.;
      fStavPerLay[j] = 0;
      fModPerStav[j] = 0;
      fStaWidth[j] = 0.;
      fDetThick[j] = 0.;
      fChipTypeID[j] = 0;
      fBuildLevel[j] = 0;
      fUpGeom[j] = 0;     
    }
  }
}

//______________________________________________________________________
AliITSUv0::~AliITSUv0() {
  //    Standard destructor
  // Inputs:
  //   none.
  // Outputs:
  //   none.
  // Return:
  //   none.
  delete [] fLayTurbo;
  delete [] fLayPhi0;
  delete [] fLayRadii;
  delete [] fLayZLength;
  delete [] fStavPerLay;
  delete [] fModPerStav;
  delete [] fStaThick;
  delete [] fStaWidth;
  delete [] fStaTilt;
  delete [] fDetThick;
  delete [] fChipTypeID;
  delete [] fBuildLevel;
  delete [] fUpGeom;
  delete [] fWrapRMin;
  delete [] fWrapRMax;
  delete [] fWrapZSpan;
  delete [] fLay2WrapV;
  //
}

//______________________________________________________________________
void AliITSUv0::AddAlignableVolumes() const{
  // Creates entries for alignable volumes associating the symbolic volume
  // name with the corresponding volume path.
  // 
  // Records in the alignable entries the transformation matrices converting
  // TGeo local coordinates (in the RS of alignable volumes) to the tracking
  // system
  // For this, this function has to run before the misalignment because we
  // are using the ideal positions in the AliITSgeom object.
  // Inputs:
  //   none.
  // Outputs:
  //   none.
  // Return:
  //   none.

  AliInfo("Add ITS alignable volumes");

  if (!gGeoManager) { AliFatal("TGeoManager doesn't exist !"); return;  }
  TString pth;
  //
  pth = Form("ALIC_1/%s_2",AliITSUGeomTGeo::GetITSVolPattern());
  // RS: to be checked with MS
  if( !gGeoManager->SetAlignableEntry(AliITSUGeomTGeo::ComposeSymNameITS(),pth.Data()) )
    AliFatal(Form("Unable to set alignable entry ! %s :: %s","ITS",pth.Data()));    
  //
  int chipNum = 0;
  //
  for (int lr=0; lr<fNLayers; lr++) {
    //
    TString wrpV = fLay2WrapV[lr]!=-1 ? Form("%s%d_1/",AliITSUGeomTGeo::GetITSWrapVolPattern(),fLay2WrapV[lr]) : "";
    pth = Form("ALIC_1/%s_2/%s%s%d_1",AliITSUGeomTGeo::GetITSVolPattern(),wrpV.Data(),AliITSUGeomTGeo::GetITSLayerPattern(),lr);
    //printf("SetAlignable: %s %s\n",snm.Data(),pth.Data());
    if( !gGeoManager->SetAlignableEntry(AliITSUGeomTGeo::ComposeSymNameLayer(lr),pth.Data()) ) {
      AliFatal(Form("Unable to set alignable entry ! %s :: %s",AliITSUGeomTGeo::ComposeSymNameLayer(lr),pth.Data()));
    }
    //
    for (int ld=0; ld<fStavPerLay[lr]; ld++) {
      //
      TString pthL = Form("%s/%s%d_%d",pth.Data(),AliITSUGeomTGeo::GetITSStavePattern(),lr,ld);
      //printf("SetAlignable: %s %s\n",snmL.Data(),pthL.Data());
      if ( !gGeoManager->SetAlignableEntry(AliITSUGeomTGeo::ComposeSymNameStave(lr,ld),pthL.Data()) ) {
	AliFatal(Form("Unable to set alignable entry ! %s :: %s",AliITSUGeomTGeo::ComposeSymNameStave(lr,ld),pthL.Data()));
      }
      //
      for (int md=0; md<fModPerStav[lr]; md++) {
	//
	TString pthM = Form("%s/%s%d_%d",pthL.Data(),AliITSUGeomTGeo::GetITSChipPattern(),lr,md);
	//
	// RS: Attention, this is a hack: AliGeomManager cannot accomodate all ITSU chips w/o
	// conflicts with TPC. For this reason we define the UID of the chip to be simply its ID
	//	int modUID = AliGeomManager::LayerToVolUID(lr+1,chipNum++); // here chipNum would be chip within the layer
	int modUID = AliITSUGeomTGeo::ChipVolUID( chipNum++ );
	// 
	if ( !gGeoManager->SetAlignableEntry(AliITSUGeomTGeo::ComposeSymNameChip(lr,ld,-1,-1,md),pthM.Data(),modUID) ) {
	  AliFatal(Form("Unable to set alignable entry ! %s :: %s",AliITSUGeomTGeo::ComposeSymNameChip(lr,ld,-1,-1,md),pthM.Data()));
	}
	//
      }
    }
  }
  //
}

//______________________________________________________________________
void AliITSUv0::SetNWrapVolumes(Int_t n)
{
  // book arrays for wrapper volumes
  if (fNWrapVol) AliFatal(Form("%d wrapper volumes already defined",fNWrapVol));
  if (n<1) return;
  fNWrapVol = n;
  fWrapRMin = new Double_t[fNWrapVol];
  fWrapRMax = new Double_t[fNWrapVol];
  fWrapZSpan= new Double_t[fNWrapVol];
  for (int i=fNWrapVol;i--;) fWrapRMin[i]=fWrapRMax[i]=fWrapZSpan[i]=-1;
  //
}

//______________________________________________________________________
void AliITSUv0::DefineWrapVolume(Int_t id, Double_t rmin,Double_t rmax, Double_t zspan)
{
  // set parameters of id-th wrapper volume
  if (id>=fNWrapVol||id<0) AliFatal(Form("id=%d of wrapper volume is not in 0-%d range",id,fNWrapVol-1));
  fWrapRMin[id] = rmin;
  fWrapRMax[id] = rmax;
  fWrapZSpan[id] = zspan;
}

//______________________________________________________________________
void AliITSUv0::CreateGeometry() {

  // Create the geometry and insert it in the mother volume ITSV
  TGeoManager *geoManager = gGeoManager;

  TGeoVolume *vALIC = geoManager->GetVolume("ALIC");

  new TGeoVolumeAssembly(AliITSUGeomTGeo::GetITSVolPattern());
  TGeoVolume *vITSV = geoManager->GetVolume(AliITSUGeomTGeo::GetITSVolPattern());
  vITSV->SetUniqueID(AliITSUGeomTGeo::GetUIDShift()); // store modID -> midUUID bitshift
  vALIC->AddNode(vITSV, 2, 0);  // Copy number is 2 to cheat AliGeoManager::CheckSymNamesLUT
  //
  const Int_t kLength=100;
  Char_t vstrng[kLength] = "xxxRS"; //?
  vITSV->SetTitle(vstrng);
  //
  // Check that we have all needed parameters
  if (fNLayers <= 0) AliFatal(Form("Wrong number of layers (%d)",fNLayers));
  //
  for (Int_t j=0; j<fNLayers; j++) {
    if (fLayRadii[j] <= 0)                 AliFatal(Form("Wrong layer radius for layer %d (%f)",j,fLayRadii[j]));
    if (fLayZLength[j] <= 0)               AliFatal(Form("Wrong layer length for layer %d (%f)",j,fLayZLength[j]));
    if (fStavPerLay[j] <= 0)               AliFatal(Form("Wrong number of staves for layer %d (%d)",j,fStavPerLay[j]));
    if (fModPerStav[j] <= 0)               AliFatal(Form("Wrong number of chips for layer %d (%d)",j,fModPerStav[j]));
    if (fStaThick[j] < 0)                  AliFatal(Form("Wrong stave thickness for layer %d (%f)",j,fStaThick[j]));
    if (fLayTurbo[j] && fStaWidth[j] <= 0) AliFatal(Form("Wrong stave width for layer %d (%f)",j,fStaWidth[j]));
    if (fDetThick[j] < 0)                  AliFatal(Form("Wrong chip thickness for layer %d (%f)",j,fDetThick[j]));
    //
    if (j > 0) {
      if (fLayRadii[j]<=fLayRadii[j-1])    AliFatal(Form("Layer %d radius (%f) is smaller than layer %d radius (%f)",
							 j,fLayRadii[j],j-1,fLayRadii[j-1]));
    } // if (j > 0)

    if (fStaThick[j] == 0) AliInfo(Form("Stave thickness for layer %d not set, using default",j));
    if (fDetThick[j] == 0) AliInfo(Form("Chip thickness for layer %d not set, using default",j));

  } // for (Int_t j=0; j<fNLayers; j++)

  // Create the wrapper volumes
  TGeoVolume **wrapVols = 0;
  if (fNWrapVol) {
    wrapVols = new TGeoVolume*[fNWrapVol];
    for (int id=0;id<fNWrapVol;id++) {
      wrapVols[id] = CreateWrapperVolume(id);
      vITSV->AddNode(wrapVols[id], 1, 0);
    }
  }
  //
  fLay2WrapV = new Int_t[fNLayers];

  // Now create the actual geometry
  for (Int_t j=0; j<fNLayers; j++) {
    TGeoVolume* dest = vITSV;
    //
    if (fLayTurbo[j]) {
      fUpGeom[j] = new AliITSUv0Layer(j,kTRUE,kFALSE);
      fUpGeom[j]->SetStaveWidth(fStaWidth[j]);
      fUpGeom[j]->SetStaveTilt(fStaTilt[j]);
    }
    else fUpGeom[j] = new AliITSUv0Layer(j,kFALSE);
    //
    fUpGeom[j]->SetPhi0(fLayPhi0[j]);
    fUpGeom[j]->SetRadius(fLayRadii[j]);
    fUpGeom[j]->SetZLength(fLayZLength[j]);
    fUpGeom[j]->SetNStaves(fStavPerLay[j]);
    fUpGeom[j]->SetNChips(fModPerStav[j]);
    fUpGeom[j]->SetChipType(fChipTypeID[j]);
    fUpGeom[j]->SetBuildLevel(fBuildLevel[j]);
    fUpGeom[j]->SetStaveModel(fStaveModel);
    AliDebug(1,Form("fBuildLevel: %d\n",fBuildLevel[j]));
    //
    if (fStaThick[j] != 0) fUpGeom[j]->SetStaveThick(fStaThick[j]);
    if (fDetThick[j] != 0) fUpGeom[j]->SetSensorThick(fDetThick[j]);
    //
    for (int iw=0;iw<fNWrapVol;iw++) {
      if (fLayRadii[j]>fWrapRMin[iw] && fLayRadii[j]<fWrapRMax[iw]) {
	AliInfo(Form("Will embed layer %d in wrapper volume %d",j,iw));
	if (fLayZLength[j]>=fWrapZSpan[iw]) AliFatal(Form("ZSpan %.3f of wrapper volume %d is less than ZSpan %.3f of layer %d",
							  fWrapZSpan[iw],iw,fLayZLength[j],j));
	dest = wrapVols[iw];
	fLay2WrapV[j] = iw;
	break;
      }
    }
    fUpGeom[j]->CreateLayer(dest);
  }
  delete[] wrapVols; // delete pointer only, not the volumes
  //
}

//______________________________________________________________________
void AliITSUv0::CreateMaterials() {
  // Create ITS materials
  //     This function defines the default materials used in the Geant
  // Monte Carlo simulations for the geometries AliITSv1, AliITSv3,
  // AliITSv11Hybrid.
  // In general it is automatically replaced by
  // the CreateMaterials routine defined in AliITSv?. Should the function
  // CreateMaterials not exist for the geometry version you are using this
  // one is used. See the definition found in AliITSv5 or the other routine
  // for a complete definition.
  // Inputs:
  //   none.
  // Outputs:
  //   none.
  // Return:
  //   none.

  Int_t   ifield = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
  Float_t fieldm = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->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"

  // AIR
  Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
  Float_t zAir[4]={6.,7.,8.,18.};
  Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
  Float_t dAir = 1.20479E-3;

  // Water
  Float_t aWater[2]={1.00794,15.9994};
  Float_t zWater[2]={1.,8.};
  Float_t wWater[2]={0.111894,0.888106};
  Float_t dWater   = 1.0;


  // Kapton
  Float_t aKapton[4]={1.00794,12.0107, 14.010,15.9994};
  Float_t zKapton[4]={1.,6.,7.,8.};
  Float_t wKapton[4]={0.026362,0.69113,0.07327,0.209235};
  Float_t dKapton   = 1.42;
 
  AliMixture(1,"AIR$",aAir,zAir,dAir,4,wAir);
  AliMedium(1, "AIR$",1,0,ifield,fieldm,tmaxfdAir,stemaxAir,deemaxAir,epsilAir,stminAir);

  AliMixture(2,"WATER$",aWater,zWater,dWater,2,wWater);
  AliMedium(2, "WATER$",2,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);

  AliMaterial(3,"SI$",0.28086E+02,0.14000E+02,0.23300E+01,0.93600E+01,0.99900E+03);
  AliMedium(3,  "SI$",3,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);

  AliMaterial(4,"BERILLIUM$",9.01, 4., 1.848, 35.3, 36.7);// From AliPIPEv3
  AliMedium(4,  "BERILLIUM$",4,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);

  AliMaterial(5,"COPPER$",0.63546E+02,0.29000E+02,0.89600E+01,0.14300E+01,0.99900E+03);
  AliMedium(5,  "COPPER$",5,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);

    
  // needed for STAVE , Carbon, kapton, Epoxy, flexcable

  //AliMaterial(6,"CARBON$",12.0107,6,2.210,999,999);
  AliMaterial(6,"CARBON$",12.0107,6,2.210/1.3,999,999);
  AliMedium(6,  "CARBON$",6,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);

  AliMixture(7,"KAPTON(POLYCH2)$", aKapton, zKapton, dKapton, 4, wKapton);
  AliMedium(7, "KAPTON(POLYCH2)$",7,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);


 
  // values below modified as compared to source AliITSv11 !

  //AliMaterial(7,"GLUE$",0.12011E+02,0.60000E+01,0.1930E+01/2.015,999,999); // original
  AliMaterial(7,"GLUE$",12.011,6,1.93/2.015,999,999);  // conform with ATLAS, Corrado, Stefan
  AliMedium(7,  "GLUE$",7,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);

  // All types of carbon
  // Unidirectional prepreg
 AliMaterial(8,"K13D2U2k$",12.0107,6,1.643,999,999);
 AliMedium(8,  "K13D2U2k$",8,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
 //Impregnated thread
 AliMaterial(9,"M60J3K$",12.0107,6,2.21,999,999);
 AliMedium(9,  "M60J3K$",9,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
 //Impregnated thread
 AliMaterial(10,"M55J6K$",12.0107,6,1.63,999,999);
 AliMedium(10,  "M55J6K$",10,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
 // Fabric(0/90)
 AliMaterial(11,"T300$",12.0107,6,1.725,999,999);
 AliMedium(11,  "T300$",11,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
 //AMEC Thermasol
 AliMaterial(12,"FGS003$",12.0107,6,1.6,999,999);
 AliMedium(12,  "FGS003$",12,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);
 // Carbon fleece
 AliMaterial(13,"CarbonFleece$",12.0107,6,0.4,999,999);
 AliMedium(13,  "CarbonFleece$",13,0,ifield,fieldm,tmaxfdSi,stemaxSi,deemaxSi,epsilSi,stminSi);

  // Flex cable
  Float_t aFCm[5]={12.0107,1.00794,14.0067,15.9994,26.981538};
  Float_t zFCm[5]={6.,1.,7.,8.,13.};
  Float_t wFCm[5]={0.520088819984,0.01983871336,0.0551367996,0.157399667056, 0.247536};
  //Float_t dFCm = 1.6087;  // original
  //Float_t dFCm = 2.55;   // conform with STAR
   Float_t dFCm = 2.595;   // conform with Corrado

  AliMixture(14,"FLEXCABLE$",aFCm,zFCm,dFCm,5,wFCm);
  AliMedium(14, "FLEXCABLE$",14,0,ifield,fieldm,tmaxfd,stemax,deemax,epsil,stmin);

}

//______________________________________________________________________
void AliITSUv0::DefineLayer(Int_t nlay, double phi0, Double_t r,
			    Double_t zlen, Int_t nstav,
			    Int_t nmod, Double_t lthick,
			    Double_t dthick, UInt_t dettypeID)
{
  //     Sets the layer parameters
  // Inputs:
  //          nlay    layer number
  //          phi0    layer phi0
  //          r       layer radius
  //          zlen    layer length
  //          nstav   number of staves
  //          nmod    number of chips per stave
  //          lthick  stave thickness (if omitted, defaults to 0)
  //          dthick  detector thickness (if omitted, defaults to 0)
  // Outputs:
  //   none.
  // Return:
  //   none.
  
  if (nlay >= fNLayers || nlay < 0) {
    AliError(Form("Wrong layer number (%d)",nlay));
    return;
  }
  
  fLayTurbo[nlay] = kFALSE;
  fLayPhi0[nlay] = phi0;
  fLayRadii[nlay] = r;
  fLayZLength[nlay] = zlen;
  fStavPerLay[nlay] = nstav;
  fModPerStav[nlay] = nmod;
  fStaThick[nlay] = lthick;
  fDetThick[nlay] = dthick;
  fChipTypeID[nlay] = dettypeID;
    
}

//______________________________________________________________________
void AliITSUv0::DefineLayerTurbo(Int_t nlay, Double_t phi0, Double_t r, Double_t zlen, Int_t nstav,
				 Int_t nmod, Double_t width, Double_t tilt,
				 Double_t lthick,Double_t dthick,
				 UInt_t dettypeID, Int_t buildLevel)
{
  //     Sets the layer parameters for a "turbo" layer
  //     (i.e. a layer whose staves overlap in phi)
  // Inputs:
  //          nlay    layer number
  //          phi0    phi of 1st stave
  //          r       layer radius
  //          zlen    layer length
  //          nstav   number of staves
  //          nmod    number of chips per stave
  //          width   stave width
  //          tilt    layer tilt angle (degrees)
  //          lthick  stave thickness (if omitted, defaults to 0)
  //          dthick  detector thickness (if omitted, defaults to 0)
  //          dettypeID  ??
  //          buildLevel (if 0, all geometry is build, used for material budget studies)
  // Outputs:
  //   none.
  // Return:
  //   none.

  if (nlay >= fNLayers || nlay < 0) {
    AliError(Form("Wrong layer number (%d)",nlay));
    return;
  }

  fLayTurbo[nlay] = kTRUE;
  fLayPhi0[nlay] = phi0;
  fLayRadii[nlay] = r;
  fLayZLength[nlay] = zlen;
  fStavPerLay[nlay] = nstav;
  fModPerStav[nlay] = nmod;
  fStaThick[nlay] = lthick;
  fStaWidth[nlay] = width;
  fStaTilt[nlay] = tilt;
  fDetThick[nlay] = dthick;
  fChipTypeID[nlay] = dettypeID;
  fBuildLevel[nlay] = buildLevel;

}

//______________________________________________________________________
void AliITSUv0::GetLayerParameters(Int_t nlay, Double_t &phi0,
				   Double_t &r, Double_t &zlen,
				   Int_t &nstav, Int_t &nmod,
				   Double_t &width, Double_t &tilt,
				   Double_t &lthick, Double_t &dthick,
				   UInt_t &dettype) const
{
  //     Gets the layer parameters
  // Inputs:
  //          nlay    layer number
  // Outputs:
  //          phi0    phi of 1st stave
  //          r       layer radius
  //          zlen    layer length
  //          nstav   number of staves
  //          nmod    number of chips per stave
  //          width   stave width
  //          tilt    stave tilt angle
  //          lthick  stave thickness
  //          dthick  detector thickness
  //          dettype detector type
  // Return:
  //   none.

  if (nlay >= fNLayers || nlay < 0) {
    AliError(Form("Wrong layer number (%d)",nlay));
    return;
  }
  
  phi0   = fLayPhi0[nlay];
  r      = fLayRadii[nlay];
  zlen   = fLayZLength[nlay];
  nstav  = fStavPerLay[nlay];
  nmod   = fModPerStav[nlay];
  width  = fStaWidth[nlay];
  tilt   = fStaTilt[nlay];
  lthick = fStaThick[nlay];
  dthick = fDetThick[nlay];
  dettype= fChipTypeID[nlay];
}

//______________________________________________________________________
TGeoVolume* AliITSUv0::CreateWrapperVolume(Int_t id)
{
  //     Creates an air-filled wrapper cylindrical volume 
  // Inputs:
  //          volume id
  // Outputs:
  //          the wrapper volume

  if (fWrapRMin[id]<0 || fWrapRMax[id]<0 || fWrapZSpan[id]<0) AliFatal(Form("Wrapper volume %d was requested but not defined",id));
  // Now create the actual shape and volume
  //
  TGeoTube *tube = new TGeoTube(fWrapRMin[id], fWrapRMax[id], fWrapZSpan[id]/2.);

  TGeoMedium *medAir = gGeoManager->GetMedium("ITS_AIR$");

  char volnam[30];
  snprintf(volnam, 29, "%s%d", AliITSUGeomTGeo::GetITSWrapVolPattern(),id);

  TGeoVolume *wrapper = new TGeoVolume(volnam, tube, medAir);

  return wrapper;
}

//______________________________________________________________________
void AliITSUv0::Init()
{
  //     Initialise the ITS after it has been created.
  UpdateInternalGeometry();
  AliITSU::Init();
  //  
}

//______________________________________________________________________
Bool_t AliITSUv0::IsLayerTurbo(Int_t nlay)
{
  //     Returns true if the layer is a "turbo" layer
  if ( nlay < 0 || nlay > fNLayers ) {
    AliError(Form("Wrong layer number %d",nlay));
    return kFALSE;
  } 
  else return fUpGeom[nlay]->IsTurbo();
}

//______________________________________________________________________
void AliITSUv0::SetDefaults()
{
  // sets the default segmentation, response, digit and raw cluster classes
}

//______________________________________________________________________
void AliITSUv0::StepManager()
{
  //    Called for every step in the ITS, then calles the AliITSUHit 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.
  // Inputs:
  //   none.
  // Outputs:
  //   none.
  // Return:
  //   none.
  if(!(this->IsActive())) return;
  if(!(TVirtualMC::GetMC()->TrackCharge())) return;
  //
  Int_t copy, lay = 0;
  Int_t id = TVirtualMC::GetMC()->CurrentVolID(copy);

  Bool_t notSens = kFALSE;
  while ((lay<fNLayers)  && (notSens = (id!=fIdSens[lay]))) ++lay;
  //printf("R: %.1f | Lay: %d  NotSens: %d\n",positionRS.Pt(), lay, notSens);
	   
  if (notSens) return;

  if(TVirtualMC::GetMC()->IsTrackExiting()) {
    AddTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber(), AliTrackReference::kITS);
  } // if Outer ITS mother Volume

  static TLorentzVector position, momentum; // Saves on calls to construtors
  static AliITSUHit hit;// Saves on calls to constructors

  TClonesArray &lhits = *(Hits());
  Int_t   cpn0, cpn1, mod, status = 0;
  //
  // Track status
  if(TVirtualMC::GetMC()->IsTrackInside())      status +=  1;
  if(TVirtualMC::GetMC()->IsTrackEntering())    status +=  2;
  if(TVirtualMC::GetMC()->IsTrackExiting())     status +=  4;
  if(TVirtualMC::GetMC()->IsTrackOut())         status +=  8;
  if(TVirtualMC::GetMC()->IsTrackDisappeared()) status += 16;
  if(TVirtualMC::GetMC()->IsTrackStop())        status += 32;
  if(TVirtualMC::GetMC()->IsTrackAlive())       status += 64;

  //
  // retrieve the indices with the volume path
  //
  if (lay < 0 || lay >= fNLayers) {
    AliError(Form("Invalid value: lay=%d. Not an ITS sensitive volume",lay));
    return; // not an ITS sensitive volume.
  } else {
    copy = 1;
    TVirtualMC::GetMC()->CurrentVolOffID(1,cpn1);
    TVirtualMC::GetMC()->CurrentVolOffID(2,cpn0);
  } //

  mod = fGeomTGeo->GetChipIndex(lay,cpn0,cpn1);
  //RS2DEL  fInitGeom.DecodeDetector(mod,lay+1,cpn0,cpn1,copy);
  //
  // Fill hit structure.
  //
  hit.SetChip(mod);
  hit.SetTrack(gAlice->GetMCApp()->GetCurrentTrackNumber());
  TVirtualMC::GetMC()->TrackPosition(position);
  TVirtualMC::GetMC()->TrackMomentum(momentum);
  hit.SetPosition(position);
  hit.SetTime(TVirtualMC::GetMC()->TrackTime());
  hit.SetMomentum(momentum);
  hit.SetStatus(status);
  hit.SetEdep(TVirtualMC::GetMC()->Edep());
  hit.SetShunt(GetIshunt());
  if(TVirtualMC::GetMC()->IsTrackEntering()){
    hit.SetStartPosition(position);
    hit.SetStartTime(TVirtualMC::GetMC()->TrackTime());
    hit.SetStartStatus(status);
    return; // don't save entering hit.
  } // end if IsEntering
    // Fill hit structure with this new hit.
    //Info("StepManager","Calling Copy Constructor");
  new(lhits[fNhits++]) AliITSUHit(hit); // Use Copy Construtor.
  // Save old position... for next hit.
  hit.SetStartPosition(position);
  hit.SetStartTime(TVirtualMC::GetMC()->TrackTime());
  hit.SetStartStatus(status);

  return;
}

//______________________________________________________________________
void AliITSUv0::SetLayerChipTypeID(Int_t lr, UInt_t id)
{
  // set det type
  if (!fChipTypeID || fNLayers<=lr) AliFatal(Form("Number of layers %d, %d is manipulated",fNLayers,lr));
  fChipTypeID[lr] = id;
}

//______________________________________________________________________
Int_t AliITSUv0::GetLayerChipTypeID(Int_t lr)
{
  // set det type
  if (!fChipTypeID || fNLayers<=lr) AliFatal(Form("Number of layers %d, %d is manipulated",fNLayers,lr));
  return fChipTypeID[lr];
}
 AliITSUv0.cxx:1
 AliITSUv0.cxx:2
 AliITSUv0.cxx:3
 AliITSUv0.cxx:4
 AliITSUv0.cxx:5
 AliITSUv0.cxx:6
 AliITSUv0.cxx:7
 AliITSUv0.cxx:8
 AliITSUv0.cxx:9
 AliITSUv0.cxx:10
 AliITSUv0.cxx:11
 AliITSUv0.cxx:12
 AliITSUv0.cxx:13
 AliITSUv0.cxx:14
 AliITSUv0.cxx:15
 AliITSUv0.cxx:16
 AliITSUv0.cxx:17
 AliITSUv0.cxx:18
 AliITSUv0.cxx:19
 AliITSUv0.cxx:20
 AliITSUv0.cxx:21
 AliITSUv0.cxx:22
 AliITSUv0.cxx:23
 AliITSUv0.cxx:24
 AliITSUv0.cxx:25
 AliITSUv0.cxx:26
 AliITSUv0.cxx:27
 AliITSUv0.cxx:28
 AliITSUv0.cxx:29
 AliITSUv0.cxx:30
 AliITSUv0.cxx:31
 AliITSUv0.cxx:32
 AliITSUv0.cxx:33
 AliITSUv0.cxx:34
 AliITSUv0.cxx:35
 AliITSUv0.cxx:36
 AliITSUv0.cxx:37
 AliITSUv0.cxx:38
 AliITSUv0.cxx:39
 AliITSUv0.cxx:40
 AliITSUv0.cxx:41
 AliITSUv0.cxx:42
 AliITSUv0.cxx:43
 AliITSUv0.cxx:44
 AliITSUv0.cxx:45
 AliITSUv0.cxx:46
 AliITSUv0.cxx:47
 AliITSUv0.cxx:48
 AliITSUv0.cxx:49
 AliITSUv0.cxx:50
 AliITSUv0.cxx:51
 AliITSUv0.cxx:52
 AliITSUv0.cxx:53
 AliITSUv0.cxx:54
 AliITSUv0.cxx:55
 AliITSUv0.cxx:56
 AliITSUv0.cxx:57
 AliITSUv0.cxx:58
 AliITSUv0.cxx:59
 AliITSUv0.cxx:60
 AliITSUv0.cxx:61
 AliITSUv0.cxx:62
 AliITSUv0.cxx:63
 AliITSUv0.cxx:64
 AliITSUv0.cxx:65
 AliITSUv0.cxx:66
 AliITSUv0.cxx:67
 AliITSUv0.cxx:68
 AliITSUv0.cxx:69
 AliITSUv0.cxx:70
 AliITSUv0.cxx:71
 AliITSUv0.cxx:72
 AliITSUv0.cxx:73
 AliITSUv0.cxx:74
 AliITSUv0.cxx:75
 AliITSUv0.cxx:76
 AliITSUv0.cxx:77
 AliITSUv0.cxx:78
 AliITSUv0.cxx:79
 AliITSUv0.cxx:80
 AliITSUv0.cxx:81
 AliITSUv0.cxx:82
 AliITSUv0.cxx:83
 AliITSUv0.cxx:84
 AliITSUv0.cxx:85
 AliITSUv0.cxx:86
 AliITSUv0.cxx:87
 AliITSUv0.cxx:88
 AliITSUv0.cxx:89
 AliITSUv0.cxx:90
 AliITSUv0.cxx:91
 AliITSUv0.cxx:92
 AliITSUv0.cxx:93
 AliITSUv0.cxx:94
 AliITSUv0.cxx:95
 AliITSUv0.cxx:96
 AliITSUv0.cxx:97
 AliITSUv0.cxx:98
 AliITSUv0.cxx:99
 AliITSUv0.cxx:100
 AliITSUv0.cxx:101
 AliITSUv0.cxx:102
 AliITSUv0.cxx:103
 AliITSUv0.cxx:104
 AliITSUv0.cxx:105
 AliITSUv0.cxx:106
 AliITSUv0.cxx:107
 AliITSUv0.cxx:108
 AliITSUv0.cxx:109
 AliITSUv0.cxx:110
 AliITSUv0.cxx:111
 AliITSUv0.cxx:112
 AliITSUv0.cxx:113
 AliITSUv0.cxx:114
 AliITSUv0.cxx:115
 AliITSUv0.cxx:116
 AliITSUv0.cxx:117
 AliITSUv0.cxx:118
 AliITSUv0.cxx:119
 AliITSUv0.cxx:120
 AliITSUv0.cxx:121
 AliITSUv0.cxx:122
 AliITSUv0.cxx:123
 AliITSUv0.cxx:124
 AliITSUv0.cxx:125
 AliITSUv0.cxx:126
 AliITSUv0.cxx:127
 AliITSUv0.cxx:128
 AliITSUv0.cxx:129
 AliITSUv0.cxx:130
 AliITSUv0.cxx:131
 AliITSUv0.cxx:132
 AliITSUv0.cxx:133
 AliITSUv0.cxx:134
 AliITSUv0.cxx:135
 AliITSUv0.cxx:136
 AliITSUv0.cxx:137
 AliITSUv0.cxx:138
 AliITSUv0.cxx:139
 AliITSUv0.cxx:140
 AliITSUv0.cxx:141
 AliITSUv0.cxx:142
 AliITSUv0.cxx:143
 AliITSUv0.cxx:144
 AliITSUv0.cxx:145
 AliITSUv0.cxx:146
 AliITSUv0.cxx:147
 AliITSUv0.cxx:148
 AliITSUv0.cxx:149
 AliITSUv0.cxx:150
 AliITSUv0.cxx:151
 AliITSUv0.cxx:152
 AliITSUv0.cxx:153
 AliITSUv0.cxx:154
 AliITSUv0.cxx:155
 AliITSUv0.cxx:156
 AliITSUv0.cxx:157
 AliITSUv0.cxx:158
 AliITSUv0.cxx:159
 AliITSUv0.cxx:160
 AliITSUv0.cxx:161
 AliITSUv0.cxx:162
 AliITSUv0.cxx:163
 AliITSUv0.cxx:164
 AliITSUv0.cxx:165
 AliITSUv0.cxx:166
 AliITSUv0.cxx:167
 AliITSUv0.cxx:168
 AliITSUv0.cxx:169
 AliITSUv0.cxx:170
 AliITSUv0.cxx:171
 AliITSUv0.cxx:172
 AliITSUv0.cxx:173
 AliITSUv0.cxx:174
 AliITSUv0.cxx:175
 AliITSUv0.cxx:176
 AliITSUv0.cxx:177
 AliITSUv0.cxx:178
 AliITSUv0.cxx:179
 AliITSUv0.cxx:180
 AliITSUv0.cxx:181
 AliITSUv0.cxx:182
 AliITSUv0.cxx:183
 AliITSUv0.cxx:184
 AliITSUv0.cxx:185
 AliITSUv0.cxx:186
 AliITSUv0.cxx:187
 AliITSUv0.cxx:188
 AliITSUv0.cxx:189
 AliITSUv0.cxx:190
 AliITSUv0.cxx:191
 AliITSUv0.cxx:192
 AliITSUv0.cxx:193
 AliITSUv0.cxx:194
 AliITSUv0.cxx:195
 AliITSUv0.cxx:196
 AliITSUv0.cxx:197
 AliITSUv0.cxx:198
 AliITSUv0.cxx:199
 AliITSUv0.cxx:200
 AliITSUv0.cxx:201
 AliITSUv0.cxx:202
 AliITSUv0.cxx:203
 AliITSUv0.cxx:204
 AliITSUv0.cxx:205
 AliITSUv0.cxx:206
 AliITSUv0.cxx:207
 AliITSUv0.cxx:208
 AliITSUv0.cxx:209
 AliITSUv0.cxx:210
 AliITSUv0.cxx:211
 AliITSUv0.cxx:212
 AliITSUv0.cxx:213
 AliITSUv0.cxx:214
 AliITSUv0.cxx:215
 AliITSUv0.cxx:216
 AliITSUv0.cxx:217
 AliITSUv0.cxx:218
 AliITSUv0.cxx:219
 AliITSUv0.cxx:220
 AliITSUv0.cxx:221
 AliITSUv0.cxx:222
 AliITSUv0.cxx:223
 AliITSUv0.cxx:224
 AliITSUv0.cxx:225
 AliITSUv0.cxx:226
 AliITSUv0.cxx:227
 AliITSUv0.cxx:228
 AliITSUv0.cxx:229
 AliITSUv0.cxx:230
 AliITSUv0.cxx:231
 AliITSUv0.cxx:232
 AliITSUv0.cxx:233
 AliITSUv0.cxx:234
 AliITSUv0.cxx:235
 AliITSUv0.cxx:236
 AliITSUv0.cxx:237
 AliITSUv0.cxx:238
 AliITSUv0.cxx:239
 AliITSUv0.cxx:240
 AliITSUv0.cxx:241
 AliITSUv0.cxx:242
 AliITSUv0.cxx:243
 AliITSUv0.cxx:244
 AliITSUv0.cxx:245
 AliITSUv0.cxx:246
 AliITSUv0.cxx:247
 AliITSUv0.cxx:248
 AliITSUv0.cxx:249
 AliITSUv0.cxx:250
 AliITSUv0.cxx:251
 AliITSUv0.cxx:252
 AliITSUv0.cxx:253
 AliITSUv0.cxx:254
 AliITSUv0.cxx:255
 AliITSUv0.cxx:256
 AliITSUv0.cxx:257
 AliITSUv0.cxx:258
 AliITSUv0.cxx:259
 AliITSUv0.cxx:260
 AliITSUv0.cxx:261
 AliITSUv0.cxx:262
 AliITSUv0.cxx:263
 AliITSUv0.cxx:264
 AliITSUv0.cxx:265
 AliITSUv0.cxx:266
 AliITSUv0.cxx:267
 AliITSUv0.cxx:268
 AliITSUv0.cxx:269
 AliITSUv0.cxx:270
 AliITSUv0.cxx:271
 AliITSUv0.cxx:272
 AliITSUv0.cxx:273
 AliITSUv0.cxx:274
 AliITSUv0.cxx:275
 AliITSUv0.cxx:276
 AliITSUv0.cxx:277
 AliITSUv0.cxx:278
 AliITSUv0.cxx:279
 AliITSUv0.cxx:280
 AliITSUv0.cxx:281
 AliITSUv0.cxx:282
 AliITSUv0.cxx:283
 AliITSUv0.cxx:284
 AliITSUv0.cxx:285
 AliITSUv0.cxx:286
 AliITSUv0.cxx:287
 AliITSUv0.cxx:288
 AliITSUv0.cxx:289
 AliITSUv0.cxx:290
 AliITSUv0.cxx:291
 AliITSUv0.cxx:292
 AliITSUv0.cxx:293
 AliITSUv0.cxx:294
 AliITSUv0.cxx:295
 AliITSUv0.cxx:296
 AliITSUv0.cxx:297
 AliITSUv0.cxx:298
 AliITSUv0.cxx:299
 AliITSUv0.cxx:300
 AliITSUv0.cxx:301
 AliITSUv0.cxx:302
 AliITSUv0.cxx:303
 AliITSUv0.cxx:304
 AliITSUv0.cxx:305
 AliITSUv0.cxx:306
 AliITSUv0.cxx:307
 AliITSUv0.cxx:308
 AliITSUv0.cxx:309
 AliITSUv0.cxx:310
 AliITSUv0.cxx:311
 AliITSUv0.cxx:312
 AliITSUv0.cxx:313
 AliITSUv0.cxx:314
 AliITSUv0.cxx:315
 AliITSUv0.cxx:316
 AliITSUv0.cxx:317
 AliITSUv0.cxx:318
 AliITSUv0.cxx:319
 AliITSUv0.cxx:320
 AliITSUv0.cxx:321
 AliITSUv0.cxx:322
 AliITSUv0.cxx:323
 AliITSUv0.cxx:324
 AliITSUv0.cxx:325
 AliITSUv0.cxx:326
 AliITSUv0.cxx:327
 AliITSUv0.cxx:328
 AliITSUv0.cxx:329
 AliITSUv0.cxx:330
 AliITSUv0.cxx:331
 AliITSUv0.cxx:332
 AliITSUv0.cxx:333
 AliITSUv0.cxx:334
 AliITSUv0.cxx:335
 AliITSUv0.cxx:336
 AliITSUv0.cxx:337
 AliITSUv0.cxx:338
 AliITSUv0.cxx:339
 AliITSUv0.cxx:340
 AliITSUv0.cxx:341
 AliITSUv0.cxx:342
 AliITSUv0.cxx:343
 AliITSUv0.cxx:344
 AliITSUv0.cxx:345
 AliITSUv0.cxx:346
 AliITSUv0.cxx:347
 AliITSUv0.cxx:348
 AliITSUv0.cxx:349
 AliITSUv0.cxx:350
 AliITSUv0.cxx:351
 AliITSUv0.cxx:352
 AliITSUv0.cxx:353
 AliITSUv0.cxx:354
 AliITSUv0.cxx:355
 AliITSUv0.cxx:356
 AliITSUv0.cxx:357
 AliITSUv0.cxx:358
 AliITSUv0.cxx:359
 AliITSUv0.cxx:360
 AliITSUv0.cxx:361
 AliITSUv0.cxx:362
 AliITSUv0.cxx:363
 AliITSUv0.cxx:364
 AliITSUv0.cxx:365
 AliITSUv0.cxx:366
 AliITSUv0.cxx:367
 AliITSUv0.cxx:368
 AliITSUv0.cxx:369
 AliITSUv0.cxx:370
 AliITSUv0.cxx:371
 AliITSUv0.cxx:372
 AliITSUv0.cxx:373
 AliITSUv0.cxx:374
 AliITSUv0.cxx:375
 AliITSUv0.cxx:376
 AliITSUv0.cxx:377
 AliITSUv0.cxx:378
 AliITSUv0.cxx:379
 AliITSUv0.cxx:380
 AliITSUv0.cxx:381
 AliITSUv0.cxx:382
 AliITSUv0.cxx:383
 AliITSUv0.cxx:384
 AliITSUv0.cxx:385
 AliITSUv0.cxx:386
 AliITSUv0.cxx:387
 AliITSUv0.cxx:388
 AliITSUv0.cxx:389
 AliITSUv0.cxx:390
 AliITSUv0.cxx:391
 AliITSUv0.cxx:392
 AliITSUv0.cxx:393
 AliITSUv0.cxx:394
 AliITSUv0.cxx:395
 AliITSUv0.cxx:396
 AliITSUv0.cxx:397
 AliITSUv0.cxx:398
 AliITSUv0.cxx:399
 AliITSUv0.cxx:400
 AliITSUv0.cxx:401
 AliITSUv0.cxx:402
 AliITSUv0.cxx:403
 AliITSUv0.cxx:404
 AliITSUv0.cxx:405
 AliITSUv0.cxx:406
 AliITSUv0.cxx:407
 AliITSUv0.cxx:408
 AliITSUv0.cxx:409
 AliITSUv0.cxx:410
 AliITSUv0.cxx:411
 AliITSUv0.cxx:412
 AliITSUv0.cxx:413
 AliITSUv0.cxx:414
 AliITSUv0.cxx:415
 AliITSUv0.cxx:416
 AliITSUv0.cxx:417
 AliITSUv0.cxx:418
 AliITSUv0.cxx:419
 AliITSUv0.cxx:420
 AliITSUv0.cxx:421
 AliITSUv0.cxx:422
 AliITSUv0.cxx:423
 AliITSUv0.cxx:424
 AliITSUv0.cxx:425
 AliITSUv0.cxx:426
 AliITSUv0.cxx:427
 AliITSUv0.cxx:428
 AliITSUv0.cxx:429
 AliITSUv0.cxx:430
 AliITSUv0.cxx:431
 AliITSUv0.cxx:432
 AliITSUv0.cxx:433
 AliITSUv0.cxx:434
 AliITSUv0.cxx:435
 AliITSUv0.cxx:436
 AliITSUv0.cxx:437
 AliITSUv0.cxx:438
 AliITSUv0.cxx:439
 AliITSUv0.cxx:440
 AliITSUv0.cxx:441
 AliITSUv0.cxx:442
 AliITSUv0.cxx:443
 AliITSUv0.cxx:444
 AliITSUv0.cxx:445
 AliITSUv0.cxx:446
 AliITSUv0.cxx:447
 AliITSUv0.cxx:448
 AliITSUv0.cxx:449
 AliITSUv0.cxx:450
 AliITSUv0.cxx:451
 AliITSUv0.cxx:452
 AliITSUv0.cxx:453
 AliITSUv0.cxx:454
 AliITSUv0.cxx:455
 AliITSUv0.cxx:456
 AliITSUv0.cxx:457
 AliITSUv0.cxx:458
 AliITSUv0.cxx:459
 AliITSUv0.cxx:460
 AliITSUv0.cxx:461
 AliITSUv0.cxx:462
 AliITSUv0.cxx:463
 AliITSUv0.cxx:464
 AliITSUv0.cxx:465
 AliITSUv0.cxx:466
 AliITSUv0.cxx:467
 AliITSUv0.cxx:468
 AliITSUv0.cxx:469
 AliITSUv0.cxx:470
 AliITSUv0.cxx:471
 AliITSUv0.cxx:472
 AliITSUv0.cxx:473
 AliITSUv0.cxx:474
 AliITSUv0.cxx:475
 AliITSUv0.cxx:476
 AliITSUv0.cxx:477
 AliITSUv0.cxx:478
 AliITSUv0.cxx:479
 AliITSUv0.cxx:480
 AliITSUv0.cxx:481
 AliITSUv0.cxx:482
 AliITSUv0.cxx:483
 AliITSUv0.cxx:484
 AliITSUv0.cxx:485
 AliITSUv0.cxx:486
 AliITSUv0.cxx:487
 AliITSUv0.cxx:488
 AliITSUv0.cxx:489
 AliITSUv0.cxx:490
 AliITSUv0.cxx:491
 AliITSUv0.cxx:492
 AliITSUv0.cxx:493
 AliITSUv0.cxx:494
 AliITSUv0.cxx:495
 AliITSUv0.cxx:496
 AliITSUv0.cxx:497
 AliITSUv0.cxx:498
 AliITSUv0.cxx:499
 AliITSUv0.cxx:500
 AliITSUv0.cxx:501
 AliITSUv0.cxx:502
 AliITSUv0.cxx:503
 AliITSUv0.cxx:504
 AliITSUv0.cxx:505
 AliITSUv0.cxx:506
 AliITSUv0.cxx:507
 AliITSUv0.cxx:508
 AliITSUv0.cxx:509
 AliITSUv0.cxx:510
 AliITSUv0.cxx:511
 AliITSUv0.cxx:512
 AliITSUv0.cxx:513
 AliITSUv0.cxx:514
 AliITSUv0.cxx:515
 AliITSUv0.cxx:516
 AliITSUv0.cxx:517
 AliITSUv0.cxx:518
 AliITSUv0.cxx:519
 AliITSUv0.cxx:520
 AliITSUv0.cxx:521
 AliITSUv0.cxx:522
 AliITSUv0.cxx:523
 AliITSUv0.cxx:524
 AliITSUv0.cxx:525
 AliITSUv0.cxx:526
 AliITSUv0.cxx:527
 AliITSUv0.cxx:528
 AliITSUv0.cxx:529
 AliITSUv0.cxx:530
 AliITSUv0.cxx:531
 AliITSUv0.cxx:532
 AliITSUv0.cxx:533
 AliITSUv0.cxx:534
 AliITSUv0.cxx:535
 AliITSUv0.cxx:536
 AliITSUv0.cxx:537
 AliITSUv0.cxx:538
 AliITSUv0.cxx:539
 AliITSUv0.cxx:540
 AliITSUv0.cxx:541
 AliITSUv0.cxx:542
 AliITSUv0.cxx:543
 AliITSUv0.cxx:544
 AliITSUv0.cxx:545
 AliITSUv0.cxx:546
 AliITSUv0.cxx:547
 AliITSUv0.cxx:548
 AliITSUv0.cxx:549
 AliITSUv0.cxx:550
 AliITSUv0.cxx:551
 AliITSUv0.cxx:552
 AliITSUv0.cxx:553
 AliITSUv0.cxx:554
 AliITSUv0.cxx:555
 AliITSUv0.cxx:556
 AliITSUv0.cxx:557
 AliITSUv0.cxx:558
 AliITSUv0.cxx:559
 AliITSUv0.cxx:560
 AliITSUv0.cxx:561
 AliITSUv0.cxx:562
 AliITSUv0.cxx:563
 AliITSUv0.cxx:564
 AliITSUv0.cxx:565
 AliITSUv0.cxx:566
 AliITSUv0.cxx:567
 AliITSUv0.cxx:568
 AliITSUv0.cxx:569
 AliITSUv0.cxx:570
 AliITSUv0.cxx:571
 AliITSUv0.cxx:572
 AliITSUv0.cxx:573
 AliITSUv0.cxx:574
 AliITSUv0.cxx:575
 AliITSUv0.cxx:576
 AliITSUv0.cxx:577
 AliITSUv0.cxx:578
 AliITSUv0.cxx:579
 AliITSUv0.cxx:580
 AliITSUv0.cxx:581
 AliITSUv0.cxx:582
 AliITSUv0.cxx:583
 AliITSUv0.cxx:584
 AliITSUv0.cxx:585
 AliITSUv0.cxx:586
 AliITSUv0.cxx:587
 AliITSUv0.cxx:588
 AliITSUv0.cxx:589
 AliITSUv0.cxx:590
 AliITSUv0.cxx:591
 AliITSUv0.cxx:592
 AliITSUv0.cxx:593
 AliITSUv0.cxx:594
 AliITSUv0.cxx:595
 AliITSUv0.cxx:596
 AliITSUv0.cxx:597
 AliITSUv0.cxx:598
 AliITSUv0.cxx:599
 AliITSUv0.cxx:600
 AliITSUv0.cxx:601
 AliITSUv0.cxx:602
 AliITSUv0.cxx:603
 AliITSUv0.cxx:604
 AliITSUv0.cxx:605
 AliITSUv0.cxx:606
 AliITSUv0.cxx:607
 AliITSUv0.cxx:608
 AliITSUv0.cxx:609
 AliITSUv0.cxx:610
 AliITSUv0.cxx:611
 AliITSUv0.cxx:612
 AliITSUv0.cxx:613
 AliITSUv0.cxx:614
 AliITSUv0.cxx:615
 AliITSUv0.cxx:616
 AliITSUv0.cxx:617
 AliITSUv0.cxx:618
 AliITSUv0.cxx:619
 AliITSUv0.cxx:620
 AliITSUv0.cxx:621
 AliITSUv0.cxx:622
 AliITSUv0.cxx:623
 AliITSUv0.cxx:624
 AliITSUv0.cxx:625
 AliITSUv0.cxx:626
 AliITSUv0.cxx:627
 AliITSUv0.cxx:628
 AliITSUv0.cxx:629
 AliITSUv0.cxx:630
 AliITSUv0.cxx:631
 AliITSUv0.cxx:632
 AliITSUv0.cxx:633
 AliITSUv0.cxx:634
 AliITSUv0.cxx:635
 AliITSUv0.cxx:636
 AliITSUv0.cxx:637
 AliITSUv0.cxx:638
 AliITSUv0.cxx:639
 AliITSUv0.cxx:640
 AliITSUv0.cxx:641
 AliITSUv0.cxx:642
 AliITSUv0.cxx:643
 AliITSUv0.cxx:644
 AliITSUv0.cxx:645
 AliITSUv0.cxx:646
 AliITSUv0.cxx:647
 AliITSUv0.cxx:648
 AliITSUv0.cxx:649
 AliITSUv0.cxx:650
 AliITSUv0.cxx:651
 AliITSUv0.cxx:652
 AliITSUv0.cxx:653
 AliITSUv0.cxx:654
 AliITSUv0.cxx:655
 AliITSUv0.cxx:656
 AliITSUv0.cxx:657
 AliITSUv0.cxx:658
 AliITSUv0.cxx:659
 AliITSUv0.cxx:660
 AliITSUv0.cxx:661
 AliITSUv0.cxx:662
 AliITSUv0.cxx:663
 AliITSUv0.cxx:664
 AliITSUv0.cxx:665
 AliITSUv0.cxx:666
 AliITSUv0.cxx:667
 AliITSUv0.cxx:668
 AliITSUv0.cxx:669
 AliITSUv0.cxx:670
 AliITSUv0.cxx:671
 AliITSUv0.cxx:672
 AliITSUv0.cxx:673
 AliITSUv0.cxx:674
 AliITSUv0.cxx:675
 AliITSUv0.cxx:676
 AliITSUv0.cxx:677
 AliITSUv0.cxx:678
 AliITSUv0.cxx:679
 AliITSUv0.cxx:680
 AliITSUv0.cxx:681
 AliITSUv0.cxx:682
 AliITSUv0.cxx:683
 AliITSUv0.cxx:684
 AliITSUv0.cxx:685
 AliITSUv0.cxx:686
 AliITSUv0.cxx:687
 AliITSUv0.cxx:688
 AliITSUv0.cxx:689
 AliITSUv0.cxx:690
 AliITSUv0.cxx:691
 AliITSUv0.cxx:692
 AliITSUv0.cxx:693
 AliITSUv0.cxx:694
 AliITSUv0.cxx:695
 AliITSUv0.cxx:696
 AliITSUv0.cxx:697
 AliITSUv0.cxx:698
 AliITSUv0.cxx:699
 AliITSUv0.cxx:700
 AliITSUv0.cxx:701
 AliITSUv0.cxx:702
 AliITSUv0.cxx:703
 AliITSUv0.cxx:704
 AliITSUv0.cxx:705
 AliITSUv0.cxx:706
 AliITSUv0.cxx:707
 AliITSUv0.cxx:708
 AliITSUv0.cxx:709
 AliITSUv0.cxx:710
 AliITSUv0.cxx:711
 AliITSUv0.cxx:712
 AliITSUv0.cxx:713
 AliITSUv0.cxx:714
 AliITSUv0.cxx:715
 AliITSUv0.cxx:716
 AliITSUv0.cxx:717
 AliITSUv0.cxx:718
 AliITSUv0.cxx:719
 AliITSUv0.cxx:720
 AliITSUv0.cxx:721
 AliITSUv0.cxx:722
 AliITSUv0.cxx:723
 AliITSUv0.cxx:724
 AliITSUv0.cxx:725
 AliITSUv0.cxx:726
 AliITSUv0.cxx:727
 AliITSUv0.cxx:728
 AliITSUv0.cxx:729
 AliITSUv0.cxx:730
 AliITSUv0.cxx:731
 AliITSUv0.cxx:732
 AliITSUv0.cxx:733
 AliITSUv0.cxx:734
 AliITSUv0.cxx:735
 AliITSUv0.cxx:736
 AliITSUv0.cxx:737
 AliITSUv0.cxx:738
 AliITSUv0.cxx:739
 AliITSUv0.cxx:740
 AliITSUv0.cxx:741
 AliITSUv0.cxx:742
 AliITSUv0.cxx:743
 AliITSUv0.cxx:744
 AliITSUv0.cxx:745
 AliITSUv0.cxx:746
 AliITSUv0.cxx:747
 AliITSUv0.cxx:748
 AliITSUv0.cxx:749
 AliITSUv0.cxx:750
 AliITSUv0.cxx:751
 AliITSUv0.cxx:752
 AliITSUv0.cxx:753
 AliITSUv0.cxx:754
 AliITSUv0.cxx:755
 AliITSUv0.cxx:756
 AliITSUv0.cxx:757
 AliITSUv0.cxx:758
 AliITSUv0.cxx:759
 AliITSUv0.cxx:760
 AliITSUv0.cxx:761
 AliITSUv0.cxx:762
 AliITSUv0.cxx:763
 AliITSUv0.cxx:764
 AliITSUv0.cxx:765
 AliITSUv0.cxx:766
 AliITSUv0.cxx:767
 AliITSUv0.cxx:768
 AliITSUv0.cxx:769
 AliITSUv0.cxx:770
 AliITSUv0.cxx:771
 AliITSUv0.cxx:772
 AliITSUv0.cxx:773