C opens an rz file and books histograms C*************************************************************************** subroutine hist_book implicit none include 'cwn5.inc' c include 'ghib$data:fidfiles.inc' integer lrecl, ierr, lun, icycle, nid,i,j,k character*64 chdir, chopt, chfile, chtop, chpath,fname REAL GENMAX(21),GENMIN(21),JETSMAX(0:20,31),JETSMIN(0:20,31) REAL JT1SMAX(0:20,12),JT1SMIN(0:20,12),METMAX(19),METMIN(19) REAL ELEMAX(0:10,37),ELEMIN(0:10,37) REAL MUMAX(0:10,31),MUMIN(0:10,31) REAL TAUMAX(0:10,19),TAUMIN(0:10,19),DILMAX(13),DILMIN(13) REAL PHOMAX(0:10,62),PHOMIN(0:10,62) REAL TAGMAX(0:4,6),TAGMIN(0:4,6) REAL VERMAX(0:20,8),VERMIN(0:20,8) REAL TRKMAX(0:100,24),TRKMIN(0:100,24) COMMON /LIMITS1/ GENMAX,GENMIN,JETSMAX,JETSMIN,JT1SMAX,JT1SMIN COMMON /LIMITS2/ METMAX,METMIN,ELEMAX,ELEMIN,MUMAX,MUMIN COMMON /LIMITS3/ TAUMAX,TAUMIN,DILMAX,DILMIN,PHOMAX,PHOMIN COMMON /LIMITS4/ TAGMAX,TAGMIN,VERMAX,VERMIN COMMON /LIMITS5/ TRKMAX,TRKMIN INTEGER IQUEST COMMON/QUEST/IQUEST(100) INTEGER NWDS_HBOOK PARAMETER (NWDS_HBOOK=6000000) REAL HMEM COMMON /PAWC/ HMEM(NWDS_HBOOK) C create rz file IQUEST(10) = 65536 LRECL = 4095 LUN = 3 CHDIR = 'hist' CHFILE='pg_sampl.hist' CHOPT = 'NQ' CALL HROPEN(LUN,CHDIR,CHFILE,CHOPT,LRECL,IERR) if (ierr.eq.0) then CHPATH = '//hist/stntuple' CHOPT = ' ' CALL HMDIR(CHPATH,CHOPT) CALL HCDIR(CHPATH,CHOPT) else write (6,*) 'File ',CHFILE,' not opened.' STOP endif C book histos C call hbook1(id,chtitl,nx,xmi,vmx) C id histogram identifier C chtitl histogram title C nx number of channels C xmi lower edge of first channel C xma upper edge of last channel C vmx upper limit of single channel content C General Variable Limits C call hbook1(10 ,'Run ',88,Genmin(2),Genmax(2),0.) call hbook1(11 ,'Event ',88,Genmin(3),Genmax(3),0.) call hbook1(12 ,'ILUM ',88,Genmin(4),Genmax(4),0.) call hbook1(13 ,'Trigbits ',88,Genmin(5),Genmax(5),0.) c call hbook1(14 ,'Goodrun ',88,Genmin(6),Genmax(6),0.) call hbook1(15 ,'Brcode ',88,Genmin(7),Genmax(7),0.) c call hbook1(16 ,'Goodtrig ',88,Genmin(8),Genmax(8),0.) call hbook1(17 ,'Trigword ',88,Genmin(9),Genmax(9),0.) call hbook1(18 ,'Njet10 ',88,Genmin(10),Genmax(10),0.) call hbook1(19 ,'Njet15 ',88,Genmin(11),Genmax(11),0.) call hbook1(20 ,'Njet20 ',88,Genmin(12),Genmax(12),0.) call hbook1(21 ,'Njet25 ',88,Genmin(13),Genmax(13),0.) call hbook1(22 ,'Nhpte ',88,Genmin(14),Genmax(14),0.) call hbook1(23 ,'Nlcent ',88,Genmin(15),Genmax(15),0.) call hbook1(24 ,'Nlplug ',88,Genmin(16),Genmax(16),0.) call hbook1(25 ,'Nhpte ',88,Genmin(17),Genmax(17),0.) call hbook1(26 ,'Ncmio ',88,Genmin(18),Genmax(18),0.) call hbook1(27 ,'Ngtau ',88,Genmin(19),Genmax(19),0.) call hbook1(28 ,'Nlepton ',88,Genmin(20),Genmax(20),0.) call hbook1(29 ,'Etout ',88,Genmin(21),Genmax(21),0.) C INITIALIZE THE JET PLOTS call hbook1(35 ,'Njt ',20,Jetsmin(1,1),Jetsmax(1,1),0.) do i=1,20 call hbook1(35+30*(i-1)+1,'JTET',88,Jetsmin(i,2),Jetsmax(i,2),0.) call hbook1(35+30*(i-1)+2,'JTM',88,Jetsmin(i,3),Jetsmax(i,3),0.) call hbook1(35+30*(i-1)+3,'JTP4 x',88,Jetsmin(i,4),Jetsmax(i,4),0.) call hbook1(35+30*(i-1)+4,'JTEVETA',88,Jetsmin(i,5),Jetsmax(i,5),0.) call hbook1(35+30*(i-1)+5,'JTDTETA',88,Jetsmin(i,6),Jetsmax(i,6),0.) call hbook1(35+30*(i-1)+6,'JTPHI',88,Jetsmin(i,7),Jetsmax(i,7),0.) call hbook1(35+30*(i-1)+7,'JTEMF',88,Jetsmin(i,8),Jetsmax(i,8),0.) call hbook1(35+30*(i-1)+8,'JTWORD',88,Jetsmin(i,9),Jetsmax(i,9),0.) call hbook1(35+30*(i-1)+9,'JTCORFM',88,Jetsmin(i,10),Jetsmax(i,10),0.) call hbook1(35+30*(i-1)+10,'JTCRFD',88,Jetsmin(i,11),Jetsmax(i,11),0.) call hbook1(35+30*(i-1)+11,'JTCRFA',88,Jetsmin(i,12),Jetsmax(i,12),0.) call hbook1(35+30*(i-1)+12,'JTNTR',88,Jetsmin(i,13),Jetsmax(i,13),0.) call hbook1(35+30*(i-1)+13,'JTNTRP',88,Jetsmin(i,14),Jetsmax(i,14),0.) call hbook1(35+30*(i-1)+14,'JTEP',88,Jetsmin(i,15),Jetsmax(i,15),0.) call hbook1(35+30*(i-1)+15,'JTNJTR',88,Jetsmin(i,16),Jetsmax(i,16),0.) call hbook1(35+30*(i-1)+16,'JTTAG',88,Jetsmin(i,17),Jetsmax(i,17),0.) call hbook1(35+30*(i-1)+17,'JTEGRD',88,Jetsmin(i,18),Jetsmax(i,18),0.) call hbook1(35+30*(i-1)+18,'JVBPB',88,Jetsmin(i,19),Jetsmax(i,19),0.) call hbook1(35+30*(i-1)+19,'JPBPB',88,Jetsmin(i,20),Jetsmax(i,20),0.) call hbook1(35+30*(i-1)+20,'SCBPB',88,Jetsmin(i,21),Jetsmax(i,21),0.) call hbook1(35+30*(i-1)+21,'SCNPB',88,Jetsmin(i,22),Jetsmax(i,22),0.) call hbook1(35+30*(i-1)+22,'JPPPB',88,Jetsmin(i,23),Jetsmax(i,23),0.) call hbook1(35+30*(i-1)+23,'JPNPB',88,Jetsmin(i,24),Jetsmax(i,24),0.) call hbook1(35+30*(i-1)+24,'JTAU1',88,Jetsmin(i,25),Jetsmax(i,25),0.) call hbook1(35+30*(i-1)+25,'JTAU2',88,Jetsmin(i,26),Jetsmax(i,26),0.) call hbook1(35+30*(i-1)+26,'JTZVTX',88,Jetsmin(i,27),Jetsmax(i,27),0.) call hbook1(35+30*(i-1)+27,'JTDZVX',88,Jetsmin(i,28),Jetsmax(i,28),0.) call hbook1(35+30*(i-1)+28,'JT_POU',88,Jetsmin(i,29),Jetsmax(i,29),0.) call hbook1(35+30*(i-1)+29,'JPVPT',88,Jetsmin(i,30),Jetsmax(i,30),0.) call hbook1(35+30*(i-1)+30,'JNOVPT',88,Jetsmin(i,31),Jetsmax(i,31),0.) ENDDO C INITIALIZE THE JT1S VARIABLES C call hbook1(640 ,'Nj1 ',20,Jt1smin(1,1),Jt1smax(1,1),0.) do i=1,20 call hbook1(640+11*(i-1)+1,'J1ET',88,Jt1smin(i,2),Jt1smax(i,2),0.) call hbook1(640+11*(i-1)+2,'J1M',88,Jt1smin(i,3),Jt1smax(i,3),0.) call hbook1(640+11*(i-1)+3,'J1P4 x',88,Jt1smin(i,4),Jt1smax(i,4),0.) call hbook1(640+11*(i-1)+4,'J1EVETA',88,Jt1smin(i,5),Jt1smax(i,5),0.) call hbook1(640+11*(i-1)+5,'J1DTETA',88,Jt1smin(i,6),Jt1smax(i,6),0.) call hbook1(640+11*(i-1)+6,'J1PHI',88,Jt1smin(i,7),Jt1smax(i,7),0.) call hbook1(640+11*(i-1)+7,'J1WORD',88,Jt1smin(i,8),Jt1smax(i,8),0.) call hbook1(640+11*(i-1)+8,'J1JTS',88,Jt1smin(i,9),Jt1smax(i,9),0.) call hbook1(640+11*(i-1)+9,'J1CRFM',88,Jt1smin(i,10),Jt1smax(i,10),0.) call hbook1(640+11*(i-1)+10,'J1CRF',88,Jt1smin(i,11),Jt1smax(i,11),0.) call hbook1(640+11*(i-1)+11,'J1EGR',88,Jt1smin(i,12),Jt1smax(i,12),0.) Enddo C INITIALIZE THE MET PLOTS C call hbook1(865 ,'MET(5) ',88,Metmin(1),Metmax(1),0.) call hbook1(866 ,'METPHI(5) ',88,Metmin(2),Metmax(2),0.) call hbook1(867 ,'METJT ',88,Metmin(3),Metmax(3),0.) call hbook1(868 ,'METJTPHI ',88,Metmin(4),Metmax(4),0.) call hbook1(869 ,'METJ1 ',88,Metmin(5),Metmax(5),0.) call hbook1(870 ,'METJ1PHI ',88,Metmin(6),Metmax(6),0.) call hbook1(871 ,'METC ',88,Metmin(7),Metmax(7),0.) call hbook1(872 ,'METCPHI ',88,Metmin(8),Metmax(8),0.) call hbook1(873 ,'HTC ',88,Metmin(9),Metmax(9),0.) call hbook1(874 ,'METCF ',88,Metmin(10),Metmax(10),0.) call hbook1(875 ,'METCFPHI ',88,Metmin(11),Metmax(11),0.) call hbook1(876 ,'HTCF ',88,Metmin(12),Metmax(12),0.) call hbook1(877 ,'ETATRKS ',88,Metmin(13),Metmax(13),0.) call hbook1(878 ,'PHITRKS ',88,Metmin(14),Metmax(14),0.) call hbook1(879 ,'PTTRKS ',88,Metmin(15),Metmax(15),0.) call hbook1(880 ,'SEEDPT ',88,Metmin(16),Metmax(16),0.) call hbook1(881 ,'SUMET ',88,Metmin(17),Metmax(17),0.) call hbook1(882 ,'SUMETJET ',88,Metmin(18),Metmax(18),0.) call hbook1(883 ,'METSIG ',88,Metmin(19),Metmax(19),0.) C INITIALIZE THE ELECTRON PLOTS C call hbook1(890 ,'NEL ',20,ELEmin(1,1),ELEmax(1,1),0.) do i=1,10 call hbook1(890+36*(i-1)+1,'EP4 X',88,ELEMIN(i,2),ELEMAX(i,2),0.) call hbook1(890+36*(i-1)+2,'EETCOR',88,ELEMIN(i,3),ELEMAX(i,3),0.) call hbook1(890+36*(i-1)+3,'EET',88,ELEMIN(i,4),ELEMAX(i,4),0.) call hbook1(890+36*(i-1)+4,'ECHARGE',88,ELEMIN(i,5),ELEMAX(i,5),0.) call hbook1(890+36*(i-1)+5,'EDTETA',88,ELEMIN(i,6),ELEMAX(i,6),0.) call hbook1(890+36*(i-1)+6,'EEVETA',88,ELEMIN(i,7),ELEMAX(i,7),0.) call hbook1(890+36*(i-1)+7,'EPHI',88,ELEMIN(i,8),ELEMAX(i,8),0.) call hbook1(890+36*(i-1)+8,'EPT',88,ELEMIN(i,9),ELEMAX(i,9),0.) call hbook1(890+36*(i-1)+9,'ESTAT',88,ELEMIN(i,10),ELEMAX(i,10),0.) call hbook1(890+36*(i-1)+10,'EZV',88,ELEMIN(i,11),ELEMAX(i,11),0.) call hbook1(890+36*(i-1)+11,'EEP',88,ELEMIN(i,12),ELEMAX(i,12),0.) call hbook1(890+36*(i-1)+12,'EHADEM',88,ELEMIN(i,13),ELEMAX(i,13),0.) call hbook1(890+36*(i-1)+13,'EXCES',88,ELEMIN(i,14),ELEMAX(i,14),0.) call hbook1(890+36*(i-1)+14,'EDELX',88,ELEMIN(i,15),ELEMAX(i,15),0.) call hbook1(890+36*(i-1)+15,'EZCES',88,ELEMIN(i,16),ELEMAX(i,16),0.) call hbook1(890+36*(i-1)+16,'EDELZ',88,ELEMIN(i,17),ELEMAX(i,17),0.) call hbook1(890+36*(i-1)+17,'ECHIS',88,ELEMIN(i,18),ELEMAX(i,18),0.) call hbook1(890+36*(i-1)+18,'ECHIW',88,ELEMIN(i,19),ELEMAX(i,19),0.) call hbook1(890+36*(i-1)+19,'EISO',88,ELEMIN(i,20),ELEMAX(i,20),0.) call hbook1(890+36*(i-1)+20,'ETISO',88,ELEMIN(i,21),ELEMAX(i,21),0.) call hbook1(890+36*(i-1)+21,'EIDWRD',88,ELEMIN(i,22),ELEMAX(i,22),0.) call hbook1(890+36*(i-1)+22,'EISWRD',88,ELEMIN(i,23),ELEMAX(i,23),0.) call hbook1(890+36*(i-1)+23,'ECONWRD',88,ELEMIN(i,24),ELEMAX(i,24),0.) call hbook1(890+36*(i-1)+24,'EZTRK',88,ELEMIN(i,25),ELEMAX(i,25),0.) call hbook1(890+36*(i-1)+25,'EFID',88,ELEMIN(i,26),ELEMAX(i,26),0.) call hbook1(890+36*(i-1)+26,'EDET',88,ELEMIN(i,27),ELEMAX(i,27),0.) call hbook1(890+36*(i-1)+27,'EBTE',88,ELEMIN(i,28),ELEMAX(i,28),0.) call hbook1(890+36*(i-1)+28,'ENASL',88,ELEMIN(i,29),ELEMAX(i,29),0.) call hbook1(890+36*(i-1)+29,'ENSSL',88,ELEMIN(i,30),ELEMAX(i,30),0.) call hbook1(890+36*(i-1)+30,'EPTRK',88,ELEMIN(i,31),ELEMAX(i,31),0.) call hbook1(890+36*(i-1)+31,'ELSHR2',88,ELEMIN(i,32),ELEMAX(i,32),0.) call hbook1(890+36*(i-1)+32,'ECHI3',88,ELEMIN(i,33),ELEMAX(i,33),0.) call hbook1(890+36*(i-1)+33,'ECHID',88,ELEMIN(i,34),ELEMAX(i,34),0.) call hbook1(890+36*(i-1)+34,'EVTXOCC',88,ELEMIN(i,35),ELEMAX(i,35),0.) call hbook1(890+36*(i-1)+35,'ETIME',88,ELEMIN(i,36),ELEMAX(i,36),0.) call hbook1(890+36*(i-1)+36,'ETRIND',88,ELEMIN(i,37),ELEMAX(i,37),0.) ENDDO C INITIALIZE THE MUON PLOTS C call hbook1(1255 ,'NMU ',20,MUmin(1,1),MUmax(1,1),0.) do i=1,15 call hbook1(1255+30*(i-1)+1,'MURAWPT',88,MUMIN(i,2),MUMAX(i,2),0.) call hbook1(1255+30*(i-1)+2,'MUP4 X',88,MUMIN(i,3),MUMAX(i,3),0.) call hbook1(1255+30*(i-1)+3,'MUPT',88,MUMIN(i,4),MUMAX(i,4),0.) call hbook1(1255+30*(i-1)+4,'MUQ',88,MUMIN(i,5),MUMAX(i,5),0.) call hbook1(1255+30*(i-1)+5,'MUHAD',88,MUMIN(i,6),MUMAX(i,6),0.) call hbook1(1255+30*(i-1)+6,'MUEM',88,MUMIN(i,7),MUMAX(i,7),0.) call hbook1(1255+30*(i-1)+7,'MUEVETA',88,MUMIN(i,8),MUMAX(i,8),0.) call hbook1(1255+30*(i-1)+8,'MUPHI',88,MUMIN(i,9),MUMAX(i,9),0.) call hbook1(1255+30*(i-1)+9,'MUSTAT',88,MUMIN(i,10),MUMAX(i,10),0.) call hbook1(1255+30*(i-1)+10,'MUZV',88,MUMIN(i,11),MUMAX(i,11),0.) call hbook1(1255+30*(i-1)+11,'MUD0',88,MUMIN(i,12),MUMAX(i,12),0.) call hbook1(1255+30*(i-1)+12,'MUDELZ',88,MUMIN(i,13),MUMAX(i,13),0.) call hbook1(1255+30*(i-1)+13,'MUDELX',88,MUMIN(i,14),MUMAX(i,14),0.) call hbook1(1255+30*(i-1)+14,'MUISO',88,MUMIN(i,15),MUMAX(i,15),0.) call hbook1(1255+30*(i-1)+15,'MUTISO',88,MUMIN(i,16),MUMAX(i,16),0.) call hbook1(1255+30*(i-1)+16,'MUCRAY',88,MUMIN(i,17),MUMAX(i,17),0.) call hbook1(1255+30*(i-1)+17,'MUIDWRD',88,MUMIN(i,18),MUMAX(i,18),0.) call hbook1(1255+30*(i-1)+18,'MUISWRD',88,MUMIN(i,19),MUMAX(i,19),0.) call hbook1(1255+30*(i-1)+19,'MUZTRK',88,MUMIN(i,20),MUMAX(i,20),0.) call hbook1(1255+30*(i-1)+20,'MUFID',88,MUMIN(i,21),MUMAX(i,21),0.) call hbook1(1255+30*(i-1)+21,'MUDET',88,MUMIN(i,22),MUMAX(i,22),0.) call hbook1(1255+30*(i-1)+22,'MUBTE',88,MUMIN(i,23),MUMAX(i,23),0.) call hbook1(1255+30*(i-1)+23,'MUNASL',88,MUMIN(i,24),MUMAX(i,24),0.) call hbook1(1255+30*(i-1)+24,'MUNSSL',88,MUMIN(i,25),MUMAX(i,25),0.) call hbook1(1255+30*(i-1)+25,'MUCHZ0',88,MUMIN(i,26),MUMAX(i,26),0.) call hbook1(1255+30*(i-1)+26,'MUCHD0',88,MUMIN(i,27),MUMAX(i,27),0.) call hbook1(1255+30*(i-1)+27,'MUTIME',88,MUMIN(i,28),MUMAX(i,28),0.) call hbook1(1255+30*(i-1)+28,'MUVTXOCC',88,MUMIN(i,29),MUMAX(i,29),0.) call hbook1(1255+30*(i-1)+29,'MUDTETA',88,MUMIN(i,30),MUMAX(i,30),0.) call hbook1(1255+30*(i-1)+30,'MUTRIND',88,MUMIN(i,31),MUMAX(i,31),0.) ENDDO C INITIALIZE THE TAU PLOTS C call hbook1(1795 ,'NTau ',20,Taumin(1,1),Taumax(1,1),0.) do i=1,10 call hbook1(1795+18*(i-1)+1,'TP4 X',88,TAUMIN(i,2),TAUMAX(i,2),0.) call hbook1(1795+18*(i-1)+2,'TET',88,TAUMIN(i,3),TAUMAX(i,3),0.) call hbook1(1795+18*(i-1)+3,'TQ',88,TAUMIN(i,4),TAUMAX(i,4),0.) call hbook1(1795+18*(i-1)+4,'TZV',88,TAUMIN(i,5),TAUMAX(i,5),0.) call hbook1(1795+18*(i-1)+5,'TDTETA',88,TAUMIN(i,6),TAUMAX(i,6),0.) call hbook1(1795+18*(i-1)+6,'TEVETA',88,TAUMIN(i,7),TAUMAX(i,7),0.) call hbook1(1795+18*(i-1)+7,'TPHI',88,TAUMIN(i,8),TAUMAX(i,8),0.) call hbook1(1795+18*(i-1)+8,'TPT',88,TAUMIN(i,9),TAUMAX(i,9),0.) call hbook1(1795+18*(i-1)+9,'TPTVIS',88,TAUMIN(i,10),TAUMAX(i,10),0.) call hbook1(1795+18*(i-1)+10,'TSTAT',88,TAUMIN(i,11),TAUMAX(i,11),0.) call hbook1(1795+18*(i-1)+11,'TDELR',88,TAUMIN(i,12),TAUMAX(i,12),0.) call hbook1(1795+18*(i-1)+12,'TEMF',88,TAUMIN(i,13),TAUMAX(i,13),0.) call hbook1(1795+18*(i-1)+13,'TMASS',88,TAUMIN(i,14),TAUMAX(i,14),0.) call hbook1(1795+18*(i-1)+14,'TPTR 1',88,TAUMIN(i,15),TAUMAX(i,15),0.) call hbook1(1795+18*(i-1)+15,'TNTR10',88,TAUMIN(i,16),TAUMAX(i,16),0.) call hbook1(1795+18*(i-1)+16,'TNPI0',88,TAUMIN(i,17),TAUMAX(i,17),0.) call hbook1(1795+18*(i-1)+17,'TMUSTB',88,TAUMIN(i,18),TAUMAX(i,18),0.) call hbook1(1795+18*(i-1)+18,'TTRIND 1',88,TAUMIN(i,19),TAUMAX(i,19),0.) ENDDO C INITIALIZE THE DILEPTON PLOTS C call hbook1(1980 ,'NDIL',20,DILMIN(1),DILMAX(1),0.) c call hbook1(1981 ,'TZFLAG',20,DILMIN(2),DILMAX(2),0.) call hbook1(1982 ,'DILWRD',20,DILMIN(3),DILMAX(3),0.) call hbook1(1983 ,'METDILEP',20,DILMIN(4),DILMAX(4),0.) call hbook1(1984 ,'NJTDILEP',20,DILMIN(5),DILMAX(5),0.) call hbook1(1985 ,'DPHIMINJ',20,DILMIN(6),DILMAX(6),0.) call hbook1(1986 ,'DPHIMINL',20,DILMIN(7),DILMAX(7),0.) call hbook1(1987 ,'DILMODE',20,DILMIN(8),DILMAX(8),0.) call hbook1(1988 ,'DILTYPE1',20,DILMIN(9),DILMAX(9),0.) call hbook1(1989 ,'EEMASS',20,DILMIN(10),DILMAX(10),0.) call hbook1(1990 ,'MUMUMASS',20,DILMIN(11),DILMAX(11),0.) call hbook1(1991 ,'TAUMASS',20,DILMIN(12),DILMAX(12),0.) call hbook1(1992 ,'PHOMASS',20,DILMIN(13),DILMAX(13),0.) C INITIALIZE THE PHOTON PLOTS C call hbook1(1995 ,'NPHO ',20,PHOmin(1,1),PHOmax(1,1),0.) call hbook1(1996,'NCENCL',88,PHOMIN(1,2),PHOMAX(1,2),0.) call hbook1(1997,'NPLUGCL',88,PHOMIN(1,3),PHOMAX(1,3),0.) call hbook1(1998,'NPHOIDL',88,PHOMIN(1,4),PHOMAX(1,4),0.) call hbook1(1999,'NPHOISOL',88,PHOMIN(1,5),PHOMAX(1,5),0.) call hbook1(2000,'NCMUS',88,PHOMIN(1,6),PHOMAX(1,6),0.) call hbook1(2001,'COSTUB',88,PHOMIN(1,7),PHOMAX(1,7),0.) do i=1,10 call hbook1(2001+53*(i-1)+1,'PDETECT',88,PHOMIN(i,10),PHOMAX(i,10),0.) call hbook1(2001+53*(i-1)+2,'PZV',88,PHOMIN(i,11),PHOMAX(i,11),0.) call hbook1(2001+53*(i-1)+3,'PET',88,PHOMIN(i,12),PHOMAX(i,12),0.) call hbook1(2001+53*(i-1)+4,'PETC',88,PHOMIN(i,13),PHOMAX(i,13),0.) call hbook1(2001+53*(i-1)+5,'PP4 X',88,PHOMIN(i,14),PHOMAX(i,14),0.) call hbook1(2001+53*(i-1)+6,'PCO4',88,PHOMIN(i,15),PHOMAX(i,15),0.) call hbook1(2001+53*(i-1)+7,'PCESX',88,PHOMIN(i,16),PHOMAX(i,16),0.) call hbook1(2001+53*(i-1)+8,'PCESZ',88,PHOMIN(i,17),PHOMAX(i,17),0.) call hbook1(2001+53*(i-1)+9,'PCPR5PH',88,PHOMIN(i,18),PHOMAX(i,18),0.) call hbook1(2001+53*(i-1)+10,'PCPR5PS',88,PHOMIN(i,19),PHOMAX(i,19),0.) call hbook1(2001+53*(i-1)+11,'PDTETA',88,PHOMIN(i,20),PHOMAX(i,20),0.) call hbook1(2001+53*(i-1)+12,'PEVETA',88,PHOMIN(i,21),PHOMAX(i,21),0.) call hbook1(2001+53*(i-1)+13,'PCESE',88,PHOMIN(i,22),PHOMAX(i,22),0.) call hbook1(2001+53*(i-1)+14,'PPHI',88,PHOMIN(i,23),PHOMAX(i,23),0.) call hbook1(2001+53*(i-1)+15,'PCO7',88,PHOMIN(i,24),PHOMAX(i,24),0.) call hbook1(2001+53*(i-1)+16,'PHADEM',88,PHOMIN(i,25),PHOMAX(i,25),0.) call hbook1(2001+53*(i-1)+17,'PLSHR',88,PHOMIN(i,26),PHOMAX(i,26),0.) call hbook1(2001+53*(i-1)+18,'PSTR2',88,PHOMIN(i,27),PHOMAX(i,27),0.) call hbook1(2001+53*(i-1)+19,'PSTRE2',88,PHOMIN(i,28),PHOMAX(i,28),0.) call hbook1(2001+53*(i-1)+20,'PWIR2',88,PHOMIN(i,29),PHOMAX(i,29),0.) call hbook1(2001+53*(i-1)+21,'PWIRE2',88,PHOMIN(i,30),PHOMAX(i,30),0.) call hbook1(2001+53*(i-1)+22,'PSUMPT4',88,PHOMIN(i,31),PHOMAX(i,31),0.) call hbook1(2001+53*(i-1)+23,'PVTXOCC',88,PHOMIN(i,32),PHOMAX(i,32),0.) call hbook1(2001+53*(i-1)+24,'PPT',88,PHOMIN(i,33),PHOMAX(i,33),0.) call hbook1(2001+53*(i-1)+25,'PN3D',88,PHOMIN(i,34),PHOMAX(i,34),0.) call hbook1(2001+53*(i-1)+26,'PCHI',88,PHOMIN(i,35),PHOMAX(i,35),0.) call hbook1(2001+53*(i-1)+27,'PCHI3X3',88,PHOMIN(i,36),PHOMAX(i,36),0.) call hbook1(2001+53*(i-1)+28,'PCHIETA',88,PHOMIN(i,37),PHOMAX(i,37),0.) call hbook1(2001+53*(i-1)+29,'PCHIPHI',88,PHOMIN(i,38),PHOMAX(i,38),0.) call hbook1(2001+53*(i-1)+30,'PSTAT',88,PHOMIN(i,39),PHOMAX(i,39),0.) call hbook1(2001+53*(i-1)+31,'PBOXISO',88,PHOMIN(i,40),PHOMAX(i,40),0.) c call hbook1(2001+53*(i-1)+32,'PCESFLG',88,PHOMIN(i,41),PHOMAX(i,41),0.) call hbook1(2001+53*(i-1)+33,'PE',88,PHOMIN(i,42),PHOMAX(i,42),0.) call hbook1(2001+53*(i-1)+34,'PCESWHT',88,PHOMIN(i,43),PHOMAX(i,43),0.) call hbook1(2001+53*(i-1)+35,'PCPRWHT',88,PHOMIN(i,44),PHOMAX(i,44),0.) call hbook1(2001+53*(i-1)+36,'PCESSLID',88,PHOMIN(i,45),PHOMAX(i,45),0.) call hbook1(2001+53*(i-1)+37,'PWRD',88,PHOMIN(i,46),PHOMAX(i,46),0.) call hbook1(2001+53*(i-1)+38,'PSTH',88,PHOMIN(i,47),PHOMAX(i,47),0.) call hbook1(2001+53*(i-1)+39,'PELIND',88,PHOMIN(i,48),PHOMAX(i,48),0.) call hbook1(2001+53*(i-1)+40,'PCHISTR',88,PHOMIN(i,49),PHOMAX(i,49),0.) call hbook1(2001+53*(i-1)+41,'PCHIWIR',88,PHOMIN(i,50),PHOMAX(i,50),0.) call hbook1(2001+53*(i-1)+42,'PCESSIG',88,PHOMIN(i,51),PHOMAX(i,51),0.) call hbook1(2001+53*(i-1)+43,'PCESPG',88,PHOMIN(i,52),PHOMAX(i,52),0.) call hbook1(2001+53*(i-1)+44,'PCESPB',88,PHOMIN(i,53),PHOMAX(i,53),0.) call hbook1(2001+53*(i-1)+45,'PCPRPG',88,PHOMIN(i,54),PHOMAX(i,54),0.) call hbook1(2001+53*(i-1)+46,'PCPRPB',88,PHOMIN(i,55),PHOMAX(i,55),0.) call hbook1(2001+53*(i-1)+47,'PPT2',88,PHOMIN(i,56),PHOMAX(i,56),0.) call hbook1(2001+53*(i-1)+48,'PNTRWD3',88,PHOMIN(i,57),PHOMAX(i,57),0.) call hbook1(2001+53*(i-1)+49,'PCESCPRX',88,PHOMIN(i,58),PHOMAX(i,58),0.) call hbook1(2001+53*(i-1)+50,'PTKCPRX',88,PHOMIN(i,59),PHOMAX(i,59),0.) call hbook1(2001+53*(i-1)+51,'PTIME',88,PHOMIN(i,60),PHOMAX(i,60),0.) call hbook1(2001+53*(i-1)+52,'PVCOR',88,PHOMIN(i,61),PHOMAX(i,61),0.) call hbook1(2001+53*(i-1)+53,'PLCOR',88,PHOMIN(i,62),PHOMAX(i,62),0.) ENDDO C INITIALIZE THE TAGGING PLOTS C call hbook1(2535,'Ntag',20,Tagmin(1,1),Tagmax(1,1),0.) call hbook1(2536,'Btgwrd',20,Tagmin(1,2),Tagmax(1,2),0.) do i=1,4 call hbook1(2536+4*(i-1)+1,'Lxy',88,TAGMIN(i,3),TAGMAX(i,3),0.) call hbook1(2536+4*(i-1)+2,'Lxyerr',88,TAGMIN(i,4),TAGMAX(i,4),0.) call hbook1(2536+4*(i-1)+3,'Ctau',88,TAGMIN(i,5),TAGMAX(i,5),0.) call hbook1(2536+4*(i-1)+4,'Tagmin',88,TAGMIN(i,6),TAGMAX(i,6),0.) enddo C INITIALIZE THE Vertex PLOTS C call hbook1(2555,'Npvert',20,Vermin(1,1),Vermax(1,1),0.) call hbook1(2556,'Vxprim x',20,Vermin(1,2),Vermax(1,2),0.) do i=1,20 call hbook1(2556+3*(i-1)+1,'Vclass',88,VERMIN(i,3),VERMAX(i,3),0.) call hbook1(2556+3*(i-1)+2,'Vz',88,VERMIN(i,4),VERMAX(i,4),0.) call hbook1(2556+3*(i-1)+3,'Vq',88,VERMIN(i,5),VERMAX(i,5),0.) enddo call hbook1(2617,'ptnov',88,Vermin(1,6),vermax(1,6),0.) call hbook1(2618,'vzev',88,Vermin(1,7),vermax(1,7),0.) call hbook1(2619,'vclassev',88,Vermin(1,8),vermax(1,8),0.) C INITIALIZE THE TRACK PLOTS C call hbook1(2625,'Ntr',20,Trkmin(1,1),Trkmax(1,1),0.) do i=1,20 !only doing the first 20 for now call hbook1(2625+23*(i-1)+1,'Trvind',88,VERMIN(i,2),VERMAX(i,2),0.) call hbook1(2625+23*(i-1)+2,'Trnsvx',88,VERMIN(i,3),VERMAX(i,3),0.) call hbook1(2625+23*(i-1)+3,'Trnasl',88,VERMIN(i,4),VERMAX(i,4),0.) call hbook1(2625+23*(i-1)+4,'Trnssl',88,VERMIN(i,5),VERMAX(i,5),0.) call hbook1(2625+23*(i-1)+5,'Trcot',88,VERMIN(i,6),VERMAX(i,6),0.) call hbook1(2625+23*(i-1)+6,'Trcurv',88,VERMIN(i,7),VERMAX(i,7),0.) call hbook1(2625+23*(i-1)+7,'Trz0',88,VERMIN(i,8),VERMAX(i,8),0.) call hbook1(2625+23*(i-1)+8,'Trd0',88,VERMIN(i,9),VERMAX(i,9),0.) call hbook1(2625+23*(i-1)+9,'Trq',88,VERMIN(i,10),VERMAX(i,10),0.) call hbook1(2625+23*(i-1)+10,'Trphi0',88,VERMIN(i,11),VERMAX(i,11),0.) call hbook1(2625+23*(i-1)+11,'Trp4 x',88,VERMIN(i,12),VERMAX(i,12),0.) call hbook1(2625+23*(i-1)+12,'Trpt',88,VERMIN(i,13),VERMAX(i,13),0.) call hbook1(2625+23*(i-1)+13,'Treta',88,VERMIN(i,14),VERMAX(i,14),0.) call hbook1(2625+23*(i-1)+14,'Trexeta',88,VERMIN(i,15),VERMAX(i,15),0.) call hbook1(2625+23*(i-1)+15,'Trexphi',88,VERMIN(i,16),VERMAX(i,16),0.) call hbook1(2625+23*(i-1)+16,'Trbcvind',88,VERMIN(i,17),VERMAX(i,17),0.) call hbook1(2625+23*(i-1)+17,'Trbcz0',88,VERMIN(i,18),VERMAX(i,18),0.) call hbook1(2625+23*(i-1)+18,'Trbcd0',88,VERMIN(i,19),VERMAX(i,19),0.) call hbook1(2625+23*(i-1)+19,'Trbcq',88,VERMIN(i,20),VERMAX(i,20),0.) call hbook1(2625+23*(i-1)+20,'Trbcp4 1',88,VERMIN(i,21),VERMAX(i,21),0.) call hbook1(2625+23*(i-1)+21,'Trbcpt',88,VERMIN(i,22),VERMAX(i,22),0.) call hbook1(2625+23*(i-1)+22,'Trbceta',88,VERMIN(i,23),VERMAX(i,23),0.) call hbook1(2625+23*(i-1)+23,'Triso4',88,VERMIN(i,24),VERMAX(i,24),0.) enddo C chopt = ' ' C call hbprof(400,'Photon Et Raw', C & 100,55242.,71023.,0.0,330.0,chopt) C call hbprof(401,'Photon CES x', C & 100,55242.,71023.,-22.0,22.0,chopt) return end C a routine for writing histograms to an rz file C*************************************************************************** subroutine hist_write implicit none include 'cwn5.inc' integer icycle call hcdir('//hist/stntuple',' ') call hrout(0,icycle,'N') call hrend('hist') return end C an event selection filter before filling histograms C*************************************************************************** subroutine hist_fill implicit none include 'cwn5.inc' c include 'ghib$data:runsum.inc' integer i,j,k real frun,fnpvert,fpn3d,pwht integer nphotight,njethad,phoindex(10),jetindex(20) real coriso,phet real phopb call hcdir('//hist/stntuple',' ') C Fill General Variable Plots C call hfill(10 , REAL(Run) , 0., 1.) call hfill(11 , REAL(Event) , 0., 1.) call hfill(12 , ILUM , 0., 1.) call hfill(13 , REAL(Trigbits) , 0., 1.) c call hfill(14 , Goodrun , 0., 1.) call hfill(15 , real(Brcode) , 0., 1.) c call hfill(16 , Goodtrig , 0., 1.) call hfill(17 , real(Trigword) , 0., 1.) call hfill(18 , real(Njet10) , 0., 1.) call hfill(19 , real(Njet15) , 0., 1.) call hfill(20 , real(Njet20) , 0., 1.) call hfill(21 , real(Njet25) , 0., 1.) call hfill(22 , real(Nhpte) , 0., 1.) call hfill(23 , real(Nlcent) , 0., 1.) call hfill(24 , real(Nlplug) , 0., 1.) call hfill(25 , real(Nhpte) , 0., 1.) call hfill(26 , real(Ncmio) , 0., 1.) call hfill(27 , real(Ngtau) , 0., 1.) call hfill(28 , real(Nlepton) , 0., 1.) call hfill(29 , real(Etout) , 0., 1.) C Fill JET PLOTS C call hfill(35 ,real(Njt) ,0., 1.) do i=1,Njt call hfill(35+30*(i-1)+1, JTET(i) ,0., 1.) call hfill(35+30*(i-1)+2, JTM(i) ,0., 1.) call hfill(35+30*(i-1)+3, JTP4(1,i) ,0., 1.) call hfill(35+30*(i-1)+4, JTEVETA(i) ,0., 1.) call hfill(35+30*(i-1)+5, JTDTETA(i) ,0., 1.) call hfill(35+30*(i-1)+6, JTPHI(i) ,0., 1.) call hfill(35+30*(i-1)+7, JTEMF(i) ,0., 1.) call hfill(35+30*(i-1)+8, real(JTWORD(i)) ,0., 1.) call hfill(35+30*(i-1)+9, JTCORFM(i) ,0., 1.) call hfill(35+30*(i-1)+10, JTCORFD(i) ,0., 1.) call hfill(35+30*(i-1)+11, JTCORFA(i) ,0., 1.) call hfill(35+30*(i-1)+12, real(JTNTR(i)) ,0., 1.) call hfill(35+30*(i-1)+13, real(JTNTRP(i)) ,0., 1.) call hfill(35+30*(i-1)+14, JTEP(i) ,0., 1.) call hfill(35+30*(i-1)+15, real(JTNJTRK0(i)),0., 1.) call hfill(35+30*(i-1)+16, real(JTTAG(i)) ,0., 1.) call hfill(35+30*(i-1)+17, JTEGRD(i) ,0., 1.) call hfill(35+30*(i-1)+18, JVBPB(i) ,0., 1.) call hfill(35+30*(i-1)+19, JPBPB(i) ,0., 1.) call hfill(35+30*(i-1)+20, SCBPB(i) ,0., 1.) call hfill(35+30*(i-1)+21, SCNPB(i) ,0., 1.) call hfill(35+30*(i-1)+22, JPPPB(i) ,0., 1.) call hfill(35+30*(i-1)+23, JPNPB(i) ,0., 1.) call hfill(35+30*(i-1)+24, JTAU1(i) ,0., 1.) call hfill(35+30*(i-1)+25, JTAU2(i) ,0., 1.) call hfill(35+30*(i-1)+26, JTZVTX(i) ,0., 1.) call hfill(35+30*(i-1)+27, JTDZVTX(i) ,0., 1.) call hfill(35+30*(i-1)+28, JT_PTOUT(i),0., 1.) call hfill(35+30*(i-1)+29, JTVPT(1,i) ,0., 1.) call hfill(35+30*(i-1)+30, JTNOVPT(i) ,0., 1.) ENDDO C Fill JT1S Plots C call hfill(640, real(Nj1) , 0., 1.) do i=1,Nj1 call hfill(640+11*(i-1)+1, J1ET(i) ,0. ,1.) call hfill(640+11*(i-1)+2, J1M(i) ,0. ,1.) call hfill(640+11*(i-1)+3, J1P4(1,i) ,0. ,1.) call hfill(640+11*(i-1)+4, J1EVETA(i) ,0. ,1.) call hfill(640+11*(i-1)+5, J1DTETA(i) ,0. ,1.) call hfill(640+11*(i-1)+6, J1PHI(i) ,0. ,1.) call hfill(640+11*(i-1)+7, real(J1WORD(i)) ,0. ,1.) call hfill(640+11*(i-1)+8, real(J1JTS(i)) ,0. ,1.) call hfill(640+11*(i-1)+9, J1CORFM(i) ,0. ,1.) call hfill(640+11*(i-1)+10, J1CORFD(i),0. ,1.) call hfill(640+11*(i-1)+11, J1EGRD(i) ,0. ,1.) Enddo C Fill THE MET PLOTS C call hfill(865 , MET(5) , 0., 1.) call hfill(866 , METPHI(5) , 0., 1.) call hfill(867 , METJT , 0., 1.) call hfill(868 , METJTPHI , 0., 1.) call hfill(869 , METJ1 , 0., 1.) call hfill(870 , METJ1PHI , 0., 1.) call hfill(871 , METC , 0., 1.) call hfill(872 , METCPHI , 0., 1.) call hfill(873 , HTC , 0., 1.) call hfill(874 , METCF , 0., 1.) call hfill(875 , METCFPHI , 0., 1.) call hfill(876 , HTCF , 0., 1.) call hfill(877 , ETATRKS , 0., 1.) call hfill(878 , PHITRKS , 0., 1.) call hfill(879 , PTTRKS , 0., 1.) call hfill(880 , SEEDPT , 0., 1.) call hfill(881 , SUMET , 0., 1.) call hfill(882 , SUMETJET , 0., 1.) call hfill(883 , METSIG , 0., 1.) C Fill THE ELECTRON PLOTS C call hfill(890 , real(NEL), 0., 1.) do i=1,NEL call hfill(890+36*(i-1)+1, EP4(1,i) , 0., 1.) call hfill(890+36*(i-1)+2, EETCOR(i) , 0., 1.) call hfill(890+36*(i-1)+3, EET(i) , 0., 1.) call hfill(890+36*(i-1)+4, ECHARGE(i) , 0., 1.) call hfill(890+36*(i-1)+5, EDTETA(i) , 0., 1.) call hfill(890+36*(i-1)+6, EEVETA(i) , 0., 1.) call hfill(890+36*(i-1)+7, EPHI(i) , 0., 1.) call hfill(890+36*(i-1)+8, EPT(i) , 0., 1.) call hfill(890+36*(i-1)+9, ESTAT(i) , 0., 1.) call hfill(890+36*(i-1)+10, EZV(i) , 0., 1.) call hfill(890+36*(i-1)+11, EEP(i) , 0., 1.) call hfill(890+36*(i-1)+12, EHADEM(i) , 0., 1.) call hfill(890+36*(i-1)+13, EXCES(i) , 0., 1.) call hfill(890+36*(i-1)+14, EDELX(i) , 0., 1.) call hfill(890+36*(i-1)+15, EZCES(i) , 0., 1.) call hfill(890+36*(i-1)+16, EDELZ(i) , 0., 1.) call hfill(890+36*(i-1)+17, ECHIS(i) , 0., 1.) call hfill(890+36*(i-1)+18, ECHIW(i) , 0., 1.) call hfill(890+36*(i-1)+19, EISO(i) , 0., 1.) call hfill(890+36*(i-1)+20, ETISO(i) , 0., 1.) call hfill(890+36*(i-1)+21, real(EIDWRD(i)) , 0., 1.) call hfill(890+36*(i-1)+22, real(EISWRD(i)) , 0., 1.) call hfill(890+36*(i-1)+23, real(ECONWRD(i)), 0., 1.) call hfill(890+36*(i-1)+24, EZTRK(i) , 0., 1.) call hfill(890+36*(i-1)+25, real(EFID(i)) , 0., 1.) call hfill(890+36*(i-1)+26, real(EDET(i)) , 0., 1.) call hfill(890+36*(i-1)+27, EBTE(i) , 0., 1.) call hfill(890+36*(i-1)+28, real(ENASL(i)) , 0., 1.) call hfill(890+36*(i-1)+29, real(ENSSL(i)) , 0., 1.) call hfill(890+36*(i-1)+30, EPTRK(i) , 0., 1.) call hfill(890+36*(i-1)+31, ELSHR2(i) , 0., 1.) call hfill(890+36*(i-1)+32, ECHI3(i) , 0., 1.) call hfill(890+36*(i-1)+33, ECHID(i) , 0., 1.) call hfill(890+36*(i-1)+34, EVTXOCC(i), 0., 1.) call hfill(890+36*(i-1)+35, ETIME(i) , 0., 1.) call hfill(890+36*(i-1)+36, real(ETRIND(i)) , 0., 1.) ENDDO C Fill THE MUON PLOTS C call hfill(1255 , real(NMU), 0., 1.) do i=1,nmu call hfill(1255+30*(i-1)+1, MURAWPT(i) , 0., 1.) call hfill(1255+30*(i-1)+2, MUP4(1,i) , 0., 1.) call hfill(1255+30*(i-1)+3, MUPT(i) , 0., 1.) call hfill(1255+30*(i-1)+4, MUQ(i) , 0., 1.) call hfill(1255+30*(i-1)+5, MUHAD(i) , 0., 1.) call hfill(1255+30*(i-1)+6, MUEM(i) , 0., 1.) call hfill(1255+30*(i-1)+7, MUEVETA(i) , 0., 1.) call hfill(1255+30*(i-1)+8, MUPHI(i) , 0., 1.) call hfill(1255+30*(i-1)+9, MUSTAT(i) , 0., 1.) call hfill(1255+30*(i-1)+10, MUZV(i) , 0., 1.) call hfill(1255+30*(i-1)+11, MUD0(i) , 0., 1.) call hfill(1255+30*(i-1)+12, MUDELZ(i) , 0., 1.) call hfill(1255+30*(i-1)+13, MUDELX(i) , 0., 1.) call hfill(1255+30*(i-1)+14, MUISO(i) , 0., 1.) call hfill(1255+30*(i-1)+15, MUTISO(i) , 0., 1.) call hfill(1255+30*(i-1)+16, real(MUCRAY(i)) , 0., 1.) call hfill(1255+30*(i-1)+17, real(MUIDWRD(i)) , 0., 1.) call hfill(1255+30*(i-1)+18, real(MUISWRD(i)) , 0., 1.) call hfill(1255+30*(i-1)+19, MUZTRK(i) , 0., 1.) call hfill(1255+30*(i-1)+20, real(MUFID(i)) , 0., 1.) call hfill(1255+30*(i-1)+21, real(MUDET(i)) , 0., 1.) call hfill(1255+30*(i-1)+22, MUBTE(i) , 0., 1.) call hfill(1255+30*(i-1)+23, real(MUNASL(i)) , 0., 1.) call hfill(1255+30*(i-1)+24, real(MUNSSL(i)) , 0., 1.) call hfill(1255+30*(i-1)+25, MUCHZ0(i) , 0., 1.) call hfill(1255+30*(i-1)+26, MUCHD0(i) , 0., 1.) call hfill(1255+30*(i-1)+27, MUTIME(i) , 0., 1.) call hfill(1255+30*(i-1)+28, MUVTXOCC(i), 0., 1.) call hfill(1255+30*(i-1)+29, MUDTETA(i) , 0., 1.) call hfill(1255+30*(i-1)+30, real(MUTRIND(i)) , 0., 1.) ENDDO C INITIALIZE THE TAU PLOTS C call hfill(1795 , real(NTAU), 0., 1.) do i=1,Ntau call hfill(1795+18*(i-1)+1, TP4(1,i) , 0., 1.) call hfill(1795+18*(i-1)+2, TET(i) , 0., 1.) call hfill(1795+18*(i-1)+3, real(TQ(i)) , 0., 1.) call hfill(1795+18*(i-1)+4, TZV(i) , 0., 1.) call hfill(1795+18*(i-1)+5, TDTETA(i) , 0., 1.) call hfill(1795+18*(i-1)+6, TEVETA(i) , 0., 1.) call hfill(1795+18*(i-1)+7, TPHI(i) , 0., 1.) call hfill(1795+18*(i-1)+8, TPT(i) , 0., 1.) call hfill(1795+18*(i-1)+9, TPTVIS(i) , 0., 1.) call hfill(1795+18*(i-1)+10, TSTAT(i) , 0., 1.) call hfill(1795+18*(i-1)+11, TDELR(i) , 0., 1.) call hfill(1795+18*(i-1)+12, TEMF(i) , 0., 1.) call hfill(1795+18*(i-1)+13, TMASS(i) , 0., 1.) call hfill(1795+18*(i-1)+14, TPTTR(1,i) , 0., 1.) call hfill(1795+18*(i-1)+15, real(TNTR10(i)) , 0., 1.) call hfill(1795+18*(i-1)+16, real(TNPI0(i)) , 0., 1.) call hfill(1795+18*(i-1)+17, real(TMUSTUB(i)) , 0., 1.) call hfill(1795+18*(i-1)+18, real(TTRIND(1,i)), 0., 1.) ENDDO C INITIALIZE THE DILEPTON PLOTS C call hfill(1980, real(NDIL), 0., 1.) c call hfill(1981, TZFLAG, 0., 1.) call hfill(1982, real(DILWRD), 0., 1.) call hfill(1983, METDILEP, 0., 1.) call hfill(1984, real(NJTDILEP), 0., 1.) call hfill(1985, DPHIMINJ, 0., 1.) call hfill(1986, DPHIMINL, 0., 1.) call hfill(1987, real(DILMODE), 0., 1.) call hfill(1988, real(DILTYPE1), 0., 1.) call hfill(1989, EEMASS, 0., 1.) call hfill(1990, MUMUMASS, 0., 1.) call hfill(1991, TAUMASS, 0., 1.) call hfill(1992, PHOMASS, 0., 1.) C INITIALIZE THE PHOTON PLOTS C call hfill(1995, real(NPHO) , 0., 1.) call hfill(1996, real(NCENCL) , 0., 1.) call hfill(1997, real(NPLUGCL) , 0., 1.) call hfill(1998, real(NPHOIDL) , 0., 1.) call hfill(1999, real(NPHOISOL) , 0., 1.) call hfill(2000, real(NCMUS) , 0., 1.) call hfill(2001, real(COSTUB) , 0., 1.) do i=1,NPHO call hfill(2001+53*(i-1)+1, real(PDETECT(i)) , 0., 1.) call hfill(2001+53*(i-1)+2, PZV(i) , 0., 1.) call hfill(2001+53*(i-1)+3, PET(i) , 0., 1.) call hfill(2001+53*(i-1)+4, PETC(i) , 0., 1.) call hfill(2001+53*(i-1)+5, PP4(1,i) , 0., 1.) call hfill(2001+53*(i-1)+6, PCO4(i) , 0., 1.) call hfill(2001+53*(i-1)+7, PCESX(i) , 0., 1.) call hfill(2001+53*(i-1)+8, PCESZ(i) , 0., 1.) call hfill(2001+53*(i-1)+9, PCPR5PH(i) , 0., 1.) call hfill(2001+53*(i-1)+10, PCPR5PS(i) , 0., 1.) call hfill(2001+53*(i-1)+11, PDTETA(i) , 0., 1.) call hfill(2001+53*(i-1)+12, PEVETA(i) , 0., 1.) call hfill(2001+53*(i-1)+13, PCESE(i) , 0., 1.) call hfill(2001+53*(i-1)+14, PPHI(i) , 0., 1.) call hfill(2001+53*(i-1)+15, PCO7(i) , 0., 1.) call hfill(2001+53*(i-1)+16, PHADEM(i) , 0., 1.) call hfill(2001+53*(i-1)+17, PLSHR(i) , 0., 1.) call hfill(2001+53*(i-1)+18, PSTR2(i) , 0., 1.) call hfill(2001+53*(i-1)+19, PSTRE2(i) , 0., 1.) call hfill(2001+53*(i-1)+20, PWIR2(i) , 0., 1.) call hfill(2001+53*(i-1)+21, PWIRE2(i) , 0., 1.) call hfill(2001+53*(i-1)+22, PSUMPT4(i) , 0., 1.) call hfill(2001+53*(i-1)+23, PVTXOCC(i) , 0., 1.) call hfill(2001+53*(i-1)+24, PPT(i) , 0., 1.) call hfill(2001+53*(i-1)+25, real(PN3D(i)) , 0., 1.) call hfill(2001+53*(i-1)+26, PCHI(i) , 0., 1.) call hfill(2001+53*(i-1)+27, PCHI3X3(i) , 0., 1.) call hfill(2001+53*(i-1)+28, PCHIETA(i) , 0., 1.) call hfill(2001+53*(i-1)+29, PCHIPHI(i) , 0., 1.) call hfill(2001+53*(i-1)+30, PSTAT(i) , 0., 1.) call hfill(2001+53*(i-1)+31, PBOXISO(i) , 0., 1.) c call hfill(2001+53*(i-1)+32, PCESFLG(i) , 0., 1.) call hfill(2001+53*(i-1)+33, PE(i) , 0., 1.) call hfill(2001+53*(i-1)+34, PCESWHT(i) , 0., 1.) call hfill(2001+53*(i-1)+35, PCPRWHT(i) , 0., 1.) call hfill(2001+53*(i-1)+36, PCESSLIDE(i), 0., 1.) call hfill(2001+53*(i-1)+37, real(PWRD(i)) , 0., 1.) call hfill(2001+53*(i-1)+38, PSTH(i) , 0., 1.) call hfill(2001+53*(i-1)+39, real(PELIND(i)) , 0., 1.) call hfill(2001+53*(i-1)+40, PCHISTR(i) , 0., 1.) call hfill(2001+53*(i-1)+41, PCHIWIR(i) , 0., 1.) call hfill(2001+53*(i-1)+42, PCESSIG(i) , 0., 1.) call hfill(2001+53*(i-1)+43, PCESPG(i) , 0., 1.) call hfill(2001+53*(i-1)+44, PCESPB(i) , 0., 1.) call hfill(2001+53*(i-1)+45, PCPRPG(i) , 0., 1.) call hfill(2001+53*(i-1)+46, PCPRPB(i) , 0., 1.) call hfill(2001+53*(i-1)+47, PPT2(i) , 0., 1.) call hfill(2001+53*(i-1)+48, real(PNTRWD3(i)) , 0., 1.) call hfill(2001+53*(i-1)+49, PCESCPRX(i), 0., 1.) call hfill(2001+53*(i-1)+50, PTKCPRX(i) , 0., 1.) call hfill(2001+53*(i-1)+51, PTIME(i) , 0., 1.) call hfill(2001+53*(i-1)+52, PVCOR(i) , 0., 1.) call hfill(2001+53*(i-1)+53, PLCOR(i) , 0., 1.) ENDDO C Fill THE TAGGING PLOTS C call hfill(2535, real(Ntag), 0., 1.) call hfill(2536, real(Btgwrd), 0., 1.) do i=1,Ntag call hfill(2536+4*(i-1)+1, Lxy(i), 0., 1.) call hfill(2536+4*(i-1)+2, Lxyerr(i), 0., 1.) call hfill(2536+4*(i-1)+3, Ctau(i), 0., 1.) call hfill(2536+4*(i-1)+4, TagCHI(i), 0., 1.) enddo C Fill THE Vertex PLOTS C call hfill(2555, real(Npvert), 0., 1.) call hfill(2556, Vxprim(1), 0., 1.) do i=1,Npvert call hfill(2556+3*(i-1)+1, real(Vclass(i)), 0., 1.) call hfill(2556+3*(i-1)+2, Vz(i), 0., 1.) call hfill(2556+3*(i-1)+3, Vq(i), 0., 1.) enddo call hfill(2617, ptnov, 0., 1.) call hfill(2618, vzev, 0., 1.) call hfill(2619, real(vclassev),0., 1.) C INITIALIZE THE TRACK PLOTS C call hfill(2625, real(Ntr), 0., 1.) do i=1,Ntr !only doing the first 20 for now call hfill(2625+23*(i-1)+1, real(Trvind(i)), 0., 1.) call hfill(2625+23*(i-1)+2, real(Trnsvx(i)), 0., 1.) call hfill(2625+23*(i-1)+3, real(Trnasl(i)), 0., 1.) call hfill(2625+23*(i-1)+4, real(Trnssl(i)), 0., 1.) call hfill(2625+23*(i-1)+5, Trcot(i), 0., 1.) call hfill(2625+23*(i-1)+6, Trcurv(i), 0., 1.) call hfill(2625+23*(i-1)+7, Trz0(i), 0., 1.) call hfill(2625+23*(i-1)+8, Trd0(i), 0., 1.) call hfill(2625+23*(i-1)+9, Trq(i), 0., 1.) call hfill(2625+23*(i-1)+10, Trphi0(i), 0., 1.) call hfill(2625+23*(i-1)+11, Trp4(1,i), 0., 1.) call hfill(2625+23*(i-1)+12, Trpt(i), 0., 1.) call hfill(2625+23*(i-1)+13, Treta(i), 0., 1.) call hfill(2625+23*(i-1)+14, Trexeta(i), 0., 1.) call hfill(2625+23*(i-1)+15, Trexphi(i), 0., 1.) call hfill(2625+23*(i-1)+16, real(Trbcvind(i)),0., 1.) call hfill(2625+23*(i-1)+17, Trbcz0(i), 0., 1.) call hfill(2625+23*(i-1)+18, Trbcd0(i), 0., 1.) call hfill(2625+23*(i-1)+19, Trbcq(i), 0., 1.) call hfill(2625+23*(i-1)+20, Trbcp4(1,i),0., 1.) call hfill(2625+23*(i-1)+21, Trbcpt(i), 0., 1.) call hfill(2625+23*(i-1)+22, Trbceta(i), 0., 1.) call hfill(2625+23*(i-1)+23, Triso4(i), 0., 1.) enddo return end Subroutine ENDVAL CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C Subroutine used to determine the MIN and Max Variable Values In the C Ntuple CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC Implicit None Include 'cwn5.inc' c Include 'ghib$data:runsum.inc' Integer i,j,k REAL GENMAX(21),GENMIN(21),JETSMAX(0:20,31),JETSMIN(0:20,31) REAL JT1SMAX(0:20,12),JT1SMIN(0:20,12),METMAX(19),METMIN(19) REAL ELEMAX(0:10,37),ELEMIN(0:10,37) REAL MUMAX(0:10,31),MUMIN(0:10,31) REAL TAUMAX(0:10,19),TAUMIN(0:10,19),DILMAX(13),DILMIN(13) REAL PHOMAX(0:10,62),PHOMIN(0:10,62) REAL TAGMAX(0:4,6),TAGMIN(0:4,6) REAL VERMAX(0:20,8),VERMIN(0:20,8) REAL TRKMAX(0:100,24),TRKMIN(0:100,24) COMMON /LIMITS1/ GENMAX,GENMIN,JETSMAX,JETSMIN,JT1SMAX,JT1SMIN COMMON /LIMITS2/ METMAX,METMIN,ELEMAX,ELEMIN,MUMAX,MUMIN COMMON /LIMITS3/ TAUMAX,TAUMIN,DILMAX,DILMIN,PHOMAX,PHOMIN COMMON /LIMITS4/ TAGMAX,TAGMIN,VERMAX,VERMIN COMMON /LIMITS5/ TRKMAX,TRKMIN CC C General Variable Limits CC If (REAL(Version).gt.Genmax(1)) Genmax(1)=REAL(Version) If (REAL(Version).lt.Genmin(1)) Genmin(1)=REAL(Version) If (REAL(Run).gt.Genmax(2)) Genmax(2)=REAL(Run) If (REAL(Run).lt.Genmin(2)) Genmin(2)=REAL(Run) If (REAL(Event).gt.Genmax(3)) Genmax(3)=REAL(Event) If (REAL(Event).lt.Genmin(3)) Genmin(3)=REAL(Event) If (Ilum.gt.Genmax(4)) Genmax(4)=Ilum If (Ilum.lt.Genmin(4)) Genmin(4)=Ilum If (REAL(Trigbits).gt.Genmax(5)) Genmax(5)=REAL(Trigbits) If (REAL(Trigbits).lt.Genmin(5)) Genmin(5)=REAL(Trigbits) c If (Goodrun.gt.Genmax(6)) Genmax(6)=Goodrun c If (Goodrun.lt.Genmin(6)) Genmin(6)=Goodrun If (REAL(Brcode).gt.Genmax(7)) Genmax(7)=REAL(Brcode) If (REAL(Brcode).lt.Genmin(7)) Genmin(7)=REAL(Brcode) c If (Goodtrig.gt.Genmax(8)) Genmax(8)=Goodtrig c If (Goodtrig.lt.Genmin(8)) Genmin(8)=Goodtrig If (REAL(Trigword).gt.Genmax(9)) Genmax(9)=REAL(Trigword) If (REAL(Trigword).lt.Genmin(9)) Genmin(9)=REAL(Trigword) If (REAL(Njet10).gt.Genmax(10)) Genmax(10)=REAL(Njet10) If (REAL(Njet10).lt.Genmin(10)) Genmin(10)=REAL(Njet10) If (REAL(Njet15).gt.Genmax(11)) Genmax(11)=REAL(Njet15) If (REAL(Njet15).lt.Genmin(11)) Genmin(11)=REAL(Njet15) If (REAL(Njet20).gt.Genmax(12)) Genmax(12)=REAL(Njet20) If (REAL(Njet20).lt.Genmin(12)) Genmin(12)=REAL(Njet20) If (REAL(Njet25).gt.Genmax(13)) Genmax(13)=REAL(Njet25) If (REAL(Njet25).lt.Genmin(13)) Genmin(13)=REAL(Njet25) If (REAL(Nhpte).gt.Genmax(14)) Genmax(14)=REAL(Nhpte) If (REAL(Nhpte).lt.Genmin(14)) Genmin(14)=REAL(Nhpte) If (REAL(Nlcent).gt.Genmax(15)) Genmax(15)=REAL(Nlcent) If (REAL(Nlcent).lt.Genmin(15)) Genmin(15)=REAL(Nlcent) If (REAL(Nlplug).gt.Genmax(16)) Genmax(15)=REAL(Nlplug) If (REAL(Nlplug).lt.Genmin(16)) Genmin(15)=REAL(Nlplug) If (REAL(Nhptm).gt.Genmax(17)) Genmax(15)=REAL(Nhptm) If (REAL(Nhptm).lt.Genmin(17)) Genmin(15)=REAL(Nhptm) If (REAL(Ncmio).gt.Genmax(18)) Genmax(15)=REAL(Ncmio) If (REAL(Ncmio).lt.Genmin(18)) Genmin(15)=REAL(Ncmio) If (REAL(Ngtau).gt.Genmax(19)) Genmax(15)=REAL(Ngtau) If (REAL(Ngtau).lt.Genmin(19)) Genmin(15)=REAL(Ngtau) If (REAL(Nlepton).gt.Genmax(20)) Genmax(15)=REAL(Nlepton) If (REAL(Nlepton).lt.Genmin(20)) Genmin(15)=REAL(Nlepton) If (Etout.gt.Genmax(21)) Genmax(15)=Etout If (Etout.lt.Genmin(21)) Genmin(15)=Etout CC C Jet Specific Variables CC If (REAL(Njt).gt.Jetsmax(1,1)) Jetsmax(1,1)=REAL(Njt) If (REAL(Njt).lt.Jetsmin(1,1)) Jetsmin(1,1)=REAL(Njt) Do i=1,njt If (Jtet(i).gt.Jetsmax(i,2)) Jetsmax(i,2)=Jtet(i) If (Jtet(i).lt.Jetsmin(i,2)) Jetsmin(i,2)=Jtet(i) If (Jtm(i).gt.Jetsmax(i,3)) Jetsmax(i,3)=Jtm(i) If (Jtm(i).lt.Jetsmin(i,3)) Jetsmin(i,3)=Jtm(i) If (Jtp4(1,i).gt.Jetsmax(i,4)) Jetsmax(i,4)=Jtp4(1,i) !x mom If (Jtp4(1,i).lt.Jetsmin(i,4)) Jetsmin(i,4)=Jtp4(1,i) !x mom If (Jteveta(i).gt.Jetsmax(i,5)) Jetsmax(i,5)=Jteveta(i) If (Jteveta(i).lt.Jetsmin(i,5)) Jetsmin(i,5)=Jteveta(i) If (Jtdteta(i).gt.Jetsmax(i,6)) Jetsmax(i,6)=Jtdteta(i) If (Jtdteta(i).lt.Jetsmin(i,6)) Jetsmin(i,6)=Jtdteta(i) If (Jtphi(i).gt.Jetsmax(i,7)) Jetsmax(i,7)=Jtphi(i) If (Jtphi(i).lt.Jetsmin(i,7)) Jetsmin(i,7)=Jtphi(i) If (Jtemf(i).gt.Jetsmax(i,8)) Jetsmax(i,8)=Jtemf(i) If (Jtemf(i).lt.Jetsmin(i,8)) Jetsmin(i,8)=Jtemf(i) If (REAL(Jtword(i)).gt.Jetsmax(i,9)) Jetsmax(i,9)=REAL(Jtword(i)) If (REAL(Jtword(i)).lt.Jetsmin(i,9)) Jetsmin(i,9)=REAL(Jtword(i)) If (Jtcorfm(i).gt.Jetsmax(i,10)) Jetsmax(i,10)=Jtcorfm(i) If (Jtcorfm(i).lt.Jetsmin(i,10)) Jetsmin(i,10)=Jtcorfm(i) If (Jtcorfd(i).gt.Jetsmax(i,11)) Jetsmax(i,11)=Jtcorfd(i) If (Jtcorfd(i).lt.Jetsmin(i,11)) Jetsmin(i,11)=Jtcorfd(i) If (Jtcorfa(i).gt.Jetsmax(i,12)) Jetsmax(i,12)=Jtcorfa(i) If (Jtcorfa(i).lt.Jetsmin(i,12)) Jetsmin(i,12)=Jtcorfa(i) If (REAL(Jtntr(i)).gt.Jetsmax(i,13)) Jetsmax(i,13)=REAL(Jtntr(i)) If (REAL(Jtntr(i)).lt.Jetsmin(i,13)) Jetsmin(i,13)=REAL(Jtntr(i)) If (REAL(Jtntrp(i)).gt.Jetsmax(i,14)) Jetsmax(i,14)=REAL(Jtntrp(i)) If (REAL(Jtntrp(i)).lt.Jetsmin(i,14)) Jetsmin(i,14)=REAL(Jtntrp(i)) If (Jtep(i).gt.Jetsmax(i,15)) Jetsmax(i,15)=Jtep(i) If (Jtep(i).lt.Jetsmin(i,15)) Jetsmin(i,15)=Jtep(i) If (REAL(Jtnjtrk0(i)).gt.Jetsmax(i,16)) Jetsmax(i,16)=REAL(Jtnjtrk0(i)) If (REAL(Jtnjtrk0(i)).lt.Jetsmin(i,16)) Jetsmin(i,16)=REAL(Jtnjtrk0(i)) If (REAL(Jttag(i)).gt.Jetsmax(i,17)) Jetsmax(i,17)=REAL(Jttag(i)) If (REAL(Jttag(i)).lt.Jetsmin(i,17)) Jetsmin(i,17)=REAL(Jttag(i)) If (Jtegrd(i).gt.Jetsmax(i,18)) Jetsmax(i,18)=Jtegrd(i) If (Jtegrd(i).lt.Jetsmin(i,18)) Jetsmin(i,18)=Jtegrd(i) If (Jvbpb(i).gt.Jetsmax(i,19)) Jetsmax(i,19)=Jvbpb(i) If (Jvbpb(i).lt.Jetsmin(i,19)) Jetsmin(i,19)=Jvbpb(i) If (Jpbpb(i).gt.Jetsmax(i,20)) Jetsmax(i,20)=Jpbpb(i) If (Jpbpb(i).lt.Jetsmin(i,20)) Jetsmin(i,20)=Jpbpb(i) If (Scbpb(i).gt.Jetsmax(i,21)) Jetsmax(i,21)=Scbpb(i) If (Scbpb(i).lt.Jetsmin(i,21)) Jetsmin(i,21)=Scbpb(i) If (Scnpb(i).gt.Jetsmax(i,22)) Jetsmax(i,22)=Scnpb(i) If (Scnpb(i).lt.Jetsmin(i,22)) Jetsmin(i,22)=Scnpb(i) If (Jpppb(i).gt.Jetsmax(i,23)) Jetsmax(i,23)=Jpppb(i) If (Jpppb(i).lt.Jetsmin(i,23)) Jetsmin(i,23)=Jpppb(i) If (Jpnpb(i).gt.Jetsmax(i,24)) Jetsmax(i,24)=Jpnpb(i) If (Jpnpb(i).lt.Jetsmin(i,24)) Jetsmin(i,24)=Jpnpb(i) If (Jtau1(i).gt.Jetsmax(i,25)) Jetsmax(i,25)=Jtau1(i) If (Jtau1(i).lt.Jetsmin(i,25)) Jetsmin(i,25)=Jtau1(i) If (Jtau2(i).gt.Jetsmax(i,26)) Jetsmax(i,26)=Jtau2(i) If (Jtau2(i).lt.Jetsmin(i,26)) Jetsmin(i,26)=Jtau2(i) If (Jtzvtx(i).gt.Jetsmax(i,27)) Jetsmax(i,27)=Jtzvtx(i) If (Jtzvtx(i).lt.Jetsmin(i,27)) Jetsmin(i,27)=Jtzvtx(i) If (Jtdzvtx(i).gt.Jetsmax(i,28)) Jetsmax(i,28)=Jtdzvtx(i) If (Jtdzvtx(i).lt.Jetsmin(i,28)) Jetsmin(i,28)=Jtdzvtx(i) If (Jt_ptout(i).gt.Jetsmax(i,29)) Jetsmax(i,29)=Jt_ptout(i) If (Jt_ptout(i).lt.Jetsmin(i,29)) Jetsmin(i,29)=Jt_ptout(i) If (Jtvpt(1,i).gt.Jetsmax(i,30)) Jetsmax(i,30)=Jtvpt(1,i) If (Jtvpt(1,i).lt.Jetsmin(i,30)) Jetsmin(i,30)=Jtvpt(1,i) If (Jtnovpt(i).gt.Jetsmax(i,31)) Jetsmax(i,31)=Jtnovpt(i) If (Jtnovpt(i).lt.Jetsmin(i,31)) Jetsmin(i,31)=Jtnovpt(i) Enddo CC C Jt1s Specific Variables CC If (REAL(Nj1).gt.Jt1smax(1,1)) Jt1smax(1,1)=REAL(Nj1) If (REAL(Nj1).lt.Jt1smin(1,1)) Jt1smin(1,1)=REAL(Nj1) Do i=1,nj1 If (J1et(i).gt.Jt1smax(i,2)) Jt1smax(i,2)= J1et(i) If (J1et(i).lt.Jt1smin(i,2)) Jt1smin(i,2)= J1et(i) If (J1m(i).gt.Jt1smax(i,3)) Jt1smax(i,3)= J1m(i) If (J1m(i).lt.Jt1smin(i,3)) Jt1smin(i,3)= J1m(i) If (J1p4(1,i).gt.Jt1smax(i,4)) Jt1smax(i,4)= J1p4(1,i) !x mom If (J1p4(1,i).lt.Jt1smin(i,4)) Jt1smin(i,4)= J1p4(1,i) !x mom If (J1eveta(i).gt.Jt1smax(i,5)) Jt1smax(i,5)= J1eveta(i) If (J1eveta(i).lt.Jt1smin(i,5)) Jt1smin(i,5)= J1eveta(i) If (J1dteta(i).gt.Jt1smax(i,6)) Jt1smax(i,6)= J1dteta(i) If (J1dteta(i).lt.Jt1smin(i,6)) Jt1smin(i,6)= J1dteta(i) If (J1phi(i).gt.Jt1smax(i,7)) Jt1smax(i,7)= J1phi(i) If (J1phi(i).lt.Jt1smin(i,7)) Jt1smin(i,7)= J1phi(i) If (REAL(J1word(i)).gt.Jt1smax(i,8)) Jt1smax(i,8)= REAL(J1word(i)) If (REAL(J1word(i)).lt.Jt1smin(i,8)) Jt1smin(i,8)= REAL(J1word(i)) If (REAL(J1jts(i)).gt.Jt1smax(i,9)) Jt1smax(i,9) = REAL(J1jts(i)) If (REAL(J1jts(i)).lt.Jt1smin(i,9)) Jt1smin(i,9) = REAL(J1jts(i)) If (J1corfm(i).gt.Jt1smax(i,10)) Jt1smax(i,10)= J1corfm(i) If (J1corfm(i).lt.Jt1smin(i,10)) Jt1smin(i,10)= J1corfm(i) If (J1corfd(i).gt.Jt1smax(i,11)) Jt1smax(i,11)=J1corfd(i) If (J1corfd(i).lt.Jt1smin(i,11)) Jt1smin(i,11)=J1corfd(i) If (J1egrd(i).gt.Jt1smax(i,12)) Jt1smax(i,12)= J1egrd(i) If (J1egrd(i).lt.Jt1smin(i,12)) Jt1smin(i,12)= J1egrd(i) Enddo CC C Met Specific Variables CC If (Met(5).gt.Metmax(1)) Metmax(1)= Met(5) If (Met(5).lt.Metmin(1)) Metmin(1)= Met(5) If (Metphi(5).gt.Metmax(2)) Metmax(2)=Metphi(5) If (Metphi(5).lt.Metmin(2)) Metmin(2)=Metphi(5) If (Metjt.gt.Metmax(3)) Metmax(3)= Metjt If (Metjt.lt.Metmin(3)) Metmin(3)= Metjt If (Metjtphi.gt.Metmax(4)) Metmax(4)= Metjtphi If (Metjtphi.lt.Metmin(4)) Metmin(4)= Metjtphi If (Metj1.gt.Metmax(5)) Metmax(5)= Metj1 If (Metj1.lt.Metmin(5)) Metmin(5)= Metj1 If (Metj1phi.gt.Metmax(6)) Metmax(6)= Metj1phi If (Metj1phi.lt.Metmin(6)) Metmin(6)= Metj1phi If (Metc.gt.Metmax(7)) Metmax(7)= Metc If (Metc.lt.Metmin(7)) Metmin(7)= Metc If (Metcphi.gt.Metmax(8)) Metmax(8)= Metcphi If (Metcphi.lt.Metmin(8)) Metmin(8)= Metcphi If (Htc.gt.Metmax(9)) Metmax(9)= Htc If (Htc.lt.Metmin(9)) Metmin(9)= Htc If (Metcf.gt.Metmax(10)) Metmax(10)= Metcf If (Metcf.lt.Metmin(10)) Metmin(10)= Metcf If (Metcfphi.gt.Metmax(11)) Metmax(11)=Metcfphi If (Metcfphi.lt.Metmin(11)) Metmin(11)=Metcfphi If (Htcf.gt.Metmax(12)) Metmax(12)= Htcf If (Htcf.lt.Metmin(12)) Metmin(12)= Htcf If (Etatrks.gt.Metmax(13)) Metmax(13)= Etatrks If (Etatrks.lt.Metmin(13)) Metmin(13)= Etatrks If (Phitrks.gt.Metmax(14)) Metmax(14)= Phitrks If (Phitrks.lt.Metmin(14)) Metmin(14)= Phitrks If (Pttrks.gt.Metmax(15)) Metmax(15)= Pttrks If (Pttrks.lt.Metmin(15)) Metmin(15)= Pttrks If (Seedpt.gt.Metmax(16)) Metmax(16)= Seedpt If (Seedpt.lt.Metmin(16)) Metmin(16)= Seedpt If (Sumet.gt.Metmax(17)) Metmax(17)= Sumet If (Sumet.lt.Metmin(17)) Metmin(17)= Sumet If (Sumetjet.gt.Metmax(18)) Metmax(18)=Sumetjet If (Sumetjet.lt.Metmin(18)) Metmin(18)=Sumetjet If (Metsig.gt.Metmax(19)) Metmax(19)= Metsig If (Metsig.lt.Metmin(19)) Metmin(19)= Metsig CC C Electron Specific Variables CC If (Real(Nel).gt.Elemax(1,1)) Elemax(1,1)=Real(Nel) If (Real(Nel).lt.Elemin(1,1)) Elemin(1,1)=Real(Nel) Do i=1,nel If (Ep4(1,i).gt.Elemax(i,2)) Elemax(1,2)= Ep4(1,i) !x mom If (Ep4(1,i).lt.Elemin(i,2)) Elemin(1,2)= Ep4(1,i) !x mom If (Eetcor(i).gt.Elemax(i,3)) Elemax(i,3)= Eetcor(i) If (Eetcor(i).lt.Elemin(i,3)) Elemin(i,3)= Eetcor(i) If (Eet(i).gt.Elemax(i,4)) Elemax(i,4)= Eet(i) If (Eet(i).lt.Elemin(i,4)) Elemin(i,4)= Eet(i) If (Echarge(i).gt.Elemax(i,5)) Elemax(i,5)= Echarge(i) If (Echarge(i).lt.Elemin(i,5)) Elemin(i,5)= Echarge(i) If (Edteta(i).gt.Elemax(i,6)) Elemax(i,6)= Edteta(i) If (Edteta(i).lt.Elemin(i,6)) Elemin(i,6)= Edteta(i) If (Eeveta(i).gt.Elemax(i,7)) Elemax(i,7)= Eeveta(i) If (Eeveta(i).lt.Elemin(i,7)) Elemin(i,7)= Eeveta(i) If (Ephi(i).gt.Elemax(i,8)) Elemax(i,8)= Ephi(i) If (Ephi(i).lt.Elemin(i,8)) Elemin(i,8)= Ephi(i) If (Ept(i).gt.Elemax(i,9)) Elemax(i,9)= Ept(i) If (Ept(i).lt.Elemin(i,9)) Elemin(i,9)= Ept(i) If (Estat(i).gt.Elemax(i,10)) Elemax(i,10)= Estat(i) If (Estat(i).lt.Elemin(i,10)) Elemin(i,10)= Estat(i) If (Ezv(i).gt.Elemax(i,11)) Elemax(i,11)= Ezv(i) If (Ezv(i).lt.Elemin(i,11)) Elemin(i,11)= Ezv(i) If (Eep(i).gt.Elemax(i,12)) Elemax(i,12)= Eep(i) If (Eep(i).lt.Elemin(i,12)) Elemin(i,12)= Eep(i) If (Ehadem(i).gt.Elemax(i,13)) Elemax(i,13)= Ehadem(i) If (Ehadem(i).lt.Elemin(i,13)) Elemin(i,13)= Ehadem(i) If (Exces(i).gt.Elemax(i,14)) Elemax(i,14)= Exces(i) If (Exces(i).lt.Elemin(i,14)) Elemin(i,14)= Exces(i) If (Edelx(i).gt.Elemax(i,15)) Elemax(i,15)= Edelx(i) If (Edelx(i).lt.Elemin(i,15)) Elemin(i,15)= Edelx(i) If (Ezces(i).gt.Elemax(i,16)) Elemax(i,16)= Ezces(i) If (Ezces(i).lt.Elemin(i,16)) Elemin(i,16)= Ezces(i) If (Edelz(i).gt.Elemax(i,17)) Elemax(i,17)= Edelz(i) If (Edelz(i).lt.Elemin(i,17)) Elemin(i,17)= Edelz(i) If (Echis(i).gt.Elemax(i,18)) Elemax(i,18)= Echis(i) If (Echis(i).lt.Elemin(i,18)) Elemin(i,18)= Echis(i) If (Echiw(i).gt.Elemax(i,19)) Elemax(i,19)= Echiw(i) If (Echiw(i).lt.Elemin(i,19)) Elemin(i,19)= Echiw(i) If (Eiso(i).gt.Elemax(i,20)) Elemax(i,20)= Eiso(i) If (Eiso(i).lt.Elemin(i,20)) Elemin(i,20)= Eiso(i) If (Etiso(i).gt.Elemax(i,21)) Elemax(i,21)= Etiso(i) If (Etiso(i).lt.Elemin(i,21)) Elemin(i,21)= Etiso(i) If (REAL(Eidwrd(i)).gt.Elemax(i,22)) Elemax(i,22)= REAL(Eidwrd(i)) If (REAL(Eidwrd(i)).lt.Elemin(i,22)) Elemin(i,22)= REAL(Eidwrd(i)) If (REAL(Eiswrd(i)).gt.Elemax(i,23)) Elemax(i,23)= REAL(Eiswrd(i)) If (REAL(Eiswrd(i)).lt.Elemin(i,23)) Elemin(i,23)= REAL(Eiswrd(i)) If (REAL(Econwrd(i)).gt.Elemax(i,24)) Elemax(i,24)= REAL(Econwrd(i)) If (REAL(Econwrd(i)).lt.Elemin(i,24)) Elemin(i,24)= REAL(Econwrd(i)) If (Eztrk(i).gt.Elemax(i,25)) Elemax(i,25)= Eztrk(i) If (Eztrk(i).lt.Elemin(i,25)) Elemin(i,25)= Eztrk(i) If (REAL(Efid(i)).gt.Elemax(i,26)) Elemax(i,26)= REAL(Efid(i)) If (REAL(Efid(i)).lt.Elemin(i,26)) Elemin(i,26)= REAL(Efid(i)) If (REAL(Edet(i)).gt.Elemax(i,27)) Elemax(i,27)= REAL(Edet(i)) If (REAL(Edet(i)).lt.Elemin(i,27)) Elemin(i,27)= REAL(Edet(i)) If (Ebte(i).gt.Elemax(i,28)) Elemax(i,28)= Ebte(i) If (Ebte(i).lt.Elemin(i,28)) Elemin(i,28)= Ebte(i) If (REAL(Enasl(i)).gt.Elemax(i,29)) Elemax(i,29)= REAL(Enasl(i)) If (REAL(Enasl(i)).lt.Elemin(i,29)) Elemin(i,29)= REAL(Enasl(i)) If (REAL(Enssl(i)).gt.Elemax(i,30)) Elemax(i,30)= REAL(Enssl(i)) If (REAL(Enssl(i)).lt.Elemin(i,30)) Elemin(i,30)= REAL(Enssl(i)) If (Eptrk(i).gt.Elemax(i,31)) Elemax(i,31)= Eptrk(i) If (Eptrk(i).lt.Elemin(i,31)) Elemin(i,31)= Eptrk(i) If (Elshr2(i).gt.Elemax(i,32)) Elemax(i,32)= Elshr2(i) If (Elshr2(i).lt.Elemin(i,32)) Elemin(i,32)= Elshr2(i) If (Echi3(i).gt.Elemax(i,33)) Elemax(i,33)= Echi3(i) If (Echi3(i).lt.Elemin(i,33)) Elemin(i,33)= Echi3(i) If (Echid(i).gt.Elemax(i,34)) Elemax(i,34)= Echid(i) If (Echid(i).lt.Elemin(i,34)) Elemin(i,34)= Echid(i) If (Evtxocc(i).gt.Elemax(i,35)) Elemax(i,35)= Evtxocc(i) If (Evtxocc(i).lt.Elemin(i,35)) Elemin(i,35)= Evtxocc(i) If (Etime(i).gt.Elemax(i,36)) Elemax(i,36)= Etime(i) If (Etime(i).lt.Elemin(i,36)) Elemin(i,36)= Etime(i) If (REAL(Etrind(i)).gt.Elemax(i,37)) Elemax(i,37)= REAL(Etrind(i)) If (REAL(Etrind(i)).lt.Elemin(i,37)) Elemin(i,37)= REAL(Etrind(i)) Enddo CC C Muon Specific Variables CC If (Real(Nmu).gt.Mumax(1,1)) Mumax(1,1)=Real(Nmu) If (Real(Nmu).lt.Mumin(1,1)) Mumin(1,1)=Real(Nmu) Do i=1,nmu If (Murawpt(i).gt.Mumax(i,2)) Mumax(i,2)= Murawpt(i) If (Murawpt(i).lt.Mumin(i,2)) Mumin(i,2)= Murawpt(i) If (Mup4(1,i).gt.Mumax(i,3)) Mumax(i,3)= Mup4(1,i) !x mom If (Mup4(1,i).lt.Mumin(i,3)) Mumin(i,3)= Mup4(1,i) !x mom If (Mupt(i).gt.Mumax(i,4)) Mumax(i,4)= Mupt(i) If (Mupt(i).lt.Mumin(i,4)) Mumin(i,4)= Mupt(i) If (Muq(i).gt.Mumax(i,5)) Mumax(i,5)= Muq(i) If (Muq(i).lt.Mumin(i,5)) Mumin(i,5)= Muq(i) If (Muhad(i).gt.Mumax(i,6)) Mumax(i,6)= Muhad(i) If (Muhad(i).lt.Mumin(i,6)) Mumin(i,6)= Muhad(i) If (Muem(i).gt.Mumax(i,7)) Mumax(i,7)= Muem(i) If (Muem(i).lt.Mumin(i,7)) Mumin(i,7)= Muem(i) If (Mueveta(i).gt.Mumax(i,8)) Mumax(i,8)= Mueveta(i) If (Mueveta(i).lt.Mumin(i,8)) Mumin(i,8)= Mueveta(i) If (Muphi(i).gt.Mumax(i,9)) Mumax(i,9)= Muphi(i) If (Muphi(i).lt.Mumin(i,9)) Mumin(i,9)= Muphi(i) If (Mustat(i).gt.Mumax(i,10)) Mumax(i,10)= Mustat(i) If (Mustat(i).lt.Mumin(i,10)) Mumin(i,10)= Mustat(i) If (Muzv(i).gt.Mumax(i,11)) Mumax(i,11)= Muzv(i) If (Muzv(i).lt.Mumin(i,11)) Mumin(i,11)= Muzv(i) If (Mud0(i).gt.Mumax(i,12)) Mumax(i,12)= Mud0(i) If (Mud0(i).lt.Mumin(i,12)) Mumin(i,12)= Mud0(i) If (Mudelz(i).gt.Mumax(i,13)) Mumax(i,13)= Mudelz(i) If (Mudelz(i).lt.Mumin(i,13)) Mumin(i,13)= Mudelz(i) If (Mudelx(i).gt.Mumax(i,14)) Mumax(i,14)= Mudelx(i) If (Mudelx(i).lt.Mumin(i,14)) Mumin(i,14)= Mudelx(i) If (Muiso(i).gt.Mumax(i,15)) Mumax(i,15)= Muiso(i) If (Muiso(i).lt.Mumin(i,15)) Mumin(i,15)= Muiso(i) If (Mutiso(i).gt.Mumax(i,16)) Mumax(i,16)= Mutiso(i) If (Mutiso(i).lt.Mumin(i,16)) Mumin(i,16)= Mutiso(i) If (REAL(Mucray(i)).gt.Mumax(i,17)) Mumax(i,17)= REAL(Mucray(i)) If (REAL(Mucray(i)).lt.Mumin(i,17)) Mumin(i,17)= REAL(Mucray(i)) If (REAL(Muidwrd(i)).gt.Mumax(i,18)) Mumax(i,18)= REAL(Muidwrd(i)) If (REAL(Muidwrd(i)).lt.Mumin(i,18)) Mumin(i,18)= REAL(Muidwrd(i)) If (REAL(Muiswrd(i)).gt.Mumax(i,19)) Mumax(i,19)= REAL(Muiswrd(i)) If (REAL(Muiswrd(i)).lt.Mumin(i,19)) Mumin(i,19)= REAL(Muiswrd(i)) If (Muztrk(i).gt.Mumax(i,20)) Mumax(i,20)= Muztrk(i) If (Muztrk(i).lt.Mumin(i,20)) Mumin(i,20)= Muztrk(i) If (REAL(Mufid(i)).gt.Mumax(i,21)) Mumax(i,21)= REAL(Mufid(i)) If (REAL(Mufid(i)).lt.Mumin(i,21)) Mumin(i,21)= REAL(Mufid(i)) If (REAL(Mudet(i)).gt.Mumax(i,22)) Mumax(i,22)= REAL(Mudet(i)) If (REAL(Mudet(i)).lt.Mumin(i,22)) Mumin(i,22)= REAL(Mudet(i)) If (Mubte(i).gt.Mumax(i,23)) Mumax(i,23)= Mubte(i) If (Mubte(i).lt.Mumin(i,23)) Mumin(i,23)= Mubte(i) If (REAL(Munasl(i)).gt.Mumax(i,24)) Mumax(i,24)= REAL(Munasl(i)) If (REAL(Munasl(i)).lt.Mumin(i,24)) Mumin(i,24)= REAL(Munasl(i)) If (REAL(Munssl(i)).gt.Mumax(i,25)) Mumax(i,25)= REAL(Munssl(i)) If (REAL(Munssl(i)).lt.Mumin(i,25)) Mumin(i,25)= REAL(Munssl(i)) If (Muchz0(i).gt.Mumax(i,26)) Mumax(i,26)= Muchz0(i) If (Muchz0(i).lt.Mumin(i,26)) Mumin(i,26)= Muchz0(i) If (Muchd0(i).gt.Mumax(i,27)) Mumax(i,27)= Muchd0(i) If (Muchd0(i).lt.Mumin(i,27)) Mumin(i,27)= Muchd0(i) If (Mutime(i).gt.Mumax(i,28)) Mumax(i,28)= Mutime(i) If (Mutime(i).lt.Mumin(i,28)) Mumin(i,28)= Mutime(i) If (Muvtxocc(i).gt.Mumax(i,29)) Mumax(i,29)= Muvtxocc(i) If (Muvtxocc(i).lt.Mumin(i,29)) Mumin(i,29)= Muvtxocc(i) If (Mudteta(i).gt.Mumax(i,30)) Mumax(i,30)= Mudteta(i) If (Mudteta(i).lt.Mumin(i,30)) Mumin(i,30)= Mudteta(i) If (REAL(Mutrind(i)).gt.Mumax(i,31)) Mumax(i,31)= REAL(Mutrind(i)) If (REAL(Mutrind(i)).lt.Mumin(i,31)) Mumin(i,31)= REAL(Mutrind(i)) Enddo CC C Tau Specific Variables CC If (Real(Ntau).gt.Taumax(1,1)) Taumax(1,1)=Real(Ntau) If (Real(Ntau).lt.Taumin(1,1)) Taumin(1,1)=Real(Ntau) Do i=1,ntau If (Tp4(1,i).gt.Taumax(i,2)) Taumax(i,2)= Tp4(1,i) !x mom If (Tp4(1,i).lt.Taumin(i,2)) Taumin(i,2)= Tp4(1,i) !x mom If (Tet(i).gt.Taumax(i,3)) Taumax(i,3)= Tet(i) If (Tet(i).lt.Taumin(i,3)) Taumin(i,3)= Tet(i) If (REAL(Tq(i)).gt.Taumax(i,4)) Taumax(i,4)= REAL(Tq(i)) If (REAL(Tq(i)).lt.Taumin(i,4)) Taumin(i,4)= REAL(Tq(i)) If (Tzv(i).gt.Taumax(i,5)) Taumax(i,5)= Tzv(i) If (Tzv(i).lt.Taumin(i,5)) Taumin(i,5)= Tzv(i) If (Tdteta(i).gt.Taumax(i,6)) Taumax(i,6)= Tdteta(i) If (Tdteta(i).lt.Taumin(i,6)) Taumin(i,6)= Tdteta(i) If (Teveta(i).gt.Taumax(i,7)) Taumax(i,7)= Teveta(i) If (Teveta(i).lt.Taumin(i,7)) Taumin(i,7)= Teveta(i) If (Tphi(i).gt.Taumax(i,8)) Taumax(i,8)= Tphi(i) If (Tphi(i).lt.Taumin(i,8)) Taumin(i,8)= Tphi(i) If (Tpt(i).gt.Taumax(i,9)) Taumax(i,9)= Tpt(i) If (Tpt(i).lt.Taumin(i,9)) Taumin(i,9)= Tpt(i) If (Tptvis(i).gt.Taumax(i,10)) Taumax(i,10)= Tptvis(i) If (Tptvis(i).lt.Taumin(i,10)) Taumin(i,10)= Tptvis(i) If (Tstat(i).gt.Taumax(i,11)) Taumax(i,11)= Tstat(i) If (Tstat(i).lt.Taumin(i,11)) Taumin(i,11)= Tstat(i) If (Tdelr(i).gt.Taumax(i,12)) Taumax(i,12)= Tdelr(i) If (Tdelr(i).lt.Taumin(i,12)) Taumin(i,12)= Tdelr(i) If (Temf(i).gt.Taumax(i,13)) Taumax(i,13)= Temf(i) If (Temf(i).lt.Taumin(i,13)) Taumin(i,13)= Temf(i) If (Tmass(i).gt.Taumax(i,14)) Taumax(i,14)= Tmass(i) If (Tmass(i).lt.Taumin(i,14)) Taumin(i,14)= Tmass(i) If (Tpttr(1,i).gt.Taumax(i,15)) Taumax(i,15)=Tpttr(1,i) !pt 1trk If (Tpttr(1,i).lt.Taumin(i,15)) Taumin(i,15)=Tpttr(1,i) !pt 1trk If (REAL(Tntr10(i)).gt.Taumax(i,16)) Taumax(i,16)= REAL(Tntr10(i)) If (REAL(Tntr10(i)).lt.Taumin(i,16)) Taumin(i,16)= REAL(Tntr10(i)) If (REAL(Tnpi0(i)).gt.Taumax(i,17)) Taumax(i,17)= REAL(Tnpi0(i)) If (REAL(Tnpi0(i)).lt.Taumin(i,17)) Taumin(i,17)= REAL(Tnpi0(i)) If (REAL(Tmustub(i)).gt.Taumax(i,18)) Taumax(i,18)= REAL(Tmustub(i)) If (REAL(Tmustub(i)).lt.Taumin(i,18)) Taumin(i,18)= REAL(Tmustub(i)) If (REAL(Ttrind(1,i)).gt.Taumax(i,19)) Taumax(i,19)= REAL(Ttrind(1,i)) !1 If (REAL(Ttrind(1,i)).lt.Taumin(i,19)) Taumin(i,19)= REAL(Ttrind(1,i)) !1 Enddo CC C Dilepton Specific Variables CC If (Real(Ndil).gt.Dilmax(1)) Dilmax(1)=Real(Ndil) If (Real(Ndil).lt.Dilmin(1)) Dilmin(1)=Real(Ndil) C If (Tzflag.gt.Dilmax(2)) Dilmax(2)= Tzflag C If (Tzflag.lt.Dilmin(2)) Dilmin(2)= Tzflag If (REAL(Dilwrd).gt.Dilmax(3)) Dilmax(3)= REAL(Dilwrd) If (REAL(Dilwrd).lt.Dilmin(3)) Dilmin(3)= REAL(Dilwrd) If (Metdilep.gt.Dilmax(4)) Dilmax(4)= Metdilep If (Metdilep.lt.Dilmin(4)) Dilmin(4)= Metdilep If (REAL(Njtdilep).gt.Dilmax(5)) Dilmax(5)= REAL(Njtdilep) If (REAL(Njtdilep).lt.Dilmin(5)) Dilmin(5)= REAL(Njtdilep) If (Dphiminj.gt.Dilmax(6)) Dilmax(6)= Dphiminj If (Dphiminj.lt.Dilmin(6)) Dilmin(6)= Dphiminj If (Dphiminl.gt.Dilmax(7)) Dilmax(7)= Dphiminl If (Dphiminl.lt.Dilmin(7)) Dilmin(7)= Dphiminl If (REAL(Dilmode).gt.Dilmax(8)) Dilmax(8)= REAL(Dilmode) If (REAL(Dilmode).lt.Dilmin(8)) Dilmin(8)= REAL(Dilmode) If (REAL(Diltype1).gt.Dilmax(9)) Dilmax(9)= REAL(Diltype1) If (REAL(Diltype1).lt.Dilmin(9)) Dilmin(9)= REAL(Diltype1) If (Eemass.gt.Dilmax(10)) Dilmax(10)= Eemass If (Eemass.lt.Dilmin(10)) Dilmin(10)= Eemass If (Mumumass.gt.Dilmax(11)) Dilmax(11)=Mumumass If (Mumumass.lt.Dilmin(11)) Dilmin(11)=Mumumass If (Taumass.gt.Dilmax(12)) Dilmax(12)= Taumass If (Taumass.lt.Dilmin(12)) Dilmin(12)= Taumass If (Phomass.gt.Dilmax(13)) Dilmax(13)= Phomass If (Phomass.lt.Dilmin(13)) Dilmin(13)= Phomass CC C Photon Specific Variables CC If (Real(Npho).gt.Phomax(1,1)) Phomax(1,1)=Real(Npho) If (Real(Npho).lt.Phomin(1,1)) Phomin(1,1)=Real(Npho) If (REAL(Ncencl).gt.Phomax(1,2)) Phomax(1,2)= REAL(Ncencl) If (REAL(Ncencl).lt.Phomin(1,2)) Phomin(1,2)= REAL(Ncencl) If (REAL(Nplugcl).gt.Phomax(1,3)) Phomax(1,3)= REAL(Nplugcl) If (REAL(Nplugcl).lt.Phomin(1,3)) Phomin(1,3)= REAL(Nplugcl) If (REAL(Nphoidl).gt.Phomax(1,4)) Phomax(1,4)= REAL(Nphoidl) If (REAL(Nphoidl).lt.Phomin(1,4)) Phomin(1,4)= REAL(Nphoidl) If (REAL(Nphoisol).gt.Phomax(1,5)) Phomax(1,5)= REAL(Nphoisol) If (REAL(Nphoisol).lt.Phomin(1,5)) Phomin(1,5)= REAL(Nphoisol) If (REAL(Nphoid).gt.Phomax(1,6)) Phomax(1,6)= REAL(Nphoid) If (REAL(Nphoid).lt.Phomin(1,6)) Phomin(1,6)= REAL(Nphoid) If (REAL(Nphoisol).gt.Phomax(1,7)) Phomax(1,7)= REAL(Nphoisol) If (REAL(Nphoisol).lt.Phomin(1,7)) Phomin(1,7)= REAL(Nphoisol) If (REAL(Ncmus).gt.Phomax(1,8)) Phomax(1,8)= REAL(Ncmus) If (REAL(Ncmus).lt.Phomin(1,8)) Phomin(1,8)= REAL(Ncmus) If (REAL(Costub).gt.Phomax(1,9)) Phomax(1,9)= REAL(Costub) If (REAL(Costub).lt.Phomin(1,9)) Phomin(1,9)= REAL(Costub) Do i=1,npho If (REAL(Pdetect(i)).gt.Phomax(i,10)) Phomax(i,10)= REAL(Pdetect(i)) If (REAL(Pdetect(i)).lt.Phomin(i,10)) Phomin(i,10)= REAL(Pdetect(i)) If (Pzv(i).gt.Phomax(i,11)) Phomax(i,11)= Pzv(i) If (Pzv(i).lt.Phomin(i,11)) Phomin(i,11)= Pzv(i) If (Pet(i).gt.Phomax(i,12)) Phomax(i,12)= Pet(i) If (Pet(i).lt.Phomin(i,12)) Phomin(i,12)= Pet(i) If (Petc(i).gt.Phomax(i,13)) Phomax(i,13)= Petc(i) If (Petc(i).lt.Phomin(i,13)) Phomin(i,13)= Petc(i) If (Pp4(1,i).gt.Phomax(i,14)) Phomax(i,14)= Pp4(1,i) !x mom If (Pp4(1,i).lt.Phomin(i,14)) Phomin(i,14)= Pp4(i,i) If (Pco4(i).gt.Phomax(i,15)) Phomax(i,15)= Pco4(i) If (Pco4(i).lt.Phomin(i,15)) Phomin(i,15)= Pco4(i) If (Pcesx(i).gt.Phomax(i,16)) Phomax(i,16)= Pcesx(i) If (Pcesx(i).lt.Phomin(i,16)) Phomin(i,16)= Pcesx(i) If (Pcesz(i).gt.Phomax(i,17)) Phomax(i,17)= Pcesz(i) If (Pcesz(i).lt.Phomin(i,17)) Phomin(i,17)= Pcesz(i) If (Pcpr5ph(i).gt.Phomax(i,18)) Phomax(i,18)= Pcpr5ph(i) If (Pcpr5ph(i).lt.Phomin(i,18)) Phomin(i,18)= Pcpr5ph(i) If (Pcpr5ps(i).gt.Phomax(i,19)) Phomax(i,19)= Pcpr5ps(i) If (Pcpr5ps(i).lt.Phomin(i,19)) Phomin(i,19)= Pcpr5ps(i) If (Pdteta(i).gt.Phomax(i,20)) Phomax(i,20)= Pdteta(i) If (Pdteta(i).lt.Phomin(i,20)) Phomin(i,20)= Pdteta(i) If (Peveta(i).gt.Phomax(i,21)) Phomax(i,21)= Peveta(i) If (Peveta(i).lt.Phomin(i,21)) Phomin(i,21)= Peveta(i) If (Pcese(i).gt.Phomax(i,22)) Phomax(i,22)= Pcese(i) If (Pcese(i).lt.Phomin(i,22)) Phomin(i,22)= Pcese(i) If (Pphi(i).gt.Phomax(i,23)) Phomax(i,23)= Pphi(i) If (Pphi(i).lt.Phomin(i,23)) Phomin(i,23)= Pphi(i) If (Pco7(i).gt.Phomax(i,24)) Phomax(i,24)= Pco7(i) If (Pco7(i).lt.Phomin(i,24)) Phomin(i,24)= Pco7(i) If (Phadem(i).gt.Phomax(i,25)) Phomax(i,25)= Phadem(i) If (Phadem(i).lt.Phomin(i,25)) Phomin(i,25)= Phadem(i) If (Plshr(i).gt.Phomax(i,26)) Phomax(i,26)= Plshr(i) If (Plshr(i).lt.Phomin(i,26)) Phomin(i,26)= Plshr(i) If (Pstr2(i).gt.Phomax(i,27)) Phomax(i,27)= Pstr2(i) If (Pstr2(i).lt.Phomin(i,27)) Phomin(i,27)= Pstr2(i) If (Pstre2(i).gt.Phomax(i,28)) Phomax(i,28)= Pstre2(i) If (Pstre2(i).lt.Phomin(i,28)) Phomin(i,28)= Pstre2(i) If (Pwir2(i).gt.Phomax(i,29)) Phomax(i,29)= Pwir2(i) If (Pwir2(i).lt.Phomin(i,29)) Phomin(i,29)= Pwir2(i) If (Pwire2(i).gt.Phomax(i,30)) Phomax(i,30)= Pwire2(i) If (Pwire2(i).lt.Phomin(i,30)) Phomin(i,30)= Pwire2(i) If (Psumpt4(i).gt.Phomax(i,31)) Phomax(i,31)= Psumpt4(i) If (Psumpt4(i).lt.Phomin(i,31)) Phomin(i,31)= Psumpt4(i) If (Pvtxocc(i).gt.Phomax(i,32)) Phomax(i,32)= Pvtxocc(i) If (Pvtxocc(i).lt.Phomin(i,32)) Phomin(i,32)= Pvtxocc(i) If (Ppt(i).gt.Phomax(i,33)) Phomax(i,33)= Ppt(i) If (Ppt(i).lt.Phomin(i,33)) Phomin(i,33)= Ppt(i) If (REAL(Pn3d(i)).gt.Phomax(i,34)) Phomax(i,34)= REAL(Pn3d(i)) If (REAL(Pn3d(i)).lt.Phomin(i,34)) Phomin(i,34)= REAL(Pn3d(i)) If (Pchi(i).gt.Phomax(i,35)) Phomax(i,35)= Pchi(i) If (Pchi(i).lt.Phomin(i,35)) Phomin(i,35)= Pchi(i) If (Pchi3x3(i).gt.Phomax(i,36)) Phomax(i,36)= Pchi3x3(i) If (Pchi3x3(i).lt.Phomin(i,36)) Phomin(i,36)= Pchi3x3(i) If (Pchieta(i).gt.Phomax(i,37)) Phomax(i,37)= Pchieta(i) If (Pchieta(i).lt.Phomin(i,37)) Phomin(i,37)= Pchieta(i) If (Pchiphi(i).gt.Phomax(i,38)) Phomax(i,38)= Pchiphi(i) If (Pchiphi(i).lt.Phomin(i,38)) Phomin(i,38)= Pchiphi(i) If (Pstat(i).gt.Phomax(i,39)) Phomax(i,39)= Pstat(i) If (Pstat(i).lt.Phomin(i,39)) Phomin(i,39)= Pstat(i) If (Pboxiso(i).gt.Phomax(i,40)) Phomax(i,40)= Pboxiso(i) If (Pboxiso(i).lt.Phomin(i,40)) Phomin(i,40)= Pboxiso(i) C If (Pcesflg(i).gt.Phomax(i,41)) Phomax(i,41)= Pcesflg(i) C If (Pcesflg(i).lt.Phomin(i,41)) Phomin(i,41)= Pcesflg(i) If (Pe(i).gt.Phomax(i,42)) Phomax(i,42)= Pe(i) If (Pe(i).lt.Phomin(i,42)) Phomin(i,42)= Pe(i) If (Pceswht(i).gt.Phomax(i,43)) Phomax(i,43)= Pceswht(i) If (Pceswht(i).lt.Phomin(i,43)) Phomin(i,43)= Pceswht(i) If (Pcprwht(i).gt.Phomax(i,44)) Phomax(i,44)= Pcprwht(i) If (Pcprwht(i).lt.Phomin(i,44)) Phomin(i,44)= Pcprwht(i) If (Pcesslide(i).gt.Phomax(i,45)) Phomax(i,45)= Pcesslide(i) If (Pcesslide(i).lt.Phomin(i,45)) Phomin(i,45)= Pcesslide(i) If (REAL(Pwrd(i)).gt.Phomax(i,46)) Phomax(i,46)= REAL(Pwrd(i)) If (REAL(Pwrd(i)).lt.Phomin(i,46)) Phomin(i,46)= REAL(Pwrd(i)) If (Psth(i).gt.Phomax(i,47)) Phomax(i,47)= Psth(i) If (Psth(i).lt.Phomin(i,47)) Phomin(i,47)= Psth(i) If (REAL(Pelind(i)).gt.Phomax(i,48)) Phomax(i,48)= REAL(Pelind(i)) If (REAL(Pelind(i)).lt.Phomin(i,48)) Phomin(i,48)= REAL(Pelind(i)) If (Pchistr(i).gt.Phomax(i,49)) Phomax(i,49)= Pchistr(i) If (Pchistr(i).lt.Phomin(i,49)) Phomin(i,49)= Pchistr(i) If (Pchiwir(i).gt.Phomax(i,50)) Phomax(i,50)= Pchiwir(i) If (Pchiwir(i).lt.Phomin(i,50)) Phomin(i,50)= Pchiwir(i) If (Pcessig(i).gt.Phomax(i,51)) Phomax(i,51)= Pcessig(i) If (Pcessig(i).lt.Phomin(i,51)) Phomin(i,51)= Pcessig(i) If (Pcespg(i).gt.Phomax(i,52)) Phomax(i,52)= Pcespg(i) If (Pcespg(i).lt.Phomin(i,52)) Phomin(i,52)= Pcespg(i) If (Pcespb(i).gt.Phomax(i,53)) Phomax(i,53)= Pcespb(i) If (Pcespb(i).lt.Phomin(i,53)) Phomin(i,53)= Pcespb(i) If (Pcprpg(i).gt.Phomax(i,54)) Phomax(i,54)= Pcprpg(i) If (Pcprpg(i).lt.Phomin(i,54)) Phomin(i,54)= Pcprpg(i) If (Pcprpb(i).gt.Phomax(i,55)) Phomax(i,55)= Pcprpb(i) If (Pcprpb(i).lt.Phomin(i,55)) Phomin(i,55)= Pcprpb(i) If (Ppt2(i).gt.Phomax(i,56)) Phomax(i,56)= Ppt2(i) If (Ppt2(i).lt.Phomin(i,56)) Phomin(i,56)= Ppt2(i) If (REAL(Pntrwd3(i)).gt.Phomax(i,57)) Phomax(i,57)= REAL(Pntrwd3(i)) If (REAL(Pntrwd3(i)).lt.Phomin(i,57)) Phomin(i,57)= REAL(Pntrwd3(i)) If (Pcescprx(i).gt.Phomax(i,58)) Phomax(i,58)= Pcescprx(i) If (Pcescprx(i).lt.Phomin(i,58)) Phomin(i,58)= Pcescprx(i) If (Ptkcprx(i).gt.Phomax(i,59)) Phomax(i,59)= Ptkcprx(i) If (Ptkcprx(i).lt.Phomin(i,59)) Phomin(i,59)= Ptkcprx(i) If (Ptime(i).gt.Phomax(i,60)) Phomax(i,60)= Ptime(i) If (Ptime(i).lt.Phomin(i,60)) Phomin(i,60)= Ptime(i) If (Pvcor(i).gt.Phomax(i,61)) Phomax(i,61)= Pvcor(i) If (Pvcor(i).lt.Phomin(i,61)) Phomin(i,61)= Pvcor(i) If (Plcor(i).gt.Phomax(i,62)) Phomax(i,62)= Plcor(i) If (Plcor(i).lt.Phomin(i,62)) Phomin(i,62)= Plcor(i) Enddo CC C Tagging Specific Variables CC If (Real(Ntag).gt.Tagmax(1,1)) Tagmax(1,1)=Real(Ntag) If (Real(Ntag).lt.Tagmin(1,1)) Tagmin(1,1)=Real(Ntag) If (Real(Btgwrd).gt.Tagmax(1,2)) Tagmax(1,2)=Real(Btgwrd) If (Real(Btgwrd).lt.Tagmin(1,2)) Tagmin(1,2)=Real(Btgwrd) Do i=1,ntag If (Lxy(i).gt.Tagmax(i,3)) Tagmax(i,3)= Lxy(i) If (Lxy(i).lt.Tagmin(i,3)) Tagmin(i,3)= Lxy(i) If (Lxyerr(i).gt.Tagmax(i,4)) Tagmax(i,4)=Lxyerr(i) If (Lxyerr(i).lt.Tagmin(i,4)) Tagmin(i,4)=Lxyerr(i) If (Ctau(i).gt.Tagmax(i,5)) Tagmax(i,5)= Ctau(i) If (Ctau(i).lt.Tagmin(i,5)) Tagmin(i,5)= Ctau(i) If (Tagchi(i).gt.Tagmax(i,6)) Tagmax(i,6)=Tagchi(i) If (Tagchi(i).lt.Tagmin(i,6)) Tagmin(i,6)=Tagchi(i) Enddo CC C Vertex Specific Variables CC If (Real(Npvert).gt.Vermax(1,1)) Vermax(1,1)=Real(Npvert) If (Real(Npvert).lt.Vermin(1,1)) Vermin(1,1)=Real(Npvert) If (Vxprim(1).gt.Vermax(1,2)) Vermax(1,2)= Vxprim(1) !x prim vert If (Vxprim(1).lt.Vermin(1,2)) Vermin(1,2)= Vxprim(1) !x prim vert Do i=1,npvert If (REAL(Vclass(i)).gt.Vermax(i,3)) Vermax(i,3)= REAL(Vclass(i)) If (REAL(Vclass(i)).lt.Vermin(i,3)) Vermin(i,3)= REAL(Vclass(i)) If (Vz(i).gt.Vermax(i,4)) Vermax(i,4)= Vz(i) If (Vz(i).lt.Vermin(i,4)) Vermin(i,4)= Vz(i) If (Vq(i).gt.Vermax(i,5)) Vermax(i,5)= Vq(i) If (Vq(i).lt.Vermin(i,5)) Vermin(i,5)= Vq(i) Enddo If (Ptnov.gt.Vermax(i,6)) Vermax(1,6)= Ptnov If (Ptnov.lt.Vermin(i,6)) Vermin(1,6)= Ptnov If (Vzev.gt.Vermax(i,7)) Vermax(1,7)= Vzev If (Vzev.lt.Vermin(i,7)) Vermin(1,7)= Vzev If (REAL(Vclassev).gt.Vermax(i,8)) Vermax(1,8)=REAL(Vclassev) If (REAL(Vclassev).lt.Vermin(i,8)) Vermin(1,8)=REAL(Vclassev) CC C Track Specific Variables CC If (Real(Ntr).gt.Trkmax(1,1)) Trkmax(1,1)=Real(Ntr) If (Real(Ntr).lt.Trkmin(1,1)) Trkmin(1,1)=Real(Ntr) Do i=1,ntr If (REAL(Trvind(i)).gt.Trkmax(i,2)) Trkmax(i,2)= REAL(Trvind(i)) If (REAL(Trvind(i)).lt.Trkmin(i,2)) Trkmin(i,2)= REAL(Trvind(i)) If (REAL(Trnsvx(i)).gt.Trkmax(i,3)) Trkmax(i,3)= REAL(Trnsvx(i)) If (REAL(Trnsvx(i)).lt.Trkmin(i,3)) Trkmin(i,3)= REAL(Trnsvx(i)) If (REAL(Trnasl(i)).gt.Trkmax(i,4)) Trkmax(i,4)= REAL(Trnasl(i)) If (REAL(Trnasl(i)).lt.Trkmin(i,4)) Trkmin(i,4)= REAL(Trnasl(i)) If (REAL(Trnssl(i)).gt.Trkmax(i,5)) Trkmax(i,5)= REAL(Trnssl(i)) If (REAL(Trnssl(i)).lt.Trkmin(i,5)) Trkmin(i,5)= REAL(Trnssl(i)) If (Trcot(i).gt.Trkmax(i,6)) Trkmax(i,6)= Trcot(i) If (Trcot(i).lt.Trkmin(i,6)) Trkmin(i,6)= Trcot(i) If (Trcurv(i).gt.Trkmax(i,7)) Trkmax(i,7)= Trcurv(i) If (Trcurv(i).lt.Trkmin(i,7)) Trkmin(i,7)= Trcurv(i) If (Trz0(i).gt.Trkmax(i,8)) Trkmax(i,8)= Trz0(i) If (Trz0(i).lt.Trkmin(i,8)) Trkmin(i,8)= Trz0(i) If (Trd0(i).gt.Trkmax(i,9)) Trkmax(i,9)= Trd0(i) If (Trd0(i).lt.Trkmin(i,9)) Trkmin(i,9)= Trd0(i) If (Trq(i).gt.Trkmax(i,10)) Trkmax(i,10)= Trq(i) If (Trq(i).lt.Trkmin(i,10)) Trkmin(i,10)= Trq(i) If (Trphi0(i).gt.Trkmax(i,11)) Trkmax(i,11)= Trphi0(i) If (Trphi0(i).lt.Trkmin(i,11)) Trkmin(i,11)= Trphi0(i) If (Trp4(1,i).gt.Trkmax(i,12)) Trkmax(i,12)= Trp4(1,i) !x only If (Trp4(1,i).lt.Trkmin(i,12)) Trkmin(i,12)= Trp4(1,i) !x only If (Trpt(i).gt.Trkmax(i,13)) Trkmax(i,13)= Trpt(i) If (Trpt(i).lt.Trkmin(i,13)) Trkmin(i,13)= Trpt(i) If (Treta(i).gt.Trkmax(i,14)) Trkmax(i,14)= Treta(i) If (Treta(i).lt.Trkmin(i,14)) Trkmin(i,14)= Treta(i) If (Trexeta(i).gt.Trkmax(i,15)) Trkmax(i,15)=Trexeta(i) If (Trexeta(i).lt.Trkmin(i,15)) Trkmin(i,15)=Trexeta(i) If (Trexphi(i).gt.Trkmax(i,16)) Trkmax(i,16)=Trexphi(i) If (Trexphi(i).lt.Trkmin(i,16)) Trkmin(i,16)=Trexphi(i) If (REAL(Trbcvind(i)).gt.Trkmax(i,17)) Trkmax(i,17)=REAL(Trbcvind(i)) If (REAL(Trbcvind(i)).lt.Trkmin(i,17)) Trkmin(i,17)=REAL(Trbcvind(i)) If (Trbcz0(i).gt.Trkmax(i,18)) Trkmax(i,18)= Trbcz0(i) If (Trbcz0(i).lt.Trkmin(i,18)) Trkmin(i,18)= Trbcz0(i) If (Trbcd0(i).gt.Trkmax(i,19)) Trkmax(i,19)= Trbcd0(i) If (Trbcd0(i).lt.Trkmin(i,19)) Trkmin(i,19)= Trbcd0(i) If (Trbcq(i).gt.Trkmax(i,20)) Trkmax(i,20)= Trbcq(i) If (Trbcq(i).lt.Trkmin(i,20)) Trkmin(i,20)= Trbcq(i) If (Trbcp4(1,i).gt.Trkmax(i,21)) Trkmax(i,21)=Trbcp4(1,i) !1 only If (Trbcp4(1,i).lt.Trkmin(i,21)) Trkmin(i,21)=Trbcp4(1,i) !1 only If (Trbcpt(i).gt.Trkmax(i,22)) Trkmax(i,22)= Trbcpt(i) If (Trbcpt(i).lt.Trkmin(i,22)) Trkmin(i,22)= Trbcpt(i) If (Trbceta(i).gt.Trkmax(i,23)) Trkmax(i,23)= Trbceta(i) If (Trbceta(i).lt.Trkmin(i,23)) Trkmin(i,23)= Trbceta(i) If (Triso4(i).gt.Trkmax(i,24)) Trkmax(i,24)= Triso4(i) If (Triso4(i).lt.Trkmin(i,24)) Trkmin(i,24)= Triso4(i) Enddo return end C an event selection filter before filling histograms C*************************************************************************** real function histfilt(xdummy) implicit none include 'cwn5.inc' c include 'c$inc:pivalu.inc' integer xdummy integer i,j,k integer nphotight,njethad,phoindex(10),jetindex(20) real deltar,etadiff,phidiff real cone_size,ht,phometjt,metx,mety,phometphi logical jethad integer vindf real vqf,vzf real phet C currently a dummy routine that passes all events histfilt = 1.0 C find tight photons return end