*-- Author : Stephen Burke SUBROUTINE FKANAL(JPL,RES,CHISQ,NPS,NRS) *-----------------------------------------Updates 06/08/92------- **: FKANAL.......SB. Minor bug fixed. *-----------------------------------------Updates 24/01/92------- **: FKANAL 30205.SB. Digi with wrong drift sign counts as a bad hit. **: FKANAL 30205.SB. Phi histograms added. *-----------------------------------------Updates---------------- ********************************************************************** * * * Analyse the KF output * * * ********************************************************************** DOUBLE PRECISION RES(2),CHISQ,FKDPHI *KEEP,FKNPL. * * Per-track values can go in H1WORK; note that LTRUE and LFIRST must * be set at least per event. * * This is about 36k words long; the remaining common blocks are * about 3.6k in total. Some of this could be in /H1WORK/, but the * blocks would have to be reorganised. * * /FKPROJ/ * /FKFILT/ * /FKSMTH/ * /FKINT/ * /FKRSID/ * /FKTRUE/ * /FKDBG/ *KEEP,FFSTEE. *KEEP,FKCNTL. *KEEP,FKMEAS. *KEEP,FKTRUE. *KEEP,FKPROJ. *KEEP,FKSMTH. *KEEP,FKRSID. *KEEP,FKFLAG. *KEEP,FFDBG. *KEEP,FKHIST. *KEEP,FKDBG. *KEND. ********************************************************************** CALL HCDIR('//PAWC/'//CKDBG,' ') CALL HFILL(IBASE+1+5*IRP(JPL),SNGL(RES(1)),0.,1.) IF (MES(JPL).EQ.2) CALL HFILL(IBASE+12,SNGL(RES(2)),0.,1.) CALL HFILL(IBASE+5+5*IRP(JPL),SNGL(CHISQ),0.,1.) * Decode steering & CALL HFILL(101,SNGL((SSMT(1,JSTART)-TRUE(1,JSTART))/ & CALL HFILL(102,SNGL((SSMT(2,JSTART)-TRUE(2,JSTART))/ & CALL HFILL(103,SNGL((SSMT(3,JSTART)-TRUE(3,JSTART))/ & CALL HFILL(104,SNGL((SSMT(4,JSTART)-TRUE(4,JSTART))/ & CALL HFILL(105,SNGL((SSMT(5,JSTART)-TRUE(5,JSTART))/ C D1OP = DABS(SSMT(3,JSTART)) - DABS(TRUE(3,JSTART)) DPHI = FKDPHI(SSMT(5,JSTART),TRUE(5,JSTART)) CALL HFILL(230,FLOAT(NPS+NRS),D1OP,1.) CALL HFILL(231,FLOAT(NPS+NRS),D1OP,1.) CALL HFILL(232,FLOAT(NPS+NRS),D1OP,1.) CALL HFILL(233,PLOG,D1OP,1.) CALL HFILL(234,PLOG,DTHETA,1.) CALL HFILL(235,PLOG,DPHI,1.) CALL HFILL(243,THETA,D1OP,1.) CALL HFILL(244,THETA,DTHETA,1.) CALL HFILL(245,THETA,DPHI,1.) CALL HFILL(253,PHI,D1OP,1.) CALL HFILL(254,PHI,DTHETA,1.) CALL HFILL(255,PHI,DPHI,1.) CALL HFILL(IBASE+1+5*IRP(KPL),SNGL(RSMT(1,KPL)),0.,1.) & CALL HFILL(IBASE+12,SNGL(RSMT(2,KPL)),0.,1.) CALL HFILL(IBASE+5+5*IRP(KPL),SNGL(CHISMT(KPL)),0.,1.) CALL HFILL(KPL+500,CHPROB,0.,1.) CALL HFILL(KPL+600,CHPROB,0.,1.) CALL HFILL(KPL+700,CHPROB,0.,1.) CALL HFILL(KPL+800,SNGL(RSMT(1,KPL)),0.,1.) & CALL HFILL(KPL+900,SNGL(RSMT(2,KPL)),0.,1.) *