Function for change in x', Lab frame

Function for the change in x' in the detector frame for change in [math]\phi[/math] and constant [math]\theta[/math] in the lab frame

[math]D2P=\sqrt{(x_{D2}-x_P)^2+(y_{D2}-y_P)^2+(z_{D2}-z_P)^2}[/math]

[math]D1P=\sqrt{(x_P-x_{D1})^2+(y_P-y_{D1})^2+(z_P-z_{D1})^2}[/math]

[math]x_1^'=\frac{((x_{D2}-x_P)^2+(y_{D2}-y_P)^2+(z_{D2}-z_P)^2)-((x_P-x_{D1})^2+(y_P-y_{D1})^2+(z_P-z_{D1})^2)}{4ae}-ae[/math]

[math]x_{D1}=r_{D1}\ cos(\phi)\qquad y_{D1}=r_{D1}cos(\phi)\qquad z_{D1}=r_{D1} cot(\theta)[/math]

[math]x_{D2}=r_{D2} cos(\phi)\qquad y_{D2}=r_{D2} sin(\phi)\qquad z_{D2}=r_{D2} cot(\theta)[/math]

[math]x_P=\frac{2.53cos(\phi)}{(cot(\theta)+cos(\phi)cot(65^{\circ})}[/math]

[math]y_P=\frac{2.53sin(\phi)}{(cot(\theta)+cos(\phi)cot(65^{\circ})}[/math]

[math]z_P=\frac{2.53cot(\theta)}{(cot(\theta)+cos(\phi)cot(65^{\circ})}[/math]

[math]x_1^'=\frac{((x_{D2}-x_P)^2+(y_{D2}-y_P)^2+(z_{D2}-z_P)^2)-((x_P-x_{D1})^2+(y_P-y_{D1})^2+(z_P-z_{D1})^2)}{4ae}-ae[/math]

[math]x_1^'=\frac{x_{D2}^2-2x_Px_{D2}+x_P^2+y_{D2}^2-2y_Py_{D2}+y_P^2+z_{D2}^2-2z_Pz_{D2}+z_P^2-x_P^2+2x_Px_{D1}-x_{D1}^2-y_P^2+2y_Py_{D1}-y_{D1}^2-z_P^2+2z_Pz_{D1}-z_{D1}^2}{4ae}-ae[/math]

[math]x_1^'=\frac{(x_{D2}^2+y_{D2}^2)-(x_{D1}^2+y_{D1}^2)+z_{D2}^2-z_{D1}^2-2x_P(x_{D2}-x_{D1})-2y_P(y_{D2}-y_{D1})-2z_P(z_{D2}-z_{D1})}{4ae}-ae[/math]

[math]x_1^'=\frac{(r_{D2}^2)-(r_{D1}^2)+cot^2(\theta)(r_{D2}^2-r_{D1}^2)-2x_P(x_{D2}-x_{D1})-2y_P(y_{D2}-y_{D1})-2z_P(z_{D2}-z_{D1})}{4ae}-ae[/math]

Expressing this as functions of [math]\phi[/math] and non-differentiable constants

[math]x_1^'=\frac{c_1+c_2-2x_P(\phi)x_{D2}(\phi)+2x_P(\phi)x_{D1}(\phi)-2y_P(\phi)y_{D2}(\phi)+2y_P(\phi)y_{D1}(\phi)-2z_P(\phi)c_3}{4c_4}-c_4[/math]

Differentiating with respect to [math]\phi[/math]

[math]x_{D1}=r_{D1} cos(\phi)\Rightarrow \dot x_{D1}=-r_{D1} sin(\phi)[/math]

[math]y_{D1}=r_{D1}sin(\phi)\Rightarrow \dot y_{D1}=r_{D1}cos(\phi)[/math]

[math]x_{D2}=r_{D2} cos(\phi)\Rightarrow \dot x_{D2}=-r_{D2} sin(\phi)[/math]

[math]y_{D2}=r_{D2}sin(\phi)\Rightarrow \dot y_{D2}=r_{D2}cos(\phi)[/math]

[math]x_P=\frac{2.52934271645cos(\phi)}{cot(\theta)+cos(\phi)cot(65^{\circ})}\Rightarrow \dot x_P=\frac{-2.52934271645cot(\theta)sin(\phi)}{(cos(\phi)cot(65^{\circ}+cot(\theta))^2}[/math]

[math]y_P=\frac{2.52934271645sin(\phi)}{cot(\theta)+cos(\phi)cot(65^{\circ})}\Rightarrow \dot y_P=\frac{-1.7206+2.52934271645 cos(\phi) cot(\theta)}{(cos(\phi) cot(65^{\circ}) + cot(\theta))^2}[/math]

[math]z_P=\frac{2.52934271645cot(\theta)}{cot(\theta)+cos(\phi)cot(65^{\circ})}\Rightarrow \dot z_P=\frac{-1.7206 cot(\theta)sin(\phi))}{(cos(\phi) cot(65) + cot(\theta))^2}[/math]

[math]\frac{dx_1^1}{d\phi}=\frac{-2}{4c_4}\frac{d}{d\phi}(x_P(\phi)x_{D2}(\phi))+\frac{2}{4c_4}\frac{d}{d\phi}(x_P(\phi)x_{D1}(\phi))-\frac{2}{4c_4}\frac{d}{d\phi}(y_P(\phi)y_{D2}(\phi))+\frac{2}{4c_4}\frac{d}{d\phi}(y_P(\phi)y_{D1}(\phi))-\frac{2c_3}{4c_4}\frac{d}{d\phi}z_P(\phi)[/math]

[math]\frac{dx_1^1}{d\phi}=\frac{-2}{4c_4} \left ( (\dot x_P(\phi)x_{D2}(\phi)+x_P(\phi)\dot x_{D2}(\phi))-(\dot x_P(\phi)x_{D1}(\phi)+x_P(\phi)\dot x_{D1}(\phi))+(\dot y_P(\phi)y_{D2}(\phi)+y_P(\phi)\dot y_{D2}(\phi))-(\dot y_P(\phi)y_{D1}(\phi)+y_P(\phi)\dot y_{D1}(\phi))+c_3\dot z_P(\phi) \right )[/math]