X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=MUON%2FAliMUONSt2GeometryBuilderV2.cxx;h=8264c822e0dbb4801085a9ebb26bc85d86207486;hp=bcf50ec5ef0314f48360a5e3d87caf1f6f9566f3;hb=71d5ee506ed0a3db4d9ef0db9fcdc842dfc24bbe;hpb=c0404a6cca997b27f75916bc052d754a59fea6a6 diff --git a/MUON/AliMUONSt2GeometryBuilderV2.cxx b/MUON/AliMUONSt2GeometryBuilderV2.cxx index bcf50ec5ef0..8264c822e0d 100644 --- a/MUON/AliMUONSt2GeometryBuilderV2.cxx +++ b/MUON/AliMUONSt2GeometryBuilderV2.cxx @@ -14,19 +14,20 @@ **************************************************************************/ // $Id$ -// + +//----------------------------------------------------------------------------- // Class AliMUONSt2GeometryBuilderV2 -// ------------------------------- +// --------------------------------- // MUON Station2 coarse geometry construction class. -//******************************************************************** // Author: SANJOY PAL ,Prof. SUKALYAN CHATTOPADHAYAY [SINP, KOLKATA] // & Dr.SHAKEEL AHMAD (AMU), INDIA -//******************************************************************** +//----------------------------------------------------------------------------- #include #include #include +#include #include "AliRun.h" #include "AliLog.h" @@ -39,14 +40,17 @@ #include "AliMUONConstants.h" #define PI 3.14159 + +/// \cond CLASSIMP ClassImp(AliMUONSt2GeometryBuilderV2) +/// \endcond //______________________________________________________________________________ AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2(AliMUON* muon) - : AliMUONVGeometryBuilder(2, 3), + : AliMUONVGeometryBuilder(2, 2), fMUON(muon) { -// Standard constructor +/// Standard constructor } @@ -55,36 +59,13 @@ AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2() : AliMUONVGeometryBuilder(), fMUON(0) { -// Default constructor -} - - -//______________________________________________________________________________ -AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2(const AliMUONSt2GeometryBuilderV2& rhs) - : AliMUONVGeometryBuilder(rhs) -{ -// Protected copy constructor - - AliFatal("Copy constructor is not implemented."); -} - -//______________________________________________________________________________ -AliMUONSt2GeometryBuilderV2::~AliMUONSt2GeometryBuilderV2() { -// +/// Default constructor } //______________________________________________________________________________ -AliMUONSt2GeometryBuilderV2& -AliMUONSt2GeometryBuilderV2::operator = (const AliMUONSt2GeometryBuilderV2& rhs) +AliMUONSt2GeometryBuilderV2::~AliMUONSt2GeometryBuilderV2() { -// Protected assignement operator - - // check assignement to self - if (this == &rhs) return *this; - - AliFatal("Assignment operator is not implemented."); - - return *this; +/// Destructor } // @@ -94,6 +75,7 @@ AliMUONSt2GeometryBuilderV2::operator = (const AliMUONSt2GeometryBuilderV2& rhs) //______________________________________________________________________________ void AliMUONSt2GeometryBuilderV2::CreateGeometry() { +/// Geometry construction // //******************************************************************** @@ -130,604 +112,610 @@ void AliMUONSt2GeometryBuilderV2::CreateGeometry() gMC->Gsvolu("SQM4","TUBS", idAir, tpar1, 5); - //================================================================================== -// Plane +// Plane //================================================================================== //Thickness of variour parts - Float_t zCbb = 0.04; //cathode pcb - Float_t zCu = 0.004; // eff. cu in cathode pcb - Float_t zRoha = 2.5; // Rhocell - Float_t zMeb = 0.08; //Mech. exit board - Float_t zEeb = 0.02; //Effective electronic exit board - //Z-positions of various parts--- in Plane-1 - - Float_t zposCbb = 0.25; // 2.5 mm => gap between anode & chatode plane - Float_t zposCu =zposCbb + zCbb; - Float_t zposRoha =zposCu + zCu; - Float_t zposMeb =zposRoha + zRoha; - Float_t zposEeb =zposMeb + zMeb; - - Float_t zposCbbBar = zposCbb + zCbb/2.; //for segment 0 & 6 - Float_t zposCuBar = zCbb/2. + zCu/2.; - Float_t zposRohaBar = zCu/2. + zRoha/2.; - Float_t zposMebBar = zRoha/2. + zMeb/2.; - Float_t zposEebBar = zMeb/2. + zEeb/2.; + Float_t zCu = 0.005; // eff. cu in cathode pcb + Float_t zCbb = 0.04; //cathode pcb + Float_t zRoha = 2.5; // Rhocell + Float_t zMeb = 0.04; //Mech. exit board //0.08 + Float_t zEeb = 0.04; //Effective electronic readout board //0.02 + + //starting Z-positions of various parts--- in Plane-1 + Float_t zposCu = 0.25; // 2.5 mm => gap between anode & chatode plane + Float_t zposCbb = zposCu + zCu; + Float_t zposRoha = zposCbb + zCbb; + Float_t zposMeb = zposRoha + zRoha; + Float_t zposEeb = zposMeb + zMeb; + + Float_t zposCuBar = zposCu + zCu/2.; //for segment 0 & 6 + Float_t zposCbbBar = zposCbb + zCbb/2.; + Float_t zposRohaBar = zposRoha + zRoha/2.; + Float_t zposMebBar = zposMeb + zMeb/2.; + Float_t zposEebBar = zposEeb + zEeb/2.; //Cathode PCB + Copper sheet + Rohacell + mech exit board + eff. electronic exit board //Segment-0 ~~~Horizantal box ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - Float_t bparH[3]; - bparH[0] = 94.5/2.; // extension beyond 0 deg in x direction - bparH[1] = 1.2/2.; // extension beyond 0 deg in y direction (3.7[total extn] - 2.5[frame dim]) - bparH[2] = zCbb/2.; // thickness of cathode sheet in z direction - gMC->Gsvolu("CB0L", "BOX", idPCB, bparH, 3); - - bparH[2] = zCu/2.; //Thickness of Copper sheet - gMC->Gsvolu("CU0L", "BOX", idCU, bparH, 3); - - bparH[2] = zRoha/2.; //Thickness of Roha cell - gMC->Gsvolu("RH0L", "BOX", idRoha, bparH, 3); - - bparH[2] = zMeb/2; //Thickness of mechanical exit board - gMC->Gsvolu("MB0L", "BOX", idPCB, bparH, 3); - - bparH[2] = zEeb/2; //Thickness of effective electronic exit board - gMC->Gsvolu("EB0L", "BOX", idCU, bparH, 3); + Float_t bparH[3]; + bparH[0] = 95.5/2.; // extension beyond 0 deg in x direction // 94.5 + bparH[1] = 1.2/2.; // extension beyond 0 deg in y direction (3.7[total extn] - 2.5[frame dim]) + bparH[2] = zCu/2.; //Thickness of Copper sheet + gMC->Gsvolu("SCU0L", "BOX", idCU, bparH, 3); + + bparH[2] = zCbb/2.; // thickness of cathode sheet in z direction + gMC->Gsvolu("SCB0L", "BOX", idPCB, bparH, 3); + + bparH[2] = zRoha/2.; //Thickness of Roha cell + gMC->Gsvolu("SRH0L", "BOX", idRoha, bparH, 3); - //Segment-1 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - Float_t pgpar[10]; // polygon - pgpar[0] = 0.; // initial angle - pgpar[1] = 90.; // increment in angle starting from initial angle - pgpar[2] = 5; // number of side - pgpar[3] = 2.; // number of plane - pgpar[4] = 0.; // z-position of the first plane - pgpar[5] = 23.1; // innner radius first plane - pgpar[6] = 117.6; // outer radious first plane - pgpar[7] = pgpar[4] + zCbb; // z-position of the second plane - pgpar[8] = pgpar[5]; // innner radius of second plane - pgpar[9] = pgpar[6]; // outer radious of second plane - gMC->Gsvolu("CB1L", "PGON", idPCB, pgpar, 10); - - pgpar[7] = pgpar[4] + zCu; // Thickness of copper-sheet - gMC->Gsvolu("CU1L", "PGON", idCU, pgpar, 10); - - pgpar[7] = pgpar[4] + zRoha; // Thickness of Roha cell - gMC->Gsvolu("RH1L", "PGON", idRoha, pgpar, 10); - - pgpar[7] = pgpar[4] + zMeb; // Thickness of mechanical exit board - gMC->Gsvolu("MB1L", "PGON", idPCB, pgpar, 10); + bparH[2] = zMeb/2; //Thickness of mechanical exit board + gMC->Gsvolu("SMB0L", "BOX", idPCB, bparH, 3); + + bparH[2] = zEeb/2; //Thickness of effective electronic exit board + gMC->Gsvolu("SEB0L", "BOX", idCU, bparH, 3); + + //Segment-1 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + Float_t pgpar[10]; // polygon + pgpar[0] = 0.; // initial angle + pgpar[1] = 90.; // increment in angle starting from initial angle + pgpar[2] = 5; // number of side + pgpar[3] = 2.; // number of plane + pgpar[4] = -zCu/2.; // z-position of the first plane + pgpar[5] = 23.1; // innner radius first plane + pgpar[6] = 117.6; // outer radious first plane + pgpar[7] = zCu/2.; // z-position of the second plane + pgpar[8] = pgpar[5]; // innner radius of second plane + pgpar[9] = pgpar[6]; // outer radious of second plane + gMC->Gsvolu("SCU1L", "PGON", idCU, pgpar, 10); + + pgpar[4] = -zCbb/2.; // z-position of the first plane + pgpar[7] = zCbb/2.; // Thickness of copper-sheet + gMC->Gsvolu("SCB1L", "PGON", idPCB, pgpar, 10); + + pgpar[4] = -zRoha/2.; // Thickness of Roha cell + pgpar[7] = zRoha/2.; // Thickness of Roha cell + gMC->Gsvolu("SRH1L", "PGON", idRoha, pgpar, 10); - pgpar[7] = pgpar[4] + zEeb; // Thickness of electronic exit board - gMC->Gsvolu("EB1L", "PGON", idCU, pgpar, 10); + pgpar[4] = -zMeb/2.; // Thickness of mechanical exit board + pgpar[7] = zMeb/2.; // Thickness of mechanical exit board + gMC->Gsvolu("SMB1L", "PGON", idPCB, pgpar, 10); + + pgpar[4] = -zEeb/2.; // Thickness of electronic readout board + pgpar[7] = zEeb/2.; // Thickness of electronic readout board + gMC->Gsvolu("SEB1L", "PGON", idCU, pgpar, 10); //Segment-2 - vertical box (simalar to horizontal bar as in Segment 0)~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - Float_t bparV[3]; - - bparV[0] = 1.0/2.; - bparV[1] = 94.5/2.; - bparV[2] = zCbb/2.; - gMC->Gsvolu("CB2L", "BOX", idPCB, bparV, 3); - - bparV[2] = zCu/2.; - gMC->Gsvolu("CU2L", "BOX", idCU, bparV, 3); - - bparV[2] = zRoha/2.; - gMC->Gsvolu("RH2L", "BOX", idRoha, bparV, 3); - - bparV[2] = zMeb/2; - gMC->Gsvolu("MB2L", "BOX", idPCB, bparV, 3); - - bparV[2] = zEeb/2; - gMC->Gsvolu("EB2L", "BOX", idCU, bparV, 3); - -//....(Setting posion of Segment 0,1,2).................................................................. - - Float_t xposHorBox = bparH[0] + 23.1; // 23.1 = 20.6(inner radius of qrd) + 2.5 (width of frame) - Float_t yposHorBox = -bparH[1]; + Float_t bparV[3]; + + bparV[0] = 1.0/2.; + bparV[1] = 95.5/2.; // 94.5 + bparV[2] = zCu/2.; + gMC->Gsvolu("SCU2L", "BOX", idCU, bparV, 3); + + bparV[2] = zCbb/2.; + gMC->Gsvolu("SCB2L", "BOX", idPCB, bparV, 3); + + + bparV[2] = zRoha/2.; + gMC->Gsvolu("SRH2L", "BOX", idRoha, bparV, 3); + + bparV[2] = zMeb/2; + gMC->Gsvolu("SMB2L", "BOX", idPCB, bparV, 3); + + bparV[2] = zEeb/2; + gMC->Gsvolu("SEB2L", "BOX", idCU, bparV, 3); - Float_t xposVerBox = -bparV[0]; - Float_t yposVerBox = bparV[1] + 23.1; - - - -//Positioning the PCB -// - - // chamber 3 - gMC->Gspos("CB0L",1, "SQM3",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY");// box horizontal - gMC->Gspos("CB1L",1, "SQM3", 0.0,0.0,zposCbb,0, "ONLY"); - gMC->Gspos("CB2L",1, "SQM3",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY");// box vertical - - gMC->Gspos("CB0L",2, "SQM3",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY");// box horizontal - gMC->Gspos("CB1L",2, "SQM3", 0.0,0.0,-(zposCbb+ zCbb),0, "ONLY"); - gMC->Gspos("CB2L",2, "SQM3",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY");// box vertical - - // chamber 4 - gMC->Gspos("CB0L",3, "SQM4",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY"); - gMC->Gspos("CB1L",3, "SQM4", 0.0,0.0,zposCbb,0, "ONLY"); - gMC->Gspos("CB2L",3, "SQM4",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY"); - - gMC->Gspos("CB0L",4, "SQM4",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY"); - gMC->Gspos("CB1L",4, "SQM4", 0.0,0.0,-(zposCbb+ zCbb),0, "ONLY"); - gMC->Gspos("CB2L",4, "SQM4",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY"); - - -//Positioning Copper sheet - - // chamber 3 - - gMC->Gspos("CU0L",1, "SQM3",xposHorBox,yposHorBox,zposCuBar,0, "ONLY"); - gMC->Gspos("CU1L",1, "SQM3", 0.0,0.0,zposCu,0, "ONLY"); - gMC->Gspos("CU2L",1, "SQM3",xposVerBox,yposVerBox,zposCuBar,0, "ONLY"); - - gMC->Gspos("CU0L",2, "SQM3",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY"); - gMC->Gspos("CU1L",2, "SQM3", 0.0,0.0,-(zposCu+ zCu),0, "ONLY"); - gMC->Gspos("CU2L",2, "SQM3",xposVerBox,yposVerBox,zposCuBar,0, "ONLY"); - // chamber 4 - - gMC->Gspos("CU0L",3, "SQM4",xposHorBox,yposHorBox,zposCuBar,0, "ONLY"); - gMC->Gspos("CU1L",3, "SQM4", 0.0,0.0,zposCu,0, "ONLY"); - gMC->Gspos("CU2L",3, "SQM4",xposVerBox,yposVerBox,zposCuBar,0, "ONLY"); - - gMC->Gspos("CU0L",4, "SQM4",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY"); - gMC->Gspos("CU1L",4, "SQM4", 0.0,0.0,-(zposCu+ zCu),0, "ONLY"); - gMC->Gspos("CU2L",4, "SQM4",xposVerBox,yposVerBox,zposCuBar,0, "ONLY"); - -//Positioning Roha cell - - // chamber 3 - - gMC->Gspos("RH0L",1, "SQM3",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY"); - gMC->Gspos("RH1L",1, "SQM3", 0.0,0.0,zposRoha,0, "ONLY"); - gMC->Gspos("RH2L",1, "SQM3",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY"); - - gMC->Gspos("RH0L",2, "SQM3",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY"); - gMC->Gspos("RH1L",2, "SQM3", 0.0,0.0,-(zposRoha+ zRoha),0, "ONLY"); - gMC->Gspos("RH2L",2, "SQM3",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY"); - - // chamber 4 - - gMC->Gspos("RH0L",3, "SQM4",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY"); - gMC->Gspos("RH1L",3, "SQM4", 0.0,0.0,zposRoha,0, "ONLY"); - gMC->Gspos("RH2L",3, "SQM4",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY"); - - gMC->Gspos("RH0L",4, "SQM4",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY"); - gMC->Gspos("RH1L",4, "SQM4", 0.0,0.0,-(zposRoha+ zRoha),0, "ONLY"); - gMC->Gspos("RH2L",4, "SQM4",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY"); - -//Positioning of Mech. exit board - - // chamber 3 - - gMC->Gspos("MB0L",1, "SQM3",xposHorBox,yposHorBox,zposMebBar,0, "ONLY"); - gMC->Gspos("MB1L",1, "SQM3", 0.0,0.0,zposMeb,0, "ONLY"); - gMC->Gspos("MB2L",1, "SQM3",xposVerBox,yposVerBox,zposMebBar,0, "ONLY"); - - gMC->Gspos("MB0L",2, "SQM3",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY"); - gMC->Gspos("MB1L",2, "SQM3", 0.0,0.0,-(zposMeb+ zMeb),0, "ONLY"); - gMC->Gspos("MB2L",2, "SQM3",xposVerBox,yposVerBox,zposMebBar,0, "ONLY"); - // chamber 4 - - gMC->Gspos("MB0L",3, "SQM4",xposHorBox,yposHorBox,zposMebBar,0, "ONLY"); - gMC->Gspos("MB1L",3, "SQM4", 0.0,0.0,zposMeb,0, "ONLY"); - gMC->Gspos("MB2L",3, "SQM4",xposVerBox,yposVerBox,zposMebBar,0, "ONLY"); - - gMC->Gspos("MB0L",4, "SQM4",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY"); - gMC->Gspos("MB1L",4, "SQM4", 0.0,0.0,-(zposMeb+ zMeb),0, "ONLY"); - gMC->Gspos("MB2L",4, "SQM4",xposVerBox,yposVerBox,zposMebBar,0, "ONLY"); - -//Positioning of Electronic exit board - - - // chamber 3 - - gMC->Gspos("EB0L",1, "SQM3",xposHorBox,yposHorBox,zposEebBar,0, "ONLY"); - gMC->Gspos("EB1L",1, "SQM3", 0.0,0.0,zposEeb,0, "ONLY"); - gMC->Gspos("EB2L",1, "SQM3",xposVerBox,yposVerBox,zposEebBar,0, "ONLY"); - - gMC->Gspos("EB0L",2, "SQM3",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY"); - gMC->Gspos("EB1L",2, "SQM3", 0.0,0.0,-(zposEeb+ zEeb),0, "ONLY"); - gMC->Gspos("EB2L",2, "SQM3",xposVerBox,yposVerBox,zposEebBar,0, "ONLY"); - // chamber 4 - - gMC->Gspos("EB0L",3, "SQM4",xposHorBox,yposHorBox,zposEebBar,0, "ONLY"); - gMC->Gspos("EB1L",3, "SQM4", 0.0,0.0,zposEeb,0, "ONLY"); - gMC->Gspos("EB2L",3, "SQM4",xposVerBox,yposVerBox,zposEebBar,0, "ONLY"); - - gMC->Gspos("EB0L",4, "SQM4",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY"); - gMC->Gspos("EB1L",4, "SQM4", 0.0,0.0,-(zposEeb+ zEeb),0, "ONLY"); - gMC->Gspos("EB2L",4, "SQM4",xposVerBox,yposVerBox,zposEebBar,0, "ONLY"); - - -//---------------------------------------------------------------------- -// Frames -//---------------------------------------------------------------------- -//Frame-1 - Float_t frame1[3] ; - frame1[0] = 101.0/2.; //100.6 = 94.5 + 2.5 + 3.6 - frame1[1] = 2.5/2.; - frame1[2] = 5.0/2.; + //....(Setting posion of Segment 0,1,2).................................................................. + + Float_t xposHorBox = bparH[0] + 23.1; // 23.1 = 20.6(inner radius of qrd) + 2.5 (width of frame) + Float_t yposHorBox = -bparH[1]; + + Float_t xposVerBox = -bparV[0]; + Float_t yposVerBox = bparV[1] + 23.1; + + + //Positioning the PCB + + + // chamber 3 + gMC->Gspos("SCB0L",1, "SQM3",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY"); + gMC->Gspos("SCB1L",1, "SQM3", 0.0,0.0,zposCbbBar,0, "ONLY"); + gMC->Gspos("SCB2L",1, "SQM3",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY"); + + gMC->Gspos("SCB0L",2, "SQM3",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY"); + gMC->Gspos("SCB1L",2, "SQM3", 0.0,0.0,-zposCbbBar,0, "ONLY"); + gMC->Gspos("SCB2L",2, "SQM3",xposVerBox,yposVerBox,-zposCbbBar,0, "ONLY"); + + // chamber 4 + gMC->Gspos("SCB0L",3, "SQM4",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY"); + gMC->Gspos("SCB1L",3, "SQM4", 0.0,0.0,zposCbbBar,0, "ONLY"); + gMC->Gspos("SCB2L",3, "SQM4",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY"); + + gMC->Gspos("SCB0L",4, "SQM4",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY"); + gMC->Gspos("SCB1L",4, "SQM4", 0.0,0.0,-zposCbbBar,0, "ONLY"); + gMC->Gspos("SCB2L",4, "SQM4",xposVerBox,yposVerBox,-zposCbbBar,0, "ONLY"); + - gMC->Gsvolu("FRM1", "BOX", idPGF30, frame1, 3); //Frame - 1 // fill with pkk GF30 - - Float_t arib1[3]; - arib1[0] = frame1[0]; - arib1[1] = 0.9/2.; - arib1[2] =(frame1[2]-0.95)/2.0; - - gMC->Gsvolu("FRA1", "BOX", idAir, arib1, 3); // fill with air - - Float_t xposarib1 = 0; - Float_t yposarib1 = -frame1[1] + arib1[1]; - Float_t zposarib1 = frame1[2] - arib1[2]; - - gMC->Gspos("FRA1",1, "FRM1", xposarib1, yposarib1, zposarib1,0, "ONLY"); //replace pkk GF30 with air(b) - gMC->Gspos("FRA1",2, "FRM1", xposarib1, yposarib1, -zposarib1,0, "ONLY"); //replace pkk GF30 with air(nb) - - Float_t rrib1[3]; - rrib1[0] = frame1[0]; - rrib1[1] = 0.6/2.; - rrib1[2] =(frame1[2]-0.95)/2.0; - - gMC->Gsvolu("FRR1", "BOX", idRoha, rrib1, 3); // fill with rohacell - - Float_t xposrrib1 = 0.0; - Float_t yposrrib1 = frame1[1] - rrib1[1]; - Float_t zposrrib1 = frame1[2] - rrib1[2]; - - gMC->Gspos("FRR1",1, "FRM1", xposrrib1, yposrrib1, zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell - gMC->Gspos("FRR1",2, "FRM1", xposrrib1, yposrrib1, -zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell + //Positioning Copper sheet + + // chamber 3 + gMC->Gspos("SCU0L",1, "SQM3",xposHorBox,yposHorBox,zposCuBar,0, "ONLY"); + gMC->Gspos("SCU1L",1, "SQM3", 0.0,0.0,zposCuBar,0, "ONLY"); + gMC->Gspos("SCU2L",1, "SQM3",xposVerBox,yposVerBox,zposCuBar,0, "ONLY"); + + gMC->Gspos("SCU0L",2, "SQM3",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY"); + gMC->Gspos("SCU1L",2, "SQM3", 0.0,0.0,-zposCuBar,0, "ONLY"); + gMC->Gspos("SCU2L",2, "SQM3",xposVerBox,yposVerBox,-zposCuBar,0, "ONLY"); + + // chamber 4 + gMC->Gspos("SCU0L",3, "SQM4",xposHorBox,yposHorBox,zposCuBar,0, "ONLY"); + gMC->Gspos("SCU1L",3, "SQM4", 0.0,0.0,zposCuBar,0, "ONLY"); + gMC->Gspos("SCU2L",3, "SQM4",xposVerBox,yposVerBox,zposCuBar,0, "ONLY"); + + gMC->Gspos("SCU0L",4, "SQM4",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY"); + gMC->Gspos("SCU1L",4, "SQM4", 0.0,0.0,-zposCuBar,0, "ONLY"); + gMC->Gspos("SCU2L",4, "SQM4",xposVerBox,yposVerBox,-zposCuBar,0, "ONLY"); - Float_t xposFr1 = frame1[0] + 20.6; - Float_t yposFr1 = -3.7 + frame1[1] ; - Float_t zposFr1 = 0.0; - - gMC->Gspos("FRM1",1, "SQM3", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1 - gMC->Gspos("FRM1",2, "SQM4", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1 - - -//...................................................................................... -//Frame-2 - - Float_t frame2[3]; - frame2[0] = 4.0/2.; - frame2[1] = 1.2/2.; - frame2[2] = 5.0/2; - - gMC->Gsvolu("FRM2", "BOX", idPGF30, frame2, 3); //Frame - 2 - - Float_t rrib2[3]; - rrib2[0] = frame2[0]-1.0/2.0; - rrib2[1] = frame2[1]; - rrib2[2] =(frame2[2]-0.95)/2.0; - - gMC->Gsvolu("FRR2", "BOX", idRoha, rrib2, 3); - - Float_t xposrrib2 = -1.0/2.0; - Float_t yposrrib2 = 0.0; - Float_t zposrrib2 = frame2[2] - rrib2[2]; - - gMC->Gspos("FRR2",1, "FRM2", xposrrib2, yposrrib2, zposrrib2,0, "ONLY");//replace pkk GF30 with rohacell - gMC->Gspos("FRR2",2, "FRM2", xposrrib2, yposrrib2, -zposrrib2,0, "ONLY");//replace pkk GF30 with roha - - + //Positioning Roha cell + + // chamber 3 + + gMC->Gspos("SRH0L",1, "SQM3",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY"); + gMC->Gspos("SRH1L",1, "SQM3", 0.0,0.0,zposRohaBar,0, "ONLY"); + gMC->Gspos("SRH2L",1, "SQM3",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY"); + + gMC->Gspos("SRH0L",2, "SQM3",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY"); + gMC->Gspos("SRH1L",2, "SQM3", 0.0,0.0,-zposRohaBar,0, "ONLY"); + gMC->Gspos("SRH2L",2, "SQM3",xposVerBox,yposVerBox,-zposRohaBar,0, "ONLY"); + + // chamber 4 + + gMC->Gspos("SRH0L",3, "SQM4",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY"); + gMC->Gspos("SRH1L",3, "SQM4", 0.0,0.0,zposRohaBar,0, "ONLY"); + gMC->Gspos("SRH2L",3, "SQM4",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY"); + + gMC->Gspos("SRH0L",4, "SQM4",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY"); + gMC->Gspos("SRH1L",4, "SQM4", 0.0,0.0,-zposRohaBar,0, "ONLY"); + gMC->Gspos("SRH2L",4, "SQM4",xposVerBox,yposVerBox,-zposRohaBar,0, "ONLY"); - Float_t xposFr2 = frame2[0] + 117.6; - Float_t yposFr2 = -frame2[1]; - Float_t zposFr2 = 0.0; + //Positioning of Mech. exit board + + // chamber 3 + + gMC->Gspos("SMB0L",1, "SQM3",xposHorBox,yposHorBox,zposMebBar,0, "ONLY"); + gMC->Gspos("SMB1L",1, "SQM3", 0.0,0.0,zposMebBar,0, "ONLY"); + gMC->Gspos("SMB2L",1, "SQM3",xposVerBox,yposVerBox,zposMebBar,0, "ONLY"); + + gMC->Gspos("SMB0L",2, "SQM3",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY"); + gMC->Gspos("SMB1L",2, "SQM3", 0.0,0.0,-zposMebBar,0, "ONLY"); + gMC->Gspos("SMB2L",2, "SQM3",xposVerBox,yposVerBox,-zposMebBar,0, "ONLY"); + // chamber 4 + + gMC->Gspos("SMB0L",3, "SQM4",xposHorBox,yposHorBox,zposMebBar,0, "ONLY"); + gMC->Gspos("SMB1L",3, "SQM4", 0.0,0.0,zposMebBar,0, "ONLY"); + gMC->Gspos("SMB2L",3, "SQM4",xposVerBox,yposVerBox,zposMebBar,0, "ONLY"); + + gMC->Gspos("SMB0L",4, "SQM4",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY"); + gMC->Gspos("SMB1L",4, "SQM4", 0.0,0.0,-zposMebBar,0, "ONLY"); + gMC->Gspos("SMB2L",4, "SQM4",xposVerBox,yposVerBox,-zposMebBar,0, "ONLY"); + + //Positioning of Electronic exit board + + + // chamber 3 - gMC->Gspos("FRM2",1, "SQM3", xposFr2, yposFr2, zposFr2,0, "ONLY");//global positing of frame in SQM3 - gMC->Gspos("FRM2",2, "SQM4", xposFr2, yposFr2, zposFr2,0, "ONLY");//global positing of frame in SQM4 - + gMC->Gspos("SEB0L",1, "SQM3",xposHorBox,yposHorBox,zposEebBar,0, "ONLY"); + gMC->Gspos("SEB1L",1, "SQM3", 0.0,0.0,zposEebBar,0, "ONLY"); + gMC->Gspos("SEB2L",1, "SQM3",xposVerBox,yposVerBox,zposEebBar,0, "ONLY"); -//...................................................................................... -//Frame-3 - - Float_t pgparFr3[10]; - pgparFr3[0] = 0.; - pgparFr3[1] = 90.; - pgparFr3[2] = 5; - pgparFr3[3] = 2.; - pgparFr3[4] = 0.; - pgparFr3[5] = 117.6; - pgparFr3[6] = 121.6; - pgparFr3[7] = pgparFr3[4] + 5.0; - pgparFr3[8] = pgparFr3[5]; - pgparFr3[9] = pgparFr3[6]; - - gMC->Gsvolu("FRM3", "PGON", idPGF30, pgparFr3, 10); - - Float_t pgparRrib3[10]; - pgparRrib3[0] = 0.; - pgparRrib3[1] = 90.; - pgparRrib3[2] = 5; - pgparRrib3[3] = 2.; - pgparRrib3[4] = 0.; - pgparRrib3[5] = 117.6; - pgparRrib3[6] = 120.6; - pgparRrib3[7] = pgparRrib3[4] +1.55 ; - pgparRrib3[8] = pgparRrib3[5]; - pgparRrib3[9] = pgparRrib3[6]; - - gMC->Gsvolu("FRR3", "PGON", idRoha, pgparRrib3, 10); - - Float_t xposrrib3 = 0.0; - Float_t yposrrib3 = 0.0; - Float_t zposrrib3 = 0.0; + gMC->Gspos("SEB0L",2, "SQM3",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY"); + gMC->Gspos("SEB1L",2, "SQM3", 0.0,0.0,-zposEebBar,0, "ONLY"); + gMC->Gspos("SEB2L",2, "SQM3",xposVerBox,yposVerBox,-zposEebBar,0, "ONLY"); + // chamber 4 + + gMC->Gspos("SEB0L",3, "SQM4",xposHorBox,yposHorBox,zposEebBar,0, "ONLY"); + gMC->Gspos("SEB1L",3, "SQM4", 0.0,0.0,zposEebBar,0, "ONLY"); + gMC->Gspos("SEB2L",3, "SQM4",xposVerBox,yposVerBox,zposEebBar,0, "ONLY"); + + gMC->Gspos("SEB0L",4, "SQM4",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY"); + gMC->Gspos("SEB1L",4, "SQM4", 0.0,0.0,-zposEebBar,0, "ONLY"); + gMC->Gspos("SEB2L",4, "SQM4",xposVerBox,yposVerBox,-zposEebBar,0, "ONLY"); + + + //---------------------------------------------------------------------- + // Frames + //---------------------------------------------------------------------- + //Frame-1 + + Float_t frame1[3] ; + frame1[0] = 101.0/2.; //100.6 = 94.5 + 2.5 + 3.6 + frame1[1] = 2.5/2.; + frame1[2] = 5.0/2.; + + gMC->Gsvolu("SFRM1", "BOX", idPGF30, frame1, 3); //Frame - 1 // fill with pkk GF30 + + Float_t arib1[3]; + arib1[0] = frame1[0]; + arib1[1] = 0.9/2.; + arib1[2] =(frame1[2]-0.95)/2.0; + + gMC->Gsvolu("SFRA1", "BOX", idAir, arib1, 3); // fill with air + + Float_t xposarib1 = 0; + Float_t yposarib1 = -frame1[1] + arib1[1]; + Float_t zposarib1 = frame1[2] - arib1[2]; + + gMC->Gspos("SFRA1",1, "SFRM1", xposarib1, yposarib1, zposarib1,0, "ONLY"); //replace pkk GF30 with air(b) + gMC->Gspos("SFRA1",2, "SFRM1", xposarib1, yposarib1, -zposarib1,0, "ONLY"); //replace pkk GF30 with air(nb) + + Float_t rrib1[3]; + rrib1[0] = frame1[0]; + rrib1[1] = 0.6/2.; + rrib1[2] =(frame1[2]-0.95)/2.0; + + gMC->Gsvolu("SFRR1", "BOX", idRoha, rrib1, 3); // fill with rohacell + + Float_t xposrrib1 = 0.0; + Float_t yposrrib1 = frame1[1] - rrib1[1]; + Float_t zposrrib1 = frame1[2] - rrib1[2]; + + gMC->Gspos("SFRR1",1, "SFRM1", xposrrib1, yposrrib1, zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell + gMC->Gspos("SFRR1",2, "SFRM1", xposrrib1, yposrrib1, -zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell + + + Float_t xposFr1 = frame1[0] + 20.6; + Float_t yposFr1 = -3.7 + frame1[1] ; + Float_t zposFr1 = 0.0; + + gMC->Gspos("SFRM1",1, "SQM3", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1 + gMC->Gspos("SFRM1",2, "SQM4", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1 + + + //...................................................................................... + //Frame-2 - gMC->Gspos("FRR3",1, "FRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY"); + Float_t frame2[3]; + frame2[0] = 4.0/2.; + frame2[1] = 1.2/2.; + frame2[2] = 5.0/2; + + gMC->Gsvolu("SFRM2", "BOX", idPGF30, frame2, 3); //Frame - 2 + + Float_t rrib2[3]; + rrib2[0] = frame2[0]-1.0/2.0; + rrib2[1] = frame2[1]; + rrib2[2] =(frame2[2]-0.95)/2.0; - zposrrib3 = 3.45; - - gMC->Gspos("FRR3",2, "FRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY"); - + gMC->Gsvolu("SFRR2", "BOX", idRoha, rrib2, 3); - - Float_t xposFr3 = 0.0; - Float_t yposFr3 = 0.0; - Float_t zposFr3 = -frame1[2]; - - gMC->Gspos("FRM3",1, "SQM3", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1 - gMC->Gspos("FRM3",2, "SQM4", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1 - - -//...................................................................................... -//Frame-4 - - Float_t frame4[3]; - frame4[0] = 1.0/2.; - frame4[1] = 4.0/2.; - frame4[2] = frame1[2]; - - gMC->Gsvolu("FRM4", "BOX", idPGF30, frame4, 3); + Float_t xposrrib2 = -1.0/2.0; + Float_t yposrrib2 = 0.0; + Float_t zposrrib2 = frame2[2] - rrib2[2]; + + gMC->Gspos("SFRR2",1, "SFRM2", xposrrib2, yposrrib2, zposrrib2,0, "ONLY");//replace pkk GF30 with rohacell + gMC->Gspos("SFRR2",2, "SFRM2", xposrrib2, yposrrib2, -zposrrib2,0, "ONLY");//replace pkk GF30 with roha + - Float_t rrib4[3]; - rrib4[0] = frame4[0]; - rrib4[1] = frame4[1]-1.0/2; - rrib4[2] =(frame4[2]-0.95)/2.0; + + Float_t xposFr2 = frame2[0] + 117.6; + Float_t yposFr2 = -frame2[1]; + Float_t zposFr2 = 0.0; + + gMC->Gspos("SFRM2",1, "SQM3", xposFr2, yposFr2, zposFr2,0, "MANY");//global positing of frame in SQM3 + gMC->Gspos("SFRM2",2, "SQM4", xposFr2, yposFr2, zposFr2,0, "MANY");//global positing of frame in SQM4 + + + //...................................................................................... + + //Frame-3 + + Float_t pgparFr3[10]; + pgparFr3[0] = 0.; + pgparFr3[1] = 90.; + pgparFr3[2] = 5; + pgparFr3[3] = 2.; + pgparFr3[4] = 0.; + pgparFr3[5] = 117.6; + pgparFr3[6] = 121.6; + pgparFr3[7] = pgparFr3[4] + 5.0; + pgparFr3[8] = pgparFr3[5]; + pgparFr3[9] = pgparFr3[6]; + + gMC->Gsvolu("SFRM3", "PGON", idPGF30, pgparFr3, 10); + + Float_t pgparRrib3[10]; + pgparRrib3[0] = 0.; + pgparRrib3[1] = 90.; + pgparRrib3[2] = 5; + pgparRrib3[3] = 2.; + pgparRrib3[4] = 0.; + pgparRrib3[5] = 117.6; + pgparRrib3[6] = 120.6; + pgparRrib3[7] = pgparRrib3[4] +1.55 ; + pgparRrib3[8] = pgparRrib3[5]; + pgparRrib3[9] = pgparRrib3[6]; + + gMC->Gsvolu("SFRR3", "PGON", idRoha, pgparRrib3, 10); + + Float_t xposrrib3 = 0.0; + Float_t yposrrib3 = 0.0; + Float_t zposrrib3 = 0.0; - gMC->Gsvolu("FRR4", "BOX", idRoha, rrib4, 3); + gMC->Gspos("SFRR3",1, "SFRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY"); + + zposrrib3 = 3.45; + + gMC->Gspos("SFRR3",2, "SFRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY"); + - Float_t xposrrib4 = 0.0; - Float_t yposrrib4 = -1.0/2; - Float_t zposrrib4 = frame4[2] - rrib4[2]; - - gMC->Gspos("FRR4",1, "FRM4", xposrrib4, yposrrib4, zposrrib4,0, "ONLY"); - gMC->Gspos("FRR4",2, "FRM4", xposrrib4, yposrrib4, -zposrrib4,0, "ONLY"); - - - - Float_t xposFr4 = -frame4[0]; - Float_t yposFr4 = -frame4[1] + 117.6; - Float_t zposFr4 = 0.0; - - gMC->Gspos("FRM4",1, "SQM3", xposFr4, yposFr4, zposFr4,0, "ONLY"); - gMC->Gspos("FRM4",2, "SQM4", xposFr4, yposFr4, zposFr4,0, "ONLY"); + + Float_t xposFr3 = 0.0; + Float_t yposFr3 = 0.0; + Float_t zposFr3 = -frame1[2]; + gMC->Gspos("SFRM3",1, "SQM3", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1 + gMC->Gspos("SFRM3",2, "SQM4", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1 + + + //...................................................................................... + //Frame-4 + + Float_t frame4[3]; + frame4[0] = 1.0/2.; + frame4[1] = 4.0/2.; + frame4[2] = frame1[2]; + + gMC->Gsvolu("SFRM4", "BOX", idPGF30, frame4, 3); + + Float_t rrib4[3]; + rrib4[0] = frame4[0]; + rrib4[1] = frame4[1]-1.0/2; + rrib4[2] =(frame4[2]-0.95)/2.0; -//...................................................................................... -//Frame-5 - Float_t frame5[3] ; - frame5[0] = 2.7/2.; - frame5[1] = 101.0/2.; - frame5[2] = 5.0/2.; + gMC->Gsvolu("SFRR4", "BOX", idRoha, rrib4, 3); + + Float_t xposrrib4 = 0.0; + Float_t yposrrib4 = -1.0/2; + Float_t zposrrib4 = frame4[2] - rrib4[2]; + + gMC->Gspos("SFRR4",1, "SFRM4", xposrrib4, yposrrib4, zposrrib4,0, "ONLY"); + gMC->Gspos("SFRR4",2, "SFRM4", xposrrib4, yposrrib4, -zposrrib4,0, "ONLY"); + - gMC->Gsvolu("FRM5", "BOX", idPGF30, frame5, 3); //Frame - 1 + + Float_t xposFr4 = -frame4[0]; + Float_t yposFr4 = -frame4[1] + 117.6; + Float_t zposFr4 = 0.0; + + gMC->Gspos("SFRM4",1, "SQM3", xposFr4, yposFr4, zposFr4,0, "MANY"); + gMC->Gspos("SFRM4",2, "SQM4", xposFr4, yposFr4, zposFr4,0, "MANY"); + + + //...................................................................................... + //Frame-5 + Float_t frame5[3] ; + frame5[0] = 2.7/2.; + frame5[1] = 101.0/2.; + frame5[2] = 5.0/2.; + + gMC->Gsvolu("SFRM5", "BOX", idPGF30, frame5, 3); //Frame - 1 - Float_t arib5[3]; - arib5[0] = 0.9/2.0; - arib5[1] = frame5[1]; - arib5[2] = (frame5[2]-0.95)/2.0; - - gMC->Gsvolu("FRA5", "BOX", idAir, arib5, 3); - - Float_t xposarib5 = -frame5[0] + arib5[0]; - Float_t yposarib5 = 0.0; - Float_t zposarib5 = frame5[2] - arib5[2]; - - gMC->Gspos("FRA5",1, "FRM5", xposarib5, yposarib5, zposarib5,0, "ONLY"); - gMC->Gspos("FRA5",2, "FRM5", xposarib5, yposarib5, -zposarib5,0, "ONLY"); - - Float_t rrib5[3]; - rrib5[0] = 0.8/2.0; - rrib5[1] = frame5[1]; - rrib5[2] = (frame5[2]-0.95)/2.0; - - gMC->Gsvolu("FRR5", "BOX", idRoha, rrib5, 3); - - Float_t xposrrib5 = frame5[0] - rrib5[0]; - Float_t yposrrib5 = 0.0; - Float_t zposrrib5 = frame5[2] - rrib5[2]; - - gMC->Gspos("FRR5",1, "FRM5", xposrrib5, yposrrib5, zposrrib5,0, "ONLY"); - gMC->Gspos("FRR5",2, "FRM5", xposrrib5, yposrrib5, -zposrrib5,0, "ONLY"); - - Float_t xposFr5 = -3.7 + frame5[0]; - Float_t yposFr5 = frame5[1] + 20.6; - Float_t zposFr5 = 0.0; - - gMC->Gspos("FRM5",1, "SQM3", xposFr5, yposFr5, zposFr5,0, "ONLY"); - gMC->Gspos("FRM5",2, "SQM4", xposFr5, yposFr5, zposFr5,0, "ONLY"); - -//...................................................................................... -//Frame -6 - - Float_t frame6[3]; - frame6[0] = 1.0/2.; - frame6[1] = 2.5/2.; - frame6[2] = frame1[2]; - - gMC->Gsvolu("FRM6", "BOX", idPGF30, frame6, 3); - - Float_t rrib6[3]; - rrib6[0] = frame6[0]; - rrib6[1] = 1.5/2.; - rrib6[2] =(frame2[2]-0.95)/2.0; - - gMC->Gsvolu("FRR6", "BOX", idRoha, rrib6, 3); + Float_t arib5[3]; + arib5[0] = 0.9/2.0; + arib5[1] = frame5[1]; + arib5[2] = (frame5[2]-0.95)/2.0; + + gMC->Gsvolu("SFRA5", "BOX", idAir, arib5, 3); + + Float_t xposarib5 = -frame5[0] + arib5[0]; + Float_t yposarib5 = 0.0; + Float_t zposarib5 = frame5[2] - arib5[2]; - Float_t xposrrib6 = 0.0; - Float_t yposrrib6 = 1.0/2.0; - Float_t zposrrib6 = frame6[2] - rrib6[2]; - - gMC->Gspos("FRR6",1, "FRM6", xposrrib6, yposrrib6, zposrrib6,0, "ONLY"); - gMC->Gspos("FRR6",2, "FRM6", xposrrib6, yposrrib6, -zposrrib6,0, "ONLY"); - - + gMC->Gspos("SFRA5",1, "SFRM5", xposarib5, yposarib5, zposarib5,0, "ONLY"); + gMC->Gspos("SFRA5",2, "SFRM5", xposarib5, yposarib5, -zposarib5,0, "ONLY"); + + Float_t rrib5[3]; + rrib5[0] = 0.8/2.0; + rrib5[1] = frame5[1]; + rrib5[2] = (frame5[2]-0.95)/2.0; + + gMC->Gsvolu("SFRR5", "BOX", idRoha, rrib5, 3); + + Float_t xposrrib5 = frame5[0] - rrib5[0]; + Float_t yposrrib5 = 0.0; + Float_t zposrrib5 = frame5[2] - rrib5[2]; + + gMC->Gspos("SFRR5",1, "SFRM5", xposrrib5, yposrrib5, zposrrib5,0, "ONLY"); + gMC->Gspos("SFRR5",2, "SFRM5", xposrrib5, yposrrib5, -zposrrib5,0, "ONLY"); + + Float_t xposFr5 = -3.7 + frame5[0]; + Float_t yposFr5 = frame5[1] + 20.6; + Float_t zposFr5 = 0.0; + + gMC->Gspos("SFRM5",1, "SQM3", xposFr5, yposFr5, zposFr5,0, "ONLY"); + gMC->Gspos("SFRM5",2, "SQM4", xposFr5, yposFr5, zposFr5,0, "ONLY"); + + //...................................................................................... + //Frame -6 + + Float_t frame6[3]; + frame6[0] = 1.0/2.; + frame6[1] = 2.5/2.; + frame6[2] = frame1[2]; + + gMC->Gsvolu("SFRM6", "BOX", idPGF30, frame6, 3); + + Float_t rrib6[3]; + rrib6[0] = frame6[0]; + rrib6[1] = 1.5/2.; + rrib6[2] =(frame2[2]-0.95)/2.0; - Float_t xposFr6 = -frame6[0]; - Float_t yposFr6 = frame6[1] + 20.6; - Float_t zposFr6 = 0.0; + gMC->Gsvolu("SFRR6", "BOX", idRoha, rrib6, 3); + + Float_t xposrrib6 = 0.0; + Float_t yposrrib6 = 1.0/2.0; + Float_t zposrrib6 = frame6[2] - rrib6[2]; - gMC->Gspos("FRM6",1, "SQM3", xposFr6, yposFr6, zposFr6,0, "ONLY"); - gMC->Gspos("FRM6",2, "SQM4", xposFr6, yposFr6, zposFr6,0, "ONLY"); - + gMC->Gspos("SFRR6",1, "SFRM6", xposrrib6, yposrrib6, zposrrib6,0, "ONLY"); + gMC->Gspos("SFRR6",2, "SFRM6", xposrrib6, yposrrib6, -zposrrib6,0, "ONLY"); + + + + Float_t xposFr6 = -frame6[0]; + Float_t yposFr6 = frame6[1] + 20.6; + Float_t zposFr6 = 0.0; + gMC->Gspos("SFRM6",1, "SQM3", xposFr6, yposFr6, zposFr6,0, "ONLY"); + gMC->Gspos("SFRM6",2, "SQM4", xposFr6, yposFr6, zposFr6,0, "ONLY"); + + //...................................................................................... //Frame - 7 inner pgon + + Float_t pgparFr7[10]; + pgparFr7[0] = 0.; + pgparFr7[1] = 90.; + pgparFr7[2] = 5; + pgparFr7[3] = 2.; + pgparFr7[4] = 0.; + pgparFr7[5] = 20.6; + pgparFr7[6] = 23.1; + pgparFr7[7] = pgparFr7[4] + 5.0; + pgparFr7[8] = pgparFr7[5]; + pgparFr7[9] = pgparFr7[6]; + + gMC->Gsvolu("SFRM7", "PGON", idPGF30, pgparFr7, 10); + + Float_t pgparRrib7[10]; + pgparRrib7[0] = 0.; + pgparRrib7[1] = 90.; + pgparRrib7[2] = 5; + pgparRrib7[3] = 2.; + pgparRrib7[4] = 0.; + pgparRrib7[5] = 21.6; + pgparRrib7[6] = 23.1; + pgparRrib7[7] = pgparRrib7[4] + 1.55; + pgparRrib7[8] = pgparRrib7[5]; + pgparRrib7[9] = pgparRrib7[6]; + + gMC->Gsvolu("SFRR7", "PGON", idRoha, pgparRrib7, 10); + + Float_t xposrrib7 = 0.0; + Float_t yposrrib7 = 0.0; + Float_t zposrrib7 = 0.0; + + gMC->Gspos("SFRR7",1, "SFRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY"); + + zposrrib7 = 3.45; + + gMC->Gspos("SFRR7",2, "SFRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY"); + + Float_t xposFr7 = 0.0; + Float_t yposFr7 = 0.0; + Float_t zposFr7 = -frame1[2]; + + gMC->Gspos("SFRM7",1, "SQM3", xposFr7, yposFr7, zposFr7,0, "ONLY"); + gMC->Gspos("SFRM7",2, "SQM4", xposFr7, yposFr7, zposFr7,0, "ONLY"); + + + //...................................................................................... + //Frame - 8 - Float_t pgparFr7[10]; - pgparFr7[0] = 0.; - pgparFr7[1] = 90.; - pgparFr7[2] = 5; - pgparFr7[3] = 2.; - pgparFr7[4] = 0.; - pgparFr7[5] = 20.6; - pgparFr7[6] = 23.1; - pgparFr7[7] = pgparFr7[4] + 5.0; - pgparFr7[8] = pgparFr7[5]; - pgparFr7[9] = pgparFr7[6]; - - gMC->Gsvolu("FRM7", "PGON", idPGF30, pgparFr7, 10); - - Float_t pgparRrib7[10]; - pgparRrib7[0] = 0.; - pgparRrib7[1] = 90.; - pgparRrib7[2] = 5; - pgparRrib7[3] = 2.; - pgparRrib7[4] = 0.; - pgparRrib7[5] = 21.6; - pgparRrib7[6] = 23.1; - pgparRrib7[7] = pgparRrib7[4] + 1.55; - pgparRrib7[8] = pgparRrib7[5]; - pgparRrib7[9] = pgparRrib7[6]; - - gMC->Gsvolu("FRR7", "PGON", idRoha, pgparRrib7, 10); - - Float_t xposrrib7 = 0.0; - Float_t yposrrib7 = 0.0; - Float_t zposrrib7 = 0.0; - - gMC->Gspos("FRR7",1, "FRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY"); - - zposrrib7 = 3.45; - - gMC->Gspos("FRR7",2, "FRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY"); - - Float_t xposFr7 = 0.0; - Float_t yposFr7 = 0.0; - Float_t zposFr7 = -frame1[2]; - - gMC->Gspos("FRM7",1, "SQM3", xposFr7, yposFr7, zposFr7,0, "ONLY"); - gMC->Gspos("FRM7",2, "SQM4", xposFr7, yposFr7, zposFr7,0, "ONLY"); - - -//...................................................................................... -//Frame - 8 - - - Float_t frame8[3] ; - frame8[0] = 2.5/2.0; - frame8[1] = 1.2/2.0; - frame8[2] = frame1[2]; - - gMC->Gsvolu("FRM8", "BOX", idPGF30, frame8, 3); //Frame - 2 - - Float_t rrib8[3]; - rrib8[0] = frame8[0]-1.0/2; - rrib8[1] = frame8[1]; - rrib8[2] =(frame8[2]-0.95)/2.0; - - gMC->Gsvolu("FRR8", "BOX", idRoha, rrib8, 3); - - Float_t xposrrib8 = -1.0/2; - Float_t yposrrib8 = 0.0; - Float_t zposrrib8 = frame8[2] - rrib8[2]; - - gMC->Gspos("FRR8",1, "FRM8", xposrrib8, yposrrib8, zposrrib8,0, "ONLY"); - gMC->Gspos("FRR8",2, "FRM8", xposrrib8, yposrrib8, -zposrrib8,0, "ONLY"); - - - - Float_t xposFr8 = frame8[0] + 20.6; - Float_t yposFr8 = -frame8[1]; - Float_t zposFr8 = 0.0; - - gMC->Gspos("FRM8",1, "SQM3", xposFr8, yposFr8, zposFr8,0, "ONLY"); - gMC->Gspos("FRM8",2, "SQM4", xposFr8, yposFr8, zposFr8,0, "ONLY"); + + Float_t frame8[3] ; + frame8[0] = 2.5/2.0; + frame8[1] = 1.2/2.0; + frame8[2] = frame1[2]; + + gMC->Gsvolu("SFRM8", "BOX", idPGF30, frame8, 3); //Frame - 2 - - - -//^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ - - Float_t zsenv = 0.5; // distance between two cathode plane - - //Segment-0 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - bparH[0] = 94.5/2.; - bparH[1] = 1.2/2.; - bparH[2] = zsenv/2.; - gMC->Gsvolu("C3G0", "BOX", idGas, bparH, 3); - gMC->Gsvolu("C4G0", "BOX", idGas, bparH, 3); + Float_t rrib8[3]; + rrib8[0] = frame8[0]-1.0/2; + rrib8[1] = frame8[1]; + rrib8[2] =(frame8[2]-0.95)/2.0; + + gMC->Gsvolu("SFRR8", "BOX", idRoha, rrib8, 3); + + Float_t xposrrib8 = -1.0/2; + Float_t yposrrib8 = 0.0; + Float_t zposrrib8 = frame8[2] - rrib8[2]; + + gMC->Gspos("SFRR8",1, "SFRM8", xposrrib8, yposrrib8, zposrrib8,0, "ONLY"); + gMC->Gspos("SFRR8",2, "SFRM8", xposrrib8, yposrrib8, -zposrrib8,0, "ONLY"); + + + + Float_t xposFr8 = frame8[0] + 20.6; + Float_t yposFr8 = -frame8[1]; + Float_t zposFr8 = 0.0; + + gMC->Gspos("SFRM8",1, "SQM3", xposFr8, yposFr8, zposFr8,0, "ONLY"); + gMC->Gspos("SFRM8",2, "SQM4", xposFr8, yposFr8, zposFr8,0, "ONLY"); + + + + //^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ + + Float_t zsenv = 0.5; // distance between two cathode plane + + //Segment-0 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + bparH[0] = 94.5/2.; + bparH[1] = 1.2/2.; + bparH[2] = zsenv/2.; + gMC->Gsvolu("SC3G0", "BOX", idGas, bparH, 3); + gMC->Gsvolu("SC4G0", "BOX", idGas, bparH, 3); + //Segment-1 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - pgpar[0] = 0.; - pgpar[1] = 90.; - pgpar[2] = 5; - pgpar[3] = 2.; - pgpar[4] = -zsenv/2.; - pgpar[5] = 23.1; - pgpar[6] = 117.6; - pgpar[7] = zsenv/2.; - pgpar[8] = pgpar[5]; - pgpar[9] = pgpar[6]; - gMC->Gsvolu("C3G1", "PGON", idGas, pgpar, 10); - gMC->Gsvolu("C4G1", "PGON", idGas, pgpar, 10); - - -//Segment-2 - vertical box ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - - bparV[0] = 1.0/2.; - bparV[1] = 95.5/2.; - bparV[2] = zsenv/2.; - gMC->Gsvolu("C3G2", "BOX", idGas, bparV, 3); - gMC->Gsvolu("C4G2", "BOX", idGas, bparV, 3); - -//........................................................................................... - - xposHorBox = bparH[0] + 23.1; - yposHorBox = -bparH[1]; - - xposVerBox = -bparV[0]; - yposVerBox = bparV[1] + 23.1; - - gMC->Gspos("C3G0",1, "SQM3", xposHorBox,yposHorBox,0.,0, "ONLY"); - gMC->Gspos("C3G1",1, "SQM3", 0.,0.,0.,0, "ONLY"); - gMC->Gspos("C3G2",1, "SQM3", xposVerBox,yposVerBox,0.,0, "ONLY"); - - - gMC->Gspos("C4G0",1, "SQM4", xposHorBox,yposHorBox,0.,0, "ONLY"); - gMC->Gspos("C4G1",1, "SQM4", 0.,0.,0.,0, "ONLY"); - gMC->Gspos("C4G2",1, "SQM4", xposVerBox,yposVerBox,0.,0, "ONLY"); + pgpar[0] = 0.; + pgpar[1] = 90.; + pgpar[2] = 5; + pgpar[3] = 2.; + pgpar[4] = -zsenv/2.; + pgpar[5] = 23.1; + pgpar[6] = 117.6; + pgpar[7] = zsenv/2.; + pgpar[8] = pgpar[5]; + pgpar[9] = pgpar[6]; + gMC->Gsvolu("SC3G1", "PGON", idGas, pgpar, 10); + gMC->Gsvolu("SC4G1", "PGON", idGas, pgpar, 10); + + //Segment-2 - vertical box ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + + bparV[0] = 1.0/2.; + bparV[1] = 95.5/2.; + bparV[2] = zsenv/2.; + gMC->Gsvolu("SC3G2", "BOX", idGas, bparV, 3); + gMC->Gsvolu("SC4G2", "BOX", idGas, bparV, 3); + + //........................................................................................... + + xposHorBox = bparH[0] + 23.1; + yposHorBox = -bparH[1]; + + xposVerBox = -bparV[0]; + yposVerBox = bparV[1] + 23.1; + + gMC->Gspos("SC3G0",1, "SQM3", xposHorBox,yposHorBox,0.,0, "ONLY"); + gMC->Gspos("SC3G1",1, "SQM3", 0.,0.,0.,0, "ONLY"); + gMC->Gspos("SC3G2",1, "SQM3", xposVerBox,yposVerBox,0.,0, "ONLY"); + + + gMC->Gspos("SC4G0",1, "SQM4", xposHorBox,yposHorBox,0.,0, "ONLY"); + gMC->Gspos("SC4G1",1, "SQM4", 0.,0.,0.,0, "ONLY"); + gMC->Gspos("SC4G2",1, "SQM4", xposVerBox,yposVerBox,0.,0, "ONLY"); + + + //^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ -//^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ //################################################################################################## @@ -780,10 +768,9 @@ void AliMUONSt2GeometryBuilderV2::CreateGeometry() } //______________________________________________________________________________ -void AliMUONSt2GeometryBuilderV2::SetTransformations() +void AliMUONSt2GeometryBuilderV2::SetVolumes() { -// Defines the transformations for the station2 chambers. -// --- +/// Defines the volumes for the station2 chambers. if (gAlice->GetModule("SHIL")) { SetMotherVolume(2, "YOUT1"); @@ -793,6 +780,12 @@ void AliMUONSt2GeometryBuilderV2::SetTransformations() // Define chamber volumes as virtual SetVolume(2, "SC03", true); SetVolume(3, "SC04", true); +} + +//______________________________________________________________________________ +void AliMUONSt2GeometryBuilderV2::SetTransformations() +{ +/// Defines the transformations for the station2 chambers. Double_t zpos1 = - AliMUONConstants::DefaultChamberZ(2); SetTranslation(2, TGeoTranslation(0., 0., zpos1)); @@ -804,15 +797,14 @@ void AliMUONSt2GeometryBuilderV2::SetTransformations() //______________________________________________________________________________ void AliMUONSt2GeometryBuilderV2::SetSensitiveVolumes() { -// Defines the sensitive volumes for station2 chambers. -// --- +/// Defines the sensitive volumes for station2 chambers. - GetGeometry(2)->SetSensitiveVolume("C3G0"); - GetGeometry(2)->SetSensitiveVolume("C3G1"); - GetGeometry(2)->SetSensitiveVolume("C3G2"); + GetGeometry(2)->SetSensitiveVolume("SC3G0"); + GetGeometry(2)->SetSensitiveVolume("SC3G1"); + GetGeometry(2)->SetSensitiveVolume("SC3G2"); - GetGeometry(3)->SetSensitiveVolume("C4G0"); - GetGeometry(3)->SetSensitiveVolume("C4G1"); - GetGeometry(3)->SetSensitiveVolume("C4G2"); + GetGeometry(3)->SetSensitiveVolume("SC4G0"); + GetGeometry(3)->SetSensitiveVolume("SC4G1"); + GetGeometry(3)->SetSensitiveVolume("SC4G2"); }