MUON/AliMUONSurveyChamber.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 /// \class AliMUONSurveyChamber
  18 /// Class for the survey processing of the ALICE DiMuon spectrometer
  19 ///
  20 /// This object provides the methods specific to the chambers (frames)
  21 ///
  22 /// \author Javier Castillo
  23 //-----------------------------------------------------------------------------
  24
  25
  26 #include <TClonesArray.h>
  27 #include <TObjString.h>
  28 #include <TString.h>
  29 #include <TH2.h>
  30
  31 #include "AliLog.h"
  32 #include "AliSurveyObj.h"
  33 #include "AliSurveyPoint.h"
  34
  35 #include "AliMUONSurveyChamber.h"
  36 #include "AliMUONSurveyDetElem.h"
  37
  38 /// \cond CLASSIMP
  39 ClassImp(AliMUONSurveyChamber)
  40 /// \endcond
  41
  42 AliMUONSurveyChamber::AliMUONSurveyChamber(Int_t lChamberId)
  43   : AliMUONSurveyObj() /// Constructor with chamber id
  44   , fChamberId(lChamberId)
  45   , fNDetElem(0)
  46   , fSurveyObj(0x0)
  47   , fSurveyDetElem(0x0)
  48 {
  49   /// Constructor with the chamber id
  50   fSurveyObj = new AliSurveyObj();
  51   fSurveyDetElem = new TClonesArray("AliMUONSurveyDetElem",4);
  52 }
  53
  54 AliMUONSurveyChamber::~AliMUONSurveyChamber() {
  55   /// Destructor
  56   if (fSurveyObj) delete fSurveyObj;
  57   if (fSurveyDetElem) fSurveyDetElem->Delete();
  58 }
  59
  60 Int_t AliMUONSurveyChamber::AddSurveyDetElem(Int_t lDetElemId) {
  61   /// Add a surveyed detection element to this chamber
  62   if (!fSurveyDetElem) {
  63     fSurveyDetElem = new TClonesArray("AliMUONSurveyDetElem",4);
  64     fNDetElem=0;
  65   }
  66   new((*fSurveyDetElem)[fNDetElem++]) AliMUONSurveyDetElem(lDetElemId,this);
  67
  68   return fNDetElem;
  69 }
  70
  71 AliMUONSurveyDetElem* AliMUONSurveyChamber::GetDetElem(Int_t lDetElemIndex) {
  72   /// Return AluMUONSurveyDetElem at lDetElemIndex
  73   if (lDetElemIndex<0||(lDetElemIndex>=fNDetElem)) {
  74     return 0x0;
  75   }
  76   else {
  77     return (AliMUONSurveyDetElem*)fSurveyDetElem->At(lDetElemIndex);
  78   }
  79 }
  80
  81 Int_t AliMUONSurveyChamber::AddStickerTargets(TObjArray *pArray, TString stBaseName, Int_t lTargetMax) {
  82   /// Add a maximum of lTargetMax sticker targets with stBaseName from the pArray of targets
  83   return AliMUONSurveyObj::AddStickerTargets(pArray, stBaseName, lTargetMax);
  84 }
  85
  86 Int_t AliMUONSurveyChamber::AddStickerTargets(TString stBaseName, Int_t lTargetMax) {
  87   /// Add a maximum of lTargetMax sticker targets with stBaseName from internal SurveyObj
  88   if (!fSurveyObj || !fSurveyObj->GetData()) {
  89     AliError("Survey data is missing!");
  90     return 0;
  91   }
  92   return AddStickerTargets(fSurveyObj->GetData(),stBaseName,lTargetMax);
  93 }
  94
  95 Int_t AliMUONSurveyChamber::AddGButtonTargets(TObjArray *pArray, TString stBaseName, Int_t lTargetMax) {
  96   /// Add a maximum of lTargetMax global targets with stBaseName from the pArray of targets
  97   return AliMUONSurveyObj::AddGButtonTargets(pArray, stBaseName, lTargetMax);
  98 }
  99
 100 Int_t AliMUONSurveyChamber::AddGButtonTargets(TString stBaseName, Int_t lTargetMax) {
 101   /// Add a maximum of lTargetMax global button targets with stBaseName from internal SurveyObj
 102   if (!fSurveyObj || !fSurveyObj->GetData()) {
 103     AliError("Survey data is missing!");
 104     return 0;
 105   }
 106   return AddGButtonTargets(fSurveyObj->GetData(),stBaseName,lTargetMax);
 107 }
 108
 109 void AliMUONSurveyChamber::SetLocalTransformation(TGeoCombiTrans *localTrf, Bool_t ownerLocalTrf) {
 110   AliMUONSurveyObj::SetLocalTransformation(localTrf,ownerLocalTrf);
 111   AliMUONSurveyObj::SetBaseTransformation(localTrf,ownerLocalTrf);
 112 }
 113
 114 void AliMUONSurveyChamber::PrintSurveyReport() {
 115   /// Print the survey report information and data
 116   printf("--> %d\n", fSurveyObj->GetEntries());
 117
 118   printf("Title: \"%s\"\n", fSurveyObj->GetReportTitle().Data());
 119   printf("Date: \"%s\"\n", fSurveyObj->GetReportDate().Data());
 120   printf("Detector: \"%s\"\n", fSurveyObj->GetDetector().Data());
 121   printf("URL: \"%s\"\n", fSurveyObj->GetURL().Data());
 122   printf("Number: \"%d\"\n", fSurveyObj->GetReportNumber());
 123   printf("Version: \"%d\"\n", fSurveyObj->GetReportVersion());
 124   printf("Observations: \"%s\"\n", fSurveyObj->GetObservations().Data());
 125   printf("Coordinate System: \"%s\"\n", fSurveyObj->GetCoordSys().Data());
 126   printf("Measurement Units: \"%s\"\n", fSurveyObj->GetUnits().Data());
 127   printf("Nr Columns: \"%d\"\n", fSurveyObj->GetNrColumns());
 128
 129   TObjArray *colNames = fSurveyObj->GetColumnNames();
 130   for (Int_t i = 0; i < colNames->GetEntries(); ++i)
 131     printf("  Column %d --> \"%s\"\n", i, ((TObjString *) colNames->At(i))->GetString().Data());
 132
 133   // Get Array of surveyed points
 134   printf("Points:\n");
 135   TObjArray *points = fSurveyObj->GetData();
 136
 137   for (Int_t i = 0; i < points->GetEntries(); ++i)
 138     printf("  Point %d --> \"%s\"  %s \n", i, ((AliSurveyPoint *) points->At(i))->GetPointName().Data(), points->At(i)->GetName());
 139
 140 }
 141
 142 void AliMUONSurveyChamber::FillCPSTHistograms(TString baseNameC, TH2 *hCPSTc, TString baseNameA, TH2 *hCPSTa) {
 143   /// Fill Chamber Plane Sticker Targest histograms for monitoring
 144   if(baseNameC.IsNull()||!hCPSTc){
 145     AliError("Need base name for points on side C and/or a histogram for them!");
 146     return;
 147   }
 148   AliMUONSurveyObj::FillSTHistograms(baseNameC,hCPSTc,baseNameA,hCPSTa);
 149 }
 150
 151 void AliMUONSurveyChamber::FillDESTHistograms(TString baseNameC, TH2 *hDESTc, TString baseNameA, TH2 *hDESTa) {
 152   /// Fill Detection Element Sticker Targest histograms for monitoring
 153   for (Int_t iDE=0; iDE<GetNDetElem(); iDE++){
 154     GetDetElem(iDE)->FillSTHistograms(baseNameC,hDESTc,baseNameA,hDESTa);
 155   }
 156 }
 157
 158 Double_t AliMUONSurveyChamber::GetMeanDetElemAlignResX() {
 159   /// Return the average uncertainty of the det. elem. translations along x parameter
 160   Double_t alignResX = 0.;
 161   for (int iDE=0; iDE<GetNDetElem(); iDE++){
 162     alignResX += GetDetElem(iDE)->GetAlignResX();
 163   }
 164   if (GetNDetElem()==0){
 165     AliError("This Chamber has 0 detection elements!");
 166     return 0.;
 167   }
 168   return alignResX/GetNDetElem();
 169 }
 170
 171 Double_t AliMUONSurveyChamber::GetMeanDetElemAlignResY() {
 172   /// Return the average uncertainty of the det. elem. translations along y parameter
 173   Double_t alignResY = 0.;
 174   for (int iDE=0; iDE<GetNDetElem(); iDE++){
 175     alignResY += GetDetElem(iDE)->GetAlignResY();
 176   }
 177   if (GetNDetElem()==0){
 178     AliError("This Chamber has 0 detection elements!");
 179     return 0.;
 180   }
 181   return alignResY/GetNDetElem();
 182 }