FMD/AliFMDv0.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 // Forward Multiplicity detector based on Silicon version 0        //
  18 //
  19 //Begin Html
  20 /*
  21 <img src="gif/AliFMDv0Class.gif">
  22 */
  23 //End Html
  24 //                                                                  //
  25 //                                                                  //
  26 //////////////////////////////////////////////////////////////////////
  27
  28 #include <TMath.h>
  29 #include <TGeometry.h>
  30 #include <TTUBE.h>
  31 #include <TNode.h>
  32 #include <TBranch.h>
  33 #include <TLorentzVector.h>
  34 #include "AliFMDv0.h"
  35 #include "AliRun.h"
  36 #include <Riostream.h>
  37 #include "AliRndm.h"
  38 #include "AliMagF.h"
  39 #include "AliFMDhit.h"
  40 #include <stdlib.h>
  41
  42 //class TGeant3;
  43 ClassImp(AliFMDv0)
  44
  45 //--------------------------------------------------------------------
  46 AliFMDv0::AliFMDv0(const char *name, const char *title):
  47  AliFMD(name,title)
  48 {
  49   //
  50   // Standart constructor for Forward Multiplicity Detector version 0
  51   //
  52   fIdSens1=0;
  53   fIdSens2=0;
  54   fIdSens3=0;
  55   fIdSens4=0;
  56   fIdSens5=0;
  57 //  setBufferSize(128000);
  58 }
  59 //-------------------------------------------------------------------------
  60
  61 void AliFMDv0::CreateGeometry()
  62 {
  63  //
  64   // Create the geometry of Forward Multiplicity Detector version 0
  65   //Detector consists of 5 volumes:
  66   // 1st covered pseudorapidity interval from 3.58 to 2.03
  67   // and placed on 62.8cm in Z-direction;
  68   // 2nd - from 1.94 to 1.51 and Z=75.2 cm;
  69   // 3d  - the same pseudorapidity interval as the 1st
  70   // but on the other side from the interaction point z=-62.8cm;
  71   // 4th - simmetricaly with the 2nd :
  72   // pseudorapidity from 1.94 to 1.51, Z=-75.2cm
  73   // 5th - from 3.71 to 5.28, Z=-345cm
  74   // Each part has 300mkm Si (sensetive area, detector itself),
  75   // 0.75cm of plastic simulated electronics material,
  76   // Al support ring 2cm thickness and 1cm width placed on
  77   // the outer radius of each Si disk;
  78    //
  79   // begin Html
  80   /*
  81    <img src="gif/AliFMDv0.gif">
  82    */
  83   //
  84
  85   Int_t *idtmed = fIdtmed->GetArray();
  86   Int_t ifmd;
  87   Int_t idrotm[999];
  88   Float_t zfmd,par[3];
  89   char name[5], nameSi[5], nameSector[5], nameRing[5];
  90
  91   Float_t rin[6], rout[6],zpos;
  92
  93   Float_t etain[5]= {3.40, 2.29, 3.68, 2.29, 5.09};
  94   Float_t etaout[6]={2.01, 1.70, 2.28, 1.70, 3.68};
  95   //  Float_t z[6]={64., 85., -64., -85., -270., -630};
  96   Float_t z[6]={62.8, 75.2, -83.4, -75.2, -340.};
  97   Float_t zDet=0.03;
  98   Float_t zElectronic=0.1;
  99   Float_t zSupport=1.;
 100
 101   Float_t zFMD=1.;
 102 //-------------------------------------------------------------------
 103  //  FMD
 104  //------------------------------------------------------------------
 105         cout<<" !!!!!!!!!!!New FMD geometry !!!!!!!!!"<<endl;
 106
 107   AliMatrix(idrotm[901], 90, 0, 90, 90, 180, 0);
 108
 109   //  gMC->Gsvolu("GSI","TUBE", idtmed[1], par, 0);
 110   gMC->Gsvolu("GEL ","TUBE", idtmed[4], par, 0);
 111   gMC->Gsvolu("GSUP","TUBE", idtmed[2], par, 0);
 112
 113   for (ifmd =0; ifmd < 5; ifmd++){
 114
 115     sprintf(name,"FMD%d",ifmd+1);
 116     sprintf(nameSi,"GSI%d",ifmd+1);
 117     sprintf(nameSector,"GSC%d",ifmd+1);
 118     sprintf(nameRing,"GRN%d",ifmd+1);
 119     printf(name,nameSi);
 120
 121     zfmd=TMath::Abs(z[ifmd]);
 122     AliFMD::Eta2Radius(etain[ifmd],zfmd,&rin[ifmd]);
 123     AliFMD::Eta2Radius(etaout[ifmd],zfmd,&rout[ifmd]);
 124
 125     par[0]=rin[ifmd]; // pipe size
 126     par[1]=rout[ifmd];
 127     par[2]=zFMD/2;
 128     gMC->Gsvolu(name,"TUBE", idtmed[3], par, 3);
 129     gMC->Gsvolu(nameSi,"TUBE", idtmed[1], par, 0);
 130
 131      if (z[ifmd] < 0){
 132       gMC->Gspos(name,1,"ALIC",0,0,z[ifmd],0, "ONLY");}
 133     else {
 134       gMC->Gspos(name,1,"ALIC",0,0,z[ifmd],idrotm[901], "ONLY");}
 135   //Silicon detector
 136     par[2]=zDet/2;
 137     zpos=zFMD/2 -par[2];
 138     gMC->Gsposp(nameSi,ifmd+1,name,0,0,zpos,0, "ONLY",par,3);
 139      //Granularity
 140     if(ifmd==1||ifmd==3)
 141       {
 142         gMC->Gsdvn(nameSector, nameSi , fSectorsSi2, 2);
 143         gMC->Gsdvn(nameRing, nameSector, fRingsSi2, 1);
 144       }
 145     else
 146       {
 147         gMC->Gsdvn(nameSector, nameSi , fSectorsSi1, 2);
 148         gMC->Gsdvn(nameRing, nameSector , fRingsSi1, 1);
 149       }
 150
 151     //Plastic slice for electronics
 152     par[2]=zElectronic/2;
 153     zpos=zpos-zDet/2-par[2];
 154     gMC->Gsposp("GEL ",ifmd+1,name,0,0,zpos,0, "ONLY",par,3);
 155
 156    //Simple Al support
 157    par[1]=rout[ifmd];
 158    par[0]=rout[ifmd]-2;
 159    par[2]=zSupport/2;
 160    zpos=zpos-zElectronic/2-par[2];
 161    //   gMC->Gsposp("GSUP",ifmd+1,name,0,0,zpos,0, "ONLY",par,3);
 162
 163
 164   }
 165
 166 }
 167
 168
 169 //------------------------------------------------------------------------
 170 void AliFMDv0::CreateMaterials()
 171 {
 172  Int_t isxfld   = gAlice->Field()->Integ();
 173  Float_t sxmgmx = gAlice->Field()->Max();
 174
 175  // Plastic CH
 176  Float_t aPlastic[2]={1.01,12.01};
 177  Float_t zPlastic[2]={1,6};
 178  Float_t wPlastic[2]={1,1};
 179  Float_t denPlastic=1.03;
 180    //
 181
 182  //*** Definition Of avaible FMD materials ***
 183  AliMaterial(0, "Si chip$", 28.0855,14.,2.33,9.36,999);
 184  AliMaterial(1, "Al supprt$", 26.980,13.,2.70,8.9,999);
 185  AliMaterial(2, "FMD Air$", 14.61, 7.3, .001205, 30423.,999);
 186  AliMixture( 5, "Plastic$",aPlastic,zPlastic,denPlastic,-2,wPlastic);
 187
 188
 189 //**
 190  AliMedium(1, "Si chip$", 0, 1, isxfld, sxmgmx, 1., .001, 1., .001, .001);
 191  AliMedium(2, "Al support$", 1, 0, isxfld, sxmgmx, 1., .001, 1., .001, .001);
 192  AliMedium(3, "FMD air$", 2, 0, isxfld, sxmgmx, 1., .001, 1., .001, .001);
 193  AliMedium(4, "Plastic$", 5, 0,isxfld, sxmgmx,  10., .01, 1., .003, .003);
 194
 195
 196
 197 }
 198 //---------------------------------------------------------------------
 199 void AliFMDv0::DrawDetector()
 200 {
 201 //
 202 // Draw a shaded view of the Forward multiplicity detector version 0
 203 //
 204
 205
 206 //Set ALIC mother transparent
 207 gMC->Gsatt("ALIC","SEEN",0);
 208 //
 209 //Set volumes visible
 210 gMC->Gsatt("FMD0","SEEN",1);
 211 gMC->Gsatt("FMD1","SEEN",1);
 212 gMC->Gsatt("FMD2","SEEN",1);
 213 gMC->Gsatt("FMD3","SEEN",1);
 214 gMC->Gsatt("FMD4","SEEN",1);
 215 gMC->Gsatt("FMD5","SEEN",1);
 216
 217 //
 218 gMC->Gdopt("hide","on");
 219 gMC->Gdopt("shad","on");
 220 gMC->SetClipBox(".");
 221 gMC->SetClipBox("*",0,1000,-1000,1000,-1000,1000);
 222 gMC->DefaultRange();
 223 gMC->Gdraw("alic",40,30,0,12,9.5,.2,0.2);
 224 gMC->Gdhead(1111,"Forward multiplicity detector");
 225 gMC->Gdopt("hide","off");
 226 }
 227 //-------------------------------------------------------------------
 228 void AliFMDv0::Init()
 229 {
 230 // Initialises version 0 of the Forward Multiplicity Detector
 231 //
 232 AliFMD::Init();
 233 fIdSens1=gMC->VolId("GRN1");
 234 fIdSens2=gMC->VolId("GRN2");
 235 fIdSens3=gMC->VolId("GRN3");
 236 fIdSens4=gMC->VolId("GRN4");
 237 fIdSens5=gMC->VolId("GRN5");
 238 if (fDebug) printf("*** FMD version 0 initialized ***\n");
 239 }
 240
 241 //-------------------------------------------------------------------
 242
 243 void AliFMDv0::StepManager()
 244 {
 245   //
 246   // Called for every step in the Forward Multiplicity Detector
 247   //
 248 }
 249