|
version 1.4, 1994/09/01 12:29:07
|
version 1.7, 1995/01/27 20:26:50
|
|
|
|
| * remove minuit. Make fit linear | * remove minuit. Make fit linear |
| * still does not do errors properly | * still does not do errors properly |
| * $Log$ | * $Log$ |
| |
* Revision 1.7 1995/01/27 20:26:50 cdaq |
| |
* (JRA) Remove Mack's personal focalplane diceamatic (z slicer) code |
| |
* |
| |
* Revision 1.6 1994/12/06 15:45:27 cdaq |
| |
* (DJM) Take slices in Z to look for best focus |
| |
* |
| |
* Revision 1.5 1994/10/12 18:52:06 cdaq |
| |
* (DJM) Initialize some variables |
| |
* (SAW) Prettify indentation |
| |
* |
| * Revision 1.4 1994/09/01 12:29:07 cdaq | * Revision 1.4 1994/09/01 12:29:07 cdaq |
| * (DJM) Make registered versions of residuals | * (DJM) Make registered versions of residuals |
| * | * |
|
|
|
| integer*4 itrack ! track loop index | integer*4 itrack ! track loop index |
| integer*4 ihit,ierr | integer*4 ihit,ierr |
| integer*4 hit,plane | integer*4 hit,plane |
| integer*4 i,j ! loop index |
integer*4 i,j,k ! loop index |
| |
* real*4 z_slice |
| |
|
| real*8 H_DPSIFUN | real*8 H_DPSIFUN |
| real*8 pos | real*8 pos |
| real*4 ray(hnum_fpray_param) | real*4 ray(hnum_fpray_param) |
|
|
|
| do plane=1,HMAX_NUM_DC_PLANES | do plane=1,HMAX_NUM_DC_PLANES |
| hdc_residual(itrack,plane)=1000 | hdc_residual(itrack,plane)=1000 |
| hdc_sing_res(itrack,plane)=1000 | hdc_sing_res(itrack,plane)=1000 |
| |
hdc1_sing_res(plane)=1000 |
| |
hdc2_sing_res(plane)=1000 |
| |
hdc1_dbl_res(plane)=1000 |
| |
hdc2_dbl_res(plane)=1000 |
| |
|
| enddo | enddo |
| enddo | enddo |
| | |
|
|
|
| do i=1,hnum_fpray_param | do i=1,hnum_fpray_param |
| do j=1,hnum_fpray_param | do j=1,hnum_fpray_param |
| AA(i,j)=0. | AA(i,j)=0. |
| |
if(j.lt.i)then |
| |
AA(i,j)=AA(j,i) |
| |
else |
| do ihit=2,HNTRACK_HITS(itrack,1)+1 | do ihit=2,HNTRACK_HITS(itrack,1)+1 |
| hit=HNTRACK_HITS(itrack,ihit) | hit=HNTRACK_HITS(itrack,ihit) |
| plane=HDC_PLANE_NUM(hit) | plane=HDC_PLANE_NUM(hit) |
| AA(i,j)=AA(i,j) + ( | AA(i,j)=AA(i,j) + ( |
| & hplane_coeff(remap(i),plane)*hplane_coeff(remap(j),plane) |
& hplane_coeff(remap(i),plane)*hplane_coeff(remap(j) |
| & /(hdc_sigma(plane)*hdc_sigma(plane))) |
$ ,plane)/(hdc_sigma(plane)*hdc_sigma(plane))) |
| enddo ! end loop on ihit | enddo ! end loop on ihit |
| |
endif ! end test on j .lt. i |
| enddo ! end loop on j | enddo ! end loop on j |
| enddo ! end loop on i | enddo ! end loop on i |
| * | * |
|
|
|
| enddo ! end loop over tracks | enddo ! end loop over tracks |
| endif | endif |
| | |
| |
** A reasonable selection of slices is presently -80,-60,-40,-20,0,20,40 |
| |
** ,60,80. Zero is the nominal midplane between the chambers, -80 |
| |
** corresponds closely to the exit flange position. This slice pattern |
| |
** is created with hz_wild=-80., hdelta_z_wild=20., and hnum_zslice=9 |
| |
* if(hz_slice_enable.ne.0)then |
| |
* do k=1,hnum_zslice |
| |
* z_slice = hz_wild + (k-1)*hdelta_z_wild |
| |
* hx_fp_wild(k) = hx_fp(1) + hxp_fp(1)*z_slice |
| |
* hy_fp_wild(k) = hy_fp(1) + hyp_fp(1)*z_slice |
| |
* enddo |
| |
* endif |
| |
|
| * calculate residuals for each chamber if in single stub mode | * calculate residuals for each chamber if in single stub mode |
| * and there were 2 tracks found one in first chanber and one in the second | * and there were 2 tracks found one in first chanber and one in the second |
| | |
Return to
h_track_fit.f
CVS log
Up to [HallC] / Analyzer / HTRACKING