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