| 2042 |
* tg(th_yz) |
* tg(th_yz) |
| 2043 |
alfayz2(ndblt)=(ym1-ym2)/(zm1-zm2) |
alfayz2(ndblt)=(ym1-ym2)/(zm1-zm2) |
| 2044 |
* y0 (cm) |
* y0 (cm) |
| 2045 |
alfayz1(ndblt)=alfayz2(ndblt)*(REAL(zini)-zm1)+ym1 ! EM GCC4.7 zm1, ym1 and alfayz1/2 are REAL |
alfayz1(ndblt)=alfayz2(ndblt)*(REAL(zini)-zm1)+ym1! EM GCC4.7 zm1, ym1 and alfayz1/2 are REAL |
| 2046 |
|
|
| 2047 |
**** -----------------------------------------------**** |
**** -----------------------------------------------**** |
| 2048 |
**** reject non phisical couples **** |
**** reject non phisical couples **** |
| 3399 |
dedxmmx = sgnl(icx)/mip(VIEW(icx),LADDER(icx)) !(1)(2) |
dedxmmx = sgnl(icx)/mip(VIEW(icx),LADDER(icx)) !(1)(2) |
| 3400 |
dedxmmy = sgnl(icy)/mip(VIEW(icy),LADDER(icy)) !(1)(2) |
dedxmmy = sgnl(icy)/mip(VIEW(icy),LADDER(icy)) !(1)(2) |
| 3401 |
clincnewc=10.*dsqrt(rymm**2+rxmm**2 |
clincnewc=10.*dsqrt(rymm**2+rxmm**2 |
| 3402 |
$ +DBLE(RCHI2_STORE(ibest)*k(ip)*(cov(1,1)+cov(2,2)))) ! EM GCC4.7 |
$ +DBLE(RCHI2_STORE(ibest)*k(ip)*(cov(1,1)+cov(2,2))))! EM GCC4.7 |
| 3403 |
endif |
endif |
| 3404 |
1188 continue |
1188 continue |
| 3405 |
enddo !end loop on couples on plane icp |
enddo !end loop on couples on plane icp |