Received: from d0sgi0.fnal.gov (d0sgi0.fnal.gov [131.225.221.44]) by d0sgim.fnal.gov (980427.SGI.8.8.8/970903.SGI.AUTOCF) via SMTP id KAA13964 for <\meena@d0sgim.fnal.gov>; Tue, 22 Dec 1998 10:42:24 -0600 (CST) Received: from FNAL.FNAL.Gov (fnal.fnal.gov [131.225.9.8]) by d0sgi0.fnal.gov (950413.SGI.8.6.12/950213.SGI.AUTOCF) via ESMTP id KAA08476 for ; Tue, 22 Dec 1998 10:42:28 -0600 Received: from d0tng.fnal.gov ("port 2632"@d0tng.fnal.gov) by FNAL.FNAL.GOV (PMDF V5.1-12 #3998) with SMTP id <01J5MW2SOAOK0004QR@FNAL.FNAL.GOV> for meena@d0sgi0.fnal.gov; Tue, 22 Dec 1998 10:43:04 -0600 CDT Date: Tue, 22 Dec 1998 10:43:02 -0600 From: bornali Subject: Please take a look in the Fortran file To: meena@fnal.gov Cc: BORNALI@dasb01 Message-id: <981222104302.22600404@DASB01.FNAL.GOV> Hi Meena, I made the changes for the calculation of the pt of the STT tracks. Please take a look at the fortran file. Please let me know if this is the right way to do . Thanks, Bornali _______________________________________________________________________________ REAL FUNCTION full_trigmenu(xdummy) * * implicit none * LOGICAL chain CHARACTER*128 cfile INTEGER idnevt,nchevt,ichevt REAL vidn1,vidn2,vidn3,vidn(10) * COMMON /pawidn/ idnevt,vidn1,vidn2,vidn3,vidn COMMON /pawchn/ chain, nchevt, ichevt COMMON /pawchc/ cfile * *-- Ntuple Variable Declarations * INTEGER run,evt,max_isa_id,nmct,id(500),nttrk(500),ifib1(500) + ,ifib16(500),idvert(500),isatrk_id(500),ntrig,itfcft1(500) + ,itfcft4(500),itfcft5(500),itfcft8(500),itfcft9(500) + ,itfcft12(500),itfcft13(500),itfcft16(500),itpatset(500) + ,itpat(500),itf1st(500),nmctrk(500),cft_trkid(500) + ,cft_trkmul(500),idtrack(500),svx_tracknum,svx_clusters(100) + ,svx_hits_used(100),svx_status(100),svx_barrels(100) + ,svx_nclust(100),svx_trackid(100),svx_truehit(100) + ,svx_samehit(100),svx_mult(100),svx_idrop(100),ntcps,ntfps + ,phiseg(1000),seedu(1000),seedv(1000),ps_ieta(1000), + ps_iphi(1000) + ,nt_nelec,nt_neclus,mcntele,mcidele(10),idparent(10),nt_njet + ,nt_njclus,nt_njseeds,nt_rjseeds(50),nt_ajseeds(50),nt_neseeds + ,nt_reseeds(50),nt_aeseeds(50),mu_version,mtm_ct1,mtm_ct2, + mtm_ct3 + ,mtm_ct4,mtm_cm1,mtm_cm2,mtm_cm3,mtm_cm4,mtm_cl1,mtm_cl2,mtm_cl3 + ,mtm_cl4,mtm_at1,mtm_at2,mtm_at3,mtm_at4,mtm_am1,mtm_am2,mtm_am3 + ,mtm_am4,mtm_al1,mtm_al2,mtm_al3,mtm_al4,mtm_xt1,mtm_xt2,mtm_xt3 + ,mtm_xt4,mtm_xm1,mtm_xm2,mtm_xm3,mtm_xm4,mtm_xl1,mtm_xl2,mtm_xl3 + ,mtm_xl4,nalg,idalg(20),maxitalg(20),njetsalg(20),npjet + ,flavpjet(100) REAL weight,xvert(500),yvert(500),zvert(500),eta(500),phi(500) + ,pt(500),zz1(500),zz16(500),ttmom1(500),ttmom2(500),hitfrac(16) + ,svx_kc(100),svx_phi0c(100),svx_sipt(100),svx_cftpt(100) + ,svx_b(100),svx_phi0(100),svx_tk(100),svx_chisqr(100) + ,svx_sichsqr(100),svx_cftchsqr(100),svx_hitr(8,100) + ,svx_hitphi(8,100),svx_chsqr_new(100),svx_b_new(100) + ,svx_tk_new(100),svx_phi0_new(100),phipos(1280),resppos(1280) + ,phineg(1280),respneg(1280),e_seedu(1000),e_seedv(1000) + ,e_isou(1000),e_isov(1000),e_confu(1000),e_confv(1000) + ,e_caltwr(1000),nt_elet(30),nt_eleta(30),nt_elphi(30), + nt_emfr(30) + ,nt_isofr(30),nt_el_time,mcetele(10),mcetaele(10),mcphiele(10) + ,nt_jtet(30),nt_jteta(30),nt_jtphi(30),nt_jt_time,nt_met,nt_metx + ,nt_mety,nt_metphi,nt_met_time,nt_ejseeds(50),nt_eeseeds(50) + ,clus_set,clus_met,clu_metphi,dralg(20),etalg(20),flagalg(20) + ,pjet1alg(20),etpjet(100),epjet(100),phipjet(100),etapjet(100) + ,flagpjet(100) INTEGER wmodes(2),zmodes(2),hmode * COMMON /pawcr4/ run,evt,weight,max_isa_id,nmct,id,nttrk,ifib1 + ,ifib16,xvert,yvert,zvert,eta,phi,pt,zz1,zz16,idvert,isatrk_id + ,ntrig,itfcft1,itfcft4,itfcft5,itfcft8,itfcft9,itfcft12,itfcft13 + ,itfcft16,itpatset,itpat,itf1st,ttmom1,ttmom2,nmctrk,cft_trkid + ,cft_trkmul,idtrack,hitfrac,svx_tracknum,svx_clusters,svx_kc + ,svx_phi0c,svx_sipt,svx_cftpt,svx_b,svx_phi0,svx_tk,svx_chisqr + ,svx_hits_used,svx_status,svx_barrels,svx_sichsqr,svx_cftchsqr + ,svx_hitr,svx_hitphi,svx_nclust,svx_trackid,svx_truehit + ,svx_samehit,svx_mult,svx_chsqr_new,svx_b_new,svx_tk_new + ,svx_phi0_new,svx_idrop,ntcps,phipos,resppos,phineg,respneg, + ntfps + ,phiseg,seedu,seedv,ps_ieta,ps_iphi,e_seedu,e_seedv,e_isou, + e_isov + ,e_confu,e_confv,e_caltwr,nt_nelec,nt_elet,nt_eleta,nt_elphi + ,nt_emfr,nt_isofr,nt_neclus,nt_el_time,mcntele,mcetele,mcetaele + ,mcphiele,mcidele,idparent,nt_njet,nt_jtet,nt_jteta,nt_jtphi + ,nt_njclus,nt_jt_time,nt_met,nt_metx,nt_mety,nt_metphi + ,nt_met_time,nt_njseeds,nt_ejseeds,nt_rjseeds,nt_ajseeds + ,nt_neseeds,nt_eeseeds,nt_reseeds,nt_aeseeds,clus_set,clus_met + ,clu_metphi,mu_version,mtm_ct1,mtm_ct2,mtm_ct3,mtm_ct4,mtm_cm1 + ,mtm_cm2,mtm_cm3,mtm_cm4,mtm_cl1,mtm_cl2,mtm_cl3,mtm_cl4,mtm_at1 + ,mtm_at2,mtm_at3,mtm_at4,mtm_am1,mtm_am2,mtm_am3,mtm_am4,mtm_al1 + ,mtm_al2,mtm_al3,mtm_al4,mtm_xt1,mtm_xt2,mtm_xt3,mtm_xt4,mtm_xm1 + ,mtm_xm2,mtm_xm3,mtm_xm4,mtm_xl1,mtm_xl2,mtm_xl3,mtm_xl4,nalg + ,idalg,dralg,etalg,maxitalg,flagalg,njetsalg,pjet1alg,npjet + ,etpjet,epjet,phipjet,etapjet,flagpjet,flavpjet,wmodes,zmodes + ,hmode * * * vector vtrig(80) * real vtrig(80) LOGICAL first DATA first /.true./ SAVE first SAVE IFE,ife2,ifm INTEGER i,j,k,ieta,jeta,jphi Integer ife,ife2,ifm INTEGER nl1muALL(4),nl1jt(7),nl1lt(7),nl1em(5),nl1cft(4), & nl1tagcft(4),nl1lj(7), NL1MUCENTRAL(4) REAL thl1jt(7),thl1lt(7),thl1em(5),thl1cft(4) DATA thl1jt/3.,5.,7.,10.,12.,14.,16./ DATA thl1em/3., 5., 7., 10., 12. / DATA thl1cft / 1.5, 3., 6., 8. / DATA thl1lt /5., 10., 15., 20., 25., 35., 45./ REAL cftphi,thl1cps,delphipos,delphineg DATA thl1cps/6.0/ real svx_sipt_new(4) REAL SIPT(4) DATA SIPT /1.5,3.0,6.0,8.0/ INTEGER NTRKSIPT(4) REAL rtemp,rft16,rft1,trakpt(500),en(10,4) REAL EMSCALE,HSCALE,MISSCALE LOGICAL tagcft(500),nl2elec(30),l2elec_can(20),tagcftp,tagcftn * **** EM_JET **** LOGICAL l1em_jet,l2em_jet REAL nl1jt_et(7), nl1em_et(5) REAL thl2emjet_elec,thl2emjet_jet INTEGER nl2emjet_elec,nl2emjet_jet DATA thl2emjet_elec,thl2emjet_jet/15.,10./ **** EM_MU **** LOGICAL l1em_mu,l2em_mu REAL thl2emmu_elec INTEGER nl2emmu_elec DATA thl2emmu_elec/7./ **** MISSING_ET **** LOGICAL l1missing_et,l2missing_et REAL thl1miss_ET,thl2miss_met,thl2miss_jet INTEGER nl2miss_et DATA thl1miss_ET /40./ DATA thl2miss_met /40./ DATA thl2miss_jet /10./ **** JET_3_MISS **** LOGICAL l1jet_3_miss, l2jet_3_miss REAL thl1jet_3_miss_et,thl2jet_3_miss,thl2met_3_miss INTEGER nl2jet_3_miss DATA thl2jet_3_miss /20./ DATA thl2met_3_miss /25./ DATA thl1jet_3_miss_et/20./ **** MU_MAX **** LOGICAL l1mu_1_high,l2mu_1_high **** MU_2_MED **** LOGICAL l1mu_2_med,l2mu_2_med **** MU_JET **** LOGICAL l1mu_jet_high,l2mu_jet_high INTEGER nl1mu_jet_high, nl2mu_jet_high REAL thl2mu_jet_high DATA thl2mu_jet_high /10./ **** JET_MULTI **** LOGICAL l1jet_multi,l2jet_multi REAL thl2jt_multi INTEGER nl2jt_multi DATA thl2jt_multi / 10. / **** JET_MAX **** LOGICAL l1jet_max, l2jet_max REAL thl2jet_max INTEGER nl2jmax DATA thl2jet_max/85./ **** JET_2_MED **** LOGICAL L1JET_2_MED,L2jet_2_med REAL THL2JET_2_MED INTEGER NL2JET_2_MED DATA THL2JET_2_MED/20./ **** JET_2_MIN_CFT3 **** LOGICAL L1JeT_2_MIN_CFT3,L2JeT_2_MIN_CFT3 REAL THL2JET_2_MIN_CFT3 INTEGER NL2JET_2_MIN_CFT3 DATA THL2JET_2_min_CFT3/10./ **** JET_3_MED **** LOGICAL L1JET_3_MED,L2jet_3_med REAL THL2JET_3_MED INTEGER NL2JET_3_MED DATA THL2JET_3_MED/10./ **** TEST **** LOGICAL l1em_1_high LOGICAL l1em_2_med ******* **** MU_Central(LOW P_T) **** LOGICAL L1MU_CENTRAL,L2MU_CENTRAL,l2mu_central_stt **** **** MU_All(HI P_T) **** LOGICAL L1MU_ALL, L2MU_ALL **** **** DIMU_Central **** LOGICAL L1DIMU_CENTRAL, L2DIMU_CENTRAL,L2DIMU_CENTRAL_STT **** **** DIMU_ALL *** LOGICAL L1DIMU_ALL, L2DIMU_ALL **** **** E_Central *** LOGICAL L1E_CENTRAL,L2E_CENTRAL,L2E_CENTRAL_STT REAL THL2CELE, thl1cele INTEGER NL1CELE,NL2CELE,NL1CFTE,NL2CELE_STT DATA THL1CELE/2./ DATA THL2CELE/2./ **** **** E_ALL *** LOGICAL L1E_FORWARD,L2E_FORWARD REAL THL2FELE INTEGER NL1FELE,NL2FELE DATA THL2FELE/7./ ***dimu_el(wino zino) logical l1dimu_el, l2dimu_el, l2dimu_el_stt integer nl1dimu_el, nl2dimu_el, nl2dimu_el_stt real thl2dimu_el data thl2dimu_el /7./ ****diel_mu(wino zino) logical l1diel_mu, l2diel_mu, l2diel_mu_stt integer nl1diel_mu, nl2diel_mu real thl2diel_mu data thl2diel_mu /5./ **** DIELECTRON **** LOGICAL l1diel,l2diel REAL thl2diel INTEGER nl2diel DATA thl2diel/5./ ***** CENTRAL DIELECTRON logical l1diel_central, l2diel_central, l2diel_central_stt integer nl1diel_central, nl2diel_central, nl2diel_central_stt real thl2diel_central, thl1diel_central data thl2diel_central/2.0/ data thl1diel_central/2.0/ ***** FORWARD DIELECTRON logical l1diel_forward, l2diel_forward, l2diel_forward_stt integer nl1diel_forward, nl2diel_forward, nl2diel_forward_stt real thl2diel_forward data thl2diel_forward/7./ ***** ***DIELECTRON(ONE FORWARD ONE CENTRAL) logical l1diel_CE_FO, l2diel_CE_FO, l2diel_stt_CE_FO integer nl1diel_ce, nl1diel_fo, nl2diel_fo, nl2diel_ce, & nl2diel_ce_stt, nl1diel_CE_FO, nl2diel_CE_FO real thl2diel_CE_FO, thl2diel_stt_CE_FO data thl2diel_CE_FO, thl2diel_stt_CE_FO /3., 3./ **** END MENU VARIABLES **** C C Initialize C IF(first) THEN first=.false. * * Rescale thresholds * EMSCALE= 1.0 HSCALE= 0.6 misscale= 1.0 DO i=1,7 IF(i.LE.5) thl1em(i) = EMscale*thl1em(i) thl1jt(i)=hscale*thl1jt(i) thl1lt(i)=hscale*thl1lt(i) ENDDO thl2emjet_elec=emscale*thl2emjet_elec thl2emjet_jet=hscale*thl2emjet_jet thl2emmu_elec=emscale*thl2emmu_elec thl2miss_met =hscale*thl2miss_met thl1miss_ET =hscale*thl1miss_ET thl2miss_jet =hscale*thl2miss_jet thl1jet_3_miss_et =hscale*thl1jet_3_miss_et thl2jet_3_miss =hscale*thl2jet_3_miss thl2met_3_miss =hscale*thl2met_3_miss thl2mu_jet_high =hscale*thl2mu_jet_high thl2jt_multi =hscale*thl2jt_multi thl2jet_max=hscale*thl2jet_max thl2cele=emscale*thl2cele thl2fele=emscale*thl2fele thl2diel=emscale*thl2diel THL2JET_2_MED = hSCALE*THL2JET_2_MED THL2JET_2_min_CFT3 = hSCALE*THL2JET_2_min_CFT3 THL2JET_3_MED = hSCALE*THL2JET_3_MED * * * ENDIF IF(nt_njseeds.GT.20) nt_njseeds=20 * * Count L1 jet towers * DO i=1,7 nl1jt(i) = 0 ENDDO DO i=1,nt_njseeds DO j=1,7 IF(nt_ejseeds(i).GT.thl1jt(j)) nl1jt(j)=nl1jt(j)+1 ENDDO ENDDO * * count L1 eM towers * DO i=1,5 nl1em(i) = 0 ENDDO DO i=1,nt_neseeds DO j=1,5 IF(nt_eeseeds(i).GT.thl1em(j)) nl1em(j)=nl1em(j)+1 ENDDO ENDDO * C -------calculate fiber pt # tracks above l1 threshold------------------ * DO i=1,4 nl1cft(i)=0 nl1tagcft(i) = 0 ENDDO DO i=1,ntrig tagcft(i)=.false. rft16=itfcft16(i) rft16=rft16/3520. rft1=itfcft1(i) rft1=rft1/1280. rtemp=rft16-rft1 IF(rtemp.EQ.0) THEN rtemp=100. ELSE rtemp = 0.0143239*abs((1.0/rtemp)) ENDIF trakpt(i)=rtemp DO j=1,4 IF(rtemp.GT.thl1cft(j)) nl1cft(j)=nl1cft(j)+1 ENDDO cftphi=itfcft16(i)*.00178499 DO k=1,ntcps tagcftp=.false. tagcftn=.false. IF( (resppos(k).NE.999.AND.resppos(k).GT.thl1cps) .OR. & (respneg(k).NE.999.AND.respneg(k).GT.thl1cps) ) THEN delphipos=abs(cftphi-phipos(k)) delphineg=abs(cftphi-phineg(k)) IF( delphipos.LT.0.02) THEN IF(resppos(k).NE.999.AND.resppos(k).GT.thl1cps) THEN tagcftp=.true. END IF END IF IF( delphineg.LT.0.02) THEN IF(respneg(k).NE.999.AND.respneg(k).GT.thl1cps) THEN tagcftn=.true. END IF END IF END IF IF(tagcftp.OR.tagcftn) THEN tagcft(i)=.true. DO j=1,4 IF(rtemp.GT.thl1cft(j))nl1tagcft(j)=nl1tagcft(j)+1 ENDDO ENDIF ENDDO ENDDO * * TALLY STT TRACKS * * DO I=1,4 * NTRKSIPT(I) = 0 * ENDDO * DO I=1,SVX_TRACKNUM * DO J=1,4 * IF(ABS(SVX_SIPT(I)) .GT. SIPT(J)) NTRKSIPT(J)=NTRKSIPT(J)+1 * ENDDO * ENDDO *** ** * TALLY STT TRACKS AND CALCULATE PT OF THE TRACKS do i=1,4 ntrksipt(i)= 0 enddo do i=1,svx_tracknum DO J=1,4 if (svx_chsqr_new(i) .le. 10) then svx_sipt_new(J) = 0.003/svx_tk_new(i) else svx_sipt_new(J) = svx_cftpt(i) endif ENDDO ENDDO DO J=1,4 DO K=1,4 IF(ABS(SVX_SIPT_new(J)) .GT. SIPT(K)) NTRKSIPT(K) & =NTRKSIPT(K)+1 ENDDO enddo * * Tally Muons * DO I=1,4 NL1MUALL(I) = 0 ENDDO NL1MUALL(1)=MTM_XM1 NL1MUALL(2)=MTM_XM2 NL1MUALL(3)=MTM_XM3 NL1MUALL(4)=MTM_XM4 DO I=1,4 NL1MUCENTRAL(I)=0 ENDDO NL1MUCENTRAL(1)=MTM_CM1 NL1MUCENTRAL(2)=MTM_CM2 NL1MUCENTRAL(3)=MTM_CM3 NL1MUCENTRAL(4)=MTM_CM4 * * Test L2 electron candidate * DO k=1,nt_nelec l2elec_can(k) = .false. IF(nt_emfr(k).GT.0.90) THEN IF(nt_isofr(k).GT.0.80) THEN l2elec_can(k) = .true. ENDIF ENDIF ENDDO * * Construct large tiles (10 in eta & 4 in phi) * DO j=1,10 DO k=1,4 en(j,k) = 0. ENDDO ENDDO DO i=1,nt_njseeds j=(nt_rjseeds(i)-1)/4 + 1 k=(nt_ajseeds(i)-1)/10 + 1 en(j,k)=en(j,k)+nt_ejseeds(i) END DO DO i=1,7 nl1lt(i)=0 DO j=1,10 DO k=1,4 IF( en(j,k).gt.thl1lt(i) )nl1lt(i)=nl1lt(i)+1 ENDDO ENDDO ENDDO * * Tally CC Tiles * DO j=1,10 DO k=1,4 IF (j.LE.7 .AND. j.GT.3) THEN DO i=1,7 IF( en(j,k).gt.thl1lt(i) )nl1lj(i)=nl1lj(i)+1 ENDDO ENDIF ENDDO ENDDO * * Build Triggers * C--------------------------------------------------------------------- l1em_1_high = .false. IF (nl1em(4) .GE. 1 .AND. NL1JT(1) .GE.1.0) l1em_1_high=.true. l1em_2_med = .false. IF (nl1em(2) .GE. 2 .and. nl1em(3).ge.1) l1em_2_med=.true. C________________________________________________________________________ C Ask for EM(1,12,2.6) and JT(2,5,2.0) ! EM_JET DO I =1,5 NL1EM_ET(I) =0 ENDDO DO I=1,NT_NESEEDS IF (NT_RESEEDS(I).LE.21) IETA=NT_RESEEDS(I)-22 IF (NT_RESEEDS(I).GE.22) IETA=NT_RESEEDS(I)-21 IF (ABS(IETA) .LT.13) THEN DO J=1,5 IF (NT_EESEEDS(I) .GT. THL1EM(J)) NL1EM_ET(J) = NL1EM_ET(J) & +1 ENDDO ENDIF ENDDO DO I=1,7 NL1JT_ET(I) = 0 ENDDO DO I=1,NT_NJSEEDS IF (NT_RJSEEDS(I).LE.21) JETA=NT_RJSEEDS(I)-22 IF (NT_RJSEEDS(I).GE.22) JETA=NT_RJSEEDS(I)-21 IF (ABS(JETA) .LT. 10.) THEN DO J=1,7 IF (NT_EJSEEDS(I).GT.THL1JT(J)) NL1JT_ET(J) = NL1JT_ET(J) + & 1 ENDDO ENDIF ENDDO C----------------------------------------------------------- C EM_JET C test trig. C----------------------------------------------------------- l1em_jet = .FALSE. L2EM_JET = .FALSE. IF (NL1EM_ET(4) .GE. 1 .AND. NL1JT_ET(2) .GE. 2) THEN l1em_jet =.true. ENDIF C Ask for L2EM(1,15,0.95,0.80) and L2JT(2,10) nl2emjet_elec=0 l2em_jet = .false. IF (l1em_jet) THEN DO k=1,nt_nelec IF(nt_elet(k).GT.thl2emjet_elec .AND. l2elec_can(k)) THEN nl2emjet_elec=nl2emjet_elec+1 END IF END DO DO k=1,nt_njet IF(nt_jtet(k).GT.thl2emjet_jet) nl2emjet_jet=nl2emjet_jet+1 END DO IF(nl2emjet_elec.GE.1 .AND. nl2emjet_jet.GE.2) l2em_jet=.true. ENDIF C________________________________________________________________________ C ASK MU(1,3) AND EM(1,7) !EM_MU nl2emmu_elec = 0 l1em_mu = .false. IF (nl1em(3).GE.1 .AND. nl1muALL(2) .GE.1) l1em_mu = .true. C Ask for L2EM(1,7,0.95,0.80) l2em_mu = .false. IF (l1em_mu) THEN DO k=1,nt_nelec IF(nt_elet(k).GT.thl2emmu_elec .AND. l2elec_can(k) ) THEN nl2emmu_elec = nl2emmu_elec+1 ENDIF ENDDO IF (nl2emmu_elec.GE.1)l2em_mu = .true. ENDIF C________________________________________________________________________ C Ask JT(1,5) and MET(40) !MISSING_ET nl2miss_et = 0. l1missing_et = .false. IF (nl1jt(2).GE.1 .AND. nt_met.GE.thl1miss_Et) & l1missing_et = .true. l2missing_et = .false. C Ask L2JT(1,10) and MET(40) IF (l1missing_et) THEN DO k=1,nt_njet IF(nt_jtet(k).GT.thl2miss_jet) nl2miss_et=nl2miss_et+1 ENDDO IF (nl2miss_et.GE.1 .AND. nt_met.GE.thl2miss_met)l2missing_et = & .true. ENDIF C________________________________________________________________________ C Ask JT(3,5) and MET(20) !JET_3_MISS nl2jet_3_miss = 0. l1jet_3_miss = .false. l2jet_3_miss = .false. IF ((nl1jt(3).GE.3) & .AND. nt_met.GE.thl1jet_3_miss_et) & l1jet_3_miss = .true. C Ask L2JT(3,20) and MET(25) IF (l1jet_3_miss) THEN DO k=1,nt_njet IF(nt_jtet(k).GT.thl2jet_3_miss) & nl2jet_3_miss=nl2jet_3_miss+1 ENDDO IF (nl2jet_3_miss.GE.3 .AND. nt_met.GE.thl2met_3_miss) & l2jet_3_miss = .true. ENDIF C________________________________________________________________________ C Ask MU(1,3) !MU_MAX l1mu_1_high = .false. l2mu_1_high = .false. IF (nl1muALL(3) .GE. 1) THEN l1mu_1_high = .true. ENDIF if(l1mu_1_high) then l2mu_1_high=.true. end if **** CENTRAL MUON(LOW P_T) ************************************* *** MU(1,3) **************************************************** L1MU_CENTRAL = .FALSE. L2MU_CENTRAL = .FALSE. L2MU_CENTRAL_STT = .FALSE. IF (NL1MUCENTRAL(2) .GE. 1) L1MU_CENTRAL = .TRUE. IF(L1MU_CENTRAL) THEN L2MU_CENTRAL=.TRUE. IF (NTRKSIPT(2) .GE.1) THEN L2MU_CENTRAL_STT=.TRUE. END IF END IF C________________________________________________________________________ C Ask MU(2,3) !MU_2_MED l1mu_2_med = .false. l2mu_2_med = .false. IF (nl1muALL(2) .GE. 2.) THEN l1mu_2_med = .true. l2mu_2_med = .true. ENDIF C________________________________________________________________________ C Ask MU(1,1) and JT(2,7) !MU_JET nl1mu_jet_high = 0. nl2mu_jet_high= 0. l1mu_jet_high = .false. l2mu_jet_high = .false. IF (nl1muALL(1).GE.1. .AND. & (nl1jt(3).GE. 2. .and. nl1jt(4).ge.1) ) THEN l1mu_jet_high = .true. ENDIF C--------L2MU_JET_HIGH C Ask MU(1,1) and L2JT(1,10) nl2mu_jet_high=0 IF(l1mu_jet_high) THEN DO k=1,nt_njet IF(nt_jtet(k).GT.thl2mu_jet_high) & nl2mu_jet_high=nl2mu_jet_high+1 ENDDO IF(nl1muALL(2).GE.1 .AND. NL2mu_jet_High.GE.1) l2mu_jet_high= & .true. END IF C_________________________________________________________________________ C Ask JT(3,7) and LJ(3,15) !JET_MULTI nl2jt_multi = 0 l1jet_multi = .false. l2jet_multi = .false. IF ((nl1lt(2) .GE. 3 .and. nl1lt(3).ge.2) .AND. & (nl1jt(2) .GE. 3 .and. nl1jt(3).ge.2)) THEN l1jet_multi = .true. ENDIF C Ask L2JT(5,10) IF (l1jet_multi) THEN DO k=1,nt_njet IF (nt_jtet(k) .GT. thl2jt_multi ) nl2jt_multi= & nl2jt_multi+1 ENDDO IF (nl2jt_multi .GE. 5) l2jet_multi = .true. ENDIF C_________________________________________________________________________ C Ask LJ(1,45) !JET_MAX nl2jmax = 0 l1jet_max = .false. l2jet_max = .false. IF (nl1lt(7) .GE. 1) l1jet_max= .true. C Ask for L2JT(85) IF (l1jet_max) THEN DO i=1,nt_njet IF(nt_jtet(i) .GT. thl2jet_max) nl2jmax=nl2jmax+1 ENDDO IF(nl2jmax.GE.1) l2jet_max=.true. ENDIF c_________________________________________________________________________ C Ask JT(2,7) !JET_2_MED nl2jET_2_mED = 0 l1jet_2_MED = .false. l2jet_2_MED = .false. IF (nl1jt(4).ge.1 .and. nl1jt(3) .GE. 2 ) THEN l1jet_2_MED = .true. ENDIF C Ask L2JT(2,20) IF (l1jet_2_MED) THEN DO k=1,nt_njet IF (nt_jtet(k) .GT. thl2jET_2_MED ) nl2jET_2_MED= & nl2jET_2_MED+1 ENDDO IF (nl2jET_2_MED .GE. 2) l2jet_2_MED = .true. ENDIF c_________________________________________________________________________ C Ask JT(2,5GeV) and cft(3,3GEv) !JET_2_min_cft3 nl2jET_2_min_cft3 = 0 l1jet_2_Min_cft3 = .false. l2jet_2_Min_cft3 = .false. IF ((nl1jt(2).ge.2) & .AND. nl1cft(2).GE.3 ) THEN l1jet_2_min_cft3 = .true. ENDIF C Ask L2JT(2,10) IF (l1jet_2_min_cft3) THEN DO k=1,nt_njet IF (nt_jtet(k) .GT. thl2jet_2_min_cft3 ) nl2jET_2_Min_cft3= & nl2jET_2_min_cft3+1 ENDDO IF (nl2jET_2_min_cft3 .GE. 2) l2jet_2_min_cft3 = .true. ENDIF C c_________________________________________________________________________ C Ask JT(3,5GeV) and JT(1,7) !JET_3_MED nl2jET_3_med = 0 l1jet_3_med = .false. l2jet_3_med = .false. IF (nl1jt(3).ge. 1 .and. nl1jt(2).ge.3) THEN l1jet_3_med = .true. ENDIF C Ask L2JT(3,10) IF (l1jet_3_med) THEN DO k=1,nt_njet IF (nt_jtet(k) .GT. thl2jet_3_med ) nl2jET_3_med= & nl2jET_3_med+1 ENDDO IF (nl2jET_3_med .GE. 3) l2jet_3_med = .true. ENDIF c_________________________________________________________________________ * DI MUON(central WITHOUT and WITH STT)!LOW PT c MU(2,3) C FOR L2 STT TRACK > 3 GEV L1DIMU_CENTRAL = .FALSE. L2DIMU_CENTRAL = .FALSE. L2DIMU_CENTRAL_STT = .FALSE. IF(NL1MUCENTRAL(2).GE.2) & L1DIMU_CENTRAL=.TRUE. IF(L1DIMU_CENTRAL) THEN L2DIMU_CENTRAL=.TRUE. IF(NTRKSIPT(1).GE.2) THEN L2DIMU_CENTRAL_STT=.TRUE. END IF END IF * * DI MUON (all) HIGH P_T * L1DIMU_ALL = .FALSE. L2DIMU_ALL = .FALSE. IF(NL1MUALL(3).GE.2) L1DIMU_ALL=.TRUE. IF(L1DIMU_ALL) THEN L2DIMU_ALL=.TRUE. END IF c-------------------------------------------------------------- * Electron(central WITHOUT and WITH STT) * EM(1,3,1.2) and CTP(1,3) * L2EM(1,5,0.95,0.80,eta=1.2,CPTmatch) * L1E_CENTRAL = .FALSE. L2E_CENTRAL = .FALSE. L2E_CENTRAL_STT = .FALSE. NL1CELE = 0 NL2CELE = 0 NL2CELE_STT=0 DO I=1,NT_NESEEDS IF(NT_RESEEDS(I).LE.21) IETA=NT_RESEEDS(I)-22 IF(NT_RESEEDS(I).GE.22) IETA=NT_RESEEDS(I)-21 IF(ABS(IETA).LE.6 .AND. NT_EESEEDS(I).GT.THL1cele) THEN NL1CELE=NL1CELE+1 END IF END DO IF(NL1CELE.GE.1 .AND. NL1TAGCFT(2).GE.1) L1E_CENTRAL=.TRUE. IF(L1E_CENTRAL) THEN DO K=1,NT_NELEC IF(ABS(NT_ELETA(K)).LT.1.2) THEN IF(NT_ELET(K).GT.THL2CELE .AND. L2ELEC_CAN(K)) THEN DO I=1,NTRIG CFTPHI=ITFCFT16(I)*.00178499 RTEMP= ABS(CFTPHI-(NT_ELPHI(K)-0.06)) IF( RTEMP.LT.0.15 .AND. & (TAGCFT(I).AND.TRAKPT(I).GT.3.0) ) THEN NL2CELE=NL2CELE+1 END IF END DO DO J=1,SVX_TRACKNUM RTEMP= ABS(SVX_PHI0C(J)-NT_ELPHI(K)) IF( RTEMP.LT.0.15.AND.SVX_SIPT(J).GT.3.0) THEN NL2CELE_STT=NL2CELE_STT+1 END IF END DO END IF END IF END DO IF(NL2CELE.GE.1) THEN L2E_CENTRAL=.TRUE. END IF IF(NL2CELE_STT.GE.1) THEN L2E_CENTRAL_STT=.TRUE. END IF END IF c--------------------------------------------------------------- * * Electron(Forward) * EM(1,7,gt1.2) * L2EM(1,7,0.95,0.80,etagt1.2) * L1E_FORWARD = .FALSE. L2E_FORWARD = .FALSE. NL1FELE = 0 NL2FELE = 0 DO I=1,NT_NESEEDS IF(NT_RESEEDS(I).LE.21) IETA=NT_RESEEDS(I)-22 IF(NT_RESEEDS(I).GE.22) IETA=NT_RESEEDS(I)-21 IF(ABS(IETA).GT.6 .AND. NT_EESEEDS(I).GT.THL1EM(3)) THEN NL1FELE=NL1FELE+1 END IF END DO IF(NL1FELE.GE.1) L1E_FORWARD=.TRUE. IF(L1E_FORWARD) THEN DO K=1,NT_NELEC IF(ABS(NT_ELETA(K)).GT.1.2) THEN IF(NT_ELET(K).GT.THL2FELE .AND. L2ELEC_CAN(K)) THEN NL2FELE=NL2FELE+1 END IF END IF END DO IF(NL2FELE.GE.1) L2E_FORWARD=.TRUE. END IF c----------------------------------------------------------- * ****WINI ZINO ( ONE CENTRAL MUON AND TWO ELECTRONS) NL2diel_mu = 0 L1diel_mu =.FALSE. L2diel_mu =.FALSE. l2diel_mu_stt=.false. IF (NL1MUCENTRAL(2).GE.1 .AND. NL1em(2).GE.2) L1diel_mu=.TRUE. IF(L1diel_mu) THEN DO K=1,NT_NELEC IF(NT_ELET(K).GT. THL2diel_mu .AND. L2ELEC_CAN(K)) THEN NL2diel_mu = NL2diel_mu+1 END IF END DO IF(NL2diel_mu.GE.2) L2diel_mu=.TRUE. IF(NL2diel_mu.GE.2 .and.ntrksipt(2).ge.1) L2diel_mu_stt=.TRUE. END IF c----------------------------------------------------------------- * wino_zino (dimuon+one electron)with and without stt * nl1dimu_el=0 NL2dimu_el = 0 NL2dimu_el_stt = 0 L1dimu_el = .FALSE. L2dimu_el =.FALSE. L2dimu_el_stt =.FALSE. DO I=1,NT_NESEEDS IF(NT_RESEEDS(I).LE.21) IETA=NT_RESEEDS(I)-22 IF(NT_RESEEDS(I).GE.22) IETA=NT_RESEEDS(I)-21 IF(ABS(IETA).LE.6 .AND. NT_EESEEDS(I).GT.THL1EM(1)) THEN nl1dimu_el=nl1dimu_el+1 endif enddo if (nl1dimu_el .ge. 1 .and. (NL1TAGCFT(2).GE.1 .and. & nl1muall(1) .ge. 2)) then L1dimu_el=.TRUE. endif ** IF(L1dimu_el) THEN DO K=1,NT_NELEC IF(ABS(NT_ELETA(K)).LT.1.2) THEN IF (NT_ELET(K).GT.THL2DIMU_EL .AND. L2ELEC_CAN(K)) THEN DO I=1,NTRIG CFTPHI=ITFCFT16(I)*.00178499 RTEMP= ABS(CFTPHI-(NT_ELPHI(K)-0.06)) IF( RTEMP.LT.0.15 .AND. (TAGCFT(I) .and. trakpt(i) & .gt. 3 )) THEN NL2dimu_el=NL2dimu_el+1 endif enddo do j=1,svx_tracknum rtemp= abs(svx_phi0c(j)-nt_elphi(k)) if (rtemp .lt. 0.15 .and. svx_sipt(j) .gt. 3.0) then nl2dimu_el_stt=nl2dimu_el_stt+1 endif enddo endif endif enddo IF(NL2DIMU_EL.GE.1) L2DIMU_El=.TRUE. IF(NL2DIMU_EL_STT.GE.1) L2DIMU_EL_STT=.TRUE. END IF c----------------------------------------------------------------- ********** DIELECTRON C Ask EM(2,7) !DI-ELECtron nl2diel = 0 l1diel = .false. l2diel = .false. IF(nl1em(2).GE.2 .and. nl1em(3).ge.1) l1diel=.true. C Ask for L2EM(2,5,0.95,0.80) IF(l1diel) THEN DO k=1,nt_nelec IF(nt_elet(k).GT.thl2diel .AND. l2elec_can(k)) THEN nl2diel=nl2diel+1 END IF END DO IF(nl2diel.GE.2) l2diel=.true. ENDIF c--------------------------------------------------------------- **central dieletron nl1diel_central=0 nl2diel_central=0 nl2diel_central_stt=0 l1diel_central =.false. l2diel_central =.false. l2diel_central_stt =.false. do i=1,nt_neseeds IF(NT_RESEEDS(I).LE.21) IETA=NT_RESEEDS(I)-22 IF(NT_RESEEDS(I).GE.22) IETA=NT_RESEEDS(I)-21 if (abs(ieta) .le. 6 .and. nt_eeseeds(i) .ge. thl1diel_central) & then nl1diel_central=nl1diel_central+1 endif enddo if (nl1diel_central .ge.2 .and. nl1tagcft(2).ge.2) & l1diel_central=.true. if(l1diel_central) then do k=1,nt_nelec if (abs(nt_eleta(k)) .lt.1.2) then if (nt_elet(k) .gt. thl2diel_central .and. l2elec_can(k)) & then do i=1,ntrig cftphi=itfcft16(i)*.00178499 rtemp=abs(cftphi-(nt_elphi(k))-0.06) IF (RTEMP .lt.0.15 .and. (tagcft(i) .and. trakpt(i).gt. & 3.0 )) then nl2diel_central=nl2diel_central+1 endif enddo do j=1,svx_tracknum rtemp=abs(svx_phi0c(j)-nt_elphi(k)) if (rtemp .lt. 0.15 .and. svx_sipt(j) .gt. 3.0) then nl2diel_central_stt=nl2diel_central_stt+1 endif enddo endif endif enddo if (nl2diel_central .ge. 2) l2diel_central=.true. if (nl2diel_central_stt .ge. 2) l2diel_central_stt=.true. endif c------------------------------------------------------------ ***forward dielectron C L1EM(2,7,GT1.2) C L2EM(2,7,0.95,0.80,ETA>1.2) nl1diel_forward=0 nl2diel_forward =0 l1diel_forward =.false. l2diel_forward =.false. l2diel_forward_stt =.false. do i=1,nt_neseeds IF(NT_RESEEDS(I).LE.21) IETA=NT_RESEEDS(I)-22 IF(NT_RESEEDS(I).GE.22) IETA=NT_RESEEDS(I)-21 if (abs(ieta) .gt. 6 .and. nt_eeseeds(i) .ge. thl1em(3)) & then nl1diel_forward=nl1diel_forward+1 endif enddo if (nl1diel_forward .ge. 2) then l1diel_forward=.true. endif if (l1diel_forward) then do k=1,nt_nelec if (abs(nt_eleta(k)) .gt. 1.2) then if (nt_elet(k) .gt. thl2diel_forward .and. l2elec_can(k)) & then nl2diel_forward=nl2diel_forward+1 endif endif enddo if (nl1diel_forward .ge. 2) l2diel_forward=.true. endif c-------------------------------------------------------------- **dielectron( one central and one forward) with and without stt nl1diel_ce=0 nl1diel_fo=0 nl2diel_CE_FO=0 nl2diel_ce=0 nl2diel_fo=0 nl2diel_ce_stt=0 l1diel_CE_FO=.false. l2diel_CE_FO=.false. l2diel_stt_CE_FO=.false. do i=1,nt_neseeds IF(NT_RESEEDS(I).LE.21) IETA=NT_RESEEDS(I)-22 IF(NT_RESEEDS(I).GE.22) IETA=NT_RESEEDS(I)-21 if (abs(ieta) .le. 6 .and. nt_eeseeds(i) .gt. thl1em(1)) then nl1diel_ce=nl1diel_ce+1 endif if (abs(ieta) .ge. 6 .and. nt_eeseeds(i) .gt. thl1em(1))then nl1diel_fo=nl1diel_fo+1 endif enddo if (nl1diel_ce .ge. 1 .and. nl1diel_fo .ge. 1.and. nl1tagcft(2) & .ge. 1.) then l1diel_CE_FO=.true. endif if (l1diel_CE_FO) then do k=1,nt_nelec if (abs(nt_eleta(k)) .lt. 1.2) then if (nt_elet(k) .gt. thl2diel_CE_FO .and. l2elec_can(k)) & then do i=1,ntrig cftphi=itfcft16(i)*.00178499 rtemp=abs(cftphi-(nt_elphi(k)-0.06)) if (rtemp .lt. 0.15 .and. (tagcft(i) .and. trakpt(i) & .gt.3.0)) then nl2diel_ce=nl2diel_ce+1 endif enddo do j=1,svx_tracknum rtemp=abs(svx_phi0c(j)-nt_elphi(k)) if (rtemp .lt. 0.15 .and. svx_sipt(j) .gt. 3.0) then nl2diel_ce_stt=nl2diel_ce_stt+1 endif enddo endif endif if (abs(nt_eleta(k)) .gt. 1.2) then if (nt_elet(k).gt.thl2diel_CE_FO .and. l2elec_can(k)) & then nl2diel_fo=nl2diel_fo+1 endif endif enddo if (nl2diel_ce .ge. 1 .and. nl2diel_fo .ge.1) l2diel_CE_FO=.true. if (nl2diel_ce_stt .ge. 1 .and.nl2diel_fo .ge.1) & l2diel_stt_CE_FO=.true. endif C C--- collect terms------------------- C vtrig(1) = vtrig(1)+1 IF (l1em_1_high) vtrig(2) = vtrig(2) +1 IF (l1em_2_med) vtrig(3) = vtrig(3) +1 IF (l1em_jet) vtrig(4) = vtrig(4) +1 IF (l2em_jet ) vtrig(5) = vtrig(5) +1 IF (l1em_mu) vtrig(6) = vtrig(6) +1 IF (l2em_mu) vtrig(7) = vtrig(7) +1 IF (l1missing_et) vtrig(8) = vtrig(8) +1 IF (l2missing_et) vtrig(9) = vtrig(9) +1 IF (l1jet_3_miss) vtrig(10) = vtrig(10) +1 IF (l2jet_3_miss) vtrig(11) = vtrig(11) +1 IF (l1mu_1_HIGH) vtrig(12) = vtrig(12) +1 IF (l2mu_1_high) vtrig(13) = vtrig(13)+1 IF (L1MU_CENTRAL) VTRIG(14) = VTRIG(14) +1 IF (L2MU_CENTRAL) VTRIG(15) = VTRIG(15) +1 IF (L2MU_CENTRAL_STT) VTRIG(16) = VTRIG(16) +1 IF (l1mu_2_med) vtrig(17) = vtrig(17) +1 IF (l2mu_2_med) vtrig(18) = vtrig(18) +1 IF (l1mu_jet_high) vtrig(19) = vtrig(19) +1 IF (l2mu_jet_high) vtrig(20) = vtrig(20) +1 IF (l1jet_multi) vtrig(21) = vtrig(21) +1 IF (l2jet_multi) vtrig(22) = vtrig(22) +1 IF (l1jet_max) vtrig(23) = vtrig(23) +1 IF (l2jet_max) vtrig(24) = vtrig(24) +1 IF (L1JET_2_MED ) vtrig(25) = vtrig(25) +1 IF (L2JET_2_MED ) vtrig(26) = vtrig(26) +1 IF (L1JET_2_Min_cft3 ) vtrig(27) = vtrig(27) +1 IF (L2JET_2_Min_cft3 ) vtrig(28) = vtrig(28) +1 IF (L1JET_3_Med ) vtrig(29) = vtrig(29) +1 IF (L2JET_3_Med ) vtrig(30) = vtrig(30) +1 IF (L1DIMU_CENTRAL) VTRIG(31) = VTRIG(31) +1 IF (L2DIMU_CENTRAL) VTRIG(32) = VTRIG(32) +1 IF (L2DIMU_CENTRAL_STT) VTRIG(33) = VTRIG(33)+1 IF (L1DIMU_ALL) VTRIG(34) = VTRIG(34) +1 IF (L2DIMU_ALL) VTRIG(35) = VTRIG(35) +1 IF (L1E_CENTRAL) VTRIG(36) = VTRIG(36) +1 IF (L2E_CENTRAL) VTRIG(37) = VTRIG(37) +1 IF (L2E_CENTRAL_STT) VTRIG(38) = VTRIG(38) +1 IF (L1E_FORWARD) VTRIG(39) = VTRIG(39) +1 IF (L2E_FORWARD) VTRIG(40) = VTRIG(40) +1 IF(L1E_CENTRAL.OR.L1E_FORWARD) VTRIG(41) = VTRIG(41)+1 IF(L2E_CENTRAL.OR.L2E_FORWARD) VTRIG(42) = VTRIG(42)+1 IF(L2E_CENTRAL_STT.OR.L2E_FORWARD) VTRIG(43) = VTRIG(43)+1 if (l1diel_mu) vtrig(44) = vtrig(44) +1 if (l2diel_mu) vtrig(45) = vtrig(45)+1 if (l2diel_mu_stt) vtrig(46) = vtrig(46) +1 if (l1dimu_el) vtrig(47) = vtrig(47) +1 if (l2dimu_el) vtrig(48) = vtrig(48) +1 if (l2dimu_el_stt) vtrig(49) = vtrig(49) +1 IF (l1diel ) vtrig(50) = vtrig(50) +1 IF (l2diel ) vtrig(51) = vtrig(51) +1 if (l1diel_central) vtrig(52) = vtrig(52) +1 if (l2diel_central) vtrig(53) = vtrig(53) +1 if (l2diel_central_stt) vtrig(54) = vtrig(54) +1 if (l1diel_forward) vtrig(55) = vtrig(55) +1 if (l2diel_forward) vtrig(56) = vtrig(56) +1 if (l1diel_CE_FO) vtrig(57) = vtrig(57) +1 if (l2diel_CE_FO) vtrig(58) = vtrig(58) +1 if (l2diel_stt_CE_FO) vtrig(59) = vtrig(59) +1 if (l1diel_central .or. l1diel_forward .or. l1diel_CE_FO) & vtrig(60)=vtrig(60)+1 if (l2diel_central .or. l2diel_forward .or. l2diel_CE_FO) & vtrig(61)=vtrig(61)+1 if (l2diel_central_stt .or. l2diel_forward .or. l2diel_stt_CE_FO) & vtrig(62)=vtrig(62)+1 c-------------------------------------------------------- c highpt only(without stt) c_________________________________________________________ IF ( l1em_jet .OR. l1em_mu .OR. l1missing_et & .OR. l1jet_3_miss .OR. l1mu_1_high .OR. l1mu_2_med & .OR. l1mu_jet_high .OR. l1jet_multi .OR. l1jet_max & .OR. l1diel .OR. L1DIMU_ALL .or. L1DIMU_CENTRAL .or. & L1E_CENTRAL .or. L1E_FORWARD .or. l1diel_mu & .or. l1dimu_el .or. l1diel_central .or. l1diel_forward .or. & l1diel_CE_FO ) & vtrig(73) = vtrig(73) +1 IF ( l2em_jet .OR. l2em_mu .OR. l2missing_et & .OR. l2jet_3_miss .OR. l2mu_1_high .OR. l2mu_2_med & .OR. l2mu_jet_high .OR. l2jet_multi .OR. l2jet_max & .OR. l2diel .OR. L2DIMU_CENTRAL .or. L2DIMU_ALL .or. & L2E_CENTRAL .or. L2E_FORWARD .or. l2diel_mu .or. l2diel_CE_FO) & vtrig(74) = vtrig(74) +1 c____________________________________________________ c highpt only (with stt) c____________________________________________________ if (l2em_jet .or.l2em_mu .or. l2missing_et .or. l2jet_3_miss .or. & l2mu_1_high .or. l2mu_2_med .or. l2mu_jet_high .or. l2jet_multi & .or. l2jet_max .or. l2diel .or. L2DIMU_ALL .or. & L2DIMU_CENTRAL_STT .or. L2E_CENTRAL_STT .or. l2diel_mu_stt & .or. l2dimu_el_stt .or. l2diel_central_stt .or. & l2diel_stt_CE_FO) & vtrig(75) =vtrig(75)+1 c_____________________________________________________ c highpt only (intendent to use stt) c_____________________________________________________ if (l1JET_2_MED .or. L1JET_2_Min_cft3 .or. L1JET_3_Med) vtrig(76) & =vtrig(76) +1 if (l2JET_2_MED .or. L2JET_2_Min_cft3 .or. L2JET_3_Med) vtrig(77) & =vtrig(77) +1 c----------------------------------------------------- C HIGH PT + SOFT STUFF c----------------------------------------------------- IF ( l1em_jet .OR. l1em_mu .OR. l1missing_et & .OR. l1jet_3_miss .OR. l1mu_1_high .OR. L1MU_CENTRAL .OR. & l1mu_2_med .OR. l1mu_jet_high .OR. l1jet_multi .OR. l1jet_max & .OR. L1JET_2_MED .OR. L1JET_2_Min_cft3 & .OR. L1JET_3_Med .OR. L1DIMU_CENTRAL .OR. L1DIMU_ALL .OR. & L1E_CENTRAL .OR. L1E_FORWARD .OR. l1diel_mu .OR. l1dimu_el .OR. & L1DIEL .OR. l1diel_centraL .OR. l1diel_forward .OR. & l1diel_CE_FO) & vtrig(78) = vtrig(78) +1 IF ( l2em_jet .OR. l2em_mu .OR. l2missing_et & .OR. l2jet_3_miss .OR. l2mu_1_high .OR.L1MU_CENTRAL .OR. & l2mu_2_med .OR. l2mu_jet_high .OR. l2jet_multi .OR. l2jet_max & .OR. L2JET_2_MED .OR. L2JET_2_Min_cft3 & .OR. L2JET_3_Med .OR. L2DIMU_CENTRAL & .OR. L2DIMU_ALL .OR. L2E_CENTRAL .OR. & L2E_FORWARD .OR. l2diel_mu .OR. l2dimu_eL & .OR. l2diel_central & .OR. l2diel_forward .OR. l2diel .OR. l2diel_CE_FO) & vtrig(79) = vtrig(79) +1 c-------------------------------------- c HIGH PT + SOFT STUFF(with STT) c-------------------------------------- IF ( l2em_jet .OR. l2em_mu .OR. l2missing_et & .OR. l2jet_3_miss .OR. l2mu_1_high .OR.L1MU_CENTRAL .OR. & l2mu_2_med .OR. l2mu_jet_high .OR. l2jet_multi .OR. l2jet_max & .OR. L2JET_2_MED .OR. L2JET_2_Min_cft3 & .OR. L2JET_3_Med .OR. L2DIMU_CENTRAL_stt & .OR. L2DIMU_ALL .OR. L2E_CENTRAL_stt .OR. & L2E_FORWARD .OR. l2diel_mu_stt .OR. l2dimu_eL_stt & .OR. l2diel_central_stt & .OR. l2diel_forward .OR. l2diel .OR. l2diel_stt_CE_FO) & vtrig(80) = vtrig(80) +1 full_trigmenu = 1. * END