--- DarthVader/TrackerLevel2/src/F77/analysissubroutines.f 2009/08/04 14:01:35 1.40 +++ DarthVader/TrackerLevel2/src/F77/analysissubroutines.f 2009/11/10 15:50:16 1.42 @@ -303,7 +303,7 @@ * 1st) decreasing n.points * 2nd) increasing chi**2 * ------------------------------------------------------- - rchi2best=1000000000. + rchi2best=1000000000. ndofbest=0 do i=1,ntracks ndof=0 @@ -706,7 +706,8 @@ stripx = stripx + corr resxPAM = res - 10 endif + 10 continue + endif * ----------------- * CLUSTER Y @@ -752,7 +753,8 @@ stripy = stripy + corr resyPAM = res - 20 endif + 20 continue + endif c=========================================================== @@ -3882,7 +3884,20 @@ if(ip.eq.6.and.ssensor.ne.0)ssensor = 3 - ssensor !notazione paolo x align LS(IP,ntr) = ssensor+10*sladder - if(id.ne.0)then +c if(id.ne.0)then +CCCCCC(10 novembre 2009) PATCH X NUCLEI +C non ho capito perche', ma durante il ritracciamento dei nuclei +C (quando una traccia viene trovata ma non e' buona) c'e' qualche variabile +C che non viene reinizializzata correttamente e i cluster esclusi +C dal fit risultano ancora inclusi... +C + cltrx(ip,ntr) = 0 + cltry(ip,ntr) = 0 + if( + $ xgood_nt(ip,ntr).eq.1.and.ygood_nt(ip,ntr).eq.1 + $ .and. + $ id.ne.0)then + c >>> is a couple cltrx(ip,ntr) = clx(nplanes-ip+1,icp_cp(id)) cltry(ip,ntr) = cly(nplanes-ip+1,icp_cp(id))