[u/mrichter/AliRoot.git] / TRD / AliTRDgeometryDetail.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.6  2002/10/14 14:57:43  hristov
Merging the VirtualMC branch to the main development branch (HEAD)

Revision 1.3.6.2  2002/07/24 10:09:30  alibrary
Updating VirtualMC

Revision 1.5  2002/06/12 09:54:35  cblume
Update of tracking code provided by Sergei

Revision 1.4  2002/03/28 14:59:07  cblume
Coding conventions

Revision 1.3  2002/02/11 14:21:16  cblume
Update of the geometry. Get rid of MANY

Revision 1.2  2001/11/08 13:13:08  cblume
Change to MANY for UCFI/M/O and UAFI/M/O

Revision 1.1  2001/11/06 17:19:41  cblume
Add detailed geometry and simple simulator

*/

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  Detailed TRD geometry for the spaceframe without holes                   //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "TVirtualMC.h"

#include "AliTRDgeometryDetail.h"
#include "AliTRDparameter.h"

ClassImp(AliTRDgeometryDetail)

//_____________________________________________________________________________
AliTRDgeometryDetail::AliTRDgeometryDetail():AliTRDgeometryFull()
{
  //
  // AliTRDgeometryDetail default constructor
  //

  Init();

}

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

}

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

  AliTRDgeometryFull::Init();

}

//_____________________________________________________________________________
void AliTRDgeometryDetail::CreateGeometry(Int_t *idtmed)
{
  //
  // Create the detailed TRD geometry without hole
  // including the MCMs and the cooling pipes
  //
  //
  // Names of the TRD volumina (xx = detector number):
  //
  //      Lower part of the readout chambers (gas volume + radiator)
  //
  //        UAxx    Aluminum frames             (Al)
  //        UBxx    G10 frames                  (C)
  //        UCxx    Inner volumes               (Air)
  //
  //      Upper part of the readout chambers (readout plane + fee)
  //
  //        UDxx    G10 frames                  (C)
  //        UExx    Inner volumes of the G10    (Air)
  //        UFxx    Aluminum frames             (Al)
  //        UGxx    Inner volumes of the Al     (Air)
  //
  //      Inner material layers
  //
  //        UHxx    Radiator                    (Rohacell)
  //        UIxx    Entrance window             (Mylar)
  //        UJxx    Drift volume                (Xe/CO2)
  //        UKxx    Amplification volume        (Xe/CO2)
  //        ULxx    Pad plane                   (Cu)
  //        UMxx    Support structure           (Rohacell)
  //        UNxx    FEE + signal lines          (Cu)
  //

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

  Float_t xpos, ypos, zpos;

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

  Char_t  cTagV[5];
  Char_t  cTagM[5];

  Int_t   idrotm;

  // Rotation matrix
  gMC->Matrix(idrotm,  0.0,  0.0, 90.0, 90.0, 90.0,  0.0);

  // The TRD mother volume for one sector (Air), full length in z-direction
  parTrd[0] = fgkSwidth1/2.;
  parTrd[1] = fgkSwidth2/2.;
  parTrd[2] = fgkSlenTR1/2.;
  parTrd[3] = fgkSheight/2.;
  gMC->Gsvolu("UTR1","TRD1",idtmed[1302-1],parTrd,kNparTrd);

  // Create the readout volumina
  CreateReadout(idtmed);

  // Create the volumina for the cooling
  CreateCooling(idtmed);

  for (Int_t icham = 0; icham < kNcham; icham++) {
    for (Int_t iplan = 0; iplan < kNplan; iplan++) {  

      Int_t iDet = GetDetectorSec(iplan,icham);

      // The lower part of the readout chambers (gas volume + radiator) 
      // The aluminum frames 
      sprintf(cTagV,"UA%02d",iDet);
      parCha[0] = fCwidth[iplan]/2.;
      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
      parCha[2] = fgkCraH/2. + fgkCdrH/2.;
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
      // The G10 frames 
      sprintf(cTagV,"UB%02d",iDet);
      parCha[0] = fCwidth[iplan]/2. - fgkCalT; 
      parCha[1] = -1.;
      parCha[2] = -1.;
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
      // The inner part (air)
      sprintf(cTagV,"UC%02d",iDet);
      parCha[0] = fCwidth[iplan]/2. - fgkCalT - fgkCclsT; 
      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCclfT;
      parCha[2] = -1.;
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);

      // The upper part of the readout chambers (readout plane + fee)
      // The G10 frames
      sprintf(cTagV,"UD%02d",iDet);
      parCha[0] = fCwidth[iplan]/2. + fgkCroW;
      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
      parCha[2] = fgkCamH/2.;
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
      // The inner part of the G10 frame (air)
      sprintf(cTagV,"UE%02d",iDet);
      parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCcuT; 
      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCcuT;
      parCha[2] = -1.;
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
      // The aluminum frames
      sprintf(cTagV,"UF%02d",iDet);
      parCha[0] = fCwidth[iplan]/2. + fgkCroW;
      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
      parCha[2] = fgkCroH/2.;
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
      // The inner part of the aluminum frames
      sprintf(cTagV,"UG%02d",iDet);
      parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCauT; 
      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCauT;
      parCha[2] = -1.;
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);

      // The material layers inside the chambers
      parCha[0] = -1.;
      parCha[1] = -1.;
      // Rohacell layer (radiator)
      parCha[2] = fgkRaThick/2;
      sprintf(cTagV,"UH%02d",iDet);
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha);
      // Mylar layer (entrance window + HV cathode) 
      parCha[2] = fgkMyThick/2;
      sprintf(cTagV,"UI%02d",iDet);
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1308-1],parCha,kNparCha);
      // Xe/Isobutane layer (drift volume) 
      parCha[2] = fgkDrThick/2.;
      sprintf(cTagV,"UJ%02d",iDet);
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
      // Xe/Isobutane layer (amplification volume)
      parCha[2] = fgkAmThick/2.;
      sprintf(cTagV,"UK%02d",iDet);
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);  
      // Cu layer (pad plane)
      parCha[2] = fgkCuThick/2;
      sprintf(cTagV,"UL%02d",iDet);
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
      // G10 layer (support structure / honeycomb)
      parCha[2] = fgkSuThick/2;
      sprintf(cTagV,"UM%02d",iDet);
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha);
      // Cu layer (FEE + signal lines)
      parCha[2] = fgkFeThick/2;
      sprintf(cTagV,"UN%02d",iDet);
      gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);

      // Position the layers in the chambers
      xpos = 0;
      ypos = 0;
      // Lower part
      // Rohacell layer (radiator)
      zpos = fgkRaZpos;
      sprintf(cTagV,"UH%02d",iDet);
      sprintf(cTagM,"UC%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
      // Mylar layer (entrance window + HV cathode)   
      zpos = fgkMyZpos;
      sprintf(cTagV,"UI%02d",iDet);
      sprintf(cTagM,"UC%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
      // Xe/Isobutane layer (drift volume) 
      zpos = fgkDrZpos;
      sprintf(cTagV,"UJ%02d",iDet);
      sprintf(cTagM,"UC%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
      // Upper part
      // Xe/Isobutane layer (amplification volume)
      zpos = fgkAmZpos;
      sprintf(cTagV,"UK%02d",iDet);
      sprintf(cTagM,"UE%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
      // Readout part
      // Cu layer (pad plane)
      zpos = fgkCuZpos; 
      sprintf(cTagV,"UL%02d",iDet);
      sprintf(cTagM,"UG%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
      // G10 layer (support structure)
      zpos = fgkSuZpos;
      sprintf(cTagV,"UM%02d",iDet);
      sprintf(cTagM,"UG%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
      // Cu layer (FEE + signal lines)
      zpos = fgkFeZpos; 
      sprintf(cTagV,"UN%02d",iDet);
      sprintf(cTagM,"UG%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");

      // Position the inner volumes of the chambers in the frames
      xpos      = 0.0;
      ypos      = 0.0;
      zpos      = 0.0;
      // The inside of the lower G10 frame
      sprintf(cTagV,"UC%02d",iDet);
      sprintf(cTagM,"UB%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
      // The lower G10 frame inside the aluminum frame
      sprintf(cTagV,"UB%02d",iDet);
      sprintf(cTagM,"UA%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
      // The inside of the upper G10 frame
      sprintf(cTagV,"UE%02d",iDet);
      sprintf(cTagM,"UD%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
      // The inside of the upper aluminum frame
      sprintf(cTagV,"UG%02d",iDet);
      sprintf(cTagM,"UF%02d",iDet);
      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");      

      // Position the frames of the chambers in the TRD mother volume
      xpos  = 0.;
      ypos  = - fClength[iplan][0] - fClength[iplan][1] - fClength[iplan][2]/2.;
      for (Int_t ic = 0; ic < icham; ic++) {
        ypos += fClength[iplan][ic];        
      }
      ypos += fClength[iplan][icham]/2.;
      zpos  = fgkCraH/2. + fgkCdrH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace);
      // The lower aluminum frame, radiator + drift region
      sprintf(cTagV,"UA%02d",iDet);
      gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
      // The upper G10 frame, amplification region
      sprintf(cTagV,"UD%02d",iDet);
      zpos += fgkCamH/2. + fgkCraH/2. + fgkCdrH/2.;
      gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
      // The upper aluminum frame
      sprintf(cTagV,"UF%02d",iDet);
      zpos += fgkCroH/2. + fgkCamH/2.;
      gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
 
      // Position the MCM volumina
      PositionReadout(iplan,icham);

      // Position the volumina for the cooling
      PositionCooling(iplan,icham,idrotm);

    }
  }

  xpos = 0.;
  ypos = 0.;
  zpos = 0.;
  gMC->Gspos("UTR1",1,"BTR1",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UTR1",2,"BTR2",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UTR1",3,"BTR3",xpos,ypos,zpos,0,"ONLY");

}

//_____________________________________________________________________________
void AliTRDgeometryDetail::CreateReadout(Int_t *idtmed) const
{
  //
  // Create the volumina of the readout electronics
  //

  const Int_t   kNparBox    = 3;

  Float_t parBox[kNparBox];
  Float_t xpos = 0.0;
  Float_t ypos = 0.0;
  Float_t zpos = 0.0;

  // The mother volume for the MCMs + connectors (air)
  parBox[0] = 3.0/2.;
  parBox[1] = 3.4/2.;
  parBox[2] = 0.5/2.;
  gMC->Gsvolu("UMCM","BOX",idtmed[1302-1],parBox,kNparBox);

  // The MCM carrier G10 layer
  parBox[0] = 3.0/2.;
  parBox[1] = 3.0/2.;
  parBox[2] = 0.1/2.;
  gMC->Gsvolu("UMC1","BOX",idtmed[1319-1],parBox,kNparBox);
  // The MCM carrier Cu layer
  parBox[0] = 3.0/2.;
  parBox[1] = 3.0/2.;
  parBox[2] = 0.0034/2.;
  gMC->Gsvolu("UMC2","BOX",idtmed[1318-1],parBox,kNparBox);
  // The MCM carrier Sn layer
  parBox[0] = 3.0/2.;
  parBox[1] = 3.0/2.;
  parBox[2] = 0.004/2.;
  gMC->Gsvolu("UMC3","BOX",idtmed[1317-1],parBox,kNparBox);
  // The MCM carrier Al layer
  parBox[0] = 3.0/2.;
  parBox[1] = 3.0/2.;
  parBox[2] = 0.05/2.;
  gMC->Gsvolu("UMC4","BOX",idtmed[1316-1],parBox,kNparBox);

  // The epoxy of chip no.1
  parBox[0] = 0.548/2.;
  parBox[1] = 0.548/2.;
  parBox[2] = 0.1/2.;
  gMC->Gsvolu("UCE1","BOX",idtmed[1321-1],parBox,kNparBox);
  // The silicon of chip no.1
  parBox[0] = 0.316/2.;
  parBox[1] = 0.316/2.;
  parBox[2] = 0.03/2.;
  gMC->Gsvolu("UCS1","BOX",idtmed[1320-1],parBox,kNparBox);

  // The epoxy of chip no.2
  parBox[0] = 1.549/2.;
  parBox[1] = 1.549/2.;
  parBox[2] = 0.1/2.;
  gMC->Gsvolu("UCE2","BOX",idtmed[1321-1],parBox,kNparBox);
  // The silicon of chip no.2
  parBox[0] = 0.894/2.;
  parBox[1] = 0.894/2.;
  parBox[2] = 0.03/2.;
  gMC->Gsvolu("UCS2","BOX",idtmed[1320-1],parBox,kNparBox);

  // The PE of the connector
  parBox[0] = 2.25/2.;
  parBox[1] = 0.4/2.;
  parBox[2] = 0.3/2.;
  gMC->Gsvolu("UCN1","BOX",idtmed[1322-1],parBox,kNparBox);
  // The Cu of the connector
  parBox[0] = 2.25/2.;
  parBox[1] = 0.4/2.;
  parBox[2] = 0.005/2.;
  gMC->Gsvolu("UCN2","BOX",idtmed[1323-1],parBox,kNparBox);

  xpos  =  0.0;
  ypos  = -0.4/2.;
  zpos  = -0.25      + 0.1/2.;
  gMC->Gspos("UMC1",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
  zpos +=  0.1/2.   + 0.0034/2.;
  gMC->Gspos("UMC2",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
  zpos +=  0.0034/2 + 0.004/2.;
  gMC->Gspos("UMC3",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
  zpos +=  0.004/2  + 0.05/2.;
  gMC->Gspos("UMC4",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
  zpos +=  0.05/2.  + 0.1/2.;
  xpos  =  1.0;
  gMC->Gspos("UCE1",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
  xpos  = -0.5;
  gMC->Gspos("UCE2",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
  zpos +=  0.1/2.   + 0.03/2.;
  xpos  =  1.0;
  gMC->Gspos("UCS1",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
  xpos  = -0.5;
  gMC->Gspos("UCS2",1,"UMCM",xpos,ypos,zpos,0,"ONLY");  
  xpos  =  0.0;
  ypos  =  3.4/2.   - 0.4/2.;
  zpos  = -0.25     + 0.3/2.;
  gMC->Gspos("UCN1",1,"UMCM",xpos,ypos,zpos,0,"ONLY");  
  zpos +=  0.3/2.   + 0.005/2.;
  gMC->Gspos("UCN2",1,"UMCM",xpos,ypos,zpos,0,"ONLY");  

}

//_____________________________________________________________________________
void AliTRDgeometryDetail::PositionReadout(Int_t ipla, Int_t icha)
{
  //
  // Position the volumina inside the readout mother volume
  //

  const Int_t   kNmcmChannel = 18;

  AliTRDparameter *parameter = new AliTRDparameter();

  Int_t nMCMrow = parameter->GetRowMax(ipla,icha,0);
  Int_t nMCMcol = parameter->GetColMax(ipla) / kNmcmChannel;

  Float_t xSize = (GetChamberWidth(ipla)       - 2.*fgkCpadW) 
                / ((Float_t) nMCMcol);
  Float_t ySize = (GetChamberLength(ipla,icha) - 2.*fgkRpadW) 
                / ((Float_t) nMCMrow);
  Float_t x0    = parameter->GetCol0(ipla);
  Float_t y0    = parameter->GetRow0(ipla,icha,0);

  Int_t iCopy = GetDetector(ipla,icha,0) * 1000;
  for (Int_t iMCMrow = 0; iMCMrow < nMCMrow; iMCMrow++) {
    for (Int_t iMCMcol = 0; iMCMcol < nMCMcol; iMCMcol++) {
      iCopy++;
      Float_t xpos = (0.5 + iMCMcol) * xSize + x0; 
      Float_t ypos = (0.5 + iMCMrow) * ySize + y0;
      Float_t zpos = fgkCH - fgkSheight/2. + 0.5/2.
                   + ipla * (fgkCH + fgkVspace);
      gMC->Gspos("UMCM",iCopy,"UTR1",xpos,ypos,zpos,0,"ONLY");    
    }
  }

  delete parameter;

}

//_____________________________________________________________________________
void AliTRDgeometryDetail::CreateCooling(Int_t *idtmed) const
{
  //
  // Create the volumina of the cooling
  //

  const Int_t kNparTube = 3;

  Float_t parTube[kNparTube];
  Float_t xpos;
  Float_t ypos;
  Float_t zpos;

  // The aluminum pipe for the cooling
  parTube[0] = 0.0;
  parTube[1] = 0.0;
  parTube[2] = 0.0;
  gMC->Gsvolu("UCOA","TUBE",idtmed[1324-1],parTube,0);

  // The cooling water
  parTube[0] =  0.0;
  parTube[1] =  0.2/2.;
  parTube[2] = -1.;
  gMC->Gsvolu("UCOW","TUBE",idtmed[1314-1],parTube,kNparTube);

  // Water inside the cooling pipe
  xpos = 0.0;
  ypos = 0.0;
  zpos = 0.0;
  gMC->Gspos("UCOW",1,"UCOA",xpos,ypos,zpos,0,"ONLY");

}

//_____________________________________________________________________________
void AliTRDgeometryDetail::PositionCooling(Int_t ipla, Int_t icha, Int_t idrotm)
{
  //
  // Position the volumina of the cooling
  //

  const Int_t kNpar = 3;

  Float_t par[kNpar];
  Float_t xpos;
  Float_t ypos;
  Float_t zpos;

  AliTRDparameter *parameter = new AliTRDparameter();

  Int_t   iCopy   = GetDetector(ipla,icha,0) * 100;
  Int_t   nMCMrow = parameter->GetRowMax(ipla,icha,0);

  Float_t ySize   = (GetChamberLength(ipla,icha) - 2.*fgkRpadW) 
                  / ((Float_t) nMCMrow);
  Float_t y0      = parameter->GetRow0(ipla,icha,0);

  // Position the cooling pipes
  for (Int_t iMCMrow = 0; iMCMrow < nMCMrow; iMCMrow++) {

    xpos   = 0.0;
    ypos   = (0.5 + iMCMrow) * ySize + y0 - 1.9;
    zpos   = fgkCH - fgkSheight/2. + 0.5/2.
                   + ipla * (fgkCH + fgkVspace);
    par[0] = 0.0;
    par[1] = 0.3/2.;
    par[2] = GetChamberWidth(ipla)/2.+ fgkCroW;
    gMC->Gsposp("UCOA",iCopy+iMCMrow,"UTR1",xpos,ypos,zpos
                      ,idrotm,"ONLY",par,kNpar);

  }

  delete parameter;

}
Commit	Line	Data
16bf9884	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$
4c475d27	18	Revision 1.6 2002/10/14 14:57:43 hristov
	19	Merging the VirtualMC branch to the main development branch (HEAD)
	20
b9d0a01d	21	Revision 1.3.6.2 2002/07/24 10:09:30 alibrary
	22	Updating VirtualMC
	23
	24	Revision 1.5 2002/06/12 09:54:35 cblume
	25	Update of tracking code provided by Sergei
	26
5443e65e	27	Revision 1.4 2002/03/28 14:59:07 cblume
	28	Coding conventions
	29
0a29d0f1	30	Revision 1.3 2002/02/11 14:21:16 cblume
	31	Update of the geometry. Get rid of MANY
	32
0a770ac9	33	Revision 1.2 2001/11/08 13:13:08 cblume
	34	Change to MANY for UCFI/M/O and UAFI/M/O
	35
a0bc1a05	36	Revision 1.1 2001/11/06 17:19:41 cblume
	37	Add detailed geometry and simple simulator
	38
16bf9884	39	*/
	40
	41	///////////////////////////////////////////////////////////////////////////////
	42	// //
0a770ac9	43	// Detailed TRD geometry for the spaceframe without holes //
16bf9884	44	// //
	45	///////////////////////////////////////////////////////////////////////////////
	46
4c475d27	47	#include "TVirtualMC.h"
16bf9884	48
16bf9884	49	#include "AliTRDgeometryDetail.h"
5443e65e	50	#include "AliTRDparameter.h"
16bf9884	51
	52	ClassImp(AliTRDgeometryDetail)
	53
	54	//_____________________________________________________________________________
	55	AliTRDgeometryDetail::AliTRDgeometryDetail():AliTRDgeometryFull()
	56	{
	57	//
	58	// AliTRDgeometryDetail default constructor
	59	//
	60
	61	Init();
	62
	63	}
	64
	65	//_____________________________________________________________________________
	66	AliTRDgeometryDetail::~AliTRDgeometryDetail()
	67	{
	68	//
	69	// AliTRDgeometryDetail destructor
	70	//
	71
	72	}
	73
	74	//_____________________________________________________________________________
	75	void AliTRDgeometryDetail::Init()
	76	{
	77	//
	78	// Initializes the geometry parameter
	79	//
	80
	81	AliTRDgeometryFull::Init();
	82
	83	}
	84
	85	//_____________________________________________________________________________
	86	void AliTRDgeometryDetail::CreateGeometry(Int_t *idtmed)
	87	{
	88	//
	89	// Create the detailed TRD geometry without hole
0a770ac9	90	// including the MCMs and the cooling pipes
	91	//
	92	//
	93	// Names of the TRD volumina (xx = detector number):
	94	//
	95	// Lower part of the readout chambers (gas volume + radiator)
	96	//
	97	// UAxx Aluminum frames (Al)
	98	// UBxx G10 frames (C)
	99	// UCxx Inner volumes (Air)
	100	//
	101	// Upper part of the readout chambers (readout plane + fee)
	102	//
	103	// UDxx G10 frames (C)
	104	// UExx Inner volumes of the G10 (Air)
	105	// UFxx Aluminum frames (Al)
	106	// UGxx Inner volumes of the Al (Air)
	107	//
	108	// Inner material layers
	109	//
	110	// UHxx Radiator (Rohacell)
	111	// UIxx Entrance window (Mylar)
	112	// UJxx Drift volume (Xe/CO2)
	113	// UKxx Amplification volume (Xe/CO2)
	114	// ULxx Pad plane (Cu)
	115	// UMxx Support structure (Rohacell)
	116	// UNxx FEE + signal lines (Cu)
16bf9884	117	//
16bf9884	118
	119	const Int_t kNparTrd = 4;
	120	const Int_t kNparCha = 3;
16bf9884	121
0a770ac9	122	Float_t xpos, ypos, zpos;
0a770ac9	123
16bf9884	124	Float_t parTrd[kNparTrd];
	125	Float_t parCha[kNparCha];
	126
0a770ac9	127	Char_t cTagV[5];
0a770ac9	128	Char_t cTagM[5];
16bf9884	129
0a770ac9	130	Int_t idrotm;
16bf9884	131
0a770ac9	132	// Rotation matrix
0a770ac9	133	gMC->Matrix(idrotm, 0.0, 0.0, 90.0, 90.0, 90.0, 0.0);
16bf9884	134
	135	// The TRD mother volume for one sector (Air), full length in z-direction
	136	parTrd[0] = fgkSwidth1/2.;
	137	parTrd[1] = fgkSwidth2/2.;
	138	parTrd[2] = fgkSlenTR1/2.;
	139	parTrd[3] = fgkSheight/2.;
0a770ac9	140	gMC->Gsvolu("UTR1","TRD1",idtmed[1302-1],parTrd,kNparTrd);
	141
	142	// Create the readout volumina
	143	CreateReadout(idtmed);
	144
	145	// Create the volumina for the cooling
	146	CreateCooling(idtmed);
16bf9884	147
0a770ac9	148	for (Int_t icham = 0; icham < kNcham; icham++) {
	149	for (Int_t iplan = 0; iplan < kNplan; iplan++) {
	150
	151	Int_t iDet = GetDetectorSec(iplan,icham);
	152
	153	// The lower part of the readout chambers (gas volume + radiator)
	154	// The aluminum frames
	155	sprintf(cTagV,"UA%02d",iDet);
	156	parCha[0] = fCwidth[iplan]/2.;
	157	parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
	158	parCha[2] = fgkCraH/2. + fgkCdrH/2.;
	159	gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
	160	// The G10 frames
	161	sprintf(cTagV,"UB%02d",iDet);
	162	parCha[0] = fCwidth[iplan]/2. - fgkCalT;
	163	parCha[1] = -1.;
	164	parCha[2] = -1.;
	165	gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
	166	// The inner part (air)
	167	sprintf(cTagV,"UC%02d",iDet);
	168	parCha[0] = fCwidth[iplan]/2. - fgkCalT - fgkCclsT;
	169	parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCclfT;
	170	parCha[2] = -1.;
	171	gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
	172
	173	// The upper part of the readout chambers (readout plane + fee)
	174	// The G10 frames
	175	sprintf(cTagV,"UD%02d",iDet);
	176	parCha[0] = fCwidth[iplan]/2. + fgkCroW;
	177	parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
	178	parCha[2] = fgkCamH/2.;
	179	gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
	180	// The inner part of the G10 frame (air)
	181	sprintf(cTagV,"UE%02d",iDet);
	182	parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCcuT;
	183	parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCcuT;
	184	parCha[2] = -1.;
	185	gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
	186	// The aluminum frames
	187	sprintf(cTagV,"UF%02d",iDet);
	188	parCha[0] = fCwidth[iplan]/2. + fgkCroW;
	189	parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
	190	parCha[2] = fgkCroH/2.;
	191	gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
	192	// The inner part of the aluminum frames
	193	sprintf(cTagV,"UG%02d",iDet);
	194	parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCauT;
	195	parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCauT;
	196	parCha[2] = -1.;
	197	gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
	198
	199	// The material layers inside the chambers
	200	parCha[0] = -1.;
	201	parCha[1] = -1.;
	202	// Rohacell layer (radiator)
	203	parCha[2] = fgkRaThick/2;
	204	sprintf(cTagV,"UH%02d",iDet);
	205	gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha);
	206	// Mylar layer (entrance window + HV cathode)
	207	parCha[2] = fgkMyThick/2;
	208	sprintf(cTagV,"UI%02d",iDet);
	209	gMC->Gsvolu(cTagV,"BOX ",idtmed[1308-1],parCha,kNparCha);
	210	// Xe/Isobutane layer (drift volume)
	211	parCha[2] = fgkDrThick/2.;
212	sprintf(cTagV,"UJ%02d",iDet);
213	gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
214	// Xe/Isobutane layer (amplification volume)
215	parCha[2] = fgkAmThick/2.;
216	sprintf(cTagV,"UK%02d",iDet);
217	gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
218	// Cu layer (pad plane)
219	parCha[2] = fgkCuThick/2;
220	sprintf(cTagV,"UL%02d",iDet);
221	gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
222	// G10 layer (support structure / honeycomb)
223	parCha[2] = fgkSuThick/2;
224	sprintf(cTagV,"UM%02d",iDet);
225	gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha);
226	// Cu layer (FEE + signal lines)
227	parCha[2] = fgkFeThick/2;
228	sprintf(cTagV,"UN%02d",iDet);
229	gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
230
231	// Position the layers in the chambers
232	xpos = 0;
233	ypos = 0;
234	// Lower part
235	// Rohacell layer (radiator)
236	zpos = fgkRaZpos;
237	sprintf(cTagV,"UH%02d",iDet);
238	sprintf(cTagM,"UC%02d",iDet);
239	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
240	// Mylar layer (entrance window + HV cathode)
241	zpos = fgkMyZpos;
242	sprintf(cTagV,"UI%02d",iDet);
243	sprintf(cTagM,"UC%02d",iDet);
244	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
245	// Xe/Isobutane layer (drift volume)
246	zpos = fgkDrZpos;
247	sprintf(cTagV,"UJ%02d",iDet);
248	sprintf(cTagM,"UC%02d",iDet);
249	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
250	// Upper part
251	// Xe/Isobutane layer (amplification volume)
252	zpos = fgkAmZpos;
253	sprintf(cTagV,"UK%02d",iDet);
254	sprintf(cTagM,"UE%02d",iDet);
255	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
256	// Readout part
257	// Cu layer (pad plane)
258	zpos = fgkCuZpos;
259	sprintf(cTagV,"UL%02d",iDet);
260	sprintf(cTagM,"UG%02d",iDet);
261	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
262	// G10 layer (support structure)
263	zpos = fgkSuZpos;
264	sprintf(cTagV,"UM%02d",iDet);
265	sprintf(cTagM,"UG%02d",iDet);
266	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
267	// Cu layer (FEE + signal lines)
268	zpos = fgkFeZpos;
269	sprintf(cTagV,"UN%02d",iDet);
270	sprintf(cTagM,"UG%02d",iDet);
271	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
272
273	// Position the inner volumes of the chambers in the frames
274	xpos = 0.0;
275	ypos = 0.0;
276	zpos = 0.0;
277	// The inside of the lower G10 frame
278	sprintf(cTagV,"UC%02d",iDet);
279	sprintf(cTagM,"UB%02d",iDet);
280	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
281	// The lower G10 frame inside the aluminum frame
282	sprintf(cTagV,"UB%02d",iDet);
283	sprintf(cTagM,"UA%02d",iDet);
284	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
285	// The inside of the upper G10 frame
286	sprintf(cTagV,"UE%02d",iDet);
287	sprintf(cTagM,"UD%02d",iDet);
288	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
289	// The inside of the upper aluminum frame
290	sprintf(cTagV,"UG%02d",iDet);
291	sprintf(cTagM,"UF%02d",iDet);
292	gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
293
294	// Position the frames of the chambers in the TRD mother volume
295	xpos = 0.;
296	ypos = - fClength[iplan][0] - fClength[iplan][1] - fClength[iplan][2]/2.;
297	for (Int_t ic = 0; ic < icham; ic++) {
298	ypos += fClength[iplan][ic];
299	}
300	ypos += fClength[iplan][icham]/2.;
301	zpos = fgkCraH/2. + fgkCdrH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace);
302	// The lower aluminum frame, radiator + drift region
303	sprintf(cTagV,"UA%02d",iDet);
304	gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
305	// The upper G10 frame, amplification region
306	sprintf(cTagV,"UD%02d",iDet);
307	zpos += fgkCamH/2. + fgkCraH/2. + fgkCdrH/2.;
308	gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
309	// The upper aluminum frame
310	sprintf(cTagV,"UF%02d",iDet);
311	zpos += fgkCroH/2. + fgkCamH/2.;
312	gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
313
314	// Position the MCM volumina
315	PositionReadout(iplan,icham);
316
317	// Position the volumina for the cooling
318	PositionCooling(iplan,icham,idrotm);
319
320	}
16bf9884	321	}
	322
	323	xpos = 0.;
	324	ypos = 0.;
	325	zpos = 0.;
0a770ac9	326	gMC->Gspos("UTR1",1,"BTR1",xpos,ypos,zpos,0,"ONLY");
	327	gMC->Gspos("UTR1",2,"BTR2",xpos,ypos,zpos,0,"ONLY");
	328	gMC->Gspos("UTR1",3,"BTR3",xpos,ypos,zpos,0,"ONLY");
16bf9884	329
	330	}
	331
	332	//_____________________________________________________________________________
0a29d0f1	333	void AliTRDgeometryDetail::CreateReadout(Int_t *idtmed) const
16bf9884	334	{
	335	//
	336	// Create the volumina of the readout electronics
	337	//
	338
	339	const Int_t kNparBox = 3;
	340
	341	Float_t parBox[kNparBox];
	342	Float_t xpos = 0.0;
	343	Float_t ypos = 0.0;
	344	Float_t zpos = 0.0;
	345
	346	// The mother volume for the MCMs + connectors (air)
	347	parBox[0] = 3.0/2.;
	348	parBox[1] = 3.4/2.;
	349	parBox[2] = 0.5/2.;
	350	gMC->Gsvolu("UMCM","BOX",idtmed[1302-1],parBox,kNparBox);
	351
	352	// The MCM carrier G10 layer
	353	parBox[0] = 3.0/2.;
	354	parBox[1] = 3.0/2.;
	355	parBox[2] = 0.1/2.;
	356	gMC->Gsvolu("UMC1","BOX",idtmed[1319-1],parBox,kNparBox);
	357	// The MCM carrier Cu layer
	358	parBox[0] = 3.0/2.;
	359	parBox[1] = 3.0/2.;
	360	parBox[2] = 0.0034/2.;
	361	gMC->Gsvolu("UMC2","BOX",idtmed[1318-1],parBox,kNparBox);
	362	// The MCM carrier Sn layer
	363	parBox[0] = 3.0/2.;
	364	parBox[1] = 3.0/2.;
	365	parBox[2] = 0.004/2.;
	366	gMC->Gsvolu("UMC3","BOX",idtmed[1317-1],parBox,kNparBox);
	367	// The MCM carrier Al layer
	368	parBox[0] = 3.0/2.;
	369	parBox[1] = 3.0/2.;
	370	parBox[2] = 0.05/2.;
	371	gMC->Gsvolu("UMC4","BOX",idtmed[1316-1],parBox,kNparBox);
	372
	373	// The epoxy of chip no.1
	374	parBox[0] = 0.548/2.;
	375	parBox[1] = 0.548/2.;
	376	parBox[2] = 0.1/2.;
	377	gMC->Gsvolu("UCE1","BOX",idtmed[1321-1],parBox,kNparBox);
	378	// The silicon of chip no.1
	379	parBox[0] = 0.316/2.;
	380	parBox[1] = 0.316/2.;
	381	parBox[2] = 0.03/2.;
	382	gMC->Gsvolu("UCS1","BOX",idtmed[1320-1],parBox,kNparBox);
	383
	384	// The epoxy of chip no.2
	385	parBox[0] = 1.549/2.;
	386	parBox[1] = 1.549/2.;
	387	parBox[2] = 0.1/2.;
	388	gMC->Gsvolu("UCE2","BOX",idtmed[1321-1],parBox,kNparBox);
	389	// The silicon of chip no.2
	390	parBox[0] = 0.894/2.;
	391	parBox[1] = 0.894/2.;
	392	parBox[2] = 0.03/2.;
	393	gMC->Gsvolu("UCS2","BOX",idtmed[1320-1],parBox,kNparBox);
	394
	395	// The PE of the connector
	396	parBox[0] = 2.25/2.;
	397	parBox[1] = 0.4/2.;
398	parBox[2] = 0.3/2.;
399	gMC->Gsvolu("UCN1","BOX",idtmed[1322-1],parBox,kNparBox);
400	// The Cu of the connector
401	parBox[0] = 2.25/2.;
402	parBox[1] = 0.4/2.;
403	parBox[2] = 0.005/2.;
404	gMC->Gsvolu("UCN2","BOX",idtmed[1323-1],parBox,kNparBox);
405
406	xpos = 0.0;
407	ypos = -0.4/2.;
408	zpos = -0.25 + 0.1/2.;
409	gMC->Gspos("UMC1",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
410	zpos += 0.1/2. + 0.0034/2.;
411	gMC->Gspos("UMC2",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
412	zpos += 0.0034/2 + 0.004/2.;
413	gMC->Gspos("UMC3",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
414	zpos += 0.004/2 + 0.05/2.;
415	gMC->Gspos("UMC4",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
416	zpos += 0.05/2. + 0.1/2.;
417	xpos = 1.0;
418	gMC->Gspos("UCE1",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
419	xpos = -0.5;
420	gMC->Gspos("UCE2",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
421	zpos += 0.1/2. + 0.03/2.;
422	xpos = 1.0;
423	gMC->Gspos("UCS1",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
424	xpos = -0.5;
425	gMC->Gspos("UCS2",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
426	xpos = 0.0;
427	ypos = 3.4/2. - 0.4/2.;
428	zpos = -0.25 + 0.3/2.;
429	gMC->Gspos("UCN1",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
430	zpos += 0.3/2. + 0.005/2.;
431	gMC->Gspos("UCN2",1,"UMCM",xpos,ypos,zpos,0,"ONLY");
432
433	}
434
435	//_____________________________________________________________________________
436	void AliTRDgeometryDetail::PositionReadout(Int_t ipla, Int_t icha)
437	{
438	//
439	// Position the volumina inside the readout mother volume
440	//
441
0a770ac9	442	const Int_t kNmcmChannel = 18;
16bf9884	443
5443e65e	444	AliTRDparameter *parameter = new AliTRDparameter();
	445
	446	Int_t nMCMrow = parameter->GetRowMax(ipla,icha,0);
	447	Int_t nMCMcol = parameter->GetColMax(ipla) / kNmcmChannel;
0a770ac9	448
	449	Float_t xSize = (GetChamberWidth(ipla) - 2.*fgkCpadW)
	450	/ ((Float_t) nMCMcol);
	451	Float_t ySize = (GetChamberLength(ipla,icha) - 2.*fgkRpadW)
	452	/ ((Float_t) nMCMrow);
5443e65e	453	Float_t x0 = parameter->GetCol0(ipla);
5443e65e	454	Float_t y0 = parameter->GetRow0(ipla,icha,0);
16bf9884	455
	456	Int_t iCopy = GetDetector(ipla,icha,0) * 1000;
	457	for (Int_t iMCMrow = 0; iMCMrow < nMCMrow; iMCMrow++) {
	458	for (Int_t iMCMcol = 0; iMCMcol < nMCMcol; iMCMcol++) {
	459	iCopy++;
	460	Float_t xpos = (0.5 + iMCMcol) * xSize + x0;
	461	Float_t ypos = (0.5 + iMCMrow) * ySize + y0;
0a770ac9	462	Float_t zpos = fgkCH - fgkSheight/2. + 0.5/2.
	463	+ ipla * (fgkCH + fgkVspace);
	464	gMC->Gspos("UMCM",iCopy,"UTR1",xpos,ypos,zpos,0,"ONLY");
16bf9884	465	}
	466	}
	467
5443e65e	468	delete parameter;
5443e65e	469
16bf9884	470	}
	471
	472	//_____________________________________________________________________________
0a29d0f1	473	void AliTRDgeometryDetail::CreateCooling(Int_t *idtmed) const
16bf9884	474	{
	475	//
	476	// Create the volumina of the cooling
	477	//
	478
0a770ac9	479	const Int_t kNparTube = 3;
16bf9884	480
0a770ac9	481	Float_t parTube[kNparTube];
	482	Float_t xpos;
	483	Float_t ypos;
	484	Float_t zpos;
16bf9884	485
16bf9884	486	// The aluminum pipe for the cooling
0a770ac9	487	parTube[0] = 0.0;
	488	parTube[1] = 0.0;
	489	parTube[2] = 0.0;
	490	gMC->Gsvolu("UCOA","TUBE",idtmed[1324-1],parTube,0);
16bf9884	491
16bf9884	492	// The cooling water
0a770ac9	493	parTube[0] = 0.0;
	494	parTube[1] = 0.2/2.;
	495	parTube[2] = -1.;
	496	gMC->Gsvolu("UCOW","TUBE",idtmed[1314-1],parTube,kNparTube);
	497
	498	// Water inside the cooling pipe
	499	xpos = 0.0;
	500	ypos = 0.0;
	501	zpos = 0.0;
16bf9884	502	gMC->Gspos("UCOW",1,"UCOA",xpos,ypos,zpos,0,"ONLY");
	503
	504	}
	505
	506	//_____________________________________________________________________________
0a770ac9	507	void AliTRDgeometryDetail::PositionCooling(Int_t ipla, Int_t icha, Int_t idrotm)
16bf9884	508	{
	509	//
	510	// Position the volumina of the cooling
	511	//
	512
0a770ac9	513	const Int_t kNpar = 3;
16bf9884	514
16bf9884	515	Float_t par[kNpar];
0a770ac9	516	Float_t xpos;
	517	Float_t ypos;
	518	Float_t zpos;
16bf9884	519
5443e65e	520	AliTRDparameter *parameter = new AliTRDparameter();
5443e65e	521
16bf9884	522	Int_t iCopy = GetDetector(ipla,icha,0) * 100;
5443e65e	523	Int_t nMCMrow = parameter->GetRowMax(ipla,icha,0);
16bf9884	524
0a770ac9	525	Float_t ySize = (GetChamberLength(ipla,icha) - 2.*fgkRpadW)
0a770ac9	526	/ ((Float_t) nMCMrow);
5443e65e	527	Float_t y0 = parameter->GetRow0(ipla,icha,0);
16bf9884	528
	529	// Position the cooling pipes
	530	for (Int_t iMCMrow = 0; iMCMrow < nMCMrow; iMCMrow++) {
	531
	532	xpos = 0.0;
0a770ac9	533	ypos = (0.5 + iMCMrow) * ySize + y0 - 1.9;
	534	zpos = fgkCH - fgkSheight/2. + 0.5/2.
	535	+ ipla * (fgkCH + fgkVspace);
	536	par[0] = 0.0;
16bf9884	537	par[1] = 0.3/2.;
0a770ac9	538	par[2] = GetChamberWidth(ipla)/2.+ fgkCroW;
	539	gMC->Gsposp("UCOA",iCopy+iMCMrow,"UTR1",xpos,ypos,zpos
	540	,idrotm,"ONLY",par,kNpar);
16bf9884	541
	542	}
	543
5443e65e	544	delete parameter;
5443e65e	545
16bf9884	546	}