MUON/AliMUONv0.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 /* $Id$ */
  17
  18 /////////////////////////////////////////////////////////
  19 //  Manager and hits classes for set:MUON version 0    //
  20 /////////////////////////////////////////////////////////
  21
  22 #include <TLorentzVector.h>
  23 #include <TVirtualMC.h>
  24
  25 #include "AliConst.h"
  26 #include "AliMUONChamber.h"
  27 #include "AliMUONChamberGeometry.h"
  28 #include "AliMUONConstants.h"
  29 #include "AliMUONFactory.h"
  30 #include "AliMUONHit.h"
  31 #include "AliMUONv0.h"
  32 #include "AliMagF.h"
  33 #include "AliRun.h"
  34 #include "AliMC.h"
  35
  36 ClassImp(AliMUONv0)
  37
  38 //___________________________________________
  39 AliMUONv0::AliMUONv0()
  40   : AliMUON()
  41 {
  42 // Constructor
  43     fChambers = 0;
  44 }
  45
  46 //___________________________________________
  47 AliMUONv0::AliMUONv0(const char *name, const char *title)
  48   : AliMUON(name,title)
  49 {
  50 // Constructor
  51     AliMUONFactory factory;
  52     factory.Build(this, title);
  53 }
  54
  55 //___________________________________________
  56 void AliMUONv0::CreateGeometry()
  57 {
  58 // Creates coarse geometry for hit density simulations
  59     Int_t *idtmed = fIdtmed->GetArray()-1099;
  60 //
  61      Float_t zpos, dAlu, tpar[3];
  62      Int_t idAir=idtmed[1100];
  63      Int_t idAlu=idtmed[1103];
  64
  65      AliMUONChamber *iChamber;
  66      // Loop over all chambers (tracking and trigger)
  67      for (Int_t ch = 0; ch < AliMUONConstants::NCh(); ch++) {
  68          char alu[8];
  69          char gas[8];
  70
  71          iChamber=(AliMUONChamber*) (*fChambers)[ch];
  72          // Z of the chamber
  73          zpos=iChamber->Z();
  74          dAlu=iChamber->DAlu();
  75          if (ch < AliMUONConstants::NTrackingCh()) {
  76            // tracking chambers
  77              sprintf(alu,"SA0%1d",ch);
  78              sprintf(gas,"SG0%1d",ch);
  79          } else {
  80            // trigger chambers
  81              sprintf(alu,"SA%2d",ch);
  82              sprintf(gas,"SG%2d",ch);
  83          }
  84 //
  85          tpar[0] = iChamber->RInner();
  86          tpar[1] = iChamber->ROuter();
  87          tpar[2] = (dAlu+0.2)/2.;
  88          if (ch !=4 && ch !=5) {
  89              gMC->Gsvolu(alu, "TUBE", idAlu, tpar, 3);
  90              tpar[2] = 0.1;
  91              gMC->Gsvolu(gas, "TUBE", idAir, tpar, 3);
  92          } else {
  93              gMC->Gsvolu(alu, "TUBE", idAlu, tpar, 3);
  94              tpar[2] = 0.1;
  95              gMC->Gsvolu(gas, "TUBE", idAir, tpar, 3);
  96          }
  97          gMC->Gspos(gas, 1, alu,  0., 0., 0., 0, "ONLY");
  98          if (ch == 4 || ch ==5) {
  99              if (gMC->VolId("DDIP")) {
 100                  gMC->Gspos(alu, 1, "DDIP", 0., 0., zpos, 0, "ONLY");
 101              } else {
 102                  gMC->Gspos(alu, 1, "ALIC", 0., 0., zpos, 0, "ONLY");
 103              }
 104          } else {
 105              gMC->Gspos(alu, 1, "ALIC", 0., 0., zpos, 0, "ONLY");
 106          }
 107      }
 108 }
 109
 110 //___________________________________________
 111 void AliMUONv0::CreateMaterials()
 112 {
 113 // Creates materials for coarse geometry
 114 // Air
 115   Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
 116   Float_t zAir[4]={6.,7.,8.,18.};
 117   Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
 118   Float_t dAir = 1.20479E-3;
 119   AliMixture(15, "AIR$      ", aAir,  zAir, dAir,4, wAir);
 120   //  AliMaterial(15, "AIR$      ", 14.61,  7.3, .001205, 30423.24, 67500);
 121   AliMaterial( 9, "ALUMINIUM$", 26.98, 13. , 2.7, 8.9, 37.2);
 122
 123   Float_t epsil  = .001; // Tracking precision,
 124   Float_t stemax = -1.;  // Maximum displacement for multiple scat
 125   Float_t tmaxfd = -20.; // Maximum angle due to field deflection
 126   Float_t deemax = -.3;  // Maximum fractional energy loss, DLS
 127   Float_t stmin  = -.8;
 128   Int_t isxfld   = gAlice->Field()->Integ();
 129   Float_t sxmgmx = gAlice->Field()->Max();
 130
 131   //
 132   //    Air
 133   AliMedium(1, "AIR_CH_US         ", 15, 1, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 134   AliMedium(4, "ALU_CH_US         ",  9, 0, isxfld, sxmgmx, tmaxfd, fMaxStepAlu,
 135             fMaxDestepAlu, epsil, stmin);
 136
 137 }
 138
 139 void AliMUONv0::Init()
 140 {
 141    //
 142    // Initialize Tracking Chambers
 143    //
 144     char vName[8];
 145     printf("\n\n\n Start Init for version 0 - CPC chamber type\n\n\n");
 146     for (Int_t i=0; i<AliMUONConstants::NCh(); i++) {
 147 // Initialise chamber
 148         ((AliMUONChamber*) (*fChambers)[i])->Init();
 149 // Set sensitive volume Id
 150         if (i < AliMUONConstants::NTrackingCh()) {
 151             // tracking chambers
 152             sprintf(vName,"SG0%1d",i);
 153         } else {
 154             // trigger chambers
 155             sprintf(vName,"SG%2d",i);
 156         }
 157         //((AliMUONChamber*) (*fChambers)[i])->SetGid(gMC->VolId(vName));
 158         ((AliMUONChamber*) (*fChambers)[i])
 159            ->GetGeometry()->SetSensitiveVolume(gMC->VolId(vName));
 160     }
 161 }
 162
 163 void AliMUONv0::StepManager()
 164 {
 165 //
 166 // Step manager for hit density simulations
 167   Int_t          copy, id;
 168   static Int_t   idvol;
 169   static Int_t   vol[2];
 170   Int_t          ipart;
 171   TLorentzVector pos;
 172   TLorentzVector mom;
 173   Float_t        theta,phi;
 174
 175   //  modifs perso
 176   static Float_t hits[15];
 177
 178   TClonesArray &lhits = *fHits;
 179   //
 180   // Only gas gap inside chamber
 181   // Tag chambers and record hits when track enters
 182   idvol=-1;
 183   id=gMC->CurrentVolID(copy);
 184
 185     for (Int_t i=1; i<=AliMUONConstants::NCh(); i++) {
 186       //if(id==((AliMUONChamber*)(*fChambers)[i-1])->GetGid()){
 187       if ( ((AliMUONChamber*)(*fChambers)[i-1])->IsSensId(id) ) {
 188           vol[0]=i;
 189           idvol=i-1;
 190       }
 191     }
 192     if (idvol == -1) return;
 193   //
 194   // Get current particle id (ipart), track position (pos)  and momentum (mom)
 195   gMC->TrackPosition(pos);
 196   gMC->TrackMomentum(mom);
 197
 198   ipart  = gMC->TrackPid();
 199   //
 200   // record hits when track enters ...
 201   if( !(gMC->TrackCharge()) ) return;
 202   if( gMC->IsTrackEntering()) {
 203 //      printf("\n Particle entering %f %f %f", pos[0], pos[1], pos[2] );
 204
 205       Double_t tc = mom[0]*mom[0]+mom[1]*mom[1];
 206       Double_t rt = TMath::Sqrt(tc);
 207       theta   = Float_t(TMath::ATan2(rt,Double_t(mom[2])))*kRaddeg;
 208       phi     = Float_t(TMath::ATan2(Double_t(mom[1]),Double_t(mom[0])))*kRaddeg;
 209       hits[0] = Float_t(ipart);             // Geant3 particle type
 210       hits[1] = pos[0];                     // X-position for hit
 211       hits[2] = pos[1];                     // Y-position for hit
 212       hits[3] = pos[2];                     // Z-position for hit
 213       hits[4] = theta;                      // theta angle of incidence
 214       hits[5] = phi;                        // phi angle of incidence
 215       hits[8] = -1;                         // first padhit
 216       hits[9] = -1;                         // last pad hit
 217       hits[10] = mom[3];                    // hit Energy
 218       hits[11] = mom[0];                    // Px
 219       hits[12] = mom[1];                    // Py
 220       hits[13] = mom[2];                    // Pz
 221       hits[14] = gMC->TrackTime();          // time of flight
 222       new(lhits[fNhits++])
 223           AliMUONHit(fIshunt,gAlice->GetMCApp()->GetCurrentTrackNumber(),vol,hits);
 224
 225   }
 226 //  if( gMC->IsTrackExiting()) gMC->StopTrack();
 227 }
 228
 229
 230
 231
 232
 233
 234
 235
 236