source: trunk/examples/advanced/Rich/src/RichTbGeometryParameters.cc @ 1282

Last change on this file since 1282 was 807, checked in by garnier, 16 years ago

update
File size: 11.2 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer                                           *
4// *                                                                  *
5// * The  Geant4 software  is  copyright of the Copyright Holders  of *
6// * the Geant4 Collaboration.  It is provided  under  the terms  and *
7// * conditions of the Geant4 Software License,  included in the file *
8// * LICENSE and available at  http://cern.ch/geant4/license .  These *
9// * include a list of copyright holders.                             *
10// *                                                                  *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work  make  any representation or  warranty, express or implied, *
14// * regarding  this  software system or assume any liability for its *
15// * use.  Please see the license in the file  LICENSE  and URL above *
16// * for the full disclaimer and the limitation of liability.         *
17// *                                                                  *
18// * This  code  implementation is the result of  the  scientific and *
19// * technical work of the GEANT4 collaboration.                      *
20// * By using,  copying,  modifying or  distributing the software (or *
21// * any work based  on the software)  you  agree  to acknowledge its *
22// * use  in  resulting  scientific  publications,  and indicate your *
23// * acceptance of all terms of the Geant4 Software license.          *
24// ********************************************************************
25//
26// Rich advanced example for Geant4
27// RichTbAnalysisManager.cc for Rich of LHCb
28// History:
29// Created: Sajan Easo (Sajan.Easo@cern.ch)
30// Revision and changes: Patricia Mendez (Patricia.Mendez@cern.ch)
31/////////////////////////////////////////////////////////////////////////////
32#include "globals.hh"
33#include "RichTbGeometryParameters.hh"
34
35G4double PixRowNumSect[NumberOfPadHpdSiPixels];
36G4double PixColNumSect[NumberOfPadHpdSiPixels];
37G4bool PixelAtSectEdge[NumberOfPadHpdSiPixels];
38
39G4double GetCurAerogelLength(G4int Aerognum ) {
40  return 2.0*AgelHalfZ[Aerognum];
41}
42RichTbPadHpdSiPixPos::RichTbPadHpdSiPixPos(G4int ipixelnum) {
43  //coord system such that Origin at the top of a sector;
44  // which is near the center of the anode.
45  //The X is along the bottom and Y goes radially outward (from top to
46  //     bottom).
47 
48
49
50  icurpixel=ipixelnum;
51  PadHpdSiPixPosX=PixColNumSect[ipixelnum]*XsizePix;
52   if(ipixelnum ==  BigPixelNum ) {
53  PadHpdSiPixPosY=PixRowNumSect[ipixelnum]*YsizePix+RowInitPointBigPixel;
54  }else {
55  PadHpdSiPixPosY=PixRowNumSect[ipixelnum]*YsizePix+RowInitPoint;
56
57  }
58}
59RichTbPadHpdSiPixPos::~RichTbPadHpdSiPixPos(){ ; }
60void InitializeRichTbGeometry(){
61  PixRowNumSect[59-1]=0.0;
62  PixRowNumSect[70-1]=1.0;
63  PixRowNumSect[60-1]=2.0;
64  PixRowNumSect[82-1]=3.0;
65  PixRowNumSect[47-1]=3.0;
66  PixRowNumSect[81-1]=4.0;
67  PixRowNumSect[48-1]=4.0;
68 
69  PixRowNumSect[91-1]=5.0;
70  PixRowNumSect[69-1]=5.0;
71  PixRowNumSect[38-1]=5.0;
72 
73  PixRowNumSect[101-1]=6.0;
74  PixRowNumSect[80-1]=6.0;
75  PixRowNumSect[49-1]=6.0;
76  PixRowNumSect[28-1]=6.0;
77 
78  PixRowNumSect[100-1]=7.0;
79  PixRowNumSect[79-1]=7.0;
80  PixRowNumSect[50-1]=7.0;
81  PixRowNumSect[29-1]=7.0;
82 
83  PixRowNumSect[99-1]=8.0;
84  PixRowNumSect[78-1]=8.0;
85  PixRowNumSect[51-1]=8.0;
86  PixRowNumSect[30-1]=8.0;
87 
88  PixRowNumSect[109-1]=9.0;
89  PixRowNumSect[90-1]=9.0;
90  PixRowNumSect[61-1]=9.0;
91  PixRowNumSect[39-1]=9.0;
92  PixRowNumSect[20-1]=9.0;
93 
94  PixRowNumSect[108-1]=10.0;
95  PixRowNumSect[89-1]=10.0;
96  PixRowNumSect[68-1]=10.0;
97  PixRowNumSect[40-1]=10.0;
98  PixRowNumSect[21-1]=10.0;
99 
100  PixRowNumSect[107-1]=11.0;
101  PixRowNumSect[88-1]=11.0;
102  PixRowNumSect[62-1]=11.0;
103  PixRowNumSect[41-1]=11.0;
104  PixRowNumSect[22-1]=11.0;
105 
106  PixRowNumSect[116-1]=12.0;
107  PixRowNumSect[98-1]=12.0;
108  PixRowNumSect[77-1]=12.0;
109  PixRowNumSect[52-1]=12.0;
110  PixRowNumSect[31-1]=12.0;
111  PixRowNumSect[13-1]=12.0;
112 
113  PixRowNumSect[115-1]=13.0;
114  PixRowNumSect[97-1]=13.0;
115  PixRowNumSect[76-1]=13.0;
116  PixRowNumSect[53-1]=13.0;
117  PixRowNumSect[32-1]=13.0;
118  PixRowNumSect[14-1]=13.0;
119 
120  PixRowNumSect[114-1]=14.0;
121  PixRowNumSect[96-1]=14.0;
122  PixRowNumSect[75-1]=14.0;
123  PixRowNumSect[54-1]=14.0;
124  PixRowNumSect[33-1]=14.0;
125  PixRowNumSect[15-1]=14.0;
126 
127  PixRowNumSect[121-1]=15.0;
128  PixRowNumSect[106-1]=15.0;
129  PixRowNumSect[87-1]=15.0;
130  PixRowNumSect[67-1]=15.0;
131  PixRowNumSect[42-1]=15.0;
132  PixRowNumSect[23-1]=15.0;
133  PixRowNumSect[8-1]=15.0;
134 
135  PixRowNumSect[120-1]=16.0;
136  PixRowNumSect[105-1]=16.0;
137  PixRowNumSect[86-1]=16.0;
138  PixRowNumSect[63-1]=16.0;
139  PixRowNumSect[43-1]=16.0;
140  PixRowNumSect[24-1]=16.0;
141  PixRowNumSect[9-1]=16.0;
142 
143  PixRowNumSect[119-1]=17.0;
144  PixRowNumSect[104-1]=17.0;
145  PixRowNumSect[85-1]=17.0;
146  PixRowNumSect[66-1]=17.0;
147  PixRowNumSect[44-1]=17.0;
148  PixRowNumSect[25-1]=17.0;
149  PixRowNumSect[10-1]=17.0;
150 
151  PixRowNumSect[125-1]=18.0;
152  PixRowNumSect[113-1]=18.0;
153  PixRowNumSect[95-1]=18.0;
154  PixRowNumSect[74-1]=18.0;
155  PixRowNumSect[55-1]=18.0;
156  PixRowNumSect[34-1]=18.0;
157  PixRowNumSect[16-1]=18.0;
158  PixRowNumSect[4-1]=18.0;
159 
160  PixRowNumSect[124-1]=19.0;
161  PixRowNumSect[112-1]=19.0;
162  PixRowNumSect[94-1]=19.0;
163  PixRowNumSect[73-1]=19.0;
164  PixRowNumSect[56-1]=19.0;
165  PixRowNumSect[35-1]=19.0;
166  PixRowNumSect[17-1]=19.0;
167  PixRowNumSect[5-1]=19.0;
168 
169  PixRowNumSect[123-1]=20.0;
170  PixRowNumSect[111-1]=20.0;
171  PixRowNumSect[93-1]=20.0;
172  PixRowNumSect[72-1]=20.0;
173  PixRowNumSect[57-1]=20.0;
174  PixRowNumSect[36-1]=20.0;
175  PixRowNumSect[18-1]=20.0;
176  PixRowNumSect[6-1]=20.0;
177 
178  PixRowNumSect[127-1]=21.0;
179  PixRowNumSect[118-1]=21.0;
180  PixRowNumSect[103-1]=21.0;
181  PixRowNumSect[84-1]=21.0;
182  PixRowNumSect[65-1]=21.0;
183  PixRowNumSect[45-1]=21.0;
184  PixRowNumSect[26-1]=21.0;
185  PixRowNumSect[11-1]=21.0;
186  PixRowNumSect[2-1]=21.0;
187 
188  PixRowNumSect[126-1]=22.0;
189  PixRowNumSect[117-1]=22.0;
190  PixRowNumSect[102-1]=22.0;
191  PixRowNumSect[83-1]=22.0;
192  PixRowNumSect[64-1]=22.0;
193  PixRowNumSect[46-1]=22.0;
194  PixRowNumSect[27-1]=22.0;
195  PixRowNumSect[12-1]=22.0;
196  PixRowNumSect[3-1]=22.0;
197 
198  PixRowNumSect[128-1]=23.0;
199  PixRowNumSect[122-1]=23.0;
200  PixRowNumSect[110-1]=23.0;
201  PixRowNumSect[92-1]=23.0;
202  PixRowNumSect[71-1]=23.0;
203  PixRowNumSect[58-1]=23.0;
204  PixRowNumSect[37-1]=23.0;
205  PixRowNumSect[19-1]=23.0;
206  PixRowNumSect[7-1]=23.0;
207  PixRowNumSect[1-1]=23.0;
208
209
210
211  PixColNumSect[59-1]=0.0;
212  PixColNumSect[70-1]=0.0;
213  PixColNumSect[60-1]=0.0;
214  PixColNumSect[82-1]=-0.5;
215  PixColNumSect[47-1]=0.5;
216  PixColNumSect[81-1]=-0.5;
217  PixColNumSect[48-1]=0.5;
218
219  PixColNumSect[91-1]=-1.0;
220  PixColNumSect[69-1]=0.0;
221  PixColNumSect[38-1]=1.0;
222
223  PixColNumSect[101-1]=-1.5;
224  PixColNumSect[80-1]=-0.5;
225  PixColNumSect[49-1]=0.5;
226  PixColNumSect[28-1]=1.5;
227
228  PixColNumSect[100-1]=-1.5;
229  PixColNumSect[79-1]=-0.5;
230  PixColNumSect[50-1]=0.5;
231  PixColNumSect[29-1]=1.5;
232
233  PixColNumSect[99-1]=-1.5;
234  PixColNumSect[78-1]=-0.5;
235  PixColNumSect[51-1]=0.5;
236  PixColNumSect[30-1]=1.5;
237
238  PixColNumSect[109-1]=-2.0;
239  PixColNumSect[90-1]=-1.0;
240  PixColNumSect[61-1]=0.0;
241  PixColNumSect[39-1]=1.0;
242  PixColNumSect[20-1]=2.0;
243
244  PixColNumSect[108-1]=-2.0;
245  PixColNumSect[89-1]=-1.0;
246  PixColNumSect[68-1]=0.0;
247  PixColNumSect[40-1]=1.0;
248  PixColNumSect[21-1]=2.0;
249
250  PixColNumSect[107-1]=-2.0;
251  PixColNumSect[88-1]=-1.0;
252  PixColNumSect[62-1]=0.0;
253  PixColNumSect[41-1]=1.0;
254  PixColNumSect[22-1]=2.0;
255
256  PixColNumSect[116-1]=-2.5;
257  PixColNumSect[98-1]=-1.5;
258  PixColNumSect[77-1]=-0.5;
259  PixColNumSect[52-1]=0.5;
260  PixColNumSect[31-1]=1.5;
261  PixColNumSect[13-1]=2.5;
262
263  PixColNumSect[115-1]=-2.5;
264  PixColNumSect[97-1]=-1.5;
265  PixColNumSect[76-1]=-0.5;
266  PixColNumSect[53-1]=0.5;
267  PixColNumSect[32-1]=1.5;
268  PixColNumSect[14-1]=2.5;
269
270  PixColNumSect[114-1]=-2.5;
271  PixColNumSect[96-1]=-1.5;
272  PixColNumSect[75-1]=-0.5;
273  PixColNumSect[54-1]=0.5;
274  PixColNumSect[33-1]=1.5;
275  PixColNumSect[15-1]=2.5;
276
277  PixColNumSect[121-1]=-3.0;
278  PixColNumSect[106-1]=-2.0;
279  PixColNumSect[87-1]=-1.0;
280  PixColNumSect[67-1]=0.0;
281  PixColNumSect[42-1]=1.0;
282  PixColNumSect[23-1]=2.0;
283  PixColNumSect[8-1]=3.0;
284
285  PixColNumSect[120-1]=-3.0;
286  PixColNumSect[105-1]=-2.0;
287  PixColNumSect[86-1]=-1.0;
288  PixColNumSect[63-1]=0.0;
289  PixColNumSect[43-1]=1.0;
290  PixColNumSect[24-1]=2.0;
291  PixColNumSect[9-1]=3.0;
292
293  PixColNumSect[119-1]=-3.0;
294  PixColNumSect[104-1]=-2.0;
295  PixColNumSect[85-1]=-1.0;
296  PixColNumSect[66-1]=0.0;
297  PixColNumSect[44-1]=1.0;
298  PixColNumSect[25-1]=2.0;
299  PixColNumSect[10-1]=3.0;
300
301  PixColNumSect[125-1]=-3.5;
302  PixColNumSect[113-1]=-2.5;
303  PixColNumSect[95-1]=-1.5;
304  PixColNumSect[74-1]=-0.5;
305  PixColNumSect[55-1]=0.5;
306  PixColNumSect[34-1]=1.5;
307  PixColNumSect[16-1]=2.5;
308  PixColNumSect[4-1]=3.5;
309
310  PixColNumSect[124-1]=-3.5;
311  PixColNumSect[112-1]=-2.5;
312  PixColNumSect[94-1]=-1.5;
313  PixColNumSect[73-1]=-0.5;
314  PixColNumSect[56-1]=0.5;
315  PixColNumSect[35-1]=1.5;
316  PixColNumSect[17-1]=2.5;
317  PixColNumSect[5-1]=3.5;
318
319  PixColNumSect[123-1]=-3.5;
320  PixColNumSect[111-1]=-2.5;
321  PixColNumSect[93-1]=-1.5;
322  PixColNumSect[72-1]=-0.5;
323  PixColNumSect[57-1]=0.5;
324  PixColNumSect[36-1]=1.5;
325  PixColNumSect[18-1]=2.5;
326  PixColNumSect[6-1]=3.5;
327
328  PixColNumSect[127-1]=-4.0;
329  PixColNumSect[118-1]=-3.0;
330  PixColNumSect[103-1]=-2.0;
331  PixColNumSect[84-1]=-1.0;
332  PixColNumSect[65-1]=0.0;
333  PixColNumSect[45-1]=1.0;
334  PixColNumSect[26-1]=2.0;
335  PixColNumSect[11-1]=3.0;
336  PixColNumSect[2-1]=4.0;
337
338  PixColNumSect[126-1]=-4.0;
339  PixColNumSect[117-1]=-3.0;
340  PixColNumSect[102-1]=-2.0;
341  PixColNumSect[83-1]=-1.0;
342  PixColNumSect[64-1]=0.0;
343  PixColNumSect[46-1]=1.0;
344  PixColNumSect[27-1]=2.0;
345  PixColNumSect[12-1]=3.0;
346  PixColNumSect[3-1]=4.0;
347
348  PixColNumSect[128-1]=-4.5;
349  PixColNumSect[122-1]=-3.5;
350  PixColNumSect[110-1]=-2.5;
351  PixColNumSect[92-1]=-1.5;
352  PixColNumSect[71-1]=-0.5;
353  PixColNumSect[58-1]=0.5;
354  PixColNumSect[37-1]=1.5;
355  PixColNumSect[19-1]=2.5;
356  PixColNumSect[7-1]=3.5;
357  PixColNumSect[1-1]=4.5;
358
359for (G4int ipixel=0;ipixel<NumberOfPadHpdSiPixels; ipixel++){
360  PixelAtSectEdge[ipixel]=false;
361  }
362 
363PixelAtSectEdge[59-1]=true;
364PixelAtSectEdge[70-1]=true;
365PixelAtSectEdge[59-1]=true;
366PixelAtSectEdge[60-1]=true;
367PixelAtSectEdge[82-1]=true;
368PixelAtSectEdge[47-1]=true;
369PixelAtSectEdge[81-1]=true;
370PixelAtSectEdge[48-1]=true;
371PixelAtSectEdge[91-1]=true;
372PixelAtSectEdge[38-1]=true;
373PixelAtSectEdge[101-1]=true;
374PixelAtSectEdge[28-1]=true;
375
376PixelAtSectEdge[100-1]=true;
377PixelAtSectEdge[29-1]=true;
378
379PixelAtSectEdge[99-1]=true;
380PixelAtSectEdge[30-1]=true;
381
382PixelAtSectEdge[109-1]=true;
383PixelAtSectEdge[20-1]=true;
384
385PixelAtSectEdge[108-1]=true;
386PixelAtSectEdge[21-1]=true;
387
388PixelAtSectEdge[107-1]=true;
389PixelAtSectEdge[22-1]=true;
390
391PixelAtSectEdge[116-1]=true;
392PixelAtSectEdge[13-1]=true;
393
394PixelAtSectEdge[115-1]=true;
395PixelAtSectEdge[14-1]=true;
396
397PixelAtSectEdge[114-1]=true;
398PixelAtSectEdge[15-1]=true;
399
400PixelAtSectEdge[121-1]=true;
401PixelAtSectEdge[8-1]=true;
402
403PixelAtSectEdge[120-1]=true;
404PixelAtSectEdge[9-1]=true;
405
406PixelAtSectEdge[119-1]=true;
407PixelAtSectEdge[10-1]=true;
408
409PixelAtSectEdge[125-1]=true;
410PixelAtSectEdge[4-1]=true;
411
412PixelAtSectEdge[124-1]=true;
413PixelAtSectEdge[5-1]=true;
414
415PixelAtSectEdge[123-1]=true;
416PixelAtSectEdge[6-1]=true;
417
418PixelAtSectEdge[127-1]=true;
419PixelAtSectEdge[2-1]=true;
420
421PixelAtSectEdge[126-1]=true;
422PixelAtSectEdge[3-1]=true;
423
424PixelAtSectEdge[128-1]=true;
425PixelAtSectEdge[1-1]=true;
426   
427   
428}
429
430
431
432
433
434
435
Note: See TracBrowser for help on using the repository browser.