[u/mrichter/AliRoot.git] / TRD / AliTRDgeometryHole.cxx

/**************************************************************************
 * 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.1  2000/02/28 19:01:42  cblume
Add new TRD classes

*/

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  TRD geometry with holes                                                  //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "AliTRDgeometryHole.h"

ClassImp(AliTRDgeometryHole)

//_____________________________________________________________________________
AliTRDgeometryHole::AliTRDgeometryHole():AliTRDgeometry()
{
  //
  // AliTRDgeometryHole default constructor
  //

  Init();

}

//_____________________________________________________________________________
AliTRDgeometryHole::~AliTRDgeometryHole()
{
  //
  // AliTRDgeometryHole destructor
  //

}

//_____________________________________________________________________________
void AliTRDgeometryHole::Init()
{
  //
  // Initializes the geometry parameter
  //

  Int_t iplan;

  // The length of the inner chambers
  for (iplan = 0; iplan < kNplan; iplan++) 
    fClengthI[iplan] = 110.0;

  // The length of the middle chambers
  fClengthM1[0] = 123.5;
  fClengthM1[1] = 131.0;
  fClengthM1[2] = 138.5;
  fClengthM1[3] = 146.0;
  fClengthM1[4] = 153.0;
  fClengthM1[5] = 160.5;

  fClengthM2[0] = 123.5 - 7.0;
  fClengthM2[1] = 131.0 - 7.0;
  fClengthM2[2] = 138.5 - 7.0;
  fClengthM2[3] = 146.0 - 7.0;
  fClengthM2[4] = 153.0 - 7.0;
  fClengthM2[5] = 160.4 - 7.0;

  // The length of the outer chambers
  fClengthO1[0] = 123.5;
  fClengthO1[1] = 131.0;
  fClengthO1[2] = 134.5;
  fClengthO1[3] = 142.0;
  fClengthO1[4] = 142.0;
  fClengthO1[5] = 134.5;

  fClengthO2[0] = 123.5;
  fClengthO2[1] = 131.0;
  fClengthO2[2] = 134.5;
  fClengthO2[3] = 142.0;
  fClengthO2[4] = 142.0;
  fClengthO2[5] = 134.5;

  fClengthO3[0] =  86.5;
  fClengthO3[1] = 101.5;
  fClengthO3[2] = 112.5;
  fClengthO3[3] = 127.5;
  fClengthO3[4] = 134.5;
  fClengthO3[5] = 134.5;

  // The maximum number of pads
  // and the position of pad 0,0,0 
  // 
  // chambers seen from the top:
  //     +----------------------------+
  //     |                            |
  //     |                            |     ^
  //     |                            | rphi|
  //     |                            |     |
  //     |0                           |     | 
  //     +----------------------------+     +------>
  //                                             z 
  // chambers seen from the side:           ^
  //     +----------------------------+ time|
  //     |                            |     |
  //     |0                           |     |
  //     +----------------------------+     +------>
  //                                             z
  //                                             

  // The pad row (z-direction)
  for (iplan = 0; iplan < kNplan; iplan++) {

    for (Int_t isect = 0; isect < kNsect; isect++) {
      Float_t clengthI = fClengthI[iplan];
      Float_t clengthM = fClengthM1[iplan];
      Float_t clengthO = fClengthO1[iplan];
      switch (isect) {
      case 12:
      case 13:
      case 14:
      case 15:
      case 16:
        clengthM = fClengthM2[iplan];
        clengthO = fClengthO2[iplan];
        break;
      case 4:
      case 5:
      case 6:
        clengthO = fClengthO3[iplan];
        break;
      };
      fRowMax[iplan][0][isect] = 1 + TMath::Nint((clengthO - 2. * kCcthick) 
                                                           / fRowPadSize - 0.5);
      fRowMax[iplan][1][isect] = 1 + TMath::Nint((clengthM - 2. * kCcthick) 
                                                           / fRowPadSize - 0.5);
      fRowMax[iplan][2][isect] = 1 + TMath::Nint((clengthI - 2. * kCcthick) 
                                                           / fRowPadSize - 0.5);
      fRowMax[iplan][3][isect] = 1 + TMath::Nint((clengthM - 2. * kCcthick) 
                                                           / fRowPadSize - 0.5);
      fRowMax[iplan][4][isect] = 1 + TMath::Nint((clengthO - 2. * kCcthick) 
                                                           / fRowPadSize - 0.5);
      fRow0[iplan][0][isect]   = -clengthI/2. - clengthM - clengthO + kCcthick; 
      fRow0[iplan][1][isect]   = -clengthI/2. - clengthM            + kCcthick;
      fRow0[iplan][2][isect]   = -clengthI/2.                       + kCcthick;
      fRow0[iplan][3][isect]   =  clengthI/2.                       + kCcthick; 
      fRow0[iplan][4][isect]   =  clengthI/2. + clengthM            + kCcthick; 
    }

  }

}

//_____________________________________________________________________________
void AliTRDgeometryHole::CreateGeometry(Int_t *idtmed)
{
  //
  // Create the TRD geometry with hole
  //

  Int_t iplan;

  const Int_t kNparTrd = 4;
  const Int_t kNparCha = 3;

  Float_t parTrd[kNparTrd];
  Float_t parCha[kNparCha];

  Float_t xpos, ypos, zpos;

  AliTRDgeometry::CreateGeometry(idtmed);

  // The TRD mother volume for one sector (Air) (dimensions identical to BTR1)
  parTrd[0] = kSwidth1/2.;
  parTrd[1] = kSwidth2/2.;
  parTrd[2] = kSlenTR1/2.;
  parTrd[3] = kSheight/2.;
  gMC->Gsvolu("TRD1","TRD1",idtmed[1302-1],parTrd,kNparTrd);
  
  // The TRD mother volume for one sector (Air) (dimensions identical to BTR2)
  parTrd[0] = kSwidth1/2.;
  parTrd[1] = kSwidth2/2.;
  parTrd[2] = kSlenTR2/2.;
  parTrd[3] = kSheight/2.;
  gMC->Gsvolu("TRD2","TRD1",idtmed[1302-1],parTrd,kNparTrd);

  // The TRD mother volume for one sector (Air) (dimensions identical to BTR3)
  parTrd[0] = kSwidth1/2.;
  parTrd[1] = kSwidth2/2.;
  parTrd[2] = kSlenTR3/2.;
  parTrd[3] = kSheight/2.;
  gMC->Gsvolu("TRD3","TRD1",idtmed[1302-1],parTrd,kNparTrd);

  // Position the chambers in the TRD mother volume
  for (iplan = 1; iplan <= kNplan; iplan++) {

    // The inner chambers ---------------------------------------------------------------

    // the aluminum frame
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthI[iplan-1]/2.;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = 0.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAFI",iplan       ,"TRD1",xpos,ypos,zpos,0,"MANY",parCha,kNparCha);

    // the inner part of the aluminum frame
    parCha[0] = fCwidth[iplan-1]/2.   - kCathick;
    parCha[1] = fClengthI[iplan-1]/2. - kCathick;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = 0.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAII",iplan       ,"TRD1",xpos,ypos,zpos,0,"ONLY",parCha,kNparCha);

    // the carbon frame
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthI[iplan-1]/2.;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = 0.;
    zpos       = kCcframe/2.            - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCFI",iplan       ,"TRD1",xpos,ypos,zpos,0,"MANY",parCha,kNparCha);

    // the inner part of the carbon frame
    parCha[0] = fCwidth[iplan-1]/2.   - kCcthick;
    parCha[1] = fClengthI[iplan-1]/2. - kCcthick;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = 0.;
    zpos       = kCcframe/2.            - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCII",iplan       ,"TRD1",xpos,ypos,zpos,0,"ONLY",parCha,kNparCha);

    // The middle chambers --------------------------------------------------------------

    // the aluminum frame
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthM1[iplan-1]/2.;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthI[iplan-1]/2.  + fClengthM1[iplan-1]/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAFM",iplan         ,"TRD1",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
    gMC->Gsposp("UAFM",iplan+  kNplan,"TRD1",xpos,-ypos,zpos,0,"MANY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthM2[iplan-1]/2.;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthM2[iplan-1]/2. - kSlenTR2/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAFM",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);

    // the inner part of the aluminum frame
    parCha[0] = fCwidth[iplan-1]/2.    - kCathick;
    parCha[1] = fClengthM1[iplan-1]/2. - kCathick;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthI[iplan-1]/2.  + fClengthM1[iplan-1]/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAIM",iplan         ,"TRD1",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
    gMC->Gsposp("UAIM",iplan+  kNplan,"TRD1",xpos,-ypos,zpos,0,"ONLY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.    - kCathick;
    parCha[1] = fClengthM2[iplan-1]/2. - kCathick;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthM2[iplan-1]/2. - kSlenTR2/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAIM",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);

    // the carbon frame
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthM1[iplan-1]/2.;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]/2.;
    zpos       = kCcframe/2.           - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCFM",iplan         ,"TRD1",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
    gMC->Gsposp("UCFM",iplan+  kNplan,"TRD1",xpos,-ypos,zpos,0,"MANY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthM2[iplan-1]/2.;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthM2[iplan-1]/2. - kSlenTR2/2.;
    zpos       = kCcframe/2.            - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCFM",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);

    // the inner part of the carbon frame
    parCha[0] = fCwidth[iplan-1]/2.    - kCcthick;
    parCha[1] = fClengthM1[iplan-1]/2. - kCcthick;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]/2.;
    zpos       = kCcframe/2.           - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCIM",iplan         ,"TRD1",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
    gMC->Gsposp("UCIM",iplan+  kNplan,"TRD1",xpos,-ypos,zpos,0,"ONLY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.    - kCcthick;
    parCha[1] = fClengthM2[iplan-1]/2. - kCcthick;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthM2[iplan-1]/2. - kSlenTR2/2.;
    zpos       = kCcframe/2.            - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCIM",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);

    // The outer chambers ---------------------------------------------------------------

    // the aluminum frame
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthO1[iplan-1]/2.;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]    + fClengthO1[iplan-1]/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAFO",iplan         ,"TRD1",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
    gMC->Gsposp("UAFO",iplan+  kNplan,"TRD1",xpos,-ypos,zpos,0,"MANY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthO2[iplan-1]/2.;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthM2[iplan-1]    + fClengthO2[iplan-1]/2. - kSlenTR2/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAFO",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthO3[iplan-1]/2.;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthO3[iplan-1]/2. - kSlenTR3/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAFO",iplan+4*kNplan,"TRD3",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);

    // the inner part of the aluminum frame
    parCha[0] = fCwidth[iplan-1]/2.    - kCathick;
    parCha[1] = fClengthO1[iplan-1]/2. - kCathick;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]    + fClengthO1[iplan-1]/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAIO",iplan         ,"TRD1",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
    gMC->Gsposp("UAIO",iplan+  kNplan,"TRD1",xpos,-ypos,zpos,0,"ONLY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.    - kCathick;
    parCha[1] = fClengthO2[iplan-1]/2. - kCathick;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthM2[iplan-1]    + fClengthO2[iplan-1]/2. - kSlenTR2/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAIO",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.    - kCathick;
    parCha[1] = fClengthO3[iplan-1]/2. - kCathick;
    parCha[2] = kCaframe/2.;
    xpos       = 0.;
    ypos       = fClengthO3[iplan-1]/2. - kSlenTR3/2.;
    zpos       = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UAIO",iplan+4*kNplan,"TRD3",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);

    // the carbon frame
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthO1[iplan-1]/2.;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]    + fClengthO1[iplan-1]/2.;
    zpos       = kCcframe/2.           - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCFO",iplan,         "TRD1",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
    gMC->Gsposp("UCFO",iplan+  kNplan,"TRD1",xpos,-ypos,zpos,0,"MANY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthO2[iplan-1]/2.;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthM2[iplan-1]    + fClengthO2[iplan-1]/2. - kSlenTR2/2.;
    zpos       = kCcframe/2.            - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCFO",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.;
    parCha[1] = fClengthO3[iplan-1]/2.;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthO3[iplan-1]/2. - kSlenTR3/2.;
    zpos       = kCcframe/2.            - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCFO",iplan+4*kNplan,"TRD3",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);

    // the inner part of the carbon frame
    parCha[0] = fCwidth[iplan-1]/2.    - kCcthick;
    parCha[1] = fClengthO1[iplan-1]/2. - kCcthick;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]    + fClengthO1[iplan-1]/2.;
    zpos       = kCcframe/2.           - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCIO",iplan         ,"TRD1",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
    gMC->Gsposp("UCIO",iplan+  kNplan,"TRD1",xpos,-ypos,zpos,0,"ONLY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.    - kCcthick;
    parCha[1] = fClengthO2[iplan-1]/2. - kCcthick;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthM2[iplan-1]    + fClengthO2[iplan-1]/2. - kSlenTR2/2.;
    zpos       = kCcframe/2.            - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCIO",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
    parCha[0] = fCwidth[iplan-1]/2.    - kCcthick;
    parCha[1] = fClengthO3[iplan-1]/2. - kCcthick;
    parCha[2] = kCcframe/2.;
    xpos       = 0.;
    ypos       = fClengthO3[iplan-1]/2. - kSlenTR3/2.;
    zpos       = kCcframe/2.            - kSheight/2. + (iplan-1) * (kCheight + kCspace);
    gMC->Gsposp("UCIO",iplan+4*kNplan,"TRD3",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);

  }

  xpos     = 0.;
  ypos     = 0.;
  zpos     = 0.;
  gMC->Gspos("TRD1",1,"BTR1",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("TRD2",1,"BTR2",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("TRD3",1,"BTR3",xpos,ypos,zpos,0,"ONLY");

}
Commit	Line	Data
f7336fa3	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/*
	17	$Log$
8230f242	18	Revision 1.1 2000/02/28 19:01:42 cblume
	19	Add new TRD classes
	20
f7336fa3	21	*/
	22
	23	///////////////////////////////////////////////////////////////////////////////
	24	// //
	25	// TRD geometry with holes //
	26	// //
	27	///////////////////////////////////////////////////////////////////////////////
	28
	29	#include "AliTRDgeometryHole.h"
	30
	31	ClassImp(AliTRDgeometryHole)
	32
	33	//_____________________________________________________________________________
	34	AliTRDgeometryHole::AliTRDgeometryHole():AliTRDgeometry()
	35	{
	36	//
	37	// AliTRDgeometryHole default constructor
	38	//
	39
	40	Init();
	41
	42	}
	43
	44	//_____________________________________________________________________________
	45	AliTRDgeometryHole::~AliTRDgeometryHole()
	46	{
8230f242	47	//
	48	// AliTRDgeometryHole destructor
	49	//
f7336fa3	50
	51	}
	52
	53	//_____________________________________________________________________________
	54	void AliTRDgeometryHole::Init()
	55	{
	56	//
	57	// Initializes the geometry parameter
	58	//
	59
	60	Int_t iplan;
	61
	62	// The length of the inner chambers
	63	for (iplan = 0; iplan < kNplan; iplan++)
	64	fClengthI[iplan] = 110.0;
	65
	66	// The length of the middle chambers
	67	fClengthM1[0] = 123.5;
	68	fClengthM1[1] = 131.0;
	69	fClengthM1[2] = 138.5;
	70	fClengthM1[3] = 146.0;
	71	fClengthM1[4] = 153.0;
	72	fClengthM1[5] = 160.5;
	73
	74	fClengthM2[0] = 123.5 - 7.0;
	75	fClengthM2[1] = 131.0 - 7.0;
	76	fClengthM2[2] = 138.5 - 7.0;
	77	fClengthM2[3] = 146.0 - 7.0;
	78	fClengthM2[4] = 153.0 - 7.0;
	79	fClengthM2[5] = 160.4 - 7.0;
	80
	81	// The length of the outer chambers
	82	fClengthO1[0] = 123.5;
	83	fClengthO1[1] = 131.0;
	84	fClengthO1[2] = 134.5;
	85	fClengthO1[3] = 142.0;
	86	fClengthO1[4] = 142.0;
	87	fClengthO1[5] = 134.5;
	88
	89	fClengthO2[0] = 123.5;
	90	fClengthO2[1] = 131.0;
	91	fClengthO2[2] = 134.5;
	92	fClengthO2[3] = 142.0;
	93	fClengthO2[4] = 142.0;
	94	fClengthO2[5] = 134.5;
	95
	96	fClengthO3[0] = 86.5;
	97	fClengthO3[1] = 101.5;
	98	fClengthO3[2] = 112.5;
	99	fClengthO3[3] = 127.5;
	100	fClengthO3[4] = 134.5;
	101	fClengthO3[5] = 134.5;
	102
	103	// The maximum number of pads
	104	// and the position of pad 0,0,0
	105	//
	106	// chambers seen from the top:
	107	// +----------------------------+
	108	// \| \|
	109	// \| \| ^
	110	// \| \| rphi\|
	111	// \| \| \|
	112	// \|0 \| \|
	113	// +----------------------------+ +------>
114	// z
115	// chambers seen from the side: ^
116	// +----------------------------+ time\|
117	// \| \| \|
118	// \|0 \| \|
119	// +----------------------------+ +------>
120	// z
121	//
122
123	// The pad row (z-direction)
124	for (iplan = 0; iplan < kNplan; iplan++) {
125
126	for (Int_t isect = 0; isect < kNsect; isect++) {
127	Float_t clengthI = fClengthI[iplan];
128	Float_t clengthM = fClengthM1[iplan];
129	Float_t clengthO = fClengthO1[iplan];
130	switch (isect) {
131	case 12:
132	case 13:
133	case 14:
134	case 15:
135	case 16:
136	clengthM = fClengthM2[iplan];
137	clengthO = fClengthO2[iplan];
138	break;
139	case 4:
140	case 5:
141	case 6:
142	clengthO = fClengthO3[iplan];
143	break;
144	};
145	fRowMax[iplan][0][isect] = 1 + TMath::Nint((clengthO - 2. * kCcthick)
146	/ fRowPadSize - 0.5);
147	fRowMax[iplan][1][isect] = 1 + TMath::Nint((clengthM - 2. * kCcthick)
148	/ fRowPadSize - 0.5);
149	fRowMax[iplan][2][isect] = 1 + TMath::Nint((clengthI - 2. * kCcthick)
150	/ fRowPadSize - 0.5);
151	fRowMax[iplan][3][isect] = 1 + TMath::Nint((clengthM - 2. * kCcthick)
152	/ fRowPadSize - 0.5);
153	fRowMax[iplan][4][isect] = 1 + TMath::Nint((clengthO - 2. * kCcthick)
154	/ fRowPadSize - 0.5);
155	fRow0[iplan][0][isect] = -clengthI/2. - clengthM - clengthO + kCcthick;
156	fRow0[iplan][1][isect] = -clengthI/2. - clengthM + kCcthick;
157	fRow0[iplan][2][isect] = -clengthI/2. + kCcthick;
158	fRow0[iplan][3][isect] = clengthI/2. + kCcthick;
159	fRow0[iplan][4][isect] = clengthI/2. + clengthM + kCcthick;
160	}
161
162	}
163
164	}
165
166	//_____________________________________________________________________________
167	void AliTRDgeometryHole::CreateGeometry(Int_t *idtmed)
168	{
169	//
170	// Create the TRD geometry with hole
171	//
172
173	Int_t iplan;
174
8230f242	175	const Int_t kNparTrd = 4;
8230f242	176	const Int_t kNparCha = 3;
f7336fa3	177
8230f242	178	Float_t parTrd[kNparTrd];
8230f242	179	Float_t parCha[kNparCha];
f7336fa3	180
	181	Float_t xpos, ypos, zpos;
	182
	183	AliTRDgeometry::CreateGeometry(idtmed);
	184
	185	// The TRD mother volume for one sector (Air) (dimensions identical to BTR1)
8230f242	186	parTrd[0] = kSwidth1/2.;
	187	parTrd[1] = kSwidth2/2.;
	188	parTrd[2] = kSlenTR1/2.;
	189	parTrd[3] = kSheight/2.;
	190	gMC->Gsvolu("TRD1","TRD1",idtmed[1302-1],parTrd,kNparTrd);
f7336fa3	191
f7336fa3	192	// The TRD mother volume for one sector (Air) (dimensions identical to BTR2)
8230f242	193	parTrd[0] = kSwidth1/2.;
	194	parTrd[1] = kSwidth2/2.;
	195	parTrd[2] = kSlenTR2/2.;
	196	parTrd[3] = kSheight/2.;
	197	gMC->Gsvolu("TRD2","TRD1",idtmed[1302-1],parTrd,kNparTrd);
f7336fa3	198
f7336fa3	199	// The TRD mother volume for one sector (Air) (dimensions identical to BTR3)
8230f242	200	parTrd[0] = kSwidth1/2.;
	201	parTrd[1] = kSwidth2/2.;
	202	parTrd[2] = kSlenTR3/2.;
	203	parTrd[3] = kSheight/2.;
	204	gMC->Gsvolu("TRD3","TRD1",idtmed[1302-1],parTrd,kNparTrd);
f7336fa3	205
	206	// Position the chambers in the TRD mother volume
	207	for (iplan = 1; iplan <= kNplan; iplan++) {
	208
	209	// The inner chambers ---------------------------------------------------------------
	210
	211	// the aluminum frame
8230f242	212	parCha[0] = fCwidth[iplan-1]/2.;
	213	parCha[1] = fClengthI[iplan-1]/2.;
	214	parCha[2] = kCaframe/2.;
f7336fa3	215	xpos = 0.;
	216	ypos = 0.;
	217	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	218	gMC->Gsposp("UAFI",iplan ,"TRD1",xpos,ypos,zpos,0,"MANY",parCha,kNparCha);
f7336fa3	219
f7336fa3	220	// the inner part of the aluminum frame
8230f242	221	parCha[0] = fCwidth[iplan-1]/2. - kCathick;
	222	parCha[1] = fClengthI[iplan-1]/2. - kCathick;
	223	parCha[2] = kCaframe/2.;
f7336fa3	224	xpos = 0.;
	225	ypos = 0.;
	226	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	227	gMC->Gsposp("UAII",iplan ,"TRD1",xpos,ypos,zpos,0,"ONLY",parCha,kNparCha);
f7336fa3	228
f7336fa3	229	// the carbon frame
8230f242	230	parCha[0] = fCwidth[iplan-1]/2.;
	231	parCha[1] = fClengthI[iplan-1]/2.;
	232	parCha[2] = kCcframe/2.;
f7336fa3	233	xpos = 0.;
	234	ypos = 0.;
	235	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	236	gMC->Gsposp("UCFI",iplan ,"TRD1",xpos,ypos,zpos,0,"MANY",parCha,kNparCha);
f7336fa3	237
f7336fa3	238	// the inner part of the carbon frame
8230f242	239	parCha[0] = fCwidth[iplan-1]/2. - kCcthick;
	240	parCha[1] = fClengthI[iplan-1]/2. - kCcthick;
	241	parCha[2] = kCcframe/2.;
f7336fa3	242	xpos = 0.;
	243	ypos = 0.;
	244	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	245	gMC->Gsposp("UCII",iplan ,"TRD1",xpos,ypos,zpos,0,"ONLY",parCha,kNparCha);
f7336fa3	246
	247	// The middle chambers --------------------------------------------------------------
	248
	249	// the aluminum frame
8230f242	250	parCha[0] = fCwidth[iplan-1]/2.;
	251	parCha[1] = fClengthM1[iplan-1]/2.;
	252	parCha[2] = kCaframe/2.;
f7336fa3	253	xpos = 0.;
	254	ypos = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]/2.;
	255	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	256	gMC->Gsposp("UAFM",iplan ,"TRD1",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
	257	gMC->Gsposp("UAFM",iplan+ kNplan,"TRD1",xpos,-ypos,zpos,0,"MANY",parCha,kNparCha);
	258	parCha[0] = fCwidth[iplan-1]/2.;
	259	parCha[1] = fClengthM2[iplan-1]/2.;
	260	parCha[2] = kCaframe/2.;
f7336fa3	261	xpos = 0.;
	262	ypos = fClengthM2[iplan-1]/2. - kSlenTR2/2.;
	263	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	264	gMC->Gsposp("UAFM",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
f7336fa3	265
f7336fa3	266	// the inner part of the aluminum frame
8230f242	267	parCha[0] = fCwidth[iplan-1]/2. - kCathick;
	268	parCha[1] = fClengthM1[iplan-1]/2. - kCathick;
	269	parCha[2] = kCaframe/2.;
f7336fa3	270	xpos = 0.;
	271	ypos = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]/2.;
	272	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	273	gMC->Gsposp("UAIM",iplan ,"TRD1",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
	274	gMC->Gsposp("UAIM",iplan+ kNplan,"TRD1",xpos,-ypos,zpos,0,"ONLY",parCha,kNparCha);
	275	parCha[0] = fCwidth[iplan-1]/2. - kCathick;
	276	parCha[1] = fClengthM2[iplan-1]/2. - kCathick;
	277	parCha[2] = kCaframe/2.;
f7336fa3	278	xpos = 0.;
	279	ypos = fClengthM2[iplan-1]/2. - kSlenTR2/2.;
	280	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	281	gMC->Gsposp("UAIM",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
f7336fa3	282
f7336fa3	283	// the carbon frame
8230f242	284	parCha[0] = fCwidth[iplan-1]/2.;
	285	parCha[1] = fClengthM1[iplan-1]/2.;
	286	parCha[2] = kCcframe/2.;
f7336fa3	287	xpos = 0.;
	288	ypos = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]/2.;
	289	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	290	gMC->Gsposp("UCFM",iplan ,"TRD1",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
	291	gMC->Gsposp("UCFM",iplan+ kNplan,"TRD1",xpos,-ypos,zpos,0,"MANY",parCha,kNparCha);
	292	parCha[0] = fCwidth[iplan-1]/2.;
	293	parCha[1] = fClengthM2[iplan-1]/2.;
	294	parCha[2] = kCcframe/2.;
f7336fa3	295	xpos = 0.;
	296	ypos = fClengthM2[iplan-1]/2. - kSlenTR2/2.;
	297	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	298	gMC->Gsposp("UCFM",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
f7336fa3	299
f7336fa3	300	// the inner part of the carbon frame
8230f242	301	parCha[0] = fCwidth[iplan-1]/2. - kCcthick;
	302	parCha[1] = fClengthM1[iplan-1]/2. - kCcthick;
	303	parCha[2] = kCcframe/2.;
f7336fa3	304	xpos = 0.;
	305	ypos = fClengthI[iplan-1]/2. + fClengthM1[iplan-1]/2.;
	306	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	307	gMC->Gsposp("UCIM",iplan ,"TRD1",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
	308	gMC->Gsposp("UCIM",iplan+ kNplan,"TRD1",xpos,-ypos,zpos,0,"ONLY",parCha,kNparCha);
	309	parCha[0] = fCwidth[iplan-1]/2. - kCcthick;
	310	parCha[1] = fClengthM2[iplan-1]/2. - kCcthick;
	311	parCha[2] = kCcframe/2.;
f7336fa3	312	xpos = 0.;
	313	ypos = fClengthM2[iplan-1]/2. - kSlenTR2/2.;
	314	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	315	gMC->Gsposp("UCIM",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
f7336fa3	316
	317	// The outer chambers ---------------------------------------------------------------
	318
	319	// the aluminum frame
8230f242	320	parCha[0] = fCwidth[iplan-1]/2.;
	321	parCha[1] = fClengthO1[iplan-1]/2.;
	322	parCha[2] = kCaframe/2.;
f7336fa3	323	xpos = 0.;
	324	ypos = fClengthI[iplan-1]/2. + fClengthM1[iplan-1] + fClengthO1[iplan-1]/2.;
	325	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	326	gMC->Gsposp("UAFO",iplan ,"TRD1",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
	327	gMC->Gsposp("UAFO",iplan+ kNplan,"TRD1",xpos,-ypos,zpos,0,"MANY",parCha,kNparCha);
	328	parCha[0] = fCwidth[iplan-1]/2.;
	329	parCha[1] = fClengthO2[iplan-1]/2.;
	330	parCha[2] = kCaframe/2.;
f7336fa3	331	xpos = 0.;
	332	ypos = fClengthM2[iplan-1] + fClengthO2[iplan-1]/2. - kSlenTR2/2.;
	333	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	334	gMC->Gsposp("UAFO",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
	335	parCha[0] = fCwidth[iplan-1]/2.;
	336	parCha[1] = fClengthO3[iplan-1]/2.;
	337	parCha[2] = kCaframe/2.;
f7336fa3	338	xpos = 0.;
	339	ypos = fClengthO3[iplan-1]/2. - kSlenTR3/2.;
	340	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	341	gMC->Gsposp("UAFO",iplan+4*kNplan,"TRD3",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
f7336fa3	342
f7336fa3	343	// the inner part of the aluminum frame
8230f242	344	parCha[0] = fCwidth[iplan-1]/2. - kCathick;
	345	parCha[1] = fClengthO1[iplan-1]/2. - kCathick;
	346	parCha[2] = kCaframe/2.;
f7336fa3	347	xpos = 0.;
	348	ypos = fClengthI[iplan-1]/2. + fClengthM1[iplan-1] + fClengthO1[iplan-1]/2.;
	349	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	350	gMC->Gsposp("UAIO",iplan ,"TRD1",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
	351	gMC->Gsposp("UAIO",iplan+ kNplan,"TRD1",xpos,-ypos,zpos,0,"ONLY",parCha,kNparCha);
	352	parCha[0] = fCwidth[iplan-1]/2. - kCathick;
	353	parCha[1] = fClengthO2[iplan-1]/2. - kCathick;
	354	parCha[2] = kCaframe/2.;
f7336fa3	355	xpos = 0.;
	356	ypos = fClengthM2[iplan-1] + fClengthO2[iplan-1]/2. - kSlenTR2/2.;
	357	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	358	gMC->Gsposp("UAIO",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
	359	parCha[0] = fCwidth[iplan-1]/2. - kCathick;
	360	parCha[1] = fClengthO3[iplan-1]/2. - kCathick;
	361	parCha[2] = kCaframe/2.;
f7336fa3	362	xpos = 0.;
	363	ypos = fClengthO3[iplan-1]/2. - kSlenTR3/2.;
	364	zpos = kCheight - kCaframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	365	gMC->Gsposp("UAIO",iplan+4*kNplan,"TRD3",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
f7336fa3	366
f7336fa3	367	// the carbon frame
8230f242	368	parCha[0] = fCwidth[iplan-1]/2.;
	369	parCha[1] = fClengthO1[iplan-1]/2.;
	370	parCha[2] = kCcframe/2.;
f7336fa3	371	xpos = 0.;
	372	ypos = fClengthI[iplan-1]/2. + fClengthM1[iplan-1] + fClengthO1[iplan-1]/2.;
	373	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	374	gMC->Gsposp("UCFO",iplan, "TRD1",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
	375	gMC->Gsposp("UCFO",iplan+ kNplan,"TRD1",xpos,-ypos,zpos,0,"MANY",parCha,kNparCha);
	376	parCha[0] = fCwidth[iplan-1]/2.;
	377	parCha[1] = fClengthO2[iplan-1]/2.;
	378	parCha[2] = kCcframe/2.;
f7336fa3	379	xpos = 0.;
	380	ypos = fClengthM2[iplan-1] + fClengthO2[iplan-1]/2. - kSlenTR2/2.;
	381	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	382	gMC->Gsposp("UCFO",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
	383	parCha[0] = fCwidth[iplan-1]/2.;
	384	parCha[1] = fClengthO3[iplan-1]/2.;
	385	parCha[2] = kCcframe/2.;
f7336fa3	386	xpos = 0.;
	387	ypos = fClengthO3[iplan-1]/2. - kSlenTR3/2.;
	388	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	389	gMC->Gsposp("UCFO",iplan+4*kNplan,"TRD3",xpos, ypos,zpos,0,"MANY",parCha,kNparCha);
f7336fa3	390
f7336fa3	391	// the inner part of the carbon frame
8230f242	392	parCha[0] = fCwidth[iplan-1]/2. - kCcthick;
	393	parCha[1] = fClengthO1[iplan-1]/2. - kCcthick;
	394	parCha[2] = kCcframe/2.;
f7336fa3	395	xpos = 0.;
	396	ypos = fClengthI[iplan-1]/2. + fClengthM1[iplan-1] + fClengthO1[iplan-1]/2.;
	397	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	398	gMC->Gsposp("UCIO",iplan ,"TRD1",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
	399	gMC->Gsposp("UCIO",iplan+ kNplan,"TRD1",xpos,-ypos,zpos,0,"ONLY",parCha,kNparCha);
	400	parCha[0] = fCwidth[iplan-1]/2. - kCcthick;
	401	parCha[1] = fClengthO2[iplan-1]/2. - kCcthick;
	402	parCha[2] = kCcframe/2.;
f7336fa3	403	xpos = 0.;
	404	ypos = fClengthM2[iplan-1] + fClengthO2[iplan-1]/2. - kSlenTR2/2.;
	405	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	406	gMC->Gsposp("UCIO",iplan+2*kNplan,"TRD2",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
	407	parCha[0] = fCwidth[iplan-1]/2. - kCcthick;
	408	parCha[1] = fClengthO3[iplan-1]/2. - kCcthick;
	409	parCha[2] = kCcframe/2.;
f7336fa3	410	xpos = 0.;
	411	ypos = fClengthO3[iplan-1]/2. - kSlenTR3/2.;
	412	zpos = kCcframe/2. - kSheight/2. + (iplan-1) * (kCheight + kCspace);
8230f242	413	gMC->Gsposp("UCIO",iplan+4*kNplan,"TRD3",xpos, ypos,zpos,0,"ONLY",parCha,kNparCha);
f7336fa3	414
	415	}
	416
	417	xpos = 0.;
	418	ypos = 0.;
	419	zpos = 0.;
	420	gMC->Gspos("TRD1",1,"BTR1",xpos,ypos,zpos,0,"ONLY");
	421	gMC->Gspos("TRD2",1,"BTR2",xpos,ypos,zpos,0,"ONLY");
	422	gMC->Gspos("TRD3",1,"BTR3",xpos,ypos,zpos,0,"ONLY");
	423
	424	}