TRD/AliTRDv0.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 */
  19
  20 ///////////////////////////////////////////////////////////////////////////////
  21 //                                                                           //
  22 //  Transition Radiation Detector version 0 -- coarse simulation             //
  23 //                                                                           //
  24 //Begin_Html
  25 /*
  26 <img src="picts/AliTRDv0Class.gif">
  27 */
  28 //End_Html
  29 //                                                                           //
  30 //                                                                           //
  31 ///////////////////////////////////////////////////////////////////////////////
  32
  33 #include <TMath.h>
  34 #include <TRandom.h>
  35 #include <TVector.h>
  36
  37 #include "AliTRDv0.h"
  38 #include "AliRun.h"
  39 #include "AliMC.h"
  40 #include "AliConst.h"
  41
  42 ClassImp(AliTRDv0)
  43
  44 //_____________________________________________________________________________
  45 AliTRDv0::AliTRDv0(const char *name, const char *title)
  46          :AliTRD(name, title)
  47 {
  48   //
  49   // Standard constructor for Transition Radiation Detector version 0
  50   //
  51
  52   fIdSens     = 0;
  53   fHitsOn     = 0;
  54
  55   fIdSpace1   = 0;
  56   fIdSpace2   = 0;
  57   fIdSpace3   = 0;
  58
  59   fIdChamber1 = 0;
  60   fIdChamber2 = 0;
  61   fIdChamber3 = 0;
  62
  63 }
  64
  65 //_____________________________________________________________________________
  66 void AliTRDv0::CreateGeometry()
  67 {
  68   //
  69   // Create the GEANT geometry for the Transition Radiation Detector - Version 0
  70   // This version covers the full azimuth.
  71   //
  72   // Author:  Christoph Blume (C.Blume@gsi.de) 20/07/99
  73   //
  74
  75   Float_t xpos, ypos, zpos;
  76
  77   // Check that FRAME is there otherwise we have no place where to put the TRD
  78   AliModule* FRAME = gAlice->GetModule("FRAME");
  79   if (!FRAME) return;
  80
  81   // Define the chambers
  82   AliTRD::CreateGeometry();
  83
  84   // Position the the TRD-sectors in all TRD-volumes in the spaceframe
  85   xpos     = 0.;
  86   ypos     = 0.;
  87   zpos     = 0.;
  88   gMC->Gspos("TRD ",1,"BTR1",xpos,ypos,zpos,0,"ONLY");
  89   gMC->Gspos("TRD ",2,"BTR2",xpos,ypos,zpos,0,"ONLY");
  90   gMC->Gspos("TRD ",3,"BTR3",xpos,ypos,zpos,0,"ONLY");
  91
  92 }
  93
  94 //_____________________________________________________________________________
  95 void AliTRDv0::CreateMaterials()
  96 {
  97   //
  98   // Create materials for the Transition Radiation Detector
  99   //
 100
 101   AliTRD::CreateMaterials();
 102
 103 }
 104
 105 //_____________________________________________________________________________
 106 void AliTRDv0::Init()
 107 {
 108   //
 109   // Initialise Transition Radiation Detector after geometry is built
 110   //
 111
 112   AliTRD::Init();
 113
 114   for (Int_t i = 0; i < 80; i++) printf("*");
 115   printf("\n");
 116
 117   // Identifier of the sensitive volume (amplification region)
 118   fIdSens     = gMC->VolId("UL06");
 119
 120   // Identifier of the TRD-spaceframe volumina
 121   fIdSpace1   = gMC->VolId("B028");
 122   fIdSpace2   = gMC->VolId("B029");
 123   fIdSpace3   = gMC->VolId("B030");
 124
 125   // Identifier of the TRD-driftchambers
 126   fIdChamber1 = gMC->VolId("UCIO");
 127   fIdChamber2 = gMC->VolId("UCIM");
 128   fIdChamber3 = gMC->VolId("UCII");
 129
 130 }
 131
 132 //_____________________________________________________________________________
 133 void AliTRDv0::StepManager()
 134 {
 135   //
 136   // Procedure called at every step in the TRD
 137   // Fast simulator. If switched on, a hit is produced when a track
 138   // crosses the border between amplification region and pad plane.
 139   //
 140
 141   Int_t   vol[3];
 142   Int_t   iIdSens, icSens;
 143   Int_t   iIdSpace, icSpace;
 144   Int_t   iIdChamber, icChamber;
 145
 146   Int_t   secMap1[10] = {  3,  7,  8,  9, 10, 11,  2,  1, 18, 17 };
 147   Int_t   secMap2[ 5] = { 16, 15, 14, 13, 12 };
 148   Int_t   secMap3[ 3] = {  5,  6,  4 };
 149
 150   Float_t hits[4];
 151
 152   TLorentzVector p;
 153   TClonesArray  &lhits = *fHits;
 154
 155   // Writing out hits enabled?
 156   if (!(fHitsOn)) return;
 157
 158   // Use only charged tracks and count them only once per volume
 159   if (gMC->TrackCharge()    &&
 160       gMC->IsTrackExiting()) {
 161
 162     // Check on sensitive volume
 163     iIdSens = gMC->CurrentVolID(icSens);
 164     if (iIdSens == fIdSens) {
 165
 166       gMC->TrackPosition(p);
 167       for (Int_t i = 0; i < 3; i++) hits[i] = p[i];
 168       // No charge created
 169       hits[3] = 0;
 170
 171       iIdSpace   = gMC->CurrentVolOffID(4,icSpace  );
 172       iIdChamber = gMC->CurrentVolOffID(1,icChamber);
 173
 174       // The sector number
 175       if      (iIdSpace == fIdSpace1)
 176         vol[0] = secMap1[icSpace-1];
 177       else if (iIdSpace == fIdSpace2)
 178         vol[0] = secMap2[icSpace-1];
 179       else if (iIdSpace == fIdSpace3)
 180         vol[0] = secMap3[icSpace-1];
 181
 182       // The chamber number
 183       //   1: outer left
 184       //   2: middle left
 185       //   3: inner
 186       //   4: middle right
 187       //   5: outer right
 188       if      (iIdChamber == fIdChamber1)
 189         vol[1] = (hits[2] < 0 ? 1 : 5);
 190       else if (iIdChamber == fIdChamber2)
 191         vol[1] = (hits[2] < 0 ? 2 : 4);
 192       else if (iIdChamber == fIdChamber3)
 193         vol[1] = 3;
 194
 195       // The plane number
 196       vol[2] = icChamber - TMath::Nint((Float_t) (icChamber / 7)) * 6;
 197
 198       new(lhits[fNhits++]) AliTRDhit(fIshunt,gAlice->CurrentTrack(),vol,hits);
 199
 200     }
 201
 202   }
 203
 204 }