MUON/AliMUONv1.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 /*
  17 $Log$
  18 Revision 1.3  2000/06/22 14:10:05  morsch
  19 HP scope problems corrected (PH)
  20
  21 Revision 1.2  2000/06/15 07:58:49  morsch
  22 Code from MUON-dev joined
  23
  24 Revision 1.1.2.14  2000/06/14 14:37:25  morsch
  25 Initialization of TriggerCircuit added (PC)
  26
  27 Revision 1.1.2.13  2000/06/09 21:55:47  morsch
  28 Most coding rule violations corrected.
  29
  30 Revision 1.1.2.12  2000/05/05 11:34:29  morsch
  31 Log inside comments.
  32
  33 Revision 1.1.2.11  2000/05/05 10:06:48  morsch
  34 Coding Rule violations regarding trigger section corrected (CP)
  35 Log messages included.
  36 */
  37
  38 /////////////////////////////////////////////////////////
  39 //  Manager and hits classes for set:MUON version 0    //
  40 /////////////////////////////////////////////////////////
  41
  42 #include <TTUBE.h>
  43 #include <TNode.h>
  44 #include <TRandom.h>
  45 #include <TLorentzVector.h>
  46 #include <iostream.h>
  47
  48 #include "AliMUONv1.h"
  49 #include "AliRun.h"
  50 #include "AliMC.h"
  51 #include "AliCallf77.h"
  52 #include "AliConst.h"
  53 #include "AliMUONChamber.h"
  54 #include "AliMUONHit.h"
  55 #include "AliMUONPadHit.h"
  56 #include "AliMUONConstants.h"
  57
  58 ClassImp(AliMUONv1)
  59
  60 //___________________________________________
  61 AliMUONv1::AliMUONv1() : AliMUON()
  62 {
  63 // Constructor
  64     fChambers = 0;
  65 }
  66
  67 //___________________________________________
  68 AliMUONv1::AliMUONv1(const char *name, const char *title)
  69        : AliMUON(name,title)
  70 {
  71 // Constructor
  72 }
  73
  74 //___________________________________________
  75 void AliMUONv1::CreateGeometry()
  76 {
  77 //
  78 //   Note: all chambers have the same structure, which could be
  79 //   easily parameterised. This was intentionally not done in order
  80 //   to give a starting point for the implementation of the actual
  81 //   design of each station.
  82   Int_t *idtmed = fIdtmed->GetArray()-1099;
  83
  84 //   Distance between Stations
  85 //
  86      Float_t bpar[3];
  87      Float_t tpar[3];
  88      Float_t pgpar[10];
  89      Float_t zpos1, zpos2, zfpos;
  90      Float_t dframep=.001; // Value for station 3 should be 6 ...
  91      Float_t dframep1=.001;
  92 //     Bool_t frames=kTRUE;
  93      Bool_t frames=kFALSE;
  94
  95      Float_t dframez=0.9;
  96      Float_t dr;
  97      Float_t dstation;
  98
  99 //
 100 //   Rotation matrices in the x-y plane
 101      Int_t idrotm[1199];
 102 //   phi=   0 deg
 103      AliMatrix(idrotm[1100],  90.,   0., 90.,  90., 0., 0.);
 104 //   phi=  90 deg
 105      AliMatrix(idrotm[1101],  90.,  90., 90., 180., 0., 0.);
 106 //   phi= 180 deg
 107      AliMatrix(idrotm[1102],  90., 180., 90., 270., 0., 0.);
 108 //   phi= 270 deg
 109      AliMatrix(idrotm[1103],  90., 270., 90.,   0., 0., 0.);
 110 //
 111      Float_t phi=2*TMath::Pi()/12/2;
 112
 113 //
 114 //   pointer to the current chamber
 115 //   pointer to the current chamber
 116      Int_t idAlu1=idtmed[1103];
 117      Int_t idAlu2=idtmed[1104];
 118 //     Int_t idAlu1=idtmed[1100];
 119 //     Int_t idAlu2=idtmed[1100];
 120      Int_t idAir=idtmed[1100];
 121      Int_t idGas=idtmed[1105];
 122
 123
 124      AliMUONChamber *iChamber, *iChamber1, *iChamber2;
 125 //********************************************************************
 126 //                            Station 1                             **
 127 //********************************************************************
 128 //  CONCENTRIC
 129      // indices 1 and 2 for first and second chambers in the station
 130      // iChamber (first chamber) kept for other quanties than Z,
 131      // assumed to be the same in both chambers
 132      iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[0];
 133      iChamber2 =(AliMUONChamber*) (*fChambers)[1];
 134      zpos1=iChamber1->Z();
 135      zpos2=iChamber2->Z();
 136      dstation = zpos2 - zpos1;
 137      zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
 138
 139 //
 140 //   Mother volume
 141      tpar[0] = iChamber->RInner()-dframep1;
 142      tpar[1] = (iChamber->ROuter()+dframep1)/TMath::Cos(phi);
 143      tpar[2] = dstation/4;
 144
 145      gMC->Gsvolu("C01M", "TUBE", idAir, tpar, 3);
 146      gMC->Gsvolu("C02M", "TUBE", idAir, tpar, 3);
 147      gMC->Gspos("C01M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
 148      gMC->Gspos("C02M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");
 149 // Aluminium frames
 150 // Outer frames
 151      pgpar[0] = 360/12/2;
 152      pgpar[1] = 360.;
 153      pgpar[2] = 12.;
 154      pgpar[3] =   2;
 155      pgpar[4] = -dframez/2;
 156      pgpar[5] = iChamber->ROuter();
 157      pgpar[6] = pgpar[5]+dframep1;
 158      pgpar[7] = +dframez/2;
 159      pgpar[8] = pgpar[5];
 160      pgpar[9] = pgpar[6];
 161      gMC->Gsvolu("C01O", "PGON", idAlu1, pgpar, 10);
 162      gMC->Gsvolu("C02O", "PGON", idAlu1, pgpar, 10);
 163      gMC->Gspos("C01O",1,"C01M", 0.,0.,-zfpos,  0,"ONLY");
 164      gMC->Gspos("C01O",2,"C01M", 0.,0.,+zfpos,  0,"ONLY");
 165      gMC->Gspos("C02O",1,"C02M", 0.,0.,-zfpos,  0,"ONLY");
 166      gMC->Gspos("C02O",2,"C02M", 0.,0.,+zfpos,  0,"ONLY");
 167 //
 168 // Inner frame
 169      tpar[0]= iChamber->RInner()-dframep1;
 170      tpar[1]= iChamber->RInner();
 171      tpar[2]= dframez/2;
 172      gMC->Gsvolu("C01I", "TUBE", idAlu1, tpar, 3);
 173      gMC->Gsvolu("C02I", "TUBE", idAlu1, tpar, 3);
 174
 175      gMC->Gspos("C01I",1,"C01M", 0.,0.,-zfpos,  0,"ONLY");
 176      gMC->Gspos("C01I",2,"C01M", 0.,0.,+zfpos,  0,"ONLY");
 177      gMC->Gspos("C02I",1,"C02M", 0.,0.,-zfpos,  0,"ONLY");
 178      gMC->Gspos("C02I",2,"C02M", 0.,0.,+zfpos,  0,"ONLY");
 179 //
 180 // Frame Crosses
 181      if (frames) {
 182
 183          bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2;
 184          bpar[1] = dframep1/2;
 185          bpar[2] = dframez/2;
 186          gMC->Gsvolu("C01B", "BOX", idAlu1, bpar, 3);
 187          gMC->Gsvolu("C02B", "BOX", idAlu1, bpar, 3);
 188
 189          gMC->Gspos("C01B",1,"C01M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 190                     idrotm[1100],"ONLY");
 191          gMC->Gspos("C01B",2,"C01M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 192                     idrotm[1100],"ONLY");
 193          gMC->Gspos("C01B",3,"C01M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 194                     idrotm[1101],"ONLY");
 195          gMC->Gspos("C01B",4,"C01M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 196                     idrotm[1101],"ONLY");
 197          gMC->Gspos("C01B",5,"C01M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 198                     idrotm[1100],"ONLY");
 199          gMC->Gspos("C01B",6,"C01M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 200                     idrotm[1100],"ONLY");
 201          gMC->Gspos("C01B",7,"C01M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 202                     idrotm[1101],"ONLY");
 203          gMC->Gspos("C01B",8,"C01M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 204                     idrotm[1101],"ONLY");
 205
 206          gMC->Gspos("C02B",1,"C02M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 207                     idrotm[1100],"ONLY");
 208          gMC->Gspos("C02B",2,"C02M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 209                     idrotm[1100],"ONLY");
 210          gMC->Gspos("C02B",3,"C02M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 211                     idrotm[1101],"ONLY");
 212          gMC->Gspos("C02B",4,"C02M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 213                     idrotm[1101],"ONLY");
 214          gMC->Gspos("C02B",5,"C02M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 215                     idrotm[1100],"ONLY");
 216          gMC->Gspos("C02B",6,"C02M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 217                     idrotm[1100],"ONLY");
 218          gMC->Gspos("C02B",7,"C02M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 219                     idrotm[1101],"ONLY");
 220          gMC->Gspos("C02B",8,"C02M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 221                     idrotm[1101],"ONLY");
 222      }
 223 //
 224 //   Chamber Material represented by Alu sheet
 225      tpar[0]= iChamber->RInner();
 226      tpar[1]= iChamber->ROuter();
 227      tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2;
 228      gMC->Gsvolu("C01A", "TUBE",  idAlu2, tpar, 3);
 229      gMC->Gsvolu("C02A", "TUBE",idAlu2, tpar, 3);
 230      gMC->Gspos("C01A", 1, "C01M", 0., 0., 0.,  0, "ONLY");
 231      gMC->Gspos("C02A", 1, "C02M", 0., 0., 0.,  0, "ONLY");
 232 //
 233 //   Sensitive volumes
 234      // tpar[2] = iChamber->DGas();
 235      tpar[2] = iChamber->DGas()/2;
 236      gMC->Gsvolu("C01G", "TUBE", idtmed[1108], tpar, 3);
 237      gMC->Gsvolu("C02G", "TUBE", idtmed[1108], tpar, 3);
 238      gMC->Gspos("C01G", 1, "C01A", 0., 0., 0.,  0, "ONLY");
 239      gMC->Gspos("C02G", 1, "C02A", 0., 0., 0.,  0, "ONLY");
 240 //
 241 // Frame Crosses to be placed inside gas
 242      if (frames) {
 243
 244          dr = (iChamber->ROuter() - iChamber->RInner());
 245          bpar[0] = TMath::Sqrt(dr*dr-dframep1*dframep1/4)/2;
 246          bpar[1] = dframep1/2;
 247          bpar[2] = iChamber->DGas()/2;
 248          gMC->Gsvolu("C01F", "BOX", idAlu1, bpar, 3);
 249          gMC->Gsvolu("C02F", "BOX", idAlu1, bpar, 3);
 250
 251          gMC->Gspos("C01F",1,"C01G", +iChamber->RInner()+bpar[0] , 0, 0,
 252                     idrotm[1100],"ONLY");
 253          gMC->Gspos("C01F",2,"C01G", -iChamber->RInner()-bpar[0] , 0, 0,
 254                     idrotm[1100],"ONLY");
 255          gMC->Gspos("C01F",3,"C01G", 0, +iChamber->RInner()+bpar[0] , 0,
 256                     idrotm[1101],"ONLY");
 257          gMC->Gspos("C01F",4,"C01G", 0, -iChamber->RInner()-bpar[0] , 0,
 258                     idrotm[1101],"ONLY");
 259
 260          gMC->Gspos("C02F",1,"C02G", +iChamber->RInner()+bpar[0] , 0, 0,
 261                     idrotm[1100],"ONLY");
 262          gMC->Gspos("C02F",2,"C02G", -iChamber->RInner()-bpar[0] , 0, 0,
 263                     idrotm[1100],"ONLY");
 264          gMC->Gspos("C02F",3,"C02G", 0, +iChamber->RInner()+bpar[0] , 0,
 265                     idrotm[1101],"ONLY");
 266          gMC->Gspos("C02F",4,"C02G", 0, -iChamber->RInner()-bpar[0] , 0,
 267                     idrotm[1101],"ONLY");
 268      }
 269
 270 //
 271 //
 272 //********************************************************************
 273 //                            Station 2                             **
 274 //********************************************************************
 275      // indices 1 and 2 for first and second chambers in the station
 276      // iChamber (first chamber) kept for other quanties than Z,
 277      // assumed to be the same in both chambers
 278      iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[2];
 279      iChamber2 =(AliMUONChamber*) (*fChambers)[3];
 280      zpos1=iChamber1->Z();
 281      zpos2=iChamber2->Z();
 282      dstation = zpos2 - zpos1;
 283      zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
 284
 285 //
 286 //   Mother volume
 287      tpar[0] = iChamber->RInner()-dframep;
 288      tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi);
 289      tpar[2] = dstation/4;
 290
 291      gMC->Gsvolu("C03M", "TUBE", idAir, tpar, 3);
 292      gMC->Gsvolu("C04M", "TUBE", idAir, tpar, 3);
 293      gMC->Gspos("C03M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
 294      gMC->Gspos("C04M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");
 295 // Aluminium frames
 296 // Outer frames
 297      pgpar[0] = 360/12/2;
 298      pgpar[1] = 360.;
 299      pgpar[2] = 12.;
 300      pgpar[3] =   2;
 301      pgpar[4] = -dframez/2;
 302      pgpar[5] = iChamber->ROuter();
 303      pgpar[6] = pgpar[5]+dframep;
 304      pgpar[7] = +dframez/2;
 305      pgpar[8] = pgpar[5];
 306      pgpar[9] = pgpar[6];
 307      gMC->Gsvolu("C03O", "PGON", idAlu1, pgpar, 10);
 308      gMC->Gsvolu("C04O", "PGON", idAlu1, pgpar, 10);
 309      gMC->Gspos("C03O",1,"C03M", 0.,0.,-zfpos,  0,"ONLY");
 310      gMC->Gspos("C03O",2,"C03M", 0.,0.,+zfpos,  0,"ONLY");
 311      gMC->Gspos("C04O",1,"C04M", 0.,0.,-zfpos,  0,"ONLY");
 312      gMC->Gspos("C04O",2,"C04M", 0.,0.,+zfpos,  0,"ONLY");
 313 //
 314 // Inner frame
 315      tpar[0]= iChamber->RInner()-dframep;
 316      tpar[1]= iChamber->RInner();
 317      tpar[2]= dframez/2;
 318      gMC->Gsvolu("C03I", "TUBE", idAlu1, tpar, 3);
 319      gMC->Gsvolu("C04I", "TUBE", idAlu1, tpar, 3);
 320
 321      gMC->Gspos("C03I",1,"C03M", 0.,0.,-zfpos,  0,"ONLY");
 322      gMC->Gspos("C03I",2,"C03M", 0.,0.,+zfpos,  0,"ONLY");
 323      gMC->Gspos("C04I",1,"C04M", 0.,0.,-zfpos,  0,"ONLY");
 324      gMC->Gspos("C04I",2,"C04M", 0.,0.,+zfpos,  0,"ONLY");
 325 //
 326 // Frame Crosses
 327      if (frames) {
 328
 329          bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2;
 330          bpar[1] = dframep/2;
 331          bpar[2] = dframez/2;
 332          gMC->Gsvolu("C03B", "BOX", idAlu1, bpar, 3);
 333          gMC->Gsvolu("C04B", "BOX", idAlu1, bpar, 3);
 334
 335          gMC->Gspos("C03B",1,"C03M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 336                     idrotm[1100],"ONLY");
 337          gMC->Gspos("C03B",2,"C03M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 338                     idrotm[1100],"ONLY");
 339          gMC->Gspos("C03B",3,"C03M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 340                     idrotm[1101],"ONLY");
 341          gMC->Gspos("C03B",4,"C03M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 342                     idrotm[1101],"ONLY");
 343          gMC->Gspos("C03B",5,"C03M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 344                     idrotm[1100],"ONLY");
 345          gMC->Gspos("C03B",6,"C03M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 346                     idrotm[1100],"ONLY");
 347          gMC->Gspos("C03B",7,"C03M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 348                     idrotm[1101],"ONLY");
 349          gMC->Gspos("C03B",8,"C03M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 350                     idrotm[1101],"ONLY");
 351
 352          gMC->Gspos("C04B",1,"C04M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 353                     idrotm[1100],"ONLY");
 354          gMC->Gspos("C04B",2,"C04M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 355                     idrotm[1100],"ONLY");
 356          gMC->Gspos("C04B",3,"C04M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 357                     idrotm[1101],"ONLY");
 358          gMC->Gspos("C04B",4,"C04M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 359                     idrotm[1101],"ONLY");
 360          gMC->Gspos("C04B",5,"C04M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 361                     idrotm[1100],"ONLY");
 362          gMC->Gspos("C04B",6,"C04M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 363                     idrotm[1100],"ONLY");
 364          gMC->Gspos("C04B",7,"C04M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 365                     idrotm[1101],"ONLY");
 366          gMC->Gspos("C04B",8,"C04M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 367                     idrotm[1101],"ONLY");
 368      }
 369 //
 370 //   Chamber Material represented by Alu sheet
 371      tpar[0]= iChamber->RInner();
 372      tpar[1]= iChamber->ROuter();
 373      tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2;
 374      gMC->Gsvolu("C03A", "TUBE", idAlu2, tpar, 3);
 375      gMC->Gsvolu("C04A", "TUBE", idAlu2, tpar, 3);
 376      gMC->Gspos("C03A", 1, "C03M", 0., 0., 0.,  0, "ONLY");
 377      gMC->Gspos("C04A", 1, "C04M", 0., 0., 0.,  0, "ONLY");
 378 //
 379 //   Sensitive volumes
 380      // tpar[2] = iChamber->DGas();
 381      tpar[2] = iChamber->DGas()/2;
 382      gMC->Gsvolu("C03G", "TUBE", idGas, tpar, 3);
 383      gMC->Gsvolu("C04G", "TUBE", idGas, tpar, 3);
 384      gMC->Gspos("C03G", 1, "C03A", 0., 0., 0.,  0, "ONLY");
 385      gMC->Gspos("C04G", 1, "C04A", 0., 0., 0.,  0, "ONLY");
 386
 387      if (frames) {
 388 //
 389 // Frame Crosses to be placed inside gas
 390          dr = (iChamber->ROuter() - iChamber->RInner());
 391          bpar[0] = TMath::Sqrt(dr*dr-dframep*dframep/4)/2;
 392          bpar[1] = dframep/2;
 393          bpar[2] = iChamber->DGas()/2;
 394          gMC->Gsvolu("C03F", "BOX", idAlu1, bpar, 3);
 395          gMC->Gsvolu("C04F", "BOX", idAlu1, bpar, 3);
 396
 397          gMC->Gspos("C03F",1,"C03G", +iChamber->RInner()+bpar[0] , 0, 0,
 398                     idrotm[1100],"ONLY");
 399          gMC->Gspos("C03F",2,"C03G", -iChamber->RInner()-bpar[0] , 0, 0,
 400                     idrotm[1100],"ONLY");
 401          gMC->Gspos("C03F",3,"C03G", 0, +iChamber->RInner()+bpar[0] , 0,
 402                     idrotm[1101],"ONLY");
 403          gMC->Gspos("C03F",4,"C03G", 0, -iChamber->RInner()-bpar[0] , 0,
 404                     idrotm[1101],"ONLY");
 405
 406          gMC->Gspos("C04F",1,"C04G", +iChamber->RInner()+bpar[0] , 0, 0,
 407                     idrotm[1100],"ONLY");
 408          gMC->Gspos("C04F",2,"C04G", -iChamber->RInner()-bpar[0] , 0, 0,
 409                     idrotm[1100],"ONLY");
 410          gMC->Gspos("C04F",3,"C04G", 0, +iChamber->RInner()+bpar[0] , 0,
 411                     idrotm[1101],"ONLY");
 412          gMC->Gspos("C04F",4,"C04G", 0, -iChamber->RInner()-bpar[0] , 0,
 413                     idrotm[1101],"ONLY");
 414      }
 415
 416 //********************************************************************
 417 //                            Station 3                             **
 418 //********************************************************************
 419      // indices 1 and 2 for first and second chambers in the station
 420      // iChamber (first chamber) kept for other quanties than Z,
 421      // assumed to be the same in both chambers
 422      iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[4];
 423      iChamber2 =(AliMUONChamber*) (*fChambers)[5];
 424      zpos1=iChamber1->Z();
 425      zpos2=iChamber2->Z();
 426      dstation = zpos2 - zpos1;
 427
 428      zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
 429 //
 430 //   Mother volume
 431      tpar[0] = iChamber->RInner()-dframep;
 432      tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi);
 433      tpar[2] = dstation/4;
 434      gMC->Gsvolu("C05M", "TUBE", idAir, tpar, 3);
 435      gMC->Gsvolu("C06M", "TUBE", idAir, tpar, 3);
 436      gMC->Gspos("C05M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
 437      gMC->Gspos("C06M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");
 438 // Aluminium frames
 439 // Outer frames
 440      pgpar[0] = 360/12/2;
 441      pgpar[1] = 360.;
 442      pgpar[2] = 12.;
 443      pgpar[3] =   2;
 444      pgpar[4] = -dframez/2;
 445      pgpar[5] = iChamber->ROuter();
 446      pgpar[6] = pgpar[5]+dframep;
 447      pgpar[7] = +dframez/2;
 448      pgpar[8] = pgpar[5];
 449      pgpar[9] = pgpar[6];
 450      gMC->Gsvolu("C05O", "PGON", idAlu1, pgpar, 10);
 451      gMC->Gsvolu("C06O", "PGON", idAlu1, pgpar, 10);
 452      gMC->Gspos("C05O",1,"C05M", 0.,0.,-zfpos,  0,"ONLY");
 453      gMC->Gspos("C05O",2,"C05M", 0.,0.,+zfpos,  0,"ONLY");
 454      gMC->Gspos("C06O",1,"C06M", 0.,0.,-zfpos,  0,"ONLY");
 455      gMC->Gspos("C06O",2,"C06M", 0.,0.,+zfpos,  0,"ONLY");
 456 //
 457 // Inner frame
 458      tpar[0]= iChamber->RInner()-dframep;
 459      tpar[1]= iChamber->RInner();
 460      tpar[2]= dframez/2;
 461      gMC->Gsvolu("C05I", "TUBE", idAlu1, tpar, 3);
 462      gMC->Gsvolu("C06I", "TUBE", idAlu1, tpar, 3);
 463
 464      gMC->Gspos("C05I",1,"C05M", 0.,0.,-zfpos,  0,"ONLY");
 465      gMC->Gspos("C05I",2,"C05M", 0.,0.,+zfpos,  0,"ONLY");
 466      gMC->Gspos("C06I",1,"C06M", 0.,0.,-zfpos,  0,"ONLY");
 467      gMC->Gspos("C06I",2,"C06M", 0.,0.,+zfpos,  0,"ONLY");
 468 //
 469 // Frame Crosses
 470      if (frames) {
 471          bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2;
 472          bpar[1] = dframep/2;
 473          bpar[2] = dframez/2;
 474          gMC->Gsvolu("C05B", "BOX", idAlu1, bpar, 3);
 475          gMC->Gsvolu("C06B", "BOX", idAlu1, bpar, 3);
 476
 477          gMC->Gspos("C05B",1,"C05M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 478                     idrotm[1100],"ONLY");
 479          gMC->Gspos("C05B",2,"C05M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 480                     idrotm[1100],"ONLY");
 481          gMC->Gspos("C05B",3,"C05M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 482                     idrotm[1101],"ONLY");
 483          gMC->Gspos("C05B",4,"C05M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 484                     idrotm[1101],"ONLY");
 485          gMC->Gspos("C05B",5,"C05M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 486                     idrotm[1100],"ONLY");
 487          gMC->Gspos("C05B",6,"C05M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 488                     idrotm[1100],"ONLY");
 489          gMC->Gspos("C05B",7,"C05M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 490                     idrotm[1101],"ONLY");
 491          gMC->Gspos("C05B",8,"C05M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 492                     idrotm[1101],"ONLY");
 493
 494          gMC->Gspos("C06B",1,"C06M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 495                     idrotm[1100],"ONLY");
 496          gMC->Gspos("C06B",2,"C06M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 497                     idrotm[1100],"ONLY");
 498          gMC->Gspos("C06B",3,"C06M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 499                     idrotm[1101],"ONLY");
 500          gMC->Gspos("C06B",4,"C06M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 501                     idrotm[1101],"ONLY");
 502          gMC->Gspos("C06B",5,"C06M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 503                     idrotm[1100],"ONLY");
 504          gMC->Gspos("C06B",6,"C06M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 505                     idrotm[1100],"ONLY");
 506          gMC->Gspos("C06B",7,"C06M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 507                     idrotm[1101],"ONLY");
 508          gMC->Gspos("C06B",8,"C06M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 509                     idrotm[1101],"ONLY");
 510      }
 511
 512
 513 //
 514 //   Chamber Material represented by Alu sheet
 515      tpar[0]= iChamber->RInner();
 516      tpar[1]= iChamber->ROuter();
 517      tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2;
 518      gMC->Gsvolu("C05A", "TUBE", idAlu2, tpar, 3);
 519      gMC->Gsvolu("C06A", "TUBE", idAlu2, tpar, 3);
 520      gMC->Gspos("C05A", 1, "C05M", 0., 0., 0.,  0, "ONLY");
 521      gMC->Gspos("C06A", 1, "C06M", 0., 0., 0.,  0, "ONLY");
 522 //
 523 //   Sensitive volumes
 524      tpar[2] = iChamber->DGas()/2.;
 525      gMC->Gsvolu("C05G", "TUBE", idGas, tpar, 3);
 526      gMC->Gsvolu("C06G", "TUBE", idGas, tpar, 3);
 527      gMC->Gspos("C05G", 1, "C05A", 0., 0., 0.,  0, "ONLY");
 528      gMC->Gspos("C06G", 1, "C06A", 0., 0., 0.,  0, "ONLY");
 529 //
 530 // Frame Crosses to be placed inside gas
 531      if (frames) {
 532          dr = (iChamber->ROuter() - iChamber->RInner());
 533          bpar[0] = TMath::Sqrt(dr*dr-dframep*dframep/4)/2;
 534          bpar[1] = dframep/2;
 535          bpar[2] = iChamber->DGas()/2;
 536          gMC->Gsvolu("C05F", "BOX", idAlu1, bpar, 3);
 537          gMC->Gsvolu("C06F", "BOX", idAlu1, bpar, 3);
 538
 539          gMC->Gspos("C05F",1,"C05G", +iChamber->RInner()+bpar[0] , 0, 0,
 540                     idrotm[1100],"ONLY");
 541          gMC->Gspos("C05F",2,"C05G", -iChamber->RInner()-bpar[0] , 0, 0,
 542                     idrotm[1100],"ONLY");
 543          gMC->Gspos("C05F",3,"C05G", 0, +iChamber->RInner()+bpar[0] , 0,
 544                     idrotm[1101],"ONLY");
 545          gMC->Gspos("C05F",4,"C05G", 0, -iChamber->RInner()-bpar[0] , 0,
 546                     idrotm[1101],"ONLY");
 547
 548          gMC->Gspos("C06F",1,"C06G", +iChamber->RInner()+bpar[0] , 0, 0,
 549                     idrotm[1100],"ONLY");
 550          gMC->Gspos("C06F",2,"C06G", -iChamber->RInner()-bpar[0] , 0, 0,
 551                     idrotm[1100],"ONLY");
 552          gMC->Gspos("C06F",3,"C06G", 0, +iChamber->RInner()+bpar[0] , 0,
 553                     idrotm[1101],"ONLY");
 554          gMC->Gspos("C06F",4,"C06G", 0, -iChamber->RInner()-bpar[0] , 0,
 555                     idrotm[1101],"ONLY");
 556 }
 557
 558 //********************************************************************
 559 //                            Station 4                             **
 560 //********************************************************************
 561      // indices 1 and 2 for first and second chambers in the station
 562      // iChamber (first chamber) kept for other quanties than Z,
 563      // assumed to be the same in both chambers
 564      iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[6];
 565      iChamber2 =(AliMUONChamber*) (*fChambers)[7];
 566      zpos1=iChamber1->Z();
 567      zpos2=iChamber2->Z();
 568      dstation = zpos2 - zpos1;
 569      zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
 570
 571 //
 572 //   Mother volume
 573      tpar[0] = iChamber->RInner()-dframep;
 574      tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi);
 575      tpar[2] = dstation/4;
 576
 577      gMC->Gsvolu("C07M", "TUBE", idAir, tpar, 3);
 578      gMC->Gsvolu("C08M", "TUBE", idAir, tpar, 3);
 579      gMC->Gspos("C07M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
 580      gMC->Gspos("C08M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");
 581 // Aluminium frames
 582 // Outer frames
 583      pgpar[0] = 360/12/2;
 584      pgpar[1] = 360.;
 585      pgpar[2] = 12.;
 586      pgpar[3] =   2;
 587      pgpar[4] = -dframez/2;
 588      pgpar[5] = iChamber->ROuter();
 589      pgpar[6] = pgpar[5]+dframep;
 590      pgpar[7] = +dframez/2;
 591      pgpar[8] = pgpar[5];
 592      pgpar[9] = pgpar[6];
 593      gMC->Gsvolu("C07O", "PGON", idAlu1, pgpar, 10);
 594      gMC->Gsvolu("C08O", "PGON", idAlu1, pgpar, 10);
 595      gMC->Gspos("C07O",1,"C07M", 0.,0.,-zfpos,  0,"ONLY");
 596      gMC->Gspos("C07O",2,"C07M", 0.,0.,+zfpos,  0,"ONLY");
 597      gMC->Gspos("C08O",1,"C08M", 0.,0.,-zfpos,  0,"ONLY");
 598      gMC->Gspos("C08O",2,"C08M", 0.,0.,+zfpos,  0,"ONLY");
 599 //
 600 // Inner frame
 601      tpar[0]= iChamber->RInner()-dframep;
 602      tpar[1]= iChamber->RInner();
 603      tpar[2]= dframez/2;
 604      gMC->Gsvolu("C07I", "TUBE", idAlu1, tpar, 3);
 605      gMC->Gsvolu("C08I", "TUBE", idAlu1, tpar, 3);
 606
 607      gMC->Gspos("C07I",1,"C07M", 0.,0.,-zfpos,  0,"ONLY");
 608      gMC->Gspos("C07I",2,"C07M", 0.,0.,+zfpos,  0,"ONLY");
 609      gMC->Gspos("C08I",1,"C08M", 0.,0.,-zfpos,  0,"ONLY");
 610      gMC->Gspos("C08I",2,"C08M", 0.,0.,+zfpos,  0,"ONLY");
 611 //
 612 // Frame Crosses
 613      if (frames) {
 614          bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2;
 615          bpar[1] = dframep/2;
 616          bpar[2] = dframez/2;
 617          gMC->Gsvolu("C07B", "BOX", idAlu1, bpar, 3);
 618          gMC->Gsvolu("C08B", "BOX", idAlu1, bpar, 3);
 619
 620          gMC->Gspos("C07B",1,"C07M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 621                     idrotm[1100],"ONLY");
 622          gMC->Gspos("C07B",2,"C07M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 623                     idrotm[1100],"ONLY");
 624          gMC->Gspos("C07B",3,"C07M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 625                     idrotm[1101],"ONLY");
 626          gMC->Gspos("C07B",4,"C07M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 627                     idrotm[1101],"ONLY");
 628          gMC->Gspos("C07B",5,"C07M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 629                     idrotm[1100],"ONLY");
 630          gMC->Gspos("C07B",6,"C07M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 631                     idrotm[1100],"ONLY");
 632          gMC->Gspos("C07B",7,"C07M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 633                     idrotm[1101],"ONLY");
 634          gMC->Gspos("C07B",8,"C07M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 635                     idrotm[1101],"ONLY");
 636
 637          gMC->Gspos("C08B",1,"C08M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 638                     idrotm[1100],"ONLY");
 639          gMC->Gspos("C08B",2,"C08M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 640                     idrotm[1100],"ONLY");
 641          gMC->Gspos("C08B",3,"C08M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 642                     idrotm[1101],"ONLY");
 643          gMC->Gspos("C08B",4,"C08M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 644                     idrotm[1101],"ONLY");
 645          gMC->Gspos("C08B",5,"C08M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 646                     idrotm[1100],"ONLY");
 647          gMC->Gspos("C08B",6,"C08M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 648                     idrotm[1100],"ONLY");
 649          gMC->Gspos("C08B",7,"C08M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 650                     idrotm[1101],"ONLY");
 651          gMC->Gspos("C08B",8,"C08M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 652                     idrotm[1101],"ONLY");
 653      }
 654
 655
 656 //
 657 //   Chamber Material represented by Alu sheet
 658      tpar[0]= iChamber->RInner();
 659      tpar[1]= iChamber->ROuter();
 660      tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2;
 661      gMC->Gsvolu("C07A", "TUBE", idAlu2, tpar, 3);
 662      gMC->Gsvolu("C08A", "TUBE", idAlu2, tpar, 3);
 663      gMC->Gspos("C07A", 1, "C07M", 0., 0., 0.,  0, "ONLY");
 664      gMC->Gspos("C08A", 1, "C08M", 0., 0., 0.,  0, "ONLY");
 665 //
 666 //   Sensitive volumes
 667      // tpar[2] = iChamber->DGas();
 668      tpar[2] = iChamber->DGas()/2;
 669      gMC->Gsvolu("C07G", "TUBE", idGas, tpar, 3);
 670      gMC->Gsvolu("C08G", "TUBE", idGas, tpar, 3);
 671      gMC->Gspos("C07G", 1, "C07A", 0., 0., 0.,  0, "ONLY");
 672      gMC->Gspos("C08G", 1, "C08A", 0., 0., 0.,  0, "ONLY");
 673 //
 674 // Frame Crosses to be placed inside gas
 675      if (frames) {
 676          dr = (iChamber->ROuter() - iChamber->RInner());
 677          bpar[0] = TMath::Sqrt(dr*dr-dframep*dframep/4)/2;
 678          bpar[1] = dframep/2;
 679          bpar[2] = iChamber->DGas()/2;
 680          gMC->Gsvolu("C07F", "BOX", idAlu1, bpar, 3);
 681          gMC->Gsvolu("C08F", "BOX", idAlu1, bpar, 3);
 682
 683          gMC->Gspos("C07F",1,"C07G", +iChamber->RInner()+bpar[0] , 0, 0,
 684                     idrotm[1100],"ONLY");
 685          gMC->Gspos("C07F",2,"C07G", -iChamber->RInner()-bpar[0] , 0, 0,
 686                     idrotm[1100],"ONLY");
 687          gMC->Gspos("C07F",3,"C07G", 0, +iChamber->RInner()+bpar[0] , 0,
 688                     idrotm[1101],"ONLY");
 689          gMC->Gspos("C07F",4,"C07G", 0, -iChamber->RInner()-bpar[0] , 0,
 690                     idrotm[1101],"ONLY");
 691
 692          gMC->Gspos("C08F",1,"C08G", +iChamber->RInner()+bpar[0] , 0, 0,
 693                     idrotm[1100],"ONLY");
 694          gMC->Gspos("C08F",2,"C08G", -iChamber->RInner()-bpar[0] , 0, 0,
 695                     idrotm[1100],"ONLY");
 696          gMC->Gspos("C08F",3,"C08G", 0, +iChamber->RInner()+bpar[0] , 0,
 697                     idrotm[1101],"ONLY");
 698          gMC->Gspos("C08F",4,"C08G", 0, -iChamber->RInner()-bpar[0] , 0,
 699                     idrotm[1101],"ONLY");
 700      }
 701 //********************************************************************
 702 //                            Station 5                             **
 703 //********************************************************************
 704      // indices 1 and 2 for first and second chambers in the station
 705      // iChamber (first chamber) kept for other quanties than Z,
 706      // assumed to be the same in both chambers
 707      iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[8];
 708      iChamber2 =(AliMUONChamber*) (*fChambers)[9];
 709      zpos1=iChamber1->Z();
 710      zpos2=iChamber2->Z();
 711      dstation = zpos2 - zpos1;
 712      zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
 713
 714 //
 715 //   Mother volume
 716      tpar[0] = iChamber->RInner()-dframep;
 717      tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi);
 718      tpar[2] = dstation/4;
 719
 720      gMC->Gsvolu("C09M", "TUBE", idAir, tpar, 3);
 721      gMC->Gsvolu("C10M", "TUBE", idAir, tpar, 3);
 722      gMC->Gspos("C09M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
 723      gMC->Gspos("C10M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");
 724 // Aluminium frames
 725 // Outer frames
 726      pgpar[0] = 360/12/2;
 727      pgpar[1] = 360.;
 728      pgpar[2] = 12.;
 729      pgpar[3] =   2;
 730      pgpar[4] = -dframez/2;
 731      pgpar[5] = iChamber->ROuter();
 732      pgpar[6] = pgpar[5]+dframep;
 733      pgpar[7] = +dframez/2;
 734      pgpar[8] = pgpar[5];
 735      pgpar[9] = pgpar[6];
 736      gMC->Gsvolu("C09O", "PGON", idAlu1, pgpar, 10);
 737      gMC->Gsvolu("C10O", "PGON", idAlu1, pgpar, 10);
 738      gMC->Gspos("C09O",1,"C09M", 0.,0.,-zfpos,  0,"ONLY");
 739      gMC->Gspos("C09O",2,"C09M", 0.,0.,+zfpos,  0,"ONLY");
 740      gMC->Gspos("C10O",1,"C10M", 0.,0.,-zfpos,  0,"ONLY");
 741      gMC->Gspos("C10O",2,"C10M", 0.,0.,+zfpos,  0,"ONLY");
 742 //
 743 // Inner frame
 744      tpar[0]= iChamber->RInner()-dframep;
 745      tpar[1]= iChamber->RInner();
 746      tpar[2]= dframez/2;
 747      gMC->Gsvolu("C09I", "TUBE", idAlu1, tpar, 3);
 748      gMC->Gsvolu("C10I", "TUBE", idAlu1, tpar, 3);
 749
 750      gMC->Gspos("C09I",1,"C09M", 0.,0.,-zfpos,  0,"ONLY");
 751      gMC->Gspos("C09I",2,"C09M", 0.,0.,+zfpos,  0,"ONLY");
 752      gMC->Gspos("C10I",1,"C10M", 0.,0.,-zfpos,  0,"ONLY");
 753      gMC->Gspos("C10I",2,"C10M", 0.,0.,+zfpos,  0,"ONLY");
 754
 755      if (frames) {
 756 //
 757 // Frame Crosses
 758
 759          bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2;
 760          bpar[1] = dframep/2;
 761          bpar[2] = dframez/2;
 762          gMC->Gsvolu("C09B", "BOX", idAlu1, bpar, 3);
 763          gMC->Gsvolu("C10B", "BOX", idAlu1, bpar, 3);
 764
 765          gMC->Gspos("C09B",1,"C09M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 766                     idrotm[1100],"ONLY");
 767          gMC->Gspos("C09B",2,"C09M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 768                     idrotm[1100],"ONLY");
 769          gMC->Gspos("C09B",3,"C09M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 770                     idrotm[1101],"ONLY");
 771          gMC->Gspos("C09B",4,"C09M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 772                     idrotm[1101],"ONLY");
 773          gMC->Gspos("C09B",5,"C09M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 774                     idrotm[1100],"ONLY");
 775          gMC->Gspos("C09B",6,"C09M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 776                     idrotm[1100],"ONLY");
 777          gMC->Gspos("C09B",7,"C09M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 778                     idrotm[1101],"ONLY");
 779          gMC->Gspos("C09B",8,"C09M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 780                     idrotm[1101],"ONLY");
 781
 782          gMC->Gspos("C10B",1,"C10M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
 783                     idrotm[1100],"ONLY");
 784          gMC->Gspos("C10B",2,"C10M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
 785                     idrotm[1100],"ONLY");
 786          gMC->Gspos("C10B",3,"C10M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
 787                     idrotm[1101],"ONLY");
 788          gMC->Gspos("C10B",4,"C10M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
 789                     idrotm[1101],"ONLY");
 790          gMC->Gspos("C10B",5,"C10M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
 791                     idrotm[1100],"ONLY");
 792          gMC->Gspos("C10B",6,"C10M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
 793                     idrotm[1100],"ONLY");
 794          gMC->Gspos("C10B",7,"C10M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
 795                     idrotm[1101],"ONLY");
 796          gMC->Gspos("C10B",8,"C10M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
 797                     idrotm[1101],"ONLY");
 798      }
 799
 800
 801 //
 802 //   Chamber Material represented by Alu sheet
 803      tpar[0]= iChamber->RInner();
 804      tpar[1]= iChamber->ROuter();
 805      tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2;
 806      gMC->Gsvolu("C09A", "TUBE", idAlu2, tpar, 3);
 807      gMC->Gsvolu("C10A", "TUBE", idAlu2, tpar, 3);
 808      gMC->Gspos("C09A", 1, "C09M", 0., 0., 0.,  0, "ONLY");
 809      gMC->Gspos("C10A", 1, "C10M", 0., 0., 0.,  0, "ONLY");
 810 //
 811 //   Sensitive volumes
 812      // tpar[2] = iChamber->DGas();
 813      tpar[2] = iChamber->DGas()/2;
 814      gMC->Gsvolu("C09G", "TUBE", idGas, tpar, 3);
 815      gMC->Gsvolu("C10G", "TUBE", idGas, tpar, 3);
 816      gMC->Gspos("C09G", 1, "C09A", 0., 0., 0.,  0, "ONLY");
 817      gMC->Gspos("C10G", 1, "C10A", 0., 0., 0.,  0, "ONLY");
 818 //
 819 // Frame Crosses to be placed inside gas
 820      if (frames) {
 821          dr = (iChamber->ROuter() - iChamber->RInner());
 822          bpar[0] = TMath::Sqrt(dr*dr-dframep*dframep/4)/2;
 823          bpar[1] = dframep/2;
 824          bpar[2] = iChamber->DGas()/2;
 825          gMC->Gsvolu("C09F", "BOX", idAlu1, bpar, 3);
 826          gMC->Gsvolu("C10F", "BOX", idAlu1, bpar, 3);
 827
 828          gMC->Gspos("C09F",1,"C09G", +iChamber->RInner()+bpar[0] , 0, 0,
 829                     idrotm[1100],"ONLY");
 830          gMC->Gspos("C09F",2,"C09G", -iChamber->RInner()-bpar[0] , 0, 0,
 831                     idrotm[1100],"ONLY");
 832          gMC->Gspos("C09F",3,"C09G", 0, +iChamber->RInner()+bpar[0] , 0,
 833                     idrotm[1101],"ONLY");
 834          gMC->Gspos("C09F",4,"C09G", 0, -iChamber->RInner()-bpar[0] , 0,
 835                     idrotm[1101],"ONLY");
 836
 837          gMC->Gspos("C10F",1,"C10G", +iChamber->RInner()+bpar[0] , 0, 0,
 838                     idrotm[1100],"ONLY");
 839          gMC->Gspos("C10F",2,"C10G", -iChamber->RInner()-bpar[0] , 0, 0,
 840                     idrotm[1100],"ONLY");
 841          gMC->Gspos("C10F",3,"C10G", 0, +iChamber->RInner()+bpar[0] , 0,
 842                     idrotm[1101],"ONLY");
 843          gMC->Gspos("C10F",4,"C10G", 0, -iChamber->RInner()-bpar[0] , 0,
 844                     idrotm[1101],"ONLY");
 845      }
 846
 847 ///////////////////////////////////////
 848 // GEOMETRY FOR THE TRIGGER CHAMBERS //
 849 ///////////////////////////////////////
 850
 851 // 03/00 P. Dupieux : introduce a slighly more realistic
 852 //                    geom. of the trigger readout planes with
 853 //                    2 Zpos per trigger plane (alternate
 854 //                    between left and right of the trigger)
 855
 856 //  Parameters of the Trigger Chambers
 857
 858
 859      const Float_t kXMC1MIN=34.;
 860      const Float_t kXMC1MED=51.;
 861      const Float_t kXMC1MAX=272.;
 862      const Float_t kYMC1MIN=34.;
 863      const Float_t kYMC1MAX=51.;
 864      const Float_t kRMIN1=50.;
 865      const Float_t kRMAX1=62.;
 866      const Float_t kRMIN2=50.;
 867      const Float_t kRMAX2=66.;
 868
 869 //   zposition of the middle of the gas gap in mother vol
 870      const Float_t kZMCm=-3.6;
 871      const Float_t kZMCp=+3.6;
 872
 873
 874 // TRIGGER STATION 1 - TRIGGER STATION 1 - TRIGGER STATION 1
 875
 876      // iChamber 1 and 2 for first and second chambers in the station
 877      // iChamber (first chamber) kept for other quanties than Z,
 878      // assumed to be the same in both chambers
 879      iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[10];
 880      iChamber2 =(AliMUONChamber*) (*fChambers)[11];
 881
 882      // 03/00
 883      // zpos1 and zpos2 are now the middle of the first and second
 884      // plane of station 1 :
 885      // zpos1=(16075+15995)/2=16035 mm, thick/2=40 mm
 886      // zpos2=(16225+16145)/2=16185 mm, thick/2=40 mm
 887      //
 888      // zpos1m=15999 mm , zpos1p=16071 mm (middles of gas gaps)
 889      // zpos2m=16149 mm , zpos2p=16221 mm (middles of gas gaps)
 890      // rem : the total thickness accounts for 1 mm of al on both
 891      // side of the RPCs (see zpos1 and zpos2), as previously
 892
 893      zpos1=iChamber1->Z();
 894      zpos2=iChamber2->Z();
 895
 896
 897 // Mother volume definition
 898      tpar[0] = iChamber->RInner();
 899      tpar[1] = iChamber->ROuter();
 900      tpar[2] = 4.0;
 901      gMC->Gsvolu("CM11", "TUBE", idAir, tpar, 3);
 902      gMC->Gsvolu("CM12", "TUBE", idAir, tpar, 3);
 903
 904 // Definition of the flange between the beam shielding and the RPC
 905      tpar[0]= kRMIN1;
 906      tpar[1]= kRMAX1;
 907      tpar[2]= 4.0;
 908
 909      gMC->Gsvolu("CF1A", "TUBE", idAlu1, tpar, 3);     //Al
 910      gMC->Gspos("CF1A", 1, "CM11", 0., 0., 0., 0, "MANY");
 911      gMC->Gspos("CF1A", 2, "CM12", 0., 0., 0., 0, "MANY");
 912
 913
 914 // FIRST PLANE OF STATION 1
 915
 916 //   ratios of zpos1m/zpos1p and inverse for first plane
 917      Float_t zmp=(zpos1-3.6)/(zpos1+3.6);
 918      Float_t zpm=1./zmp;
 919
 920
 921 // Definition of prototype for chambers in the first plane
 922
 923      tpar[0]= 0.;
 924      tpar[1]= 0.;
 925      tpar[2]= 0.;
 926
 927      gMC->Gsvolu("CC1A", "BOX ", idAlu1, tpar, 0);           //Al
 928      gMC->Gsvolu("CB1A", "BOX ", idtmed[1107], tpar, 0);     //Bakelite
 929      gMC->Gsvolu("CG1A", "BOX ", idtmed[1106], tpar, 0);     //Gas streamer
 930
 931 // chamber type A
 932      tpar[0] = -1.;
 933      tpar[1] = -1.;
 934
 935      const Float_t kXMC1A=kXMC1MED+(kXMC1MAX-kXMC1MED)/2.;
 936      const Float_t kYMC1Am=0.;
 937      const Float_t kYMC1Ap=0.;
 938
 939      tpar[2] = 0.1;
 940      gMC->Gsposp("CG1A", 1, "CB1A", 0., 0., 0., 0, "ONLY",tpar,3);
 941      tpar[2] = 0.3;
 942      gMC->Gsposp("CB1A", 1, "CC1A", 0., 0., 0., 0, "ONLY",tpar,3);
 943
 944      tpar[2] = 0.4;
 945      tpar[0] = (kXMC1MAX-kXMC1MED)/2.;
 946      tpar[1] = kYMC1MIN;
 947
 948      gMC->Gsposp("CC1A", 1, "CM11",kXMC1A,kYMC1Am,kZMCm, 0, "ONLY", tpar, 3);
 949      gMC->Gsposp("CC1A", 2, "CM11",-kXMC1A,kYMC1Ap,kZMCp, 0, "ONLY", tpar, 3);
 950
 951 //  chamber type B
 952      Float_t tpar1save=tpar[1];
 953      Float_t y1msave=kYMC1Am;
 954      Float_t y1psave=kYMC1Ap;
 955
 956      tpar[0] = (kXMC1MAX-kXMC1MIN)/2.;
 957      tpar[1] = (kYMC1MAX-kYMC1MIN)/2.;
 958
 959      const Float_t kXMC1B=kXMC1MIN+tpar[0];
 960      const Float_t kYMC1Bp=(y1msave+tpar1save)*zpm+tpar[1];
 961      const Float_t kYMC1Bm=(y1psave+tpar1save)*zmp+tpar[1];
 962
 963      gMC->Gsposp("CC1A", 3, "CM11",kXMC1B,kYMC1Bp,kZMCp, 0, "ONLY", tpar, 3);
 964      gMC->Gsposp("CC1A", 4, "CM11",-kXMC1B,kYMC1Bm,kZMCm, 0, "ONLY", tpar, 3);
 965      gMC->Gsposp("CC1A", 5, "CM11",kXMC1B,-kYMC1Bp,kZMCp, 0, "ONLY", tpar, 3);
 966      gMC->Gsposp("CC1A", 6, "CM11",-kXMC1B,-kYMC1Bm,kZMCm, 0, "ONLY", tpar, 3);
 967
 968 //  chamber type C  (end of type B !!)
 969      tpar1save=tpar[1];
 970      y1msave=kYMC1Bm;
 971      y1psave=kYMC1Bp;
 972
 973      tpar[0] = kXMC1MAX/2;
 974      tpar[1] = kYMC1MAX/2;
 975
 976      const Float_t kXMC1C=tpar[0];
 977 // warning : same Z than type B
 978      const Float_t kYMC1Cp=(y1psave+tpar1save)*1.+tpar[1];
 979      const Float_t kYMC1Cm=(y1msave+tpar1save)*1.+tpar[1];
 980
 981      gMC->Gsposp("CC1A", 7, "CM11",kXMC1C,kYMC1Cp,kZMCp, 0, "ONLY", tpar, 3);
 982      gMC->Gsposp("CC1A", 8, "CM11",-kXMC1C,kYMC1Cm,kZMCm, 0, "ONLY", tpar, 3);
 983      gMC->Gsposp("CC1A", 9, "CM11",kXMC1C,-kYMC1Cp,kZMCp, 0, "ONLY", tpar, 3);
 984      gMC->Gsposp("CC1A", 10, "CM11",-kXMC1C,-kYMC1Cm,kZMCm, 0, "ONLY", tpar, 3);
 985
 986 //  chamber type D, E and F (same size)
 987      tpar1save=tpar[1];
 988      y1msave=kYMC1Cm;
 989      y1psave=kYMC1Cp;
 990
 991      tpar[0] = kXMC1MAX/2.;
 992      tpar[1] = kYMC1MIN;
 993
 994      const Float_t kXMC1D=tpar[0];
 995      const Float_t kYMC1Dp=(y1msave+tpar1save)*zpm+tpar[1];
 996      const Float_t kYMC1Dm=(y1psave+tpar1save)*zmp+tpar[1];
 997
 998      gMC->Gsposp("CC1A", 11, "CM11",kXMC1D,kYMC1Dm,kZMCm, 0, "ONLY", tpar, 3);
 999      gMC->Gsposp("CC1A", 12, "CM11",-kXMC1D,kYMC1Dp,kZMCp, 0, "ONLY", tpar, 3);
1000      gMC->Gsposp("CC1A", 13, "CM11",kXMC1D,-kYMC1Dm,kZMCm, 0, "ONLY", tpar, 3);
1001      gMC->Gsposp("CC1A", 14, "CM11",-kXMC1D,-kYMC1Dp,kZMCp, 0, "ONLY", tpar, 3);
1002
1003
1004      tpar1save=tpar[1];
1005      y1msave=kYMC1Dm;
1006      y1psave=kYMC1Dp;
1007      const Float_t kYMC1Ep=(y1msave+tpar1save)*zpm+tpar[1];
1008      const Float_t kYMC1Em=(y1psave+tpar1save)*zmp+tpar[1];
1009
1010      gMC->Gsposp("CC1A", 15, "CM11",kXMC1D,kYMC1Ep,kZMCp, 0, "ONLY", tpar, 3);
1011      gMC->Gsposp("CC1A", 16, "CM11",-kXMC1D,kYMC1Em,kZMCm, 0, "ONLY", tpar, 3);
1012      gMC->Gsposp("CC1A", 17, "CM11",kXMC1D,-kYMC1Ep,kZMCp, 0, "ONLY", tpar, 3);
1013      gMC->Gsposp("CC1A", 18, "CM11",-kXMC1D,-kYMC1Em,kZMCm, 0, "ONLY", tpar, 3);
1014
1015      tpar1save=tpar[1];
1016      y1msave=kYMC1Em;
1017      y1psave=kYMC1Ep;
1018      const Float_t kYMC1Fp=(y1msave+tpar1save)*zpm+tpar[1];
1019      const Float_t kYMC1Fm=(y1psave+tpar1save)*zmp+tpar[1];
1020
1021      gMC->Gsposp("CC1A", 19, "CM11",kXMC1D,kYMC1Fm,kZMCm, 0, "ONLY", tpar, 3);
1022      gMC->Gsposp("CC1A", 20, "CM11",-kXMC1D,kYMC1Fp,kZMCp, 0, "ONLY", tpar, 3);
1023      gMC->Gsposp("CC1A", 21, "CM11",kXMC1D,-kYMC1Fm,kZMCm, 0, "ONLY", tpar, 3);
1024      gMC->Gsposp("CC1A", 22, "CM11",-kXMC1D,-kYMC1Fp,kZMCp, 0, "ONLY", tpar, 3);
1025
1026 // Positioning first plane in ALICE
1027      gMC->Gspos("CM11", 1, "ALIC", 0., 0., zpos1, 0, "ONLY");
1028
1029 // End of geometry definition for the first plane of station 1
1030
1031
1032
1033 // SECOND PLANE OF STATION 1 : proj ratio = zpos2/zpos1
1034
1035      const Float_t kZ12=zpos2/zpos1;
1036
1037 // Definition of prototype for chambers in the second plane of station 1
1038
1039      tpar[0]= 0.;
1040      tpar[1]= 0.;
1041      tpar[2]= 0.;
1042
1043      gMC->Gsvolu("CC2A", "BOX ", idAlu1, tpar, 0);           //Al
1044      gMC->Gsvolu("CB2A", "BOX ", idtmed[1107], tpar, 0);     //Bakelite
1045      gMC->Gsvolu("CG2A", "BOX ", idtmed[1106], tpar, 0);     //Gas streamer
1046
1047 // chamber type A
1048      tpar[0] = -1.;
1049      tpar[1] = -1.;
1050
1051      const Float_t kXMC2A=kXMC1A*kZ12;
1052      const Float_t kYMC2Am=0.;
1053      const Float_t kYMC2Ap=0.;
1054
1055      tpar[2] = 0.1;
1056      gMC->Gsposp("CG2A", 1, "CB2A", 0., 0., 0., 0, "ONLY",tpar,3);
1057      tpar[2] = 0.3;
1058      gMC->Gsposp("CB2A", 1, "CC2A", 0., 0., 0., 0, "ONLY",tpar,3);
1059
1060      tpar[2] = 0.4;
1061      tpar[0] = ((kXMC1MAX-kXMC1MED)/2.)*kZ12;
1062      tpar[1] = kYMC1MIN*kZ12;
1063
1064      gMC->Gsposp("CC2A", 1, "CM12",kXMC2A,kYMC2Am,kZMCm, 0, "ONLY", tpar, 3);
1065      gMC->Gsposp("CC2A", 2, "CM12",-kXMC2A,kYMC2Ap,kZMCp, 0, "ONLY", tpar, 3);
1066
1067
1068 //  chamber type B
1069
1070      tpar[0] = ((kXMC1MAX-kXMC1MIN)/2.)*kZ12;
1071      tpar[1] = ((kYMC1MAX-kYMC1MIN)/2.)*kZ12;
1072
1073      const Float_t kXMC2B=kXMC1B*kZ12;
1074      const Float_t kYMC2Bp=kYMC1Bp*kZ12;
1075      const Float_t kYMC2Bm=kYMC1Bm*kZ12;
1076      gMC->Gsposp("CC2A", 3, "CM12",kXMC2B,kYMC2Bp,kZMCp, 0, "ONLY", tpar, 3);
1077      gMC->Gsposp("CC2A", 4, "CM12",-kXMC2B,kYMC2Bm,kZMCm, 0, "ONLY", tpar, 3);
1078      gMC->Gsposp("CC2A", 5, "CM12",kXMC2B,-kYMC2Bp,kZMCp, 0, "ONLY", tpar, 3);
1079      gMC->Gsposp("CC2A", 6, "CM12",-kXMC2B,-kYMC2Bm,kZMCm, 0, "ONLY", tpar, 3);
1080
1081
1082 //  chamber type C   (end of type B !!)
1083
1084      tpar[0] = (kXMC1MAX/2)*kZ12;
1085      tpar[1] = (kYMC1MAX/2)*kZ12;
1086
1087      const Float_t kXMC2C=kXMC1C*kZ12;
1088      const Float_t kYMC2Cp=kYMC1Cp*kZ12;
1089      const Float_t kYMC2Cm=kYMC1Cm*kZ12;
1090      gMC->Gsposp("CC2A", 7, "CM12",kXMC2C,kYMC2Cp,kZMCp, 0, "ONLY", tpar, 3);
1091      gMC->Gsposp("CC2A", 8, "CM12",-kXMC2C,kYMC2Cm,kZMCm, 0, "ONLY", tpar, 3);
1092      gMC->Gsposp("CC2A", 9, "CM12",kXMC2C,-kYMC2Cp,kZMCp, 0, "ONLY", tpar, 3);
1093      gMC->Gsposp("CC2A", 10, "CM12",-kXMC2C,-kYMC2Cm,kZMCm, 0, "ONLY", tpar, 3);
1094
1095 //  chamber type D, E and F (same size)
1096
1097      tpar[0] = (kXMC1MAX/2.)*kZ12;
1098      tpar[1] = kYMC1MIN*kZ12;
1099
1100      const Float_t kXMC2D=kXMC1D*kZ12;
1101      const Float_t kYMC2Dp=kYMC1Dp*kZ12;
1102      const Float_t kYMC2Dm=kYMC1Dm*kZ12;
1103      gMC->Gsposp("CC2A", 11, "CM12",kXMC2D,kYMC2Dm,kZMCm, 0, "ONLY", tpar, 3);
1104      gMC->Gsposp("CC2A", 12, "CM12",-kXMC2D,kYMC2Dp,kZMCp, 0, "ONLY", tpar, 3);
1105      gMC->Gsposp("CC2A", 13, "CM12",kXMC2D,-kYMC2Dm,kZMCm, 0, "ONLY", tpar, 3);
1106      gMC->Gsposp("CC2A", 14, "CM12",-kXMC2D,-kYMC2Dp,kZMCp, 0, "ONLY", tpar, 3);
1107
1108      const Float_t kYMC2Ep=kYMC1Ep*kZ12;
1109      const Float_t kYMC2Em=kYMC1Em*kZ12;
1110      gMC->Gsposp("CC2A", 15, "CM12",kXMC2D,kYMC2Ep,kZMCp, 0, "ONLY", tpar, 3);
1111      gMC->Gsposp("CC2A", 16, "CM12",-kXMC2D,kYMC2Em,kZMCm, 0, "ONLY", tpar, 3);
1112      gMC->Gsposp("CC2A", 17, "CM12",kXMC2D,-kYMC2Ep,kZMCp, 0, "ONLY", tpar, 3);
1113      gMC->Gsposp("CC2A", 18, "CM12",-kXMC2D,-kYMC2Em,kZMCm, 0, "ONLY", tpar, 3);
1114
1115
1116      const Float_t kYMC2Fp=kYMC1Fp*kZ12;
1117      const Float_t kYMC2Fm=kYMC1Fm*kZ12;
1118      gMC->Gsposp("CC2A", 19, "CM12",kXMC2D,kYMC2Fm,kZMCm, 0, "ONLY", tpar, 3);
1119      gMC->Gsposp("CC2A", 20, "CM12",-kXMC2D,kYMC2Fp,kZMCp, 0, "ONLY", tpar, 3);
1120      gMC->Gsposp("CC2A", 21, "CM12",kXMC2D,-kYMC2Fm,kZMCm, 0, "ONLY", tpar, 3);
1121      gMC->Gsposp("CC2A", 22, "CM12",-kXMC2D,-kYMC2Fp,kZMCp, 0, "ONLY", tpar, 3);
1122
1123 // Positioning second plane of station 1 in ALICE
1124
1125      gMC->Gspos("CM12", 1, "ALIC", 0., 0., zpos2, 0, "ONLY");
1126
1127 // End of geometry definition for the second plane of station 1
1128
1129
1130
1131 // TRIGGER STATION 2 - TRIGGER STATION 2 - TRIGGER STATION 2
1132
1133      // 03/00
1134      // zpos3 and zpos4 are now the middle of the first and second
1135      // plane of station 2 :
1136      // zpos3=(17075+16995)/2=17035 mm, thick/2=40 mm
1137      // zpos4=(17225+17145)/2=17185 mm, thick/2=40 mm
1138      //
1139      // zpos3m=16999 mm , zpos3p=17071 mm (middles of gas gaps)
1140      // zpos4m=17149 mm , zpos4p=17221 mm (middles of gas gaps)
1141      // rem : the total thickness accounts for 1 mm of al on both
1142      // side of the RPCs (see zpos3 and zpos4), as previously
1143      iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[12];
1144      iChamber2 =(AliMUONChamber*) (*fChambers)[13];
1145      Float_t zpos3=iChamber1->Z();
1146      Float_t zpos4=iChamber2->Z();
1147
1148
1149 // Mother volume definition
1150      tpar[0] = iChamber->RInner();
1151      tpar[1] = iChamber->ROuter();
1152      tpar[2] = 4.0;
1153
1154      gMC->Gsvolu("CM21", "TUBE", idAir, tpar, 3);
1155      gMC->Gsvolu("CM22", "TUBE", idAir, tpar, 3);
1156
1157 // Definition of the flange between the beam shielding and the RPC
1158 //  ???? interface shielding
1159
1160      tpar[0]= kRMIN2;
1161      tpar[1]= kRMAX2;
1162      tpar[2]= 4.0;
1163
1164      gMC->Gsvolu("CF2A", "TUBE", idAlu1, tpar, 3);            //Al
1165      gMC->Gspos("CF2A", 1, "CM21", 0., 0., 0., 0, "MANY");
1166      gMC->Gspos("CF2A", 2, "CM22", 0., 0., 0., 0, "MANY");
1167
1168
1169
1170 // FIRST PLANE OF STATION 2 : proj ratio = zpos3/zpos1
1171
1172      const Float_t kZ13=zpos3/zpos1;
1173
1174 // Definition of prototype for chambers in the first plane of station 2
1175      tpar[0]= 0.;
1176      tpar[1]= 0.;
1177      tpar[2]= 0.;
1178
1179      gMC->Gsvolu("CC3A", "BOX ", idAlu1, tpar, 0);           //Al
1180      gMC->Gsvolu("CB3A", "BOX ", idtmed[1107], tpar, 0);     //Bakelite
1181      gMC->Gsvolu("CG3A", "BOX ", idtmed[1106], tpar, 0);     //Gas streamer
1182
1183
1184 // chamber type A
1185      tpar[0] = -1.;
1186      tpar[1] = -1.;
1187
1188      const Float_t kXMC3A=kXMC1A*kZ13;
1189      const Float_t kYMC3Am=0.;
1190      const Float_t kYMC3Ap=0.;
1191
1192      tpar[2] = 0.1;
1193      gMC->Gsposp("CG3A", 1, "CB3A", 0., 0., 0., 0, "ONLY",tpar,3);
1194      tpar[2] = 0.3;
1195      gMC->Gsposp("CB3A", 1, "CC3A", 0., 0., 0., 0, "ONLY",tpar,3);
1196
1197      tpar[2] = 0.4;
1198      tpar[0] = ((kXMC1MAX-kXMC1MED)/2.)*kZ13;
1199      tpar[1] = kYMC1MIN*kZ13;
1200      gMC->Gsposp("CC3A", 1, "CM21",kXMC3A,kYMC3Am,kZMCm, 0, "ONLY", tpar, 3);
1201      gMC->Gsposp("CC3A", 2, "CM21",-kXMC3A,kYMC3Ap,kZMCp, 0, "ONLY", tpar, 3);
1202
1203
1204 //  chamber type B
1205      tpar[0] = ((kXMC1MAX-kXMC1MIN)/2.)*kZ13;
1206      tpar[1] = ((kYMC1MAX-kYMC1MIN)/2.)*kZ13;
1207
1208      const Float_t kXMC3B=kXMC1B*kZ13;
1209      const Float_t kYMC3Bp=kYMC1Bp*kZ13;
1210      const Float_t kYMC3Bm=kYMC1Bm*kZ13;
1211      gMC->Gsposp("CC3A", 3, "CM21",kXMC3B,kYMC3Bp,kZMCp, 0, "ONLY", tpar, 3);
1212      gMC->Gsposp("CC3A", 4, "CM21",-kXMC3B,kYMC3Bm,kZMCm, 0, "ONLY", tpar, 3);
1213      gMC->Gsposp("CC3A", 5, "CM21",kXMC3B,-kYMC3Bp,kZMCp, 0, "ONLY", tpar, 3);
1214      gMC->Gsposp("CC3A", 6, "CM21",-kXMC3B,-kYMC3Bm,kZMCm, 0, "ONLY", tpar, 3);
1215
1216
1217 //  chamber type C  (end of type B !!)
1218      tpar[0] = (kXMC1MAX/2)*kZ13;
1219      tpar[1] = (kYMC1MAX/2)*kZ13;
1220
1221      const Float_t kXMC3C=kXMC1C*kZ13;
1222      const Float_t kYMC3Cp=kYMC1Cp*kZ13;
1223      const Float_t kYMC3Cm=kYMC1Cm*kZ13;
1224      gMC->Gsposp("CC3A", 7, "CM21",kXMC3C,kYMC3Cp,kZMCp, 0, "ONLY", tpar, 3);
1225      gMC->Gsposp("CC3A", 8, "CM21",-kXMC3C,kYMC3Cm,kZMCm, 0, "ONLY", tpar, 3);
1226      gMC->Gsposp("CC3A", 9, "CM21",kXMC3C,-kYMC3Cp,kZMCp, 0, "ONLY", tpar, 3);
1227      gMC->Gsposp("CC3A", 10, "CM21",-kXMC3C,-kYMC3Cm,kZMCm, 0, "ONLY", tpar, 3);
1228
1229
1230 //  chamber type D, E and F (same size)
1231
1232      tpar[0] = (kXMC1MAX/2.)*kZ13;
1233      tpar[1] = kYMC1MIN*kZ13;
1234
1235      const Float_t kXMC3D=kXMC1D*kZ13;
1236      const Float_t kYMC3Dp=kYMC1Dp*kZ13;
1237      const Float_t kYMC3Dm=kYMC1Dm*kZ13;
1238      gMC->Gsposp("CC3A", 11, "CM21",kXMC3D,kYMC3Dm,kZMCm, 0, "ONLY", tpar, 3);
1239      gMC->Gsposp("CC3A", 12, "CM21",-kXMC3D,kYMC3Dp,kZMCp, 0, "ONLY", tpar, 3);
1240      gMC->Gsposp("CC3A", 13, "CM21",kXMC3D,-kYMC3Dm,kZMCm, 0, "ONLY", tpar, 3);
1241      gMC->Gsposp("CC3A", 14, "CM21",-kXMC3D,-kYMC3Dp,kZMCp, 0, "ONLY", tpar, 3);
1242
1243      const Float_t kYMC3Ep=kYMC1Ep*kZ13;
1244      const Float_t kYMC3Em=kYMC1Em*kZ13;
1245      gMC->Gsposp("CC3A", 15, "CM21",kXMC3D,kYMC3Ep,kZMCp, 0, "ONLY", tpar, 3);
1246      gMC->Gsposp("CC3A", 16, "CM21",-kXMC3D,kYMC3Em,kZMCm, 0, "ONLY", tpar, 3);
1247      gMC->Gsposp("CC3A", 17, "CM21",kXMC3D,-kYMC3Ep,kZMCp, 0, "ONLY", tpar, 3);
1248      gMC->Gsposp("CC3A", 18, "CM21",-kXMC3D,-kYMC3Em,kZMCm, 0, "ONLY", tpar, 3);
1249
1250      const Float_t kYMC3Fp=kYMC1Fp*kZ13;
1251      const Float_t kYMC3Fm=kYMC1Fm*kZ13;
1252      gMC->Gsposp("CC3A", 19, "CM21",kXMC3D,kYMC3Fm,kZMCm, 0, "ONLY", tpar, 3);
1253      gMC->Gsposp("CC3A", 20, "CM21",-kXMC3D,kYMC3Fp,kZMCp, 0, "ONLY", tpar, 3);
1254      gMC->Gsposp("CC3A", 21, "CM21",kXMC3D,-kYMC3Fm,kZMCm, 0, "ONLY", tpar, 3);
1255      gMC->Gsposp("CC3A", 22, "CM21",-kXMC3D,-kYMC3Fp,kZMCp, 0, "ONLY", tpar, 3);
1256
1257
1258 // Positioning first plane of station 2 in ALICE
1259
1260      gMC->Gspos("CM21", 1, "ALIC", 0., 0., zpos3, 0, "ONLY");
1261
1262 // End of geometry definition for the first plane of station 2
1263
1264
1265
1266
1267 // SECOND PLANE OF STATION 2 : proj ratio = zpos4/zpos1
1268
1269      const Float_t kZ14=zpos4/zpos1;
1270
1271 // Definition of prototype for chambers in the second plane of station 2
1272
1273      tpar[0]= 0.;
1274      tpar[1]= 0.;
1275      tpar[2]= 0.;
1276
1277      gMC->Gsvolu("CC4A", "BOX ", idAlu1, tpar, 0);           //Al
1278      gMC->Gsvolu("CB4A", "BOX ", idtmed[1107], tpar, 0);     //Bakelite
1279      gMC->Gsvolu("CG4A", "BOX ", idtmed[1106], tpar, 0);     //Gas streamer
1280
1281 // chamber type A
1282      tpar[0] = -1.;
1283      tpar[1] = -1.;
1284
1285      const Float_t kXMC4A=kXMC1A*kZ14;
1286      const Float_t kYMC4Am=0.;
1287      const Float_t kYMC4Ap=0.;
1288
1289      tpar[2] = 0.1;
1290      gMC->Gsposp("CG4A", 1, "CB4A", 0., 0., 0., 0, "ONLY",tpar,3);
1291      tpar[2] = 0.3;
1292      gMC->Gsposp("CB4A", 1, "CC4A", 0., 0., 0., 0, "ONLY",tpar,3);
1293
1294      tpar[2] = 0.4;
1295      tpar[0] = ((kXMC1MAX-kXMC1MED)/2.)*kZ14;
1296      tpar[1] = kYMC1MIN*kZ14;
1297      gMC->Gsposp("CC4A", 1, "CM22",kXMC4A,kYMC4Am,kZMCm, 0, "ONLY", tpar, 3);
1298      gMC->Gsposp("CC4A", 2, "CM22",-kXMC4A,kYMC4Ap,kZMCp, 0, "ONLY", tpar, 3);
1299
1300
1301 //  chamber type B
1302      tpar[0] = ((kXMC1MAX-kXMC1MIN)/2.)*kZ14;
1303      tpar[1] = ((kYMC1MAX-kYMC1MIN)/2.)*kZ14;
1304
1305      const Float_t kXMC4B=kXMC1B*kZ14;
1306      const Float_t kYMC4Bp=kYMC1Bp*kZ14;
1307      const Float_t kYMC4Bm=kYMC1Bm*kZ14;
1308      gMC->Gsposp("CC4A", 3, "CM22",kXMC4B,kYMC4Bp,kZMCp, 0, "ONLY", tpar, 3);
1309      gMC->Gsposp("CC4A", 4, "CM22",-kXMC4B,kYMC4Bm,kZMCm, 0, "ONLY", tpar, 3);
1310      gMC->Gsposp("CC4A", 5, "CM22",kXMC4B,-kYMC4Bp,kZMCp, 0, "ONLY", tpar, 3);
1311      gMC->Gsposp("CC4A", 6, "CM22",-kXMC4B,-kYMC4Bm,kZMCm, 0, "ONLY", tpar, 3);
1312
1313
1314 //  chamber type C   (end of type B !!)
1315      tpar[0] =(kXMC1MAX/2)*kZ14;
1316      tpar[1] =  (kYMC1MAX/2)*kZ14;
1317
1318      const Float_t kXMC4C=kXMC1C*kZ14;
1319      const Float_t kYMC4Cp=kYMC1Cp*kZ14;
1320      const Float_t kYMC4Cm=kYMC1Cm*kZ14;
1321      gMC->Gsposp("CC4A", 7, "CM22",kXMC4C,kYMC4Cp,kZMCp, 0, "ONLY", tpar, 3);
1322      gMC->Gsposp("CC4A", 8, "CM22",-kXMC4C,kYMC4Cm,kZMCm, 0, "ONLY", tpar, 3);
1323      gMC->Gsposp("CC4A", 9, "CM22",kXMC4C,-kYMC4Cp,kZMCp, 0, "ONLY", tpar, 3);
1324      gMC->Gsposp("CC4A", 10, "CM22",-kXMC4C,-kYMC4Cm,kZMCm, 0, "ONLY", tpar, 3);
1325
1326
1327 //  chamber type D, E and F (same size)
1328      tpar[0] = (kXMC1MAX/2.)*kZ14;
1329      tpar[1] =  kYMC1MIN*kZ14;
1330
1331      const Float_t kXMC4D=kXMC1D*kZ14;
1332      const Float_t kYMC4Dp=kYMC1Dp*kZ14;
1333      const Float_t kYMC4Dm=kYMC1Dm*kZ14;
1334      gMC->Gsposp("CC4A", 11, "CM22",kXMC4D,kYMC4Dm,kZMCm, 0, "ONLY", tpar, 3);
1335      gMC->Gsposp("CC4A", 12, "CM22",-kXMC4D,kYMC4Dp,kZMCp, 0, "ONLY", tpar, 3);
1336      gMC->Gsposp("CC4A", 13, "CM22",kXMC4D,-kYMC4Dm,kZMCm, 0, "ONLY", tpar, 3);
1337      gMC->Gsposp("CC4A", 14, "CM22",-kXMC4D,-kYMC4Dp,kZMCp, 0, "ONLY", tpar, 3);
1338
1339      const Float_t kYMC4Ep=kYMC1Ep*kZ14;
1340      const Float_t kYMC4Em=kYMC1Em*kZ14;
1341      gMC->Gsposp("CC4A", 15, "CM22",kXMC4D,kYMC4Ep,kZMCp, 0, "ONLY", tpar, 3);
1342      gMC->Gsposp("CC4A", 16, "CM22",-kXMC4D,kYMC4Em,kZMCm, 0, "ONLY", tpar, 3);
1343      gMC->Gsposp("CC4A", 17, "CM22",kXMC4D,-kYMC4Ep,kZMCp, 0, "ONLY", tpar, 3);
1344      gMC->Gsposp("CC4A", 18, "CM22",-kXMC4D,-kYMC4Em,kZMCm, 0, "ONLY", tpar, 3);
1345
1346      const Float_t kYMC4Fp=kYMC1Fp*kZ14;
1347      const Float_t kYMC4Fm=kYMC1Fm*kZ14;
1348      gMC->Gsposp("CC4A", 19, "CM22",kXMC4D,kYMC4Fm,kZMCm, 0, "ONLY", tpar, 3);
1349      gMC->Gsposp("CC4A", 20, "CM22",-kXMC4D,kYMC4Fp,kZMCp, 0, "ONLY", tpar, 3);
1350      gMC->Gsposp("CC4A", 21, "CM22",kXMC4D,-kYMC4Fm,kZMCm, 0, "ONLY", tpar, 3);
1351      gMC->Gsposp("CC4A", 22, "CM22",-kXMC4D,-kYMC4Fp,kZMCp, 0, "ONLY", tpar, 3);
1352
1353
1354 // Positioning second plane of station 2 in ALICE
1355
1356      gMC->Gspos("CM22", 1, "ALIC", 0., 0., zpos4, 0, "ONLY");
1357
1358 // End of geometry definition for the second plane of station 2
1359
1360 // End of trigger geometry definition
1361
1362 }
1363
1364
1365
1366 //___________________________________________
1367 void AliMUONv1::CreateMaterials()
1368 {
1369   // *** DEFINITION OF AVAILABLE MUON MATERIALS ***
1370   //
1371   //     Ar-CO2 gas
1372     Float_t ag1[3]   = { 39.95,12.01,16. };
1373     Float_t zg1[3]   = { 18.,6.,8. };
1374     Float_t wg1[3]   = { .8,.0667,.13333 };
1375     Float_t dg1      = .001821;
1376     //
1377     //     Ar-buthane-freon gas -- trigger chambers
1378     Float_t atr1[4]  = { 39.95,12.01,1.01,19. };
1379     Float_t ztr1[4]  = { 18.,6.,1.,9. };
1380     Float_t wtr1[4]  = { .56,.1262857,.2857143,.028 };
1381     Float_t dtr1     = .002599;
1382     //
1383     //     Ar-CO2 gas
1384     Float_t agas[3]  = { 39.95,12.01,16. };
1385     Float_t zgas[3]  = { 18.,6.,8. };
1386     Float_t wgas[3]  = { .74,.086684,.173316 };
1387     Float_t dgas     = .0018327;
1388     //
1389     //     Ar-Isobutane gas (80%+20%) -- tracking
1390     Float_t ag[3]    = { 39.95,12.01,1.01 };
1391     Float_t zg[3]    = { 18.,6.,1. };
1392     Float_t wg[3]    = { .8,.057,.143 };
1393     Float_t dg       = .0019596;
1394     //
1395     //     Ar-Isobutane-Forane-SF6 gas (49%+7%+40%+4%) -- trigger
1396     Float_t atrig[5] = { 39.95,12.01,1.01,19.,32.066 };
1397     Float_t ztrig[5] = { 18.,6.,1.,9.,16. };
1398     Float_t wtrig[5] = { .49,1.08,1.5,1.84,0.04 };
1399     Float_t dtrig    = .0031463;
1400     //
1401     //     bakelite
1402
1403     Float_t abak[3] = {12.01 , 1.01 , 16.};
1404     Float_t zbak[3] = {6.     , 1.   , 8.};
1405     Float_t wbak[3] = {6.     , 6.   , 1.};
1406     Float_t dbak = 1.4;
1407
1408     Float_t epsil, stmin, deemax, tmaxfd, stemax;
1409
1410     Int_t iSXFLD   = gAlice->Field()->Integ();
1411     Float_t sXMGMX = gAlice->Field()->Max();
1412     //
1413     // --- Define the various materials for GEANT ---
1414     AliMaterial(9, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
1415     AliMaterial(10, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
1416     AliMaterial(15, "AIR$      ", 14.61, 7.3, .001205, 30423.24, 67500);
1417     AliMixture(19, "Bakelite$", abak, zbak, dbak, -3, wbak);
1418     AliMixture(20, "ArC4H10 GAS$", ag, zg, dg, 3, wg);
1419     AliMixture(21, "TRIG GAS$", atrig, ztrig, dtrig, -5, wtrig);
1420     AliMixture(22, "ArCO2 80%$", ag1, zg1, dg1, 3, wg1);
1421     AliMixture(23, "Ar-freon $", atr1, ztr1, dtr1, 4, wtr1);
1422     AliMixture(24, "ArCO2 GAS$", agas, zgas, dgas, 3, wgas);
1423
1424     epsil  = .001; // Tracking precision,
1425     stemax = -1.;  // Maximum displacement for multiple scat
1426     tmaxfd = -20.; // Maximum angle due to field deflection
1427     deemax = -.3;  // Maximum fractional energy loss, DLS
1428     stmin  = -.8;
1429     //
1430     //    Air
1431     AliMedium(1, "AIR_CH_US         ", 15, 1, iSXFLD, sXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
1432     //
1433     //    Aluminum
1434
1435     AliMedium(4, "ALU_CH_US          ", 9, 0, iSXFLD, sXMGMX, tmaxfd, fMaxStepAlu,
1436             fMaxDestepAlu, epsil, stmin);
1437     AliMedium(5, "ALU_CH_US          ", 10, 0, iSXFLD, sXMGMX, tmaxfd, fMaxStepAlu,
1438             fMaxDestepAlu, epsil, stmin);
1439     //
1440     //    Ar-isoC4H10 gas
1441
1442     AliMedium(6, "AR_CH_US          ", 20, 1, iSXFLD, sXMGMX, tmaxfd, fMaxStepGas,
1443             fMaxDestepGas, epsil, stmin);
1444 //
1445     //    Ar-Isobuthane-Forane-SF6 gas
1446
1447     AliMedium(7, "GAS_CH_TRIGGER    ", 21, 1, iSXFLD, sXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
1448
1449     AliMedium(8, "BAKE_CH_TRIGGER   ", 19, 0, iSXFLD, sXMGMX, tmaxfd, fMaxStepAlu,
1450             fMaxDestepAlu, epsil, stmin);
1451
1452     AliMedium(9, "ARG_CO2   ", 22, 1, iSXFLD, sXMGMX, tmaxfd, fMaxStepGas,
1453             fMaxDestepAlu, epsil, stmin);
1454 }
1455
1456 //___________________________________________
1457
1458 void AliMUONv1::Init()
1459 {
1460    //
1461    // Initialize Tracking Chambers
1462    //
1463
1464    printf("\n\n\n Start Init for version 1 - CPC chamber type\n\n\n");
1465    Int_t i;
1466    for (i=0; i<AliMUONConstants::NCh(); i++) {
1467        ( (AliMUONChamber*) (*fChambers)[i])->Init();
1468    }
1469
1470    //
1471    // Set the chamber (sensitive region) GEANT identifier
1472    AliMC* gMC = AliMC::GetMC();
1473    ((AliMUONChamber*)(*fChambers)[0])->SetGid(gMC->VolId("C01G"));
1474    ((AliMUONChamber*)(*fChambers)[1])->SetGid(gMC->VolId("C02G"));
1475    ((AliMUONChamber*)(*fChambers)[2])->SetGid(gMC->VolId("C03G"));
1476    ((AliMUONChamber*)(*fChambers)[3])->SetGid(gMC->VolId("C04G"));
1477    ((AliMUONChamber*)(*fChambers)[4])->SetGid(gMC->VolId("C05G"));
1478    ((AliMUONChamber*)(*fChambers)[5])->SetGid(gMC->VolId("C06G"));
1479    ((AliMUONChamber*)(*fChambers)[6])->SetGid(gMC->VolId("C07G"));
1480    ((AliMUONChamber*)(*fChambers)[7])->SetGid(gMC->VolId("C08G"));
1481    ((AliMUONChamber*)(*fChambers)[8])->SetGid(gMC->VolId("C09G"));
1482    ((AliMUONChamber*)(*fChambers)[9])->SetGid(gMC->VolId("C10G"));
1483    ((AliMUONChamber*)(*fChambers)[10])->SetGid(gMC->VolId("CG1A"));
1484    ((AliMUONChamber*)(*fChambers)[11])->SetGid(gMC->VolId("CG2A"));
1485    ((AliMUONChamber*)(*fChambers)[12])->SetGid(gMC->VolId("CG3A"));
1486    ((AliMUONChamber*)(*fChambers)[13])->SetGid(gMC->VolId("CG4A"));
1487
1488    printf("\n\n\n Finished Init for version 0 - CPC chamber type\n\n\n");
1489
1490    //cp
1491    printf("\n\n\n Start Init for Trigger Circuits\n\n\n");
1492    for (i=0; i<AliMUONConstants::NTriggerCircuit(); i++) {
1493      ( (AliMUONTriggerCircuit*) (*fTriggerCircuits)[i])->Init(i);
1494    }
1495    printf(" Finished Init for Trigger Circuits\n\n\n");
1496    //cp
1497
1498 }
1499
1500 //___________________________________________
1501 void AliMUONv1::StepManager()
1502 {
1503   Int_t          copy, id;
1504   static Int_t   idvol;
1505   static Int_t   vol[2];
1506   Int_t          ipart;
1507   TLorentzVector pos;
1508   TLorentzVector mom;
1509   Float_t        theta,phi;
1510   Float_t        destep, step;
1511
1512   static Float_t eloss, eloss2, xhit, yhit, tof, tlength;
1513   const  Float_t kBig=1.e10;
1514   //  modifs perso
1515   static Float_t hits[15];
1516
1517   TClonesArray &lhits = *fHits;
1518
1519   //
1520   // Set maximum step size for gas
1521   // numed=gMC->GetMedium();
1522   //
1523   // Only charged tracks
1524   if( !(gMC->TrackCharge()) ) return;
1525   //
1526   // Only gas gap inside chamber
1527   // Tag chambers and record hits when track enters
1528   idvol=-1;
1529   id=gMC->CurrentVolID(copy);
1530
1531     for (Int_t i=1; i<=AliMUONConstants::NCh(); i++) {
1532       if(id==((AliMUONChamber*)(*fChambers)[i-1])->GetGid()){
1533           vol[0]=i;
1534           idvol=i-1;
1535       }
1536     }
1537     if (idvol == -1) return;
1538   //
1539   // Get current particle id (ipart), track position (pos)  and momentum (mom)
1540   gMC->TrackPosition(pos);
1541   gMC->TrackMomentum(mom);
1542
1543   ipart  = gMC->TrackPid();
1544   //Int_t ipart1 = gMC->IdFromPDG(ipart);
1545   //printf("ich, ipart %d %d \n",vol[0],ipart1);
1546
1547   //
1548   // momentum loss and steplength in last step
1549   destep = gMC->Edep();
1550   step   = gMC->TrackStep();
1551
1552   //
1553   // record hits when track enters ...
1554   if( gMC->IsTrackEntering()) {
1555       gMC->SetMaxStep(fMaxStepGas);
1556       Double_t tc = mom[0]*mom[0]+mom[1]*mom[1];
1557       Double_t rt = TMath::Sqrt(tc);
1558       Double_t pmom = TMath::Sqrt(tc+mom[2]*mom[2]);
1559       Double_t tx=mom[0]/pmom;
1560       Double_t ty=mom[1]/pmom;
1561       Double_t tz=mom[2]/pmom;
1562       Double_t s=((AliMUONChamber*)(*fChambers)[idvol])
1563           ->ResponseModel()
1564           ->Pitch()/tz;
1565       theta   = Float_t(TMath::ATan2(rt,Double_t(mom[2])))*kRaddeg;
1566       phi     = Float_t(TMath::ATan2(Double_t(mom[1]),Double_t(mom[0])))*kRaddeg;
1567       hits[0] = Float_t(ipart);         // Geant3 particle type
1568       hits[1] = pos[0]+s*tx;                 // X-position for hit
1569       hits[2] = pos[1]+s*ty;                 // Y-position for hit
1570       hits[3] = pos[2]+s*tz;                 // Z-position for hit
1571       hits[4] = theta;                  // theta angle of incidence
1572       hits[5] = phi;                    // phi angle of incidence
1573       hits[8] = (Float_t) fNPadHits;   // first padhit
1574       hits[9] = -1;                     // last pad hit
1575
1576       // modifs perso
1577       hits[10] = mom[3]; // hit momentum P
1578       hits[11] = mom[0]; // Px/P
1579       hits[12] = mom[1]; // Py/P
1580       hits[13] = mom[2]; // Pz/P
1581       // fin modifs perso
1582       tof=gMC->TrackTime();
1583       hits[14] = tof;    // Time of flight
1584       // phi angle of incidence
1585       tlength = 0;
1586       eloss   = 0;
1587       eloss2  = 0;
1588       xhit    = pos[0];
1589       yhit    = pos[1];
1590       // Only if not trigger chamber
1591       if(idvol<10) {
1592           //
1593           //  Initialize hit position (cursor) in the segmentation model
1594           ((AliMUONChamber*) (*fChambers)[idvol])
1595               ->SigGenInit(pos[0], pos[1], pos[2]);
1596       } else {
1597           //geant3->Gpcxyz();
1598           //printf("In the Trigger Chamber #%d\n",idvol-9);
1599       }
1600   }
1601   eloss2+=destep;
1602
1603   //
1604   // Calculate the charge induced on a pad (disintegration) in case
1605   //
1606   // Mip left chamber ...
1607   if( gMC->IsTrackExiting() || gMC->IsTrackStop() || gMC->IsTrackDisappeared()){
1608       gMC->SetMaxStep(kBig);
1609       eloss   += destep;
1610       tlength += step;
1611
1612       Float_t x0,y0;
1613
1614       if(idvol<10) {
1615 // tracking chambers
1616           x0 = 0.5*(xhit+pos[0]);
1617           y0 = 0.5*(yhit+pos[1]);
1618       } else {
1619 // trigger chambers
1620           x0=xhit;
1621           y0=yhit;
1622       }
1623
1624
1625       if (eloss >0)  MakePadHits(x0,y0,eloss,tof,idvol);
1626
1627
1628       hits[6]=tlength;
1629       hits[7]=eloss2;
1630       if (fNPadHits > (Int_t)hits[8]) {
1631           hits[8]= hits[8]+1;
1632           hits[9]= (Float_t) fNPadHits;
1633       }
1634
1635       new(lhits[fNhits++])
1636           AliMUONHit(fIshunt,gAlice->CurrentTrack(),vol,hits);
1637       eloss = 0;
1638       //
1639       // Check additional signal generation conditions
1640       // defined by the segmentation
1641       // model (boundary crossing conditions)
1642   } else if
1643       (((AliMUONChamber*) (*fChambers)[idvol])
1644        ->SigGenCond(pos[0], pos[1], pos[2]))
1645   {
1646       ((AliMUONChamber*) (*fChambers)[idvol])
1647           ->SigGenInit(pos[0], pos[1], pos[2]);
1648 //      printf("\n-> MakePadHits, reason special %d",ipart);
1649       if (eloss > 0 && idvol < 10)
1650           MakePadHits(0.5*(xhit+pos[0]),0.5*(yhit+pos[1]),eloss,tof,idvol);
1651       xhit     = pos[0];
1652       yhit     = pos[1];
1653       eloss    = destep;
1654       tlength += step ;
1655       //
1656       // nothing special  happened, add up energy loss
1657   } else {
1658       eloss   += destep;
1659       tlength += step ;
1660   }
1661 }
1662
1663