Class Reference for E1039 Core & Analysis Software
GenFitExtrapolator.cxx
Go to the documentation of this file.
1 #include "GenFitExtrapolator.h"
2 #include "GFField.h"
3 
4 #include <phool/recoConsts.h>
5 #include <phfield/PHFieldUtility.h>
6 #include <phfield/PHFieldConfig_v3.h>
7 #include <phfield/PHField.h>
8 
9 #include <GenFit/FieldManager.h>
10 #include <GenFit/MaterialEffects.h>
11 #include <GenFit/TGeoMaterialInterface.h>
12 #include <GenFit/MeasuredStateOnPlane.h>
13 #include <GenFit/SharedPlanePtr.h>
14 #include <GenFit/RKTrackRep.h>
15 
16 //ROOT
17 #include <TGeoManager.h>
18 #include <TMath.h>
19 
20 //
21 //#include <G4SystemOfUnits.hh>
22 
23 #include <memory>
24 #include <cassert>
25 
26 //#define _DEBUG_ON
27 
28 namespace
29 {
30  //static flag to indicate the initialized has been done
31  static bool inited = false;
32 
33  //Simple swimming settings
34  static int NSTEPS_TARGET = 100;
35  static int NSTEPS_SHIELDING = 50;
36  static int NSTEPS_FMAG = 100;
37 
38  static double FMAG_LENGTH;
39  static double Z_UPSTREAM;
40 
41  //initialize global variables
42  void initGlobalVariables()
43  {
44  if(!inited)
45  {
46  inited = true;
47  recoConsts* rc = recoConsts::instance();
48 
49  NSTEPS_TARGET = rc->get_IntFlag("NSTEPS_TARGET");
50  NSTEPS_SHIELDING = rc->get_IntFlag("NSTEPS_SHIELDING");
51  NSTEPS_FMAG = rc->get_IntFlag("NSTEPS_FMAG");
52 
53  FMAG_LENGTH = rc->get_DoubleFlag("FMAG_LENGTH");
54  Z_UPSTREAM = rc->get_DoubleFlag("Z_UPSTREAM");
55  }
56  }
57 }
58 
59 static const double c_light = 2.99792458e+8 * m/s;
60 
61 using namespace std;
62 
63 GenFitExtrapolator::GenFitExtrapolator():
64  pos_i(TVector3()), mom_i(TVector3()), cov_i(TMatrixDSym(5)),
65  pos_f(TVector3()), mom_f(TVector3()), cov_f(TMatrixDSym(5)),
66  jac_sd2sc(TMatrixD(5,5)), jac_sc2sd(TMatrixD(5,5)), propM(TMatrixD(5,5)),
67  jac_genfit2legacy(TMatrixD(5,5)), jac_legacy2genfit(TMatrixD(5,5))
68 {
69  initGlobalVariables();
70 }
71 
72 GenFitExtrapolator::~GenFitExtrapolator()
73  {}
74 
75 bool GenFitExtrapolator::init(const PHField* field, const TGeoManager *geom)
76 {
77  iParType = 1;
78 
79  assert(field);
80  SQGenFit::GFField *fieldMap = new SQGenFit::GFField(field);
81  genfit::FieldManager::getInstance()->init(fieldMap);
82 
83  _tgeo_manager = const_cast<TGeoManager*>(geom);
84 
85 #ifdef _DEBUG_ON
86  double z_test = 1000;
87  LogInfo("");
88  {
89  double p[4] = {0, 0, z_test*cm, 0};
90  double B[3] = {0, 0, 0};
91  field->GetFieldValue(p, B);
92  cout << "PHField (CLHEP) at Z = " << z_test << endl;
93  cout << B[0] << ", " << B[1] << ", " << B[2] << endl;
94  }
95  {
96  genfit::AbsBField *f = genfit::FieldManager::getInstance()->getField();
97  TVector3 H = f->get(TVector3(0,0,z_test));
98  H *= kilogauss/tesla;
99  cout << "genfit::AbsBField (tesla) at Z = " << z_test << endl;
100  H.Print();
101 
102  z_test = 250;
103  H = f->get(TVector3(0,0,z_test));
104  H *= kilogauss/tesla;
105  cout << "genfit::AbsBField (tesla) at Z = " << z_test << endl;
106  H.Print();
107  }
108 #endif
109 
110  _tgeo_manager->Export("GenFitExtrapolatorGeom.root");
111 
112  genfit::MaterialEffects::getInstance()->init(new genfit::TGeoMaterialInterface());
113 
114  return true;
115 }
116 
117 void GenFitExtrapolator::setInitialStateWithCov(
118  double z_in,
119  TMatrixD& state_in,
120  TMatrixD& cov_in)
121 {
122 #ifdef _DEBUG_ON
123  LogInfo("z_in: ") << z_in << endl;
124  state_in.Print();
125  cov_in.Print();
126 #endif
127 
128  convertSVtoMP(z_in, state_in, mom_i, pos_i);
129 
130  if(state_in[0][0] > 0)
131  {
132  setParticleType(1);
133  }
134  else
135  {
136  setParticleType(-1);
137  }
138 
140 /*
141  TMatrixD cov_sd(5, 5);
142  for(int i = 0; i < 5; i++)
143  {
144  for(int j = 0; j < 5; j++)
145  {
146  cov_sd[i][j] = cov_in[i][j];
147 
148  if(i == 0) cov_sd[i][j] = double(iParType)*cov_sd[i][j];
149  if(j == 0) cov_sd[i][j] = double(iParType)*cov_sd[i][j];
150  }
151  }
152 
154  TRSDSC(iParType, mom_i, pos_i);
155  TMatrixD jac_sd2sc_T = jac_sd2sc;
156  jac_sd2sc_T.T();
157 
158  TMatrixD cov_sc = jac_sd2sc*cov_sd*jac_sd2sc_T;
159  for(int i = 0; i < 5; i++)
160  {
161  for(int j = 0; j < 5; j++)
162  {
163  cov_i[i][j] = cov_sc[i][j];
164  }
165  }
166 
167 */
168 
170  TRLEGACY2GENFIT(iParType, mom_i, pos_i);
171  TMatrixD jac_legacy2genfit_T = jac_legacy2genfit;
172  jac_legacy2genfit_T.T();
173  TMatrixD cov_hold= jac_legacy2genfit*cov_in*jac_legacy2genfit_T;
174 
175  for(int i = 0; i < 5; i++)
176  {
177  for(int j = 0; j < 5; j++)
178  {
179  cov_i[i][j] = cov_hold[i][j];
180  }
181  }
182 
183 
184 }
185 
186 void GenFitExtrapolator::setParticleType(int type) {
187  iParType = type;
188 }
189 
190 void GenFitExtrapolator::getFinalStateWithCov(TMatrixD& state_out, TMatrixD& cov_out) {
191  convertMPtoSV(mom_f, pos_f, state_out);
192 
193 /*
195  TMatrixD cov_sc(5, 5);
196  for(int i = 0; i < 5; i++)
197  {
198  for(int j = 0; j < 5; j++)
199  {
200  cov_sc[i][j] = cov_f[i][j];
201 
202  if(direction == GenFitExtrapolator::Backward)
203  {
204  if(i == 1) cov_sc[i][j] = -cov_sc[i][j];
205  if(j == 1) cov_sc[i][j] = -cov_sc[i][j];
206  if(i == 3) cov_sc[i][j] = -cov_sc[i][j];
207  if(j == 3) cov_sc[i][j] = -cov_sc[i][j];
208  }
209  }
210  }
211 
212  //convert from SC to SD error matrix
213  TRSCSD(iParType, mom_f, pos_f);
214  TMatrixD jac_sc2sd_T = jac_sc2sd;
215  jac_sc2sd_T.T();
216 
217  cov_out = jac_sc2sd*cov_sc*jac_sc2sd_T;
218 
219  for(int i = 0; i < 5; i++)
220  {
221  for(int j = 0; j < 5; j++)
222  {
223  if(i == 0) cov_out[i][j] = double(iParType)*cov_out[i][j];
224  if(j == 0) cov_out[i][j] = double(iParType)*cov_out[i][j];
225  }
226  }
227 */
228 
230  TRGENFIT2LEGACY(iParType, mom_i, pos_i);
231  TMatrixD jac_genfit2legacy_T = jac_genfit2legacy;
232  jac_genfit2legacy_T.T();
233  TMatrixD cov_hold= jac_genfit2legacy*cov_f*jac_genfit2legacy_T;
234 
235  for(int i = 0; i < 5; i++)
236  {
237  for(int j = 0; j < 5; j++)
238  {
239  cov_out[i][j] = cov_hold[i][j];
240  }
241  }
242 
243 }
244 
245 void GenFitExtrapolator::getPropagator(TMatrixD& prop) {
246  if(fabs(pos_i[2] - pos_f[2]) < 1E-3)
247  {
248  prop.UnitMatrix();
249  return;
250  }
251 
253 
254 // for(int i = 0; i < 5; i++)
255 // {
256 // for(int j = 0; j < 5; j++)
257 // {
258 // prop[i][j] = g4eProp[i][j];
259 //
260 // if(i == 0) prop[i][j] = double(iParType)*prop[i][j];
261 // if(j == 0) prop[i][j] = double(iParType)*prop[i][j];
262 // if(direction == GenFitExtrapolator::Backward)
263 // {
264 // if(i == 1) prop[i][j] = -prop[i][j];
265 // if(j == 1) prop[i][j] = -prop[i][j];
266 // if(i == 3) prop[i][j] = -prop[i][j];
267 // if(j == 3) prop[i][j] = -prop[i][j];
268 // }
269 // }
270 // }
271 //
272 // prop = jac_sc2sd*prop*jac_sd2sc;
273 
274  for(int i = 0; i < 5; i++)
275  {
276  for(int j = 0; j < 5; j++)
277  {
278  prop[i][j] = propM[i][j];
279 
280  }
281 
282  }
283 
284  prop = jac_legacy2genfit*prop*jac_genfit2legacy;
285 
286 }
287 
288 bool GenFitExtrapolator::extrapolateTo(double z_out) {
289 
291  if(pos_i[2] > 2400 || pos_i[2] < -600 || z_out > 2400 || z_out < -600)
292  {
293  return false;
294  }
295 
297  if(fabs(pos_i[2] - z_out) < 1E-3)
298  {
299  mom_f = mom_i;
300  pos_f = pos_i;
301  cov_f = cov_i;
302 
303  return true;
304  }
305 
306  int pid = iParType > 0 ? -13 : 13;
307 #ifdef _DEBUG_ON
308  cout
309  << "extrapolateToPlane: "
310  << "From: " << pos_i.Z()
311  << " -> " << z_out
312  << endl;
313 #endif
314 
315  auto up_rep = std::unique_ptr<genfit::AbsTrackRep> (new genfit::RKTrackRep(pid));
316  auto rep = up_rep.get();
318  genfit::SharedPlanePtr destPlane(new genfit::DetPlane(TVector3(0, 0, z_out), TVector3(0, 0, 1)));
319 
320  std::unique_ptr<genfit::MeasuredStateOnPlane> currentState = std::unique_ptr < genfit::MeasuredStateOnPlane > (new genfit::MeasuredStateOnPlane(rep));
322  genfit::SharedPlanePtr startplane(new genfit::DetPlane(TVector3(0, 0, pos_i[2]), TVector3(0, 0, 1)));
323  currentState->setPlane(startplane);
324  //rep->setPosMomCov(*currentState, pos_i, mom_i, cov_i);
325  currentState->setPosMom(pos_i, mom_i);
326  currentState->setCov(cov_i);
327 
328  try {
329  travelLength = rep->extrapolateToPlane(*currentState, destPlane);
330  } catch (...) {
331 #ifdef _DEBUG_ON
332  cout << "extrapolateToPlane failed!" << endl;
333 #endif
334  return false;
335  }
336 
337  currentState->getPosMom(pos_f, mom_f);
338  cov_f = currentState->getCov();
339  propM = rep->getJacobian();
340 
341  return true;
342 }
343 
344 //int GenFitExtrapolator::propagate() {
345 //
346 // return -1;
347 //}
348 
349 double GenFitExtrapolator::extrapolateToIP() {
350 
351  //Store the steps on each point
352  TVector3 mom[NSTEPS_FMAG + NSTEPS_TARGET + 1];
353  TVector3 pos[NSTEPS_FMAG + NSTEPS_TARGET + 1];
354 
355  //Step size in FMAG/target area, unit is cm.
356  //FIXME Units
357  double step_fmag = FMAG_LENGTH/NSTEPS_FMAG;//*cm;
358  double step_target = fabs(Z_UPSTREAM)/NSTEPS_TARGET;//*cm;
359 
360  //Start from FMAG face downstream
361  extrapolateTo(FMAG_LENGTH);
362  pos[0] = pos_f;
363  mom[0] = mom_f;
364 
365  //Now make the real swimming
366  int iStep = 1;
367  for(; iStep <= NSTEPS_FMAG; ++iStep)
368  {
369  pos_i = pos[iStep-1];
370  mom_i = mom[iStep-1];
371 
372  extrapolateTo((pos_i[2] - step_fmag));
373 
374  pos[iStep] = pos_f;
375  mom[iStep] = mom_f;
376  }
377 
378  for(; iStep < NSTEPS_FMAG+NSTEPS_TARGET+1; ++iStep)
379  {
380  pos_i = pos[iStep-1];
381  mom_i = mom[iStep-1];
382 
383  extrapolateTo((pos_i[2] - step_target));
384 
385  pos[iStep] = pos_f;
386  mom[iStep] = mom_f;
387  }
388 
389  //Find the one step with minimum DCA
390  double dca_min = 1E6;
391  for(int i = 0; i < NSTEPS_FMAG+NSTEPS_TARGET+1; ++i)
392  {
393  double dca = sqrt(pos[i][0]*pos[i][0] + pos[i][1]*pos[i][1]);
394  if(dca < dca_min)
395  {
396  dca_min = dca;
397  iStep = i;
398  }
399  }
400 
401  return pos[iStep][2];
402 }
403 void GenFitExtrapolator::convertSVtoMP(double z, TMatrixD& state, TVector3& mom, TVector3& pos)
404 {
405  double p = fabs(1. / state[0][0]);
406  double pz = p / sqrt(1. + state[1][0] * state[1][0] + state[2][0] * state[2][0]);
407  double px = pz * state[1][0];
408  double py = pz * state[2][0];
409 
410  double x = state[3][0];
411  double y = state[4][0];
412 
413  pos.SetXYZ(x, y, z);
414  mom.SetXYZ(px, py, pz);
415 }
416 
417 void GenFitExtrapolator::convertMPtoSV(TVector3& mom, TVector3& pos, TMatrixD& state) {
418  TVector3 mom_gev = mom;
419  TVector3 pos_cm = pos;
420 
421  state[0][0] = double(iParType) / mom_gev.Mag();
422  state[1][0] = mom_gev.x() / mom_gev.z();
423  state[2][0] = mom_gev.y() / mom_gev.z();
424 
425  state[3][0] = pos_cm.x();
426  state[4][0] = pos_cm.y();
427 }
428 
429 void GenFitExtrapolator::TRLEGACY2GENFIT(int charge, TVector3 mom_input, TVector3 pos_input)
430 {
432  jac_legacy2genfit.Zero();
433  jac_legacy2genfit[0][0]=1.;
434  jac_legacy2genfit[1][1]=-1.;
435  jac_legacy2genfit[2][2]=-1.;
436  jac_legacy2genfit[3][3]=-1.;
437  jac_legacy2genfit[4][4]=-1.;
438 
439 }
440 
441 void GenFitExtrapolator::TRGENFIT2LEGACY(int charge, TVector3 mom_input, TVector3 pos_input)
442 {
444  jac_genfit2legacy.Zero();
445  jac_genfit2legacy[0][0]=1.;
446  jac_genfit2legacy[1][1]=-1.;
447  jac_genfit2legacy[2][2]=-1.;
448  jac_genfit2legacy[3][3]=-1.;
449  jac_genfit2legacy[4][4]=-1.;
450 
451 }
452 
453 void GenFitExtrapolator::TRSDSC(int charge, TVector3 mom_input, TVector3 pos_input)
454 {
456  TVector3 mom(mom_input[0], mom_input[1], mom_input[2]);
457  TVector3 pos(pos_input[0], pos_input[1], pos_input[2]);
458 
460  TVector3 DJ(1., 0., 0.);
461  TVector3 DK(0., 1., 0.);
462  TVector3 DI = DJ.Cross(DK);
463 
465  double p_inv = 1./mom.Mag();
466  double vp = mom.Dot(DJ)/mom.Dot(DI);
467  double wp = mom.Dot(DK)/mom.Dot(DI);
468  double lambda = TMath::Pi()/2. - mom.Theta();
469  //double phi = mom.Phi();
470 
471  TVector3 TVW;
472  TVW.SetX(1./sqrt(1. + vp*vp + wp*wp));
473  TVW.SetY(vp*TVW.X());
474  TVW.SetZ(wp*TVW.X());
475 
476  TVector3 TN;
477  for(int i = 0; i < 3; i++)
478  {
479  TN[i] = TVW[0]*DI[i] + TVW[1]*DJ[i] + TVW[2]*DK[i];
480  }
481 
482  double cosl = cos(lambda);
483  double cosl1 = 1./cosl;
484 
485  TVector3 UN(-TN.Y()*cosl1, TN.X()*cosl1, 0.);
486  TVector3 VN(-TN.Z()*UN.Y(), TN.Z()*UN.X(), cosl);
487 
488  double UJ = UN.Dot(DJ);
489  double UK = UN.Dot(DK);
490  double VJ = VN.Dot(DJ);
491  double VK = VN.Dot(DK);
492 
493  jac_sd2sc.Zero();
494  jac_sd2sc[0][0] = 1.;
495  jac_sd2sc[1][1] = TVW[0]*VJ;
496  jac_sd2sc[1][2] = TVW[0]*VK;
497  jac_sd2sc[2][1] = TVW[0]*UJ*cosl1;
498  jac_sd2sc[2][2] = TVW[0]*UK*cosl1;
499  jac_sd2sc[3][3] = UJ;
500  jac_sd2sc[3][4] = UK;
501  jac_sd2sc[4][3] = VJ;
502  jac_sd2sc[4][4] = VK;
503 
504  //Takes cm output kGauss (0.1*tesla)
505  genfit::AbsBField *field = genfit::FieldManager::getInstance()->getField();
506 
507  if(charge != 0 && field)
508  {
509  TVector3 H = field->get(pos);
510  H *= kilogauss/tesla;
511 #ifdef _DEBUG_ON
512  LogInfo("");
513  pos.Print();
514  H.Print();
515 #endif
516  double HA = H.Mag();
517  double HAM = HA*p_inv*tesla*GeV;
518  double HM;
519  if(HA < 1E-6)
520  {
521  HM = 0.;
522  }
523  else
524  {
525  HM = charge/HA;
526  }
527 
528  double Q = -HAM*c_light/(km/ns);
529  double sinz = -H.Dot(UN)*HM;
530  double cosz = H.Dot(VN)*HM;
531 
532  jac_sd2sc[1][3] = -Q*TVW[1]*sinz;
533  jac_sd2sc[1][4] = -Q*TVW[2]*sinz;
534  jac_sd2sc[2][3] = -Q*TVW[1]*cosz*cosl1;
535  jac_sd2sc[2][4] = -Q*TVW[2]*cosz*cosl1;
536  }
537 }
538 
539 void GenFitExtrapolator::TRSCSD(int charge, TVector3 mom_input, TVector3 pos_input)
540 {
542  TVector3 mom(mom_input[0], mom_input[1], mom_input[2]);
543  TVector3 pos(pos_input[0], pos_input[1], pos_input[2]);
544 
546  TVector3 DJ(1., 0., 0.);
547  TVector3 DK(0., 1., 0.);
548  TVector3 DI = DJ.Cross(DK);
549 
551  double p_inv = 1./mom.Mag();
552  //double vp = mom.Dot(DJ)/mom.Dot(DI);
553  //double wp = mom.Dot(DK)/mom.Dot(DI);
554  double lambda = TMath::Pi()/2. - mom.Theta();
555  double phi = mom.Phi();
556 
557  double cosl = cos(lambda);
558  double sinp = sin(phi);
559  double cosp = cos(phi);
560 
561  TVector3 TN(cosl*cosp, cosl*sinp, sin(lambda));
562  TVector3 TVW;
563  TVW.SetX(TN.Dot(DI));
564  TVW.SetY(TN.Dot(DJ));
565  TVW.SetZ(TN.Dot(DK));
566 
567  double T1R = 1./TVW[0];
568  double T2R = T1R*T1R;
569  TVector3 UN(-sinp, cosp, 0.);
570  TVector3 VN(-TN.Z()*UN.Y(), TN.Z()*UN.X(), cosl);
571 
572  double UJ = UN.Dot(DJ);
573  double UK = UN.Dot(DK);
574  double VJ = VN.Dot(DJ);
575  double VK = VN.Dot(DK);
576  double UI = UN.Dot(DI);
577  double VI = VN.Dot(DI);
578 
579  jac_sc2sd.Zero();
580  jac_sc2sd[0][0] = 1.;
581  jac_sc2sd[1][1] = -UK*T2R;
582  jac_sc2sd[1][2] = VK*cosl*T2R;
583  jac_sc2sd[2][1] = UJ*T2R;
584  jac_sc2sd[2][2] = -VJ*cosl*T2R;
585  jac_sc2sd[3][3] = VK*T1R;
586  jac_sc2sd[3][4] = -UK*T1R;
587  jac_sc2sd[4][3] = -VJ*T1R;
588  jac_sc2sd[4][4] = UJ*T1R;
589 
590  genfit::AbsBField *field = genfit::FieldManager::getInstance()->getField();
591  if(charge != 0 && field)
592  {
593  TVector3 H = field->get(pos);
594  H *= kilogauss/tesla;
595 
596  double HA = H.Mag();
597  double HAM = HA*p_inv;
598  double HM;
599  if(HA < 1E-6)
600  {
601  HM = 0.;
602  }
603  else
604  {
605  HM = charge/HA;
606  }
607 
608  double Q = -HAM*c_light/(km/ns);
609  double sinz = -H.Dot(UN)*HM;
610  double cosz = H.Dot(VN)*HM;
611  double T3R = Q*T1R*T1R*T1R;
612 
613  jac_sc2sd[1][3] = -UI*(VK*cosz - UK*sinz)*T3R;
614  jac_sc2sd[1][4] = -VI*(VK*cosz - UK*sinz)*T3R;
615  jac_sc2sd[2][3] = UI*(VJ*cosz - UJ*sinz)*T3R;
616  jac_sc2sd[2][4] = VI*(VJ*cosz - UJ*sinz)*T3R;
617  }
618 }
619 
620 void GenFitExtrapolator::print() {
621  cout << "Propagating " << iParType << ":" << endl;
622  cout << "From "
623  << "(" << pos_i.X() << ", " << pos_i.Y() << ", "<< pos_i.Z() << ") "
624  << " To "
625  << "(" << pos_f.X() << ", " << pos_f.Y() << ", "<< pos_f.Z() << ") "
626  << endl;
627 
628  cout << "Momentum change: From "
629  << "(" << mom_i.X() << ", " << mom_i.Y() << ", "<< mom_i.Z() << ") "
630  << " To "
631  << "(" << mom_f.X() << ", " << mom_f.Y() << ", "<< mom_f.Z() << ") "
632  << endl;
633 
634  cout << "Initial error matrix: " << endl;
635  for(int i = 0; i < 5; i++)
636  {
637  for(int j = 0; j < 5; j++)
638  {
639  cout << cov_i[i][j] << " ";
640  }
641  cout << endl;
642  }
643 
644  cout << "Final error matrix: " << endl;
645  for(int i = 0; i < 5; i++)
646  {
647  for(int j = 0; j < 5; j++)
648  {
649  cout << cov_f[i][j] << " ";
650  }
651  cout << endl;
652  }
653 }
LogInfo
#define LogInfo(message)
Definition: DbSvc.cc:15
c_light
static const double c_light
Definition: GenFitExtrapolator.cxx:59
GenFitExtrapolator::init
bool init(const PHField *field, const TGeoManager *geom)
Initialize geometry and physics.
Definition: GenFitExtrapolator.cxx:75
GenFitExtrapolator::convertMPtoSV
void convertMPtoSV(TVector3 &mom, TVector3 &pos, TMatrixD &state)
Definition: GenFitExtrapolator.cxx:417
GenFitExtrapolator::convertSVtoMP
void convertSVtoMP(double z, TMatrixD &state, TVector3 &mom, TVector3 &pos)
Transformation between the state vector and the mom/pos.
Definition: GenFitExtrapolator.cxx:403
GenFitExtrapolator::extrapolateToIP
double extrapolateToIP()
Extrapolate to the primary vertex.
Definition: GenFitExtrapolator.cxx:349
GenFitExtrapolator::extrapolateTo
bool extrapolateTo(double z_out)
Definition: GenFitExtrapolator.cxx:288
GenFitExtrapolator::TRSDSC
void TRSDSC(int charge, TVector3 mom_input, TVector3 pos_input)
Definition: GenFitExtrapolator.cxx:453
GenFitExtrapolator::setParticleType
void setParticleType(int type)
Set particle type.
Definition: GenFitExtrapolator.cxx:186
GenFitExtrapolator::TRLEGACY2GENFIT
void TRLEGACY2GENFIT(int charge, TVector3 mom_input, TVector3 pos_input)
Definition: GenFitExtrapolator.cxx:429
GenFitExtrapolator::getPropagator
void getPropagator(TMatrixD &prop)
Get the propagator.
Definition: GenFitExtrapolator.cxx:245
GenFitExtrapolator::print
void print()
Debug print.
Definition: GenFitExtrapolator.cxx:620
GenFitExtrapolator::TRGENFIT2LEGACY
void TRGENFIT2LEGACY(int charge, TVector3 mom_input, TVector3 pos_input)
Tranformation between GenFit and Legacy plane; Abi.
Definition: GenFitExtrapolator.cxx:441
GenFitExtrapolator::getFinalStateWithCov
void getFinalStateWithCov(TMatrixD &state_out, TMatrixD &cov_out)
Get the final state parameters and covariance.
Definition: GenFitExtrapolator.cxx:190
GenFitExtrapolator::setInitialStateWithCov
void setInitialStateWithCov(double z_in, TMatrixD &state_in, TMatrixD &cov_in)
Set input initial state parameters.
Definition: GenFitExtrapolator.cxx:117
GenFitExtrapolator::TRSCSD
void TRSCSD(int charge, TVector3 mom_input, TVector3 pos_input)
Definition: GenFitExtrapolator.cxx:539
PHField
transient DST object for field storage and access
Definition: PHField.h:14
PHField::GetFieldValue
virtual void GetFieldValue(const double Point[4], double *Bfield) const =0
PHFlag::get_DoubleFlag
virtual double get_DoubleFlag(const std::string &name) const
Definition: PHFlag.cc:49
PHFlag::get_IntFlag
virtual int get_IntFlag(const std::string &name) const
Definition: PHFlag.cc:117
recoConsts::instance
static recoConsts * instance()
Definition: recoConsts.cc:7