Particle ID and event selection
Static Public Member Functions | List of all members
ElectronCuts Class Reference

Static Public Member Functions

static boolean CC_nphe_cut (double nphe)
 
static boolean EC_outer_vs_EC_inner_cut (double pcal_energy)
 
static boolean EC_sampling_fraction_cut (double partp, int pcal_sector, double pcal_energy, double ecin_energy, double ecout_energy)
 
static boolean EC_hit_position_fiducial_cut_homogeneous (int pcal_sector, double lv, double lw, Level lvl)
 
static boolean DC_fiducial_cut_XY (int dc_sector, int region, double x, double y, int partpid, boolean isinbending)
 
static boolean DC_z_vertex_cut (int pcal_sector, double partvz, boolean isinbending)
 

Member Function Documentation

◆ CC_nphe_cut()

static boolean CC_nphe_cut ( double  nphe)
inlinestatic

Cherenkov nphe cut

Parameters
nphenumber of photoelectrons
11 {
12 double nphe_min = 2;
13 return nphe > nphe_min;
14 }

◆ DC_fiducial_cut_XY()

static boolean DC_fiducial_cut_XY ( int  dc_sector,
int  region,
double  x,
double  y,
int  partpid,
boolean  isinbending 
)
inlinestatic
Parameters
dc_sectorsector of hits in DC
regionspecify fiducial cuts for which region to use
xx for region 1 or 2 or 3 from REC::Traj
yy for region 1 or 2 or 3 from REC::Traj
partpidpid assigned to particle candidate
isinbendingTrue if magnetic field is inbending
141 {
142
143// new cut parameters for the linear cut based on x and y coordinates (inbending field):
144// replace it in the function: bool DC_fiducial_cut_XY(int j, int region)
145// (optimized for electrons, do not use it for hadrons)
146//
147
148
149 double[][][][] maxparams_in = {
150 { {{-14.563,0.60032},{-19.6768,0.58729},{-22.2531,0.544896}},{{-12.7486,0.587631},{-18.8093,0.571584},{-19.077,0.519895}},
151 {{-11.3481,0.536385},{-18.8912,0.58099},{-18.8584,0.515956}},{{-10.7248,0.52678},{-18.2058,0.559429},{-22.0058,0.53808}},
152 {{-16.9644,0.688637},{-17.1012,0.543961},{-21.3974,0.495489}},{{-13.4454,0.594051},{-19.4173,0.58875},{-22.8771,0.558029}}
153 },
154 { {{-6.2928,0.541828},{-16.7759,0.57962},{-32.5232,0.599023}},{{-6.3996,0.543619},{-16.7429,0.578472},{-32.5408,0.600826}},
155 {{-5.49712,0.53463},{-16.1294,0.576928},{-32.5171,0.597735}},{{-6.4374,0.54839},{-16.9511,0.582143},{-33.0501,0.59995}},
156 {{-5.30128,0.529377},{-16.1229,0.579019},{-30.7768,0.593861}},{{-5.89201,0.541124},{-16.1245,0.575001},{-32.2617,0.601506}}
157 },
158 { {{-6.3618,0.546384},{-17.0277,0.582344},{-34.9276,0.612875}},{{-6.36432,0.546268},{-15.8404,0.574102},{-33.0627,0.599142}},
159 {{-6.34357,0.548411},{-16.0496,0.575913},{-34.8535,0.610211}},{{-5.8568,0.541784},{-16.1124,0.576473},{-32.8547,0.599033}},
160 {{-5.91941,0.536801},{-15.726,0.575211},{-34.0964,0.606777}},{{-5.55498,0.536609},{-15.9853,0.579705},{-33.4886,0.606439}}
161 },
162 { {{-12.594,0.613062},{-18.4504,0.588136},{-16.3157,0.529461}},{{-12.3417,0.61231},{-18.1498,0.590748},{-13.8106,0.52335}},
163 {{-12.1761,0.609307},{-15.919,0.572156},{-13.0598,0.5194}},{{-12.5467,0.612645},{-16.2129,0.572974},{-12.8611,0.51252}},
164 {{-13.0976,0.615928},{-16.9233,0.580972},{-13.0906,0.519738}},{{-12.884,0.622133},{-17.2566,0.585572},{-12.1874,0.510124}}
165 },
166 { {{-6.51157,0.545763},{-16.4246,0.583603},{-32.2001,0.60425}},{{-6.21169,0.541872},{-16.8484,0.591172},{-31.7785,0.606234}},
167 {{-5.89452,0.54464},{-16.612,0.591506},{-29.9143,0.589656}},{{-6.68908,0.553374},{-16.2993,0.585165},{-30.252,0.59519}},
168 {{-6.17185,0.540496},{-16.7197,0.591664},{-31.619,0.608306}},{{-5.7526,0.541761},{-16.2054,0.587326},{-31.3653,0.604081}}
169 },
170 { {{-11.8798,0.62389},{-20.2212,0.610786},{-16.4137,0.51337}},{{-12.0817,0.631621},{-20.7511,0.610844},{-16.9407,0.522958}},
171 {{-9.72746,0.605471},{-20.4903,0.622337},{-15.3363,0.520589}},{{-12.4566,0.627481},{-20.238,0.606098},{-20.7651,0.56974}},
172 {{-11.6712,0.622265},{-18.2649,0.591062},{-19.2569,0.580894}},{{-12.0943,0.630674},{-22.4432,0.633366},{-17.2197,0.537965}}
173 }
174 };
175
176 double[][][][] minparams_in = {
177 { {{12.2692,-0.583057},{17.6233,-0.605722},{19.7018,-0.518429}},{{12.1191,-0.582662},{16.8692,-0.56719},{20.9153,-0.534871}},
178 {{11.4562,-0.53549},{19.3201,-0.590815},{20.1025,-0.511234}},{{13.202,-0.563346},{20.3542,-0.575843},{23.6495,-0.54525}},
179 {{12.0907,-0.547413},{17.1319,-0.537551},{17.861,-0.493782}},{{13.2856,-0.594915},{18.5707,-0.597428},{21.6804,-0.552287}}
180 },
181 { {{5.35616,-0.531295},{16.9702,-0.583819},{36.3388,-0.612192}},{{6.41665,-0.543249},{17.3455,-0.584322},{37.1294,-0.61791}},
182 {{6.86336,-0.550492},{17.2747,-0.575263},{39.6389,-0.625934}},{{6.82938,-0.558897},{17.8618,-0.599931},{39.3376,-0.631517}},
183 {{6.05547,-0.54347},{15.7765,-0.569165},{35.6589,-0.611349}},{{6.3468,-0.544882},{16.7144,-0.578363},{38.2501,-0.617055}}
184 },
185 { {{6.70668,-0.558853},{17.0627,-0.587751},{36.1194,-0.617417}},{{6.3848,-0.542992},{16.6355,-0.581708},{34.6781,-0.609794}},
186 {{6.36802,-0.539521},{15.9829,-0.569165},{32.5691,-0.59588}},{{5.94912,-0.546191},{18.0321,-0.601764},{36.5238,-0.619185}},
187 {{5.65108,-0.541684},{15.5009,-0.567131},{34.0489,-0.602048}},{{6.71064,-0.547956},{16.4449,-0.577051},{34.4375,-0.602515}}
188 },
189 { {{12.4734,-0.608063},{16.1064,-0.575034},{16.0751,-0.536452}},{{12.1936,-0.6034},{15.9302,-0.571271},{14.2791,-0.520157}},
190 {{12.216,-0.600017},{14.8741,-0.56304},{11.1766,-0.498955}},{{12.7941,-0.616044},{17.1516,-0.583616},{11.6077,-0.500028}},
191 {{12.7448,-0.611315},{16.2814,-0.572461},{13.1033,-0.506663}},{{12.7949,-0.612051},{16.1565,-0.569143},{12.9295,-0.504203}}
192 },
193 { {{7.19022,-0.562083},{16.5946,-0.591266},{31.9033,-0.589167}},{{7.80002,-0.571429},{17.8587,-0.595543},{36.5772,-0.630136}},
194 {{7.96121,-0.569485},{17.8085,-0.592936},{37.553,-0.632848}},{{7.52041,-0.566112},{17.3385,-0.603462},{33.7712,-0.606047}},
195 {{7.35796,-0.562782},{15.2865,-0.57433},{29.8283,-0.574685}},{{7.80003,-0.571429},{16.1751,-0.583286},{39.1972,-0.642803}}
196 },
197 { {{13.4466,-0.633911},{22.0097,-0.62205},{18.8862,-0.519652}},{{13.0534,-0.626648},{20.2994,-0.60581},{19.3973,-0.573994}},
198 {{12.547,-0.62145},{18.9322,-0.596491},{16.2331,-0.546036}},{{14.5339,-0.64585},{20.0211,-0.608462},{19.0405,-0.563914}},
199 {{12.7388,-0.617954},{21.1677,-0.621012},{15.4502,-0.525165}},{{13.4019,-0.63075},{16.6584,-0.554797},{19.0302,-0.55004}}
200 }
201 };
202
203
204 //double maxparams_out[6][6][3][2] =
205 double[][][][] maxparams_out = {
206 { {{-9.86221, 0.565985},{-16.4397, 0.569087},{-29.7787, 0.586842}},
207 {{-10.2065, 0.565541},{-16.5554, 0.571394},{-28.933, 0.582078}},
208 {{-8.48034, 0.550706},{-16.4397, 0.569087},{-27.1037, 0.563767}},
209 {{-6.77188, 0.53062},{-16.4397, 0.569087},{-30.485, 0.587534}},
210 {{-8.00705, 0.543502},{-16.4038, 0.571178},{-27.7934, 0.573472}},
211 {{-10.3328, 0.571942},{-16.69, 0.575252},{-30.8177, 0.592418}}
212 },
213 { {{-5.43811, 0.550931},{-17.1906, 0.57936},{-18.552, 0.546789}},
214 {{-5.46281, 0.549659},{-18.0351, 0.588876},{-17.6981, 0.549803}},
215 {{-3.26087, 0.531677},{-16.3762, 0.578005},{-17.6831, 0.55049}},
216 {{-4.5985, 0.542017},{-17.2735, 0.581566},{-16.7013, 0.538853}},
217 {{-6.83053, 0.561019},{-16.5082, 0.579816},{-18.0846, 0.553592}},
218 {{-5.67358, 0.5558},{-18.8196, 0.594965},{-19.4333, 0.560965}}
219 },
220 { {{-12.6317, 0.611023},{-16.5644, 0.578978},{-11.5882, 0.496324}},
221 {{-12.8886, 0.614807},{-17.0847, 0.584072},{-14.9561, 0.532125}},
222 {{-11.4504, 0.600574},{-16.3862, 0.57885},{-12.3309, 0.515431}},
223 {{-12.2256, 0.609801},{-16.2134, 0.574306},{-12.7661, 0.515787}},
224 {{-12.6311, 0.611069},{-16.2486, 0.577577},{-12.6783, 0.519597}},
225 {{-12.6937, 0.615423},{-16.1427, 0.57847},{-11.5156, 0.509458}}
226 },
227 { {{-5.95834, 0.538479},{-15.8909, 0.570164},{-30.2922, 0.586335}},
228 {{-6.15277, 0.542134},{-16.1129, 0.573794},{-31.6024, 0.592681}},
229 {{-6.12341, 0.542023},{-16.1611, 0.575971},{-29.8604, 0.581528}},
230 {{-6.37691, 0.546536},{-16.8501, 0.580239},{-30.0623, 0.580497}},
231 {{-5.96605, 0.537402},{-15.7154, 0.5704},{-31.2955, 0.594146}},
232 {{-5.86704, 0.539556},{-16.2268, 0.580945},{-31.2345, 0.590849}}
233 },
234 { {{-11.7796, 0.614043},{-19.0763, 0.595015},{-18.804, 0.559538}},
235 {{-12.4399, 0.623126},{-19.1733, 0.600646},{-17.675, 0.557016}},
236 {{-10.4158, 0.605483},{-18.0044, 0.595497},{-17.5441, 0.556504}},
237 {{-12.1552, 0.617782},{-19.7134, 0.603519},{-17.3756, 0.549676}},
238 {{-11.3901, 0.612121},{-18.2429, 0.596796},{-10.0097, 0.482578}},
239 {{-12.5004, 0.626384},{-19.9266, 0.60993},{-16.4668, 0.543148}}
240 },
241 { {{-5.60572, 0.537153},{-16.3196, 0.582537},{-32.4336, 0.601487}},
242 {{-5.52369, 0.532985},{-15.2055, 0.568935},{-31.9046, 0.600079}},
243 {{-5.78558, 0.546316},{-16.3328, 0.583765},{-36.0074, 0.617008}},
244 {{-5.82321, 0.542839},{-15.9551, 0.580441},{-31.4304, 0.597132}},
245 {{-5.36526, 0.535923},{-15.9219, 0.586886},{-30.4245, 0.599613}},
246 {{-5.14766, 0.53037},{-14.1986, 0.561504},{-31.7548, 0.60233}}
247 }
248 };
249
250
251 //double minparams_out[6][6][3][2] =
252 double[][][][] minparams_out = {
253 { {{8.07831, -0.548881},{16.4382, -0.569075},{33.7768, -0.607402}},
254 {{8.51057, -0.551773},{16.7782, -0.571381},{32.2613, -0.600686}},
255 {{8.5232, -0.552628},{16.4274, -0.56775},{31.1516, -0.584708}},
256 {{7.98845, -0.544571},{16.4381, -0.569077},{31.8093, -0.595237}},
257 {{7.46705, -0.538557},{16.7414, -0.573345},{31.1888, -0.586751}},
258 {{7.82627, -0.538957},{16.2409, -0.565872},{32.1089, -0.596846}}
259 },
260 { {{7.1519, -0.563678},{16.1038, -0.571795},{20.0449, -0.559802}},
261 {{6.38228, -0.553174},{16.4526, -0.576382},{19.3523, -0.556484}},
262 {{7.11359, -0.561586},{17.2815, -0.578095},{14.9667, -0.53314}},
263 {{5.89053, -0.556406},{17.4946, -0.585038},{17.3607, -0.545739}},
264 {{7.08253, -0.562099},{15.1516, -0.569192},{16.9665, -0.545949}},
265 {{5.53089, -0.546315},{16.4962, -0.574014},{17.9593, -0.545788}}
266 },
267 { {{12.4879, -0.610527},{16.7782, -0.575065},{11.7704, -0.511182}},
268 {{12.1931, -0.604779},{15.6443, -0.560967},{12.7304, -0.515606}},
269 {{12.206, -0.602999},{16.5979, -0.573274},{12.3971, -0.513795}},
270 {{11.5538, -0.604186},{16.6974, -0.576753},{12.7385, -0.517811}},
271 {{12.9718, -0.611968},{17.7233, -0.583943},{10.6601, -0.49233}},
272 {{12.2966, -0.607592},{15.923, -0.564133},{13.9314, -0.525363}}
273 },
274 { {{5.92493, -0.539308},{17.4444, -0.586183},{31.6974, -0.591988}},
275 {{5.467, -0.525876},{16.0649, -0.570869},{30.5937, -0.590071}},
276 {{5.67798, -0.531096},{16.5072, -0.57205},{30.7922, -0.586727}},
277 {{6.85795, -0.558336},{14.9425, -0.545596},{31.3159, -0.592865}},
278 {{6.0155, -0.545283},{16.0649, -0.570869},{30.6644, -0.587002}},
279 {{6.18343, -0.539055},{17.4516, -0.583221},{32.6264, -0.594317}}
280 },
281 { {{12.9118, -0.618907},{19.7061, -0.60171},{18.9352, -0.559461}},
282 {{13.0612, -0.618743},{19.0954, -0.595406},{19.7019, -0.568119}},
283 {{12.4007, -0.613459},{17.544, -0.581147},{12.8175, -0.511017}},
284 {{13.3144, -0.625596},{18.9225, -0.594001},{15.1524, -0.530046}},
285 {{13.101, -0.620887},{18.5616, -0.595279},{14.8807, -0.533111}},
286 {{12.2964, -0.613529},{19.0686, -0.595276},{19.2596, -0.562706}}
287 },
288 { {{5.34118, -0.530584},{16.3015, -0.585185},{38.7808, -0.641362}},
289 {{6.68051, -0.548747},{16.4236, -0.583598},{38.4718, -0.630423}},
290 {{6.87, -0.552602},{16.4285, -0.57977},{36.8889, -0.624053}},
291 {{7.15338, -0.565067},{16.9387, -0.595922},{37.2398, -0.624177}},
292 {{6.06995, -0.550001},{15.7376, -0.577755},{32.6004, -0.601595}},
293 {{6.20459, -0.543148},{14.6326, -0.561623},{39.2154, -0.631762}}
294 }
295 };
296
297
298 double[][][][] minparams = isinbending ? minparams_in : minparams_out;
299 double[][][][] maxparams = isinbending ? maxparams_in : maxparams_out;
300
301 double X = x;
302 double Y = y;
303
304 double X_new = X * Math.cos(Math.toRadians(-60*(dc_sector-1))) - Y * Math.sin(Math.toRadians(-60*(dc_sector-1)));
305 Y = X * Math.sin(Math.toRadians(-60*(dc_sector-1))) + Y * Math.cos(Math.toRadians(-60*(dc_sector-1)));
306 X = X_new;
307
308 int pid = 0;
309
310 switch (partpid)
311 {
312 case 11:
313 pid = 0;
314 break;
315 case 2212:
316 pid = 1;
317 break;
318 case 211:
319 pid = 2;
320 break;
321 case -211:
322 pid = 3;
323 break;
324 case 321:
325 pid = 4;
326 break;
327 case -321:
328 pid = 5;
329 break;
330 default:
331 return false;
332 }
333
334 //if(inbending == true) pid = 0; // use only for electrons in inbending case
335
336 double calc_min = minparams[pid][dc_sector - 1][region-1][0] + minparams[pid][dc_sector - 1][region-1][1] * X;
337 double calc_max = maxparams[pid][dc_sector - 1][region-1][0] + maxparams[pid][dc_sector - 1][region-1][1] * X;
338
339 return (Y > calc_min) && (Y < calc_max);
340 }

◆ DC_z_vertex_cut()

static boolean DC_z_vertex_cut ( int  pcal_sector,
double  partvz,
boolean  isinbending 
)
inlinestatic
Parameters
pcal_sectorsector of hits in PCAL
partvzZ in vertex
isinbendingtrue for inbending, false for outbending
348 {
349
350 double vz_min_sect_inb[] = {-13, -13, -13, -13, -13, -13};
351 double vz_max_sect_inb[] = { 12, 12, 12, 12, 12, 12};
352
353 double vz_min_sect_outb[] = {-20, -20, -20, -20, -20, -20};
354 double vz_max_sect_outb[] = { 12, 12, 12, 12, 12, 12};
355
356 double[] vz_min_sect = new double[6];
357 double[] vz_max_sect = new double[6];
358
359 for(int i = 0; i < 6; i++) {
360 if(isinbending) {
361 vz_min_sect[i] = vz_min_sect_inb[i];
362 vz_max_sect[i] = vz_max_sect_inb[i];
363 } else {
364 vz_min_sect[i] = vz_min_sect_outb[i];
365 vz_max_sect[i] = vz_max_sect_outb[i];
366 }
367 }
368
369 int isec = pcal_sector - 1;
370 double vz_min = vz_min_sect[isec];
371 double vz_max = vz_max_sect[isec];
372
373 return partvz > vz_min && partvz < vz_max;
374 }

◆ EC_hit_position_fiducial_cut_homogeneous()

static boolean EC_hit_position_fiducial_cut_homogeneous ( int  pcal_sector,
double  lv,
double  lw,
Level  lvl 
)
inlinestatic
Parameters
pcal_sectorsector of hits in PCAL
lvpcal_lv
lwpcal_lw
lvlfiducial cut level

inbending:

76 {
77
78 // Cut using the natural directions of the scintillator bars/ fibers:
79
82 //
83 double[] min_v_tight_inb = {19.0, 19.0, 19.0, 19.0, 19.0, 19.0};
84 double[] min_v_med_inb = {14.0, 14.0, 14.0, 14.0, 14.0, 14.0};
85 double[] min_v_loose_inb = {9.0, 9.0, 9.0, 9.0, 9.0, 9.0 };
86 double[] min_v_loosest_inb = {5.0, 5.0, 5.0, 5.0, 5.0, 5.0 };
87 //
88 double[] max_v_tight_inb = {400, 400, 400, 400, 400, 400};
89 double[] max_v_med_inb = {400, 400, 400, 400, 400, 400};
90 double[] max_v_loose_inb = {400, 400, 400, 400, 400, 400};
91 double[] max_v_loosest_inb = {400, 400, 400, 400, 400, 400};
92 //
93 double[] min_w_tight_inb = {19.0, 19.0, 19.0, 19.0, 19.0, 19.0};
94 double[] min_w_med_inb = {14.0, 14.0, 14.0, 14.0, 14.0, 14.0};
95 double[] min_w_loose_inb = {9.0, 9.0, 9.0, 9.0, 9.0, 9.0 };
96 double[] min_w_loosest_inb = {5.0, 5.0, 5.0, 5.0, 5.0, 5.0 };
97 //
98 double[] max_w_tight_inb = {400, 400, 400, 400, 400, 400};
99 double[] max_w_med_inb = {400, 400, 400, 400, 400, 400};
100 double[] max_w_loose_inb = {400, 400, 400, 400, 400, 400};
101 double[] max_w_loosest_inb = {400, 400, 400, 400, 400, 400};
102
103
105
106 int isec = pcal_sector - 1;
107 double min_v = min_v_loose_inb[isec];
108 double max_v = max_v_loose_inb[isec];
109 double min_w = min_w_loose_inb[isec];
110 double max_w = max_w_loose_inb[isec];
111
112 if(lvl == Level.MEDIUM) {
113 min_v = min_v_med_inb[isec];
114 max_v = max_v_med_inb[isec];
115 min_w = min_w_med_inb[isec];
116 max_w = max_w_med_inb[isec];
117 } else if(lvl == Level.TIGHT) {
118 min_v = min_v_tight_inb[isec];
119 max_v = max_v_tight_inb[isec];
120 min_w = min_w_tight_inb[isec];
121 max_w = max_w_tight_inb[isec];
122 } else if(lvl == Level.LOOSEST) {
123 min_v = min_v_loosest_inb[isec];
124 max_v = max_v_loosest_inb[isec];
125 min_w = min_w_loosest_inb[isec];
126 max_w = max_w_loosest_inb[isec];
127 }
128
129 return lv > min_v && lv < max_v && lw > min_w && lw < max_w;
130 }

◆ EC_outer_vs_EC_inner_cut()

static boolean EC_outer_vs_EC_inner_cut ( double  pcal_energy)
inlinestatic

PCAL energy

Parameters
pcal_energyenergy deposited in PCAL
20 {
21 final double edep_tight = 0.06, edep_medium = 0.07, edep_loose = 0.09;
22 return pcal_energy > edep_medium;
23 }

◆ EC_sampling_fraction_cut()

static boolean EC_sampling_fraction_cut ( double  partp,
int  pcal_sector,
double  pcal_energy,
double  ecin_energy,
double  ecout_energy 
)
inlinestatic
Parameters
partpparticle momentum
pcal_sectorsector of hits in PCAL
pcal_energyenergy in PCAL
ecin_energyenergy in ECIN
ecout_energyenergy in ECOUT
33 {
34
35 double[][] ecal_e_sampl_mu = {{ 0.2531, 0.2550, 0.2514, 0.2494, 0.2528, 0.2521 },
36 { -0.6502, -0.7472, -0.7674, -0.4913, -0.3988, -0.703 },
37 { 4.939, 5.350, 5.102, 6.440, 6.149, 4.957 }
38 };
39
40 double[][] ecal_e_sampl_sigm = {{ 2.726e-3, 4.157e-3, 5.222e-3, 5.398e-3, 8.453e-3, 6.533e-3 },
41 { 1.062, 0.859, 0.5564, 0.6576, 0.3242, 0.4423 },
42 { -4.089, -3.318, -2.078, -2.565, -0.8223, -1.274 }
43 };
44
45 double sigma_range = 3.5;
46
47 double ectotal_energy = pcal_energy + ecin_energy + ecout_energy;
48 int isec = pcal_sector-1;
49 double mean = ecal_e_sampl_mu[0][isec] + ecal_e_sampl_mu[1][isec]/1000*Math.pow(partp-ecal_e_sampl_mu[2][isec],2);
50 double sigma = ecal_e_sampl_sigm[0][isec] + ecal_e_sampl_sigm[1][isec]/(10*(partp-ecal_e_sampl_sigm[2][isec]));
51 double upper_lim_total = mean + sigma_range * sigma;
52 double lower_lim_total = mean - sigma_range * sigma;
53
54 boolean pass_band = ectotal_energy/partp <= upper_lim_total && ectotal_energy/partp >= lower_lim_total;
55 boolean pass_triangle = false;
56
57 if(partp < 4.5) {
58 pass_triangle = true;
59 }
60 else {
61 pass_triangle = ecin_energy/partp > (0.2 - pcal_energy/partp);
62 }
63
64 return pass_band && pass_triangle;
65 }
The documentation for this class was generated from the following file:
Generated by doxygen 1.9.3