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