START/AliSTART.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 //  START (T-Zero) Detector                                            //
  21 //  This class contains the base procedures for the START     //
  22 //  detector                                                                 //
  23 //                                                                           //
  24 //Begin_Html
  25 /*
  26 <img src="gif/AliSTARTClass.gif">
  27 </pre>
  28 <br clear=left>
  29 <font size=+2 color=red>
  30 <p>The responsible person for this module is
  31 <a href="mailto:Alla.Maevskaia@cern.ch">Alla Maevskaia</a>.
  32 </font>
  33 <pre>
  34 */
  35 //End_Html
  36 //                                                                           //
  37 //                                                                           //
  38 ///////////////////////////////////////////////////////////////////////////////
  39
  40 #include <Riostream.h>
  41
  42 #include <TFile.h>
  43 #include <TGeometry.h>
  44 #include <TMath.h>
  45 #include <TNode.h>
  46 #include <TParticle.h>
  47 #include <TRandom.h>
  48 #include <TTUBE.h>
  49 #include <TVirtualMC.h>
  50
  51 #include "AliLoader.h"
  52 #include "AliRun.h"
  53 #include "AliSTART.h"
  54 #include "AliSTARTdigit.h"
  55 #include "AliSTARThit.h"
  56 #include "AliSTARThitPhoton.h"
  57 #include "AliSTARTvertex.h"
  58 #include "AliMC.h"
  59 #include "AliSTARTDigitizer.h"
  60
  61 ClassImp(AliSTART)
  62
  63 static  AliSTARTdigit *digits;
  64
  65 //_____________________________________________________________________________
  66 AliSTART::AliSTART()
  67 {
  68   //
  69   // Default constructor for class AliSTART
  70   //
  71   fIshunt   = 1;
  72   fHits     = 0;
  73   fDigits   = 0;
  74   fPhotons  = 0;
  75 }
  76
  77 //_____________________________________________________________________________
  78 AliSTART::AliSTART(const char *name, const char *title)
  79        : AliDetector(name,title)
  80 {
  81   //
  82   // Standard constructor for START Detector
  83   //
  84
  85
  86   //
  87   // Initialise Hit array
  88   fHits       = new TClonesArray("AliSTARThit",  405);
  89   gAlice->GetMCApp()->AddHitList(fHits);
  90
  91   fPhotons  = new TClonesArray("AliSTARThitPhoton", 10000);
  92   gAlice->GetMCApp()->AddHitList (fPhotons);
  93
  94   fIshunt     =  1;
  95   fIdSens   =  0;
  96   fNPhotons =  0;
  97   SetMarkerColor(kRed);
  98 }
  99
 100 //_____________________________________________________________________________
 101 AliSTART::~AliSTART() {
 102   if (fHits) {
 103     fHits->Delete();
 104     delete fHits;
 105   }
 106   if (fPhotons) {
 107     fPhotons->Delete();
 108     delete fPhotons;
 109   }
 110 }
 111
 112 //_____________________________________________________________________________
 113 void AliSTART::AddHit(Int_t track, Int_t *vol, Float_t *hits)
 114 {
 115   //
 116   // Add a START hit
 117   //
 118   TClonesArray &lhits = *fHits;
 119   new(lhits[fNhits++]) AliSTARThit(fIshunt,track,vol,hits);
 120 }
 121
 122 //_____________________________________________________________________________
 123 void AliSTART::AddHitPhoton(Int_t track, Int_t *vol, Float_t *hits)
 124 {
 125   //  Add a START hit of photons
 126
 127   TClonesArray &lhits = *fPhotons;
 128   new(lhits[fNPhotons++]) AliSTARThitPhoton(fIshunt,track,vol,hits);
 129 }
 130
 131 //_____________________________________________________________________________
 132
 133 void AliSTART::AddDigit(Int_t * /*tracks*/, Int_t * /*digits*/)
 134 {
 135
 136   //  Add a START digit to the list. Dummy function.
 137
 138 }
 139
 140 //_____________________________________________________________________________
 141 void AliSTART::BuildGeometry()
 142 {
 143   //
 144   // Build simple ROOT TNode geometry for event display
 145   //
 146   TNode *node, *top;
 147   const int kColorSTART  = 19;
 148
 149   top=gAlice->GetGeometry()->GetNode("alice");
 150
 151   // START define the different volumes
 152   new TRotMatrix("rotx999","rot999",  90,0,90,90,180,0);
 153
 154   new TTUBE("S_0ST1","START  volume 1","void",5.,10.7,5.3);
 155   top->cd();
 156   node = new TNode("0ST1","0ST01","S_0ST1",0,0,-69.7,"");
 157   node->SetLineColor(kColorSTART);
 158   fNodes->Add(node);
 159
 160   new TTUBE("S_0ST2","START volume 2","void",5.,10.7,5.3);
 161   top->cd();
 162   node = new TNode("0ST2","0ST2","S_0ST2",0,0,350,"rotx999");
 163   node->SetLineColor(kColorSTART);
 164   fNodes->Add(node);
 165 }
 166
 167 //_____________________________________________________________________________
 168 Int_t AliSTART::DistanceToPrimitive(Int_t /*px*/, Int_t /*py*/)
 169 {
 170   //
 171   // Calculate the distance from the mouse to the START on the screen
 172   // Dummy routine
 173   //
 174   return 9999;
 175 }
 176
 177 //-------------------------------------------------------------------------
 178 void AliSTART::Init()
 179 {
 180   //
 181   // Initialis the START after it has been built
 182   Int_t i;
 183   //
 184   if(fDebug) {
 185     printf("\n%s: ",ClassName());
 186     for(i=0;i<35;i++) printf("*");
 187     printf(" START_INIT ");
 188     for(i=0;i<35;i++) printf("*");
 189     printf("\n%s: ",ClassName());
 190     //
 191     // Here the START initialisation code (if any!)
 192     for(i=0;i<80;i++) printf("*");
 193     printf("\n");
 194   }
 195 }
 196
 197 //---------------------------------------------------------------------------
 198 void AliSTART::MakeBranch(Option_t* option)
 199 {
 200   //
 201   // Specific START branches
 202   //
 203   // Create Tree branches for the START.
 204   Int_t buffersize = 4000;
 205   char branchname[20];
 206   sprintf(branchname,"%s",GetName());
 207
 208
 209   const char *cD = strstr(option,"D");
 210   const char *cH = strstr(option,"H");
 211
 212   if (cH && fLoader->TreeH())
 213   {
 214      if (fPhotons == 0x0) fPhotons  = new TClonesArray("AliSTARThitPhoton", 10000);
 215      sprintf (branchname, "%shitPhoton", GetName());
 216      MakeBranchInTree (fLoader->TreeH(), branchname, &fPhotons, 50000, 0);
 217      if (fHits == 0x0) fHits  = new TClonesArray("AliSTARThit",  405);
 218   }
 219
 220   AliDetector::MakeBranch(option);
 221
 222   if (cD) {
 223     digits = new AliSTARTdigit();
 224     MakeBranchInTree(fLoader->TreeD(), branchname, "AliSTARTdigit", digits, buffersize, 1, 0);
 225   }
 226 }
 227
 228 //_____________________________________________________________________________
 229 void AliSTART::ResetHits()
 230 {
 231   AliDetector::ResetHits();
 232
 233   fNPhotons = 0;
 234   if (fPhotons)  fPhotons->Clear();
 235 }
 236
 237 //_____________________________________________________________________________
 238 void AliSTART::SetTreeAddress()
 239 {
 240   TBranch  *branch;
 241   TTree    *treeH;
 242
 243
 244   treeH = TreeH();
 245
 246   if (treeH)
 247     {
 248       if (fPhotons == 0x0) fPhotons  = new TClonesArray("AliSTARThitPhoton", 10000);
 249       branch = treeH->GetBranch("STARThitPhoton");
 250       if (branch)  branch->SetAddress(&fPhotons);
 251       if (fHits == 0x0) fHits  = new TClonesArray("AliSTARThit",  405);
 252     }
 253
 254   AliDetector::SetTreeAddress();
 255
 256 }
 257
 258
 259 //_____________________________________________________________________________
 260
 261 void AliSTART::Hit2digit(Int_t /*evnum*/)
 262 {
 263 }
 264
 265 //_____________________________________________________________________________
 266 AliDigitizer* AliSTART::CreateDigitizer(AliRunDigitizer* manager)
 267 {
 268   return new AliSTARTDigitizer(manager);
 269 }