MUON/AliMUONGeometryModule.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 // Class AliMUONGeometryModule
  19 // -----------------------------
  20 // Class for definition of the detector module parameters
  21 // (the transformations of detection elements, mapping between
  22 //  sensitive volumes and detection elements).
  23 //
  24 // Author: Ivana Hrivnacova, IPN Orsay
  25
  26 #include <TVirtualMC.h>
  27 #include <TGeoMatrix.h>
  28 #include <TObjArray.h>
  29 #include <TArrayI.h>
  30 #include <Riostream.h>
  31
  32 #include "AliLog.h"
  33
  34 #include "AliMUONGeometryModule.h"
  35 #include "AliMUONGeometryEnvelope.h"
  36 #include "AliMUONGeometryEnvelopeStore.h"
  37 #include "AliMUONGeometryDetElement.h"
  38 #include "AliMUONGeometryStore.h"
  39 #include "AliMUONGeometrySVMap.h"
  40 #include "AliMUONGeometryDEIndexing.h"
  41
  42 ClassImp(AliMUONGeometryModule)
  43
  44 //______________________________________________________________________________
  45 AliMUONGeometryModule::AliMUONGeometryModule(Int_t moduleId)
  46  : TObject(),
  47    fModuleId(moduleId),
  48    fMotherVolume("ALIC"),
  49    fNofSVs(0),
  50    fSVVolumeIds(0),
  51    fTransformation(0),
  52    fEnvelopes(0),
  53    fDEIndexing(0),
  54    fDetElements(0),
  55    fSVMap(0)
  56 {
  57 // Standard constructor
  58
  59   // Chamber transformation
  60   fTransformation = new TGeoCombiTrans("");
  61
  62   // Arrays of volumes Ids
  63   fSVVolumeIds = new TArrayI(20);
  64
  65   // Sensitive volumes map
  66   fSVMap = new AliMUONGeometrySVMap(100);
  67
  68   // Get indexing
  69   fDEIndexing = new AliMUONGeometryDEIndexing(fModuleId, 0);
  70
  71   // Det elements transformation stores
  72   fDetElements = new AliMUONGeometryStore(fDEIndexing);
  73
  74   // Envelope store
  75   fEnvelopes = new AliMUONGeometryEnvelopeStore(fDetElements);
  76 }
  77
  78
  79 //______________________________________________________________________________
  80 AliMUONGeometryModule::AliMUONGeometryModule()
  81  : TObject(),
  82    fModuleId(0),
  83    fMotherVolume(),
  84    fNofSVs(0),
  85    fSVVolumeIds(0),
  86    fTransformation(0),
  87    fEnvelopes(0),
  88    fDEIndexing(0),
  89    fDetElements(0),
  90    fSVMap(0)
  91 {
  92 // Default constructor
  93 }
  94
  95
  96 //______________________________________________________________________________
  97 AliMUONGeometryModule::AliMUONGeometryModule(const AliMUONGeometryModule& rhs)
  98   : TObject(rhs)
  99 {
 100   AliFatal("Copy constructor is not implemented.");
 101 }
 102
 103 //______________________________________________________________________________
 104 AliMUONGeometryModule::~AliMUONGeometryModule() {
 105 //
 106
 107   delete fTransformation;
 108   delete fSVVolumeIds;
 109   delete fEnvelopes;
 110   delete fDEIndexing;
 111   delete fDetElements;
 112   delete fSVMap;
 113 }
 114
 115 //______________________________________________________________________________
 116 AliMUONGeometryModule&
 117 AliMUONGeometryModule::operator = (const AliMUONGeometryModule& rhs)
 118 {
 119   // check assignement to self
 120   if (this == &rhs) return *this;
 121
 122   AliFatal("Assignment operator is not implemented.");
 123
 124   return *this;
 125 }
 126
 127 //
 128 // private methods
 129 //
 130
 131 //______________________________________________________________________________
 132 Int_t AliMUONGeometryModule::GetSVIndex(Int_t svVolId) const
 133 {
 134 // Returns the index of the volume specified by volId
 135 // if it is present in the list of sensitive volumes
 136 // (or -1 if not present).
 137
 138   for (Int_t i=0; i<fNofSVs; i++) {
 139       if (fSVVolumeIds->At(i) == svVolId) return i;
 140   }
 141   return -1;
 142 }
 143
 144 //
 145 // public methods
 146 //
 147
 148 //______________________________________________________________________________
 149 void  AliMUONGeometryModule::Global2Local(Int_t detElemId,
 150                                   Float_t xg, Float_t yg, Float_t zg,
 151                                   Float_t& xl, Float_t& yl, Float_t& zl) const
 152 {
 153 // Transforms point from the global reference frame (ALIC)
 154 // to the local reference frame of the detection element specified
 155 // by detElemId.
 156 // ---
 157
 158   // Get detection element
 159   AliMUONGeometryDetElement* detElement = GetDetElement(detElemId);
 160   if (!detElement) return;
 161
 162   // Transform point
 163   detElement->Global2Local(xg, yg, zg, xl, yl, zl);
 164 }
 165
 166 //______________________________________________________________________________
 167 void  AliMUONGeometryModule::Global2Local(Int_t detElemId,
 168                                   Double_t xg, Double_t yg, Double_t zg,
 169                                   Double_t& xl, Double_t& yl, Double_t& zl) const
 170 {
 171 // Transforms point from the global reference frame (ALIC)
 172 // to the local reference frame of the detection element specified
 173 // by detElemId.
 174 // ---
 175
 176    // Get detection element
 177    AliMUONGeometryDetElement* detElement = GetDetElement(detElemId);
 178    if (!detElement) return;
 179
 180    // Transform point
 181    detElement->Global2Local(xg, yg, zg, xl, yl, zl);
 182 }
 183
 184 //______________________________________________________________________________
 185 void  AliMUONGeometryModule::Local2Global(Int_t detElemId,
 186                  Float_t xl, Float_t yl, Float_t zl,
 187                  Float_t& xg, Float_t& yg, Float_t& zg) const
 188 {
 189 // Transforms point from the local reference frame of the detection element
 190 // specified by detElemId to the global reference frame (ALIC).
 191 // ---
 192
 193   // Get detection element
 194   AliMUONGeometryDetElement* detElement = GetDetElement(detElemId);
 195   if (!detElement) return;
 196
 197    // Transform point
 198   detElement->Local2Global(xl, yl, zl, xg, yg, zg);
 199 }
 200
 201 //______________________________________________________________________________
 202 void  AliMUONGeometryModule::Local2Global(Int_t detElemId,
 203                  Double_t xl, Double_t yl, Double_t zl,
 204                  Double_t& xg, Double_t& yg, Double_t& zg) const
 205 {
 206 // Transforms point from the local reference frame of the detection element
 207 // specified by detElemId to the global reference frame (ALIC).
 208 // ---
 209
 210    // Get detection element
 211    AliMUONGeometryDetElement* detElement = GetDetElement(detElemId);
 212    if (!detElement) return;
 213
 214    // Transform point
 215    detElement->Local2Global(xl, yl, zl, xg, yg, zg);
 216 }
 217
 218 //______________________________________________________________________________
 219 void  AliMUONGeometryModule::SetTranslation(const TGeoTranslation& translation)
 220 {
 221 // Sets the module position wrt world.
 222 // ---
 223
 224   fTransformation
 225     ->SetTranslation(const_cast<Double_t*>(translation.GetTranslation()));
 226 }
 227
 228 //______________________________________________________________________________
 229 void  AliMUONGeometryModule::SetRotation(const TGeoRotation& rotation)
 230 {
 231 // Sets the module rotation wrt ALIC.
 232 // ---
 233
 234   TGeoRotation* rot = new TGeoRotation();
 235   rot->SetMatrix(const_cast<Double_t*>(rotation.GetRotationMatrix()));
 236
 237   fTransformation->SetRotation(rot);
 238 }
 239
 240 //______________________________________________________________________________
 241 void  AliMUONGeometryModule::SetSensitiveVolume(Int_t svVolId)
 242 {
 243 // Adds the volume specified by volId to the list of sensitive
 244 // volumes
 245 // ---
 246
 247   // Resize TArrayI if needed
 248   if (fSVVolumeIds->GetSize() == fNofSVs) fSVVolumeIds->Set(2*fNofSVs);
 249
 250   fSVVolumeIds->AddAt(svVolId, fNofSVs++);
 251 }
 252
 253 //______________________________________________________________________________
 254 void  AliMUONGeometryModule::SetSensitiveVolume(const TString& volName)
 255 {
 256 // Adds the volume specified by volName to the list of sensitive
 257 // volumes
 258 // ---
 259
 260   SetSensitiveVolume(gMC->VolId(volName));
 261 }
 262
 263 //______________________________________________________________________________
 264 void  AliMUONGeometryModule::SetAlign(Bool_t align)
 265 {
 266 // Sets alignement option to enevelope store.
 267 // ---
 268
 269   fEnvelopes->SetAlign(align);
 270 }
 271
 272 //______________________________________________________________________________
 273 AliMUONGeometryDetElement*
 274 AliMUONGeometryModule::FindBySensitiveVolume(const TString& sensVolume) const
 275 {
 276 // Finds TGeoCombiTrans for the detector element Id specified by aligned volume
 277 // ---
 278
 279   Int_t detElemId = fSVMap->GetDetElemId(sensVolume);
 280
 281   if (!detElemId) return 0;
 282         // The specified sensitive volume is not in the map
 283
 284   return (AliMUONGeometryDetElement*)fDetElements->Get(detElemId);
 285 }
 286
 287 //______________________________________________________________________________
 288 Bool_t AliMUONGeometryModule::IsSensitiveVolume(Int_t volId) const
 289 {
 290 // Checks if the volume specified by volId is present in the list
 291 // of sensitive volumes.
 292 // ---
 293
 294   for (Int_t i=0; i<fNofSVs; i++) {
 295       if (fSVVolumeIds->At(i) == volId) return kTRUE;
 296   }
 297   return kFALSE;
 298 }
 299
 300 //______________________________________________________________________________
 301 Bool_t AliMUONGeometryModule::IsSensitiveVolume(const TString& volName) const
 302 {
 303 // Checks if the volume specified by volName  is present in the list
 304 // of sensitive volumes.
 305 // ---
 306
 307   return IsSensitiveVolume(gMC->VolId(volName));
 308 }
 309
 310 //______________________________________________________________________________
 311 AliMUONGeometryDetElement*
 312 AliMUONGeometryModule::GetDetElement(Int_t detElemId) const
 313 {
 314 // Returns thethe detection element specified by detElemId.
 315 // Gives error if detection element is not defined.
 316 // ---
 317
 318    // Get detection element
 319    AliMUONGeometryDetElement* detElement
 320      = (AliMUONGeometryDetElement*) fDetElements->Get(detElemId);
 321
 322    if (!detElement) {
 323      AliErrorStream()
 324        << "Detection element " << detElemId
 325        << " not found in module " << fModuleId << endl;
 326      return 0;
 327    }
 328
 329    return detElement;
 330 }
 331
 332 /*
 333 //______________________________________________________________________________
 334 Int_t  AliMUONGeometryModule::GetNofDetElements() const
 335 {
 336 // Returns the number of detection elements
 337 // ---
 338
 339   return fDEIndexing->GetNofDetElements();
 340 }
 341
 342 //______________________________________________________________________________
 343 Int_t  AliMUONGeometryModule::GetDetElemId(Int_t i) const
 344 {
 345 // Returns the i-th detection element id
 346 // ---
 347
 348   return fDEIndexing->GetDetElementId(i);
 349 }
 350 */