FMD/AliFMDv1.cxx

   1 ///////////////////////////////////////////////////////////////////////////////
   2 //                                                                           //
   3 //  Forward Multiplicity Detector Version 1                                  //
   4 //                                                                           //
   5 //Begin_Html
   6 /*
   7 <img src="picts/AliFMDv1Class.gif">
   8 </pre>
   9 <br clear=left>
  10 <font size=+2 color=red>
  11 <p>The responsible person for this module is
  12 <a href="mailto:Valeri.Kondratiev@cern.ch">Valeri Kondratiev</a>.
  13 </font>
  14 <pre>
  15 */
  16 //End_Html
  17 //                                                                           //
  18 ///////////////////////////////////////////////////////////////////////////////
  19
  20 #include "AliRun.h"
  21 #include "AliFMDv1.h"
  22 #include "AliMC.h"
  23 #include "AliConst.h"
  24
  25
  26 ClassImp(AliFMDv1)
  27
  28 //_____________________________________________________________________________
  29 AliFMDv1::AliFMDv1()
  30 {
  31   //
  32   // Defautl constructor for FMD version 1
  33   //
  34 }
  35
  36 //_____________________________________________________________________________
  37 AliFMDv1::AliFMDv1(const char *name, const char *title)
  38   : AliFMD(name,title)
  39 {
  40   //
  41   // Standard constructor for FMD version 1
  42   //
  43   AliModule *start = gAlice->GetModule("START");
  44   if(start) {
  45     Error("ctor","This version of FMD is incompatible with START\n");
  46     exit(1);
  47   }
  48 }
  49
  50 //_____________________________________________________________________________
  51 void AliFMDv1::CreateGeometry()
  52 {
  53   //
  54   // Creation of the geometry of the FMD version 1
  55   //
  56   //Begin_Html
  57   /*
  58     <img src="picts/AliFMDv1Tree.gif">
  59   */
  60   //End_Html
  61   //Begin_Html
  62   /*
  63     <img src="picts/AliFMDv1.gif">
  64   */
  65   //End_Html
  66
  67
  68   Float_t rout;
  69   Float_t z;
  70   Float_t h1, h2, t0, t1, t2;
  71   Float_t tt;
  72   Float_t par[3], rin;
  73   const Float_t v1 = 42.0;
  74   const Float_t v2 = 4.5;
  75   const Float_t v3 = 5.62;
  76   const Float_t v4 = 16.0;
  77
  78   Int_t *idtmed = fIdtmed->GetArray()-899;
  79
  80   // ********************************************************
  81   //     DEFINE RIGHT DISK#3  OF FMD
  82   // ********************************************************
  83
  84   //     Define parameters
  85
  86   rin  = 4.5;
  87   rout = 10.5;
  88   z    = 85.;
  89   tt   = 2.5;
  90   h1   = TMath::Sqrt(rout * rout + v1/ (tt * kPI)) - rout;
  91   h2   = .4;
  92   t0   = v2/ (h1 * kPI * (h1 + rout * 2.));
  93   t1   = v3/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
  94   t2   = v4/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
  95
  96   // *******Inner slice***********
  97   //     Inner steel wall
  98   par[0] = rin - .02;
  99   par[1] = rin;
 100   par[2] = 1.5;
 101   gMC->Gsvolu("IWR3", "TUBE", idtmed[899], par, 3);
 102   gMC->Gspos("IWR3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
 103   //     Front steel wall
 104   par[0] = rin;
 105   par[1] = rout;
 106   par[2] = .01;
 107   gMC->Gsvolu("FWR3", "TUBE", idtmed[899], par, 3);
 108   gMC->Gspos("FWR3", 1, "ALIC", 0., 0., z + .01, 0, "ONLY");
 109   //     Rear steel wall
 110   gMC->Gsvolu("RWR3", "TUBE", idtmed[899], par, 3);
 111   gMC->Gspos("RWR3", 1, "ALIC", 0., 0., z + 2.99, 0, "ONLY");
 112   //     MCP
 113   par[0] = rin;
 114   par[1] = rout;
 115   par[2] = .07;
 116   gMC->Gsvolu("MPR3", "TUBE", idtmed[900], par, 3);
 117   gMC->Gspos("MPR3", 1, "ALIC", 0., 0., z + 1.57, 0, "ONLY");
 118   //     Silicon plate
 119   par[0] = rin;
 120   par[1] = rout;
 121   par[2] = .019;
 122   gMC->Gsvolu("SPR3", "TUBE", idtmed[901], par, 3);
 123   gMC->Gspos("SPR3", 1, "ALIC", 0., 0., z + 1.719, 0, "ONLY");
 124   //     Summator plate
 125   par[0] = rin;
 126   par[1] = rout;
 127   par[2] = .01;
 128   gMC->Gsvolu("SMR3", "TUBE", idtmed[902], par, 3);
 129   gMC->Gspos("SMR3", 1, "ALIC", 0., 0., z + 2.01, 0, "ONLY");
 130   // *******Outer slice *******
 131   //     Ceramic plate
 132   par[0] = rout;
 133   par[1] = rout + h1;
 134   par[2] = 1.25;
 135   gMC->Gsvolu("CPR3", "TUBE", idtmed[902], par, 3);
 136   gMC->Gspos("CPR3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
 137   //     Covar spring
 138   par[0] = rout;
 139   par[1] = rout + h1;
 140   par[2] = t0 / 2.;
 141   gMC->Gsvolu("C1R3", "TUBE", idtmed[903], par, 3);
 142   gMC->Gspos("C1R3", 1, "ALIC", 0., 0., z + .25 - t0 / 2., 0, "ONLY");
 143   gMC->Gsvolu("C2R3", "TUBE", idtmed[903], par, 3);
 144   gMC->Gspos("C2R3", 1, "ALIC", 0., 0., z + 2.75 + t0 / 2., 0, "ONLY");
 145   //     Getter camera
 146   par[0] = rout + h1;
 147   par[1] = rout + h1 + h2;
 148   par[2] = t1 / 2.;
 149   gMC->Gsvolu("GKR3", "TUBE", idtmed[903], par, 3);
 150   gMC->Gspos("GKR3", 1, "ALIC", 0., 0., z + .25 + t1 / 2., 0, "ONLY");
 151   //     Ceramic slice
 152   par[0] = rout + h1;
 153   par[1] = rout + h1 + h2;
 154   par[2] = t2 / 2.;
 155   gMC->Gsvolu("SCR3", "TUBE", idtmed[902], par, 3);
 156   gMC->Gspos("SCR3", 1, "ALIC", 0., 0., z + .25 + t1 + t2 / 2., 0, "ONLY");
 157   // ******Electronic slice *******
 158   //     Silicon ring
 159   par[0] = rout + h1 + h2;
 160   par[1] = rout + h1 + h2 + 5.;
 161   par[2] = .025;
 162   gMC->Gsvolu("SER3", "TUBE", idtmed[901], par, 3);
 163   gMC->Gspos("SER3", 1, "ALIC", 0., 0., z + 1.57 - .025, 0, "ONLY");
 164   //     Ceramic ring
 165   par[0] = rout + h1 + h2;
 166   par[1] = rout + h1 + h2 + 5.;
 167   par[2] = .025;
 168   gMC->Gsvolu("CER3", "TUBE", idtmed[902], par, 3);
 169   gMC->Gspos("CER3", 1, "ALIC", 0., 0., z + 1.58 + .025, 0, "ONLY");
 170   // ***********************************************************
 171   //     DEFINE LEFT DISK#3 OF FMD
 172   // ***********************************************************
 173
 174   //     Define parameters
 175
 176   rin  = 4.5;
 177   rout = 10.5;
 178   z    = -85.;
 179   tt   = 2.5;
 180   h1   = TMath::Sqrt(rout * rout + v1/ (tt * kPI)) - rout;
 181   h2   = .4;
 182   t0   = v2/ (h1 * kPI * (h1 + rout * 2.));
 183   t1   = v3/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
 184   t2   = v4/  (h2 * kPI * (h2 + (h1 + rout) * 2.));
 185
 186   // *******Inner slice***********
 187   //     Inner steel wall
 188   par[0] = rin - .02;
 189   par[1] = rin;
 190   par[2] = 1.5;
 191   gMC->Gsvolu("IWL3", "TUBE", idtmed[899], par, 3);
 192   gMC->Gspos("IWL3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
 193   //     Front steel wall
 194   par[0] = rin;
 195   par[1] = rout;
 196   par[2] = .01;
 197   gMC->Gsvolu("FWL3", "TUBE", idtmed[899], par, 3);
 198   gMC->Gspos("FWL3", 1, "ALIC", 0., 0., z - .01, 0, "ONLY");
 199   //     Rear steel wall
 200   gMC->Gsvolu("RWL3", "TUBE", idtmed[899], par, 3);
 201   gMC->Gspos("RWL3", 1, "ALIC", 0., 0., z - 2.99, 0, "ONLY");
 202   //     MCP
 203   par[0] = rin;
 204   par[1] = rout;
 205   par[2] = .07;
 206   gMC->Gsvolu("MPL3", "TUBE", idtmed[900], par, 3);
 207   gMC->Gspos("MPL3", 1, "ALIC", 0., 0., z - 1.57, 0, "ONLY");
 208   //     Silicon plate
 209   par[0] = rin;
 210   par[1] = rout;
 211   par[2] = .019;
 212   gMC->Gsvolu("SPL3", "TUBE", idtmed[901], par, 3);
 213   gMC->Gspos("SPL3", 1, "ALIC", 0., 0., z - 1.719, 0, "ONLY");
 214   //     Summator plate
 215   par[0] = rin;
 216   par[1] = rout;
 217   par[2] = .01;
 218   gMC->Gsvolu("SML3", "TUBE", idtmed[902], par, 3);
 219   gMC->Gspos("SML3", 1, "ALIC", 0., 0., z - 2.01, 0, "ONLY");
 220   // *******Outer slice *******
 221   //     Ceramic plate
 222   par[0] = rout;
 223   par[1] = rout + h1;
 224   par[2] = 1.25;
 225   gMC->Gsvolu("CPL3", "TUBE", idtmed[902], par, 3);
 226   gMC->Gspos("CPL3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
 227   //     Covar spring
 228   par[0] = rout;
 229   par[1] = rout + h1;
 230   par[2] = t0 / 2.;
 231   gMC->Gsvolu("C1L3", "TUBE", idtmed[903], par, 3);
 232   gMC->Gspos("C1L3", 1, "ALIC", 0., 0., z - .25 + t0 / 2., 0, "ONLY");
 233   gMC->Gsvolu("C2L3", "TUBE", idtmed[903], par, 3);
 234   gMC->Gspos("C2L3", 1, "ALIC", 0., 0., z - 2.75 - t0 / 2., 0, "ONLY");
 235   //     Getter camera
 236   par[0] = rout + h1;
 237   par[1] = rout + h1 + h2;
 238   par[2] = t1 / 2.;
 239   gMC->Gsvolu("GKL3", "TUBE", idtmed[903], par, 3);
 240   gMC->Gspos("GKL3", 1, "ALIC", 0., 0., z - .25 - t1 / 2., 0, "ONLY");
 241   //     Ceramic slice
 242   par[0] = rout + h1;
 243   par[1] = rout + h1 + h2;
 244   par[2] = t2 / 2.;
 245   gMC->Gsvolu("SCL3", "TUBE", idtmed[902], par, 3);
 246   gMC->Gspos("SCL3", 1, "ALIC", 0., 0., z - .25 - t1 - t2 / 2., 0, "ONLY");
 247   // ******Electronic slice *******
 248   //     Silicon ring
 249   par[0] = rout + h1 + h2;
 250   par[1] = rout + h1 + h2 + 5.;
 251   par[2] = .025;
 252   gMC->Gsvolu("SEL3", "TUBE", idtmed[901], par, 3);
 253   gMC->Gspos("SEL3", 1, "ALIC", 0., 0., z - 1.57 + .025, 0, "ONLY");
 254   //     Ceramic ring
 255   par[0] = rout + h1 + h2;
 256   par[1] = rout + h1 + h2 + 5.;
 257   par[2] = .025;
 258   gMC->Gsvolu("CEL3", "TUBE", idtmed[902], par, 3);
 259   gMC->Gspos("CEL3", 1, "ALIC", 0., 0., z - 1.58 - .025, 0, "ONLY");
 260   // ********************************************************
 261   //     DEFINE RIGHT DISK#2  OF FMD
 262   // ********************************************************
 263
 264   //     Define parameters
 265
 266   rin  = 8.;
 267   rout = 14.;
 268   z    = 69.7;
 269   tt   = 2.5;
 270   h1   = TMath::Sqrt(rout * rout + v1/ (tt * kPI)) - rout;
 271   h2   = .4;
 272   t0   = v2/ (h1 * kPI * (h1 + rout * 2.));
 273   t1   = v3/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
 274   t2   = v4/  (h2 * kPI * (h2 + (h1 + rout) * 2.));
 275
 276   // *******Inner slice***********
 277   //     Inner steel wall
 278   par[0] = rin - .02;
 279   par[1] = rin;
 280   par[2] = 1.5;
 281   gMC->Gsvolu("IWR2", "TUBE", idtmed[899], par, 3);
 282   gMC->Gspos("IWR2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
 283   //     Front steel wall
 284   par[0] = rin;
 285   par[1] = rout;
 286   par[2] = .01;
 287   gMC->Gsvolu("FWR2", "TUBE", idtmed[899], par, 3);
 288   gMC->Gspos("FWR2", 1, "ALIC", 0., 0., z + .01, 0, "ONLY");
 289   //     Rear steel wall
 290   gMC->Gsvolu("RWR2", "TUBE", idtmed[899], par, 3);
 291   gMC->Gspos("RWR2", 1, "ALIC", 0., 0., z + 2.99, 0, "ONLY");
 292   //     MCP
 293   par[0] = rin;
 294   par[1] = rout;
 295   par[2] = .07;
 296   gMC->Gsvolu("MPR2", "TUBE", idtmed[900], par, 3);
 297   gMC->Gspos("MPR2", 1, "ALIC", 0., 0., z + 1.57, 0, "ONLY");
 298   //     Silicon plate
 299   par[0] = rin;
 300   par[1] = rout;
 301   par[2] = .019;
 302   gMC->Gsvolu("SPR2", "TUBE", idtmed[901], par, 3);
 303   gMC->Gspos("SPR2", 1, "ALIC", 0., 0., z + 1.719, 0, "ONLY");
 304   //     Summator plate
 305   par[0] = rin;
 306   par[1] = rout;
 307   par[2] = .01;
 308   gMC->Gsvolu("SMR2", "TUBE", idtmed[902], par, 3);
 309   gMC->Gspos("SMR2", 1, "ALIC", 0., 0., z + 2.01, 0, "ONLY");
 310   // *******Outer slice *******
 311   //     Ceramic plate
 312   par[0] = rout;
 313   par[1] = rout + h1;
 314   par[2] = 1.25;
 315   gMC->Gsvolu("CPR2", "TUBE", idtmed[902], par, 3);
 316   gMC->Gspos("CPR2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
 317   //     Covar spring
 318   par[0] = rout;
 319   par[1] = rout + h1;
 320   par[2] = t0 / 2.;
 321   gMC->Gsvolu("C1R2", "TUBE", idtmed[903], par, 3);
 322   gMC->Gspos("C1R2", 1, "ALIC", 0., 0., z + .25 - t0 / 2., 0, "ONLY");
 323   gMC->Gsvolu("C2R2", "TUBE", idtmed[903], par, 3);
 324   gMC->Gspos("C2R2", 1, "ALIC", 0., 0., z + 2.75 + t0 / 2., 0, "ONLY");
 325   //     Getter camera
 326   par[0] = rout + h1;
 327   par[1] = rout + h1 + h2;
 328   par[2] = t1 / 2.;
 329   gMC->Gsvolu("GKR2", "TUBE", idtmed[903], par, 3);
 330   gMC->Gspos("GKR2", 1, "ALIC", 0., 0., z + .25 + t1 / 2., 0, "ONLY");
 331   //     Ceramic slice
 332   par[0] = rout + h1;
 333   par[1] = rout + h1 + h2;
 334   par[2] = t2 / 2.;
 335   gMC->Gsvolu("SCR2", "TUBE", idtmed[902], par, 3);
 336   gMC->Gspos("SCR2", 1, "ALIC", 0., 0., z + .25 + t1 + t2 / 2., 0, "ONLY");
 337   // ******Electronic slice *******
 338   //     Silicon ring
 339   par[0] = rout + h1 + h2;
 340   par[1] = rout + h1 + h2 + 5.;
 341   par[2] = .025;
 342   gMC->Gsvolu("SER2", "TUBE", idtmed[901], par, 3);
 343   gMC->Gspos("SER2", 1, "ALIC", 0., 0., z + 1.57 - .025, 0, "ONLY");
 344   //     Ceramic ring
 345   par[0] = rout + h1 + h2;
 346   par[1] = rout + h1 + h2 + 5.;
 347   par[2] = .025;
 348   gMC->Gsvolu("CER2", "TUBE", idtmed[902], par, 3);
 349   gMC->Gspos("CER2", 1, "ALIC", 0., 0., z + 1.58 + .025, 0, "ONLY");
 350   // ***********************************************************
 351   //       DEFINE LEFT DISK#2 OF FMD
 352   // ***********************************************************
 353
 354   //       Define parameters
 355
 356   rin  = 8.;
 357   rout = 14.;
 358   z    = -69.7;
 359   tt   = 2.5;
 360   h1   = TMath::Sqrt(rout * rout + v1/ (tt * kPI)) - rout;
 361   h2   = .4;
 362   t0   = v2/ (h1 * kPI * (h1 + rout * 2.));
 363   t1   = v3/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
 364   t2   = v4/  (h2 * kPI * (h2 + (h1 + rout) * 2.));
 365
 366   // *******Inner slice***********
 367   //     Inner steel wall
 368   par[0] = rin - .02;
 369   par[1] = rin;
 370   par[2] = 1.5;
 371   gMC->Gsvolu("IWL2", "TUBE", idtmed[899], par, 3);
 372   gMC->Gspos("IWL2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
 373   //     Front steel wall
 374   par[0] = rin;
 375   par[1] = rout;
 376   par[2] = .01;
 377   gMC->Gsvolu("FWL2", "TUBE", idtmed[899], par, 3);
 378   gMC->Gspos("FWL2", 1, "ALIC", 0., 0., z - .01, 0, "ONLY");
 379   //     Rear steel wall
 380   gMC->Gsvolu("RWL2", "TUBE", idtmed[899], par, 3);
 381   gMC->Gspos("RWL2", 1, "ALIC", 0., 0., z - 2.99, 0, "ONLY");
 382   //     MCP
 383   par[0] = rin;
 384   par[1] = rout;
 385   par[2] = .07;
 386   gMC->Gsvolu("MPL2", "TUBE", idtmed[900], par, 3);
 387   gMC->Gspos("MPL2", 1, "ALIC", 0., 0., z - 1.57, 0, "ONLY");
 388   //     Silicon plate
 389   par[0] = rin;
 390   par[1] = rout;
 391   par[2] = .019;
 392   gMC->Gsvolu("SPL2", "TUBE", idtmed[901], par, 3);
 393   gMC->Gspos("SPL2", 1, "ALIC", 0., 0., z - 1.719, 0, "ONLY");
 394   //     Summator plate
 395   par[0] = rin;
 396   par[1] = rout;
 397   par[2] = .01;
 398   gMC->Gsvolu("SML2", "TUBE", idtmed[902], par, 3);
 399   gMC->Gspos("SML2", 1, "ALIC", 0., 0., z - 2.01, 0, "ONLY");
 400   // *******Outer slice *******
 401   //     Ceramic plate
 402   par[0] = rout;
 403   par[1] = rout + h1;
 404   par[2] = 1.25;
 405   gMC->Gsvolu("CPL2", "TUBE", idtmed[902], par, 3);
 406   gMC->Gspos("CPL2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
 407   //     Covar spring
 408   par[0] = rout;
 409   par[1] = rout + h1;
 410   par[2] = t0 / 2.;
 411   gMC->Gsvolu("C1L2", "TUBE", idtmed[903], par, 3);
 412   gMC->Gspos("C1L2", 1, "ALIC", 0., 0., z - .25 + t0 / 2., 0, "ONLY");
 413   gMC->Gsvolu("C2L2", "TUBE", idtmed[903], par, 3);
 414   gMC->Gspos("C2L2", 1, "ALIC", 0., 0., z - 2.75 - t0 / 2., 0, "ONLY");
 415   //     Getter camera
 416   par[0] = rout + h1;
 417   par[1] = rout + h1 + h2;
 418   par[2] = t1 / 2.;
 419   gMC->Gsvolu("GKL2", "TUBE", idtmed[903], par, 3);
 420   gMC->Gspos("GKL2", 1, "ALIC", 0., 0., z - .25 - t1 / 2., 0, "ONLY");
 421   //     Ceramic slice
 422   par[0] = rout + h1;
 423   par[1] = rout + h1 + h2;
 424   par[2] = t2 / 2.;
 425   gMC->Gsvolu("SCL2", "TUBE", idtmed[902], par, 3);
 426   gMC->Gspos("SCL2", 1, "ALIC", 0., 0., z - .25 - t1 - t2 / 2., 0, "ONLY");
 427   // ******Electronic slice *******
 428   //     Silicon ring
 429   par[0] = rout + h1 + h2;
 430   par[1] = rout + h1 + h2 + 5.;
 431   par[2] = .025;
 432   gMC->Gsvolu("SEL2", "TUBE", idtmed[901], par, 3);
 433   gMC->Gspos("SEL2", 1, "ALIC", 0., 0., z - 1.57 + .025, 0, "ONLY");
 434   //     Ceramic ring
 435   par[0] = rout + h1 + h2;
 436   par[1] = rout + h1 + h2 + 5.;
 437   par[2] = .025;
 438   gMC->Gsvolu("CEL2", "TUBE", idtmed[902], par, 3);
 439   gMC->Gspos("CEL2", 1, "ALIC", 0., 0., z - 1.58 - .025, 0, "ONLY");
 440   // ********************************************************
 441   //     DEFINE RIGHT DISK#1  OF FMD
 442   // ********************************************************
 443
 444   //     Define parameters
 445
 446   rin  = 8.;
 447   rout = 17.5;
 448   z    = 42.5;
 449   tt   = 2.5;
 450   h1   = TMath::Sqrt(rout * rout + v1/ (tt * kPI)) - rout;
 451   h2   = .4;
 452   t0   = v2/ (h1 * kPI * (h1 + rout * 2.));
 453   t1   = v3/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
 454   t2   = v4/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
 455
 456   // *******Inner slice***********
 457   //     Inner steel wall
 458   par[0] = rin - .02;
 459   par[1] = rin;
 460   par[2] = 1.5;
 461   gMC->Gsvolu("IWR1", "TUBE", idtmed[899], par, 3);
 462   gMC->Gspos("IWR1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
 463   //     Front steel wall
 464   par[0] = rin;
 465   par[1] = rout;
 466   par[2] = .01;
 467   gMC->Gsvolu("FWR1", "TUBE", idtmed[899], par, 3);
 468   gMC->Gspos("FWR1", 1, "ALIC", 0., 0., z + .01, 0, "ONLY");
 469   //     Rear steel wall
 470   gMC->Gsvolu("RWR1", "TUBE", idtmed[899], par, 3);
 471   gMC->Gspos("RWR1", 1, "ALIC", 0., 0., z + 2.99, 0, "ONLY");
 472   //     MCP
 473   par[0] = rin;
 474   par[1] = rout;
 475   par[2] = .07;
 476   gMC->Gsvolu("MPR1", "TUBE", idtmed[900], par, 3);
 477   gMC->Gspos("MPR1", 1, "ALIC", 0., 0., z + 1.57, 0, "ONLY");
 478   //     Silicon plate
 479   par[0] = rin;
 480   par[1] = rout;
 481   par[2] = .019;
 482   gMC->Gsvolu("SPR1", "TUBE", idtmed[901], par, 3);
 483   gMC->Gspos("SPR1", 1, "ALIC", 0., 0., z + 1.719, 0, "ONLY");
 484   //     Summator plate
 485   par[0] = rin;
 486   par[1] = rout;
 487   par[2] = .01;
 488   gMC->Gsvolu("SMR1", "TUBE", idtmed[902], par, 3);
 489   gMC->Gspos("SMR1", 1, "ALIC", 0., 0., z + 2.01, 0, "ONLY");
 490   // *******Outer slice *******
 491   //     Ceramic plate
 492   par[0] = rout;
 493   par[1] = rout + h1;
 494   par[2] = 1.25;
 495   gMC->Gsvolu("CPR1", "TUBE", idtmed[902], par, 3);
 496   gMC->Gspos("CPR1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
 497   //     Covar spring
 498   par[0] = rout;
 499   par[1] = rout + h1;
 500   par[2] = t0 / 2.;
 501   gMC->Gsvolu("C1R1", "TUBE", idtmed[903], par, 3);
 502   gMC->Gspos("C1R1", 1, "ALIC", 0., 0., z + .25 - t0 / 2., 0, "ONLY");
 503   gMC->Gsvolu("C2R1", "TUBE", idtmed[903], par, 3);
 504   gMC->Gspos("C2R1", 1, "ALIC", 0., 0., z + 2.75 + t0 / 2., 0, "ONLY");
 505   //     Getter camera
 506   par[0] = rout + h1;
 507   par[1] = rout + h1 + h2;
 508   par[2] = t1 / 2.;
 509   gMC->Gsvolu("GKR1", "TUBE", idtmed[903], par, 3);
 510   gMC->Gspos("GKR1", 1, "ALIC", 0., 0., z + .25 + t1 / 2., 0, "ONLY");
 511   //     Ceramic slice
 512   par[0] = rout + h1;
 513   par[1] = rout + h1 + h2;
 514   par[2] = t2 / 2.;
 515   gMC->Gsvolu("SCR1", "TUBE", idtmed[902], par, 3);
 516   gMC->Gspos("SCR1", 1, "ALIC", 0., 0., z + .25 + t1 + t2 / 2., 0, "ONLY");
 517   // ******Electronic slice *******
 518   //     Silicon ring
 519   par[0] = rout + h1 + h2;
 520   par[1] = rout + h1 + h2 + 5.;
 521   par[2] = .025;
 522   gMC->Gsvolu("SER1", "TUBE", idtmed[901], par, 3);
 523   gMC->Gspos("SER1", 1, "ALIC", 0., 0., z + 1.57 - .025, 0, "ONLY");
 524   //     Ceramic ring
 525   par[0] = rout + h1 + h2;
 526   par[1] = rout + h1 + h2 + 5.;
 527   par[2] = .025;
 528   gMC->Gsvolu("CER1", "TUBE", idtmed[902], par, 3);
 529   gMC->Gspos("CER1", 1, "ALIC", 0., 0., z + 1.58 + .025, 0, "ONLY");
 530   // ***********************************************************
 531   //     DEFINE LEFT DISK#1 OF FMD
 532   // ***********************************************************
 533
 534   //     Define parameters
 535
 536   rin  = 8.;
 537   rout = 17.5;
 538   z    = -42.5;
 539   tt   = 2.5;
 540   h1   = TMath::Sqrt(rout * rout + v1/ (tt * kPI)) - rout;
 541   h2   = .4;
 542   t0   = v2/ (h1 * kPI * (h1 + rout * 2.));
 543   t1   = v3/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
 544   t2   = v4/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
 545
 546   // *******Inner slice***********
 547   //     Inner steel wall
 548   par[0] = rin - .02;
 549   par[1] = rin;
 550   par[2] = 1.5;
 551   gMC->Gsvolu("IWL1", "TUBE", idtmed[899], par, 3);
 552   gMC->Gspos("IWL1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
 553   //     Front steel wall
 554   par[0] = rin;
 555   par[1] = rout;
 556   par[2] = .01;
 557   gMC->Gsvolu("FWL1", "TUBE", idtmed[899], par, 3);
 558   gMC->Gspos("FWL1", 1, "ALIC", 0., 0., z - .01, 0, "ONLY");
 559   //     Rear steel wall
 560   gMC->Gsvolu("RWL1", "TUBE", idtmed[899], par, 3);
 561   gMC->Gspos("RWL1", 1, "ALIC", 0., 0., z - 2.99, 0, "ONLY");
 562   //     MCP
 563   par[0] = rin;
 564   par[1] = rout;
 565   par[2] = .07;
 566   gMC->Gsvolu("MPL1", "TUBE", idtmed[900], par, 3);
 567   gMC->Gspos("MPL1", 1, "ALIC", 0., 0., z - 1.57, 0, "ONLY");
 568   //     Silicon plate
 569   par[0] = rin;
 570   par[1] = rout;
 571   par[2] = .019;
 572   gMC->Gsvolu("SPL1", "TUBE", idtmed[901], par, 3);
 573   gMC->Gspos("SPL1", 1, "ALIC", 0., 0., z - 1.719, 0, "ONLY");
 574   //     Summator plate
 575   par[0] = rin;
 576   par[1] = rout;
 577   par[2] = .01;
 578   gMC->Gsvolu("SML1", "TUBE", idtmed[902], par, 3);
 579   gMC->Gspos("SML1", 1, "ALIC", 0., 0., z - 2.01, 0, "ONLY");
 580   // *******Outer slice *******
 581   //     Ceramic plate
 582   par[0] = rout;
 583   par[1] = rout + h1;
 584   par[2] = 1.25;
 585   gMC->Gsvolu("CPL1", "TUBE", idtmed[902], par, 3);
 586   gMC->Gspos("CPL1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
 587   //     Covar spring
 588   par[0] = rout;
 589   par[1] = rout + h1;
 590   par[2] = t0 / 2.;
 591   gMC->Gsvolu("C1L1", "TUBE", idtmed[903], par, 3);
 592   gMC->Gspos("C1L1", 1, "ALIC", 0., 0., z - .25 + t0 / 2., 0, "ONLY");
 593   gMC->Gsvolu("C2L1", "TUBE", idtmed[903], par, 3);
 594   gMC->Gspos("C2L1", 1, "ALIC", 0., 0., z - 2.75 - t0 / 2., 0, "ONLY");
 595   //     Getter camera
 596   par[0] = rout + h1;
 597   par[1] = rout + h1 + h2;
 598   par[2] = t1 / 2.;
 599   gMC->Gsvolu("GKL1", "TUBE", idtmed[903], par, 3);
 600   gMC->Gspos("GKL1", 1, "ALIC", 0., 0., z - .25 - t1 / 2., 0, "ONLY");
 601   //     Ceramic slice
 602   par[0] = rout + h1;
 603   par[1] = rout + h1 + h2;
 604   par[2] = t2 / 2.;
 605   gMC->Gsvolu("SCL1", "TUBE", idtmed[902], par, 3);
 606   gMC->Gspos("SCL1", 1, "ALIC", 0., 0., z - .25 - t1 - t2 / 2., 0, "ONLY");
 607   // ******Electronic slice *******
 608   //     Silicon ring
 609   par[0] = rout + h1 + h2;
 610   par[1] = rout + h1 + h2 + 5.;
 611   par[2] = .025;
 612   gMC->Gsvolu("SEL1", "TUBE", idtmed[901], par, 3);
 613   gMC->Gspos("SEL1", 1, "ALIC", 0., 0., z - 1.57 + .025, 0, "ONLY");
 614   //     Ceramic ring
 615   par[0] = rout + h1 + h2;
 616   par[1] = rout + h1 + h2 + 5.;
 617   par[2] = .025;
 618   gMC->Gsvolu("CEL1", "TUBE", idtmed[902], par, 3);
 619   gMC->Gspos("CEL1", 1, "ALIC", 0., 0., z - 1.58 - .025, 0, "ONLY");
 620   // ***********************************************************
 621   //     DEFINE LEFT DISK#4 OF FMD
 622   // ***********************************************************
 623
 624   //     Define parameters
 625
 626   rin  = 4.2;
 627   rout = 13.;
 628   z    = -229.5;
 629   tt   = 2.5;
 630   h1   = TMath::Sqrt(rout * rout + v1/ (tt * kPI)) - rout;
 631   h2   = .4;
 632   t0   = v2/ (h1 * kPI * (h1 + rout * 2.));
 633   t1   = v3/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
 634   t2   = v4/ (h2 * kPI * (h2 + (h1 + rout) * 2.));
 635
 636   // *******Inner slice***********
 637   //     Inner steel wall
 638   par[0] = rin - .02;
 639   par[1] = rin;
 640   par[2] = 1.5;
 641   gMC->Gsvolu("IWL4", "TUBE", idtmed[899], par, 3);
 642   gMC->Gspos("IWL4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
 643   //     Front steel wall
 644   par[0] = rin;
 645   par[1] = rout;
 646   par[2] = .01;
 647   gMC->Gsvolu("FWL4", "TUBE", idtmed[899], par, 3);
 648   gMC->Gspos("FWL4", 1, "ALIC", 0., 0., z - .01, 0, "ONLY");
 649   //     Rear steel wall
 650   gMC->Gsvolu("RWL4", "TUBE", idtmed[899], par, 3);
 651   gMC->Gspos("RWL4", 1, "ALIC", 0., 0., z - 2.99, 0, "ONLY");
 652   //     MCP
 653   par[0] = rin;
 654   par[1] = rout;
 655   par[2] = .07;
 656   gMC->Gsvolu("MPL4", "TUBE", idtmed[900], par, 3);
 657   gMC->Gspos("MPL4", 1, "ALIC", 0., 0., z - 1.57, 0, "ONLY");
 658   //     Silicon plate
 659   par[0] = rin;
 660   par[1] = rout;
 661   par[2] = .019;
 662   gMC->Gsvolu("SPL4", "TUBE", idtmed[901], par, 3);
 663   gMC->Gspos("SPL4", 1, "ALIC", 0., 0., z - 1.719, 0, "ONLY");
 664   //     Summator plate
 665   par[0] = rin;
 666   par[1] = rout;
 667   par[2] = .01;
 668   gMC->Gsvolu("SML4", "TUBE", idtmed[902], par, 3);
 669   gMC->Gspos("SML4", 1, "ALIC", 0., 0., z - 2.01, 0, "ONLY");
 670   // *******Outer slice *******
 671   //     Ceramic plate
 672   par[0] = rout;
 673   par[1] = rout + h1;
 674   par[2] = 1.25;
 675   gMC->Gsvolu("CPL4", "TUBE", idtmed[902], par, 3);
 676   gMC->Gspos("CPL4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
 677   //     Covar spring
 678   par[0] = rout;
 679   par[1] = rout + h1;
 680   par[2] = t0 / 2.;
 681   gMC->Gsvolu("C1L4", "TUBE", idtmed[903], par, 3);
 682   gMC->Gspos("C1L4", 1, "ALIC", 0., 0., z - .25 + t0 / 2., 0, "ONLY");
 683   gMC->Gsvolu("C2L4", "TUBE", idtmed[903], par, 3);
 684   gMC->Gspos("C2L4", 1, "ALIC", 0., 0., z - 2.75 - t0 / 2., 0, "ONLY");
 685   //       Getter camera
 686   par[0] = rout + h1;
 687   par[1] = rout + h1 + h2;
 688   par[2] = t1 / 2.;
 689   gMC->Gsvolu("GKL4", "TUBE", idtmed[903], par, 3);
 690   gMC->Gspos("GKL4", 1, "ALIC", 0., 0., z - .25 - t1 / 2., 0, "ONLY");
 691   //     Ceramic slice
 692   par[0] = rout + h1;
 693   par[1] = rout + h1 + h2;
 694   par[2] = t2 / 2.;
 695   gMC->Gsvolu("SCL4", "TUBE", idtmed[902], par, 3);
 696   gMC->Gspos("SCL4", 1, "ALIC", 0., 0., z - .25 - t1 - t2 / 2., 0, "ONLY");
 697   // ******Electronic slice *******
 698   //     Silicon ring
 699   par[0] = rout + h1 + h2;
 700   par[1] = rout + h1 + h2 + 5.;
 701   par[2] = .025;
 702   gMC->Gsvolu("SEL4", "TUBE", idtmed[901], par, 3);
 703   gMC->Gspos("SEL4", 1, "ALIC", 0., 0., z - 1.57 + .025, 0, "ONLY");
 704   //     Ceramic ring
 705   par[0] = rout + h1 + h2;
 706   par[1] = rout + h1 + h2 + 5.;
 707   par[2] = .025;
 708   gMC->Gsvolu("CEL4", "TUBE", idtmed[902], par, 3);
 709   gMC->Gspos("CEL4", 1, "ALIC", 0., 0., z - 1.58 - .025, 0, "ONLY");
 710 }
 711
 712 //_____________________________________________________________________________
 713 void AliFMDv1::DrawModule()
 714 {
 715   //
 716   // Draw a shaded view of the FMD version 1
 717   //
 718
 719
 720   // Set everything unseen
 721   gMC->Gsatt("*", "seen", -1);
 722   //
 723   // Set ALIC mother transparent
 724   gMC->Gsatt("ALIC","SEEN",0);
 725   //
 726   // Set the volumes visible
 727   gMC->Gsatt("IWR3","seen",1);
 728   gMC->Gsatt("FWR3","seen",1);
 729   gMC->Gsatt("RWR3","seen",1);
 730   gMC->Gsatt("MPR3","seen",1);
 731   gMC->Gsatt("SPR3","seen",1);
 732   gMC->Gsatt("SMR3","seen",1);
 733   gMC->Gsatt("CPR3","seen",1);
 734   gMC->Gsatt("C1R3","seen",1);
 735   gMC->Gsatt("C2R3","seen",1);
 736   gMC->Gsatt("GKR3","seen",1);
 737   gMC->Gsatt("SCR3","seen",1);
 738   gMC->Gsatt("SER3","seen",1);
 739   gMC->Gsatt("CER3","seen",1);
 740   gMC->Gsatt("IWL3","seen",1);
 741   gMC->Gsatt("FWL3","seen",1);
 742   gMC->Gsatt("RWL3","seen",1);
 743   gMC->Gsatt("MPL3","seen",1);
 744   gMC->Gsatt("SPL3","seen",1);
 745   gMC->Gsatt("SML3","seen",1);
 746   gMC->Gsatt("CPL3","seen",1);
 747   gMC->Gsatt("C1L3","seen",1);
 748   gMC->Gsatt("C2L3","seen",1);
 749   gMC->Gsatt("GKL3","seen",1);
 750   gMC->Gsatt("SCL3","seen",1);
 751   gMC->Gsatt("SEL3","seen",1);
 752   gMC->Gsatt("CEL3","seen",1);
 753   gMC->Gsatt("IWR2","seen",1);
 754   gMC->Gsatt("FWR2","seen",1);
 755   gMC->Gsatt("RWR2","seen",1);
 756   gMC->Gsatt("MPR2","seen",1);
 757   gMC->Gsatt("SPR2","seen",1);
 758   gMC->Gsatt("SMR2","seen",1);
 759   gMC->Gsatt("CPR2","seen",1);
 760   gMC->Gsatt("C1R2","seen",1);
 761   gMC->Gsatt("C2R2","seen",1);
 762   gMC->Gsatt("GKR2","seen",1);
 763   gMC->Gsatt("SCR2","seen",1);
 764   gMC->Gsatt("SER2","seen",1);
 765   gMC->Gsatt("CER2","seen",1);
 766   gMC->Gsatt("IWL2","seen",1);
 767   gMC->Gsatt("FWL2","seen",1);
 768   gMC->Gsatt("RWL2","seen",1);
 769   gMC->Gsatt("MPL2","seen",1);
 770   gMC->Gsatt("SPL2","seen",1);
 771   gMC->Gsatt("SML2","seen",1);
 772   gMC->Gsatt("CPL2","seen",1);
 773   gMC->Gsatt("C1L2","seen",1);
 774   gMC->Gsatt("C2L2","seen",1);
 775   gMC->Gsatt("GKL2","seen",1);
 776   gMC->Gsatt("SCL2","seen",1);
 777   gMC->Gsatt("SEL2","seen",1);
 778   gMC->Gsatt("CEL2","seen",1);
 779   gMC->Gsatt("IWR1","seen",1);
 780   gMC->Gsatt("FWR1","seen",1);
 781   gMC->Gsatt("RWR1","seen",1);
 782   gMC->Gsatt("MPR1","seen",1);
 783   gMC->Gsatt("SPR1","seen",1);
 784   gMC->Gsatt("SMR1","seen",1);
 785   gMC->Gsatt("CPR1","seen",1);
 786   gMC->Gsatt("C1R1","seen",1);
 787   gMC->Gsatt("C2R1","seen",1);
 788   gMC->Gsatt("GKR1","seen",1);
 789   gMC->Gsatt("SCR1","seen",1);
 790   gMC->Gsatt("SER1","seen",1);
 791   gMC->Gsatt("CER1","seen",1);
 792   gMC->Gsatt("IWL1","seen",1);
 793   gMC->Gsatt("FWL1","seen",1);
 794   gMC->Gsatt("RWL1","seen",1);
 795   gMC->Gsatt("MPL1","seen",1);
 796   gMC->Gsatt("SPL1","seen",1);
 797   gMC->Gsatt("SML1","seen",1);
 798   gMC->Gsatt("CPL1","seen",1);
 799   gMC->Gsatt("C1L1","seen",1);
 800   gMC->Gsatt("C2L1","seen",1);
 801   gMC->Gsatt("GKL1","seen",1);
 802   gMC->Gsatt("SCL1","seen",1);
 803   gMC->Gsatt("SEL1","seen",1);
 804   gMC->Gsatt("CEL1","seen",1);
 805   gMC->Gsatt("IWL4","seen",1);
 806   gMC->Gsatt("FWL4","seen",1);
 807   gMC->Gsatt("RWL4","seen",1);
 808   gMC->Gsatt("MPL4","seen",1);
 809   gMC->Gsatt("SPL4","seen",1);
 810   gMC->Gsatt("SML4","seen",1);
 811   gMC->Gsatt("CPL4","seen",1);
 812   gMC->Gsatt("C1L4","seen",1);
 813   gMC->Gsatt("C2L4","seen",1);
 814   gMC->Gsatt("GKL4","seen",1);
 815   gMC->Gsatt("SCL4","seen",1);
 816   gMC->Gsatt("SEL4","seen",1);
 817   gMC->Gsatt("CEL4","seen",1);
 818   //
 819   gMC->Gdopt("hide", "on");
 820   gMC->Gdopt("shad", "on");
 821   gMC->Gsatt("*", "fill", 7);
 822   gMC->SetClipBox(".");
 823   gMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
 824   gMC->DefaultRange();
 825   gMC->Gdraw("alic", 40, 30, 0, 6, 9, .08, .08);
 826   gMC->Gdhead(1111, "Forward Multiplicity Detector version 1");
 827   gMC->Gdman(13, 9, "MAN");
 828 }
 829
 830 //_____________________________________________________________________________
 831 void AliFMDv1::CreateMaterials()
 832 {
 833   //
 834   // Create materials for version 1 of FMD
 835   //
 836
 837   //     Steel for chamber
 838   Float_t ast[4]  = { 55.847,58.71,51.996,47.9 };
 839   Float_t zst[4]  = { 26.,28.,24.,22. };
 840   Float_t wst[4]  = { .6,.18,.12,.1 };
 841   //     Lead glass for MCP
 842   Float_t amcp[3] = { 15.9994,28.086,207.19 };
 843   Float_t zmcp[3] = { 8.,14.,82. };
 844   Float_t wmcp[3] = { .43,.375,.195 };
 845   //     Ceramic for plates
 846   Float_t acer[2] = { 26.98,15.9994 };
 847   Float_t zcer[2] = { 13.,8. };
 848   Float_t wcer[2] = { .4,.6 };
 849
 850   Float_t epsil, stmin, deemax, tmaxfd, stemax;
 851
 852   Int_t   ISXFLD = gAlice->Field()->Integ();
 853   Float_t SXMGMX = gAlice->Field()->Max();
 854
 855   //     Silicon
 856
 857   AliMaterial(2, "SI$", 28.086, 14., 2.33, 9.36, 45.5);
 858
 859   //     Covar
 860
 861   AliMaterial(4, "FE$", 55.85, 26., 8.3, 1.67, 15.9);
 862
 863   //     Steel
 864
 865   AliMixture(0, "FMD_STEEL$", ast, zst, 7.9, 4, wst);
 866
 867   //     Lead glass
 868
 869   AliMixture(1, "FMD_MCP$", amcp, zmcp, 1.55, 3, wmcp);
 870
 871   //     Ceramic
 872
 873   AliMixture(3, "FMD_CERAMIC$", acer, zcer, 3.96, -2, wcer);
 874   // *******************************************************
 875   //     Defines tracking media parameters.
 876   // *******************************************************
 877   epsil  = .001; // Tracking precision, DLS
 878   stemax = -1.;  // Maximum displacement for multiple scattering
 879   tmaxfd = -20.; // Maximum angle due to field deflection
 880   deemax = -.3;  // Maximum fractional energy loss, DLS
 881   stmin  = -.8;
 882   // ********************************************************
 883   AliMedium(0, "STEEL_L3          ", 0, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
 884   AliMedium(1, "LEAD_GLASS_L3     ", 1, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
 885   AliMedium(2, "SILICON_L3        ", 2, 1, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
 886   AliMedium(3, "CERAMIC_L3        ", 3, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
 887   AliMedium(4, "COVAR_L3          ", 4, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
 888 }
 889