[u/mrichter/AliRoot.git] / MUON / AliMUONSt1GeometryBuilder.cxx

// $Id$
//
// Class AliMUONSt1GeometryBuilder
// -------------------------------
// MUON Station1 coarse geometry construction class.
// Extracted from AliMUONv1
// by Ivana Hrivnacova, IPN Orsay
// Included in AliRoot 2004/01/23

#include <TVirtualMC.h>
#include <TGeoMatrix.h>

#include "AliMUONSt1GeometryBuilder.h"
#include "AliMUON.h"
#include "AliMUONChamber.h"
#include "AliMUONChamberGeometry.h"

ClassImp(AliMUONSt1GeometryBuilder)

//______________________________________________________________________________
AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder(AliMUON* muon)
 : AliMUONVGeometryBuilder(&muon->Chamber(0), &muon->Chamber(1)),
   fMUON(muon)
{
// Standard constructor

}

//______________________________________________________________________________
AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder()
 : AliMUONVGeometryBuilder(),
   fMUON(0)
{
// Default constructor
}


//______________________________________________________________________________
AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder(const AliMUONSt1GeometryBuilder& rhs)
  : AliMUONVGeometryBuilder(rhs)
{
  Fatal("Copy constructor", 
        "Copy constructor is not implemented.");
}

//______________________________________________________________________________
AliMUONSt1GeometryBuilder::~AliMUONSt1GeometryBuilder() {
//
}

//______________________________________________________________________________
AliMUONSt1GeometryBuilder& 
AliMUONSt1GeometryBuilder::operator = (const AliMUONSt1GeometryBuilder& rhs) 
{
  // check assignement to self
  if (this == &rhs) return *this;

  Fatal("operator=", 
        "Assignment operator is not implemented.");
    
  return *this;  
}

//
// public methods
//

//______________________________________________________________________________
void AliMUONSt1GeometryBuilder::CreateGeometry()
{
// From AliMUONv1::CreateGeometry()
// ---

//********************************************************************
//                            Station 1                             **
//********************************************************************
//  CONCENTRIC
     // indices 1 and 2 for first and second chambers in the station
     // iChamber (first chamber) kept for other quanties than Z,
     // assumed to be the same in both chambers

     // Get tracking medias Ids     
     Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
     Int_t idAir= idtmed[1100]; // medium 1
     Int_t idAlu1=idtmed[1103]; // medium 4
     Int_t idAlu2=idtmed[1104]; // medium 5
     Int_t idGas=idtmed[1108];  // medium 9 = Ar-CO2 gas (80%+20%)
     Bool_t frameCrosses=kTRUE;     

     // Rotation matrices in the x-y plane  
     // phi=   0 deg
     Int_t irot1;
     fMUON->AliMatrix(irot1,  90.,   0., 90.,  90., 0., 0.);
     // phi=  90 deg
     Int_t irot2;
     fMUON->AliMatrix(irot2,  90.,  90., 90., 180., 0., 0.);

     AliMUONChamber* iChamber1 = GetChamber(0);
     AliMUONChamber* iChamber2 = GetChamber(1);
     AliMUONChamber* iChamber = iChamber1;

     // DGas decreased from standard one (0.5)
     iChamber->SetDGas(0.4); 
     iChamber2->SetDGas(0.4);

     // DAlu increased from standard one (3% of X0),
     // because more electronics with smaller pads
     iChamber->SetDAlu(3.5 * 8.9 / 100.); 
     iChamber2->SetDAlu(3.5 * 8.9 / 100.);

     // Half of the total thickness of frame crosses (including DAlu)
     // for each chamber in stations 1 and 2:
     // 3% of X0 of composite material,
     // but taken as Aluminium here, with same thickness in number of X0
     Float_t dframez = 3. * 8.9 / 100;
     Float_t zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
             // The same parameters are defined in builder for station 2 

     // Mother volume
     // Outer excess and inner recess for mother volume radius
     // with respect to ROuter and RInner
     Float_t dframep=.001; // Value for station 3 should be 6 ...
     // Width (RdPhi) of the frame crosses for stations 1 and 2 (cm)
     // Float_t dframep1=.001;
     Float_t dframep1 = 11.0;
     Float_t phi=2*TMath::Pi()/12/2;
             // The same parameters are defined in builder for station 2 

     Float_t tpar[3];
     Float_t dstation = (-iChamber2->Z()) - (-iChamber1->Z());
     tpar[0] = iChamber->RInner()-dframep; 
     tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi);
     tpar[2] = dstation/5;

     gMC->Gsvolu("S01M", "TUBE", idAir, tpar, 3);
     gMC->Gsvolu("S02M", "TUBE", idAir, tpar, 3);

     // CHANGED
     //gMC->Gspos("S01M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
     //gMC->Gspos("S02M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");     
     GetChamber(0)->GetGeometry()->AddEnvelope("S01M", false);
     GetChamber(1)->GetGeometry()->AddEnvelope("S02M", false);

// // Aluminium frames
// // Outer frames
//      pgpar[0] = 360/12/2;
//      pgpar[1] = 360.;
//      pgpar[2] = 12.;
//      pgpar[3] =   2;
//      pgpar[4] = -dframez/2;
//      pgpar[5] = iChamber->ROuter();
//      pgpar[6] = pgpar[5]+dframep1;
//      pgpar[7] = +dframez/2;
//      pgpar[8] = pgpar[5];
//      pgpar[9] = pgpar[6];
//      gMC->Gsvolu("S01O", "PGON", idAlu1, pgpar, 10);
//      gMC->Gsvolu("S02O", "PGON", idAlu1, pgpar, 10);
//      gMC->Gspos("S01O",1,"S01M", 0.,0.,-zfpos,  0,"ONLY");
//      gMC->Gspos("S01O",2,"S01M", 0.,0.,+zfpos,  0,"ONLY");
//      gMC->Gspos("S02O",1,"S02M", 0.,0.,-zfpos,  0,"ONLY");
//      gMC->Gspos("S02O",2,"S02M", 0.,0.,+zfpos,  0,"ONLY");
// //
// // Inner frame
//      tpar[0]= iChamber->RInner()-dframep1;
//      tpar[1]= iChamber->RInner();
//      tpar[2]= dframez/2;
//      gMC->Gsvolu("S01I", "TUBE", idAlu1, tpar, 3);
//      gMC->Gsvolu("S02I", "TUBE", idAlu1, tpar, 3);

//      gMC->Gspos("S01I",1,"S01M", 0.,0.,-zfpos,  0,"ONLY");
//      gMC->Gspos("S01I",2,"S01M", 0.,0.,+zfpos,  0,"ONLY");
//      gMC->Gspos("S02I",1,"S02M", 0.,0.,-zfpos,  0,"ONLY");
//      gMC->Gspos("S02I",2,"S02M", 0.,0.,+zfpos,  0,"ONLY");
//
// Frame Crosses
     if (frameCrosses) {
         // outside gas
         // security for inside mother volume
         Float_t bpar[3];
	 bpar[0] = (iChamber->ROuter() - iChamber->RInner())
	   * TMath::Cos(TMath::ASin(dframep1 /
				   (iChamber->ROuter() - iChamber->RInner())))
	   / 2.0;
	 bpar[1] = dframep1/2;
	 // total thickness will be (4 * bpar[2]) for each chamber,
	 // which has to be equal to (2 * dframez) - DAlu
	 bpar[2] = (2.0 * dframez - iChamber->DAlu()) / 4.0;
	 gMC->Gsvolu("S01B", "BOX", idAlu1, bpar, 3);
	 gMC->Gsvolu("S02B", "BOX", idAlu1, bpar, 3);
	 
	 gMC->Gspos("S01B",1,"S01M", +iChamber->RInner()+bpar[0] , 0,-zfpos, 
		    irot1,"ONLY");
	 gMC->Gspos("S01B",2,"S01M", -iChamber->RInner()-bpar[0] , 0,-zfpos, 
		    irot1,"ONLY");
	 gMC->Gspos("S01B",3,"S01M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, 
		    irot2,"ONLY");
	 gMC->Gspos("S01B",4,"S01M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, 
		    irot2,"ONLY");
	 gMC->Gspos("S01B",5,"S01M", +iChamber->RInner()+bpar[0] , 0,+zfpos, 
		    irot1,"ONLY");
	 gMC->Gspos("S01B",6,"S01M", -iChamber->RInner()-bpar[0] , 0,+zfpos, 
		    irot1,"ONLY");
	 gMC->Gspos("S01B",7,"S01M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, 
		    irot2,"ONLY");
	 gMC->Gspos("S01B",8,"S01M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, 
		    irot2,"ONLY");
	 
	 gMC->Gspos("S02B",1,"S02M", +iChamber->RInner()+bpar[0] , 0,-zfpos, 
		    irot1,"ONLY");
	 gMC->Gspos("S02B",2,"S02M", -iChamber->RInner()-bpar[0] , 0,-zfpos, 
		    irot1,"ONLY");
	 gMC->Gspos("S02B",3,"S02M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, 
		    irot2,"ONLY");
	 gMC->Gspos("S02B",4,"S02M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, 
		    irot2,"ONLY");
	 gMC->Gspos("S02B",5,"S02M", +iChamber->RInner()+bpar[0] , 0,+zfpos, 
		    irot1,"ONLY");
	 gMC->Gspos("S02B",6,"S02M", -iChamber->RInner()-bpar[0] , 0,+zfpos, 
		    irot1,"ONLY");
	 gMC->Gspos("S02B",7,"S02M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, 
		    irot2,"ONLY");
	 gMC->Gspos("S02B",8,"S02M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, 
		    irot2,"ONLY");
     }
//
//   Chamber Material represented by Alu sheet
     tpar[0]= iChamber->RInner();
     tpar[1]= iChamber->ROuter();
     tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2;
     gMC->Gsvolu("S01A", "TUBE",  idAlu2, tpar, 3);
     gMC->Gsvolu("S02A", "TUBE",idAlu2, tpar, 3);
     gMC->Gspos("S01A", 1, "S01M", 0., 0., 0.,  0, "ONLY");
     gMC->Gspos("S02A", 1, "S02M", 0., 0., 0.,  0, "ONLY");
//     
//   Sensitive volumes
     // tpar[2] = iChamber->DGas();
     tpar[2] = iChamber->DGas()/2;
     gMC->Gsvolu("S01G", "TUBE", idGas, tpar, 3);
     gMC->Gsvolu("S02G", "TUBE", idGas, tpar, 3);
     gMC->Gspos("S01G", 1, "S01A", 0., 0., 0.,  0, "ONLY");
     gMC->Gspos("S02G", 1, "S02A", 0., 0., 0.,  0, "ONLY");
//
// Frame Crosses to be placed inside gas
     // NONE: chambers are sensitive everywhere
//      if (frameCrosses) {

// 	 dr = (iChamber->ROuter() - iChamber->RInner());
// 	 bpar[0] = TMath::Sqrt(dr*dr-dframep1*dframep1/4)/2;
// 	 bpar[1] = dframep1/2;
// 	 bpar[2] = iChamber->DGas()/2;
// 	 gMC->Gsvolu("S01F", "BOX", idAlu1, bpar, 3);
// 	 gMC->Gsvolu("S02F", "BOX", idAlu1, bpar, 3);
	 
// 	 gMC->Gspos("S01F",1,"S01G", +iChamber->RInner()+bpar[0] , 0, 0, 
// 		    irot1,"ONLY");
// 	 gMC->Gspos("S01F",2,"S01G", -iChamber->RInner()-bpar[0] , 0, 0, 
// 		    irot1,"ONLY");
// 	 gMC->Gspos("S01F",3,"S01G", 0, +iChamber->RInner()+bpar[0] , 0, 
// 		    irot2,"ONLY");
// 	 gMC->Gspos("S01F",4,"S01G", 0, -iChamber->RInner()-bpar[0] , 0, 
// 		    irot2,"ONLY");
	 
// 	 gMC->Gspos("S02F",1,"S02G", +iChamber->RInner()+bpar[0] , 0, 0, 
// 		    irot1,"ONLY");
// 	 gMC->Gspos("S02F",2,"S02G", -iChamber->RInner()-bpar[0] , 0, 0, 
// 		    irot1,"ONLY");
// 	 gMC->Gspos("S02F",3,"S02G", 0, +iChamber->RInner()+bpar[0] , 0, 
// 		    irot2,"ONLY");
// 	 gMC->Gspos("S02F",4,"S02G", 0, -iChamber->RInner()-bpar[0] , 0, 
// 		    irot2,"ONLY");
//      }
}

//______________________________________________________________________________
void AliMUONSt1GeometryBuilder::SetTransformations()
{
// Defines the transformations for the station2 chambers.
// ---

  AliMUONChamber* iChamber1 = GetChamber(0);
  Double_t zpos1= - iChamber1->Z(); 
  iChamber1->GetGeometry()
    ->SetTranslation(TGeoTranslation(0., 0., zpos1));

  AliMUONChamber* iChamber2 = GetChamber(1);
  Double_t zpos2 = - iChamber2->Z(); 
  iChamber2->GetGeometry()
    ->SetTranslation(TGeoTranslation(0., 0., zpos2));
}

//______________________________________________________________________________
void AliMUONSt1GeometryBuilder::SetSensitiveVolumes()
{
// Defines the sensitive volumes for station1 chambers.
// ---

  GetChamber(0)->GetGeometry()->SetSensitiveVolume("S01G");
  GetChamber(1)->GetGeometry()->SetSensitiveVolume("S02G");
}
Commit	Line	Data
d1cd2474	1	// $Id$
	2	//
	3	// Class AliMUONSt1GeometryBuilder
	4	// -------------------------------
	5	// MUON Station1 coarse geometry construction class.
d1cd2474	6	// Extracted from AliMUONv1
d1cd2474	7	// by Ivana Hrivnacova, IPN Orsay
5f1df83a	8	// Included in AliRoot 2004/01/23
d1cd2474	9
	10	#include <TVirtualMC.h>
	11	#include <TGeoMatrix.h>
	12
	13	#include "AliMUONSt1GeometryBuilder.h"
	14	#include "AliMUON.h"
	15	#include "AliMUONChamber.h"
	16	#include "AliMUONChamberGeometry.h"
	17
	18	ClassImp(AliMUONSt1GeometryBuilder)
	19
	20	//______________________________________________________________________________
	21	AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder(AliMUON* muon)
	22	: AliMUONVGeometryBuilder(&muon->Chamber(0), &muon->Chamber(1)),
	23	fMUON(muon)
	24	{
	25	// Standard constructor
	26
	27	}
	28
	29	//______________________________________________________________________________
	30	AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder()
	31	: AliMUONVGeometryBuilder(),
	32	fMUON(0)
	33	{
	34	// Default constructor
	35	}
	36
	37
	38	//______________________________________________________________________________
	39	AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder(const AliMUONSt1GeometryBuilder& rhs)
	40	: AliMUONVGeometryBuilder(rhs)
	41	{
	42	Fatal("Copy constructor",
	43	"Copy constructor is not implemented.");
	44	}
	45
	46	//______________________________________________________________________________
	47	AliMUONSt1GeometryBuilder::~AliMUONSt1GeometryBuilder() {
	48	//
	49	}
	50
	51	//______________________________________________________________________________
	52	AliMUONSt1GeometryBuilder&
	53	AliMUONSt1GeometryBuilder::operator = (const AliMUONSt1GeometryBuilder& rhs)
	54	{
	55	// check assignement to self
	56	if (this == &rhs) return *this;
	57
	58	Fatal("operator=",
	59	"Assignment operator is not implemented.");
	60
	61	return *this;
	62	}
	63
	64	//
	65	// public methods
	66	//
	67
	68	//______________________________________________________________________________
	69	void AliMUONSt1GeometryBuilder::CreateGeometry()
	70	{
	71	// From AliMUONv1::CreateGeometry()
	72	// ---
73
74	//********************************************************************
75	// Station 1 **
76	//********************************************************************
77	// CONCENTRIC
78	// indices 1 and 2 for first and second chambers in the station
79	// iChamber (first chamber) kept for other quanties than Z,
80	// assumed to be the same in both chambers
81
82	// Get tracking medias Ids
83	Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
84	Int_t idAir= idtmed[1100]; // medium 1
85	Int_t idAlu1=idtmed[1103]; // medium 4
86	Int_t idAlu2=idtmed[1104]; // medium 5
87	Int_t idGas=idtmed[1108]; // medium 9 = Ar-CO2 gas (80%+20%)
88	Bool_t frameCrosses=kTRUE;
89
90	// Rotation matrices in the x-y plane
91	// phi= 0 deg
92	Int_t irot1;
93	fMUON->AliMatrix(irot1, 90., 0., 90., 90., 0., 0.);
94	// phi= 90 deg
95	Int_t irot2;
96	fMUON->AliMatrix(irot2, 90., 90., 90., 180., 0., 0.);
97
98	AliMUONChamber* iChamber1 = GetChamber(0);
99	AliMUONChamber* iChamber2 = GetChamber(1);
100	AliMUONChamber* iChamber = iChamber1;
101
102	// DGas decreased from standard one (0.5)
103	iChamber->SetDGas(0.4);
104	iChamber2->SetDGas(0.4);
105
106	// DAlu increased from standard one (3% of X0),
107	// because more electronics with smaller pads
108	iChamber->SetDAlu(3.5 * 8.9 / 100.);
109	iChamber2->SetDAlu(3.5 * 8.9 / 100.);
110
111	// Half of the total thickness of frame crosses (including DAlu)
112	// for each chamber in stations 1 and 2:
113	// 3% of X0 of composite material,
114	// but taken as Aluminium here, with same thickness in number of X0
115	Float_t dframez = 3. * 8.9 / 100;
116	Float_t zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
117	// The same parameters are defined in builder for station 2
118
119	// Mother volume
120	// Outer excess and inner recess for mother volume radius
121	// with respect to ROuter and RInner
122	Float_t dframep=.001; // Value for station 3 should be 6 ...
123	// Width (RdPhi) of the frame crosses for stations 1 and 2 (cm)
124	// Float_t dframep1=.001;
125	Float_t dframep1 = 11.0;
126	Float_t phi=2*TMath::Pi()/12/2;
127	// The same parameters are defined in builder for station 2
128
129	Float_t tpar[3];
130	Float_t dstation = (-iChamber2->Z()) - (-iChamber1->Z());
131	tpar[0] = iChamber->RInner()-dframep;
132	tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi);
133	tpar[2] = dstation/5;
134
135	gMC->Gsvolu("S01M", "TUBE", idAir, tpar, 3);
136	gMC->Gsvolu("S02M", "TUBE", idAir, tpar, 3);
137
138	// CHANGED
139	//gMC->Gspos("S01M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
140	//gMC->Gspos("S02M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");
141	GetChamber(0)->GetGeometry()->AddEnvelope("S01M", false);
142	GetChamber(1)->GetGeometry()->AddEnvelope("S02M", false);
143
144	// // Aluminium frames
145	// // Outer frames
146	// pgpar[0] = 360/12/2;
147	// pgpar[1] = 360.;
148	// pgpar[2] = 12.;
149	// pgpar[3] = 2;
150	// pgpar[4] = -dframez/2;
151	// pgpar[5] = iChamber->ROuter();
152	// pgpar[6] = pgpar[5]+dframep1;
153	// pgpar[7] = +dframez/2;
154	// pgpar[8] = pgpar[5];
155	// pgpar[9] = pgpar[6];
156	// gMC->Gsvolu("S01O", "PGON", idAlu1, pgpar, 10);
157	// gMC->Gsvolu("S02O", "PGON", idAlu1, pgpar, 10);
158	// gMC->Gspos("S01O",1,"S01M", 0.,0.,-zfpos, 0,"ONLY");
159	// gMC->Gspos("S01O",2,"S01M", 0.,0.,+zfpos, 0,"ONLY");
160	// gMC->Gspos("S02O",1,"S02M", 0.,0.,-zfpos, 0,"ONLY");
161	// gMC->Gspos("S02O",2,"S02M", 0.,0.,+zfpos, 0,"ONLY");
162	// //
163	// // Inner frame
164	// tpar[0]= iChamber->RInner()-dframep1;
165	// tpar[1]= iChamber->RInner();
166	// tpar[2]= dframez/2;
167	// gMC->Gsvolu("S01I", "TUBE", idAlu1, tpar, 3);
168	// gMC->Gsvolu("S02I", "TUBE", idAlu1, tpar, 3);
169
170	// gMC->Gspos("S01I",1,"S01M", 0.,0.,-zfpos, 0,"ONLY");
171	// gMC->Gspos("S01I",2,"S01M", 0.,0.,+zfpos, 0,"ONLY");
172	// gMC->Gspos("S02I",1,"S02M", 0.,0.,-zfpos, 0,"ONLY");
173	// gMC->Gspos("S02I",2,"S02M", 0.,0.,+zfpos, 0,"ONLY");
174	//
175	// Frame Crosses
176	if (frameCrosses) {
177	// outside gas
178	// security for inside mother volume
179	Float_t bpar[3];
180	bpar[0] = (iChamber->ROuter() - iChamber->RInner())
181	* TMath::Cos(TMath::ASin(dframep1 /
182	(iChamber->ROuter() - iChamber->RInner())))
183	/ 2.0;
184	bpar[1] = dframep1/2;
185	// total thickness will be (4 * bpar[2]) for each chamber,
186	// which has to be equal to (2 * dframez) - DAlu
187	bpar[2] = (2.0 * dframez - iChamber->DAlu()) / 4.0;
188	gMC->Gsvolu("S01B", "BOX", idAlu1, bpar, 3);
189	gMC->Gsvolu("S02B", "BOX", idAlu1, bpar, 3);
190
191	gMC->Gspos("S01B",1,"S01M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
192	irot1,"ONLY");
193	gMC->Gspos("S01B",2,"S01M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
194	irot1,"ONLY");
195	gMC->Gspos("S01B",3,"S01M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
196	irot2,"ONLY");
197	gMC->Gspos("S01B",4,"S01M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
198	irot2,"ONLY");
199	gMC->Gspos("S01B",5,"S01M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
200	irot1,"ONLY");
201	gMC->Gspos("S01B",6,"S01M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
202	irot1,"ONLY");
203	gMC->Gspos("S01B",7,"S01M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
204	irot2,"ONLY");
205	gMC->Gspos("S01B",8,"S01M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
206	irot2,"ONLY");
207
208	gMC->Gspos("S02B",1,"S02M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
209	irot1,"ONLY");
210	gMC->Gspos("S02B",2,"S02M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
211	irot1,"ONLY");
212	gMC->Gspos("S02B",3,"S02M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
213	irot2,"ONLY");
214	gMC->Gspos("S02B",4,"S02M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
215	irot2,"ONLY");
216	gMC->Gspos("S02B",5,"S02M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
217	irot1,"ONLY");
218	gMC->Gspos("S02B",6,"S02M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
219	irot1,"ONLY");
220	gMC->Gspos("S02B",7,"S02M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
221	irot2,"ONLY");
222	gMC->Gspos("S02B",8,"S02M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
223	irot2,"ONLY");
224	}
225	//
226	// Chamber Material represented by Alu sheet
227	tpar[0]= iChamber->RInner();
228	tpar[1]= iChamber->ROuter();
229	tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2;
230	gMC->Gsvolu("S01A", "TUBE", idAlu2, tpar, 3);
231	gMC->Gsvolu("S02A", "TUBE",idAlu2, tpar, 3);
232	gMC->Gspos("S01A", 1, "S01M", 0., 0., 0., 0, "ONLY");
233	gMC->Gspos("S02A", 1, "S02M", 0., 0., 0., 0, "ONLY");
234	//
235	// Sensitive volumes
236	// tpar[2] = iChamber->DGas();
237	tpar[2] = iChamber->DGas()/2;
238	gMC->Gsvolu("S01G", "TUBE", idGas, tpar, 3);
239	gMC->Gsvolu("S02G", "TUBE", idGas, tpar, 3);
240	gMC->Gspos("S01G", 1, "S01A", 0., 0., 0., 0, "ONLY");
241	gMC->Gspos("S02G", 1, "S02A", 0., 0., 0., 0, "ONLY");
242	//
243	// Frame Crosses to be placed inside gas
244	// NONE: chambers are sensitive everywhere
245	// if (frameCrosses) {
246
247	// dr = (iChamber->ROuter() - iChamber->RInner());
248	// bpar[0] = TMath::Sqrt(drdr-dframep1dframep1/4)/2;
249	// bpar[1] = dframep1/2;
250	// bpar[2] = iChamber->DGas()/2;
251	// gMC->Gsvolu("S01F", "BOX", idAlu1, bpar, 3);
252	// gMC->Gsvolu("S02F", "BOX", idAlu1, bpar, 3);
253
254	// gMC->Gspos("S01F",1,"S01G", +iChamber->RInner()+bpar[0] , 0, 0,
255	// irot1,"ONLY");
256	// gMC->Gspos("S01F",2,"S01G", -iChamber->RInner()-bpar[0] , 0, 0,
257	// irot1,"ONLY");
258	// gMC->Gspos("S01F",3,"S01G", 0, +iChamber->RInner()+bpar[0] , 0,
259	// irot2,"ONLY");
260	// gMC->Gspos("S01F",4,"S01G", 0, -iChamber->RInner()-bpar[0] , 0,
261	// irot2,"ONLY");
262
263	// gMC->Gspos("S02F",1,"S02G", +iChamber->RInner()+bpar[0] , 0, 0,
264	// irot1,"ONLY");
265	// gMC->Gspos("S02F",2,"S02G", -iChamber->RInner()-bpar[0] , 0, 0,
266	// irot1,"ONLY");
267	// gMC->Gspos("S02F",3,"S02G", 0, +iChamber->RInner()+bpar[0] , 0,
268	// irot2,"ONLY");
269	// gMC->Gspos("S02F",4,"S02G", 0, -iChamber->RInner()-bpar[0] , 0,
270	// irot2,"ONLY");
271	// }
272	}
273
274	//______________________________________________________________________________
275	void AliMUONSt1GeometryBuilder::SetTransformations()
276	{
277	// Defines the transformations for the station2 chambers.
278	// ---
279
280	AliMUONChamber* iChamber1 = GetChamber(0);
281	Double_t zpos1= - iChamber1->Z();
282	iChamber1->GetGeometry()
283	->SetTranslation(TGeoTranslation(0., 0., zpos1));
284
285	AliMUONChamber* iChamber2 = GetChamber(1);
286	Double_t zpos2 = - iChamber2->Z();
287	iChamber2->GetGeometry()
288	->SetTranslation(TGeoTranslation(0., 0., zpos2));
289	}
290
291	//______________________________________________________________________________
292	void AliMUONSt1GeometryBuilder::SetSensitiveVolumes()
293	{
294	// Defines the sensitive volumes for station1 chambers.
295	// ---
296
297	GetChamber(0)->GetGeometry()->SetSensitiveVolume("S01G");
298	GetChamber(1)->GetGeometry()->SetSensitiveVolume("S02G");
299	}