Fork me on GitHub

source: git/cards/ILCgen/ILCgen_CTagging_30.tcl@ 4564cba

Last change on this file since 4564cba 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: 6.4 KB
Line 
1 # misidentification rate (uds)
2 # average efficiency: 0.0168
3 add EfficiencyFormula {0} {
4 ( abs(eta)<0.867 )*( energy<30 )*( 0.0294 )+
5 ( abs(eta)<0.867 )*( 30<=energy && energy<60 )*( 0.0285 )+
6 ( abs(eta)<0.867 )*( 60<=energy && energy<100 )*( 0.0286 )+
7 ( abs(eta)<0.867 )*( 100<=energy && energy<150 )*( 0.0301 )+
8 ( abs(eta)<0.867 )*( 150<=energy )*( 0.0311 )+
9 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( energy<30 )*( 0.0262 )+
10 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 30<=energy && energy<60 )*( 0.0239 )+
11 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 60<=energy && energy<100 )*( 0.0235 )+
12 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 100<=energy && energy<150 )*( 0.0244 )+
13 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 150<=energy )*( 0.0273 )+
14 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( energy<30 )*( 0.0197 )+
15 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 30<=energy && energy<60 )*( 0.0165 )+
16 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 60<=energy && energy<100 )*( 0.0183 )+
17 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 100<=energy && energy<150 )*( 0.0216 )+
18 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 150<=energy )*( 0.0221 )+
19 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( energy<30 )*( 0.0103 )+
20 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 30<=energy && energy<60 )*( 0.0125 )+
21 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 60<=energy && energy<100 )*( 0.0116 )+
22 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 100<=energy && energy<150 )*( 0.0138 )+
23 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 150<=energy )*( 0.0167 )+
24 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( energy<30 )*( 0.00848 )+
25 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 30<=energy && energy<60 )*( 0.00687 )+
26 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 60<=energy && energy<100 )*( 0.00718 )+
27 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 100<=energy && energy<150 )*( 0.00868 )+
28 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 150<=energy )*( 0.00706 )+
29 ( 2.99<=abs(eta) )*( energy<30 )*( 0.00685 )+
30 ( 2.99<=abs(eta) )*( 30<=energy && energy<60 )*( 0.00459 )+
31 ( 2.99<=abs(eta) )*( 60<=energy && energy<100 )*( 0.00485 )+
32 ( 2.99<=abs(eta) )*( 100<=energy && energy<150 )*( 0.00657 )+
33 ( 2.99<=abs(eta) )*( 150<=energy )*( 0.0083 )
34 }
35
36 # c-tagging efficiency
37 # average efficiency: 0.315
38 add EfficiencyFormula {4} {
39 ( abs(eta)<0.867 )*( energy<30 )*( 0.497 )+
40 ( abs(eta)<0.867 )*( 30<=energy && energy<60 )*( 0.598 )+
41 ( abs(eta)<0.867 )*( 60<=energy && energy<100 )*( 0.652 )+
42 ( abs(eta)<0.867 )*( 100<=energy && energy<150 )*( 0.686 )+
43 ( abs(eta)<0.867 )*( 150<=energy )*( 0.702 )+
44 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( energy<30 )*( 0.351 )+
45 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 30<=energy && energy<60 )*( 0.48 )+
46 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 60<=energy && energy<100 )*( 0.565 )+
47 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 100<=energy && energy<150 )*( 0.614 )+
48 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 150<=energy )*( 0.649 )+
49 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( energy<30 )*( 0.194 )+
50 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 30<=energy && energy<60 )*( 0.321 )+
51 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 60<=energy && energy<100 )*( 0.432 )+
52 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 100<=energy && energy<150 )*( 0.517 )+
53 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 150<=energy )*( 0.567 )+
54 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( energy<30 )*( 0.0942 )+
55 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 30<=energy && energy<60 )*( 0.161 )+
56 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 60<=energy && energy<100 )*( 0.239 )+
57 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 100<=energy && energy<150 )*( 0.313 )+
58 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 150<=energy )*( 0.37 )+
59 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( energy<30 )*( 0.0216 )+
60 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 30<=energy && energy<60 )*( 0.0524 )+
61 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 60<=energy && energy<100 )*( 0.0616 )+
62 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 100<=energy && energy<150 )*( 0.075 )+
63 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 150<=energy )*( 0.0988 )+
64 ( 2.99<=abs(eta) )*( energy<30 )*( 0.00502 )+
65 ( 2.99<=abs(eta) )*( 30<=energy && energy<60 )*( 0.0135 )+
66 ( 2.99<=abs(eta) )*( 60<=energy && energy<100 )*( 0.0276 )+
67 ( 2.99<=abs(eta) )*( 100<=energy && energy<150 )*( 0.0428 )+
68 ( 2.99<=abs(eta) )*( 150<=energy )*( 0.0402 )
69 }
70
71 # misidentification rate (b)
72 # average efficiency: 0.105
73 add EfficiencyFormula {5} {
74 ( abs(eta)<0.867 )*( energy<30 )*( 0.212 )+
75 ( abs(eta)<0.867 )*( 30<=energy && energy<60 )*( 0.143 )+
76 ( abs(eta)<0.867 )*( 60<=energy && energy<100 )*( 0.112 )+
77 ( abs(eta)<0.867 )*( 100<=energy && energy<150 )*( 0.0862 )+
78 ( abs(eta)<0.867 )*( 150<=energy )*( 0.0749 )+
79 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( energy<30 )*( 0.215 )+
80 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 30<=energy && energy<60 )*( 0.148 )+
81 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 60<=energy && energy<100 )*( 0.117 )+
82 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 100<=energy && energy<150 )*( 0.0992 )+
83 ( 0.867<=abs(eta) && abs(eta)<1.47 )*( 150<=energy )*( 0.0827 )+
84 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( energy<30 )*( 0.182 )+
85 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 30<=energy && energy<60 )*( 0.161 )+
86 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 60<=energy && energy<100 )*( 0.119 )+
87 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 100<=energy && energy<150 )*( 0.107 )+
88 ( 1.47<=abs(eta) && abs(eta)<1.83 )*( 150<=energy )*( 0.0884 )+
89 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( energy<30 )*( 0.137 )+
90 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 30<=energy && energy<60 )*( 0.136 )+
91 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 60<=energy && energy<100 )*( 0.133 )+
92 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 100<=energy && energy<150 )*( 0.108 )+
93 ( 1.83<=abs(eta) && abs(eta)<2.3 )*( 150<=energy )*( 0.101 )+
94 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( energy<30 )*( 0.0688 )+
95 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 30<=energy && energy<60 )*( 0.0724 )+
96 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 60<=energy && energy<100 )*( 0.0761 )+
97 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 100<=energy && energy<150 )*( 0.0751 )+
98 ( 2.3<=abs(eta) && abs(eta)<2.99 )*( 150<=energy )*( 0.064 )+
99 ( 2.99<=abs(eta) )*( energy<30 )*( 0.0264 )+
100 ( 2.99<=abs(eta) )*( 30<=energy && energy<60 )*( 0.0337 )+
101 ( 2.99<=abs(eta) )*( 60<=energy && energy<100 )*( 0.0504 )+
102 ( 2.99<=abs(eta) )*( 100<=energy && energy<150 )*( 0.0501 )+
103 ( 2.99<=abs(eta) )*( 150<=energy )*( 0.056 )
104 }
Note: See TracBrowser for help on using the repository browser.