1 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
2 * See cxx source for full Copyright notice */
7 /// \class AliMUONGeometryMisAligner
8 /// \brief Class for misalignment of geometry transformations
10 // Authors: Bruce Becker, Javier Castillo
13 #ifndef ALI_MUON_GEOMETRY_MIS_ALIGNER_H
14 #define ALI_MUON_GEOMETRY_MIS_ALIGNER_H
18 class AliMUONGeometryTransformer;
24 class AliMUONGeometryMisAligner:public TObject
27 AliMUONGeometryMisAligner(Double_t cartXMisAligM, Double_t cartXMisAligW, Double_t cartYMisAligM, Double_t cartYMisAligW, Double_t angMisAligM, Double_t angMisAligW);
28 AliMUONGeometryMisAligner(Double_t cartMisAligM, Double_t cartMisAligW, Double_t angMisAligM, Double_t angMisAligW);
29 AliMUONGeometryMisAligner(Double_t cartMisAligW, Double_t angMisAligW);
30 AliMUONGeometryMisAligner();
31 virtual ~AliMUONGeometryMisAligner();
33 //_________________________________________________________________
36 // return a misaligned geometry obtained from the existing one.
37 AliMUONGeometryTransformer* MisAlign(const AliMUONGeometryTransformer* transformer,
38 Bool_t verbose = kFALSE);
40 /// Set cartesian displacement parameters different along x, y
41 void SetCartMisAlig(Double_t xmean, Double_t xwidth, Double_t ymean, Double_t ywidth, Double_t zmean = 0., Double_t zwidth = 0.)
42 {fDetElemMisAlig[0][0] = xmean; fDetElemMisAlig[0][1] = xwidth; fDetElemMisAlig[1][0] = ymean; fDetElemMisAlig[1][1] = ywidth; fDetElemMisAlig[2][0] = zmean; fDetElemMisAlig[2][1] = zwidth; }
44 /// Set cartesian displacement parameters, the same along x, y
45 void SetCartMisAlig(Double_t mean, Double_t width)
46 {fDetElemMisAlig[0][0] = mean; fDetElemMisAlig[0][1] = width; fDetElemMisAlig[1][0] = mean; fDetElemMisAlig[1][1] = width;}
48 /// Set angular displacement
49 void SetAngMisAlig(Double_t zmean, Double_t zwidth, Double_t xmean = 0., Double_t xwidth = 0., Double_t ymean = 0., Double_t ywidth = 0.)
50 {fDetElemMisAlig[3][0] = xmean; fDetElemMisAlig[3][1] = xwidth; fDetElemMisAlig[4][0] = ymean; fDetElemMisAlig[4][1] = ywidth; fDetElemMisAlig[5][0] = zmean; fDetElemMisAlig[5][1] = zwidth;}
52 /// Set cartesian displacement (Kept for backward compatibility)
53 void SetMaxCartMisAlig(Double_t width)
54 {fDetElemMisAlig[0][0] = 0.0; fDetElemMisAlig[0][1] = width; fDetElemMisAlig[1][0] = 0.0; fDetElemMisAlig[1][1] = width;}
56 /// Set angular displacement (Kept for backward compatibility)
57 void SetMaxAngMisAlig(Double_t width)
58 {fDetElemMisAlig[5][0] = 0.0; fDetElemMisAlig[5][1] = width;}
60 void SetXYAngMisAligFactor(Double_t factor);
62 void SetZCartMisAligFactor(Double_t factor);
64 /// Set option for gaussian distribution
65 void SetUseGaus(Bool_t usegaus)
66 {fUseGaus=usegaus; fUseUni=!usegaus;}
68 /// Set option for uniform distribution
69 void SetUseUni(Bool_t useuni)
70 {fUseGaus=!useuni; fUseUni=useuni;}
72 /// Set module (half chambers) cartesian displacement parameters
73 void SetModuleCartMisAlig(Double_t xmean, Double_t xwidth, Double_t ymean, Double_t ywidth, Double_t zmean, Double_t zwidth)
74 {fModuleMisAlig[0][0] = xmean; fModuleMisAlig[0][1] = xwidth; fModuleMisAlig[1][0] = ymean; fModuleMisAlig[1][1] = ywidth; fModuleMisAlig[2][0] = zmean; fModuleMisAlig[2][1] = zwidth;}
76 /// Set module (half chambers) cartesian displacement parameters
77 void SetModuleAngMisAlig(Double_t xmean, Double_t xwidth, Double_t ymean, Double_t ywidth, Double_t zmean, Double_t zwidth)
78 {fModuleMisAlig[3][0] = xmean; fModuleMisAlig[3][1] = xwidth; fModuleMisAlig[4][0] = ymean; fModuleMisAlig[4][1] = ywidth; fModuleMisAlig[5][0] = zmean; fModuleMisAlig[5][1] = zwidth;}
80 /// Set alignment resolution to misalign objects to be stored in CDB
81 void SetAlignmentResolution(const TClonesArray* misAlignArray, Int_t chId=-1, Double_t chResX=-1., Double_t chResY=-1., Double_t deResX=-1., Double_t deResY=-1.);
85 AliMUONGeometryMisAligner(const AliMUONGeometryMisAligner & right);
87 AliMUONGeometryMisAligner & operator =(const AliMUONGeometryMisAligner &right);
91 // return a misaligned transformation
92 TGeoCombiTrans MisAlignDetElem(const TGeoCombiTrans& transform) const;
93 TGeoCombiTrans MisAlignModule(const TGeoCombiTrans& transform) const;
94 void GetUniMisAlign(Double_t cartMisAlig[3], Double_t angMisAlig[3], const Double_t lParMisAlig[6][2]) const;
95 void GetGausMisAlign(Double_t cartMisAlig[3], Double_t angMisAlig[3], const Double_t lParMisAlig[6][2]) const;
97 Bool_t fUseUni; ///< use uniform distribution for misaligmnets
98 Bool_t fUseGaus; ///< use gaussian distribution for misaligmnets
99 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)
100 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)
102 Double_t fXYAngMisAligFactor; ///< factor (<1) to apply to angular misalignment range since range of motion is restricted out of the xy plane
103 Double_t fZCartMisAligFactor; ///< factor (<1) to apply to cartetian misalignment range since range of motion is restricted in z direction
104 TRandom *fDisplacementGenerator; ///< random number generator for the displacements
106 ClassDef(AliMUONGeometryMisAligner,3) // Geometry parametrisation
109 #endif //ALI_MUON_GEOMETRY_MIS_ALIGNER_H