Fork me on GitHub

source: git/cards/ILCgen/ILCgen_HCAL_Binning.tcl@ f76741a

Last change on this file since f76741a was d192152, checked in by Aleksander Filip Zarnecki <zarnecki@…>, 4 years ago

Generic ILC detector model for Delphes. First import from iLCSoft/ILCDelphes

  • Property mode set to 100644
File size: 1.8 KB
Line 
1#P.Sopicki: based on plots from D.Jeans:
2# corrected by A.F.Zarnecki
3
4set pi [expr {acos(-1)} ]
5
6# HCAL barrel
7# 0.015 segm in eta and phi(rad) - towers (5x5 mm^2) for |eta| <= 1.1
8
9# to get dphi = 0.015rad/0.86deg -> 360/0.86 /cosh(0.55) = ~362.7 -> 360
10# to get deta = 0.015 2.2/0.015=146.(6)
11set PhiBins {}
12 for {set i -180} {$i <= 180} {incr i} {
13 add PhiBins [expr {$i * $pi/180.0} ]
14 }
15
16for {set i 1} {$i <= 147} {incr i} {
17 set eta [expr {-1.1 + $i * 2.2/147.0} ]
18 add EtaPhiBins $eta $PhiBins
19 }
20
21# HCAL endcaps
22# eta 1.1 - 2 : dphi(rad) and deta = 0.03
23# eta 2 - 2.5 : d = 0.065
24# eta 2.5 - 2.8 : d = 0.08
25
26#360/1.72 deg/cosh(1.65) : ~77.6 -> 76 divisions
27#0.9/0.03: 30
28set PhiBins {}
29 for {set i -38} {$i <= 38} {incr i} {
30 add PhiBins [expr {$i * $pi/38} ]
31 }
32
33for {set i 1} {$i <= 30} {incr i} {
34 set eta [expr {-2.0 + $i * 0.9/30.0} ]
35 add EtaPhiBins $eta $PhiBins
36 }
37
38for {set i 1} {$i <= 30} {incr i} {
39 set eta [expr {1.1 + $i * 0.9/30.0} ]
40 add EtaPhiBins $eta $PhiBins
41 }
42
43#360/3.72 deg/cosh(2.25) : ~20 divisions
44#0.5/0.065: ~7.65 -> 8
45set PhiBins {}
46 for {set i -20} {$i <= 20} {incr i} {
47 add PhiBins [expr {$i * $pi/20} ]
48 }
49
50for {set i 1} {$i <= 8} {incr i} {
51 set eta [expr {-2.5 + $i * 0.5/8.0} ]
52 add EtaPhiBins $eta $PhiBins
53 }
54
55for {set i 1} {$i <= 8} {incr i} {
56 set eta [expr {2.0 + $i * 0.5/8.0} ]
57 add EtaPhiBins $eta $PhiBins
58}
59
60#360/4.58 deg/cosh(2.65) : ~11 -> 10 divisions
61#0.3/0.08: ~4
62set PhiBins {}
63 for {set i -5} {$i <= 5} {incr i} {
64 add PhiBins [expr {$i * $pi/5} ]
65 }
66
67for {set i 0} {$i <= 4} {incr i} {
68 set eta [expr {-2.8 + $i * 0.3/4.0} ]
69 add EtaPhiBins $eta $PhiBins
70 }
71
72for {set i 1} {$i <= 4} {incr i} {
73 set eta [expr {2.5 + $i * 0.3/4.0} ]
74 add EtaPhiBins $eta $PhiBins
75 }
Note: See TracBrowser for help on using the repository browser.