]> git.uio.no Git - u/mrichter/AliRoot.git/blob - MUON/AliMUONSegmentationV01.h
fWSN->Eval(0.001) to avoid fpe.
[u/mrichter/AliRoot.git] / MUON / AliMUONSegmentationV01.h
1 #ifndef ALIMUONSEGMENTATIONV01_H
2 #define ALIMUONSEGMENTATIONV01_H
3
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6
7 /* $Id$ */
8
9 /////////////////////////////////////////////////////
10 //  Segmentation and Response classes version 01   //
11 /////////////////////////////////////////////////////
12 class AliMUON;
13 class TObjArray;
14
15 #include "AliMUONSegmentationV0.h"
16 #include "TArrayI.h" // because the object, and not the pointer,
17 #include "TArrayF.h" // belongs to the class
18
19
20 class AliMUONSegmentationV01 :
21 public AliMUONSegmentationV0 {
22  public:
23     AliMUONSegmentationV01();
24     AliMUONSegmentationV01(Int_t nsec);
25     AliMUONSegmentationV01(const AliMUONSegmentationV01 & segmentation);
26     
27     virtual ~AliMUONSegmentationV01();
28     
29     //    
30     // Set Chamber Segmentation Parameters
31     // 
32     virtual  void    SetPadDivision(Int_t ndiv[4]);
33     // Radii
34     virtual  void    SetSegRadii(Float_t  r[4]);
35     virtual  void    SetOffsetY(Float_t off) {fOffsetY = off;}
36     //
37     // Transform from pad (wire) to real coordinates and vice versa
38     //
39     // Transform from pad to real coordinates
40     virtual void    GetPadI(Float_t x ,Float_t y ,Int_t   &ix,Int_t &iy);
41     virtual void    GetPadI(Float_t x, Float_t y , Float_t z, Int_t &ix, Int_t &iy);
42     // Transform from real to pad coordinates
43     virtual void    GetPadC(Int_t   ix,Int_t   iy,Float_t &x ,Float_t &y );
44     virtual void    GetPadC(Int_t ix, Int_t iy, Float_t &x, Float_t &y, Float_t &z) 
45         {z=fZ; GetPadC(ix, iy, x , y);}
46     //
47     // Initialisation
48     virtual void Init(Int_t chamber);
49     //
50     // Get member data
51     //
52     // Pad size in x by Sector
53     virtual Float_t Dpx(Int_t isec) const;
54     // Pad size in y by Sector
55     virtual Float_t Dpy(Int_t isec) const;
56     // Max number of Pads in x
57     virtual Int_t   Npx() const {return fNpxS[fNsec-1][1]+1;}
58     //
59     virtual void    SetPad(Int_t ix,Int_t iy);
60     //
61     // Iterate over pads
62     // Initialiser
63     virtual void  FirstPad(Float_t xhit, Float_t yhit, Float_t dx, Float_t dy);
64     virtual void  FirstPad(Float_t xhit, Float_t yhit, Float_t zhit, Float_t dx, Float_t dy);
65     // Stepper
66     virtual void  NextPad();
67     // Condition
68     virtual Int_t MorePads();
69     // Get next neighbours 
70     virtual void Neighbours
71         (Int_t iX, Int_t iY, Int_t* Nlist, Int_t Xlist[10], Int_t Ylist[10]);
72     //
73     // Current Pad during Integration
74     // current sector
75     virtual Int_t ISector()  {return fSector;}
76     // calculate sector from pad coordinates
77     virtual Int_t Sector(Int_t ix, Int_t iy);
78     //
79     // Integration
80     // Current integration limits
81      virtual void IntegrationLimits
82         (Float_t& x1, Float_t& x2, Float_t& y1, Float_t& y2);
83     // Test points for auto calibration
84     void GiveTestPoints(Int_t &n, Float_t *x, Float_t *y) const;
85     //
86     // Draw segmentation zones
87     virtual void Draw(const char *opt="") const;
88     // Function for systematic corrections
89     // Set the correction function
90     virtual void SetCorrFunc(Int_t dum, TF1* func);
91     // Get the correction function
92     virtual TF1* CorrFunc(Int_t iZone) const;
93     // assignment operator
94     AliMUONSegmentationV01& operator=(const AliMUONSegmentationV01& rhs);
95     ClassDef(AliMUONSegmentationV01,1) // Segmentation approximating circular zones with different pad size
96  protected:
97     //  Geometry
98     //
99     Int_t       fNsec;           // Number of sectors
100     TArrayF*    fRSec;           // Sector outer radia
101     TArrayI*    fNDiv;           // Pad size division
102     TArrayF*    fDpxD;           // y pad width per sector
103     Float_t     fOffsetY;        // Staggering offset in y
104     // Segmentation map
105     Int_t      fNpxS[10][1000];  //  Number of pads per sector in x
106     Float_t    fCx[10][1000];    //  pad-sector contour x vs y  
107     // Chamber region consideres during disintegration
108     // (lower left and upper right corner)
109     //
110     Float_t fXmin; // ! lower left  x
111     Float_t fXmax; // ! lower left  y
112     Float_t fYmin; // ! upper right x
113     Float_t fYmax; // ! upper right y 
114
115     //
116     // Current pad during integration (cursor for disintegration)
117     Int_t   fSector;   // ! Current sector
118     //
119     TObjArray *fCorrA; // ! Array of correction functions
120 };
121 #endif
122
123
124
125
126
127
128
129