Separating run-dependent mapping data from data, which are not
[u/mrichter/AliRoot.git] / MUON / AliMUONGeometryMisAligner.h
1 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
2  * See cxx source for full Copyright notice                               */
3
4 // $Id$
5
6 /// \ingroup geometry
7 /// \class AliMUONGeometryMisAligner
8 /// \brief Class for misalignment of geometry transformations
9 //
10 //  Authors: Bruce Becker, Javier Castillo
11
12
13 #ifndef ALI_MUON_GEOMETRY_MIS_ALIGNER_H
14 #define ALI_MUON_GEOMETRY_MIS_ALIGNER_H
15
16 #include <TObject.h>
17
18 class AliMUONGeometryTransformer;
19
20 class TRandom;
21 class TGeoCombiTrans;
22
23 class AliMUONGeometryMisAligner:public TObject
24 {
25  public:
26   AliMUONGeometryMisAligner(Double_t cartXMisAligM, Double_t cartXMisAligW, Double_t cartYMisAligM, Double_t cartYMisAligW, Double_t angMisAligM, Double_t angMisAligW);
27   AliMUONGeometryMisAligner(Double_t cartMisAligM, Double_t cartMisAligW, Double_t angMisAligM, Double_t angMisAligW);
28   AliMUONGeometryMisAligner(Double_t cartMisAligW, Double_t angMisAligW);
29   AliMUONGeometryMisAligner();
30   virtual ~AliMUONGeometryMisAligner();
31   
32   //_________________________________________________________________
33   // methods
34   
35   // return a misaligned geometry obtained from the existing one.
36   AliMUONGeometryTransformer* MisAlign(const AliMUONGeometryTransformer* transformer, 
37                                        Bool_t verbose = kFALSE);
38   
39   /// Set cartesian displacement parameters different along x, y
40   void SetCartMisAlig(Double_t xmean, Double_t xwidth, Double_t ymean, Double_t ywidth, Double_t zmean = 0., Double_t zwidth = 0.)
41     {fDetElemMisAlig[0][0] = xmean; fDetElemMisAlig[0][1] = xwidth; fDetElemMisAlig[1][0] = ymean; fDetElemMisAlig[1][1] = ywidth; fDetElemMisAlig[2][0] = zmean; fDetElemMisAlig[2][1] = zwidth; }
42
43   /// Set cartesian displacement parameters, the same along x, y
44   void SetCartMisAlig(Double_t mean, Double_t width)
45     {fDetElemMisAlig[0][0] = mean; fDetElemMisAlig[0][1] = width; fDetElemMisAlig[1][0] = mean; fDetElemMisAlig[1][1] = width;}
46   
47   /// Set angular displacement
48   void SetAngMisAlig(Double_t zmean, Double_t zwidth, Double_t xmean = 0., Double_t xwidth = 0., Double_t ymean = 0., Double_t ywidth = 0.)
49     {fDetElemMisAlig[3][0] = xmean; fDetElemMisAlig[3][1] = xwidth; fDetElemMisAlig[4][0] = ymean; fDetElemMisAlig[4][1] = ywidth; fDetElemMisAlig[5][0] = zmean; fDetElemMisAlig[5][1] = zwidth;}
50   
51   /// Set cartesian displacement (Kept for backward compatibility)
52   void SetMaxCartMisAlig(Double_t width) 
53     {fDetElemMisAlig[0][0] = 0.0; fDetElemMisAlig[0][1] = width; fDetElemMisAlig[1][0] = 0.0; fDetElemMisAlig[1][1] = width;}
54   
55   /// Set angular displacement (Kept for backward compatibility)
56   void SetMaxAngMisAlig(Double_t width) 
57     {fDetElemMisAlig[5][0] = 0.0; fDetElemMisAlig[5][1] = width;}
58
59   void SetXYAngMisAligFactor(Double_t factor);
60
61   void SetZCartMisAligFactor(Double_t factor);
62
63   /// Set option for gaussian distribution 
64   void SetUseGaus(Bool_t usegaus)
65     {fUseGaus=usegaus; fUseUni=!usegaus;}
66
67   /// Set option for uniform distribution 
68   void SetUseUni(Bool_t useuni)
69     {fUseGaus=!useuni; fUseUni=useuni;}
70
71   /// Set module (half chambers) cartesian displacement parameters
72   void SetModuleCartMisAlig(Double_t xmean, Double_t xwidth, Double_t ymean, Double_t ywidth, Double_t zmean, Double_t zwidth) 
73     {fModuleMisAlig[0][0] = xmean; fModuleMisAlig[0][1] = xwidth; fModuleMisAlig[1][0] = ymean; fModuleMisAlig[1][1] = ywidth; fModuleMisAlig[2][0] = zmean; fModuleMisAlig[2][1] = zwidth;}
74
75   /// Set module (half chambers) cartesian displacement parameters
76   void SetModuleAngMisAlig(Double_t xmean, Double_t xwidth, Double_t ymean, Double_t ywidth, Double_t zmean, Double_t zwidth) 
77     {fModuleMisAlig[3][0] = xmean; fModuleMisAlig[3][1] = xwidth; fModuleMisAlig[4][0] = ymean; fModuleMisAlig[4][1] = ywidth; fModuleMisAlig[5][0] = zmean; fModuleMisAlig[5][1] = zwidth;}
78   
79  protected:
80   /// Not implemented
81   AliMUONGeometryMisAligner(const AliMUONGeometryMisAligner & right);
82   /// Not implemented
83   AliMUONGeometryMisAligner & operator =(const AliMUONGeometryMisAligner &right);
84   
85   
86  private:
87   // return a misaligned transformation
88   TGeoCombiTrans MisAlignDetElem(const TGeoCombiTrans& transform) const;
89   TGeoCombiTrans MisAlignModule(const TGeoCombiTrans& transform) const;
90   void GetUniMisAlign(Double_t cartMisAlig[3], Double_t angMisAlig[3], const Double_t lParMisAlig[6][2]) const;
91   void GetGausMisAlign(Double_t cartMisAlig[3], Double_t angMisAlig[3], const Double_t lParMisAlig[6][2]) const;
92
93   Bool_t fUseUni;            ///< use uniform distribution for misaligmnets
94   Bool_t fUseGaus;           ///< use gaussian distribution for misaligmnets
95   Double_t fDetElemMisAlig[6][2]; ///< Mean and width of the displacements of the detection elements along x,y,z (translations) and about x,y,z (rotations)
96   Double_t fModuleMisAlig[6][2];  ///< Mean and width of the displacements of the modules along x,y,z (translations) and about x,y,z (rotations)  
97
98   Double_t fXYAngMisAligFactor;  ///< factor (<1) to apply to angular misalignment range since range of motion is restricted out of the xy plane
99   Double_t fZCartMisAligFactor; ///< factor (<1) to apply to cartetian misalignment range since range of motion is restricted in z direction
100   TRandom *fDisplacementGenerator;  ///< random number generator for the displacements
101   
102   ClassDef(AliMUONGeometryMisAligner,3) // Geometry parametrisation
103 };
104
105 #endif //ALI_MUON_GEOMETRY_MIS_ALIGNER_H
106
107
108
109