ATLAS Offline Software
MuCCaFitter.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 #include "MdtCalibFitters/MuCCaFitter.h"
6 
7 #include <cmath>
8 
9 #include "AthenaKernel/getMessageSvc.h"
10 #include "GaudiKernel/MsgStream.h"
11 
12 namespace MuonCalib {
13 
14  void MuCCaFitter::printLevel(int level) {
15  if (level > 0) m_debug = true;
16  }
17 
18  bool MuCCaFitter::fit(MuonCalibSegment& seg) const {
19  // select all hits
20  HitSelection selection(seg.mdtHitsOnTrack(), 0);
21  // call fit function
22  return fit(seg, selection);
23  }
24 
25  bool MuCCaFitter::fit(MuonCalibSegment& seg, HitSelection selection) const {
26  MsgStream log(Athena::getMessageSvc(), "MuCCaFitter");
27  if (log.level() <= MSG::DEBUG) log << MSG::DEBUG << "fit() New seg:" << endmsg;
28 
29  int N = seg.mdtHitsOnTrack();
30 
31  if (N < 2) { return false; }
32 
33  if ((int)selection.size() != N) {
34  selection.clear();
35  selection.assign(N, 0);
36  } else {
37  int used(0);
38  for (int i = 0; i < N; ++i) {
39  if (selection[i] == 0) ++used;
40  }
41  if (used < 2) {
42  log << MSG::WARNING << "fit() TO FEW HITS SELECTED" << endmsg;
43  return false;
44  }
45  }
46 
47  const Amg::Vector3D& pos = seg.position();
48 
49  double S(0), Sy(0), Sz(0), Zc{0}, Yc{0};
50  std::vector<double> y(N), x(N), z(N);
51  std::vector<double> r(N), sr(N), w(N), rw(N);
52  int ii{0}, jj{0};
53  {
54  for (const MuonCalibSegment::MdtHitPtr& hit : seg.mdtHOT()) {
55  const MdtCalibHitBase& h = *hit;
56  y[ii] = getY(h.localPosition());
57  z[ii] = getZ(h.localPosition());
58  r[ii] = std::abs(h.driftRadius());
59 
60  if (h.sigma2DriftRadius() > 0) {
61  w[ii] = 1. / (h.sigma2DriftRadius());
62  sr[ii] = std::sqrt(h.sigma2DriftRadius());
63  } else {
64  w[ii] = 0.;
65  sr[ii] = 0.;
66  }
67  if (log.level() <= MSG::DEBUG)
68  log << MSG::DEBUG << "fit() MuCCaFitter: (" << x[ii] << "," << y[ii] << ") R = " << r[ii] << " W " << w[ii] << endmsg;
69  rw[ii] = r[ii] * w[ii];
70  if (selection[jj]) {
71  ++jj;
72  continue;
73  }
74  S += w[ii];
75  Sz += w[ii] * z[ii];
76  Sy += w[ii] * y[ii];
77  ++ii;
78  ++jj;
79  }
80  }
81  Zc = Sz / S;
82  Yc = Sy / S;
83 
84  std::unique_ptr<MuCCaFitterImplementation> Fitter = std::make_unique<MuCCaFitterImplementation>();
85  if (log.level() <= MSG::DEBUG) {
86  for (int i = 0; i != ii; i++) {
87  log << MSG::DEBUG << "fit() MuCCaFitter hits passed to computepam Z=" << x[i] << " phi=" << y[i] << " zzz=" << z[i]
88  << " drift=" << r[i] << " error=" << w[i] << endmsg;
89  }
90  }
91  Fitter->Computeparam3(ii, z, y, r, sr);
92  if (log.level() <= MSG::DEBUG) {
93  log << MSG::DEBUG << "fit() MuCCaFitter computed track a=" << Fitter->get_a() << " da= " << Fitter->get_da()
94  << " b= " << Fitter->get_b() << " db= " << Fitter->get_db() << " corrab= " << Fitter->get_corrab()
95  << " chi2f= " << Fitter->get_chi2f() << endmsg;
96  }
97 
98  double afit = Fitter->get_a();
99  double chi2f = Fitter->get_chi2f();
100 
101  std::vector<double> dist(N), ddist(N);
102  double R, dis;
103  double theta = std::atan(afit);
104  if (theta < 0.) theta = M_PI + theta;
105  double sinus = std::sin(theta);
106  double cosin = std::cos(theta);
107  double d = -(getZ(pos) - Zc) * sinus + (getY(pos) - Yc) * cosin;
108  if (log.level() <= MSG::DEBUG) {
109  log << MSG::DEBUG << "fit() MuCCaFitter>>> theta= " << theta << " sinus= " << sinus << " cosin= " << cosin << endmsg;
110  log << MSG::DEBUG << "fit() MuCCaFitter>>> getZ( pos )= " << getZ(pos) << " getY( pos )= " << getY(pos) << " Zc= " << Zc
111  << " Yc= " << Yc << endmsg;
112  }
113  double Syy(0), Szy(0), Syyzz(0), Att(0);
114  R = 0;
115  for (int i = 0; i < N; ++i) {
116  if (selection[i]) continue;
117  Syy += (y[i] - Yc) * (y[i] - Yc) * w[i];
118  Szy += (y[i] - Yc) * (z[i] - Zc) * w[i];
119  Syyzz += ((y[i] - Yc) - (z[i] - Zc)) * ((y[i] - Yc) + (z[i] - Zc)) * w[i];
120  dis = (y[i] - Yc) * cosin - (z[i] - Zc) * sinus;
121  if (dis > d) {
122  R -= rw[i];
123  } else {
124  R += rw[i];
125  }
126  }
127  Att = Syy + cosin * (2 * sinus * Szy - cosin * Syyzz);
128  d = R / S;
129  if (log.level() <= MSG::DEBUG)
130  log << MSG::DEBUG << "fit() MuCCaFitter>>> d= " << d << " R= " << R << " S= " << S << " Att= " << Att << endmsg;
131  for (int i = 0; i < N; ++i) {
132  if (selection[i]) continue;
133  dist[i] = cosin * (y[i] - Yc) - sinus * (z[i] - Zc) - d;
134  double dth = -(sinus * (y[i] - Yc) + cosin * (z[i] - Zc)) * (std::sqrt(1. / Att));
135  ddist[i] = std::sqrt(dth * dth + (1. / S));
136  }
137  if (log.level() <= MSG::DEBUG) log << MSG::DEBUG << "fit() Transforming back to real world" << endmsg;
138  Amg::Vector3D ndir = getVec(0., sinus, cosin);
139  Amg::Vector3D npos = getVec(0., Yc + cosin * d, Zc - sinus * d);
140 
141  if (log.level() <= MSG::DEBUG)
142  log << MSG::DEBUG << "fit() New line: position " << npos << " direction " << ndir << " chi2f " << chi2f << endmsg;
143 
144  seg.set(chi2f / (N - 2), npos, ndir);
145 
146  int i{0};
147  for (const MuonCalibSegment::MdtHitPtr& hit_ptr : seg.mdtHOT()) {
148  hit_ptr->setDistanceToTrack(dist[i], ddist[i]);
149  ++i;
150  }
151 
152  if (log.level() <= MSG::DEBUG) log << MSG::DEBUG << "fit() fit done" << endmsg;
153  return true;
154  }
155  void MuCCaFitterImplementation::Computeparam3(int number_of_hits, const std::vector<double>& x, const std::vector<double>& y, const std::vector<double>& r,
156  const std::vector<double>& sr) {
157  /***************************************/
158  /* Fit a line to n=number_of_hits hits */
159  /***************************************/
160  double xout[100], yout[100];
161  double dist[100], rsub[100];
162  double chi2outn[4], chi2min;
163  double aref[4], bref[4];
164  int fhit, i, j, icouple = 0, lhit;
165  /* compute 4 tanget lines from first and last hit */
166  fhit = 0;
167  lhit = number_of_hits - 1;
168  Computelinparnew(x[fhit], y[fhit], r[fhit], x[lhit], y[lhit], r[lhit]);
169  for (int ii = 0; ii < 4; ii++) {
170  bref[ii] = m_bfpar[ii];
171  aref[ii] = m_angularcoefficient[ii];
172  }
173  /* choose the REFERENCE LINE (aref[icouple],bref[icouple])*/
174  for (i = 0; i < 4; i++) {
175  for (j = 0; j < number_of_hits; j++) {
176  double d;
177  d = std::abs(aref[i] * x[j] + bref[i] - y[j]);
178  dist[j] = d / std::sqrt(1. + aref[i] * aref[i]);
179  rsub[j] = r[j] - dist[j];
180  }
181  double chi2out = 0;
182  for (int ii = 1; ii < number_of_hits - 1; ii++) { chi2out += rsub[ii] * rsub[ii] / (sr[ii] * sr[ii]); }
183  chi2outn[i] = chi2out;
184  }
185  chi2min = 999999999.;
186  for (i = 0; i < 4; i++) {
187  if (chi2outn[i] < chi2min) {
188  chi2min = chi2outn[i];
189  icouple = i;
190  }
191  }
192  /* Compute TRACK POINTS */
193  Computehitpsemes(number_of_hits, x, y, r, aref[icouple], bref[icouple]);
194  for (i = 0; i < number_of_hits; i++) {
195  xout[i] = m_xpoint[i];
196  yout[i] = m_ypoint[i];
197  }
198  /* Compute a & b parameters of the track from TRACK POINTS */
199  // int idim;
200  // ifail;
201  double aoutn, bout, sig2a, sig2b, corrab;
202  double temp, det;
203  double hesse[2][2];
204  double W, WX, WX2, WY, WXY;
205  // double errormatrix[2][2];
206  W = 0.;
207  WX = 0.;
208  WX2 = 0.;
209  WY = 0.;
210  WXY = 0;
211  for (int i = 0; i < number_of_hits; i++) {
212  temp = 1. / (sr[i] * sr[i]);
213  W += temp;
214  WX += xout[i] * temp;
215  WX2 += xout[i] * xout[i] * temp;
216  WY += yout[i] * temp;
217  WXY += xout[i] * yout[i] * temp;
218  }
219  det = W * WX2 - WX * WX;
220  aoutn = (W * WXY - WY * WX) / det;
221  bout = (WY * WX2 - WX * WXY) / det;
222  hesse[1][1] = W;
223  hesse[0][0] = WX2;
224  hesse[1][0] = WX;
225  hesse[0][1] = WX;
226  // idim = 2;
227  /* invert hessian matrix */
228  hesse[1][1] = hesse[0][0] / det;
229  hesse[0][0] = hesse[1][1] / det;
230  hesse[1][0] = -1. * (hesse[0][1] / det);
231  hesse[0][1] = -1. * (hesse[0][1] / det);
232  sig2a = hesse[0][0];
233  sig2b = hesse[1][1];
234  corrab = hesse[1][0];
235  double deno, btest, bs1, bs2, db1, db2;
236  btest = bout;
237  bout = 0;
238  deno = 0;
239  for (int i = 0; i < number_of_hits; i++) {
240  bs1 = (y[i] - aoutn * x[i] - r[i] * std::sqrt(1 + aoutn * aoutn));
241  bs2 = (y[i] - aoutn * x[i] + r[i] * std::sqrt(1 + aoutn * aoutn));
242  db1 = std::abs(bs1 - btest);
243  db2 = std::abs(bs2 - btest);
244  if (db1 < db2) {
245  bout += bs1 / (sr[i] * sr[i]);
246  } else {
247  bout += bs2 / (sr[i] * sr[i]);
248  }
249  deno += 1 / (sr[i] * sr[i]);
250  }
251  bout /= deno;
252  /* compute chi2 from residuals*/
253  double resd;
254  double chi2f = 0;
255  for (int i = 0; i < number_of_hits; i++) {
256  double xi, yi;
257  xi = (aoutn * y[i] - aoutn * bout + x[i]) / (1. + aoutn * aoutn);
258  yi = aoutn * xi + bout;
259  resd = (std::sqrt((xi - x[i]) * (xi - x[i]) + (yi - y[i]) * (yi - y[i])) - r[i]) / sr[i];
260  chi2f += resd * resd;
261  }
262  /* set global variables (method output)*/
263  m_aoutn = aoutn;
264  m_sig2a = sig2a;
265  m_bout = bout;
266  m_sig2b = sig2b;
267  m_corrab = corrab;
268  m_chi2f = chi2f;
269  }
270 
271  double MuCCaFitterImplementation::get_a() { return m_aoutn; }
272  double MuCCaFitterImplementation::get_b() { return m_bout; }
273  double MuCCaFitterImplementation::get_da() { return m_sig2a; }
274  double MuCCaFitterImplementation::get_db() { return m_sig2b; }
275  double MuCCaFitterImplementation::get_corrab() { return m_corrab; }
276  double MuCCaFitterImplementation::get_chi2f() { return m_chi2f; }
277 
278  void MuCCaFitterImplementation::Computelinparnew(double x1, double y1, double r1, double x2, double y2, double r2) {
279  /********************************************/
280  /* Compute the 4 lines tangent to 2 circles */
281  /********************************************/
282  double delta;
283  double averagephi, dist, dphiin, dphiex;
284  double bpar[4], phi;
285  double bfparn[2];
286  bfparn[0] = 0;
287  bfparn[1] = 0;
288  // int segnobfn[2];
289  // int segnoc[2][4],ncandid[4],ncand;
290  int ncandid[4], ncand;
291  double bcand[2][4];
292  int segnob[] = {1, -1, 1, -1};
293  int firsttime, i;
294  double angularcoefficient[4];
295  double bfpar[4];
296  // int segnobf[4];
297  double dy, dx;
298 
299  /* compute a parameters */
300  dx = x2 - x1;
301  dy = y2 - y1;
302  averagephi = std::atan2(dy, dx);
303  dist = std::sqrt((x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1));
304  delta = r2 + r1;
305 
306  dphiin = std::asin(delta / dist);
307 
308  delta = r2 - r1;
309  dphiex = std::asin(delta / dist);
310 
311  int f = 1;
312  phi = averagephi + f * dphiex;
313  if (phi < 0) { phi = 2 * M_PI + (phi); }
314  angularcoefficient[0] = std::tan(phi);
315 
316  f = -1;
317  phi = averagephi + f * dphiex;
318  if (phi < 0) { phi = 2 * M_PI + (phi); }
319  angularcoefficient[1] = std::tan(phi);
320 
321  f = 1;
322  phi = averagephi + f * dphiin;
323  if (phi < 0) { phi = 2 * M_PI + (phi); }
324  angularcoefficient[2] = std::tan(phi);
325 
326  f = -1;
327  phi = averagephi + f * dphiin;
328  if (phi < 0) { phi = 2 * M_PI + (phi); }
329  angularcoefficient[3] = std::tan(phi);
330 
331  /* Compute b parameters */
332  for (i = 0; i < 4; i++) {
333  bpar[0] = y1 - angularcoefficient[i] * x1 + segnob[0] * r1 * std::sqrt(1 + angularcoefficient[i] * angularcoefficient[i]);
334  bpar[1] = y1 - angularcoefficient[i] * x1 + segnob[1] * r1 * std::sqrt(1 + angularcoefficient[i] * angularcoefficient[i]);
335  bpar[2] = y2 - angularcoefficient[i] * x2 + segnob[2] * r2 * std::sqrt(1 + angularcoefficient[i] * angularcoefficient[i]);
336  bpar[3] = y2 - angularcoefficient[i] * x2 + segnob[3] * r2 * std::sqrt(1 + angularcoefficient[i] * angularcoefficient[i]);
337  double delta = 0.00001;
338  ncand = 0;
339  if (std::abs(bpar[0] - bpar[2]) < delta) {
340  bfparn[ncand] = bpar[0];
341  ncand = ncand + 1;
342  }
343  if (std::abs(bpar[0] - bpar[3]) < delta) {
344  bfparn[ncand] = bpar[0];
345  ncand = ncand + 1;
346  }
347  if (std::abs(bpar[1] - bpar[2]) < delta) {
348  bfparn[ncand] = bpar[1];
349  ncand = ncand + 1;
350  }
351  if (std::abs(bpar[1] - bpar[3]) < delta) {
352  bfparn[ncand] = bpar[1];
353  ncand = ncand + 1;
354  }
355  bcand[0][i] = bfparn[0];
356  bcand[1][i] = bfparn[1];
357  ncandid[i] = ncand;
358  }
359  firsttime = 0;
360  for (i = 0; i < 4; i++) {
361  if (ncandid[i] == 1) {
362  bfpar[i] = bcand[0][i];
363  } else {
364  bfpar[i] = bcand[firsttime][i];
365  firsttime++;
366  }
367  }
368  /* set global variables (method output)*/
369  for (i = 0; i < 4; i++) {
370  m_bfpar[i] = bfpar[i];
371  m_angularcoefficient[i] = angularcoefficient[i];
372  }
373  }
374 
375  void MuCCaFitterImplementation::Computelin(double x1, double y1, double r1, double x2, double y2, double r2) {
376  /********************************************/
377  /* Compute the 4 lines tangent to 2 circles */
378  /********************************************/
379  double delta{0}, averagephi{0}, dist{0}, dphiin{0}, dphiex{0};
380  double bpar[4], phi{0};
381  double bfparn[2];
382  bfparn[0] = 0;
383  bfparn[1] = 0;
384  double bcand[2][4];
385  int segnob[] = {1, -1, 1, -1};
386  int ncandid[4], ncand;
387  int i{0}, firsttime{0};
388  double angularcoefficient[4];
389  double bfpar[4];
390  double dy{0}, dx{0};
391 
392  /* compute a parameters */
393  dx = x2 - x1;
394  dy = y2 - y1;
395  averagephi = std::atan2(dy, dx);
396  dist = std::hypot(dx, dy);
397 
398  delta = r2 + r1;
399  dphiin = std::asin(delta / dist);
400 
401  delta = r2 - r1;
402  dphiex = std::asin(delta / dist);
403 
404  int f = 1;
405  phi = averagephi + f * dphiex;
406  if (phi < 0) { phi = 2 * M_PI + (phi); }
407  angularcoefficient[0] = std::tan(phi);
408 
409  f = -1;
410  phi = averagephi + f * dphiex;
411  if (phi < 0) { phi = 2 * M_PI + (phi); }
412  angularcoefficient[1] = std::tan(phi);
413 
414  f = 1;
415  phi = averagephi + f * dphiin;
416  if (phi < 0) { phi = 2 * M_PI + (phi); }
417  angularcoefficient[2] = std::tan(phi);
418 
419  f = -1;
420  phi = averagephi + f * dphiin;
421  if (phi < 0) { phi = 2 * M_PI + (phi); }
422  angularcoefficient[3] = std::tan(phi);
423 
424  /* Compute b parameters */
425  for (i = 0; i < 4; i++) {
426  bpar[0] = y1 - angularcoefficient[i] * x1 + segnob[0] * r1 * std::sqrt(1 + angularcoefficient[i] * angularcoefficient[i]);
427  bpar[1] = y1 - angularcoefficient[i] * x1 + segnob[1] * r1 * std::sqrt(1 + angularcoefficient[i] * angularcoefficient[i]);
428  bpar[2] = y2 - angularcoefficient[i] * x2 + segnob[2] * r2 * std::sqrt(1 + angularcoefficient[i] * angularcoefficient[i]);
429  bpar[3] = y2 - angularcoefficient[i] * x2 + segnob[3] * r2 * std::sqrt(1 + angularcoefficient[i] * angularcoefficient[i]);
430  double delta = 0.00001;
431  ncand = 0;
432  if (std::abs(bpar[0] - bpar[2]) < delta) {
433  bfparn[ncand] = bpar[0];
434  ncand = ncand + 1;
435  }
436  if (std::abs(bpar[0] - bpar[3]) < delta) {
437  bfparn[ncand] = bpar[0];
438  ncand = ncand + 1;
439  }
440  if (std::abs(bpar[1] - bpar[2]) < delta) {
441  bfparn[ncand] = bpar[1];
442  ncand = ncand + 1;
443  }
444  if (std::abs(bpar[1] - bpar[3]) < delta) {
445  bfparn[ncand] = bpar[1];
446  ncand = ncand + 1;
447  }
448  bcand[0][i] = bfparn[0];
449  bcand[1][i] = bfparn[1];
450  ncandid[i] = ncand;
451  }
452  firsttime = 0;
453  for (i = 0; i < 4; i++) {
454  if (ncandid[i] == 1) {
455  bfpar[i] = bcand[0][i];
456  } else {
457  bfpar[i] = bcand[firsttime][i];
458  firsttime++;
459  }
460  }
461  /* set global variables (method output)*/
462  for (int i = 0; i < 4; i++) {
463  m_bfpar[i] = bfpar[i];
464  m_angularcoefficient[i] = angularcoefficient[i];
465  }
466  }
467  void MuCCaFitterImplementation::Computehitpsemes(int nhit, const std::vector<double>& xcirc, const std::vector<double>& ycirc,
468  const std::vector<double>& rcirc, double a, double b) {
469  /************************/
470  /* Compute TRACK POINTS */
471  /************************/
472  int i, nhitc;
473  // double xpoint[nhit],ypoint[nhit];
474  /* loop over hit couples */
475  nhitc = nhit / 2;
476  for (i = 0; i < nhitc; i++) {
477  /* Compute TRACK POINTS for 2 circles */
478  computehitsfromcircles(xcirc[2 * i], ycirc[2 * i], rcirc[2 * i], xcirc[2 * i + 1], ycirc[2 * i + 1], rcirc[2 * i + 1], a, b);
479  m_xpoint[2 * i] = m_xout0;
480  m_ypoint[2 * i] = m_yout0;
481  m_xpoint[2 * i + 1] = m_xout1;
482  m_ypoint[2 * i + 1] = m_yout1;
483  }
484  if (2 * nhitc != nhit) {
485  /* Compute TRACK POINTS for 2 circles */
486  computehitsfromcircles(xcirc[nhit - 2], ycirc[nhit - 2], rcirc[nhit - 2], xcirc[nhit - 1], ycirc[nhit - 1], rcirc[nhit - 1], a,
487  b);
488  m_xpoint[nhit - 2] = m_xout0;
489  m_ypoint[nhit - 2] = m_yout0;
490  m_xpoint[nhit - 1] = m_xout1;
491  m_ypoint[nhit - 1] = m_yout1;
492  }
493  }
494 
495  void MuCCaFitterImplementation::computehitsfromcircles(double x0, double y0, double r0, double x1, double y1, double r1, double a,
496  double b) {
497  /**************************************/
498  /* Compute TRACK POINTS for 2 circles */
499  /**************************************/
500  double ang[4], bpar[4];
501  int i;
502  double xout0 = 0, yout0 = 0;
503  double xout1 = 0, yout1 = 0;
504  double dist, dist1, xint, yint, xref, yref;
505 
506  /* Compute the 4 lines tangent to 2 circles */
507  Computelin(x0, y0, r0, x1, y1, r1);
508  for (i = 0; i < 4; i++) {
509  bpar[i] = m_bfpar[i];
510  ang[i] = m_angularcoefficient[i];
511  }
512 
513  dist = 9999999;
514  /* for each of the 4 tangents : */
515  for (i = 0; i < 4; i++) {
516  /* compute tangent point */
517  xint = (x0 + ang[i] * y0 - ang[i] * bpar[i]) / (1 + ang[i] * ang[i]);
518  yint = ang[i] * (xint) + bpar[i];
519  /* compute point from reference line */
520  double bprime;
521  if (a != 0) {
522  bprime = y0 + x0 / a;
523  xref = (bprime - b) * a / (a * a + 1);
524  yref = (b + bprime * a * a) / (a * a + 1);
525  } else {
526  xref = x0;
527  yref = b;
528  }
529  /* choose tangent point closest to that of the reference lines (TRACK POINT)*/
530  dist1 = std::sqrt((xint - xref) * (xint - xref) + (yint - yref) * (yint - yref));
531  if (dist1 < dist) {
532  dist = dist1;
533  xout0 = xint;
534  yout0 = yint;
535  }
536  }
537  dist = 9999999;
538  for (i = 0; i < 4; i++) {
539  /* compute tangent point */
540  xint = (x1 + ang[i] * y1 - ang[i] * bpar[i]) / (1 + ang[i] * ang[i]);
541  yint = ang[i] * (xint) + bpar[i];
542  /* compute point from reference line */
543  double bprime;
544  if (a != 0) {
545  bprime = y1 + x1 / a;
546  xref = (bprime - b) * a / (a * a + 1);
547  yref = (b + bprime * a * a) / (a * a + 1);
548  } else {
549  xref = x1;
550  yref = b;
551  }
552  /* choose tangent point closest to that of the reference line (TRACK POINT)*/
553  dist1 = std::sqrt((xint - xref) * (xint - xref) + (yint - yref) * (yint - yref));
554  if (dist1 < dist) {
555  dist = dist1;
556  xout1 = xint;
557  yout1 = yint;
558  }
559  }
560  /* set global variables (method output)*/
561  m_xout0 = xout0;
562  m_yout0 = yout0;
563  m_xout1 = xout1;
564  m_yout1 = yout1;
565  }
566 } // namespace MuonCalib
used
MuonCalib::MuCCaFitterImplementation::Computelin
void Computelin(double x1, double y1, double r1, double x2, double y2, double r2)
Definition: MuCCaFitter.cxx:375
beamspotman.r
def r
Definition: beamspotman.py:676
plotBeamSpotCompare.x1
x1
Definition: plotBeamSpotCompare.py:216
MuonCalib::MuCCaFitter::getZ
double getZ(const Amg::Vector3D &p) const
Definition: MuCCaFitter.h:58
getMessageSvc.h
singleton-like access to IMessageSvc via open function and helper
MuonCalib::MuCCaFitterImplementation::m_angularcoefficient
double m_angularcoefficient[4]
parameters of the 4 tangent lines from first and last hit
Definition: MuCCaFitter.h:88
TRTCalib_Extractor.det
det
Definition: TRTCalib_Extractor.py:36
MuonCalib::MuCCaFitterImplementation::get_da
double get_da()
Definition: MuCCaFitter.cxx:273
JetTiledMap::W
@ W
Definition: TiledEtaPhiMap.h:44
hist_file_dump.d
d
Definition: hist_file_dump.py:137
plotBeamSpotCompare.x2
x2
Definition: plotBeamSpotCompare.py:218
MuonCalib::MuCCaFitterImplementation::m_ypoint
double m_ypoint[100]
Definition: MuCCaFitter.h:97
M_PI
#define M_PI
Definition: ActiveFraction.h:11
MuonCalib::MuonCalibSegment
Definition: MuonCalibSegment.h:39
MuonCalib::MuonCalibSegment::position
const Amg::Vector3D & position() const
retrieve local position of segment (on station level)
Definition: MuonCalibSegment.cxx:185
MuonCalib::MuCCaFitterImplementation::m_yout1
double m_yout1
Definition: MuCCaFitter.h:94
MuonCalib::MuCCaFitterImplementation::Computelinparnew
void Computelinparnew(double x1, double y1, double r1, double x2, double y2, double r2)
Definition: MuCCaFitter.cxx:278
JetTiledMap::N
@ N
Definition: TiledEtaPhiMap.h:44
MuonCalib::MuCCaFitterImplementation::get_chi2f
double get_chi2f()
Definition: MuCCaFitter.cxx:276
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
MuonCalib::MuCCaFitter::fit
bool fit(MuonCalibSegment &seg) const
fit using all hits
Definition: MuCCaFitter.cxx:18
MuonCalib::MuCCaFitterImplementation::m_bfpar
double m_bfpar[4]
Definition: MuCCaFitter.h:89
checkFileSG.rw
rw
Definition: checkFileSG.py:124
x
#define x
JetTiledMap::S
@ S
Definition: TiledEtaPhiMap.h:44
Athena::getMessageSvc
IMessageSvc * getMessageSvc(bool quiet=false)
Definition: getMessageSvc.cxx:20
makeTRTBarrelCans.y1
tuple y1
Definition: makeTRTBarrelCans.py:15
MuonCalib::MuCCaFitterImplementation::m_xpoint
double m_xpoint[100]
track points
Definition: MuCCaFitter.h:96
drawFromPickle.atan
atan
Definition: drawFromPickle.py:36
MuonCalib::MuCCaFitterImplementation::m_xout0
double m_xout0
track points from 2 circles
Definition: MuCCaFitter.h:91
MuonCalib::MuCCaFitterImplementation::m_bout
double m_bout
track constant term
Definition: MuCCaFitter.h:101
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
MuonCalib::MuCCaFitterImplementation::m_sig2a
double m_sig2a
track slope's variance
Definition: MuCCaFitter.h:100
bsCompare.db1
int db1
Definition: bsCompare.py:40
MuonCalib::MuCCaFitterImplementation::get_b
double get_b()
Definition: MuCCaFitter.cxx:272
MuonCalib::MuCCaFitterImplementation::m_xout1
double m_xout1
Definition: MuCCaFitter.h:93
python.SystemOfUnits.sr
int sr
Definition: SystemOfUnits.py:113
DeMoUpdate.db2
string db2
Definition: DeMoUpdate.py:1032
lumiFormat.i
int i
Definition: lumiFormat.py:85
MuonCalib::MuCCaFitterImplementation::computehitsfromcircles
void computehitsfromcircles(double x0, double y0, double r0, double x1, double y1, double r1, double a, double b)
Definition: MuCCaFitter.cxx:495
z
#define z
MuonCalib::MuCCaFitterImplementation::m_chi2f
double m_chi2f
chisquared
Definition: MuCCaFitter.h:104
MuonCalib::MuCCaFitterImplementation::m_sig2b
double m_sig2b
error on track constant term
Definition: MuCCaFitter.h:102
endmsg
#define endmsg
Definition: AnalysisConfig_Ntuple.cxx:63
makeTRTBarrelCans.y2
tuple y2
Definition: makeTRTBarrelCans.py:18
MuonCalib::MuCCaFitterImplementation::get_corrab
double get_corrab()
Definition: MuCCaFitter.cxx:275
MuonCalib
CscCalcPed - algorithm that finds the Cathode Strip Chamber pedestals from an RDO.
Definition: CscCalcPed.cxx:22
hist_file_dump.f
f
Definition: hist_file_dump.py:135
drawFromPickle.tan
tan
Definition: drawFromPickle.py:36
MuonCalib::MuCCaFitterImplementation::get_db
double get_db()
Definition: MuCCaFitter.cxx:274
MuonCalib::MuonCalibSegment::mdtHitsOnTrack
unsigned int mdtHitsOnTrack() const
retrieve the number of MdtCalibHitBase s assigned to this segment
Definition: MuonCalibSegment.cxx:146
selection
const std::string selection
Definition: fbtTestBasics.cxx:74
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:77
MuonCalib::MuCCaFitterImplementation::m_aoutn
double m_aoutn
results
Definition: MuCCaFitter.h:99
MuonCalib::MuCCaFitterImplementation::m_yout0
double m_yout0
Definition: MuCCaFitter.h:92
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
python.LumiBlobConversion.pos
pos
Definition: LumiBlobConversion.py:18
MuonCalib::MuCCaFitter::getVec
Amg::Vector3D getVec(double x, double y, double z) const
Definition: MuCCaFitter.h:59
makeTRTBarrelCans.dy
tuple dy
Definition: makeTRTBarrelCans.py:21
MuonCalib::MuonCalibSegment::set
void set(double chi2, const Amg::Vector3D &pos, const Amg::Vector3D &dir)
Definition: MuonCalibSegment.cxx:128
MuonCalib::MuCCaFitter::m_debug
bool m_debug
Definition: MuCCaFitter.h:46
a
TList * a
Definition: liststreamerinfos.cxx:10
y
#define y
h
MuonCalib::MdtCalibHitBase
Definition: MdtCalibHitBase.h:38
DEBUG
#define DEBUG
Definition: page_access.h:11
makeTRTBarrelCans.dx
tuple dx
Definition: makeTRTBarrelCans.py:20
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
MuonCalib::MuCCaFitterImplementation::Computeparam3
void Computeparam3(int number_of_hits, const std::vector< double > &x, const std::vector< double > &y, const std::vector< double > &r, const std::vector< double > &sr)
Definition: MuCCaFitter.cxx:155
MuonCalib::MuCCaFitterImplementation::m_corrab
double m_corrab
correlation term
Definition: MuCCaFitter.h:103
MuonCalib::MuCCaFitter::getY
double getY(const Amg::Vector3D &p) const
these methods are needed to change the reference frame between the local one of the hit and one used ...
Definition: MuCCaFitter.h:56
MuonCalib::MuonCalibSegment::mdtHOT
const MdtHitVec & mdtHOT() const
retrieve the full set of MdtCalibHitBase s assigned to this segment
Definition: MuonCalibSegment.cxx:147
MuonCalib::MuCCaFitterImplementation::get_a
double get_a()
Definition: MuCCaFitter.cxx:271
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
python.IoTestsLib.w
def w
Definition: IoTestsLib.py:200
MuonCalib::MuCCaFitter::printLevel
void printLevel(int level)
set print level
Definition: MuCCaFitter.cxx:14
MuCCaFitter.h
MuonCalib::MuonCalibSegment::MdtHitPtr
std::shared_ptr< MdtCalibHitBase > MdtHitPtr
typedef for a collection of MdtCalibHitBase s
Definition: MuonCalibSegment.h:44
MuonCalib::MuCCaFitterImplementation::Computehitpsemes
void Computehitpsemes(int nhit, const std::vector< double > &xcirc, const std::vector< double > &ycirc, const std::vector< double > &rcirc, double a, double b)
Definition: MuCCaFitter.cxx:467
MuonCalib::IMdtSegmentFitter::HitSelection
std::vector< unsigned int > HitSelection
Definition: IMdtSegmentFitter.h:32
Generated on Sun Dec 22 2024 21:14:49 for ATLAS Offline Software by   doxygen 1.8.18