MUON/AliMUONSt2GeometryBuilderV2.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 // $Id$
  17
  18 //-----------------------------------------------------------------------------
  19 // Class AliMUONSt2GeometryBuilderV2
  20 // ---------------------------------
  21 // MUON Station2 coarse geometry construction class.
  22 // Author: SANJOY PAL ,Prof. SUKALYAN CHATTOPADHAYAY  [SINP, KOLKATA]
  23 //         &  Dr.SHAKEEL AHMAD (AMU), INDIA
  24 //-----------------------------------------------------------------------------
  25
  26
  27 #include <TVirtualMC.h>
  28 #include <TGeoMatrix.h>
  29 #include <Riostream.h>
  30
  31 #include "AliRun.h"
  32 #include "AliLog.h"
  33
  34 #include "AliMUONSt2GeometryBuilderV2.h"
  35 #include "AliMUON.h"
  36 #include "AliMUONConstants.h"
  37 #include "AliMUONGeometryModule.h"
  38 #include "AliMUONGeometryEnvelopeStore.h"
  39 #include "AliMUONConstants.h"
  40
  41 #define PI 3.14159
  42
  43 /// \cond CLASSIMP
  44 ClassImp(AliMUONSt2GeometryBuilderV2)
  45 /// \endcond
  46
  47 //______________________________________________________________________________
  48 AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2(AliMUON* muon)
  49  : AliMUONVGeometryBuilder(2, 2),
  50    fMUON(muon)
  51 {
  52 /// Standard constructor
  53
  54 }
  55
  56 //______________________________________________________________________________
  57 AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2()
  58  : AliMUONVGeometryBuilder(),
  59    fMUON(0)
  60 {
  61 /// Default constructor
  62 }
  63
  64 //______________________________________________________________________________
  65 AliMUONSt2GeometryBuilderV2::~AliMUONSt2GeometryBuilderV2()
  66 {
  67 /// Destructor
  68 }
  69
  70 //
  71 // public methods
  72 //
  73
  74 //______________________________________________________________________________
  75 void AliMUONSt2GeometryBuilderV2::CreateGeometry()
  76 {
  77 /// Geometry construction
  78
  79 //
  80 //********************************************************************
  81 //                            Station 2                             **
  82 //********************************************************************
  83      // indices 1 and 2 for first and second chambers in the station
  84      // iChamber (first chamber) kept for other quanties than Z,
  85      // assumed to be the same in both chambers
  86
  87      // Get tracking medias Ids
  88      Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
  89      Int_t idAir  = idtmed[1100]; // medium 1
  90      Int_t idGas  = idtmed[1108]; // medium Ar-CO2 gas (80%+20%)
  91      Int_t idPCB  = idtmed[1122]; // medium FR4
  92      Int_t idCU   = idtmed[1110]; // medium copper
  93      Int_t idRoha = idtmed[1113]; // medium roha cell
  94      Int_t idPGF30= idtmed[1123]; // medium for Frame Eq.to Bakelite
  95      //Int_t idScru = idtmed[1128]; // screw material - Stainless Steel(18%Cr,9%Ni,Fe)
  96
  97
  98
  99 /*########################################################################################
 100     Create volume for one Quadrant
 101 ##########################################################################################*/
 102      Float_t tpar1[5];
 103      tpar1[0] = 20.55;
 104      tpar1[1] = 123.5;
 105      tpar1[2] = 6.8/2;
 106      tpar1[3] = -12.0;
 107      tpar1[4] = 102.0;
 108
 109
 110      gMC->Gsvolu("SQM3","TUBS", idAir, tpar1, 5);
 111      gMC->Gsvolu("SQM4","TUBS", idAir, tpar1, 5);
 112
 113
 114 //==================================================================================
 115 //                                Plane
 116 //==================================================================================
 117
 118 //Thickness of variour parts
 119   Float_t zCu   = 0.005;      // eff. cu in cathode pcb
 120   Float_t zCbb  = 0.04;       //cathode pcb
 121   Float_t zRoha = 2.5;        // Rhocell
 122   Float_t zMeb = 0.04;      //Mech. exit board //0.08
 123   Float_t zEeb = 0.04;     //Effective electronic readout board //0.02
 124
 125   //starting Z-positions of various parts--- in Plane-1
 126   Float_t zposCu    = 0.25;   // 2.5 mm => gap between anode & chatode plane
 127   Float_t zposCbb   = zposCu + zCu;
 128   Float_t zposRoha  = zposCbb + zCbb;
 129   Float_t zposMeb   = zposRoha + zRoha;
 130   Float_t zposEeb   = zposMeb + zMeb;
 131
 132   Float_t zposCuBar   = zposCu +  zCu/2.; //for segment 0 & 6
 133   Float_t zposCbbBar  = zposCbb + zCbb/2.;
 134   Float_t zposRohaBar = zposRoha + zRoha/2.;
 135   Float_t zposMebBar  = zposMeb + zMeb/2.;
 136   Float_t zposEebBar  = zposEeb + zEeb/2.;
 137
 138
 139  //Cathode PCB + Copper sheet + Rohacell + mech exit board + eff. electronic exit board
 140
 141  //Segment-0 ~~~Horizantal box  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 142   Float_t bparH[3];
 143   bparH[0] = 95.5/2.; // extension beyond 0 deg in x direction // 94.5
 144   bparH[1] = 1.2/2.;  // extension beyond 0 deg in y direction (3.7[total extn] - 2.5[frame dim])
 145   bparH[2] = zCu/2.;     //Thickness of Copper sheet
 146   gMC->Gsvolu("SCU0L", "BOX", idCU, bparH, 3);
 147
 148   bparH[2] = zCbb/2.; // thickness of cathode sheet in z direction
 149   gMC->Gsvolu("SCB0L", "BOX", idPCB, bparH, 3);
 150
 151   bparH[2] = zRoha/2.;     //Thickness of Roha cell
 152   gMC->Gsvolu("SRH0L", "BOX", idRoha, bparH, 3);
 153
 154   bparH[2] = zMeb/2;       //Thickness of mechanical exit board
 155   gMC->Gsvolu("SMB0L", "BOX", idPCB, bparH, 3);
 156
 157   bparH[2] = zEeb/2;          //Thickness of effective electronic  exit board
 158   gMC->Gsvolu("SEB0L", "BOX", idCU, bparH, 3);
 159
 160   //Segment-1 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 161   Float_t pgpar[10]; // polygon
 162   pgpar[0] = 0.;  // initial angle
 163   pgpar[1] = 90.; // increment in angle starting from initial angle
 164   pgpar[2] = 5;   // number of side
 165   pgpar[3] = 2.; // number of plane
 166   pgpar[4] = -zCu/2.; // z-position of the first plane
 167   pgpar[5] = 23.1; // innner radius first plane
 168   pgpar[6] = 117.6;  // outer radious first plane
 169   pgpar[7] = zCu/2.; // z-position of the second plane
 170   pgpar[8] = pgpar[5];  // innner radius of second plane
 171   pgpar[9] = pgpar[6];  // outer radious of second plane
 172   gMC->Gsvolu("SCU1L", "PGON", idCU, pgpar, 10);
 173
 174   pgpar[4] = -zCbb/2.; // z-position of the first plane
 175   pgpar[7] = zCbb/2.;  // Thickness of copper-sheet
 176   gMC->Gsvolu("SCB1L", "PGON", idPCB, pgpar, 10);
 177
 178   pgpar[4] = -zRoha/2.;  // Thickness of Roha cell
 179   pgpar[7] = zRoha/2.;  // Thickness of Roha cell
 180   gMC->Gsvolu("SRH1L", "PGON", idRoha, pgpar, 10);
 181
 182   pgpar[4] = -zMeb/2.;  // Thickness of mechanical exit board
 183   pgpar[7] = zMeb/2.;  // Thickness of mechanical exit board
 184   gMC->Gsvolu("SMB1L", "PGON", idPCB, pgpar, 10);
 185
 186   pgpar[4] = -zEeb/2.;  // Thickness of electronic readout board
 187   pgpar[7] = zEeb/2.;  // Thickness of electronic readout board
 188   gMC->Gsvolu("SEB1L", "PGON", idCU, pgpar, 10);
 189
 190
 191 //Segment-2 - vertical box (simalar to horizontal bar as in Segment 0)~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 192
 193   Float_t bparV[3];
 194
 195   bparV[0] = 1.0/2.;
 196   bparV[1] = 95.5/2.; // 94.5
 197   bparV[2] = zCu/2.;
 198   gMC->Gsvolu("SCU2L", "BOX", idCU, bparV, 3);
 199
 200   bparV[2] = zCbb/2.;
 201   gMC->Gsvolu("SCB2L", "BOX", idPCB, bparV, 3);
 202
 203
 204   bparV[2] = zRoha/2.;
 205   gMC->Gsvolu("SRH2L", "BOX", idRoha, bparV, 3);
 206
 207   bparV[2] = zMeb/2;
 208   gMC->Gsvolu("SMB2L", "BOX", idPCB, bparV, 3);
 209
 210   bparV[2] = zEeb/2;
 211   gMC->Gsvolu("SEB2L", "BOX", idCU, bparV, 3);
 212
 213   //....(Setting posion of Segment 0,1,2)..................................................................
 214
 215   Float_t xposHorBox =  bparH[0] + 23.1; // 23.1 = 20.6(inner radius of qrd) + 2.5 (width of frame)
 216   Float_t yposHorBox = -bparH[1];
 217
 218   Float_t xposVerBox = -bparV[0];
 219   Float_t yposVerBox =  bparV[1] + 23.1;
 220
 221
 222   //Positioning the PCB
 223
 224
 225   // chamber 3
 226   gMC->Gspos("SCB0L",1, "SQM3",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY");
 227   gMC->Gspos("SCB1L",1, "SQM3", 0.0,0.0,zposCbbBar,0, "ONLY");
 228   gMC->Gspos("SCB2L",1, "SQM3",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY");
 229
 230   gMC->Gspos("SCB0L",2, "SQM3",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY");
 231   gMC->Gspos("SCB1L",2, "SQM3", 0.0,0.0,-zposCbbBar,0, "ONLY");
 232   gMC->Gspos("SCB2L",2, "SQM3",xposVerBox,yposVerBox,-zposCbbBar,0, "ONLY");
 233
 234   // chamber 4
 235   gMC->Gspos("SCB0L",3, "SQM4",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY");
 236   gMC->Gspos("SCB1L",3, "SQM4", 0.0,0.0,zposCbbBar,0, "ONLY");
 237   gMC->Gspos("SCB2L",3, "SQM4",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY");
 238
 239   gMC->Gspos("SCB0L",4, "SQM4",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY");
 240   gMC->Gspos("SCB1L",4, "SQM4", 0.0,0.0,-zposCbbBar,0, "ONLY");
 241   gMC->Gspos("SCB2L",4, "SQM4",xposVerBox,yposVerBox,-zposCbbBar,0, "ONLY");
 242
 243
 244   //Positioning Copper sheet
 245
 246   // chamber 3
 247
 248   gMC->Gspos("SCU0L",1, "SQM3",xposHorBox,yposHorBox,zposCuBar,0, "ONLY");
 249   gMC->Gspos("SCU1L",1, "SQM3", 0.0,0.0,zposCuBar,0, "ONLY");
 250   gMC->Gspos("SCU2L",1, "SQM3",xposVerBox,yposVerBox,zposCuBar,0, "ONLY");
 251
 252   gMC->Gspos("SCU0L",2, "SQM3",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY");
 253   gMC->Gspos("SCU1L",2, "SQM3", 0.0,0.0,-zposCuBar,0, "ONLY");
 254   gMC->Gspos("SCU2L",2, "SQM3",xposVerBox,yposVerBox,-zposCuBar,0, "ONLY");
 255
 256   // chamber 4
 257   gMC->Gspos("SCU0L",3, "SQM4",xposHorBox,yposHorBox,zposCuBar,0, "ONLY");
 258   gMC->Gspos("SCU1L",3, "SQM4", 0.0,0.0,zposCuBar,0, "ONLY");
 259   gMC->Gspos("SCU2L",3, "SQM4",xposVerBox,yposVerBox,zposCuBar,0, "ONLY");
 260
 261   gMC->Gspos("SCU0L",4, "SQM4",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY");
 262   gMC->Gspos("SCU1L",4, "SQM4", 0.0,0.0,-zposCuBar,0, "ONLY");
 263   gMC->Gspos("SCU2L",4, "SQM4",xposVerBox,yposVerBox,-zposCuBar,0, "ONLY");
 264
 265   //Positioning Roha cell
 266
 267   // chamber 3
 268
 269   gMC->Gspos("SRH0L",1, "SQM3",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY");
 270   gMC->Gspos("SRH1L",1, "SQM3", 0.0,0.0,zposRohaBar,0, "ONLY");
 271   gMC->Gspos("SRH2L",1, "SQM3",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY");
 272
 273   gMC->Gspos("SRH0L",2, "SQM3",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY");
 274   gMC->Gspos("SRH1L",2, "SQM3", 0.0,0.0,-zposRohaBar,0, "ONLY");
 275   gMC->Gspos("SRH2L",2, "SQM3",xposVerBox,yposVerBox,-zposRohaBar,0, "ONLY");
 276
 277   // chamber 4
 278
 279   gMC->Gspos("SRH0L",3, "SQM4",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY");
 280   gMC->Gspos("SRH1L",3, "SQM4", 0.0,0.0,zposRohaBar,0, "ONLY");
 281   gMC->Gspos("SRH2L",3, "SQM4",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY");
 282
 283   gMC->Gspos("SRH0L",4, "SQM4",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY");
 284   gMC->Gspos("SRH1L",4, "SQM4", 0.0,0.0,-zposRohaBar,0, "ONLY");
 285   gMC->Gspos("SRH2L",4, "SQM4",xposVerBox,yposVerBox,-zposRohaBar,0, "ONLY");
 286
 287   //Positioning of Mech. exit board
 288
 289   // chamber 3
 290
 291   gMC->Gspos("SMB0L",1, "SQM3",xposHorBox,yposHorBox,zposMebBar,0, "ONLY");
 292   gMC->Gspos("SMB1L",1, "SQM3", 0.0,0.0,zposMebBar,0, "ONLY");
 293   gMC->Gspos("SMB2L",1, "SQM3",xposVerBox,yposVerBox,zposMebBar,0, "ONLY");
 294
 295   gMC->Gspos("SMB0L",2, "SQM3",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY");
 296   gMC->Gspos("SMB1L",2, "SQM3", 0.0,0.0,-zposMebBar,0, "ONLY");
 297   gMC->Gspos("SMB2L",2, "SQM3",xposVerBox,yposVerBox,-zposMebBar,0, "ONLY");
 298   // chamber 4
 299
 300   gMC->Gspos("SMB0L",3, "SQM4",xposHorBox,yposHorBox,zposMebBar,0, "ONLY");
 301   gMC->Gspos("SMB1L",3, "SQM4", 0.0,0.0,zposMebBar,0, "ONLY");
 302   gMC->Gspos("SMB2L",3, "SQM4",xposVerBox,yposVerBox,zposMebBar,0, "ONLY");
 303
 304   gMC->Gspos("SMB0L",4, "SQM4",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY");
 305   gMC->Gspos("SMB1L",4, "SQM4", 0.0,0.0,-zposMebBar,0, "ONLY");
 306   gMC->Gspos("SMB2L",4, "SQM4",xposVerBox,yposVerBox,-zposMebBar,0, "ONLY");
 307
 308   //Positioning of Electronic exit board
 309
 310
 311   // chamber 3
 312
 313   gMC->Gspos("SEB0L",1, "SQM3",xposHorBox,yposHorBox,zposEebBar,0, "ONLY");
 314   gMC->Gspos("SEB1L",1, "SQM3", 0.0,0.0,zposEebBar,0, "ONLY");
 315   gMC->Gspos("SEB2L",1, "SQM3",xposVerBox,yposVerBox,zposEebBar,0, "ONLY");
 316
 317   gMC->Gspos("SEB0L",2, "SQM3",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY");
 318   gMC->Gspos("SEB1L",2, "SQM3", 0.0,0.0,-zposEebBar,0, "ONLY");
 319   gMC->Gspos("SEB2L",2, "SQM3",xposVerBox,yposVerBox,-zposEebBar,0, "ONLY");
 320   // chamber 4
 321
 322   gMC->Gspos("SEB0L",3, "SQM4",xposHorBox,yposHorBox,zposEebBar,0, "ONLY");
 323   gMC->Gspos("SEB1L",3, "SQM4", 0.0,0.0,zposEebBar,0, "ONLY");
 324   gMC->Gspos("SEB2L",3, "SQM4",xposVerBox,yposVerBox,zposEebBar,0, "ONLY");
 325
 326   gMC->Gspos("SEB0L",4, "SQM4",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY");
 327   gMC->Gspos("SEB1L",4, "SQM4", 0.0,0.0,-zposEebBar,0, "ONLY");
 328   gMC->Gspos("SEB2L",4, "SQM4",xposVerBox,yposVerBox,-zposEebBar,0, "ONLY");
 329
 330
 331   //----------------------------------------------------------------------
 332   //                         Frames
 333   //----------------------------------------------------------------------
 334   //Frame-1
 335
 336   Float_t frame1[3] ;
 337   frame1[0] = 101.0/2.;             //100.6 = 94.5 + 2.5 + 3.6
 338   frame1[1] = 2.5/2.;
 339   frame1[2] = 5.0/2.;
 340
 341   gMC->Gsvolu("SFRM1", "BOX", idPGF30, frame1, 3); //Frame - 1 // fill with pkk GF30
 342
 343   Float_t arib1[3];
 344   arib1[0] = frame1[0];
 345   arib1[1] = 0.9/2.;
 346   arib1[2] =(frame1[2]-0.95)/2.0;
 347
 348   gMC->Gsvolu("SFRA1", "BOX", idAir, arib1, 3); // fill with air
 349
 350   Float_t xposarib1 = 0;
 351   Float_t yposarib1 = -frame1[1] + arib1[1];
 352   Float_t zposarib1 = frame1[2] - arib1[2];
 353
 354   gMC->Gspos("SFRA1",1, "SFRM1", xposarib1, yposarib1, zposarib1,0, "ONLY");  //replace pkk GF30 with air(b)
 355   gMC->Gspos("SFRA1",2, "SFRM1", xposarib1, yposarib1, -zposarib1,0, "ONLY"); //replace pkk GF30 with air(nb)
 356
 357   Float_t rrib1[3];
 358   rrib1[0] = frame1[0];
 359   rrib1[1] = 0.6/2.;
 360   rrib1[2] =(frame1[2]-0.95)/2.0;
 361
 362   gMC->Gsvolu("SFRR1", "BOX", idRoha, rrib1, 3); // fill with rohacell
 363
 364   Float_t xposrrib1 = 0.0;
 365   Float_t yposrrib1 = frame1[1] - rrib1[1];
 366   Float_t zposrrib1 = frame1[2] - rrib1[2];
 367
 368   gMC->Gspos("SFRR1",1, "SFRM1", xposrrib1, yposrrib1, zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell
 369   gMC->Gspos("SFRR1",2, "SFRM1", xposrrib1, yposrrib1, -zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell
 370
 371
 372   Float_t xposFr1 = frame1[0] + 20.6;
 373   Float_t yposFr1 = -3.7 + frame1[1] ;
 374   Float_t zposFr1 = 0.0;
 375
 376   gMC->Gspos("SFRM1",1, "SQM3", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1
 377   gMC->Gspos("SFRM1",2, "SQM4", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1
 378
 379
 380   //......................................................................................
 381   //Frame-2
 382
 383   Float_t frame2[3];
 384   frame2[0] = 4.0/2.;
 385   frame2[1] = 1.2/2.;
 386   frame2[2] = 5.0/2;
 387
 388   gMC->Gsvolu("SFRM2", "BOX", idPGF30, frame2, 3); //Frame - 2
 389
 390   Float_t rrib2[3];
 391   rrib2[0] = frame2[0]-1.0/2.0;
 392   rrib2[1] = frame2[1];
 393   rrib2[2] =(frame2[2]-0.95)/2.0;
 394
 395   gMC->Gsvolu("SFRR2", "BOX", idRoha, rrib2, 3);
 396
 397   Float_t xposrrib2 = -1.0/2.0;
 398   Float_t yposrrib2 = 0.0;
 399   Float_t zposrrib2 = frame2[2] - rrib2[2];
 400
 401   gMC->Gspos("SFRR2",1, "SFRM2", xposrrib2, yposrrib2, zposrrib2,0, "ONLY");//replace pkk GF30 with rohacell
 402   gMC->Gspos("SFRR2",2, "SFRM2", xposrrib2, yposrrib2, -zposrrib2,0, "ONLY");//replace pkk GF30 with roha
 403
 404
 405
 406   Float_t xposFr2 = frame2[0] + 117.6;
 407   Float_t yposFr2 = -frame2[1];
 408   Float_t zposFr2 = 0.0;
 409
 410   gMC->Gspos("SFRM2",1, "SQM3", xposFr2, yposFr2, zposFr2,0, "MANY");//global positing of frame in SQM3
 411   gMC->Gspos("SFRM2",2, "SQM4", xposFr2, yposFr2, zposFr2,0, "MANY");//global positing of frame in SQM4
 412
 413
 414   //......................................................................................
 415
 416   //Frame-3
 417
 418   Float_t pgparFr3[10];
 419   pgparFr3[0] = 0.;
 420   pgparFr3[1] = 90.;
 421   pgparFr3[2] = 5;
 422   pgparFr3[3] = 2.;
 423   pgparFr3[4] = 0.;
 424   pgparFr3[5] = 117.6;
 425   pgparFr3[6] = 121.6;
 426   pgparFr3[7] = pgparFr3[4] + 5.0;
 427   pgparFr3[8] = pgparFr3[5];
 428   pgparFr3[9] = pgparFr3[6];
 429
 430   gMC->Gsvolu("SFRM3", "PGON", idPGF30, pgparFr3, 10);
 431
 432   Float_t pgparRrib3[10];
 433   pgparRrib3[0] = 0.;
 434   pgparRrib3[1] = 90.;
 435   pgparRrib3[2] = 5;
 436   pgparRrib3[3] = 2.;
 437   pgparRrib3[4] = 0.;
 438   pgparRrib3[5] = 117.6;
 439   pgparRrib3[6] = 120.6;
 440   pgparRrib3[7] = pgparRrib3[4] +1.55 ;
 441   pgparRrib3[8] = pgparRrib3[5];
 442   pgparRrib3[9] = pgparRrib3[6];
 443
 444   gMC->Gsvolu("SFRR3", "PGON", idRoha, pgparRrib3, 10);
 445
 446   Float_t xposrrib3 = 0.0;
 447   Float_t yposrrib3 = 0.0;
 448   Float_t zposrrib3 = 0.0;
 449
 450   gMC->Gspos("SFRR3",1, "SFRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY");
 451
 452   zposrrib3 = 3.45;
 453
 454   gMC->Gspos("SFRR3",2, "SFRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY");
 455
 456
 457
 458   Float_t xposFr3 = 0.0;
 459   Float_t yposFr3 = 0.0;
 460   Float_t zposFr3 = -frame1[2];
 461
 462   gMC->Gspos("SFRM3",1, "SQM3", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1
 463   gMC->Gspos("SFRM3",2, "SQM4", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1
 464
 465
 466   //......................................................................................
 467   //Frame-4
 468
 469   Float_t frame4[3];
 470   frame4[0] = 1.0/2.;
 471   frame4[1] = 4.0/2.;
 472   frame4[2] = frame1[2];
 473
 474   gMC->Gsvolu("SFRM4", "BOX", idPGF30, frame4, 3);
 475
 476   Float_t rrib4[3];
 477   rrib4[0] = frame4[0];
 478   rrib4[1] = frame4[1]-1.0/2;
 479   rrib4[2] =(frame4[2]-0.95)/2.0;
 480
 481   gMC->Gsvolu("SFRR4", "BOX", idRoha, rrib4, 3);
 482
 483   Float_t xposrrib4 = 0.0;
 484   Float_t yposrrib4 = -1.0/2;
 485   Float_t zposrrib4 = frame4[2] - rrib4[2];
 486
 487   gMC->Gspos("SFRR4",1, "SFRM4", xposrrib4, yposrrib4, zposrrib4,0, "ONLY");
 488   gMC->Gspos("SFRR4",2, "SFRM4", xposrrib4, yposrrib4, -zposrrib4,0, "ONLY");
 489
 490
 491
 492   Float_t xposFr4 = -frame4[0];
 493   Float_t yposFr4 = -frame4[1] + 117.6;
 494   Float_t zposFr4 = 0.0;
 495
 496   gMC->Gspos("SFRM4",1, "SQM3", xposFr4, yposFr4, zposFr4,0, "MANY");
 497   gMC->Gspos("SFRM4",2, "SQM4", xposFr4, yposFr4, zposFr4,0, "MANY");
 498
 499
 500   //......................................................................................
 501   //Frame-5
 502   Float_t frame5[3] ;
 503   frame5[0] = 2.7/2.;
 504   frame5[1] = 101.0/2.;
 505   frame5[2] = 5.0/2.;
 506
 507   gMC->Gsvolu("SFRM5", "BOX", idPGF30, frame5, 3); //Frame - 1
 508
 509   Float_t arib5[3];
 510   arib5[0] = 0.9/2.0;
 511   arib5[1] = frame5[1];
 512   arib5[2] = (frame5[2]-0.95)/2.0;
 513
 514   gMC->Gsvolu("SFRA5", "BOX", idAir, arib5, 3);
 515
 516   Float_t xposarib5 = -frame5[0] + arib5[0];
 517   Float_t yposarib5 = 0.0;
 518   Float_t zposarib5 = frame5[2] - arib5[2];
 519
 520   gMC->Gspos("SFRA5",1, "SFRM5", xposarib5, yposarib5, zposarib5,0, "ONLY");
 521   gMC->Gspos("SFRA5",2, "SFRM5", xposarib5, yposarib5, -zposarib5,0, "ONLY");
 522
 523   Float_t rrib5[3];
 524   rrib5[0] = 0.8/2.0;
 525   rrib5[1] = frame5[1];
 526   rrib5[2] = (frame5[2]-0.95)/2.0;
 527
 528   gMC->Gsvolu("SFRR5", "BOX", idRoha, rrib5, 3);
 529
 530   Float_t xposrrib5 = frame5[0] - rrib5[0];
 531   Float_t yposrrib5 = 0.0;
 532   Float_t zposrrib5 = frame5[2] - rrib5[2];
 533
 534   gMC->Gspos("SFRR5",1, "SFRM5", xposrrib5, yposrrib5, zposrrib5,0, "ONLY");
 535   gMC->Gspos("SFRR5",2, "SFRM5", xposrrib5, yposrrib5, -zposrrib5,0, "ONLY");
 536
 537   Float_t xposFr5 = -3.7 + frame5[0];
 538   Float_t yposFr5 = frame5[1] + 20.6;
 539   Float_t zposFr5 = 0.0;
 540
 541   gMC->Gspos("SFRM5",1, "SQM3", xposFr5, yposFr5, zposFr5,0, "ONLY");
 542   gMC->Gspos("SFRM5",2, "SQM4", xposFr5, yposFr5, zposFr5,0, "ONLY");
 543
 544   //......................................................................................
 545   //Frame -6
 546
 547   Float_t frame6[3];
 548   frame6[0] = 1.0/2.;
 549   frame6[1] = 2.5/2.;
 550   frame6[2] = frame1[2];
 551
 552   gMC->Gsvolu("SFRM6", "BOX", idPGF30, frame6, 3);
 553
 554   Float_t rrib6[3];
 555   rrib6[0] = frame6[0];
 556   rrib6[1] = 1.5/2.;
 557   rrib6[2] =(frame2[2]-0.95)/2.0;
 558
 559   gMC->Gsvolu("SFRR6", "BOX", idRoha, rrib6, 3);
 560
 561   Float_t xposrrib6 = 0.0;
 562   Float_t yposrrib6 = 1.0/2.0;
 563   Float_t zposrrib6 = frame6[2] - rrib6[2];
 564
 565   gMC->Gspos("SFRR6",1, "SFRM6", xposrrib6, yposrrib6, zposrrib6,0, "ONLY");
 566   gMC->Gspos("SFRR6",2, "SFRM6", xposrrib6, yposrrib6, -zposrrib6,0, "ONLY");
 567
 568
 569
 570   Float_t xposFr6 = -frame6[0];
 571   Float_t yposFr6 = frame6[1] + 20.6;
 572   Float_t zposFr6 = 0.0;
 573
 574   gMC->Gspos("SFRM6",1, "SQM3", xposFr6, yposFr6, zposFr6,0, "ONLY");
 575   gMC->Gspos("SFRM6",2, "SQM4", xposFr6, yposFr6, zposFr6,0, "ONLY");
 576
 577
 578 //......................................................................................
 579 //Frame - 7 inner pgon
 580
 581   Float_t pgparFr7[10];
 582   pgparFr7[0] = 0.;
 583   pgparFr7[1] = 90.;
 584   pgparFr7[2] = 5;
 585   pgparFr7[3] = 2.;
 586   pgparFr7[4] = 0.;
 587   pgparFr7[5] = 20.6;
 588   pgparFr7[6] = 23.1;
 589   pgparFr7[7] = pgparFr7[4] + 5.0;
 590   pgparFr7[8] = pgparFr7[5];
 591   pgparFr7[9] = pgparFr7[6];
 592
 593   gMC->Gsvolu("SFRM7", "PGON", idPGF30, pgparFr7, 10);
 594
 595   Float_t pgparRrib7[10];
 596   pgparRrib7[0] = 0.;
 597   pgparRrib7[1] = 90.;
 598   pgparRrib7[2] = 5;
 599   pgparRrib7[3] = 2.;
 600   pgparRrib7[4] = 0.;
 601   pgparRrib7[5] = 21.6;
 602   pgparRrib7[6] = 23.1;
 603   pgparRrib7[7] = pgparRrib7[4] + 1.55;
 604   pgparRrib7[8] = pgparRrib7[5];
 605   pgparRrib7[9] = pgparRrib7[6];
 606
 607   gMC->Gsvolu("SFRR7", "PGON", idRoha, pgparRrib7, 10);
 608
 609   Float_t xposrrib7 = 0.0;
 610   Float_t yposrrib7 = 0.0;
 611   Float_t zposrrib7 = 0.0;
 612
 613   gMC->Gspos("SFRR7",1, "SFRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY");
 614
 615   zposrrib7 = 3.45;
 616
 617   gMC->Gspos("SFRR7",2, "SFRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY");
 618
 619   Float_t xposFr7 = 0.0;
 620   Float_t yposFr7 = 0.0;
 621   Float_t zposFr7 = -frame1[2];
 622
 623   gMC->Gspos("SFRM7",1, "SQM3", xposFr7, yposFr7, zposFr7,0, "ONLY");
 624   gMC->Gspos("SFRM7",2, "SQM4", xposFr7, yposFr7, zposFr7,0, "ONLY");
 625
 626
 627   //......................................................................................
 628   //Frame - 8
 629
 630
 631   Float_t frame8[3] ;
 632   frame8[0] = 2.5/2.0;
 633   frame8[1] = 1.2/2.0;
 634   frame8[2] = frame1[2];
 635
 636   gMC->Gsvolu("SFRM8", "BOX", idPGF30, frame8, 3); //Frame - 2
 637
 638   Float_t rrib8[3];
 639   rrib8[0] = frame8[0]-1.0/2;
 640   rrib8[1] = frame8[1];
 641   rrib8[2] =(frame8[2]-0.95)/2.0;
 642
 643   gMC->Gsvolu("SFRR8", "BOX", idRoha, rrib8, 3);
 644
 645   Float_t xposrrib8 = -1.0/2;
 646   Float_t yposrrib8 = 0.0;
 647   Float_t zposrrib8 = frame8[2] - rrib8[2];
 648
 649   gMC->Gspos("SFRR8",1, "SFRM8", xposrrib8, yposrrib8, zposrrib8,0, "ONLY");
 650   gMC->Gspos("SFRR8",2, "SFRM8", xposrrib8, yposrrib8, -zposrrib8,0, "ONLY");
 651
 652
 653
 654   Float_t xposFr8 = frame8[0] + 20.6;
 655   Float_t yposFr8 = -frame8[1];
 656   Float_t zposFr8 = 0.0;
 657
 658   gMC->Gspos("SFRM8",1, "SQM3", xposFr8, yposFr8, zposFr8,0, "ONLY");
 659   gMC->Gspos("SFRM8",2, "SQM4", xposFr8, yposFr8, zposFr8,0, "ONLY");
 660
 661
 662
 663
 664   //^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 665
 666   Float_t zsenv = 0.5; // distance between two cathode plane
 667
 668   //Segment-0 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 669   bparH[0] = 94.5/2.;
 670   bparH[1] = 1.2/2.;
 671   bparH[2] = zsenv/2.;
 672   gMC->Gsvolu("SC3G0", "BOX", idGas, bparH, 3);
 673   gMC->Gsvolu("SC4G0", "BOX", idGas, bparH, 3);
 674
 675  //Segment-1 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 676   pgpar[0] = 0.;
 677   pgpar[1] = 90.;
 678   pgpar[2] = 5;
 679   pgpar[3] = 2.;
 680   pgpar[4] = -zsenv/2.;
 681   pgpar[5] = 23.1;
 682   pgpar[6] = 117.6;
 683   pgpar[7] = zsenv/2.;
 684   pgpar[8] = pgpar[5];
 685   pgpar[9] = pgpar[6];
 686   gMC->Gsvolu("SC3G1", "PGON", idGas, pgpar, 10);
 687   gMC->Gsvolu("SC4G1", "PGON", idGas, pgpar, 10);
 688
 689
 690   //Segment-2 - vertical box ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 691
 692   bparV[0] = 1.0/2.;
 693   bparV[1] = 95.5/2.;
 694   bparV[2] = zsenv/2.;
 695   gMC->Gsvolu("SC3G2", "BOX", idGas, bparV, 3);
 696   gMC->Gsvolu("SC4G2", "BOX", idGas, bparV, 3);
 697
 698   //...........................................................................................
 699
 700   xposHorBox =  bparH[0] + 23.1;
 701   yposHorBox = -bparH[1];
 702
 703   xposVerBox = -bparV[0];
 704   yposVerBox = bparV[1] + 23.1;
 705
 706   gMC->Gspos("SC3G0",1, "SQM3", xposHorBox,yposHorBox,0.,0, "ONLY");
 707   gMC->Gspos("SC3G1",1, "SQM3", 0.,0.,0.,0, "ONLY");
 708   gMC->Gspos("SC3G2",1, "SQM3", xposVerBox,yposVerBox,0.,0, "ONLY");
 709
 710
 711   gMC->Gspos("SC4G0",1, "SQM4", xposHorBox,yposHorBox,0.,0, "ONLY");
 712   gMC->Gspos("SC4G1",1, "SQM4", 0.,0.,0.,0, "ONLY");
 713   gMC->Gspos("SC4G2",1, "SQM4", xposVerBox,yposVerBox,0.,0, "ONLY");
 714
 715
 716   //^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 717
 718
 719
 720 //##################################################################################################
 721 //   Positioning Quadrant  in chamber#3 and chamber#4
 722 //##################################################################################################
 723 /******Transformations for  Quadrant**********************************************
 724       ||        I  => Quadrant I:   no rotation
 725       ||
 726   II. || I.    II  => Quadrant II:  Reflaction of Quadrant I in XZ plane
 727       ||           => TGeoRotation("Qrot2",90.,180.,90.,90.,180.,0.);
 728 =============
 729       ||       III => Quadrant III: 180 degree rotation of Quadrant I in XY plane
 730  III. || IV.       => TGeoRotation("Qrot3",90.,180.,90.,90.,180.,0.);
 731       ||        IV => Quadrant IV:-180 degree rotation of Quadrant II in XY plane
 732                    => TGeoRotation("Qrot4",90.,0.,90.,-90.,180.,0.);
 733 **********************************************************************************************/
 734
 735  Int_t detElemId1 =  1;  // quadrant I
 736  Int_t detElemId2 =  0;  // quadrant II
 737  Int_t detElemId3 =  3;  // quadrant III
 738  Int_t detElemId4 =  2;  // quadrant IV
 739
 740  //Float_t halfChamber = zCbb + zCu + zRoha + zMeb + zEeb + zsenv/2;
 741 //   cout<<  "\n half_chamber \t" << half_chamber << endl;
 742
 743  Float_t halfChamber = tpar1[2];// tpar1[2] = 6.8/2;
 744
 745 // ------------------------------St2 Chamber3------------------------------------------------
 746
 747  //    GetEnvelopes(2)->AddEnvelope("S3M0", 300, true,TGeoTranslation(0.,0.,0.));
 748     GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId1, 1, TGeoTranslation( 0., 0., - halfChamber));
 749     GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId2, 2, TGeoTranslation( 0., 0., + halfChamber),
 750                                  TGeoRotation("Qrot3",90.,180.,90.,90.,180.,0.));
 751     GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId3, 3, TGeoTranslation( 0., 0., - halfChamber),
 752                                   TGeoRotation("Qrot3",90.,180.,90.,270.,0.,0.));
 753     GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId4, 4, TGeoTranslation( 0., 0., + halfChamber),
 754                                   TGeoRotation("Qrot4",90.,0.,90.,-90.,180.,0.));
 755
 756 //--------------------------------St2 Chamber4-------------------------------------------------
 757
 758     GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId1, 1, TGeoTranslation( 0., 0., - halfChamber));
 759     GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId2, 2, TGeoTranslation( 0., 0., + halfChamber),
 760                                  TGeoRotation("Qrot2",90.,180.,90.,90.,180.,0.));
 761     GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId3, 3, TGeoTranslation( 0., 0., - halfChamber),
 762                                   TGeoRotation("Qrot3",90.,180.,90.,270.,0.,0.));
 763     GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId4, 4, TGeoTranslation( 0., 0., + halfChamber),
 764                                   TGeoRotation("Qrot4",90.,0.,90.,-90.,180.,0.));
 765
 766 //**********************************************************************************************
 767
 768 }
 769 //______________________________________________________________________________
 770 void AliMUONSt2GeometryBuilderV2::SetVolumes()
 771 {
 772 /// Defines the volumes for the station2 chambers.
 773
 774   if (gAlice->GetModule("SHIL")) {
 775     SetMotherVolume(2, "YOUT1");
 776     SetMotherVolume(3, "YOUT1");
 777   }
 778
 779   // Define chamber volumes as virtual
 780   SetVolume(2, "SC03", true);
 781   SetVolume(3, "SC04", true);
 782 }
 783
 784 //______________________________________________________________________________
 785 void AliMUONSt2GeometryBuilderV2::SetTransformations()
 786 {
 787 /// Defines the transformations for the station2 chambers.
 788
 789   Double_t zpos1 = - AliMUONConstants::DefaultChamberZ(2);
 790   SetTranslation(2, TGeoTranslation(0., 0., zpos1));
 791
 792   Double_t zpos2 = - AliMUONConstants::DefaultChamberZ(3);
 793   SetTranslation(3, TGeoTranslation(0., 0., zpos2));
 794 }
 795
 796 //______________________________________________________________________________
 797 void AliMUONSt2GeometryBuilderV2::SetSensitiveVolumes()
 798 {
 799 /// Defines the sensitive volumes for station2 chambers.
 800
 801   GetGeometry(2)->SetSensitiveVolume("SC3G0");
 802   GetGeometry(2)->SetSensitiveVolume("SC3G1");
 803   GetGeometry(2)->SetSensitiveVolume("SC3G2");
 804
 805   GetGeometry(3)->SetSensitiveVolume("SC4G0");
 806   GetGeometry(3)->SetSensitiveVolume("SC4G1");
 807   GetGeometry(3)->SetSensitiveVolume("SC4G2");
 808 }
 809