BOSS 7.1.3
BESIII Offline Software System
Loading...
Searching...
No Matches
SimplePIDSvc.cxx
Go to the documentation of this file.
1#include "GaudiKernel/IDataProviderSvc.h"
2#include "GaudiKernel/MsgStream.h"
3#include "GaudiKernel/SmartDataPtr.h"
4#include "EventModel/EventModel.h"
5#include "EvtRecEvent/EvtRecEvent.h"
6#include "EventModel/EventHeader.h"
7#include "SimplePIDSvc/SimplePIDSvc.h"
8
9#include "DstEvent/TofHitStatus.h"
10#include "TMath.h"
11#include "TFile.h"
12#include "TMatrixD.h"
13#include "TArray.h"
14#include <fstream>
15#include <iostream>
16#include <cstdlib>
17#include <cmath>
18using namespace std;
19DECLARE_COMPONENT(SimplePIDSvc)
20
21SimplePIDSvc::SimplePIDSvc(const std::string& name, ISvcLocator* svcLoc) : base_class(name, svcLoc)
22{
23 declareProperty("DedxChiCut", m_dedx_chi_cut = 4);
24 declareProperty("TofChiCut", m_tof_chi_cut = 4);
25 declareProperty("IsTofCorr", m_tof_corr = true);
26 declareProperty("IsDedxCorr", m_dedx_corr = true);
27 declareProperty("EidRatio", m_eid_ratio = 0.80);
28}
29
30SimplePIDSvc::~SimplePIDSvc() {;}
31
32StatusCode SimplePIDSvc::initialize()
33{
34 MsgStream log(messageService(), name());
35 log << MSG::INFO << "in SimplePIDSvc initialize()" << endreq;
36
37 StatusCode sc = Service::initialize();
38
39 sc = serviceLocator()->service("EventDataSvc", eventSvc_, true);
40
41 loadHistogram();
42 loadSecondPar();//the second correct factor
43
44 return sc;
45}
46
47StatusCode SimplePIDSvc::finalize()
48{
49 MsgStream log(messageService(), name());
50 log << MSG::INFO << "in SimplePIDSvc finalize()" << endreq;
51
52 StatusCode sc = Service::finalize();
53
54 for (unsigned int i = 0; i < 2; i++)
55 {
56 for (unsigned int j = 0; j < 4; j++)
57 {
58 f_dedx[i][j]->Close();
59 f_tof_q[i][j]->Close();
60 f_tof_bgcost[i][j]->Close();
61 f_tof_wgt[i][j]->Close();
62 f_tof_final[i][j]->Close();
63 f_tofec_q[i][j]->Close();
64 f_tofec_bg[i][j]->Close();
65 f_tofec_cost[i][j]->Close();
66 }
67 }
68 for (unsigned int i = 0; i < 3; i++)
69 {
70 for (unsigned int j = 0; j < 4; j++)
71 {
72 f_emc[i][j]->Close();
73 }
74 }
75
76 return sc;
77}
78
79/*StatusCode SimplePIDSvc::queryInterface(const InterfaceID& riid, void** ppvIF)
80 {
81 if (ISimplePIDSvc::interfaceID().versionMatch(riid))
82 {
83 *ppvIF = dynamic_cast<ISimplePIDSvc*>(this);
84 }
85 else
86 {
87 return Service::queryInterface(riid, ppvIF);
88 }
89 addRef();
90 return StatusCode::SUCCESS;
91 }
92 */
93
94
95void SimplePIDSvc::preparePID(EvtRecTrack* track)
96{
97
98 SmartDataPtr<Event::EventHeader> eventHeaderpid(eventSvc_, "/Event/EventHeader");
99 m_run = eventHeaderpid->runNumber();
100
101 if (track->isMdcKalTrackValid())
102 {
103 RecMdcKalTrack *mdckalTrk = track->mdcKalTrack();
104 RecMdcKalTrack::PidType trk_type[5] = {
105 RecMdcKalTrack::electron,
106 RecMdcKalTrack::muon,
107 RecMdcKalTrack::pion,
108 RecMdcKalTrack::kaon,
109 RecMdcKalTrack::proton,
110 };
111 double mass[5] = {
112 0.000511,
113 0.105658,
114 0.13957,
115 0.493677,
116 0.938272,
117 };
118 for(unsigned int pid = 0; pid < 5; pid++)
119 {
120 mdckalTrk->setPidType(trk_type[pid]);
121 m_p[pid] = mdckalTrk->p();
122 m_betagamma[pid] = m_p[pid] / mass[pid];
123 m_charge[pid] = mdckalTrk->charge();
124 m_cost[pid] = cos(mdckalTrk->theta());
125 }
126 }
127 else
128 {
129 for(unsigned int i = 0; i < 5; i++)
130 {
131 m_p[i] = -99;
132 m_betagamma[i] = -99;
133 m_cost[i] = -99;
134 m_charge[i] = 0;
135 }
136
137 }
138
139 //dE/dx PID
140 loadDedxInfo(track);
141 if (m_dedx_corr)
142 {
143 dedxCorrection();
144 dedxSecondCorrection();
145 }
146 //TOF PID
147 loadTOFInfo(track);
148 if (m_tof_corr)
149 {
150 if (m_tof_barrel == 1)
151 {
152 tofBarrelCorrection();
153 tofBarrelSecondCorrection();
154 }
155 else if (m_tof_barrel == 0)
156 {
157 tofEndcapCorrection();
158 tofEndcapSecondCorrection();
159 }
160 }
161 //EMC
162 loadEMCInfo(track);
163
164 calprob();
165}
166
167void SimplePIDSvc::calprob()
168{
169 bool usededx = false;
170 bool usetof = false;
171
172 for (unsigned int i = 0; i < 5 ;i++)
173 {
174 if (!usededx && fabs(m_dedx_chi[i]) < m_dedx_chi_cut)
175 usededx = true;
176 if (!usetof && fabs(m_tof_chi[i]) < m_tof_chi_cut)
177 usetof = true;
178
179 m_dedx_only[i] = false;
180 }
181 if (!usededx)
182 {
183 for(unsigned int i = 0; i < 5; i++)
184 m_dedx_chi[i] = -99;
185 }
186 if (!usetof)
187 {
188 for(unsigned int i = 0; i < 5; i++)
189 m_tof_chi[i] = -99;
190 }
191
192 for (unsigned int i = 0; i < 5; i++)
193 {
194 m_prob[i] = -99;
195 double chi2 = 0;
196 int ndf = 0;
197
198 if (usededx && usetof)
199 {
200 chi2 = pow(m_dedx_chi[i], 2) + pow(m_tof_chi[i], 2);
201 ndf = 2;
202 }
203 else if (usededx && !usetof)
204 {
205 chi2 = pow(m_dedx_chi[i], 2);
206 ndf = 1;
207 m_dedx_only[i] = true;
208 }
209 else if (!usededx && usetof)
210 {
211 chi2 = pow(m_tof_chi[i],2);
212 ndf = 1;
213 }
214 if (ndf > 0 && chi2 > 0)
215 m_prob[i] = TMath::Prob(chi2, ndf);
216 }
217}
218
219int SimplePIDSvc::getRunIdx(int run_no)
220{
221 // -1: beyond correction region
222 // 0: 2010 psi(3770) data
223 // 1: 2011 psi(3770) data
224 // 2: 2010 psi(3770) mc
225 // 3: 2011 psi(3770) mc
226 const int RUN_BEGIN_DATA_10 = 11414;
227 const int RUN_END_DATA_10 = 14604;
228 const int RUN_BEGIN_MC_10 = -14604;
229 const int RUN_END_MC_10 = -11414;
230 const int RUN_BEGIN_DATA_11 = 20448;
231 const int RUN_END_DATA_11 = 23454;
232 const int RUN_BEGIN_MC_11 = -23454;
233 const int RUN_END_MC_11 = -20448;
234
235 const int RUN_BEGIN_DATA_22 = 70521;
236 const int RUN_END_DATA_22 = 73930;
237 const int RUN_BEGIN_MC_22 = -73930;
238 const int RUN_END_MC_22 = -70521;
239
240 const int RUN_BEGIN_DATA_23 = 74031;
241 const int RUN_END_DATA_23 = 78536;
242 const int RUN_BEGIN_MC_23 = -78536;
243 const int RUN_END_MC_23 = -74031;
244
245 const int RUN_BEGIN_DATA_24 = 78615;
246 const int RUN_END_DATA_24 = 81094;
247 const int RUN_BEGIN_MC_24 = -81094;
248 const int RUN_END_MC_24 = -78615;
249
250 const int RUN_BEGIN_DATA_OffRes_24 = 81095;
251 const int RUN_END_DATA_OffRes_24 = 81631;
252 const int RUN_BEGIN_MC_OffRes_24 = -81631;
253 const int RUN_END_MC_OffRes_24 = -81095;
254
255 if (run_no >= RUN_BEGIN_DATA_10 && run_no <= RUN_END_DATA_10)
256 return 0;
257 else if (run_no >= RUN_BEGIN_DATA_11 && run_no <= RUN_END_DATA_11)
258 return 1;
259 else if (run_no >= RUN_BEGIN_MC_10 && run_no <= RUN_END_MC_10)
260 return 2;
261 else if (run_no >= RUN_BEGIN_MC_11 && run_no <= RUN_END_MC_11)
262 return 3;
263
264 else if (run_no >= RUN_BEGIN_DATA_22 && run_no <= RUN_END_DATA_22)
265 return 1;
266 else if (run_no >= RUN_BEGIN_MC_22 && run_no <= RUN_END_MC_22)
267 return 3;
268 else if (run_no >= RUN_BEGIN_DATA_23 && run_no <= RUN_END_DATA_23)
269 return 1;
270 else if (run_no >= RUN_BEGIN_MC_23 && run_no <= RUN_END_MC_23)
271 return 3;
272 else if (run_no >= RUN_BEGIN_DATA_24 && run_no <= RUN_END_DATA_24)
273 return 1;
274 else if (run_no >= RUN_BEGIN_MC_24 && run_no <= RUN_END_MC_24)
275 return 3;
276 else if (run_no >= RUN_BEGIN_DATA_OffRes_24 && run_no <= RUN_END_DATA_OffRes_24)
277 return 1;
278 else if (run_no >= RUN_BEGIN_MC_OffRes_24 && run_no <= RUN_END_MC_OffRes_24)
279 return 3;
280 else
281 return -1;
282
283}
284
285void SimplePIDSvc::loadTOFInfo(EvtRecTrack *track)
286{
287 //Initialization
288 for (unsigned int i = 0; i < 8; i++)
289 {
290 for (unsigned int j = 0; j < 5; j++)
291 m_tof_dt[i][j] = -99.;
292 m_tof_ph[i] = -99.;
293 }
294 for (unsigned int i = 0; i < 2; i++)
295 {
296 m_tof_zr[i] = -9999.;
297 m_tof_counter[i] = -1;
298 }
299 for (unsigned int i = 0; i < 5; i++)
300 {
301 m_tof_chi[i] = -99.;
302 }
303 m_tof_barrel = -1;
304
305 if (!track->isExtTrackValid() || !track->isTofTrackValid()) return;
306
307 SmartRefVector<RecTofTrack> tofTrk = track->tofTrack();
308 SmartRefVector<RecTofTrack>::iterator it;
309 RecExtTrack* extTrk = track->extTrack();
310 double zrhit[2];
311 zrhit[0] = extTrk->tof1Position().z();
312 zrhit[1] = extTrk->tof2Position().z();
313
314 TofHitStatus *hitst = new TofHitStatus;
315
316 for (it = tofTrk.begin(); it != tofTrk.end(); it++)
317 {
318 unsigned int st = (*it)->status();
319 hitst->setStatus(st);
320 if (hitst->is_raw()) continue; //empty TOF hit
321 bool barrel = hitst->is_barrel();
322 bool readout = hitst->is_readout();
323 bool counter = hitst->is_counter();
324 bool cluster = hitst->is_cluster();
325 int layer = hitst->layer();
326 double tof = (*it)->tof();
327 double ph = (*it)->ph();
328 m_tof_counter[layer-1] = (*it)->tofID();
329
330 if (barrel)
331 {
332 m_tof_barrel = 1;
333 }
334 else
335 {
336 m_tof_barrel = 0;
337 zrhit[0] = extTrk->tof1Position().rho();
338 }
339 m_tof_zr[0] = zrhit[0];
340 m_tof_zr[1] = zrhit[1];
341
342 int idx = -1;
343 if (readout)
344 {
345 if (barrel)
346 idx = ((st & 0xC0) >> 5) + (((st ^ 0x20) & 0x20) >> 5) - 2;
347 else
348 idx = 7;
349 }
350 else if (counter)
351 {
352 idx = layer + 3;
353 }
354 else if (cluster)
355 {
356 idx = 6;
357 }
358 if (idx == -1) continue;
359 m_tof_ph[idx] = ph;
360 for (unsigned int i = 0; i < 5; i++)
361 {
362 double offset = (*it)->toffset(i);
363 double texp = (*it)->texp(i);
364 if (texp < 0.0) continue;
365 double dt = tof - offset - texp;
366 m_tof_dt[idx][i] = dt;
367 }
368 }
369 delete hitst;
370}
371
372void SimplePIDSvc::tofBarrelCorrection()
373{
374 const double EPS = 1e-4;
375 const double BG_LOW = 0.20;
376 const double BG_UP = 7.40;
377 const double COST_LOW = -0.81;
378 const double COST_UP = 0.81;
379 const double Q_LOW = 0.;
380 const double Q_UP = 9000.;
381 const double P_LOW = 0.2;
382 const double P_UP = 1.3;
383 const int BIN_BG = 15;
384 const int BIN_COST = 15;
385 const int BIN_P = 15;
386 double BG[BIN_BG + 1] = {0.20, 0.87, 1.11, 1.35, 1.55, 1.72, 1.91, 2.17, 2.63, 3.05, 3.47, 3.93, 4.50, 5.27, 6.00, 7.40};
387 double COST[BIN_COST + 1] = {-0.81, -0.64, -0.53, -0.43, -0.33, -0.23, -0.13, -0.04, 0.05, 0.14, 0.24, 0.34, 0.44, 0.54, 0.65, 0.81};
388 double P[BIN_P + 1] = {0.20, 0.47, 0.56, 0.65, 0.72, 0.79, 0.86, 0.92, 0.98, 1.03, 1.08, 1.13, 1.17, 1.22, 1.26, 1.30};
389 int idx = getRunIdx(m_run);
390
391 if (idx != -1)
392 {
393 for (unsigned int i = 0; i < 4; i++)
394 {// only correct e, pi, K
395 double bg;
396 int bin_bg, bin_cost, bin_wgt;
397 int pid;
398 if (i == 0)
399 {
400 bg = max(P_LOW+EPS, min(P_UP-EPS, m_p[i]));
401 bin_bg = findBin(P, BIN_P, bg);
402 pid = 0;
403 }
404 else if (i == 2 || i == 3)
405 {
406 bg = max(BG_LOW+EPS, min(BG_UP-EPS, m_betagamma[i]));
407 bin_bg = findBin(BG, BIN_BG, bg);
408 pid = 1;
409 }
410 else
411 {
412 continue;
413 }
414 double cost = m_cost[i];
415 int charge = m_charge[i];
416 double t[5], q;
417 double offset, sigma;
418 double offset_q, offset_bgcost;
419 int flag[4] = {0, 0, 0, 0, };
420 cost = max(COST_LOW+EPS, min(COST_UP-EPS, cost));
421 bin_cost = findBin(COST, BIN_COST, cost);
422 if (bin_bg == -1 || bin_cost == -1) continue;
423
424 //corrections
425 for (unsigned int j = 0; j < 4; j++)
426 {
427 t[j] = m_tof_dt[j][i];
428 if (fabs(t[j] + 99.) < EPS)//no readout
429 flag[j] = 0;
430 else
431 flag[j] = 1;
432 q = m_tof_ph[j];
433 q = max(Q_LOW+EPS, min(Q_UP-EPS, q));
434 if (charge == 1)
435 {
436 offset_q = h_tof_p_q_offset[pid][idx][j]->Interpolate( q );
437 offset_bgcost = h_tof_p_bgcost_offset[pid][idx][j]->Interpolate( bg, cost );
438 t[j] = t[j] - offset_q - offset_bgcost;
439 }
440 else
441 {
442 offset_q = h_tof_m_q_offset[pid][idx][j]->Interpolate( q );
443 offset_bgcost = h_tof_m_bgcost_offset[pid][idx][j]->Interpolate( bg, cost );
444 t[j] = t[j] - offset_q - offset_bgcost;
445 }
446 }
447 bin_wgt = flag[0]*8 + flag[1]*4 + flag[2]*2 + flag[3] - 1;
448 if (bin_wgt == -1) continue;
449 t[4] = 0;
450 for (unsigned int j = 0; j < 4; j++)
451 {
452 if (charge == 1)
453 t[4] += t[j] * h_tof_p_wgt[pid][idx][bin_wgt][j]->GetBinContent( bin_bg+1, bin_cost+1 );
454 else
455 t[4] += t[j] * h_tof_m_wgt[pid][idx][bin_wgt][j]->GetBinContent( bin_bg+1, bin_cost+1 );
456 }
457 if (charge == 1)
458 {
459 t[4] /= h_tof_p_wgt[pid][idx][bin_wgt][4]->GetBinContent( bin_bg+1, bin_cost+1 );
460 offset = h_tof_p_final_offset[pid][idx][bin_wgt]->Interpolate( bg, cost );
461 sigma = h_tof_p_final_sigma[pid][idx][bin_wgt]-> Interpolate( bg, cost );
462 m_tof_chi[i] = (t[4] - offset) / sigma;
463 }
464 else
465 {
466 t[4] /= h_tof_m_wgt[pid][idx][bin_wgt][4]->GetBinContent( bin_bg+1, bin_cost+1 );
467 offset = h_tof_m_final_offset[pid][idx][bin_wgt]->Interpolate( bg, cost );
468 sigma = h_tof_m_final_sigma[pid][idx][bin_wgt]-> Interpolate( bg, cost );
469 m_tof_chi[i] = (t[4] - offset) / sigma;
470 }
471 }
472 }
473}
474
475void SimplePIDSvc::tofEndcapCorrection()
476{
477 const double EPS = 1e-4;
478 const double BG_LOW = 0.30;
479 const double BG_UP = 7.40;
480 const double Q_LOW = 0.;
481 const double Q_UP = 6000.;
482 const double COST_EAST_LOW = 0.720;
483 const double COST_EAST_UP = 0.930;
484 const double COST_WEST_LOW = -0.930;
485 const double COST_WEST_UP = -0.720;
486 const double P_LOW = 0.2;
487 const double P_UP = 1.3;
488
489 int idx = getRunIdx(m_run);
490
491 if (idx != -1)
492 {
493 for (unsigned int i = 0; i < 4; i++)
494 {// only correct e, pi, K
495 int pid;
496 double bg;
497 if (i == 0)
498 {
499 bg = max(P_LOW+EPS, min(P_UP-EPS, m_p[i]));
500 pid = 0;
501 }
502 else if (i == 2 || i == 3)
503 {
504 bg = max(BG_LOW+EPS, min(BG_UP-EPS, m_betagamma[i]));
505 pid = 1;
506 }
507 else
508 {
509 continue;
510 }
511
512 int flag; //0:east, 1:west
513 double cost = m_cost[i];
514 int charge = m_charge[i];
515 double t = m_tof_dt[7][i];
516 double q = m_tof_ph[7];
517 double off_q, off_bg, off_cost;
518 double sg_q, sg_bg, sg_cost;
519 if (cost > 0)
520 {
521 flag = 0;
522 cost = max(COST_EAST_LOW+EPS, min(COST_EAST_UP-EPS, cost));
523 }
524 else
525 {
526 flag = 1;
527 cost = max(COST_WEST_LOW+EPS, min(COST_WEST_UP-EPS, cost));
528 }
529 q = max(Q_LOW+EPS, min(Q_UP-EPS, q));
530
531 //corrections
532 if (charge == 1)
533 {
534 off_q = h_tofec_p_q_offset[pid][idx][flag] ->Interpolate( q );
535 sg_q = h_tofec_p_q_sigma[pid][idx][flag] ->Interpolate( q );
536 off_bg = h_tofec_p_bg_offset[pid][idx][flag] ->Interpolate( bg );
537 sg_bg = h_tofec_p_bg_sigma[pid][idx][flag] ->Interpolate( bg );
538 off_cost = h_tofec_p_cost_offset[pid][idx][flag]->Interpolate( cost );
539 sg_cost = h_tofec_p_cost_sigma[pid][idx][flag] ->Interpolate( cost );
540 m_tof_chi[i] = (((t - off_q) / sg_q - off_bg) / sg_bg - off_cost) / sg_cost;
541 }
542 else
543 {
544 off_q = h_tofec_m_q_offset[pid][idx][flag] ->Interpolate( q );
545 sg_q = h_tofec_m_q_sigma[pid][idx][flag] ->Interpolate( q );
546 off_bg = h_tofec_m_bg_offset[pid][idx][flag] ->Interpolate( bg );
547 sg_bg = h_tofec_m_bg_sigma[pid][idx][flag] ->Interpolate( bg );
548 off_cost = h_tofec_m_cost_offset[pid][idx][flag]->Interpolate( cost );
549 sg_cost = h_tofec_m_cost_sigma[pid][idx][flag] ->Interpolate( cost );
550 m_tof_chi[i] = (((t - off_q) / sg_q - off_bg) / sg_bg - off_cost) / sg_cost;
551 }
552 }
553 }
554}
555
556void SimplePIDSvc::loadDedxInfo(EvtRecTrack *track)
557{
558 if (track->isMdcDedxValid())
559 {
560 RecMdcDedx* dedx_trk = track->mdcDedx();
561 for (unsigned int i = 0; i < 5; i++)
562 m_dedx_chi[i] = dedx_trk->chi(i);
563 }
564 else
565 {
566 for (unsigned int i = 0; i < 5; i++)
567 m_dedx_chi[i] = -99;
568 }
569}
570
571void SimplePIDSvc::dedxCorrection()
572{
573 int idx = getRunIdx(m_run);
574 const double EPS = 1e-4;
575 const double BG_LOW = 0.20;
576 const double BG_UP = 7.40;
577 const double COST_LOW = -0.93;
578 const double COST_UP = 0.93;
579 const double P_LOW = 0.2;
580 const double P_UP = 1.3;
581 if (idx != -1)
582 {
583 double offset, sigma;
584 for (unsigned int i = 0; i < 4; i++)
585 {// only correct e, pi, K
586 double bg;
587 int pid;
588 if (i == 0)
589 {
590 bg = max(P_LOW+EPS, min(P_UP-EPS, m_p[i]));
591 pid = 0;
592 }
593 else if (i == 2 || i == 3)
594 {
595 bg = max(BG_LOW+EPS, min(BG_UP-EPS, m_betagamma[i]));
596 pid = 1;
597 }
598 else
599 {
600 continue;
601 }
602 double cost = m_cost[i];
603 double charge = m_charge[i];
604 cost = max(COST_LOW+EPS, min(COST_UP-EPS, cost));
605 if (charge == 1)
606 {
607 offset = h_dedx_p_offset[pid][idx]->Interpolate( bg, cost );
608 sigma = h_dedx_p_sigma[pid][idx] ->Interpolate( bg, cost );
609 m_dedx_chi[i] = (m_dedx_chi[i] - offset) / sigma;
610 }
611 else
612 {
613 offset = h_dedx_m_offset[pid][idx]->Interpolate( bg, cost );
614 sigma = h_dedx_m_sigma[pid][idx] ->Interpolate( bg, cost );
615 m_dedx_chi[i] = (m_dedx_chi[i] - offset) / sigma;
616 }
617 }
618 }
619}
620
621void SimplePIDSvc::tofBarrelSecondCorrection()
622{
623
624 int idx = getRunIdx(m_run);
625 const double EPS = 1e-4;
626 const double P_LOW = 0.0;
627 const double P_UP = 1.3;
628 const int BIN_P = 10;
629
630 const int RUN_BEGIN_DATA_10 = 11414;
631 const int RUN_END_DATA_10 = 14604;
632 const int RUN_BEGIN_MC_10 = -14604;
633 const int RUN_END_MC_10 = -11414;
634 const int RUN_BEGIN_DATA_11 = 20448;
635 const int RUN_END_DATA_11 = 23454;
636 const int RUN_BEGIN_MC_11 = -23454;
637 const int RUN_END_MC_11 = -20448;
638
639 const int RUN_BEGIN_DATA_22 = 70521;
640 const int RUN_END_DATA_22 = 73930;
641 const int RUN_BEGIN_MC_22 = -73930;
642 const int RUN_END_MC_22 = -70521;
643
644 const int RUN_BEGIN_DATA_23 = 74031;
645 const int RUN_END_DATA_23 = 78536;
646 const int RUN_BEGIN_MC_23 = -78536;
647 const int RUN_END_MC_23 = -74031;
648
649 const int RUN_BEGIN_DATA_24 = 78615;
650 const int RUN_END_DATA_24 = 81094;
651 const int RUN_BEGIN_MC_24 = -81094;
652 const int RUN_END_MC_24 = -78615;
653
654 const int RUN_BEGIN_DATA_OffRes_24 = 81095;
655 const int RUN_END_DATA_OffRes_24 = 81631;
656 const int RUN_BEGIN_MC_OffRes_24 = -81631;
657 const int RUN_END_MC_OffRes_24 = -81095;
658
659 double P[BIN_P + 1] = {0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.3};
660 if (idx != -1)
661 {
662
663 for (unsigned int i = 2; i < 4; i++){// second correct for tof of pi k
664
665 int aa=99,bb=99;
666
667 int bin_p;
668 double ptk = m_p[i];
669 ptk = max(P_LOW+EPS, min(P_UP-EPS, m_p[i]));
670 bin_p = findBin(P, BIN_P, ptk);
671
672 if(i==2){// pi
673
674 if((m_run>=RUN_BEGIN_DATA_10 && m_run<=RUN_END_DATA_10 && idx == 0)||(m_run>=RUN_BEGIN_DATA_11 && m_run<=RUN_END_DATA_11 && idx == 1))
675 {aa=2; bb=2;}// round0304 data
676
677 if((m_run>=RUN_BEGIN_MC_10 && m_run<=RUN_END_MC_10 && idx == 2)||(m_run>=RUN_BEGIN_MC_11 && m_run<=RUN_END_MC_11 && idx == 3))
678 {aa=2; bb=3;}// round0304 inc
679
680 if(m_run>=RUN_BEGIN_DATA_22 && m_run<=RUN_END_DATA_22 && idx == 1)
681 {aa=3; bb=2;}// round15 data
682
683 if(m_run>=RUN_BEGIN_MC_22 && m_run<=RUN_END_MC_22 && idx == 3)
684 {aa=3; bb=3;}// round15 inc
685
686 if(m_run>=RUN_BEGIN_DATA_23 && m_run<=RUN_END_DATA_23 && idx == 1)
687 {aa=7; bb=2;}// round16 data added by Yijia Zeng
688
689 if(m_run>=RUN_BEGIN_MC_23 && m_run<=RUN_END_MC_23 && idx == 3)
690 {aa=7; bb=3;}// round16 inc added by Yijia Zeng
691
692 if(m_run>=RUN_BEGIN_DATA_24 && m_run<=RUN_END_DATA_24 && idx == 1)
693 {aa=10; bb=2;}// round17 data added by Yijia Zeng
694
695 if(m_run>=RUN_BEGIN_MC_24 && m_run<=RUN_END_MC_24 && idx == 3)
696 {aa=10; bb=3;}// round17 inc added by Yijia Zeng
697
698 if(m_run>=RUN_BEGIN_DATA_OffRes_24 && m_run<=RUN_END_DATA_OffRes_24 && idx == 1)
699 {aa=10; bb=2;}// psipp off-resonance data added by Yijia Zeng
700
701 if(m_run>=RUN_BEGIN_MC_OffRes_24 && m_run<=RUN_END_MC_OffRes_24 && idx == 3)
702 {aa=10; bb=3;}// psipp off-resonance inc added by Yijia Zeng
703
704 }
705
706 if(i==3){// k
707
708 if((m_run>=RUN_BEGIN_DATA_10 && m_run<=RUN_END_DATA_10 && idx == 0)||(m_run>=RUN_BEGIN_DATA_11 && m_run<=RUN_END_DATA_11 && idx == 1))
709 {aa=2; bb=0;}
710
711 if((m_run>=RUN_BEGIN_MC_10 && m_run<=RUN_END_MC_10 && idx == 2)||(m_run>=RUN_BEGIN_MC_11 && m_run<=RUN_END_MC_11 && idx == 3))
712 {aa=2; bb=1;}
713
714 if(m_run>=RUN_BEGIN_DATA_22 && m_run<=RUN_END_DATA_22 && idx == 1)
715 {aa=3; bb=0;}
716
717 if(m_run>=RUN_BEGIN_MC_22 && m_run<=RUN_END_MC_22 && idx == 3)
718 {aa=3; bb=1;}
719
720 if(m_run>=RUN_BEGIN_DATA_23 && m_run<=RUN_END_DATA_23 && idx == 1)
721 {aa=7; bb=0;}// round16 data added by Yijia Zeng
722
723 if(m_run>=RUN_BEGIN_MC_23 && m_run<=RUN_END_MC_23 && idx == 3)
724 {aa=7; bb=1;}// round16 inc added by Yijia Zeng
725
726 if(m_run>=RUN_BEGIN_DATA_24 && m_run<=RUN_END_DATA_24 && idx == 1)
727 {aa=10; bb=0;}// round17 data added by Yijia Zeng
728
729 if(m_run>=RUN_BEGIN_MC_24 && m_run<=RUN_END_MC_24 && idx == 3)
730 {aa=10; bb=1;}// round17 inc added by Yijia Zeng
731
732 if(m_run>=RUN_BEGIN_DATA_OffRes_24 && m_run<=RUN_END_DATA_OffRes_24 && idx == 1)
733 {aa=10; bb=0;}// psipp off-resonance data added by Yijia Zeng
734
735 if(m_run>=RUN_BEGIN_MC_OffRes_24 && m_run<=RUN_END_MC_OffRes_24 && idx == 3)
736 {aa=10; bb=1;}// psipp off-resonance inc added by Yijia Zeng
737
738 }
739
740 if(m_tof_chi[i]!=-99 && m_run>0){
741
742 //NO Correction on the higher FIVE momentum region from Yijia for round 16 & 17 data sample
743 //NO Correction on the higher SIX momentum region from Kaikai He for round03-04 and round15
744 //The Corrected Distribution is somehow worse than the Original one
745 if(!((m_run>=RUN_BEGIN_DATA_23 && m_run<=RUN_END_DATA_23)||(m_run>=RUN_BEGIN_DATA_24 && m_run<=RUN_END_DATA_24)||(m_run>=RUN_BEGIN_DATA_OffRes_24 && m_run<=RUN_END_DATA_OffRes_24))){
746 if(bin_p<4){
747 m_tof_chi[i] = (m_tof_chi[i] - m_gaussion_mean[aa][bb][bin_p]) / (m_gaussion_sigmab[aa][bb][bin_p] / m_gaussion_sigmab[aa][bb+1][bin_p]);
748 }
749 if(bin_p>=4){
750 m_tof_chi[i] = (m_tof_chi[i] - m_gaussion_mean[aa][bb][bin_p]);
751 }
752 }
753 else{
754 if(bin_p<5){
755 m_tof_chi[i] = (m_tof_chi[i] - m_gaussion_mean[aa][bb][bin_p]) / (m_gaussion_sigmab[aa][bb][bin_p] / m_gaussion_sigmab[aa][bb+1][bin_p]);
756 }
757 if(bin_p>=5){
758 m_tof_chi[i] = (m_tof_chi[i] - m_gaussion_mean[aa][bb][bin_p]);
759 }
760 }
761
762 }
763 if(m_tof_chi[i]!=-99 && m_run<0){
764 m_tof_chi[i] = (m_tof_chi[i] - m_gaussion_mean[aa][bb][bin_p]) / (m_gaussion_sigmab[aa][bb][bin_p] / m_gaussion_sigmab[aa][bb][bin_p]);
765 }
766
767 }
768
769 }//end idx!= -1
770
771}
772
773void SimplePIDSvc::tofEndcapSecondCorrection()
774{
775
776 int idx = getRunIdx(m_run);
777 const double EPS = 1e-4;
778 const double P_LOW = 0.0;
779 const double P_UP = 1.3;
780 const int BIN_P = 10;
781
782 const int RUN_BEGIN_DATA_10 = 11414;
783 const int RUN_END_DATA_10 = 14604;
784 const int RUN_BEGIN_MC_10 = -14604;
785 const int RUN_END_MC_10 = -11414;
786 const int RUN_BEGIN_DATA_11 = 20448;
787 const int RUN_END_DATA_11 = 23454;
788 const int RUN_BEGIN_MC_11 = -23454;
789 const int RUN_END_MC_11 = -20448;
790
791 const int RUN_BEGIN_DATA_22 = 70521;
792 const int RUN_END_DATA_22 = 73930;
793 const int RUN_BEGIN_MC_22 = -73930;
794 const int RUN_END_MC_22 = -70521;
795
796 const int RUN_BEGIN_DATA_23 = 74031;
797 const int RUN_END_DATA_23 = 78536;
798 const int RUN_BEGIN_MC_23 = -78536;
799 const int RUN_END_MC_23 = -74031;
800
801 const int RUN_BEGIN_DATA_24 = 78615;
802 const int RUN_END_DATA_24 = 81094;
803 const int RUN_BEGIN_MC_24 = -81094;
804 const int RUN_END_MC_24 = -78615;
805
806 const int RUN_BEGIN_DATA_OffRes_24 = 81095;
807 const int RUN_END_DATA_OffRes_24 = 81631;
808 const int RUN_BEGIN_MC_OffRes_24 = -81631;
809 const int RUN_END_MC_OffRes_24 = -81095;
810
811 double P[BIN_P + 1] = {0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.3};
812 if (idx != -1)
813 {
814
815 for (unsigned int i = 2; i < 4; i++){// second correct for tof of pi k
816
817 int aa=99,bb=99;
818
819 int bin_p;
820 double ptk = m_p[i];
821 ptk = max(P_LOW+EPS, min(P_UP-EPS, m_p[i]));
822 bin_p = findBin(P, BIN_P, ptk);
823
824 if(i==2){// pi
825
826 if((m_run>=RUN_BEGIN_DATA_10 && m_run<=RUN_END_DATA_10 && idx == 0)||(m_run>=RUN_BEGIN_DATA_11 && m_run<=RUN_END_DATA_11 && idx == 1))
827 {aa=4; bb=2;}// round0304 data endcap
828
829 if((m_run>=RUN_BEGIN_MC_10 && m_run<=RUN_END_MC_10 && idx == 2)||(m_run>=RUN_BEGIN_MC_11 && m_run<=RUN_END_MC_11 && idx == 3))
830 {aa=4; bb=3;}// round0304 inc endcap
831
832 if(m_run>=RUN_BEGIN_DATA_22 && m_run<=RUN_END_DATA_22 && idx == 1)
833 {aa=5; bb=2;}// round15 data endcap
834
835 if(m_run>=RUN_BEGIN_MC_22 && m_run<=RUN_END_MC_22 && idx == 3)
836 {aa=5; bb=3;}// round15 inc endcap
837
838 if(m_run>=RUN_BEGIN_DATA_23 && m_run<=RUN_END_DATA_23 && idx == 1)
839 {aa=8; bb=2;}// round16 data endcap added by Yijia Zeng
840
841 if(m_run>=RUN_BEGIN_MC_23 && m_run<=RUN_END_MC_23 && idx == 3)
842 {aa=8; bb=3;}// round16 inc endcap added by Yijia Zeng
843
844 if(m_run>=RUN_BEGIN_DATA_24 && m_run<=RUN_END_DATA_24 && idx == 1)
845 {aa=11; bb=2;}// round17 data endcap added by Yijia Zeng
846
847 if(m_run>=RUN_BEGIN_MC_24 && m_run<=RUN_END_MC_24 && idx == 3)
848 {aa=11; bb=3;}// round17 inc endcap added by Yijia Zeng
849
850 if(m_run>=RUN_BEGIN_DATA_OffRes_24 && m_run<=RUN_END_DATA_OffRes_24 && idx == 1)
851 {aa=11; bb=2;}// psipp off-resonance data endcap added by Yijia Zeng
852
853 if(m_run>=RUN_BEGIN_MC_OffRes_24 && m_run<=RUN_END_MC_OffRes_24 && idx == 3)
854 {aa=11; bb=3;}// psipp off-resonance inc endcap added by Yijia Zeng
855
856 }
857
858 if(i==3){// k
859
860 if((m_run>=RUN_BEGIN_DATA_10 && m_run<=RUN_END_DATA_10 && idx == 0)||(m_run>=RUN_BEGIN_DATA_11 && m_run<=RUN_END_DATA_11 && idx == 1))
861 {aa=4; bb=0;}
862
863 if((m_run>=RUN_BEGIN_MC_10 && m_run<=RUN_END_MC_10 && idx == 2)||(m_run>=RUN_BEGIN_MC_11 && m_run<=RUN_END_MC_11 && idx == 3))
864 {aa=4; bb=1;}
865
866 if(m_run>=RUN_BEGIN_DATA_22 && m_run<=RUN_END_DATA_22 && idx == 1)
867 {aa=5; bb=0;}
868
869 if(m_run>=RUN_BEGIN_MC_22 && m_run<=RUN_END_MC_22 && idx == 3)
870 {aa=5; bb=1;}
871
872 if(m_run>=RUN_BEGIN_DATA_23 && m_run<=RUN_END_DATA_23 && idx == 1)
873 {aa=8; bb=0;}// round16 data endcap added by Yijia Zeng
874
875 if(m_run>=RUN_BEGIN_MC_23 && m_run<=RUN_END_MC_23 && idx == 3)
876 {aa=8; bb=1;}// round16 inc endcap added by Yijia Zeng
877
878 if(m_run>=RUN_BEGIN_DATA_24 && m_run<=RUN_END_DATA_24 && idx == 1)
879 {aa=11; bb=0;}// round17 data endcap added by Yijia Zeng
880
881 if(m_run>=RUN_BEGIN_MC_24 && m_run<=RUN_END_MC_24 && idx == 3)
882 {aa=11; bb=1;}// round17 inc endcap added by Yijia Zeng
883
884 if(m_run>=RUN_BEGIN_DATA_OffRes_24 && m_run<=RUN_END_DATA_OffRes_24 && idx == 1)
885 {aa=11; bb=0;}// psipp off-resonance data endcap added by Yijia Zeng
886
887 if(m_run>=RUN_BEGIN_MC_OffRes_24 && m_run<=RUN_END_MC_OffRes_24 && idx == 3)
888 {aa=11; bb=1;}// psipp off-resonance inc endcap added by Yijia Zeng
889
890 }
891
892 if(m_tof_chi[i]!=-99 && ( aa==5 || aa==8 || aa==11)){// for round15 and round16 and round17 data and inc
893
894 m_tof_chi[i] = (m_tof_chi[i] - m_gaussion_mean[aa][bb][bin_p]) / (60./110.);
895
896 }
897
898 if(m_tof_chi[i]!=-99 && aa==4){// for round0304 data and inc
899
900 m_tof_chi[i] = (m_tof_chi[i] - m_gaussion_mean[aa][bb][bin_p]);
901
902 }
903
904
905 }
906
907 }//end idx!= -1
908
909
910}
911
912void SimplePIDSvc::dedxSecondCorrection()
913
914{
915
916 int idx = getRunIdx(m_run);
917 const double EPS = 1e-4;
918 const double P_LOW = 0.0;
919 const double P_UP = 1.3;
920 const int BIN_P = 10;
921
922 const int RUN_BEGIN_DATA_10 = 11414;
923 const int RUN_END_DATA_10 = 14604;
924 const int RUN_BEGIN_MC_10 = -14604;
925 const int RUN_END_MC_10 = -11414;
926 const int RUN_BEGIN_DATA_11 = 20448;
927 const int RUN_END_DATA_11 = 23454;
928 const int RUN_BEGIN_MC_11 = -23454;
929 const int RUN_END_MC_11 = -20448;
930
931 const int RUN_BEGIN_DATA_22 = 70521;
932 const int RUN_END_DATA_22 = 73930;
933 const int RUN_BEGIN_MC_22 = -73930;
934 const int RUN_END_MC_22 = -70521;
935
936 const int RUN_BEGIN_DATA_23 = 74031;
937 const int RUN_END_DATA_23 = 78536;
938 const int RUN_BEGIN_MC_23 = -78536;
939 const int RUN_END_MC_23 = -74031;
940
941 const int RUN_BEGIN_DATA_24 = 78615;
942 const int RUN_END_DATA_24 = 81094;
943 const int RUN_BEGIN_MC_24 = -81094;
944 const int RUN_END_MC_24 = -78615;
945
946 const int RUN_BEGIN_DATA_OffRes_24 = 81095;
947 const int RUN_END_DATA_OffRes_24 = 81631;
948 const int RUN_BEGIN_MC_OffRes_24 = -81631;
949 const int RUN_END_MC_OffRes_24 = -81095;
950
951 double P[BIN_P + 1] = {0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.3};
952 if (idx != -1)
953 {
954
955 for (unsigned int i = 2; i < 4; i++){// second correct for dedx of pi k
956
957 int aa=99,bb=99;
958
959 int bin_p;
960 double ptk = m_p[i];
961 ptk = max(P_LOW+EPS, min(P_UP-EPS, m_p[i]));
962 bin_p = findBin(P, BIN_P, ptk);
963
964 if(i==2){// pi
965
966 if((m_run>=RUN_BEGIN_DATA_10 && m_run<=RUN_END_DATA_10 && idx == 0)||(m_run>=RUN_BEGIN_DATA_11 && m_run<=RUN_END_DATA_11 && idx == 1))
967 {aa=0; bb=2;}// round0304 data
968
969 if((m_run>=RUN_BEGIN_MC_10 && m_run<=RUN_END_MC_10 && idx == 2)||(m_run>=RUN_BEGIN_MC_11 && m_run<=RUN_END_MC_11 && idx == 3))
970 {aa=0; bb=3;}// round0304 inc
971
972 if(m_run>=RUN_BEGIN_DATA_22 && m_run<=RUN_END_DATA_22 && idx == 1)
973 {aa=1; bb=2;}// round15 data
974
975 if(m_run>=RUN_BEGIN_MC_22 && m_run<=RUN_END_MC_22 && idx == 3)
976 {aa=1; bb=3;}// round15 inc
977
978 if(m_run>=RUN_BEGIN_DATA_23 && m_run<=RUN_END_DATA_23 && idx == 1)
979 {aa=6; bb=2;}// round16 data added by Yijia Zeng
980
981 if(m_run>=RUN_BEGIN_MC_23 && m_run<=RUN_END_MC_23 && idx == 3)
982 {aa=6; bb=3;}// round16 inc added by Yijia Zeng
983
984 if(m_run>=RUN_BEGIN_DATA_24 && m_run<=RUN_END_DATA_24 && idx == 1)
985 {aa=9; bb=2;}// round17 data added by Yijia Zeng
986
987 if(m_run>=RUN_BEGIN_MC_24 && m_run<=RUN_END_MC_24 && idx == 3)
988 {aa=9; bb=3;}// round17 inc added by Yijia Zeng
989
990 if(m_run>=RUN_BEGIN_DATA_OffRes_24 && m_run<=RUN_END_DATA_OffRes_24 && idx == 1)
991 {aa=9; bb=2;}// psipp off-resonance data added by Yijia Zeng
992
993 if(m_run>=RUN_BEGIN_MC_OffRes_24 && m_run<=RUN_END_MC_OffRes_24 && idx == 3)
994 {aa=9; bb=3;}// psipp off-resonance inc added by Yijia Zeng
995
996 }
997
998 if(i==3){// k
999
1000 if((m_run>=RUN_BEGIN_DATA_10 && m_run<=RUN_END_DATA_10 && idx == 0)||(m_run>=RUN_BEGIN_DATA_11 && m_run<=RUN_END_DATA_11 && idx == 1))
1001 {aa=0; bb=0;}
1002
1003 if((m_run>=RUN_BEGIN_MC_10 && m_run<=RUN_END_MC_10 && idx == 2)||(m_run>=RUN_BEGIN_MC_11 && m_run<=RUN_END_MC_11 && idx == 3))
1004 {aa=0; bb=1;}
1005
1006 if(m_run>=RUN_BEGIN_DATA_22 && m_run<=RUN_END_DATA_22 && idx == 1)
1007 {aa=1; bb=0;}
1008
1009 if(m_run>=RUN_BEGIN_MC_22 && m_run<=RUN_END_MC_22 && idx == 3)
1010 {aa=1; bb=1;}
1011
1012 if(m_run>=RUN_BEGIN_DATA_23 && m_run<=RUN_END_DATA_23 && idx == 1)
1013 {aa=6; bb=0;}// round16 data added by Yijia Zeng
1014
1015 if(m_run>=RUN_BEGIN_MC_23 && m_run<=RUN_END_MC_23 && idx == 3)
1016 {aa=6; bb=1;}// round16 inc added by Yijia Zeng
1017
1018 if(m_run>=RUN_BEGIN_DATA_24 && m_run<=RUN_END_DATA_24 && idx == 1)
1019 {aa=9; bb=0;}// round17 data added by Yijia Zeng
1020
1021 if(m_run>=RUN_BEGIN_MC_24 && m_run<=RUN_END_MC_24 && idx == 3)
1022 {aa=9; bb=1;}// round17 inc added by Yijia Zeng
1023
1024 if(m_run>=RUN_BEGIN_DATA_OffRes_24 && m_run<=RUN_END_DATA_OffRes_24 && idx == 1)
1025 {aa=9; bb=0;}// psipp off-resonance data added by Yijia Zeng
1026
1027 if(m_run>=RUN_BEGIN_MC_OffRes_24 && m_run<=RUN_END_MC_OffRes_24 && idx == 3)
1028 {aa=9; bb=1;}// psipp off-resonance inc added by Yijia Zeng
1029
1030 }
1031
1032 if(m_dedx_chi[i]!=-99 && m_run>0){
1033 m_dedx_chi[i] = (m_dedx_chi[i] - m_gaussion_mean[aa][bb][bin_p]) / (m_gaussion_sigmab[aa][bb][bin_p] / m_gaussion_sigmab[aa][bb+1][bin_p]);
1034 }
1035 if(m_dedx_chi[i]!=-99 && m_run<0){
1036 m_dedx_chi[i] = (m_dedx_chi[i] - m_gaussion_mean[aa][bb][bin_p]) / (m_gaussion_sigmab[aa][bb][bin_p] / m_gaussion_sigmab[aa][bb][bin_p]);
1037 }
1038
1039 }
1040
1041 }//end idx!= -1
1042
1043}
1044
1045void SimplePIDSvc::loadSecondPar()
1046{
1047 string path = getenv("SIMPLEPIDSVCROOT");
1048
1049 for(int i=0; i<12; i++){
1050
1051 const char *dir;
1052
1053 if(i == 0) dir = "round0304_dedx";
1054 else if(i == 1) dir = "round15_dedx";
1055 else if(i == 2) dir = "round0304_tof";
1056 else if(i == 3) dir = "round15_tof";
1057 else if(i == 4) dir = "round0304_tof_endcap";
1058 else if(i == 5) dir = "round15_tof_endcap";
1059 //Added by Yijia Zeng for round 16 data sample
1060 else if(i == 6) dir = "round16_dedx";
1061 else if(i == 7) dir = "round16_tof";
1062 else if(i == 8) dir = "round16_tof_endcap";
1063 //Added by Yijia Zeng for round 17 data sample
1064 else if(i == 9) dir = "round17_dedx";
1065 else if(i == 10) dir = "round17_tof";
1066 else if(i == 11) dir = "round17_tof_endcap";
1067
1068 else{
1069 cout << "Boundary Error! " << endl;
1070 exit(1);
1071 }
1072
1073 for(int j=0; j<4; j++){
1074
1075 const char *name;
1076
1077 if(j == 0) name = "data_K";
1078 else if(j == 1) name = "inc_K";
1079 else if(j == 2) name = "data_pi";
1080 else if(j == 3) name = "inc_pi";
1081 else{
1082 cout << "Boundary Error! " << endl;
1083 exit(1);
1084 }
1085
1086 ifstream second_cor( path + Form("/share/second_correct/%s/%s_%s.dat",dir,dir,name));
1087
1088 for(int m=0; m<10; m++){
1089
1090 second_cor>>m_gaussion_mean[i][j][m]>>m_gaussion_sigma[i][j][m]>>m_gaussion_sigmab[i][j][m];
1091
1092 }
1093
1094 second_cor.close();
1095
1096 }
1097 }
1098
1099
1100}
1101
1102void SimplePIDSvc::loadHistogram()
1103{
1104 string path = getenv("SIMPLEPIDSVCROOT");
1105 vector<string> filename;
1106 for (unsigned int idx = 0; idx < 2; idx++)
1107 {
1108 const char *dir;
1109 if (idx == 0)
1110 dir = "electron";
1111 else if (idx == 1)
1112 dir = "kpi";
1113 else
1114 {
1115 cout << "Boundary Error! " << endl;
1116 exit(1);
1117 }
1118
1119 //dedx
1120 filename.clear();
1121 filename.push_back( path + Form("/share/%s/dedx/dedx_d10.root", dir) );
1122 filename.push_back( path + Form("/share/%s/dedx/dedx_d11.root", dir) );
1123 filename.push_back( path + Form("/share/%s/dedx/dedx_m10.root", dir) );
1124 filename.push_back( path + Form("/share/%s/dedx/dedx_m11.root", dir) );
1125 for (unsigned int i = 0; i < filename.size(); i++)
1126 {
1127 f_dedx[idx][i] = new TFile(filename[i].c_str(), "READ");
1128 const char *name;
1129 if (i == 0)
1130 name = "d10";
1131 else if (i == 1)
1132 name = "d11";
1133 else if (i == 2)
1134 name = "m10";
1135 else if (i == 3)
1136 name = "m11";
1137 else
1138 {
1139 cout << "Boundary Error! " << endl;
1140 exit(1);
1141 }
1142 h_dedx_p_offset[idx][i] = (TH2D*)f_dedx[idx][i]->Get( Form("h_dedx_p_offset_%s", name) );
1143 h_dedx_p_sigma[idx][i] = (TH2D*)f_dedx[idx][i]->Get( Form("h_dedx_p_sigma_%s" , name) );
1144 h_dedx_m_offset[idx][i] = (TH2D*)f_dedx[idx][i]->Get( Form("h_dedx_m_offset_%s", name) );
1145 h_dedx_m_sigma[idx][i] = (TH2D*)f_dedx[idx][i]->Get( Form("h_dedx_m_sigma_%s" , name) );
1146 }
1147 //tof_barrel q
1148 filename.clear();
1149 filename.push_back( path + Form("/share/%s/tof_barrel/tof_q_d10.root", dir) );
1150 filename.push_back( path + Form("/share/%s/tof_barrel/tof_q_d11.root", dir) );
1151 filename.push_back( path + Form("/share/%s/tof_barrel/tof_q_m10.root", dir) );
1152 filename.push_back( path + Form("/share/%s/tof_barrel/tof_q_m11.root", dir) );
1153 for (unsigned int i = 0; i < filename.size(); i++)
1154 {
1155 f_tof_q[idx][i] = new TFile(filename[i].c_str(), "READ");
1156 const char *name;
1157 if (i == 0)
1158 name = "d10";
1159 else if (i == 1)
1160 name = "d11";
1161 else if (i == 2)
1162 name = "m10";
1163 else if (i == 3)
1164 name = "m11";
1165 else
1166 {
1167 cout << "Boundary Error! " << endl;
1168 exit(1);
1169 }
1170 for (unsigned int j = 0; j < 4; j++)
1171 {
1172 h_tof_p_q_offset[idx][i][j] = (TH1D*)f_tof_q[idx][i]->Get( Form("h_tof_p_q_offset_%s_%d", name, j) );
1173 h_tof_m_q_offset[idx][i][j] = (TH1D*)f_tof_q[idx][i]->Get( Form("h_tof_m_q_offset_%s_%d", name, j) );
1174 h_tof_p_q_sigma[idx][i][j] = (TH1D*)f_tof_q[idx][i]->Get( Form("h_tof_p_q_sigma_%s_%d" , name, j) );
1175 h_tof_m_q_sigma[idx][i][j] = (TH1D*)f_tof_q[idx][i]->Get( Form("h_tof_m_q_sigma_%s_%d" , name, j) );
1176 }
1177 }
1178 //tof_barrel bg&cost
1179 filename.clear();
1180 filename.push_back( path + Form("/share/%s/tof_barrel/tof_bg_cost_d10.root", dir) );
1181 filename.push_back( path + Form("/share/%s/tof_barrel/tof_bg_cost_d11.root", dir) );
1182 filename.push_back( path + Form("/share/%s/tof_barrel/tof_bg_cost_m10.root", dir) );
1183 filename.push_back( path + Form("/share/%s/tof_barrel/tof_bg_cost_m11.root", dir) );
1184 for (unsigned int i = 0; i < filename.size(); i++)
1185 {
1186 f_tof_bgcost[idx][i] = new TFile(filename[i].c_str(), "READ");
1187 const char *name;
1188 if (i == 0)
1189 name = "d10";
1190 else if (i == 1)
1191 name = "d11";
1192 else if (i == 2)
1193 name = "m10";
1194 else if (i == 3)
1195 name = "m11";
1196 else
1197 {
1198 cout << "Boundary Error! " << endl;
1199 exit(1);
1200 }
1201 for (unsigned int j = 0; j < 4; j++)
1202 {
1203 h_tof_p_bgcost_offset[idx][i][j] = (TH2D*)f_tof_bgcost[idx][i]->Get( Form("h_tof_p_bgcost_offset_%s_%d", name, j) );
1204 h_tof_m_bgcost_offset[idx][i][j] = (TH2D*)f_tof_bgcost[idx][i]->Get( Form("h_tof_m_bgcost_offset_%s_%d", name, j) );
1205 h_tof_p_bgcost_sigma[idx][i][j] = (TH2D*)f_tof_bgcost[idx][i]->Get( Form("h_tof_p_bgcost_sigma_%s_%d" , name, j) );
1206 h_tof_m_bgcost_sigma[idx][i][j] = (TH2D*)f_tof_bgcost[idx][i]->Get( Form("h_tof_m_bgcost_sigma_%s_%d" , name, j) );
1207 }
1208 }
1209 //tof_barrel wgt
1210 filename.clear();
1211 filename.push_back( path + Form("/share/%s/tof_barrel/tof_wgt_d10.root", dir) );
1212 filename.push_back( path + Form("/share/%s/tof_barrel/tof_wgt_d11.root", dir) );
1213 filename.push_back( path + Form("/share/%s/tof_barrel/tof_wgt_m10.root", dir) );
1214 filename.push_back( path + Form("/share/%s/tof_barrel/tof_wgt_m11.root", dir) );
1215 for (unsigned int i = 0; i < filename.size(); i++)
1216 {
1217 f_tof_wgt[idx][i] = new TFile(filename[i].c_str(), "READ");
1218 const char *name;
1219 if (i == 0)
1220 name = "d10";
1221 else if (i == 1)
1222 name = "d11";
1223 else if (i == 2)
1224 name = "m10";
1225 else if (i == 3)
1226 name = "m11";
1227 else
1228 {
1229 cout << "Boundary Error! " << endl;
1230 exit(1);
1231 }
1232 for (unsigned int j = 0; j < 15; j++)
1233 {
1234 for (unsigned int k = 0; k < 5; k++)
1235 {
1236 h_tof_p_wgt[idx][i][j][k] = (TH2D*)f_tof_wgt[idx][i]->Get( Form("h_tof_p_wgt_%s_%d_%d", name, j, k) );
1237 h_tof_m_wgt[idx][i][j][k] = (TH2D*)f_tof_wgt[idx][i]->Get( Form("h_tof_m_wgt_%s_%d_%d", name, j, k) );
1238 }
1239 }
1240 }
1241 //tof_barrel corr
1242 filename.clear();
1243 filename.push_back( path + Form("/share/%s/tof_barrel/tof_final_d10.root", dir) );
1244 filename.push_back( path + Form("/share/%s/tof_barrel/tof_final_d11.root", dir) );
1245 filename.push_back( path + Form("/share/%s/tof_barrel/tof_final_m10.root", dir) );
1246 filename.push_back( path + Form("/share/%s/tof_barrel/tof_final_m11.root", dir) );
1247 for (unsigned int i = 0; i < filename.size(); i++)
1248 {
1249 f_tof_final[idx][i] = new TFile(filename[i].c_str(), "READ");
1250 const char *name;
1251 if (i == 0)
1252 name = "d10";
1253 else if (i == 1)
1254 name = "d11";
1255 else if (i == 2)
1256 name = "m10";
1257 else if (i == 3)
1258 name = "m11";
1259 else
1260 {
1261 cout << "Boundary Error! " << endl;
1262 exit(1);
1263 }
1264 for (unsigned int j = 0; j < 15; j++)
1265 {
1266 h_tof_p_final_offset[idx][i][j] = (TH2D*)f_tof_final[idx][i]->Get( Form("h_tof_p_final_offset_%s_%d", name, j) );
1267 h_tof_m_final_offset[idx][i][j] = (TH2D*)f_tof_final[idx][i]->Get( Form("h_tof_m_final_offset_%s_%d", name, j) );
1268 h_tof_p_final_sigma[idx][i][j] = (TH2D*)f_tof_final[idx][i]->Get( Form("h_tof_p_final_sigma_%s_%d" , name, j) );
1269 h_tof_m_final_sigma[idx][i][j] = (TH2D*)f_tof_final[idx][i]->Get( Form("h_tof_m_final_sigma_%s_%d" , name, j) );
1270 }
1271 }
1272 //tof_endcap q
1273 filename.clear();
1274 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_q_d10.root", dir) );
1275 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_q_d11.root", dir) );
1276 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_q_m10.root", dir) );
1277 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_q_m11.root", dir) );
1278 for (unsigned int i = 0; i < filename.size(); i++)
1279 {
1280 f_tofec_q[idx][i] = new TFile(filename[i].c_str(), "READ");
1281 const char *name;
1282 if (i == 0)
1283 name = "d10";
1284 else if (i == 1)
1285 name = "d11";
1286 else if (i == 2)
1287 name = "m10";
1288 else if (i == 3)
1289 name = "m11";
1290 else
1291 {
1292 cout << "Boundary Error! " << endl;
1293 exit(1);
1294 }
1295 for (unsigned int j = 0; j < 2; j++)
1296 {
1297 h_tofec_p_q_offset[idx][i][j] = (TH1D*)f_tofec_q[idx][i]->Get( Form("h_tofec_p_q_offset_%s_%d", name, j) );
1298 h_tofec_m_q_offset[idx][i][j] = (TH1D*)f_tofec_q[idx][i]->Get( Form("h_tofec_m_q_offset_%s_%d", name, j) );
1299 h_tofec_p_q_sigma[idx][i][j] = (TH1D*)f_tofec_q[idx][i]->Get( Form("h_tofec_p_q_sigma_%s_%d" , name, j) );
1300 h_tofec_m_q_sigma[idx][i][j] = (TH1D*)f_tofec_q[idx][i]->Get( Form("h_tofec_m_q_sigma_%s_%d" , name, j) );
1301 }
1302 }
1303 //tof_endcap bg
1304 filename.clear();
1305 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_bg_d10.root", dir) );
1306 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_bg_d11.root", dir) );
1307 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_bg_m10.root", dir) );
1308 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_bg_m11.root", dir) );
1309 for (unsigned int i = 0; i < filename.size(); i++)
1310 {
1311 f_tofec_bg[idx][i] = new TFile(filename[i].c_str(), "READ");
1312 const char *name;
1313 if (i == 0)
1314 name = "d10";
1315 else if (i == 1)
1316 name = "d11";
1317 else if (i == 2)
1318 name = "m10";
1319 else if (i == 3)
1320 name = "m11";
1321 else
1322 {
1323 cout << "Boundary Error! " << endl;
1324 exit(1);
1325 }
1326 for (unsigned int j = 0; j < 2; j++)
1327 {
1328 h_tofec_p_bg_offset[idx][i][j] = (TH1D*)f_tofec_bg[idx][i]->Get( Form("h_tofec_p_bg_offset_%s_%d", name, j) );
1329 h_tofec_m_bg_offset[idx][i][j] = (TH1D*)f_tofec_bg[idx][i]->Get( Form("h_tofec_m_bg_offset_%s_%d", name, j) );
1330 h_tofec_p_bg_sigma[idx][i][j] = (TH1D*)f_tofec_bg[idx][i]->Get( Form("h_tofec_p_bg_sigma_%s_%d" , name, j) );
1331 h_tofec_m_bg_sigma[idx][i][j] = (TH1D*)f_tofec_bg[idx][i]->Get( Form("h_tofec_m_bg_sigma_%s_%d" , name, j) );
1332 }
1333 }
1334 //tof_endcap cost
1335 filename.clear();
1336 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_cost_d10.root", dir) );
1337 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_cost_d11.root", dir) );
1338 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_cost_m10.root", dir) );
1339 filename.push_back( path + Form("/share/%s/tof_endcap/tofec_cost_m11.root", dir) );
1340 for (unsigned int i = 0; i < filename.size(); i++)
1341 {
1342 f_tofec_cost[idx][i] = new TFile(filename[i].c_str(), "READ");
1343 const char *name;
1344 if (i == 0)
1345 name = "d10";
1346 else if (i == 1)
1347 name = "d11";
1348 else if (i == 2)
1349 name = "m10";
1350 else if (i == 3)
1351 name = "m11";
1352 else
1353 {
1354 cout << "Boundary Error! " << endl;
1355 exit(1);
1356 }
1357 for (unsigned int j = 0; j < 2; j++)
1358 {
1359 h_tofec_p_cost_offset[idx][i][j] = (TH1D*)f_tofec_cost[idx][i]->Get( Form("h_tofec_p_cost_offset_%s_%d", name, j) );
1360 h_tofec_m_cost_offset[idx][i][j] = (TH1D*)f_tofec_cost[idx][i]->Get( Form("h_tofec_m_cost_offset_%s_%d", name, j) );
1361 h_tofec_p_cost_sigma[idx][i][j] = (TH1D*)f_tofec_cost[idx][i]->Get( Form("h_tofec_p_cost_sigma_%s_%d" , name, j) );
1362 h_tofec_m_cost_sigma[idx][i][j] = (TH1D*)f_tofec_cost[idx][i]->Get( Form("h_tofec_m_cost_sigma_%s_%d" , name, j) );
1363 }
1364 }
1365 }
1366 for (unsigned int idx = 0; idx < 3; idx++)
1367 {
1368 const char *dir;
1369 if (idx == 0)
1370 dir = "electron/emc";
1371 else if (idx == 1)
1372 dir = "kpi/emc_pion";
1373 else if (idx == 2)
1374 dir = "kpi/emc_kaon";
1375 else
1376 {
1377 cout << "Boundary Error! " << endl;
1378 exit(1);
1379 }
1380 //emc
1381 filename.clear();
1382 filename.push_back( path + Form("/share/%s/emc_d10.root", dir) );
1383 filename.push_back( path + Form("/share/%s/emc_d11.root", dir) );
1384 filename.push_back( path + Form("/share/%s/emc_m10.root", dir) );
1385 filename.push_back( path + Form("/share/%s/emc_m11.root", dir) );
1386 for (unsigned int i = 0; i < filename.size(); i++)
1387 {
1388 f_emc[idx][i] = new TFile(filename[i].c_str(), "READ");
1389 const char *name;
1390 if (i == 0)
1391 name = "d10";
1392 else if (i == 1)
1393 name = "d11";
1394 else if (i == 2)
1395 name = "m10";
1396 else if (i == 3)
1397 name = "m11";
1398 else
1399 {
1400 cout << "Boundary Error! " << endl;
1401 exit(1);
1402 }
1403 for (unsigned int j = 0; j < 15; j++)
1404 {
1405 for (unsigned int k = 0; k < 25; k++)
1406 {
1407 h_emc_ep[idx][i][j][k] = (TH1D*)f_emc[idx][i]->Get( Form("h_ep_%s_%d_%d", name, j, k) );
1408 h_emc_e35[idx][i][j][k] = (TH1D*)f_emc[idx][i]->Get( Form("h_e35_%s_%d_%d", name, j, k) );
1409 }
1410 }
1411 }
1412 }
1413 cout << "Successfully Return from Loading Initializations by package SimplePIDSvc ... " << endl;
1414}
1415
1416int SimplePIDSvc::findBin(double *a, int length, double value)
1417{
1418 for (int i = 0; i < length; i++)
1419 {
1420 if (value > a[i] && value <= a[i+1])
1421 {
1422 return i;
1423 }
1424 }
1425 if (value < a[0])
1426 {
1427 return 0;
1428 }
1429 else
1430 {
1431 return length;
1432 }
1433}
1434
1435double SimplePIDSvc::getChi2(int i)
1436{
1437 return pow(m_dedx_chi[i], 2) + pow(m_tof_chi[i], 2);
1438}
1439
1440void SimplePIDSvc::loadEMCInfo(EvtRecTrack *track)
1441{
1442 //Initialization
1443 for (unsigned int i = 0; i < 5; i++)
1444 {
1445 m_emc_eop[i] = -99.;
1446 m_emc_likelihood[i] = -99.;
1447 }
1448 m_emc_e = -99.;
1449 m_emc_e13 = -99.;
1450 m_emc_e35 = -99.;
1451 m_emc_sec = -99.;
1452 m_emc_lat = -99.;
1453 m_lh_electron = -99.;
1454
1455 if (!track->isEmcShowerValid()) return;
1456
1457 RecEmcShower* emc_trk = track->emcShower();
1458 m_emc_e = emc_trk->energy();
1459 for (unsigned int i = 0; i < 5; i++)
1460 {
1461 m_emc_eop[i] = m_emc_e / fabs(m_p[i]);
1462 }
1463 double eseed = emc_trk->eSeed();
1464 double e3 = emc_trk->e3x3();
1465 double e5 = emc_trk->e5x5();
1466 m_emc_sec = emc_trk->secondMoment() / 1000.;
1467 m_emc_lat = emc_trk->latMoment();
1468 if (e3 != 0)
1469 {
1470 m_emc_e13 = eseed / e3;
1471 }
1472 if (e5 != 0)
1473 {
1474 m_emc_e35 = e3 / e5;
1475 }
1476}
1477
1478bool SimplePIDSvc::calEMCLikelihood()
1479{
1480 if (m_emc_eop[0] < 0)
1481 return false;
1482
1483 int idx = getRunIdx(m_run);
1484 const Int_t BIN_P = 15;
1485 const Int_t BIN_COST = 25;
1486 const Int_t BIN_PID = 3;
1487 const double EPS = 1e-4;
1488 //electron
1489 double P[BIN_PID][BIN_P + 1] = {
1490 {0.20, 0.47, 0.56, 0.65, 0.72, 0.79, 0.86, 0.92, 0.98, 1.03, 1.08, 1.13, 1.17, 1.22, 1.26, 1.30},
1491 {0.20, 0.26, 0.31, 0.35, 0.39, 0.42, 0.46, 0.49, 0.53, 0.57, 0.62, 0.67, 0.73, 0.80, 0.88, 1.05},
1492 {0.20, 0.33, 0.39, 0.43, 0.48, 0.52, 0.56, 0.61, 0.67, 0.73, 0.76, 0.81, 0.85, 0.90, 0.96, 1.05}, };
1493 double COST[BIN_PID][BIN_COST + 1] = {
1494 {-0.930, -0.910, -0.905, -0.897, -0.890, -0.881, -0.871, -0.858, -0.775, -0.732, -0.669, -0.561, -0.330, 0.199, 0.515, 0.645, 0.718, 0.766, 0.804, 0.870, 0.882, 0.891, 0.898, 0.906, 0.913, 0.930},
1495 {-0.930, -0.810, -0.728, -0.648, -0.574, -0.501, -0.431, -0.364, -0.295, -0.228, -0.161, -0.096, -0.031, 0.035, 0.100, 0.167, 0.234, 0.301, 0.370, 0.439, 0.510, 0.580, 0.655, 0.733, 0.813, 0.930},
1496 {-0.930, -0.804, -0.721, -0.643, -0.568, -0.497, -0.429, -0.362, -0.293, -0.228, -0.161, -0.096, -0.029, 0.035, 0.100, 0.166, 0.233, 0.298, 0.365, 0.432, 0.500, 0.571, 0.644, 0.722, 0.805, 0.930}, };
1497
1498 double vp, vcost;
1499 int pid;
1500 int bin_p, bin_cost;
1501 for (unsigned int i = 0; i < 4; i++)
1502 {
1503 //only e, pi ,K
1504 if (i == 0)
1505 pid = 0;
1506 else if (i == 2)
1507 pid = 1;
1508 else if (i == 3)
1509 pid = 2;
1510 else
1511 continue;
1512
1513 vp = max(P[pid][0]+EPS, min(P[pid][BIN_P]-EPS, m_p[i]));
1514 vcost = max(COST[pid][0]+EPS, min(COST[pid][BIN_COST]-EPS, m_cost[i]));
1515 bin_p = findBin(P[pid], BIN_P, vp);
1516 bin_cost = findBin(COST[pid], BIN_COST, vcost);
1517
1518 m_emc_likelihood[i] = h_emc_ep[pid][idx][bin_p][bin_cost]->Interpolate(m_emc_eop[i]) * h_emc_e35[pid][idx][bin_p][bin_cost]->Interpolate(m_emc_e35);
1519 }
1520 double a = m_prob[0] > 0 ? m_prob[0] : 0;
1521 double b = m_prob[2] > 0 ? m_prob[2] : 0;
1522 double c = m_prob[3] > 0 ? m_prob[3] : 0;
1523 double sum = a * m_emc_likelihood[0] + b * m_emc_likelihood[2] + c * m_emc_likelihood[3];
1524
1525 if (sum > 0 && m_prob[0] > 0)
1526 {
1527 m_lh_electron = m_prob[0] * m_emc_likelihood[0] / sum;
1528 return true;
1529 }
1530 else
1531 {
1532 return false;
1533 }
1534}
1535
1536bool SimplePIDSvc::ispion()
1537{
1538 if (m_prob[2] > 0.00 && m_prob[2] > m_prob[3])
1539 return true;
1540 else
1541 return false;
1542}
1543
1544bool SimplePIDSvc::iskaon()
1545{
1546 if (m_prob[3] > 0.00 && m_prob[3] > m_prob[2])
1547 return true;
1548 else
1549 return false;
1550}
1551
1552//bool SimplePIDSvc::iselectron_(bool emc)
1553//{
1554// if (m_prob[0] > 0 && m_prob[0] > m_prob[2] && m_prob[0] > m_prob[3])
1555// {
1556// if (!emc)
1557// return true;
1558// else if (fabs(m_cost[0]) < 0.7 && m_emc_eop[0] > 0 && m_emc_eop[0] < 0.8)
1559// return false;
1560// else if (fabs(m_cost[0]) >= 0.7 && fabs(m_cost[0])<0.8 && m_emc_eop[0] > 0 && m_emc_eop[0] < -7.5*fabs(m_cost[0])+6.05)
1561// return false;
1562// else if (fabs(m_cost[0]) > 0.85 && m_emc_eop[0] > 0 && m_emc_eop[0] < 0.6)
1563// return false;
1564// else
1565// return true;
1566// }
1567// else
1568// return false;
1569//}
1570
1571bool SimplePIDSvc::iselectron(bool emc)
1572{
1573 if (!emc)
1574 {
1575 if (m_prob[0] > 0 && m_prob[0] > m_prob[2] && m_prob[0] > m_prob[3])
1576 return true;
1577 else
1578 return false;
1579 }
1580 else
1581 {
1582 if (calEMCLikelihood())
1583 {
1584 if (m_lh_electron > m_eid_ratio)
1585 return true;
1586 else
1587 return false;
1588 }
1589 else
1590 {
1591 if (m_prob[0] > 0 && m_prob[0] > m_prob[2] && m_prob[0] > m_prob[3])
1592 return true;
1593 else
1594 return false;
1595 }
1596 }
1597}
cos
double cos(const BesAngle a)
Definition BesAngle.h:213
P
double P(RecMdcKalTrack *trk)
Definition CalibEventSelect.cxx:87
mass
double mass
Definition CosmicGenerator.cxx:138
sigma
TTree * sigma
Definition CurConst_dEdXCal.cxx:6
q
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
Definition KKsem.h:33
dt
TGraph2DErrors * dt
Definition McCor.cxx:45
EPS
const double EPS
Definition TRunge.cxx:43
t
TTree * t
Definition binning.cxx:23
DstEmcShower::latMoment
double latMoment() const
Definition DstEmcShower.h:52
DstEmcShower::eSeed
double eSeed() const
Definition DstEmcShower.h:47
DstEmcShower::e3x3
double e3x3() const
Definition DstEmcShower.h:48
DstEmcShower::secondMoment
double secondMoment() const
Definition DstEmcShower.h:51
DstEmcShower::e5x5
double e5x5() const
Definition DstEmcShower.h:49
DstEmcShower::energy
double energy() const
Definition DstEmcShower.h:45
DstExtTrack::tof1Position
const Hep3Vector tof1Position() const
Definition DstExtTrack.h:58
DstExtTrack::tof2Position
const Hep3Vector tof2Position() const
Definition DstExtTrack.h:94
DstMdcDedx::chi
double chi(int i) const
Definition DstMdcDedx.h:58
DstMdcKalTrack::theta
const double theta() const
Definition DstMdcKalTrack.h:70
DstMdcKalTrack::setPidType
static void setPidType(PidType pidType)
Definition DstMdcKalTrack.h:59
DstMdcKalTrack::p
const double p() const
Definition DstMdcKalTrack.h:82
DstMdcKalTrack::charge
const int charge() const
Definition DstMdcKalTrack.h:65
EvtRecTrack::isMdcDedxValid
bool isMdcDedxValid()
Definition EvtRecTrack.h:45
EvtRecTrack::mdcDedx
RecMdcDedx * mdcDedx()
Definition EvtRecTrack.h:55
EvtRecTrack::isExtTrackValid
bool isExtTrackValid()
Definition EvtRecTrack.h:49
EvtRecTrack::extTrack
RecExtTrack * extTrack()
Definition EvtRecTrack.h:56
EvtRecTrack::isMdcKalTrackValid
bool isMdcKalTrackValid()
Definition EvtRecTrack.h:44
EvtRecTrack::tofTrack
SmartRefVector< RecTofTrack > tofTrack()
Definition EvtRecTrack.h:57
EvtRecTrack::isTofTrackValid
bool isTofTrackValid()
Definition EvtRecTrack.h:46
EvtRecTrack::emcShower
RecEmcShower * emcShower()
Definition EvtRecTrack.h:58
EvtRecTrack::isEmcShowerValid
bool isEmcShowerValid()
Definition EvtRecTrack.h:47
EvtRecTrack::mdcKalTrack
RecMdcKalTrack * mdcKalTrack()
Definition EvtRecTrack.h:54
SimplePIDSvc::SimplePIDSvc
SimplePIDSvc(const std::string &name, ISvcLocator *svcLoc)
Definition SimplePIDSvc.cxx:21
SimplePIDSvc::~SimplePIDSvc
virtual ~SimplePIDSvc()
Definition SimplePIDSvc.cxx:30
SimplePIDSvc::getChi2
double getChi2(int i)
Definition SimplePIDSvc.cxx:1435
SimplePIDSvc::preparePID
void preparePID(EvtRecTrack *track)
Definition SimplePIDSvc.cxx:95
SimplePIDSvc::initialize
virtual StatusCode initialize()
Definition SimplePIDSvc.cxx:32
SimplePIDSvc::iselectron
bool iselectron(bool emc=true)
Definition SimplePIDSvc.cxx:1571
SimplePIDSvc::finalize
virtual StatusCode finalize()
Definition SimplePIDSvc.cxx:47
TofHitStatus::is_barrel
bool is_barrel() const
Definition TofHitStatus.h:26
TofHitStatus::layer
unsigned int layer() const
Definition TofHitStatus.h:28
TofHitStatus::is_cluster
bool is_cluster() const
Definition TofHitStatus.h:25
TofHitStatus::setStatus
void setStatus(unsigned int status)
Definition TofHitStatus.cxx:134
TofHitStatus::is_counter
bool is_counter() const
Definition TofHitStatus.h:24
TofHitStatus::is_readout
bool is_readout() const
Definition TofHitStatus.h:23
TofHitStatus::is_raw
bool is_raw() const
Definition TofHitStatus.h:22
const
EvtCyclic3::c_str
char * c_str(Index i)
Definition EvtCyclic3.cc:252
bpkt::ifstream
std::ifstream ifstream
Definition bpkt_streams.h:44
charge
float charge
Definition prepare_data.cxx:28
bg
float bg
Definition prepare_data.cxx:28
b
const double b
Definition slope.cxx:9