FMD/AliFMD.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 //  Forward Multiplicity Detector based on Silicon plates                    //
  18 //  This class contains the base procedures for the Forward Multiplicity     //
  19 //  detector                                                                 //
  20 //  Detector consists of 6 Si volumes covered pseudorapidity interval         //
  21 //  from 1.6 to 6.0.                                                         //
  22 //                                                                           //
  23 //Begin_Html
  24 /*
  25 <img src="gif/AliFMDClass.gif">
  26 </pre>
  27 <br clear=left>
  28 <font size=+2 color=red>
  29 <p>The responsible person for this module is
  30 <a href="mailto:Alla.Maevskaia@cern.ch">Alla Maevskaia</a>.
  31 </font>
  32 <pre>
  33 */
  34 //End_Html
  35 //                                                                           //
  36 //                                                                           //
  37 ///////////////////////////////////////////////////////////////////////////////
  38
  39 #define DEBUG
  40 #include <TMath.h>
  41 #include <TGeometry.h>
  42 #include <TTUBE.h>
  43 #include <TTree.h>
  44 #include <TNode.h>
  45
  46 #include <TClonesArray.h>
  47 #include <TLorentzVector.h>
  48 #include "AliFMDv1.h"
  49 #include "AliRun.h"
  50 #include "AliMC.h"
  51 #include <iostream.h>
  52 #include <fstream.h>
  53 #include "AliMagF.h"
  54 #include "AliFMDhit.h"
  55 #include "AliFMDdigit.h"
  56 #include <stdlib.h>
  57
  58
  59 ClassImp(AliFMD)
  60
  61   //_____________________________________________________________________________
  62   AliFMD::AliFMD(): AliDetector()
  63 {
  64   //
  65   // Default constructor for class AliFMD
  66   //
  67   fIshunt   = 0;
  68 }
  69
  70 //_____________________________________________________________________________
  71 AliFMD::AliFMD(const char *name, const char *title)
  72   : AliDetector(name,title)
  73 {
  74   //
  75   // Standard constructor for Forward Multiplicity Detector
  76   //
  77
  78   //
  79   // Initialise Hit array
  80   fHits     = new TClonesArray("AliFMDhit", 1000);
  81   // Digits for each Si disk
  82   fDigits   = new TClonesArray("AliFMDdigit", 1000);
  83   fSDigits  = new TClonesArray("AliFMDdigit", 1000) ;
  84   gAlice->AddHitList(fHits);
  85
  86   fIshunt     =  0;
  87   fIdSens1    =  0;
  88
  89   SetMarkerColor(kRed);
  90 }
  91
  92
  93 AliFMD::~AliFMD()
  94 {
  95   if (fHits) {
  96       fHits->Delete();
  97       delete fHits;
  98   }
  99   delete fDigits ;
 100   delete fSDigits ;
 101
 102 }
 103 //_____________________________________________________________________________
 104 void AliFMD::AddHit(Int_t track, Int_t *vol, Float_t *hits)
 105 {
 106   //
 107   // Add a hit to the list
 108   //
 109   TClonesArray &lhits = *fHits;
 110   new(lhits[fNhits++]) AliFMDhit(fIshunt,track,vol,hits);
 111 }
 112 //_____________________________________________________________________________
 113 void AliFMD::AddDigit( Int_t *digits)
 114 {
 115   // add a real digit - as coming from data
 116
 117   // printf("AddDigit\n");
 118
 119    TClonesArray &ldigits = *fDigits;
 120    // new((*fDigits)[fNdigits++]) AliFMDdigit(digits);
 121   new(ldigits[fNdigits++]) AliFMDdigit(digits);
 122
 123 }
 124 //_____________________________________________________________________________
 125 void AliFMD::BuildGeometry()
 126 {
 127   //
 128   // Build simple ROOT TNode geometry for event display
 129   //
 130   TNode *node, *top;
 131   const int kColorFMD  = 7;
 132   //
 133   top=gAlice->GetGeometry()->GetNode("alice");
 134
 135   // FMD define the different volumes
 136   new TRotMatrix("rot901","rot901",  90, 0, 90, 90, 180, 0);
 137
 138   new TTUBE("S_FMD0","FMD  volume 0","void",4.73,17.7,1.5);
 139   top->cd();
 140   node = new TNode("FMD0","FMD0","S_FMD0",0,0,64,"");
 141   node->SetLineColor(kColorFMD);
 142   fNodes->Add(node);
 143
 144   new TTUBE("S_FMD1","FMD  volume 1","void",23.4,36.,1.5);
 145   top->cd();
 146   node = new TNode("FMD1","FMD1","S_FMD1",0,0,85,"");
 147   node->SetLineColor(kColorFMD);
 148   fNodes->Add(node);
 149
 150   new TTUBE("S_FMD2","FMD  volume 2","void",4.73,17.7,1.5);
 151   top->cd();
 152   node = new TNode("FMD2","FMD2","S_FMD2",0,0,-64,"");
 153   node->SetLineColor(kColorFMD);
 154   fNodes->Add(node);
 155
 156   new TTUBE("S_FMD3","FMD  volume 3","void",23.4,36.,1.5);
 157   top->cd();
 158   node = new TNode("FMD3","FMD3","S_FMD3",0,0,-85,"");
 159   node->SetLineColor(kColorFMD);
 160   fNodes->Add(node);
 161
 162   new TTUBE("S_FMD4","FMD  volume 4","void",5,15,0.015);
 163   top->cd();
 164   node = new TNode("FMD4","FMD4","S_FMD4",0,0,-270,"");
 165   node->SetLineColor(kColorFMD);
 166   fNodes->Add(node);
 167
 168
 169   new TTUBE("S_FMD5","FMD  volume 5","void",5,14,0.015);
 170   top->cd();
 171   node = new TNode("FMD5","FMD5","S_FMD5",0,0,-630,"");
 172   node->SetLineColor(kColorFMD);
 173   fNodes->Add(node);
 174
 175 }
 176
 177 //_____________________________________________________________________________
 178 Int_t AliFMD::DistanceToPrimitive(Int_t px, Int_t py)
 179 {
 180   //
 181   // Calculate the distance from the mouse to the FMD on the screen
 182   // Dummy routine
 183   //
 184   return 9999;
 185 }
 186 //___________________________________________
 187 void AliFMD::ResetHits()
 188 {
 189   // Reset number of clusters and the cluster array for this detector
 190   AliDetector::ResetHits();
 191 }
 192
 193 //____________________________________________
 194 void AliFMD::ResetDigits()
 195 {
 196     //
 197     // Reset number of digits and the digits array for this detector
 198    AliDetector::ResetHits();
 199    //
 200 }
 201
 202 //-------------------------------------------------------------------------
 203 void AliFMD::Init()
 204 {
 205   //
 206   // Initialis the FMD after it has been built
 207   Int_t i;
 208   AliMC* pMC = AliMC::GetMC();
 209   //
 210   printf("\n");
 211   for(i=0;i<35;i++) printf("*");
 212   printf(" FMD_INIT ");
 213   for(i=0;i<35;i++) printf("*");
 214   printf("\n");
 215   //
 216   // Here the FMD initialisation code (if any!)
 217   for(i=0;i<80;i++) printf("*");
 218   printf("\n");
 219   //
 220   //
 221   fIdSens1=pMC->VolId("GRIN"); //Si sensetive volume
 222
 223 }
 224 //---------------------------------------------------------------------
 225 void AliFMD::MakeBranch(Option_t* option, char *file)
 226 {
 227   // Create Tree branches for the FMD.
 228   const Int_t kBufferSize = 4000;
 229   char branchname[10];
 230   sprintf(branchname,"%s",GetName());
 231
 232   AliDetector::MakeBranch(option,file);
 233
 234   const char *cD = strstr(option,"D");
 235
 236   if (cD) {
 237
 238      gAlice->MakeBranchInTree(gAlice->TreeD(),
 239                                    branchname, &fDigits, kBufferSize, file) ;
 240
 241     printf("Making Branch %s for digits\n",branchname);
 242     gAlice->TreeD()->Print();
 243   }
 244 }
 245
 246 //_____________________________________________________________________________
 247 void AliFMD::SetTreeAddress()
 248 {
 249   // Set branch address for the Hits and Digits Tree.
 250   char branchname[30];
 251   AliDetector::SetTreeAddress();
 252
 253   TBranch *branch;
 254   TTree *treeD = gAlice->TreeD();
 255
 256
 257   if (treeD) {
 258         if (fDigits) {
 259               branch = treeD->GetBranch(branchname);
 260               if (branch) branch->SetAddress(&fDigits);
 261         }
 262
 263   }
 264   if(fSDigits)
 265     fSDigits->Clear();
 266
 267   if (gAlice->TreeS()  && fSDigits ) {
 268     branch = gAlice->TreeS()->GetBranch("FMD");
 269     if (branch) branch->SetAddress(&fSDigits) ;
 270   }
 271
 272
 273 }
 274
 275 //---------------------------------------------------------------------
 276 /*
 277 void AliFMD::SDigits2Digits()
 278 {
 279      Int_t  ibg=0, bgr=0;
 280     Int_t nbgr_ev= 0;
 281
 282     if (ibg) {
 283          printf("nbgr_ev %d\n",nbgr_ev);
 284  //        Hit2Digits(nbgr_ev,"Add"," ","galice_bgr.root");
 285     }
 286        else {
 287    //      Hit2Digits(nbgr_ev,"rien","","");
 288     }
 289
 290 }
 291 */
 292 //---------------------------------------------------------------------
 293
 294
 295
 296 void AliFMD::Eta2Radius(Float_t eta, Float_t zDisk, Float_t *radius)
 297 {
 298   Float_t expEta=TMath::Exp(-eta);
 299   Float_t theta=TMath::ATan(expEta);
 300   theta=2.*theta;
 301   Float_t rad=zDisk*(TMath::Tan(theta));
 302   *radius=rad;
 303
 304   printf(" eta %f radius %f\n", eta, rad);
 305 }
 306
 307 //---------------------------------------------------------------------
 308
 309 void AliFMD::Hits2SDigits(){
 310 #ifdef DEBUG
 311   cout<<"ALiFMD::Hits2SDigits> start...\n";
 312 #endif
 313   char * fileSDigits = 0 ;
 314   char * fileHeader=0 ;
 315   AliFMDSDigitizer * sd = new AliFMDSDigitizer("mgalice.root","FMD.SDigit.root") ;
 316   sd->Exec("") ;
 317   sd->Print("");
 318
 319   delete sd ;
 320 }
 321
 322
 323