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

// $Id$
//
// Class AliMUONTriggerGeometryBuilder
// -----------------------------------
// MUON Trigger stations geometry construction class.
//
// Author: Philippe Crochette, LPC Clermont-Ferrand

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

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

ClassImp(AliMUONTriggerGeometryBuilder)

//______________________________________________________________________________
AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder(AliMUON* muon)
 : AliMUONVGeometryBuilder(&muon->Chamber(10), &muon->Chamber(11),&muon->Chamber(12),&muon->Chamber(13)),
   fMUON(muon)
{
// Standard constructor

}

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


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

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

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

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

//
// public methods
//

//______________________________________________________________________________
void AliMUONTriggerGeometryBuilder::CreateGeometry()
{
// From AliMUONv1::CreateGeometry()

    /* 
       zpos1 and zpos2 are the middle of the first and second
       planes of station 1 (+1m for second station):
       zpos1=(zpos1m+zpos1p)/2=(15999+16071)/2=16035 mm, thick/2=40 mm
       zpos2=(zpos2m+zpos2p)/2=(16169+16241)/2=16205 mm, thick/2=40 mm
       zposxm and zposxp= middles of gaz gaps within a detection plane
       rem: the total thickness accounts for 1 mm of al on both
       side of the RPCs (see zpos1 and zpos2)
    */
    
    Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
    Int_t idAir= idtmed[1100]; // medium 1
    Int_t idAlu1=idtmed[1103]; // medium 4
    Float_t tpar[3];
    Double_t dpar[3];    
    
// vertical gap between right and left chambers (kDXZERO*2=4cm)
    const Float_t kDXZERO=2.; 
// main distances for chamber definition in first plane/first station
    const Float_t kXMIN=34.;       
    const Float_t kXMED=51.;                                
    const Float_t kXMAX=272.; 
// kXMAX will become 255. in real life. segmentation to be updated accordingly
// (see fig.2-4 & 2-5 of Local Trigger Board PRR)
    const Float_t kYMIN=34.;                              
    const Float_t kYMAX=51.;                              
// inner/outer radius of flange between beam shield. and chambers (1/station)
    const Float_t kRMIN[2]={50.,50.};
    const Float_t kRMAX[2]={64.,68.};
// z position of the middle of the gas gap in mother vol 
    const Float_t kZm=-3.6;
    const Float_t kZp=+3.6;     

    AliMUONChamber *iChamber, *iChamber1;
    iChamber1 = GetChamber(10);
    Float_t zpos1=-iChamber1->Z(); 
    
// ratio of zpos1m/zpos1p and inverse for first plane
    Float_t zmp=(zpos1-3.6)/(zpos1+3.6);
    Float_t zpm=1./zmp;
    
    Int_t icount=0;  // chamber counter (0 1 2 3)
    
    for (Int_t istation=0; istation<2; istation++) { // loop on stations
	for (Int_t iplane=0; iplane<2; iplane++) { // loop on detection planes
	    
	    Int_t iVolNum=1; // counter Volume Number
	    icount = Int_t(iplane*TMath::Power(2,0))+
		Int_t(istation*TMath::Power(2,1));
	    
	    iChamber = GetChamber(10+icount);
	    Float_t zpos = - iChamber->Z();	     
	    
// Flange between beam shielding and RPC 
	    tpar[0]= kRMIN[istation];
	    tpar[1]= kRMAX[istation];
	    tpar[2]= 4.0;	    
	    char volFlange[5];
	    sprintf(volFlange,"SF%dA",icount+1);	 
	    gMC->Gsvolu(volFlange,"TUBE",idAlu1,tpar,3);     // Al
            // changed by ivana
	    //gMC->Gspos(volFlange,1,"ALIC",0.,0.,zpos,0,"MANY");
	    iChamber->GetGeometry()->AddEnvelope(volFlange, false, "MANY");
	    
// scaling factor
	    Float_t zRatio = zpos / zpos1;
	    
// envelopes (same size except line 5, all virtual)
	    char volEnv[18][5];
	    tpar[1] = kYMIN * zRatio; 
	    tpar[2] = 0.4;
	    Int_t i=0;    // counter
	    for (Int_t icolumn=0; icolumn<2; icolumn++) {
		for (Int_t iline=1; iline<10; iline++){
		    tpar[0] = (kXMAX/2.) * zRatio;
		    if (iline==5) tpar[0] = ((kXMAX-kXMED)/2.)*zRatio;
		    if (icolumn==0) 
			sprintf(volEnv[i],"S%dL%d",icount,iline);
		    else
			sprintf(volEnv[i],"S%dR%d",icount,iline);
		    gMC->Gsvolu(volEnv[i],"BOX",idAir,tpar,0); 
		    i++;
		}
	    }

// chamber prototype
	    tpar[0]= 0.;
	    tpar[1]= 0.;
	    tpar[2]= 0.;	    
	    char volAlu[5];     // Alu 
	    char volBak[5];     // Bakelite
	    char volGaz[5];     // Gas streamer	    
	    sprintf(volAlu,"SC%dA",icount+1);
	    sprintf(volBak,"SB%dA",icount+1);
	    sprintf(volGaz,"SG%dA",icount+1);
	    gMC->Gsvolu(volAlu,"BOX",idAlu1,tpar,0);         // Al
	    gMC->Gsvolu(volBak,"BOX",idtmed[1107],tpar,0);   // Bakelite
	    gMC->Gsvolu(volGaz,"BOX",idtmed[1106],tpar,0);   // Gas streamer
	    tpar[0] = -1.;
	    tpar[1] = -1.;
	    tpar[2] = 0.1;    
	    gMC->Gsposp(volGaz,1,volBak,0.,0.,0.,0,"ONLY",tpar,3);
	    tpar[2] = 0.3;
	    gMC->Gsposp(volBak,1,volAlu,0.,0.,0.,0,"ONLY",tpar,3);

// chamber type A
	    Float_t xEnv = (kDXZERO+kXMED+(kXMAX-kXMED)/2.)*zRatio;
	    Float_t yEnvM = 0.;	 // y low position of envelope in chamber
	    Float_t yEnvP = 0.;	 // y up position of envelope in chamber
	    Float_t yEnvPsave = 0.; // tmp data
	    Float_t yEnvMsave = 0.; // tmp data
	    Float_t xpos = 0.; // x position of RPC in envelope	    
	    Float_t ypos = 0.; // y position of RPC in envelope
	    dpar[2] = 0.4;	    
	    dpar[0] = ((kXMAX-kXMED)/2.)*zRatio;
	    dpar[1] = kYMIN * zRatio;

	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[4], true, TGeoTranslation(xEnv,yEnvM,kZm));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[13], true, TGeoTranslation(-xEnv,yEnvP,kZp));

	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[4],iVolNum++,3, dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[13],iVolNum++,3, dpar);	    

// chamber type B (plus envelope chambers B & C)   
	    xEnv = (kDXZERO+kXMAX/2.)*zRatio;
	    yEnvPsave = yEnvP;
	    yEnvMsave = yEnvM;
	    yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
	    yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
	    dpar[0] = ((kXMAX-kXMIN)/2.) * zRatio;
	    dpar[1] = ((kYMAX-kYMIN)/2.) * zRatio;
	    xpos = kXMIN/2. * zRatio;
	    ypos = (kYMIN - kYMIN/4.) * zRatio;

	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[3], true, TGeoTranslation( xEnv,-yEnvP,kZp));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[5], true, TGeoTranslation( xEnv, yEnvP,kZp));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[12], true, TGeoTranslation(-xEnv,-yEnvM,kZm));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[14], true, TGeoTranslation(-xEnv, yEnvM,kZm));

	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar);

// chamber type C (note: same Z than type B)
	    dpar[0] = (kXMAX/2)*zRatio;
	    dpar[1] = (kYMAX/2)*zRatio;
	    xpos = 0.;	    
	    ypos = ((kYMAX - kYMIN)/2.) * zRatio;

	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar);
    
// chamber type D, E and F (same size)
// D	    
	    yEnvPsave = yEnvP;
	    yEnvMsave = yEnvM;
	    yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
	    yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
	    dpar[0] = (kXMAX/2.)*zRatio;
	    dpar[1] =  kYMIN*zRatio;

	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[2], true, TGeoTranslation(xEnv, yEnvM,kZm));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[6], true, TGeoTranslation(-xEnv, yEnvP,kZp));	    
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[11], true, TGeoTranslation(xEnv,-yEnvM,kZm));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[15], true, TGeoTranslation(-xEnv,-yEnvP,kZp));

	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[2],iVolNum++,3, dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[6],iVolNum++,3, dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[11],iVolNum++,3, dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[15],iVolNum++,3, dpar);

// E
	    yEnvPsave = yEnvP;
	    yEnvMsave = yEnvM;
	    yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
	    yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;

	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[1], true, TGeoTranslation(xEnv, yEnvP,kZp));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[7], true, TGeoTranslation(-xEnv, yEnvM,kZm));	    
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[10], true, TGeoTranslation(xEnv,-yEnvP,kZp));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[16], true, TGeoTranslation(-xEnv,-yEnvM,kZm));

	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[1],iVolNum++,3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[7],iVolNum++,3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[10],iVolNum++,3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[16],iVolNum++,3,dpar);


// F
	    yEnvPsave = yEnvP;
	    yEnvMsave = yEnvM;
	    yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
	    yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;

	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[0], true, TGeoTranslation(xEnv, yEnvM,kZm));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[8], true, TGeoTranslation(-xEnv, yEnvP,kZp));	    
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[9], true, TGeoTranslation(xEnv,-yEnvM,kZm));
	    GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[17], true, TGeoTranslation(-xEnv,-yEnvP,kZp));

	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[0],iVolNum++,3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[8],iVolNum++,3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[9],iVolNum++,3,dpar);
	    GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[17],iVolNum++,3,dpar);

	} // end loop on detection planes
    } // end loop on stations    
}

//______________________________________________________________________________
void AliMUONTriggerGeometryBuilder::SetTransformations() 
{
// Defines the transformations for the trigger chambers.
// ---
    Double_t zpos1, zpos2;    
    AliMUONChamber *iChamber1, *iChamber2;

    iChamber1 = GetChamber(10);
    zpos1= - iChamber1->Z(); 
    iChamber1->GetGeometry()
	->SetTranslation(TGeoTranslation(0., 0., zpos1));
    
    iChamber2 = GetChamber(11);
    zpos2 = - iChamber2->Z(); 
    iChamber2->GetGeometry()
	->SetTranslation(TGeoTranslation(0., 0., zpos2));

    iChamber1 = GetChamber(12);
    zpos1 = - iChamber1->Z(); 
    iChamber1->GetGeometry()
	->SetTranslation(TGeoTranslation(0., 0., zpos1));
    
    iChamber2 = GetChamber(13);
    zpos2 = - iChamber2->Z(); 
    iChamber2->GetGeometry()
	->SetTranslation(TGeoTranslation(0., 0., zpos2));
}

//______________________________________________________________________________
void AliMUONTriggerGeometryBuilder::SetSensitiveVolumes()
{
// Defines the sensitive volumes for trigger station chambers.
// ---

  GetChamber(10)->GetGeometry()->SetSensitiveVolume("SG1A");
  GetChamber(11)->GetGeometry()->SetSensitiveVolume("SG2A");
  GetChamber(12)->GetGeometry()->SetSensitiveVolume("SG3A");
  GetChamber(13)->GetGeometry()->SetSensitiveVolume("SG4A");
}
Commit	Line	Data
d1cd2474	1	// $Id$
	2	//
	3	// Class AliMUONTriggerGeometryBuilder
	4	// -----------------------------------
	5	// MUON Trigger stations geometry construction class.
	6	//
	7	// Author: Philippe Crochette, LPC Clermont-Ferrand
	8
	9	#include <TVirtualMC.h>
	10	#include <TGeoMatrix.h>
	11
	12	#include "AliMUONTriggerGeometryBuilder.h"
	13	#include "AliMUON.h"
	14	#include "AliMUONChamber.h"
	15	#include "AliMUONChamberGeometry.h"
	16
	17	ClassImp(AliMUONTriggerGeometryBuilder)
	18
	19	//______________________________________________________________________________
	20	AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder(AliMUON* muon)
	21	: AliMUONVGeometryBuilder(&muon->Chamber(10), &muon->Chamber(11),&muon->Chamber(12),&muon->Chamber(13)),
	22	fMUON(muon)
	23	{
	24	// Standard constructor
	25
	26	}
	27
	28	//______________________________________________________________________________
	29	AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder()
	30	: AliMUONVGeometryBuilder(),
	31	fMUON(0)
	32	{
	33	// Default constructor
	34	}
	35
	36
	37	//______________________________________________________________________________
	38	AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder(const AliMUONTriggerGeometryBuilder& rhs)
	39	: AliMUONVGeometryBuilder(rhs)
	40	{
	41	Fatal("Copy constructor",
	42	"Copy constructor is not implemented.");
	43	}
	44
	45	//______________________________________________________________________________
	46	AliMUONTriggerGeometryBuilder::~AliMUONTriggerGeometryBuilder() {
	47	//
	48	}
	49
	50	//______________________________________________________________________________
	51	AliMUONTriggerGeometryBuilder&
	52	AliMUONTriggerGeometryBuilder::operator = (const AliMUONTriggerGeometryBuilder& rhs)
	53	{
	54	// check assignement to self
	55	if (this == &rhs) return *this;
	56
	57	Fatal("operator=",
	58	"Assignment operator is not implemented.");
	59
	60	return *this;
	61	}
	62
	63	//
	64	// public methods
65	//
66
67	//______________________________________________________________________________
68	void AliMUONTriggerGeometryBuilder::CreateGeometry()
69	{
70	// From AliMUONv1::CreateGeometry()
71
72	/*
73	zpos1 and zpos2 are the middle of the first and second
74	planes of station 1 (+1m for second station):
75	zpos1=(zpos1m+zpos1p)/2=(15999+16071)/2=16035 mm, thick/2=40 mm
76	zpos2=(zpos2m+zpos2p)/2=(16169+16241)/2=16205 mm, thick/2=40 mm
77	zposxm and zposxp= middles of gaz gaps within a detection plane
78	rem: the total thickness accounts for 1 mm of al on both
79	side of the RPCs (see zpos1 and zpos2)
80	*/
81
82	Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
83	Int_t idAir= idtmed[1100]; // medium 1
84	Int_t idAlu1=idtmed[1103]; // medium 4
85	Float_t tpar[3];
86	Double_t dpar[3];
87
88	// vertical gap between right and left chambers (kDXZERO*2=4cm)
89	const Float_t kDXZERO=2.;
90	// main distances for chamber definition in first plane/first station
91	const Float_t kXMIN=34.;
92	const Float_t kXMED=51.;
93	const Float_t kXMAX=272.;
94	// kXMAX will become 255. in real life. segmentation to be updated accordingly
95	// (see fig.2-4 & 2-5 of Local Trigger Board PRR)
96	const Float_t kYMIN=34.;
97	const Float_t kYMAX=51.;
98	// inner/outer radius of flange between beam shield. and chambers (1/station)
99	const Float_t kRMIN[2]={50.,50.};
100	const Float_t kRMAX[2]={64.,68.};
101	// z position of the middle of the gas gap in mother vol
102	const Float_t kZm=-3.6;
103	const Float_t kZp=+3.6;
104
105	AliMUONChamber iChamber, iChamber1;
106	iChamber1 = GetChamber(10);
107	Float_t zpos1=-iChamber1->Z();
108
109	// ratio of zpos1m/zpos1p and inverse for first plane
110	Float_t zmp=(zpos1-3.6)/(zpos1+3.6);
111	Float_t zpm=1./zmp;
112
113	Int_t icount=0; // chamber counter (0 1 2 3)
114
115	for (Int_t istation=0; istation<2; istation++) { // loop on stations
116	for (Int_t iplane=0; iplane<2; iplane++) { // loop on detection planes
117
118	Int_t iVolNum=1; // counter Volume Number
119	icount = Int_t(iplane*TMath::Power(2,0))+
120	Int_t(istation*TMath::Power(2,1));
121
122	iChamber = GetChamber(10+icount);
123	Float_t zpos = - iChamber->Z();
124
125	// Flange between beam shielding and RPC
126	tpar[0]= kRMIN[istation];
127	tpar[1]= kRMAX[istation];
128	tpar[2]= 4.0;
129	char volFlange[5];
130	sprintf(volFlange,"SF%dA",icount+1);
131	gMC->Gsvolu(volFlange,"TUBE",idAlu1,tpar,3); // Al
132	// changed by ivana
133	//gMC->Gspos(volFlange,1,"ALIC",0.,0.,zpos,0,"MANY");
6b82c1f0	134	iChamber->GetGeometry()->AddEnvelope(volFlange, false, "MANY");
d1cd2474	135
	136	// scaling factor
	137	Float_t zRatio = zpos / zpos1;
	138
	139	// envelopes (same size except line 5, all virtual)
	140	char volEnv[18][5];
	141	tpar[1] = kYMIN * zRatio;
	142	tpar[2] = 0.4;
	143	Int_t i=0; // counter
	144	for (Int_t icolumn=0; icolumn<2; icolumn++) {
	145	for (Int_t iline=1; iline<10; iline++){
	146	tpar[0] = (kXMAX/2.) * zRatio;
	147	if (iline==5) tpar[0] = ((kXMAX-kXMED)/2.)*zRatio;
	148	if (icolumn==0)
	149	sprintf(volEnv[i],"S%dL%d",icount,iline);
	150	else
	151	sprintf(volEnv[i],"S%dR%d",icount,iline);
	152	gMC->Gsvolu(volEnv[i],"BOX",idAir,tpar,0);
	153	i++;
	154	}
	155	}
	156
	157	// chamber prototype
	158	tpar[0]= 0.;
	159	tpar[1]= 0.;
	160	tpar[2]= 0.;
	161	char volAlu[5]; // Alu
	162	char volBak[5]; // Bakelite
	163	char volGaz[5]; // Gas streamer
	164	sprintf(volAlu,"SC%dA",icount+1);
	165	sprintf(volBak,"SB%dA",icount+1);
	166	sprintf(volGaz,"SG%dA",icount+1);
	167	gMC->Gsvolu(volAlu,"BOX",idAlu1,tpar,0); // Al
	168	gMC->Gsvolu(volBak,"BOX",idtmed[1107],tpar,0); // Bakelite
	169	gMC->Gsvolu(volGaz,"BOX",idtmed[1106],tpar,0); // Gas streamer
	170	tpar[0] = -1.;
	171	tpar[1] = -1.;
	172	tpar[2] = 0.1;
	173	gMC->Gsposp(volGaz,1,volBak,0.,0.,0.,0,"ONLY",tpar,3);
	174	tpar[2] = 0.3;
	175	gMC->Gsposp(volBak,1,volAlu,0.,0.,0.,0,"ONLY",tpar,3);
	176
	177	// chamber type A
	178	Float_t xEnv = (kDXZERO+kXMED+(kXMAX-kXMED)/2.)*zRatio;
	179	Float_t yEnvM = 0.; // y low position of envelope in chamber
	180	Float_t yEnvP = 0.; // y up position of envelope in chamber
	181	Float_t yEnvPsave = 0.; // tmp data
	182	Float_t yEnvMsave = 0.; // tmp data
	183	Float_t xpos = 0.; // x position of RPC in envelope
	184	Float_t ypos = 0.; // y position of RPC in envelope
	185	dpar[2] = 0.4;
	186	dpar[0] = ((kXMAX-kXMED)/2.)*zRatio;
	187	dpar[1] = kYMIN * zRatio;
	188
	189	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[4], true, TGeoTranslation(xEnv,yEnvM,kZm));
	190	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[13], true, TGeoTranslation(-xEnv,yEnvP,kZp));
	191
	192	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[4],iVolNum++,3, dpar);
	193	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[13],iVolNum++,3, dpar);
	194
	195	// chamber type B (plus envelope chambers B & C)
	196	xEnv = (kDXZERO+kXMAX/2.)*zRatio;
	197	yEnvPsave = yEnvP;
	198	yEnvMsave = yEnvM;
199	yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
200	yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
201	dpar[0] = ((kXMAX-kXMIN)/2.) * zRatio;
202	dpar[1] = ((kYMAX-kYMIN)/2.) * zRatio;
203	xpos = kXMIN/2. * zRatio;
204	ypos = (kYMIN - kYMIN/4.) * zRatio;
205
206	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[3], true, TGeoTranslation( xEnv,-yEnvP,kZp));
207	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[5], true, TGeoTranslation( xEnv, yEnvP,kZp));
208	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[12], true, TGeoTranslation(-xEnv,-yEnvM,kZm));
209	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[14], true, TGeoTranslation(-xEnv, yEnvM,kZm));
210
211	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
212	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
213	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar);
214	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar);
215
216	// chamber type C (note: same Z than type B)
217	dpar[0] = (kXMAX/2)*zRatio;
218	dpar[1] = (kYMAX/2)*zRatio;
219	xpos = 0.;
220	ypos = ((kYMAX - kYMIN)/2.) * zRatio;
221
222	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
223	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
224	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar);
225	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar);
226
227	// chamber type D, E and F (same size)
228	// D
229	yEnvPsave = yEnvP;
230	yEnvMsave = yEnvM;
231	yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
232	yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
233	dpar[0] = (kXMAX/2.)*zRatio;
234	dpar[1] = kYMIN*zRatio;
235
236	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[2], true, TGeoTranslation(xEnv, yEnvM,kZm));
237	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[6], true, TGeoTranslation(-xEnv, yEnvP,kZp));
238	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[11], true, TGeoTranslation(xEnv,-yEnvM,kZm));
239	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[15], true, TGeoTranslation(-xEnv,-yEnvP,kZp));
240
241	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[2],iVolNum++,3, dpar);
242	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[6],iVolNum++,3, dpar);
243	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[11],iVolNum++,3, dpar);
244	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[15],iVolNum++,3, dpar);
245
246	// E
247	yEnvPsave = yEnvP;
248	yEnvMsave = yEnvM;
249	yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
250	yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
251
252	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[1], true, TGeoTranslation(xEnv, yEnvP,kZp));
253	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[7], true, TGeoTranslation(-xEnv, yEnvM,kZm));
254	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[10], true, TGeoTranslation(xEnv,-yEnvP,kZp));
255	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[16], true, TGeoTranslation(-xEnv,-yEnvM,kZm));
256
257	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[1],iVolNum++,3,dpar);
258	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[7],iVolNum++,3,dpar);
259	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[10],iVolNum++,3,dpar);
260	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[16],iVolNum++,3,dpar);
261
262
263	// F
264	yEnvPsave = yEnvP;
265	yEnvMsave = yEnvM;
266	yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
267	yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
268
269	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[0], true, TGeoTranslation(xEnv, yEnvM,kZm));
270	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[8], true, TGeoTranslation(-xEnv, yEnvP,kZp));
271	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[9], true, TGeoTranslation(xEnv,-yEnvM,kZm));
272	GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[17], true, TGeoTranslation(-xEnv,-yEnvP,kZp));
273
274	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[0],iVolNum++,3,dpar);
275	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[8],iVolNum++,3,dpar);
276	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[9],iVolNum++,3,dpar);
277	GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[17],iVolNum++,3,dpar);
278
279	} // end loop on detection planes
280	} // end loop on stations
281	}
282
283	//______________________________________________________________________________
284	void AliMUONTriggerGeometryBuilder::SetTransformations()
285	{
286	// Defines the transformations for the trigger chambers.
287	// ---
288	Double_t zpos1, zpos2;
289	AliMUONChamber iChamber1, iChamber2;
290
291	iChamber1 = GetChamber(10);
292	zpos1= - iChamber1->Z();
293	iChamber1->GetGeometry()
294	->SetTranslation(TGeoTranslation(0., 0., zpos1));
295
296	iChamber2 = GetChamber(11);
297	zpos2 = - iChamber2->Z();
298	iChamber2->GetGeometry()
299	->SetTranslation(TGeoTranslation(0., 0., zpos2));
300
301	iChamber1 = GetChamber(12);
302	zpos1 = - iChamber1->Z();
303	iChamber1->GetGeometry()
304	->SetTranslation(TGeoTranslation(0., 0., zpos1));
305
306	iChamber2 = GetChamber(13);
307	zpos2 = - iChamber2->Z();
308	iChamber2->GetGeometry()
309	->SetTranslation(TGeoTranslation(0., 0., zpos2));
310	}
311
312	//______________________________________________________________________________
313	void AliMUONTriggerGeometryBuilder::SetSensitiveVolumes()
314	{
315	// Defines the sensitive volumes for trigger station chambers.
316	// ---
317
318	GetChamber(10)->GetGeometry()->SetSensitiveVolume("SG1A");
319	GetChamber(11)->GetGeometry()->SetSensitiveVolume("SG2A");
320	GetChamber(12)->GetGeometry()->SetSensitiveVolume("SG3A");
321	GetChamber(13)->GetGeometry()->SetSensitiveVolume("SG4A");
322	}
323