version 1.4, 1994/09/01 12:29:07
|
version 1.5, 1994/10/12 18:52:06
|
|
|
* 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.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 |
* | * |
|
|
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 |
* | * |