FMD/AliFMDv1.cxx

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