94 jout <<
"Printing Derivatives for Ring " <<
ring <<
" Straw " <<
straw << std::endl;
104 jout <<
" phi x,y,z " <<
phiX <<
" " <<
phiY <<
" " <<
phiZ << std::endl;
106 jout <<
"sin " <<
sin(
phiZ) <<
" cos " << cos(
phiZ) << std::endl;
107 jout <<
"udir_mag " <<
udir_mag << std::endl;
114 double sinPhiX=
sin(
phiX);
double cosPhiX=cos(
phiX);
115 double sinPhiY=
sin(
phiY);
double cosPhiY=cos(
phiY);
116 double sinPhiZ=
sin(
phiZ);
double cosPhiZ=cos(
phiZ);
119 double cosStereo=cos(
stereo);
121 double x0Shift =
x0+
r0*cosPhiStraw;
122 double y0Shift =
y0+
r0*sinPhiStraw;
128 +sinPhiX*(-
z0*cosPhiZ*sinPhiY+y0Shift*sinPhiZ);
130 +cosPhiY*y0Shift*sinPhiX-x0Shift*sinPhiY);
132 -cosPhiX*(cosPhiZ*y0Shift+
z0*sinPhiY*sinPhiZ);
141 +sinPhiY*sinPhiZ*(y0Shift*cosPhiX-
z0*sinPhiX);
144 +sinPhiX*(y0Shift*cosPhiZ*sinPhiY+
z0*sinPhiZ)
145 +cosPhiX*(
z0*cosPhiZ*sinPhiY-y0Shift*sinPhiZ);
165 +sinPhiZ*(cosPhiStraw*sinPhiX*tanStereo-cosPhiX));
168 -sinPhiZ*(cosPhiX*sinPhiY
169 +(cosPhiY*sinPhiStraw-cosPhiStraw*sinPhiX*sinPhiY)*tanStereo));
178 +cosPhiX*(cosPhiZ+cosPhiStraw*sinPhiY*sinPhiZ*tanStereo));
180 -(cosPhiStraw*cosPhiY*sinPhiX+sinPhiStraw*sinPhiY)*tanStereo);
182 +sinPhiX*(sinPhiZ-cosPhiStraw*cosPhiZ*sinPhiY*tanStereo)
183 +cosPhiX*(cosPhiZ*sinPhiY+cosPhiStraw*sinPhiZ*tanStereo));
193 +(cosPhiY*sinPhiStraw-cosPhiStraw*sinPhiX*sinPhiY)*tanStereo);
vector< double > GetCDCWireDerivatives()
vector< double > derivatives