62 #ifndef _DHELICAL_FIT_H_
63 #define _DHELICAL_FIT_H_
74 #include "JANA/jerror.h"
78 #define atan2f(x,y) atan2((double)x,(double)y)
109 jerror_t AddStereoHit(
const DCDCWire *wire);
110 jerror_t AddHit(
float r,
float phi,
float z);
111 jerror_t AddHitXYZ(
float x,
float y,
float z);
112 jerror_t AddHitXYZ(
float x,
float y,
float z,
float covx,
float covy,
113 float covxy,
bool is_axial=
false);
115 jerror_t PruneHit(
int idx);
116 jerror_t
Clear(
void);
117 jerror_t FitCircle(
void);
118 double ChisqCircle(
void);
119 jerror_t FitCircleRiemann(
float rc=0.);
120 jerror_t FitCircleRiemann(
float z_vertex,
float BeamRMS);
121 jerror_t FitLineRiemann(
void);
122 jerror_t FitCircleStraightTrack();
123 void SearchPtrans(
double ptrans_max=9.0,
double ptrans_step=0.5);
124 void QuickPtrans(
void);
125 jerror_t GuessChargeFromCircleFit(
void);
126 jerror_t FitTrack(
void);
127 jerror_t FitTrackRiemann(
float rc);
128 jerror_t FitCircleAndLineRiemann(
float rc);
129 jerror_t FitTrack_FixedZvertex(
float z_vertex);
130 jerror_t FitLine_FixedZvertex(
float z_vertex);
131 jerror_t Fill_phi_circle(vector<DHFHit_t*> hits,
float x0,
float y0);
132 inline const vector<DHFHit_t*>
GetHits()
const {
return hits;}
138 jerror_t PrintChiSqVector(
void)
const;
139 jerror_t
Print(
void)
const;
140 void FindSenseOfRotation(
void);
141 jerror_t Dump(
void)
const;
146 n.SetXYZ(N[0],N[1],N[2]);
173 jerror_t FillTrackParams(
void);
179 double xavg[3],var_avg;
184 #endif //_DHELICAL_FIT_H_
void SetMagneticFieldMap(const DMagneticFieldMap *map)
class DFDCPseudo: definition for a reconstructed point in the FDC
float chisq
chi-sq contribution of this hit
const DMagneticFieldMap * bfield
pointer to magnetic field map
ChiSqSourceType_t chisq_source
const DMagneticFieldMap * GetMagneticFieldMap() const
vector< DHFHit_t * > hits
const vector< DHFHit_t * > GetHits() const
float phi_circle
< error info in radial coordinates
float z
point in lab coordinates
double covr
< error info in RPhi coordinates
void GetPlaneParameters(double &c, DVector3 &n)