MUON/AliMUONManuContourMaker.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 AliMUONManuContourMaker
  19 ///
  20 /// Maker of manu contours.
  21 ///
  22 /// Make use of the AliMUONContourMaker class, but this one contains
  23 /// specific things for MUON (as the mapping, for instance), hence its
  24 /// separation from AliMUONContourMaker.
  25 ///
  26 /// This class requires that the mapping is loaded before anything can be done.
  27 ///
  28 /// \author Laurent Aphecetche, Subatech
  29
  30 #include "AliMUONManuContourMaker.h"
  31
  32 #include "AliCodeTimer.h"
  33 #include "AliLog.h"
  34 #include "AliMUONContour.h"
  35 #include "AliMUONContourMaker.h"
  36 #include "AliMUONPolygon.h"
  37 #include "AliMpCathodType.h"
  38 #include "AliMpConnection.h"
  39 #include "AliMpConstants.h"
  40 #include "AliMpDDLStore.h"
  41 #include "AliMpDEManager.h"
  42 #include "AliMpIntPair.h"
  43 #include "AliMpMotifPosition.h"
  44 #include "AliMpMotifType.h"
  45 #include "AliMpManuIterator.h"
  46 #include "AliMpPlaneType.h"
  47 #include "AliMpSegmentation.h"
  48 #include "AliMpUID.h"
  49 #include "AliMpVMotif.h"
  50 #include "AliMpVSegmentation.h"
  51 #include "TGeoMatrix.h"
  52 #include "TObjArray.h"
  53 #include "TObjString.h"
  54 #include "TString.h"
  55 #include "TVector2.h"
  56
  57 ///\cond CLASSIMP
  58 ClassImp(AliMUONManuContourMaker)
  59 ///\endcond
  60
  61 //_____________________________________________________________________________
  62 AliMUONManuContourMaker::AliMUONManuContourMaker(AliMpExMap* deTransformations)
  63 : TObject(), fDETransformations(deTransformations), fLocalManuContours(222,1)
  64 {
  65   fLocalManuContours.SetOwnerKeyValue(kTRUE,kTRUE);
  66 }
  67
  68 //_____________________________________________________________________________
  69 AliMUONManuContourMaker::~AliMUONManuContourMaker()
  70 {
  71 }
  72
  73 //_____________________________________________________________________________
  74 AliMUONContour*
  75 AliMUONManuContourMaker::CreateManuContour(Int_t detElemId, Int_t manuId, const char* name) const
  76 {
  77   /// Create the contour of a given manu (global coordinates)
  78
  79   AliCodeTimerAuto("");
  80
  81   TString sname(name);
  82
  83   if ( sname.Length()==0 )
  84   {
  85     sname = ManuPathName(detElemId,manuId);
  86   }
  87
  88   const AliMpVSegmentation* seg = AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
  89   const AliMpMotifPosition* motifPos = seg->MotifPosition(manuId);
  90
  91   AliMUONContour* contour = CreateMotifContour(*motifPos);
  92
  93   if (!contour)
  94   {
  95     AliError(Form("Could not build contour %s",sname.Data()));
  96     return 0x0;
  97   }
  98
  99   contour->SetName(sname.Data());
 100
 101   contour->Offset(motifPos->GetPositionX()-seg->GetPositionX(),
 102                   motifPos->GetPositionY()-seg->GetPositionY());
 103
 104   TGeoHMatrix* matrix = 0x0;
 105
 106   if ( fDETransformations )
 107   {
 108     matrix = static_cast<TGeoHMatrix*>(fDETransformations->GetValue(detElemId));
 109     if ( matrix ) contour->Transform(*matrix);
 110   }
 111
 112   return contour;
 113 }
 114
 115
 116 //_____________________________________________________________________________
 117 AliMUONContour*
 118 AliMUONManuContourMaker::CreateMotifContour(const AliMpMotifPosition& motifPosition) const
 119 {
 120   /// Create the contour of a given MOTIF (i.e. local coordinates only).
 121
 122   AliCodeTimerAuto("");
 123
 124   TString mpName(NameIt(motifPosition));
 125
 126   AliMUONContour* contour = static_cast<AliMUONContour*>(fLocalManuContours.GetValue(mpName.Data()));
 127
 128   if ( contour )
 129   {
 130     // if we have already done the job, just have to clone it and we are done
 131     return static_cast<AliMUONContour*>(contour->Clone());
 132   }
 133
 134   TObjArray polygons(AliMpConstants::ManuNofChannels()); // array of AliMUONPolygon objects
 135   polygons.SetOwner(kTRUE);
 136
 137   AliMpVMotif* motif = motifPosition.GetMotif();
 138
 139   AliMpMotifType* motifType = motif->GetMotifType();
 140
 141   if ( motifType->IsFull() )
 142   {
 143     // motif is a simple rectangle. No need to loop over pads, we can
 144     // compute the contour right here and now.
 145     polygons.Add(new AliMUONPolygon(0.0,0.0,motif->DimensionX(),motif->DimensionY()));
 146   }
 147   else
 148   {
 149     for ( Int_t i = 0; i <= AliMpConstants::ManuNofChannels(); ++i )
 150     {
 151       AliMpConnection* connection = motifType->FindConnectionByGassiNum(i);
 152
 153       if ( connection )
 154       {
 155         Int_t ix = connection->GetLocalIx();
 156         Int_t iy = connection->GetLocalIy();
 157
 158         Double_t x,y,dx,dy;
 159
 160         motif->GetPadDimensionsByIndices(ix,iy,dx,dy);
 161         motif->PadPositionLocal(ix,iy,x,y);
 162
 163         AliMUONPolygon* pol = new AliMUONPolygon(x,y,dx,dy);
 164         polygons.Add(pol);
 165       }
 166     }
 167   }
 168
 169   AliMUONContourMaker maker;
 170
 171   contour = maker.CreateContour(polygons);
 172
 173   if (!contour || !contour->IsValid() )
 174   {
 175     AliError(Form("Failed to properly create contour %s contour = %p",mpName.Data(),contour));
 176     if ( contour )
 177     {
 178       AliError(Form("nofVertices=%d area.isvalid=%d",contour->NumberOfVertices(),contour->Area().IsValid()));
 179       StdoutToAliError(contour->Area().Print(););
 180     }
 181     delete contour;
 182     return 0x0;
 183   }
 184
 185   {
 186     AliCodeTimerAuto("localmanucontour.add");
 187     fLocalManuContours.Add(new TObjString(mpName),contour);
 188   }
 189
 190   return static_cast<AliMUONContour*>(contour->Clone());
 191 }
 192
 193 //_____________________________________________________________________________
 194 TObjArray*
 195 AliMUONManuContourMaker::GenerateManuContours(Bool_t stopAtError)
 196 {
 197   /// Generate the contours for all the manus, taking into account the given transformation
 198   /// (to go from local to global). That transformation need not be the real one (i.e.
 199   /// it can be an "exploded" one to ease visualization).
 200
 201   AliCodeTimerAuto("");
 202
 203   TObjArray* manuContours = new TObjArray;
 204
 205   manuContours->SetOwner(kTRUE);
 206
 207   AliMpManuIterator it;
 208   Int_t detElemId, manuId;
 209   Int_t nmanus(0);
 210   Int_t nok(0);
 211
 212   while ( it.Next(detElemId,manuId) )
 213   {
 214     ++nmanus;
 215     AliMUONContour* contour = CreateManuContour(detElemId,manuId);
 216     if (contour)
 217     {
 218       manuContours->Add(contour);
 219     }
 220     else
 221     {
 222       if ( stopAtError )
 223       {
 224         break;
 225       }
 226     }
 227     ++nok;
 228   }
 229
 230   AliInfo(Form("%d manus. %d contours successfully created",nmanus,nok));
 231
 232   return manuContours;
 233 }
 234
 235 //_____________________________________________________________________________
 236 TString
 237 AliMUONManuContourMaker::NameIt(const AliMpMotifPosition& motifPosition) const
 238 {
 239   /// Get the name of an AliMpMotifPosition
 240
 241   AliMpVMotif* motif = motifPosition.GetMotif();
 242   TString name(Form("%s",motif->GetID().Data()));
 243
 244   for ( Int_t i = 0; i < motif->GetNofPadDimensions(); ++i )
 245   {
 246     name += Form("/%7.3f-%7.3f:",motif->GetPadDimensionX(i),motif->GetPadDimensionY(i));
 247   }
 248
 249   return name;
 250 }
 251
 252 //_____________________________________________________________________________
 253 TString
 254 AliMUONManuContourMaker::ManuPathName(Int_t detElemId, Int_t manuId, Bool_t withCathodeName)
 255 {
 256   /// Get the name of a manu
 257
 258   AliMp::PlaneType planeType;
 259   if ( manuId & AliMpConstants::ManuMask(AliMp::kNonBendingPlane) )
 260   {
 261     planeType = AliMp::kNonBendingPlane;
 262   }
 263   else
 264   {
 265     planeType = AliMp::kBendingPlane;
 266   }
 267   AliMp::CathodType cathodeType = AliMpDEManager::GetCathod(detElemId,planeType);
 268
 269   Int_t chamberId = AliMpDEManager::GetChamberId(detElemId);
 270   Int_t stationId = chamberId/2;
 271
 272   Int_t busPatchId = AliMpDDLStore::Instance()->GetBusPatchId(detElemId, manuId);
 273
 274   AliMpUID id(cathodeType,stationId,chamberId,detElemId,busPatchId,manuId);
 275
 276   if ( withCathodeName ) return id.PathName();
 277
 278   TString name(id.PathName());
 279
 280   name.ReplaceAll("Cathode0/","");
 281   name.ReplaceAll("Cathode1/","");
 282
 283   return name;
 284 }
 285
 286
 287
 288