21 vector<const DFDCIntersection*> fdcintersections;
22 loop->Get(fdcintersections);
24 for(
unsigned int i=0; i<fdcintersections.size(); i++){
28 MakePseudo(fdcintersection->
hit1, fdcintersection->
wire1, fdcintersection->
pos);
29 MakePseudo(fdcintersection->
hit2, fdcintersection->
wire2, fdcintersection->
pos);
48 pseudo->
dw = 1.0/
sqrt(12.0);
49 pseudo->
s = (R-a)*udir;
50 pseudo->
ds = 1.0/
sqrt(12.0);
52 pseudo->
time = hit->
t;
54 pseudo->
xy.Set(pos.X(),pos.Y());
58 cov(0,0) = pseudo->
dw*pseudo->
dw; cov(1,0) = 0.0;
59 cov(0,1) = 0.0; cov(1,1) = pseudo->
ds*pseudo->
ds;
63 double cos_angle = cos(wire->
angle);
65 Rot(0,0) = cos_angle; Rot(1,0) = sin_angle;
66 Rot(0,1) = -sin_angle; Rot(1,1) = cos_angle;
67 DMatrix RotT(DMatrix::kTransposed, Rot);
70 pseudo->
covxx=RotCov(0,0);
71 pseudo->
covxy=RotCov(1,0);
72 pseudo->
covyy=RotCov(1,1);
74 _data.push_back(pseudo);
DVector2 xy
rough x,y coordinates in lab coordinate system
int status
status word for pseudopoint
const DFDCWire * wire
DFDCWire for this wire.
class DFDCPseudo: definition for a reconstructed point in the FDC
double ds
local coordinate of pseudopoint in the direction along the wire and its uncertainty ...
double s
< wire position computed from cathode data, assuming the avalanche occurs at the wire ...
double covyy
Covariance terms for (x,y)
double time
time corresponding to this pseudopoint.
jerror_t evnt(JEventLoop *eventLoop, uint64_t eventnumber)
Called every event.
void MakePseudo(const DFDCHit *hit, const DFDCWire *wire, const DVector3 &pos)
double dw
local coordinate of pseudopoint in the direction perpendicular to the wires and its uncertainty ...
class DFDCHit: definition for a basic FDC hit data type.