|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
|
// H4Geom.h
//
// Class which manages geometry information of the super-module
// Author: Giovanni Franzoni
//
#ifndef H4Geom_H
#define H4Geom_H
#include <string>
class H4Geom
{
public:
// Geometry is period dependent. The different geometries are listed below
enum GeomPeriod_t {
Undef, // Cause the program to crash
Year2003, // Test beam during 2003: SM0 and SM1
Spring2004, // June-July 2004: E0' tests
Automn2004 // Tests of one full supermodule
};
// Constants defining supermodule geometry
enum SMGeom_t {
kSModulesInEcal = 36, // Number of Super Modules in whole ECAL
kModules = 4, // Number of modules per supermodule
kTriggerTowers = 68, // Number of trigger towers per supermodule
kTTandMems = 70, // Number of tt per supermodule, including MEM boxes
kTowersInPhi = 4, // Number of trigger towers in phi
kTowersInEta = 17, // Number of trigger towers in eta
kCrystals = 1700, // Number of crystals per supermodule
kCrystalsWithMems = 1750, // Number of channels per supermodule, mem included
kCrystalsM1 = 500, // Number of crystals per supermodule
kCrystalsM2 = 400, // Number of crystals per supermodule
kCrystalsM3 = 400, // Number of crystals per supermodule
kCrystalsM4 = 400, // Number of crystals per supermodule
kCrystalsInPhi = 20, // Number of crystals in phi
kCrystalsInEta = 85, // Number of crystals in eta
kCrystalsInEtaM1 = 25, // Number of crystals in eta for module 1
kCrystalsInEtaM2 = 20, // Number of crystals in eta for module 2
kCrystalsInEtaM3 = 20, // Number of crystals in eta for module 3
kCrystalsInEtaM4 = 20, // Number of crystals in eta for module 4
kCrystalsPerTower = 25, // Number of crystals per trigger tower
kCardsPerTower = 5, // Number of VFE cards per trigger tower
kChannelsPerCard = 5, // Number of channels per VFE card
kSamplesInEvent = 10, // Number of samples in event
kSamplesInPNEvent = 50, // Number of samples in event
kPNs = 10 // number of pn diode for laser monitoring (number [0,9])
};
// Default Constructor, mainly for Root
H4Geom() ;
// Destructor: Does nothing?
virtual ~H4Geom();
// Initialize geometry with config file
bool init();
// Retuns the crystal number in the super module for a given
// tower number in the super module and crystal number in the tower
int getSMCrystalNumber(int tower, int crystal) const ;
// Retuns the crystal number in the super module for a given
// tower number, strip_id number and crystal_id number
// necessary for output of the data-parser, Ecal Monitoring
int getSMCrystalNumber(int tower, int strip_id, int crystal_id) const ;
// Retuns the tower number, strip_id number and crystal_id number
// for a given crystal in the SM numbering
// This is the inverse of int getSMCrystalNumber(int tower, int strip_id, int crystal_id) const ;
void getTowerStripChannelNumber(int& tower, int& strip_id, int& crystal_id, int sm_num) const;
// Retuns the crystal number in a tower for a given
// crystal number in the super module
void getTowerCrystalNumber(int &tower, int &crystal, int smCrystal) const ;
// Returns the crystal number (readout order) in a tower
// for a given position in the tower (crystalNbGeom=0 is the
// lower-right corner and crystalNbGeom=24 is the upper-left corner)
int getTowerCrystalNumber(int smTowerNb, int crystalNbGeom) const ;
// Returns the crystal coordinates (eta, phi index) for a given
// crystal number in the super module
void getCrystalCoord(int &eta, int &phi, int smCrystal) const ;
// Retuns the crystal number in the super module for given coordinates
int getSMCrystalFromCoord(int eta, int phi) const ;
// Returns left neighbour of a sm crystal.
// Input and output are crystal numbers in the super module.
// A negative output means outside of the supermodule.
int getLeft(int smCrystal) const ;
// Returns right neighbour of a sm crystal.
// Input and output are crystal numbers in the super module.
// A negative output means outside of the supermodule.
int getRight(int smCrystal) const ;
// Returns upper neighbour of a sm crystal.
// Input and output are crystal numbers in the super module.
// A negative output means outside of the supermodule.
int getUpper(int smCrystal) const ;
// Returns lower neighbour of a sm crystal.
// Input and output are crystal numbers in the super module.
// A negative output means outside of the supermodule.
int getLower(int smCrystal) const ;
// Returns left neighbour of a crystal referenced by its coordinates.
// New coordonates overwrite the old ones. No check is done to see
// if it corresponds to a real crystal. To be used with caution.
void mvLeft(int &eta, int &phi) const ;
// Returns right neighbour of a crystal referenced by its coordinates.
// New coordonates overwrite the old ones. No check is done to see
// if it corresponds to a real crystal. To be used with caution.
void mvRight(int &eta, int &phi) const ;
// Returns upper neighbour of a crystal referenced by its coordinates.
// New coordonates overwrite the old ones. No check is done to see
// if it corresponds to a real crystal. To be used with caution.
void mvUp(int &eta, int &phi) const ;
// Returns lower neighbour of a crystal referenced by its coordinates.
// New coordonates overwrite the old ones. No check is done to see
// if it corresponds to a real crystal. To be used with caution.
void mvDown(int &eta, int &phi) const ;
// Returns the 25 crystals of tower towerNb in the super module.
// Output are crystal numbers in the super module.
// By default, the order in the output array (tower) corresponds to
// geometric order (index 0 is lower-right corner).
// if order=readout, the order in the output array (tower)
// corresponds to the readout scheme (depends on the kind of tower)
void getTower(int * tower, int towerNb, std::string order = "geom") const ;
// Returns the 5 crystals belonging to the same VFE board as smCrystal.
// Input and output are crystal numbers in the super module.
// By default, the order in the output array (VFE) corresponds to
// The geometric order (index 0 is lower-right corner).
// if order=readout, the order in the output array (VFE)
// corresponds to the readout scheme (depends on the kind of tower)
void getVFE(int * VFE, int smCrystal, std::string order = "geom") const ;
// Returns sm crystal numbers for crystals in a window of
// size widthxheight centered around a given smCrystal.
// width and height must be odd.
// The order in the output array (window) is defined
// by the geometric order (index 0 is lower-right corner).
void getWindow(int * window, int smCrystal, int width, int height) const ;
// Tests if low voltage board is on the right size of the tower.
// Readout scheme depends on that.
bool rightTower(int tower) const ;
// Tests if low voltage board is on the left size of the tower
// Readout scheme depends on that.
bool leftTower(int tower) const ;
// first half of SM, shown by one laser shot
bool isInFirstHalf(int numberInSM)
{
int eta; int phi;
getCrystalCoord(eta,phi, numberInSM);
if (phi > 19 || eta < 20 )
{return true; }
else
{return false; }
};
GeomPeriod_t GetGeomPeriod() const {return geometry_;}
static void SetGeomPeriod(GeomPeriod_t geometry);
int getHalf(int TT);
private:
bool IsGeomPeriodDefined() const;
const static int crystalChannelMap[5][5];
const static int crystalMap[25];
const static int WhichHalf[69];
static GeomPeriod_t geometry_;
};
#endif
|