MFT/AliMFTSegmentation.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 //
  18 //      Segmentation class for the planes of the ALICE Muon Forward Tracker
  19 //
  20 //      Contact author: antonio.uras@cern.ch
  21 //
  22 //====================================================================================================================================================
  23
  24 #include "TFile.h"
  25 #include "TNtuple.h"
  26 #include "TClonesArray.h"
  27 #include "TMath.h"
  28 #include "AliMFTPlane.h"
  29 #include "AliMFTSegmentation.h"
  30
  31 ClassImp(AliMFTSegmentation)
  32
  33 //====================================================================================================================================================
  34
  35 AliMFTSegmentation::AliMFTSegmentation():
  36   TObject(),
  37   fMFTPlanes(0)
  38 {
  39
  40   // default constructor
  41
  42 }
  43
  44 //====================================================================================================================================================
  45
  46 AliMFTSegmentation::AliMFTSegmentation(const Char_t *nameGeomFile):
  47   TObject(),
  48   fMFTPlanes(new TClonesArray("AliMFTPlane", fNMaxPlanes))
  49 {
  50
  51   // constructor
  52
  53   Float_t zCenter, rMin, rMax, pixelSizeX, pixelSizeY, thicknessActive, thicknessSupport, thicknessReadout;
  54   Float_t equivalentSilicon, equivalentSiliconBeforeFront, equivalentSiliconBeforeBack;
  55
  56   TFile *geomFile = new TFile(nameGeomFile);
  57   TNtuple *geomNtuple = (TNtuple*) geomFile->Get("AliMFTGeometry");
  58
  59   geomNtuple -> SetBranchAddress("zCenter", &zCenter);
  60   geomNtuple -> SetBranchAddress("rMin",    &rMin);
  61   geomNtuple -> SetBranchAddress("rMax",    &rMax);
  62   geomNtuple -> SetBranchAddress("pixelSizeX", &pixelSizeX);
  63   geomNtuple -> SetBranchAddress("pixelSizeY", &pixelSizeY);
  64   geomNtuple -> SetBranchAddress("thicknessActive",  &thicknessActive);
  65   geomNtuple -> SetBranchAddress("thicknessSupport", &thicknessSupport);
  66   geomNtuple -> SetBranchAddress("thicknessReadout", &thicknessReadout);
  67   geomNtuple -> SetBranchAddress("equivalentSilicon",            &equivalentSilicon);
  68   geomNtuple -> SetBranchAddress("equivalentSiliconBeforeFront", &equivalentSiliconBeforeFront);
  69   geomNtuple -> SetBranchAddress("equivalentSiliconBeforeBack",  &equivalentSiliconBeforeBack);
  70
  71   Int_t nPlanes = geomNtuple->GetEntries();
  72
  73   for (Int_t iPlane=0; iPlane<nPlanes; iPlane++) {
  74
  75     // Create new plane
  76
  77     printf("Setting segmentation for MFT plane #%02d\n", iPlane);
  78
  79     geomNtuple -> GetEntry(iPlane);
  80     zCenter = TMath::Abs(zCenter);
  81
  82     AliMFTPlane *plane = new AliMFTPlane(Form("MFTPlane_%02d", iPlane), Form("MFTPlane_%02d", iPlane));
  83     plane -> Init(iPlane, zCenter, rMin, rMax, pixelSizeX, pixelSizeY, thicknessActive, thicknessSupport, thicknessReadout);
  84     plane -> SetEquivalentSilicon(equivalentSilicon);
  85     plane -> SetEquivalentSiliconBeforeFront(equivalentSiliconBeforeFront);
  86     plane -> SetEquivalentSiliconBeforeBack(equivalentSiliconBeforeBack);
  87     plane -> CreateStructure();
  88
  89     new ((*fMFTPlanes)[fMFTPlanes->GetEntries()]) AliMFTPlane(*plane);
  90
  91   }
  92
  93   delete geomFile;
  94
  95   printf("MFT segmentation set!\n");
  96
  97 }
  98
  99 //====================================================================================================================================================
 100
 101 THnSparseC* AliMFTSegmentation::GetDetElem(Int_t detElemID) const {
 102
 103   // Find det elem
 104
 105   Int_t planeNb = detElemID/fNMaxDetElemPerPlane;
 106   Int_t detElemNb = detElemID - planeNb*fNMaxDetElemPerPlane;
 107
 108   THnSparseC *detElem = GetPlane(planeNb)->GetActiveElement(detElemNb);
 109
 110   return detElem;
 111
 112 }
 113
 114 //====================================================================================================================================================
 115
 116 Bool_t AliMFTSegmentation::Hit2PixelID(Double_t xHit, Double_t yHit, Int_t detElemID, Int_t &xPixel, Int_t &yPixel) {
 117
 118   THnSparseC *detElem = GetDetElem(detElemID);
 119
 120   if ( xHit<detElem->GetAxis(0)->GetXmin() ||
 121        xHit>detElem->GetAxis(0)->GetXmax() ||
 122        yHit<detElem->GetAxis(1)->GetXmin() ||
 123        yHit>detElem->GetAxis(1)->GetXmax() ) return kFALSE;
 124
 125   xPixel = detElem->GetAxis(0)->FindBin(xHit) - 1;
 126   yPixel = detElem->GetAxis(1)->FindBin(yHit) - 1;
 127
 128   return kTRUE;
 129
 130 }
 131
 132 //====================================================================================================================================================
 133