*-- Author : Girish D. Patel 07/12/93 SUBROUTINE FCHKQR ********************************************************************** * * * Fill histograms from hits in COMMON/H1WORK/, entered by FILLQR * * * * JVM 22/7/92 * ********************************************************************** *KEEP,FMOBIN. *KEEP,FMOWRK. * planar hit data... * radial hit data... *KEND. *KEEP,FMOHIS. *KEND. * check that COMMON/H1WORK/ has some data .. * cycle through the cells ... CALL SHS (IHS(12),0,CHITS) !SGI * analyse existing hits in cell NC, if there are any CALL SHD (IHS(19),0,TDIF,0.) !SGI CALL SHD (IHS(19),0,TDIF,FLOAT(NSMOD+1)) !SGI * * complete set of hits available here * set flag "LFREE" to flag Qs which have not yet been plotted in * hist ID=68 ... CALL SHD (IHS(23),0,X2,13.5) !SGI CALL SHD (IHS(23),0,X2,FLOAT(NSMOD+1)+13.5) !SGI CALL SHD (IHS(23),0,X2,FLOAT(NLAY+4)+13.5) !SGI CALL SHD (IHS(23),0,X3,13.5) !SGI CALL SHD (IHS(23),0,X3,FLOAT(NSMOD+1)+13.5) !SGI CALL SHD (IHS(23),0,X3,FLOAT(NLAY+4)+13.5) !SGI CALL SHD (IHS(23),0,X1,0.5) !SGI CALL SHD (IHS(23),0,X1,FLOAT(NSMOD+1)+0.5) !SGI CALL SHD (IHS(23),0,X1,FLOAT(NLAY+4)+0.5) !SGI * RAD = (RR(NC,2,J1)+RR(NC,1,J2)+RR(NC,4,J3)+RR(NC,3,J4))/4.0 * TTW(1) = TTR(NC,1,J1) * TTW(2) = TTR(NC,2,J2) * TTW(3) = TTR(NC,3,J3) * TTW(4) = TTR(NC,4,J4) * RRW(1) = RR(NC,1,J1) * RRW(2) = RR(NC,2,J2) * RRW(3) = RR(NC,3,J3) * RRW(4) = RR(NC,4,J4) * SRZ = -0.3*RRW(1) - 0.1*RRW(2) + 0.1*RRW(3) + 0.3*RRW(4) * CRZ = RRW(1) + 0.5*RRW(2) - 0.5*RRW(4) CALL SHD (IHS(22),0,QQWR(J),0.) !SGI CALL SHD (IHS(22),0,QQWR(J),FLOAT(NSMOD+1)) !SGI * RPRED = SRZ*FLOAT(J) + CRZ * IF(RPRED.LT.80.0) THEN * YRPRED = RPRED*FLOAT(NBINR)/80.0 * IRPRED = YRPRED * IF(X1.GT.0.0) THEN * CALL SHDW(IHS(17),0,TTW(J)-60.,YRPRED,-1.) !SGI * CALL SHDW(IHS(17),0,TTW(J)+60.,YRPRED,+1.) !SGI * ELSE * CALL SHDW(IHS(18),0,TTW(J)-60.,YRPRED,-1.) !SGI * CALL SHDW(IHS(18),0,TTW(J)+60.,YRPRED,+1.) !SGI * ENDIF * ENDIF * IF(RAD.LT.80.0.AND.ABS(X1).GT.8.0.AND.ABS(X1).LT.32.0) THEN * ENDIF CALL SHD (IHS(23),0,X2,26.5) !SGI CALL SHD (IHS(23),0,X2,FLOAT(NSMOD+1)+26.5) !SGI CALL SHD (IHS(23),0,X2,FLOAT(NLAY+4)+26.5) !SGI CALL SHD (IHS(23),0,X3,26.5) !SGI CALL SHD (IHS(23),0,X3,FLOAT(NSMOD+1)+26.5) !SGI CALL SHD (IHS(23),0,X3,FLOAT(NLAY+4)+26.5) !SGI * zero arrays ready for next event CALL VZERO( TTR(0,1,1),NW3 ) CALL VZERO( QQR(0,1,1),NW3 ) CALL VZERO( NHITR(0,1),NW2 ) *