*-- Author : Stephen Burke SUBROUTINE FFBKLK *-----------------------------------------Updates 22/11/93------- **: FFBKLK.......SB. Changes to monitoring histograms. *-----------------------------------------Updates 27/07/93------- **: FFBKLK 30907 SB. Changes to monitoring histograms. **: FFBKLK 30907 RP. Farm changes. *-----------------------------------------Updates 02/06/93------- **: FFBKLK 30907 SB. Correct bad update. *-----------------------------------------Updates 03/05/93------- **: FFBKLK 30907 SB. Change range of MC histos. *-----------------------------------------Updates 03/03/93------- **: FFBKLK 30907 SB. New monitoring histograms. *-----------------------------------------Updates 03/03/93------- **: FFBKLK 30907 SB. New monitoring histograms. *-----------------------------------------Updates 30/11/92------- **: FFBKLK 30907 SB. New monitoring histograms. *-----------------------------------------Updates 06/08/92------- **: FFBKLK 30205.SB. Remove vertex histos, add xy, change scales. *-----------------------------------------Updates 24/01/92------- **: FFBKLK 30205.SB. Add some new histograms, with new numbering. **: FFBKLK 30205.SB. Change range of histograms for data. *-----------------------------------------Updates---------------- ********************************************************************** * * * Book LOOK monitoring histograms for the Kalman filter * * * ********************************************************************** *KEEP,H1EVDT. * * IDATA type of information (HEAD bank word 6) : * * 0 - real data H1 * 1 - MC data H1SIM * 2 - real data CERN tests * 3 - MC data ARCET * * MONTE = .TRUE. if IDATA=1 * KEVENT = event processed counter for H1REC * *KEEP,BCS. *KEEP,BOSMDL. C ------BOSMDL C ------ *KEEP,FTHIST. * indices of filter farm histos *KEND. ********************************************************************** * * Histogram offset (for filter) * * * Residual histogram range for MC and data * C RANGE = 0.06 CALL BHS(IHX+1,0,60,1.,61.) CALL STEXT(IHX+1,4,'Tracks/event') CALL BHS(IHX+2,0,37,0.,37.) CALL STEXT(IHX+2,4,'Planar hits/track') CALL BHS(IHX+3,0,37,0.,37.) CALL STEXT(IHX+3,4,'Radial hits/track') CALL BHS(IHX+4,0,60,-RANGE,RANGE) CALL STEXT(IHX+4,4,'Drift residuals (planars)') CALL STEXT(IHX+4,1,'cm') CALL BHS(IHX+5,0,60,-RANGE,RANGE) CALL STEXT(IHX+5,4,'Drift residuals (radials)') CALL STEXT(IHX+5,1,'cm') CALL BHS(IHX+6,0,60,-12.0,12.0) CALL STEXT(IHX+6,4,'Radius residuals (radials)') CALL STEXT(IHX+6,1,'cm') CALL BHS(IHX+7,0,64,-3.2,3.2) CALL STEXT(IHX+7,4,'Phi') CALL BHS(IHX+8,0,60,0.,0.6) CALL STEXT(IHX+8,4,'Theta') CALL BHS(IHX+9,0,60,-0.3,0.3) CALL STEXT(IHX+9,4,'Q/Momentum') CALL STEXT(IHX+9,1,'GeV**-1') CALL BHS(IHX+10,0,60,-1.0,2.0) CALL STEXT(IHX+10,4,'LOG(Momentum)') CALL BHS(IHX+11,0,60,15.,75.) CALL STEXT(IHX+11,4,'R start') CALL STEXT(IHX+11,1,'cm') CALL BHS(IHX+12,0,64,-3.2,3.2) CALL STEXT(IHX+12,4,'PHI start') CALL BHS(IHX+13,0,60,0.,2400.) CALL STEXT(IHX+13,4,'Number of planar hits') CALL BHS(IHX+14,0,60,0.,2400.) CALL STEXT(IHX+14,4,'Number of radial hits') CALL BHS(IHX+15,0,50,0.,1.) CALL STEXT(IHX+15,4,'Fraction of planar hits used') CALL BHS(IHX+16,0,50,0.,1.) CALL STEXT(IHX+16,4,'Fraction of radial hits used') CALL BHS(IHX+17,0,37,0.,37.) CALL STEXT(IHX+17,4,'Planar hits per segment') CALL BHS(IHX+18,0,37,0.,37.) CALL STEXT(IHX+18,4,'Radial hits per segment') CALL BHS(IHX+19,0,64,0.,64.) CALL STEXT(IHX+19,4,'Supermodule mask') C CALL BHS(IHX+20,0,60,0.,240.) C CALL STEXT(IHX+20,4,'Chi-squared between start/end parameters') CALL BHS(IHX+20,0,64,0.,64.) CALL STEXT(IHX+20,4,'Secondary/tertiary mask') CALL BHS(IHX+51,0,60,0.,12.) CALL STEXT(IHX+51,4,'Chi-squared/DOF') CALL BHS(IHX+52,0,50,0.,1.) CALL STEXT(IHX+52,4,'Chi-squared probability') CALL BHS(IHX+53,0,60,0.,0.6) CALL STEXT(IHX+53,4,'Fractional error on 1/r') CALL BHS(IHX+54,0,60,0.,0.012) CALL STEXT(IHX+54,4,'Error on theta') CALL BHS(IHX+55,0,60,0.,0.15) CALL STEXT(IHX+55,4,'Error on xy') CALL STEXT(IHX+55,1,'cm') CALL BHD(IHX+71,0,32,0.,32.,36,0.,36.) CALL STEXT(IHX+71,4,'Planar used hit map (.gt.1 planar segment)') CALL BHD(IHX+72,0,48,0.,48.,36,0.,36.) CALL STEXT(IHX+72,4,'Radial used hit map (.gt.1 planar segment)') CALL BHD(IHX+73,0,32,0.,32.,36,0.,36.) CALL STEXT(IHX+73,4,'Planar unused hit map') CALL BHD(IHX+74,0,48,0.,48.,36,0.,36.) CALL STEXT(IHX+74,4,'Radial unused hit map') CALL BHS(IHX+75,0,60,-RANGE,RANGE) CALL STEXT(IHX+75,4,'Drift residuals (planars) (>1 seg)') CALL STEXT(IHX+75,1,'cm') CALL BHS(IHX+76,0,60,-RANGE,RANGE) CALL STEXT(IHX+76,4,'Drift residuals (radials) (>1 seg)') CALL STEXT(IHX+76,1,'cm') CALL BHS(IHX+77,0,60,-3.0,3.0) CALL STEXT(IHX+77,4,'Planar drift distance') CALL STEXT(IHX+77,1,'cm') CALL BHS(IHX+78,0,60,-6.0,6.0) CALL STEXT(IHX+78,4,'Radial drift distance') CALL STEXT(IHX+78,1,'cm') CALL BHS(IHX+79,0,60,-0.15,0.15) CALL STEXT(IHX+79,4,'Planar C1 checksum') CALL STEXT(IHX+79,1,'cm') CALL BHS(IHX+80,0,60,-0.15,0.15) CALL STEXT(IHX+80,4,'Planar C2 checksum') CALL STEXT(IHX+80,1,'cm') CALL BHS(IHX+81,0,60,-0.15,0.15) CALL STEXT(IHX+81,4,'Radial C1 checksum') CALL STEXT(IHX+81,1,'cm') CALL BHS(IHX+82,0,60,-0.15,0.15) CALL STEXT(IHX+82,4,'Radial C2 checksum') CALL STEXT(IHX+82,1,'cm') CALL BHS(IHX+83,0,60,-RANGE,RANGE) CALL STEXT(IHX+83,4,'Drift residuals (planars) (<5 mm)') CALL STEXT(IHX+83,1,'cm') CALL BHS(IHX+84,0,60,-RANGE,RANGE) CALL STEXT(IHX+84,4,'Drift residuals (radials) (<5 mm)') CALL STEXT(IHX+84,1,'cm') CALL BHD(IHX+85,0,32,0.,32.,36,0.,36.) CALL STEXT(IHX+85,4,'Planar used hit map') CALL BHD(IHX+86,0,48,0.,48.,36,0.,36.) CALL STEXT(IHX+86,4,'Radial used hit map') CALL BHS(IHX+87,0,60,-0.15,0.15) CALL STEXT(IHX+87,4,'Planar C2 checksum (crossing wire plane)') CALL STEXT(IHX+87,1,'cm') CALL BHS(IHX+88,0,60,-0.15,0.15) CALL STEXT(IHX+88,4,'Radial C2 checksum (crossing wire plane)') CALL STEXT(IHX+88,1,'cm') CALL BHS(IHX+89,0,60,-0.15,0.15) CALL STEXT(IHX+89,4,'Planar C2 checksum (crossing cathode)') CALL STEXT(IHX+89,1,'cm') CALL BHS(IHX+90,0,60,-0.15,0.15) CALL STEXT(IHX+90,4,'Radial C2 checksum (crossing cathode)') CALL STEXT(IHX+90,1,'cm') CALL BHD(IHX+68,0,60,-RANGE,RANGE,50,0.,2500.) CALL STEXT(IHX+68,4,'Drift residuals (planars) vs Q') CALL BHD(IHX+69,0,60,-RANGE,RANGE,50,0.,2500.) CALL STEXT(IHX+69,4,'Drift residuals (radials) vs Q') CALL BHS(IHX+70,0,60,-12.0,12.0,50,0.,2500.) CALL STEXT(IHX+70,4,'Radius residuals (radials) vs Q') * * Next lot for debug only * CALL BHS(IHX+56,0,65,0.,130.) CALL STEXT(IHX+56,4,'Track length in z') CALL BHS(IHX+57,0,60,0.,300.) CALL STEXT(IHX+57,4,'Number of planar segments') CALL BHS(IHX+58,0,60,0.,300.) CALL STEXT(IHX+58,4,'Number of radial segments') CALL BHD(IHX+59,0,40,-80.,80.,40,-80.,80.) CALL STEXT(IHX+59,4,'x-y position of starting point') CALL BHD(IHX+60,0,40,-80.,80.,40,-80.,80.) CALL STEXT(IHX+60,4,'x-y position of end point') CALL BHS(IHX+61,0,60,15.,75.) CALL STEXT(IHX+61,4,'R end') CALL BHS(IHX+62,0,64,-3.2,3.2) CALL STEXT(IHX+62,4,'phi end') CALL BHD(IHX+63,0,37,0.,37.,37,0.,37.) CALL STEXT(IHX+63,4,'Radial hits/track vs planar hits/track') CALL BHD(IHX+64,0,50,0.,1.0,73,0.,73.) CALL STEXT(IHX+64,4,'Hits/track vs theta') CALL BHD(IHX+65,0,60,-RANGE,RANGE,64,0.,64.) CALL STEXT(IHX+65,4,'Drift residuals (planars) vs IMAP') CALL BHD(IHX+66,0,60,-RANGE,RANGE,64,0.,64.) CALL STEXT(IHX+66,4,'Drift residuals (radials) vs IMAP') CALL BHS(IHX+67,0,60,-12.0,12.0,64,0.,64.) CALL STEXT(IHX+67,4,'Radius residuals (radials) vs IMAP') *