[u/mrichter/AliRoot.git] / TRD / AliTRDgeometry.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$
*/

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  TRD geometry class                                                       //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "AliTRDgeometry.h"
#include "AliTRDrecPoint.h"

ClassImp(AliTRDgeometry)

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

  Init();

}

//_____________________________________________________________________________
AliTRDgeometry::~AliTRDgeometry()
{

}

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

  Int_t iplan;

  // The width of the chambers
  fCwidth[0]    =  99.6;
  fCwidth[1]    = 104.1;
  fCwidth[2]    = 108.5;
  fCwidth[3]    = 112.9;
  fCwidth[4]    = 117.4;
  fCwidth[5]    = 121.8;

  // The default pad dimensions
  fRowPadSize  = 4.5;
  fColPadSize  = 1.0;
  fTimeBinSize = 0.1;

  // 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 column (rphi-direction)  
  for (iplan = 0; iplan < kNplan; iplan++) {
    fColMax[iplan] = 1 + TMath::Nint((fCwidth[iplan] - 2. * kCcthick) 
                                                     / fColPadSize - 0.5);
    fCol0[iplan]   = -fCwidth[iplan]/2. + kCcthick;
  }

  // The time bucket
  fTimeMax = 1 + TMath::Nint(kDrThick / fTimeBinSize - 0.5);
  for (iplan = 0; iplan < kNplan; iplan++) {
    fTime0[iplan]  = kRmin + kCcframe/2. + kDrZpos - 0.5 * kDrThick
                           + iplan * (kCheight + kCspace);
  } 

}

//_____________________________________________________________________________
void AliTRDgeometry::CreateGeometry(Int_t *idtmed)
{
  //
  // Create the TRD geometry
  //
  // Author: Christoph Blume (C.Blume@gsi.de) 20/07/99
  //
  // The volumes:
  //    TRD1-3     (Air)   --- The TRD mother volumes for one sector. 
  //                           To be placed into the spaceframe.
  //
  //    UAFI(/M/O) (Al)    --- The aluminum frame of the inner(/middle/outer) chambers (readout)
  //    UCFI(/M/O) (C)     --- The carbon frame of the inner(/middle/outer) chambers 
  //                           (driftchamber + radiator)
  //    UAII(/M/O) (Air)   --- The inner part of the readout of the inner(/middle/outer) chambers
  //    UFII(/M/O) (Air)   --- The inner part of the chamner and radiator of the 
  //                           inner(/middle/outer) chambers
  //
  // The material layers in one chamber:
  //    UL01       (G10)   --- The gas seal of the radiator
  //    UL02       (CO2)   --- The gas in the radiator
  //    UL03       (PE)    --- The foil stack
  //    UL04       (Mylar) --- Entrance window to the driftvolume and HV-cathode
  //    UL05       (Xe)    --- The driftvolume
  //    UL06       (Xe)    --- The amplification region
  //    
  //    UL07       (Cu)    --- The pad plane
  //    UL08       (G10)   --- The Nomex honeycomb support structure
  //    UL09       (Cu)    --- FEE and signal lines
  //    UL10       (PE)    --- The cooling devices
  //    UL11       (Water) --- The cooling water

  const Int_t npar_cha = 3;

  Float_t par_dum[3];
  Float_t par_cha[npar_cha];

  Float_t xpos, ypos, zpos;

  // The aluminum frames - readout + electronics (Al)
  // The inner chambers
  gMC->Gsvolu("UAFI","BOX ",idtmed[1301-1],par_dum,0);
  // The middle chambers
  gMC->Gsvolu("UAFM","BOX ",idtmed[1301-1],par_dum,0);
  // The outer chambers
  gMC->Gsvolu("UAFO","BOX ",idtmed[1301-1],par_dum,0);

  // The inner part of the aluminum frames (Air)
  // The inner chambers
  gMC->Gsvolu("UAII","BOX ",idtmed[1302-1],par_dum,0);
  // The middle chambers
  gMC->Gsvolu("UAIM","BOX ",idtmed[1302-1],par_dum,0);
  // The outer chambers
  gMC->Gsvolu("UAIO","BOX ",idtmed[1302-1],par_dum,0);

  // The carbon frames - radiator + driftchamber (C)
  // The inner chambers
  gMC->Gsvolu("UCFI","BOX ",idtmed[1307-1],par_dum,0);
  // The middle chambers
  gMC->Gsvolu("UCFM","BOX ",idtmed[1307-1],par_dum,0);
  // The outer chambers
  gMC->Gsvolu("UCFO","BOX ",idtmed[1307-1],par_dum,0);

  // The inner part of the carbon frames (Air)
  // The inner chambers
  gMC->Gsvolu("UCII","BOX ",idtmed[1302-1],par_dum,0);
  // The middle chambers
  gMC->Gsvolu("UCIM","BOX ",idtmed[1302-1],par_dum,0);
  // The outer chambers
  gMC->Gsvolu("UCIO","BOX ",idtmed[1302-1],par_dum,0);

  // The material layers inside the chambers
  par_cha[0] = -1.;
  par_cha[1] = -1.;
  // G10 layer (radiator seal)
  par_cha[2] = kSeThick/2;
  gMC->Gsvolu("UL01","BOX ",idtmed[1313-1],par_cha,npar_cha);
  // CO2 layer (radiator)
  par_cha[2] = kRaThick/2;
  gMC->Gsvolu("UL02","BOX ",idtmed[1312-1],par_cha,npar_cha);
  // PE layer (radiator)
  par_cha[2] = kPeThick/2;
  gMC->Gsvolu("UL03","BOX ",idtmed[1303-1],par_cha,npar_cha);
  // Mylar layer (entrance window + HV cathode) 
  par_cha[2] = kMyThick/2;
  gMC->Gsvolu("UL04","BOX ",idtmed[1308-1],par_cha,npar_cha);
  // Xe/Isobutane layer (drift volume, sensitive) 
  par_cha[2] = kDrThick/2.;
  gMC->Gsvolu("UL05","BOX ",idtmed[1309-1],par_cha,npar_cha);
  // Xe/Isobutane layer (amplification volume, not sensitive)
  par_cha[2] = kAmThick/2.;
  gMC->Gsvolu("UL06","BOX ",idtmed[1309-1],par_cha,npar_cha);
  
  // Cu layer (pad plane)
  par_cha[2] = kCuThick/2;
  gMC->Gsvolu("UL07","BOX ",idtmed[1305-1],par_cha,npar_cha);
  // G10 layer (support structure)
  par_cha[2] = kSuThick/2;
  gMC->Gsvolu("UL08","BOX ",idtmed[1313-1],par_cha,npar_cha);
  // Cu layer (FEE + signal lines)
  par_cha[2] = kFeThick/2;
  gMC->Gsvolu("UL09","BOX ",idtmed[1305-1],par_cha,npar_cha);
  // PE layer (cooling devices)
  par_cha[2] = kCoThick/2;
  gMC->Gsvolu("UL10","BOX ",idtmed[1303-1],par_cha,npar_cha);
  // Water layer (cooling)
  par_cha[2] = kWaThick/2;
  gMC->Gsvolu("UL11","BOX ",idtmed[1314-1],par_cha,npar_cha);

  // Position the layers in the chambers
  xpos = 0;
  ypos = 0;

  // G10 layer (radiator seal)
  zpos = kSeZpos;
  gMC->Gspos("UL01",1,"UCII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL01",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL01",3,"UCIO",xpos,ypos,zpos,0,"ONLY");
  // CO2 layer (radiator)
  zpos = kRaZpos;
  gMC->Gspos("UL02",1,"UCII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL02",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL02",3,"UCIO",xpos,ypos,zpos,0,"ONLY");
  // PE layer (radiator)
  zpos = 0;
  gMC->Gspos("UL03",1,"UL02",xpos,ypos,zpos,0,"ONLY");
  // Mylar layer (entrance window + HV cathode)   
  zpos = kMyZpos;
  gMC->Gspos("UL04",1,"UCII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL04",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL04",3,"UCIO",xpos,ypos,zpos,0,"ONLY");
  // Xe/Isobutane layer (drift volume) 
  zpos = kDrZpos;
  gMC->Gspos("UL05",1,"UCII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL05",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL05",3,"UCIO",xpos,ypos,zpos,0,"ONLY");
  // Xe/Isobutane layer (amplification volume)
  zpos = kAmZpos;
  gMC->Gspos("UL06",1,"UCII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL06",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL06",3,"UCIO",xpos,ypos,zpos,0,"ONLY");

  // Cu layer (pad plane)
  zpos = kCuZpos;
  gMC->Gspos("UL07",1,"UAII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL07",2,"UAIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL07",3,"UAIO",xpos,ypos,zpos,0,"ONLY");
  // G10 layer (support structure)
  zpos = kSuZpos;
  gMC->Gspos("UL08",1,"UAII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL08",2,"UAIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL08",3,"UAIO",xpos,ypos,zpos,0,"ONLY");
  // Cu layer (FEE + signal lines)
  zpos = kFeZpos; 
  gMC->Gspos("UL09",1,"UAII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL09",2,"UAIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL09",3,"UAIO",xpos,ypos,zpos,0,"ONLY");
  // PE layer (cooling devices)
  zpos = kCoZpos;
  gMC->Gspos("UL10",1,"UAII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL10",2,"UAIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL10",3,"UAIO",xpos,ypos,zpos,0,"ONLY");
  // Water layer (cooling)
  zpos = kWaZpos;
  gMC->Gspos("UL11",1,"UAII",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL11",1,"UAIM",xpos,ypos,zpos,0,"ONLY");
  gMC->Gspos("UL11",1,"UAIO",xpos,ypos,zpos,0,"ONLY");

}

//_____________________________________________________________________________
Bool_t AliTRDgeometry::Local2Global(Int_t idet, Float_t *local, Float_t *global)
{
  //
  // Converts local pad-coordinates (row,col,time) into 
  // global ALICE reference frame coordinates (x,y,z)
  //

  Int_t        icham     = GetChamber(idet);    // Chamber info (0-4)
  Int_t        isect     = GetSector(idet);     // Sector info  (0-17)
  Int_t        iplan     = GetPlane(idet);      // Plane info   (0-5)

  Float_t      padRow    = local[0];            // Pad Row position
  Float_t      padCol    = local[1];            // Pad Column position
  Float_t      timeSlice = local[2];            // Time "position"

  Float_t      row0      = GetRow0(iplan,icham,isect);
  Float_t      col0      = GetCol0(iplan);
  Float_t      time0     = GetTime0(iplan);

  Float_t      rot[3];

  // calculate (x,y) position in rotated chamber
  rot[1] = col0  + padCol    * fColPadSize;
  rot[0] = time0 + timeSlice * fTimeBinSize;
  // calculate z-position:
  rot[2] = row0  + padRow    * fRowPadSize;

  // Rotate back to original position
  return RotateBack(idet,rot,global);

}
 
//_____________________________________________________________________________
Bool_t AliTRDgeometry::Local2Global(Int_t iplan, Int_t icham, Int_t isect
                                  , Float_t *local, Float_t *global)
{
  //
  // Converts local pad-coordinates (row,col,time) into 
  // global ALICE reference frame coordinates (x,y,z)
  //

  Int_t        idet      = GetDetector(iplan,icham,isect); // Detector number

  Float_t      padRow    = local[0];                       // Pad Row position
  Float_t      padCol    = local[1];                       // Pad Column position
  Float_t      timeSlice = local[2];                       // Time "position"

  Float_t      row0      = GetRow0(iplan,icham,isect);
  Float_t      col0      = GetCol0(iplan);
  Float_t      time0     = GetTime0(iplan);

  Float_t      rot[3];

  // calculate (x,y,z) position in rotated chamber
  rot[1] = col0  + padCol    * fColPadSize;
  rot[0] = time0 + timeSlice * fTimeBinSize;
  rot[2] = row0  + padRow    * fRowPadSize;

  // Rotate back to original position
  return RotateBack(idet,rot,global);

}

//_____________________________________________________________________________
Bool_t AliTRDgeometry::Rotate(Int_t d, Float_t *pos, Float_t *rot)
{
  //
  // Rotates all chambers in the position of sector 0 and transforms
  // the coordinates in the ALICE restframe <pos> into the 
  // corresponding local frame <rot>.
  //

  Int_t   sector = GetSector(d);

  Float_t phi    = -2.0 * kPI /  (Float_t) kNsect * ((Float_t) sector + 0.5);

  rot[0] =  pos[0] * TMath::Cos(phi) + pos[1] * TMath::Sin(phi);
  rot[1] = -pos[0] * TMath::Sin(phi) + pos[1] * TMath::Cos(phi);
  rot[2] =  pos[2];

  return kTRUE;

}

//_____________________________________________________________________________
Bool_t AliTRDgeometry::RotateBack(Int_t d, Float_t *rot, Float_t *pos)
{
  //
  // Rotates a chambers from the position of sector 0 into its
  // original position and transforms the corresponding local frame 
  // coordinates <rot> into the coordinates of the ALICE restframe <pos>.
  //

  Int_t   sector = GetSector(d);

  Float_t phi    =  2.0 * kPI /  (Float_t) kNsect * ((Float_t) sector + 0.5);

  rot[0] =  pos[0] * TMath::Cos(phi) + pos[1] * TMath::Sin(phi);
  rot[1] = -pos[0] * TMath::Sin(phi) + pos[1] * TMath::Cos(phi);
  rot[2] =  pos[2];

  return kTRUE;

}

//_____________________________________________________________________________
Int_t AliTRDgeometry::GetDetector(Int_t p, Int_t c, Int_t s)
{
  //
  // Convert plane / chamber / sector into detector number
  //

  return (p + c * kNplan + s * kNplan * kNcham);

}

//_____________________________________________________________________________
Int_t AliTRDgeometry::GetPlane(Int_t d)
{
  //
  // Reconstruct the plane number from the detector number
  //

  return ((Int_t) (d % kNplan));

}

//_____________________________________________________________________________
Int_t AliTRDgeometry::GetChamber(Int_t d)
{
  //
  // Reconstruct the chamber number from the detector number
  //

  return ((Int_t) (d % (kNplan * kNcham)) / kNplan);

}

//_____________________________________________________________________________
Int_t AliTRDgeometry::GetSector(Int_t d)
{
  //
  // Reconstruct the sector number from the detector number
  //

  return ((Int_t) (d / (kNplan * kNcham)));

}

//_____________________________________________________________________________
void AliTRDgeometry::GetGlobal(const AliRecPoint *p, TVector3 &pos, TMatrix &mat)
{
  // 
  // Returns the global coordinate and error matrix of a AliTRDrecPoint
  //

  GetGlobal(p,pos);

}

//_____________________________________________________________________________
void AliTRDgeometry::GetGlobal(const AliRecPoint *p, TVector3 &pos)
{
  // 
  // Returns the global coordinate and error matrix of a AliTRDrecPoint
  //

  Int_t detector = ((AliTRDrecPoint *) p)->GetDetector();

  Float_t global[3];
  Float_t local[3];
  local[0] = pos.X();
  local[1] = pos.Y();
  local[2] = pos.Z();

  if (Local2Global(detector,local,global)) {
    pos.SetX(global[0]);
    pos.SetY(global[1]);
    pos.SetZ(global[2]);
  }
  else {
    pos.SetX(0.0);
    pos.SetY(0.0);
    pos.SetZ(0.0);
  }

}
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$
	18	*/
	19
	20	///////////////////////////////////////////////////////////////////////////////
	21	// //
	22	// TRD geometry class //
	23	// //
	24	///////////////////////////////////////////////////////////////////////////////
	25
	26	#include "AliTRDgeometry.h"
	27	#include "AliTRDrecPoint.h"
	28
	29	ClassImp(AliTRDgeometry)
	30
	31	//_____________________________________________________________________________
	32	AliTRDgeometry::AliTRDgeometry():AliGeometry()
	33	{
	34	//
	35	// AliTRDgeometry default constructor
	36	//
	37
	38	Init();
	39
	40	}
	41
	42	//_____________________________________________________________________________
	43	AliTRDgeometry::~AliTRDgeometry()
	44	{
	45
	46	}
	47
	48	//_____________________________________________________________________________
	49	void AliTRDgeometry::Init()
	50	{
	51	//
	52	// Initializes the geometry parameter
	53	//
	54
	55	Int_t iplan;
	56
	57	// The width of the chambers
	58	fCwidth[0] = 99.6;
	59	fCwidth[1] = 104.1;
	60	fCwidth[2] = 108.5;
	61	fCwidth[3] = 112.9;
	62	fCwidth[4] = 117.4;
	63	fCwidth[5] = 121.8;
	64
65	// The default pad dimensions
66	fRowPadSize = 4.5;
67	fColPadSize = 1.0;
68	fTimeBinSize = 0.1;
69
70	// The maximum number of pads
71	// and the position of pad 0,0,0
72	//
73	// chambers seen from the top:
74	// +----------------------------+
75	// \| \|
76	// \| \| ^
77	// \| \| rphi\|
78	// \| \| \|
79	// \|0 \| \|
80	// +----------------------------+ +------>
81	// z
82	// chambers seen from the side: ^
83	// +----------------------------+ time\|
84	// \| \| \|
85	// \|0 \| \|
86	// +----------------------------+ +------>
87	// z
88	//
89
90	// The pad column (rphi-direction)
91	for (iplan = 0; iplan < kNplan; iplan++) {
92	fColMax[iplan] = 1 + TMath::Nint((fCwidth[iplan] - 2. * kCcthick)
93	/ fColPadSize - 0.5);
94	fCol0[iplan] = -fCwidth[iplan]/2. + kCcthick;
95	}
96
97	// The time bucket
98	fTimeMax = 1 + TMath::Nint(kDrThick / fTimeBinSize - 0.5);
99	for (iplan = 0; iplan < kNplan; iplan++) {
100	fTime0[iplan] = kRmin + kCcframe/2. + kDrZpos - 0.5 * kDrThick
101	+ iplan * (kCheight + kCspace);
102	}
103
104	}
105
106	//_____________________________________________________________________________
107	void AliTRDgeometry::CreateGeometry(Int_t *idtmed)
108	{
109	//
110	// Create the TRD geometry
111	//
112	// Author: Christoph Blume (C.Blume@gsi.de) 20/07/99
113	//
114	// The volumes:
115	// TRD1-3 (Air) --- The TRD mother volumes for one sector.
116	// To be placed into the spaceframe.
117	//
118	// UAFI(/M/O) (Al) --- The aluminum frame of the inner(/middle/outer) chambers (readout)
119	// UCFI(/M/O) (C) --- The carbon frame of the inner(/middle/outer) chambers
120	// (driftchamber + radiator)
121	// UAII(/M/O) (Air) --- The inner part of the readout of the inner(/middle/outer) chambers
122	// UFII(/M/O) (Air) --- The inner part of the chamner and radiator of the
123	// inner(/middle/outer) chambers
124	//
125	// The material layers in one chamber:
126	// UL01 (G10) --- The gas seal of the radiator
127	// UL02 (CO2) --- The gas in the radiator
128	// UL03 (PE) --- The foil stack
129	// UL04 (Mylar) --- Entrance window to the driftvolume and HV-cathode
130	// UL05 (Xe) --- The driftvolume
131	// UL06 (Xe) --- The amplification region
132	//
133	// UL07 (Cu) --- The pad plane
134	// UL08 (G10) --- The Nomex honeycomb support structure
135	// UL09 (Cu) --- FEE and signal lines
136	// UL10 (PE) --- The cooling devices
137	// UL11 (Water) --- The cooling water
138
139	const Int_t npar_cha = 3;
140
141	Float_t par_dum[3];
142	Float_t par_cha[npar_cha];
143
144	Float_t xpos, ypos, zpos;
145
146	// The aluminum frames - readout + electronics (Al)
147	// The inner chambers
148	gMC->Gsvolu("UAFI","BOX ",idtmed[1301-1],par_dum,0);
149	// The middle chambers
150	gMC->Gsvolu("UAFM","BOX ",idtmed[1301-1],par_dum,0);
151	// The outer chambers
152	gMC->Gsvolu("UAFO","BOX ",idtmed[1301-1],par_dum,0);
153
154	// The inner part of the aluminum frames (Air)
155	// The inner chambers
156	gMC->Gsvolu("UAII","BOX ",idtmed[1302-1],par_dum,0);
157	// The middle chambers
158	gMC->Gsvolu("UAIM","BOX ",idtmed[1302-1],par_dum,0);
159	// The outer chambers
160	gMC->Gsvolu("UAIO","BOX ",idtmed[1302-1],par_dum,0);
161
162	// The carbon frames - radiator + driftchamber (C)
163	// The inner chambers
164	gMC->Gsvolu("UCFI","BOX ",idtmed[1307-1],par_dum,0);
165	// The middle chambers
166	gMC->Gsvolu("UCFM","BOX ",idtmed[1307-1],par_dum,0);
167	// The outer chambers
168	gMC->Gsvolu("UCFO","BOX ",idtmed[1307-1],par_dum,0);
169
170	// The inner part of the carbon frames (Air)
171	// The inner chambers
172	gMC->Gsvolu("UCII","BOX ",idtmed[1302-1],par_dum,0);
173	// The middle chambers
174	gMC->Gsvolu("UCIM","BOX ",idtmed[1302-1],par_dum,0);
175	// The outer chambers
176	gMC->Gsvolu("UCIO","BOX ",idtmed[1302-1],par_dum,0);
177
178	// The material layers inside the chambers
179	par_cha[0] = -1.;
180	par_cha[1] = -1.;
181	// G10 layer (radiator seal)
182	par_cha[2] = kSeThick/2;
183	gMC->Gsvolu("UL01","BOX ",idtmed[1313-1],par_cha,npar_cha);
184	// CO2 layer (radiator)
185	par_cha[2] = kRaThick/2;
186	gMC->Gsvolu("UL02","BOX ",idtmed[1312-1],par_cha,npar_cha);
187	// PE layer (radiator)
188	par_cha[2] = kPeThick/2;
189	gMC->Gsvolu("UL03","BOX ",idtmed[1303-1],par_cha,npar_cha);
190	// Mylar layer (entrance window + HV cathode)
191	par_cha[2] = kMyThick/2;
192	gMC->Gsvolu("UL04","BOX ",idtmed[1308-1],par_cha,npar_cha);
193	// Xe/Isobutane layer (drift volume, sensitive)
194	par_cha[2] = kDrThick/2.;
195	gMC->Gsvolu("UL05","BOX ",idtmed[1309-1],par_cha,npar_cha);
196	// Xe/Isobutane layer (amplification volume, not sensitive)
197	par_cha[2] = kAmThick/2.;
198	gMC->Gsvolu("UL06","BOX ",idtmed[1309-1],par_cha,npar_cha);
199
200	// Cu layer (pad plane)
201	par_cha[2] = kCuThick/2;
202	gMC->Gsvolu("UL07","BOX ",idtmed[1305-1],par_cha,npar_cha);
203	// G10 layer (support structure)
204	par_cha[2] = kSuThick/2;
205	gMC->Gsvolu("UL08","BOX ",idtmed[1313-1],par_cha,npar_cha);
206	// Cu layer (FEE + signal lines)
207	par_cha[2] = kFeThick/2;
208	gMC->Gsvolu("UL09","BOX ",idtmed[1305-1],par_cha,npar_cha);
209	// PE layer (cooling devices)
210	par_cha[2] = kCoThick/2;
211	gMC->Gsvolu("UL10","BOX ",idtmed[1303-1],par_cha,npar_cha);
212	// Water layer (cooling)
213	par_cha[2] = kWaThick/2;
214	gMC->Gsvolu("UL11","BOX ",idtmed[1314-1],par_cha,npar_cha);
215
216	// Position the layers in the chambers
217	xpos = 0;
218	ypos = 0;
219
220	// G10 layer (radiator seal)
221	zpos = kSeZpos;
222	gMC->Gspos("UL01",1,"UCII",xpos,ypos,zpos,0,"ONLY");
223	gMC->Gspos("UL01",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
224	gMC->Gspos("UL01",3,"UCIO",xpos,ypos,zpos,0,"ONLY");
225	// CO2 layer (radiator)
226	zpos = kRaZpos;
227	gMC->Gspos("UL02",1,"UCII",xpos,ypos,zpos,0,"ONLY");
228	gMC->Gspos("UL02",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
229	gMC->Gspos("UL02",3,"UCIO",xpos,ypos,zpos,0,"ONLY");
230	// PE layer (radiator)
231	zpos = 0;
232	gMC->Gspos("UL03",1,"UL02",xpos,ypos,zpos,0,"ONLY");
233	// Mylar layer (entrance window + HV cathode)
234	zpos = kMyZpos;
235	gMC->Gspos("UL04",1,"UCII",xpos,ypos,zpos,0,"ONLY");
236	gMC->Gspos("UL04",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
237	gMC->Gspos("UL04",3,"UCIO",xpos,ypos,zpos,0,"ONLY");
238	// Xe/Isobutane layer (drift volume)
239	zpos = kDrZpos;
240	gMC->Gspos("UL05",1,"UCII",xpos,ypos,zpos,0,"ONLY");
241	gMC->Gspos("UL05",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
242	gMC->Gspos("UL05",3,"UCIO",xpos,ypos,zpos,0,"ONLY");
243	// Xe/Isobutane layer (amplification volume)
244	zpos = kAmZpos;
245	gMC->Gspos("UL06",1,"UCII",xpos,ypos,zpos,0,"ONLY");
246	gMC->Gspos("UL06",2,"UCIM",xpos,ypos,zpos,0,"ONLY");
247	gMC->Gspos("UL06",3,"UCIO",xpos,ypos,zpos,0,"ONLY");
248
249	// Cu layer (pad plane)
250	zpos = kCuZpos;
251	gMC->Gspos("UL07",1,"UAII",xpos,ypos,zpos,0,"ONLY");
252	gMC->Gspos("UL07",2,"UAIM",xpos,ypos,zpos,0,"ONLY");
253	gMC->Gspos("UL07",3,"UAIO",xpos,ypos,zpos,0,"ONLY");
254	// G10 layer (support structure)
255	zpos = kSuZpos;
256	gMC->Gspos("UL08",1,"UAII",xpos,ypos,zpos,0,"ONLY");
257	gMC->Gspos("UL08",2,"UAIM",xpos,ypos,zpos,0,"ONLY");
258	gMC->Gspos("UL08",3,"UAIO",xpos,ypos,zpos,0,"ONLY");
259	// Cu layer (FEE + signal lines)
260	zpos = kFeZpos;
261	gMC->Gspos("UL09",1,"UAII",xpos,ypos,zpos,0,"ONLY");
262	gMC->Gspos("UL09",2,"UAIM",xpos,ypos,zpos,0,"ONLY");
263	gMC->Gspos("UL09",3,"UAIO",xpos,ypos,zpos,0,"ONLY");
264	// PE layer (cooling devices)
265	zpos = kCoZpos;
266	gMC->Gspos("UL10",1,"UAII",xpos,ypos,zpos,0,"ONLY");
267	gMC->Gspos("UL10",2,"UAIM",xpos,ypos,zpos,0,"ONLY");
268	gMC->Gspos("UL10",3,"UAIO",xpos,ypos,zpos,0,"ONLY");
269	// Water layer (cooling)
270	zpos = kWaZpos;
271	gMC->Gspos("UL11",1,"UAII",xpos,ypos,zpos,0,"ONLY");
272	gMC->Gspos("UL11",1,"UAIM",xpos,ypos,zpos,0,"ONLY");
273	gMC->Gspos("UL11",1,"UAIO",xpos,ypos,zpos,0,"ONLY");
274
275	}
276
277	//_____________________________________________________________________________
278	Bool_t AliTRDgeometry::Local2Global(Int_t idet, Float_t local, Float_t global)
279	{
280	//
281	// Converts local pad-coordinates (row,col,time) into
282	// global ALICE reference frame coordinates (x,y,z)
283	//
284
285	Int_t icham = GetChamber(idet); // Chamber info (0-4)
286	Int_t isect = GetSector(idet); // Sector info (0-17)
287	Int_t iplan = GetPlane(idet); // Plane info (0-5)
288
289	Float_t padRow = local[0]; // Pad Row position
290	Float_t padCol = local[1]; // Pad Column position
291	Float_t timeSlice = local[2]; // Time "position"
292
293	Float_t row0 = GetRow0(iplan,icham,isect);
294	Float_t col0 = GetCol0(iplan);
295	Float_t time0 = GetTime0(iplan);
296
297	Float_t rot[3];
298
299	// calculate (x,y) position in rotated chamber
300	rot[1] = col0 + padCol * fColPadSize;
301	rot[0] = time0 + timeSlice * fTimeBinSize;
302	// calculate z-position:
303	rot[2] = row0 + padRow * fRowPadSize;
304
305	// Rotate back to original position
306	return RotateBack(idet,rot,global);
307
308	}
309
310	//_____________________________________________________________________________
311	Bool_t AliTRDgeometry::Local2Global(Int_t iplan, Int_t icham, Int_t isect
312	, Float_t local, Float_t global)
313	{
314	//
315	// Converts local pad-coordinates (row,col,time) into
316	// global ALICE reference frame coordinates (x,y,z)
317	//
318
319	Int_t idet = GetDetector(iplan,icham,isect); // Detector number
320
321	Float_t padRow = local[0]; // Pad Row position
322	Float_t padCol = local[1]; // Pad Column position
323	Float_t timeSlice = local[2]; // Time "position"
324
325	Float_t row0 = GetRow0(iplan,icham,isect);
326	Float_t col0 = GetCol0(iplan);
327	Float_t time0 = GetTime0(iplan);
328
329	Float_t rot[3];
330
331	// calculate (x,y,z) position in rotated chamber
332	rot[1] = col0 + padCol * fColPadSize;
333	rot[0] = time0 + timeSlice * fTimeBinSize;
334	rot[2] = row0 + padRow * fRowPadSize;
335
336	// Rotate back to original position
337	return RotateBack(idet,rot,global);
338
339	}
340
341	//_____________________________________________________________________________
342	Bool_t AliTRDgeometry::Rotate(Int_t d, Float_t pos, Float_t rot)
343	{
344	//
345	// Rotates all chambers in the position of sector 0 and transforms
346	// the coordinates in the ALICE restframe <pos> into the
347	// corresponding local frame <rot>.
348	//
349
350	Int_t sector = GetSector(d);
351
352	Float_t phi = -2.0 * kPI / (Float_t) kNsect * ((Float_t) sector + 0.5);
353
354	rot[0] = pos[0] * TMath::Cos(phi) + pos[1] * TMath::Sin(phi);
355	rot[1] = -pos[0] * TMath::Sin(phi) + pos[1] * TMath::Cos(phi);
356	rot[2] = pos[2];
357
358	return kTRUE;
359
360	}
361
362	//_____________________________________________________________________________
363	Bool_t AliTRDgeometry::RotateBack(Int_t d, Float_t rot, Float_t pos)
364	{
365	//
366	// Rotates a chambers from the position of sector 0 into its
367	// original position and transforms the corresponding local frame
368	// coordinates <rot> into the coordinates of the ALICE restframe <pos>.
369	//
370
371	Int_t sector = GetSector(d);
372
373	Float_t phi = 2.0 * kPI / (Float_t) kNsect * ((Float_t) sector + 0.5);
374
375	rot[0] = pos[0] * TMath::Cos(phi) + pos[1] * TMath::Sin(phi);
376	rot[1] = -pos[0] * TMath::Sin(phi) + pos[1] * TMath::Cos(phi);
377	rot[2] = pos[2];
378
379	return kTRUE;
380
381	}
382
383	//_____________________________________________________________________________
384	Int_t AliTRDgeometry::GetDetector(Int_t p, Int_t c, Int_t s)
385	{
386	//
387	// Convert plane / chamber / sector into detector number
388	//
389
390	return (p + c * kNplan + s * kNplan * kNcham);
391
392	}
393
394	//_____________________________________________________________________________
395	Int_t AliTRDgeometry::GetPlane(Int_t d)
396	{
397	//
398	// Reconstruct the plane number from the detector number
399	//
400
401	return ((Int_t) (d % kNplan));
402
403	}
404
405	//_____________________________________________________________________________
406	Int_t AliTRDgeometry::GetChamber(Int_t d)
407	{
408	//
409	// Reconstruct the chamber number from the detector number
410	//
411
412	return ((Int_t) (d % (kNplan * kNcham)) / kNplan);
413
414	}
415
416	//_____________________________________________________________________________
417	Int_t AliTRDgeometry::GetSector(Int_t d)
418	{
419	//
420	// Reconstruct the sector number from the detector number
421	//
422
423	return ((Int_t) (d / (kNplan * kNcham)));
424
425	}
426
427	//_____________________________________________________________________________
428	void AliTRDgeometry::GetGlobal(const AliRecPoint *p, TVector3 &pos, TMatrix &mat)
429	{
430	//
431	// Returns the global coordinate and error matrix of a AliTRDrecPoint
432	//
433
434	GetGlobal(p,pos);
435
436	}
437
438	//_____________________________________________________________________________
439	void AliTRDgeometry::GetGlobal(const AliRecPoint *p, TVector3 &pos)
440	{
441	//
442	// Returns the global coordinate and error matrix of a AliTRDrecPoint
443	//
444
445	Int_t detector = ((AliTRDrecPoint *) p)->GetDetector();
446
447	Float_t global[3];
448	Float_t local[3];
449	local[0] = pos.X();
450	local[1] = pos.Y();
451	local[2] = pos.Z();
452
453	if (Local2Global(detector,local,global)) {
454	pos.SetX(global[0]);
455	pos.SetY(global[1]);
456	pos.SetZ(global[2]);
457	}
458	else {
459	pos.SetX(0.0);
460	pos.SetY(0.0);
461	pos.SetZ(0.0);
462	}
463
464	}