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 /* $Id$ */
  17
  18 ///////////////////////////////////////////////////////////////////////////////
  19 //                                                                           //
  20 //  Transition Radiation Detector version 0 -- fast simulator                //
  21 //                                                                           //
  22 //Begin_Html
  23 /*
  24 <img src="picts/AliTRDfullClass.gif">
  25 */
  26 //End_Html
  27 //                                                                           //
  28 //                                                                           //
  29 ///////////////////////////////////////////////////////////////////////////////
  30
  31 #include <stdlib.h>
  32
  33 #include <TLorentzVector.h>
  34 #include <TMath.h>
  35 #include <TRandom.h>
  36 #include <TVector.h>
  37 #include <TVirtualMC.h>
  38
  39 #include "AliConst.h"
  40 #include "AliRun.h"
  41 #include "AliTRDgeometry.h"
  42 #include "AliTRDhit.h"
  43 #include "AliTRDv0.h"
  44
  45 ClassImp(AliTRDv0)
  46
  47 //_____________________________________________________________________________
  48 AliTRDv0::AliTRDv0():AliTRD()
  49 {
  50   //
  51   // AliTRDv0 default constructor
  52   //
  53
  54   fHitsOn     = 0;
  55
  56 }
  57
  58 //_____________________________________________________________________________
  59 AliTRDv0::AliTRDv0(const char *name, const char *title)
  60          :AliTRD(name, title)
  61 {
  62   //
  63   // Standard constructor for Transition Radiation Detector version 0
  64   //
  65
  66   fHitsOn     = 0;
  67
  68 }
  69
  70 //_____________________________________________________________________________
  71 AliTRDv0::~AliTRDv0()
  72 {
  73   //
  74   // AliTRDv0 destructor
  75   //
  76
  77 }
  78
  79 //_____________________________________________________________________________
  80 void AliTRDv0::CreateGeometry()
  81 {
  82   //
  83   // Create the GEANT geometry for the Transition Radiation Detector - Version 0
  84   // This version covers the full azimuth.
  85   //
  86
  87   // Check that FRAME is there otherwise we have no place where to put the TRD
  88   AliModule* frame = gAlice->GetModule("FRAME");
  89   if (!frame) return;
  90
  91   // Define the chambers
  92   AliTRD::CreateGeometry();
  93
  94 }
  95
  96 //_____________________________________________________________________________
  97 void AliTRDv0::CreateMaterials()
  98 {
  99   //
 100   // Create materials for the Transition Radiation Detector
 101   //
 102
 103   AliTRD::CreateMaterials();
 104
 105 }
 106
 107 //_____________________________________________________________________________
 108 void AliTRDv0::Init()
 109 {
 110   //
 111   // Initialize Transition Radiation Detector after geometry is built
 112   //
 113
 114   AliTRD::Init();
 115
 116   printf("          Fast simulator\n\n");
 117   for (Int_t i = 0; i < 80; i++) printf("*");
 118   printf("\n");
 119
 120 }
 121
 122 //_____________________________________________________________________________
 123 void AliTRDv0::StepManager()
 124 {
 125   //
 126   // Procedure called at every step in the TRD
 127   // Fast simulator. If switched on, a hit is produced when a track
 128   // crosses the border between amplification region and pad plane.
 129   //
 130
 131   Int_t   pla = 0;
 132   Int_t   cha = 0;
 133   Int_t   sec = 0;
 134
 135   Float_t hits[3];
 136   Int_t   det;
 137
 138   TLorentzVector p;
 139
 140   // Use pad plane as sensitive volume
 141   TString  cIdSens = "L";
 142   TString  cIdCurrent;
 143   Char_t   cIdChamber[3];
 144            cIdChamber[2] = 0;
 145
 146   const Int_t kNplan = AliTRDgeometry::Nplan();
 147
 148   // Writing out hits enabled?
 149   if (!(fHitsOn)) return;
 150
 151   // Use only charged tracks and count them only once per volume
 152   if (gMC->TrackCharge()    &&
 153       gMC->IsTrackEntering()) {
 154
 155     // Check on sensitive volume
 156     cIdCurrent = gMC->CurrentVolName();
 157     if (cIdSens == cIdCurrent[1]) {
 158
 159       gMC->TrackPosition(p);
 160       for (Int_t i = 0; i < 3; i++) hits[i] = p[i];
 161
 162       // The sector number (0 - 17)
 163       // The numbering goes clockwise and starts at y = 0
 164       Float_t phi = kRaddeg*TMath::ATan2(hits[0],hits[1]);
 165       if (phi < 90.)
 166         phi = phi + 270.;
 167       else
 168         phi = phi -  90.;
 169       sec = ((Int_t) (phi / 20));
 170
 171       // The plane and chamber number
 172       cIdChamber[0] = cIdCurrent[2];
 173       cIdChamber[1] = cIdCurrent[3];
 174       Int_t idChamber = atoi(cIdChamber);
 175       cha = ((Int_t) idChamber / kNplan);
 176       pla = ((Int_t) idChamber % kNplan);
 177       det = fGeometry->GetDetector(pla,cha,sec);
 178
 179       AddHit(gAlice->CurrentTrack(),det,hits,0,kTRUE);
 180
 181     }
 182
 183   }
 184
 185 }