MUON/mapping/AliMpMotifSpecial.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 // $MpId: AliMpMotifSpecial.cxx,v 1.12 2006/05/24 13:58:41 ivana Exp $
  18 // Category: motif
  19 //
  20 // Class AliMpMotifSpecial
  21 // -----------------------
  22 // Class that defines a motif with its unique ID
  23 // and the motif type.
  24 // Included in AliRoot: 2003/05/02
  25 // Authors: David Guez, Ivana Hrivnacova; IPN Orsay
  26
  27 #include "AliMpMotifSpecial.h"
  28 #include "AliMpMotifType.h"
  29 #include "AliMpIntPair.h"
  30 #include "AliMpConstants.h"
  31
  32 #include "AliLog.h"
  33
  34 #include <TString.h>
  35
  36 /// \cond CLASSIMP
  37 ClassImp(AliMpMotifSpecial)
  38 /// \endcond
  39
  40
  41 //______________________________________________________________________________
  42 AliMpMotifSpecial::AliMpMotifSpecial():
  43   AliMpVMotif(),
  44   fDimensions(),
  45   fPadDimensionsVector(),
  46   fPadDimensionsVector2()
  47 {
  48   /// Default constructor
  49 }
  50
  51
  52 //______________________________________________________________________________
  53 AliMpMotifSpecial::AliMpMotifSpecial(const TString &id,
  54                                      AliMpMotifType *motifType)
  55   : AliMpVMotif(id,motifType),
  56     fDimensions(),
  57 #ifdef WITH_STL
  58     fPadDimensionsVector(),
  59 #endif
  60 #ifdef WITH_ROOT
  61     fPadDimensionsVector(true),
  62 #endif
  63     fPadDimensionsVector2()
  64
  65 {
  66   /// Standard constructor.
  67
  68 #ifdef WITH_STL
  69   fPadDimensionsVector.resize(motifType->GetNofPadsX()*motifType->GetNofPadsY());
  70 #endif
  71 }
  72
  73 //______________________________________________________________________________
  74 AliMpMotifSpecial::~AliMpMotifSpecial()
  75 {
  76   /// Destructor
  77 }
  78
  79
  80 //
  81 // private methods
  82 //
  83
  84 //______________________________________________________________________________
  85 Int_t AliMpMotifSpecial::VectorIndex(const AliMpIntPair& indices) const
  86 {
  87 /// Transform indices to linear vector index
  88
  89   return indices.GetFirst()*GetMotifType()->GetNofPadsY() + indices.GetSecond();
  90 }
  91
  92
  93 //
  94 // public methods
  95 //
  96
  97 #include <Riostream.h>
  98 //______________________________________________________________________________
  99 TVector2
 100 AliMpMotifSpecial::GetPadDimensions(const AliMpIntPair& localIndices) const
 101 {
 102 /// Return the dimensions of pad located at the given indices
 103
 104   if (GetMotifType()->HasPad(localIndices)) {
 105 #ifdef WITH_STL
 106     return fPadDimensionsVector[VectorIndex(localIndices)];
 107 #endif
 108 #ifdef WITH_ROOT
 109     if (!fPadDimensionsVector.GetValue(localIndices)) {
 110       Warning("GetPadDimensions","Indices outside limits");
 111       return TVector2(0.,0.);
 112     }
 113     else
 114       return  *((TVector2*)fPadDimensionsVector.GetValue(localIndices));
 115 #endif
 116   }
 117   else {
 118     Warning("GetPadDimensions","Indices outside limits");
 119     return TVector2(0.,0.);
 120   }
 121 }
 122
 123 //______________________________________________________________________________
 124 Int_t AliMpMotifSpecial::GetNofPadDimensions() const
 125 {
 126 /// Return number of different pad dimensions in this motif
 127
 128 #ifdef WITH_STL
 129   return fPadDimensionsVector2.size();
 130 #endif
 131
 132 #ifdef WITH_ROOT
 133   return fPadDimensionsVector2.GetEntriesFast();
 134 #endif
 135 }
 136
 137 //______________________________________________________________________________
 138 TVector2 AliMpMotifSpecial::GetPadDimensions(Int_t i) const
 139 {
 140 /// Returns the i-th different pad dimensions
 141
 142   if (i<0 || i>GetNofPadDimensions()) {
 143     AliFatal("Index outside limits.");
 144     return TVector2();
 145   }
 146
 147 #ifdef WITH_STL
 148   return fPadDimensionsVector2[i];
 149 #endif
 150
 151 #ifdef WITH_ROOT
 152   return *((TVector2*) fPadDimensionsVector2[i]);
 153 #endif
 154 }
 155
 156 //______________________________________________________________________________
 157 void AliMpMotifSpecial::CalculateDimensions()
 158 {
 159   /// Calculate motif dimensions and keep them in fDimensions data
 160
 161   Int_t i,j;
 162   Double_t sizeY=0.;
 163   Double_t sizeX=0.;
 164
 165   Double_t* tabSizeX = new Double_t[GetMotifType()->GetNofPadsY()];
 166   for (j=0;j<GetMotifType()->GetNofPadsY();++j) tabSizeX[j]=0.0;
 167
 168   for (i=0;i<GetMotifType()->GetNofPadsX();++i) {
 169     Double_t trSizeY=0.;
 170     for (j=0;j<GetMotifType()->GetNofPadsY();++j) {
 171       TVector2 dim = GetPadDimensions(AliMpIntPair(i,j));
 172       trSizeY+=dim.Y();
 173       tabSizeX[j]+=dim.X();
 174     }
 175     if (trSizeY>sizeY) sizeY=trSizeY;
 176   }
 177   for (j=0;j<GetMotifType()->GetNofPadsY();++j) {
 178     if (tabSizeX[j]>sizeX) sizeX = tabSizeX[j];
 179   }
 180
 181   delete [] tabSizeX;
 182
 183   fDimensions = TVector2(sizeX,sizeY);
 184 }
 185
 186 //______________________________________________________________________________
 187 TVector2 AliMpMotifSpecial::Dimensions() const
 188 {
 189   /// Give the dimension of the motif
 190
 191   return fDimensions;
 192 }
 193
 194 //______________________________________________________________________________
 195 TVector2
 196 AliMpMotifSpecial::PadPositionLocal(const AliMpIntPair& localIndices) const
 197 {
 198   /// Give the local position of the pad number (ix,iy)
 199   /// (0,0 is the center of the motif)
 200
 201   TVector2 dim = GetPadDimensions(localIndices);
 202
 203   Double_t posX= dim.X();
 204   for (Int_t i=0;i<localIndices.GetFirst();++i) {
 205     posX+=2.*GetPadDimensions(AliMpIntPair(i,localIndices.GetSecond())).X();
 206   }
 207
 208   Double_t posY= dim.Y();
 209   for (Int_t j=0;j<localIndices.GetSecond();++j) {
 210     posY+=2.*GetPadDimensions(AliMpIntPair(localIndices.GetFirst(),j)).Y();
 211   }
 212
 213   return TVector2(posX,posY)-Dimensions();
 214
 215 }
 216 //______________________________________________________________________________
 217 AliMpIntPair AliMpMotifSpecial::PadIndicesLocal(const TVector2& localPos) const
 218 {
 219   /// Return the pad indices from a given local position
 220   /// or AliMpIntPair::Invalid() if this position doesn't correspond to any valid
 221   /// connection
 222   ///
 223   /// *SOLEIL* : This code suppose that
 224   /// - 1) all cells have the same size along the Y direction
 225   /// - 2) the column 0 is entierly filled
 226
 227
 228   // First : find the j index
 229   TVector2 pos = localPos + Dimensions();
 230   Int_t j=0;
 231   Double_t y=pos.Y();
 232
 233   while (j<GetMotifType()->GetNofPadsY()) {
 234     TVector2 padDim = GetPadDimensions(AliMpIntPair(0,j));
 235     y-=2.*padDim.Y();
 236     if (y<0.) break;
 237     j++;
 238   }
 239
 240   // Test if it's outside limits
 241   if (j==GetMotifType()->GetNofPadsY()){
 242     Warning("PadIndicesLocal","The position is outside the motif");
 243     return AliMpIntPair::Invalid();
 244   }
 245
 246
 247   // now find the i index, in the j_th row
 248   Int_t i=0;
 249   Double_t x=pos.X();
 250
 251   while (i<GetMotifType()->GetNofPadsX()) {
 252     TVector2 padDim = GetPadDimensions(AliMpIntPair(i,j));
 253     x-=2.*padDim.X();
 254     if (x<0.) break;
 255     i++;
 256   }
 257
 258
 259   // Test if it's outside limits
 260
 261   if (i==GetMotifType()->GetNofPadsX()){
 262     Warning("PadIndicesLocal","The position is outside the motif");
 263     return AliMpIntPair::Invalid();
 264   }
 265
 266   // then return the found (i,j)
 267   return AliMpIntPair(i,j);
 268 }
 269 //______________________________________________________________________________
 270 void AliMpMotifSpecial::SetPadDimensions(const AliMpIntPair& localIndices,
 271                                          const TVector2& dimensions)
 272 {
 273   /// Set the dimensions of the pad located at \a localIndices to the given
 274   /// \a dimensions
 275
 276   if ( !GetMotifType()->HasPad(localIndices)){
 277     Warning("SetPadDimensions","Pad indices outside limits");
 278     return;
 279   }
 280
 281   // fill the dimensions map vector
 282 #ifdef WITH_STL
 283   fPadDimensionsVector[VectorIndex(localIndices)]=dimensions;
 284
 285   // fill the vector of different pad dimensions
 286   // only if these dimensions are not yet present
 287   Bool_t isPresent = false;
 288   for (Int_t i=0; i<GetNofPadDimensions(); i++) {
 289     if (AliMpConstants::IsEqual(fPadDimensionsVector2[i], dimensions))
 290       isPresent = true;
 291   }
 292
 293   if (!isPresent) fPadDimensionsVector2.push_back(dimensions);
 294 #endif
 295
 296 #ifdef WITH_ROOT
 297   TVector2* dimensionsObj = new TVector2(dimensions);
 298   fPadDimensionsVector.Add(localIndices, dimensionsObj);
 299
 300   // fill the vector of different pad dimensions
 301   // only if these dimensions are not yet present
 302   Bool_t isPresent = false;
 303   for (Int_t i=0; i<GetNofPadDimensions(); i++) {
 304     if (AliMpConstants::IsEqual(*((TVector2*) fPadDimensionsVector2[i]), dimensions))
 305       isPresent = true;
 306   }
 307
 308   if (!isPresent) fPadDimensionsVector2.Add(dimensionsObj);
 309 #endif
 310
 311 }