REAL FUNCTION my(XDUMMY) ********************************************************* * * * This file was generated by HUWFUN. * * * ********************************************************* * * Ntuple Id: 1 * Ntuple Title: PNN ntup * Creation: 28/06/:2 17.08.05 * ********************************************************* * LOGICAL CHAIN CHARACTER*128 CFILE INTEGER IDNEVT,NCHEVT,ICHEVT REAL OBS(13) * COMMON /PAWIDN/ IDNEVT,OBS COMMON /PAWCHN/ CHAIN, NCHEVT, ICHEVT COMMON /PAWCHC/ CFILE * *-- Ntuple Variable Declarations * REAL scal_tsp,scal_io2,scal_io3,utbfld,rsdt1,rsdt2,trs,erd,tavgz + ,tntg,rsintime,cmupi1(2),tmuon1(2),emuon1(2),emuonc1(2),dt1,emut1 + ,emutc1,prod1,tmuav1,muped1(2),piped1(2),cmuend1(2),cpiend1(2) + ,tpion1(2),epion1(2),epionc1(2),ers,rrs,cos3d,pdc,zf,elast + ,elast_cor,tk,tpi,etg,rtg,tgx,tgy,tgz,ek,eic,ric,eicest,ptot + ,ptote,rtot,etot,b4strob,tdot2t,ictime,tds,rsz1d,rsz2d,rsz1d2 + ,rsz2d2,zutout,radius,d0,phi0,tanl,chidfxy,chidfz,rsxy1d,rsxy2d + ,rscrsd(10),rsr1,rsp1,rsz1,rsr2,rsp2,rsz2,rscrsr(10),rscrsp(10) + ,eicout,eicin,ticout,ticin,tdtic_i,tdtic_o,tdaic_i,tdaic_o,PEDS1 + ,PEDD1,PEDS2,PEDD2,TDFRAW1(128),TDFPS1(128),TDFP1D1(128) + ,TDFP2D1(128),TDFRAW2(128),TDFPS2(128),TDFP1D2(128),TDFP2D2(128) + ,PR_RF,PR_RFZ1,PR_RFZ2,PR_RFZ2B,LIKE,LIKE2,XIN,YIN,DB4,DB4TIP + ,DVXTIP,DVXPI,KPIANG,DEXTIP,LOGUTCQ,LOGUTCQI(6),rs_rate(10) + ,cos3d1,radius1,d01,phi01,tanl1,chidfxy1,tgx1,tgy1,tkold,tpiold + ,DRIFT_TIME(24),RESID(24),ADCW(24),ZRES(6),ZZPOS(6),CL_EPRCOR + ,DELTARP,DELTARE,SIG_P,SIG_E,FDEDX,CL_RSDEDX,CL_RSDEDX2 + ,CHIMAX_RSDEDX,DZUTSC14,L14ACC,TK_TG(100),EK_TG(100),TPI_TG(100) + ,EPI_TG(100),TPIOP_TG(100),EPIOP_TG(100),CCDTK(100),CCDEK(100) + ,CCDHK(100) INTEGER run,event,splno,splrec,icut,lay,tlay,sec,tsec,ssec,thex + ,t2sec,rssec1,rssec2,ipiflg,layv4,secv4,ndclay,ictpi,itgqualt + ,nhxy,nhz,ncross,iicout,iicin,NRAW1,NPS1,NP1D1,NP2D1,NRAW2,NPS2 + ,NP1D2,NP2D2,UTCHNU1,UTCHNU2,UTCHNU3,nhxy1,NK_TG1,ELK_TG1(100) + ,NPI_TG1,ELPI_TG1(100),NHITXY,NHITZ,NLAYER(24),NFOIL(6) + ,LAY_CHIMAX,NPRPAT,ICODE_L14,ICODE_L14ACC,NK_TG,ELK_TG(100) + ,NPI_TG,ELPI_TG(100),NPIOP_TG,ELPIOP_TG(100),NEDGE_TG + ,ELEDGE_TG(50),NTGVET_TG,NTRIK_TG,NKCLUST,NKCLUST2,IFAIL_KIC + ,NK_TG_S,N1,N2,N3,N4,N5,N6,N7,N8,N9,N10,N11,N12,N13,N14,N15 + ,CCD1(40),CCD2(40),CCD3(40),CCD4(40),CCD5(40),CCD6(40),CCD7(40) + ,CCD8(40),CCD9(40),CCD10(40),CCD11(40),CCD12(40),CCD13(40) + ,CCD14(40),CCD15(40),CCDPF(15),CCDBIN(15),NELVET,ISEV(50),ISTEV5 + ,ISTEV5M1,IFAIL_TDFL,NHEXCT,ICODE_TDFL(25),NHIT_STHEX,NNOFID + ,CHNOFID(288),IFID_CODEX(288),NPVRD,ELEMPVRD(100),PVRDTE1(100) + ,PVRDTR1(100),PVRDTH1(100),PVRDTF1(100),PVRDTE2(100),PVRDTR2(100) + ,PVRDTH2(100),PVRDTF2(100),NPVBV,ELEMPVBV(50),NPVBL,ELEMPVBL(50) + ,NPVEC,ELEMPVEC(80),NPVIC,ELEMPVIC(40),NPVCO,ELEMPVCO(20),NPVCM + ,ELEMPVCM(20),NPVTG,ELEMPVTG(100),NPVBS,ELEMPVBS(32),NSABTRD1 + ,NSABTRD2 REAL CCDRK(100),CCDTPI(100),CCDEPI(100),CCDHPI(100),CCDRPI(100) + ,CCDTPIOP(100),CCDEPIOP(100),CCDHPIOP(100),CCDRPIOP(100) + ,TEDGE_TG(50),EEDGE_TG(50),ETGVET_TG,DIST_TARGF,EPIMAX_TG + ,TPIMAX_TG,XTGT_TP,YTGT_TP,ZTGT_TP,DTGT_TP,RTDIF_TG,EPIMAXK + ,PHIVTX,DRPMAX,DRPMIN,PLN_SW,PLNV_SW,X_B4SW,Y_B4SW,DTK,DTPI,T1MT2 + ,EELEC1,EELEC2,EELEC2C,TELEC,DTELEC,PROBTD,PROBTDN,EMUT_TDFL(25) + ,EMUTC_TDFL(25),TMUAV_TDFL(25),DT_TDFL(25),DTMU_TDFL(25) + ,CMUEND_TDFL(2,25),CPIEND_TDFL(2,25),CMUPI_TDFL(2,25) + ,PROD_TDFL(25),ESUM_STHEX,EMAX_STHEX,ESTHEX(4),TIME_FIDX(288) + ,TPVRD(100),EPVRD(100),EPVRDPR,TZPVRD(100),EZPVRD(100),TPVBV(50) + ,EPVBV(50),EPVBVPR,TZPVBV(50),EZPVBV(50),TPVBL(50),EPVBL(50) + ,EPVBLPR,TZPVBL(50),EZPVBL(50),TPVEC(80),EPVEC(80),HPVEC(80) + ,RPVEC(80),EPVECPR,EPVEUPR,EPVECGT,EPVEUGT,EPVECFERA(80) + ,TPVIC(40),EPVIC(40),EPVICPR,EPVVCPR,TPVCO(20),EPVCO(20),EPVCOPR + ,TPVCM(20),EPVCM(20),EPVCMPR,TPVTG(100),EPVTG(100),EPVTGPR + ,CCDTPV(100),CCDEPV(100),CCDHPV(100),CCDRPV(100),TPVBS(32) + ,HPVBS(32),TPVBSTD,ESABTRD1(100),ESABTRD2(100),ESASTRD1(100) + ,ESASTRD2(100),TSABTRD1(100),TSABTRD2(100),TSASTRD1(100) + ,TSASTRD2(100) REAL SSABTRD1,SSABTRD2,SSASTRD1,SSASTRD2,ESANTRD(100) + ,EBASTRD1(100),EBASTRD2(100),TSANTRD(100) INTEGER NSASTRD1,NSASTRD2,ELSABTRD1(100),ELSABTRD2(100) + ,ELSASTRD1(100),ELSASTRD2(100),NSANTRD,NNABTRD,NBASTRD1,NBASTRD2 LOGICAL skim(10),cut(100),lct(21),ext(25),lev0(16),lhex,hhex(6) + ,l0z,l11,l12,ustop_hex,QB4,LCCDK(100),LCCDPI(100),LCCDPIOP(100) + ,CCDLOW(15),FPVECHI(80),FPVECT0(80),FPVECSTR(80),FPVEC2ND(80) + ,FPVECSFT(80),FPVECDNP(80),FPVECSTRD(80),LCCDPV(100),FPVBSHI(32) + ,FPVBST0(32),FPVBSSTR(32),FPVBSSFT(32),FPVBSNOTD * COMMON /PAWCR4/ run,event,splno,splrec,scal_tsp,scal_io2,scal_io3 + ,utbfld,skim,cut,icut,lay,tlay,sec,tsec,ssec,thex,t2sec,rssec1 + ,rssec2,rsdt1,rsdt2,trs,erd,lct,ext,lev0,lhex,hhex,l0z,l11,l12 + ,tavgz,tntg,rsintime,ustop_hex,ipiflg,cmupi1,tmuon1,emuon1 + ,emuonc1,dt1,emut1,emutc1,prod1,tmuav1,muped1,piped1,cmuend1 + ,cpiend1,tpion1,epion1,epionc1,ers,rrs,layv4,secv4,cos3d,pdc + ,ndclay,zf,elast,elast_cor,tk,tpi,ictpi,itgqualt,etg,rtg,tgx,tgy + ,tgz,ek,eic,ric,eicest,ptot,ptote,rtot,etot,b4strob,tdot2t,ictime + ,tds,rsz1d,rsz2d,rsz1d2,rsz2d2,zutout,radius,d0,phi0,tanl,nhxy + ,nhz,chidfxy,chidfz,rsxy1d,rsxy2d,ncross,rscrsd,rsr1,rsp1,rsz1 + ,rsr2,rsp2,rsz2,rscrsr,rscrsp,iicout,iicin,eicout,eicin,ticout + ,ticin,tdtic_i,tdtic_o,tdaic_i,tdaic_o,NRAW1,NPS1,PEDS1,NP1D1 + ,NP2D1,PEDD1,NRAW2,NPS2,PEDS2,NP1D2,NP2D2,PEDD2,TDFRAW1,TDFPS1 + ,TDFP1D1,TDFP2D1,TDFRAW2,TDFPS2,TDFP1D2,TDFP2D2,PR_RF,PR_RFZ1 + ,PR_RFZ2,PR_RFZ2B,LIKE,LIKE2,XIN,YIN,DB4,DB4TIP,DVXTIP,QB4,DVXPI + ,KPIANG,DEXTIP,UTCHNU1,UTCHNU2,UTCHNU3,LOGUTCQ,LOGUTCQI,rs_rate + ,cos3d1,radius1,d01,phi01,tanl1,chidfxy1,nhxy1,tgx1,tgy1,tkold + ,tpiold,NK_TG1,ELK_TG1,NPI_TG1,ELPI_TG1,NHITXY,NHITZ,NLAYER + ,DRIFT_TIME,RESID,ADCW,ZRES,NFOIL,ZZPOS,CL_EPRCOR,DELTARP,DELTARE + ,SIG_P,SIG_E,FDEDX,CL_RSDEDX,CL_RSDEDX2,CHIMAX_RSDEDX,LAY_CHIMAX + ,NPRPAT,DZUTSC14,L14ACC,ICODE_L14,ICODE_L14ACC,NK_TG,TK_TG,EK_TG + ,ELK_TG,NPI_TG,TPI_TG,EPI_TG,ELPI_TG,NPIOP_TG,TPIOP_TG,EPIOP_TG + ,ELPIOP_TG,CCDTK,CCDEK,CCDHK,CCDRK,LCCDK,CCDTPI,CCDEPI,CCDHPI + ,CCDRPI,LCCDPI,CCDTPIOP,CCDEPIOP,CCDHPIOP,CCDRPIOP,LCCDPIOP + ,NEDGE_TG,TEDGE_TG,EEDGE_TG,ELEDGE_TG,ETGVET_TG,NTGVET_TG + ,DIST_TARGF,EPIMAX_TG,TPIMAX_TG,NTRIK_TG,XTGT_TP,YTGT_TP,ZTGT_TP + ,DTGT_TP,RTDIF_TG,EPIMAXK,PHIVTX,DRPMAX,DRPMIN,NKCLUST,NKCLUST2 + ,PLN_SW,PLNV_SW,X_B4SW,Y_B4SW,IFAIL_KIC,DTK,DTPI,NK_TG_S,N1,N2,N3 + ,N4,N5,N6,N7,N8,N9,N10,N11,N12,N13,N14,N15,CCD1,CCD2,CCD3,CCD4 + ,CCD5,CCD6,CCD7,CCD8,CCD9,CCD10,CCD11,CCD12,CCD13,CCD14,CCD15 + ,CCDPF,CCDBIN,CCDLOW,NELVET,ISEV,ISTEV5,ISTEV5M1,T1MT2,EELEC1 + ,EELEC2,EELEC2C,TELEC,DTELEC,PROBTD,PROBTDN,IFAIL_TDFL,NHEXCT + ,ICODE_TDFL,EMUT_TDFL,EMUTC_TDFL,TMUAV_TDFL,DT_TDFL,DTMU_TDFL + ,CMUEND_TDFL,CPIEND_TDFL,CMUPI_TDFL,PROD_TDFL,NHIT_STHEX + ,ESUM_STHEX,EMAX_STHEX,ESTHEX,NNOFID,CHNOFID,TIME_FIDX,IFID_CODEX + ,NPVRD,TPVRD,EPVRD,ELEMPVRD,EPVRDPR,TZPVRD,EZPVRD,PVRDTE1,PVRDTR1 + ,PVRDTH1,PVRDTF1,PVRDTE2,PVRDTR2,PVRDTH2,PVRDTF2,NPVBV,TPVBV + ,EPVBV,ELEMPVBV,EPVBVPR,TZPVBV,EZPVBV,NPVBL,TPVBL,EPVBL,ELEMPVBL + ,EPVBLPR,TZPVBL,EZPVBL,NPVEC,TPVEC,EPVEC,HPVEC,RPVEC,ELEMPVEC + ,EPVECPR,EPVEUPR,FPVECHI,FPVECT0,FPVECSTR,FPVEC2ND,FPVECSFT + ,FPVECDNP,FPVECSTRD,EPVECGT,EPVEUGT,EPVECFERA,NPVIC,TPVIC,EPVIC + ,ELEMPVIC,EPVICPR,EPVVCPR,NPVCO,TPVCO,EPVCO,ELEMPVCO,EPVCOPR + ,NPVCM,TPVCM,EPVCM,ELEMPVCM,EPVCMPR,NPVTG,TPVTG,EPVTG,ELEMPVTG + ,EPVTGPR,CCDTPV,CCDEPV,CCDHPV,CCDRPV,LCCDPV,NPVBS,TPVBS,HPVBS + ,ELEMPVBS,TPVBSTD,FPVBSHI,FPVBST0,FPVBSSTR,FPVBSSFT,FPVBSNOTD + ,NSABTRD1,NSABTRD2,NSASTRD1,NSASTRD2,ESABTRD1,ESABTRD2,ESASTRD1 + ,ESASTRD2,TSABTRD1,TSABTRD2,TSASTRD1,TSASTRD2,SSABTRD1,SSABTRD2 + ,SSASTRD1,SSASTRD2,ELSABTRD1,ELSABTRD2,ELSASTRD1,ELSASTRD2 + ,NSANTRD,NNABTRD,NBASTRD1,NBASTRD2,ESANTRD,EBASTRD1,EBASTRD2 + ,TSANTRD * * *-- Enter user code here * call hfill(100,ptot,rtot,1.) if (.not.cut(39)) then ! kincut1 call hfill(101,ptot,rtot,1.) if (.not.cut(40)) then ! pscut1 call hfill(102,ptot,rtot,1.) c select kp2 region if (ptot.gt.195..and.ptot.lt.215.) then do i = 1, npvbv call hfill(300,tpvbv(i),epvbv(i),1.) enddo do i = 1, npvbl call hfill(301,tpvbl(i),epvbl(i),1.) enddo do i = 1, npvec call hfill(302,tpvec(i),epvec(i),1.) enddo do i = 1, npvrd call hfill(303,tpvrd(i),epvrd(i),1.) enddo do i = 1, npvtg call hfill(304,tpvtg(i),epvtg(i),1.) enddo do i = 1, npvco c write(6,*)tpvco(i), epvco(i), elempvco(i) if (elempvco(i).le.12) then call hfill(305,tpvco(i),epvco(i),1.) ii=elempvco(i) call hfill(1100+ii,tpvco(i),0.,1.) call hfill(3100+ii,epvco(i),0.,1.) else call hfill(306,tpvco(i),epvco(i),1.) ii=elempvco(i)-12 call hfill(2100+ii,tpvco(i),0.,1.) call hfill(4100+ii,epvco(i),0.,1.) endif enddo do i = 1, npvcm call hfill(307,tpvcm(i),epvcm(i),1.) enddo endif c slect tail region if (ptot.lt.195.) then do i = 1, npvbv call hfill(400,tpvbv(i),epvbv(i),1.) enddo do i = 1, npvbl call hfill(401,tpvbl(i),epvbl(i),1.) enddo do i = 1, npvec call hfill(402,tpvec(i),epvec(i),1.) enddo do i = 1, npvrd call hfill(403,tpvrd(i),epvrd(i),1.) enddo do i = 1, npvtg call hfill(404,tpvtg(i),epvtg(i),1.) enddo do i = 1, npvco c write(6,*)tpvco(i), epvco(i), elempvco(i) if (elempvco(i).le.12) then call hfill(405,tpvco(i),epvco(i),1.) ii=elempvco(i) call hfill(1000+ii,tpvco(i),0.,1.) call hfill(3000+ii,epvco(i),0.,1.) else call hfill(406,tpvco(i),epvco(i),1.) ii=elempvco(i)-12 call hfill(2000+ii,tpvco(i),0.,1.) call hfill(4000+ii,epvco(i),0.,1.) endif enddo do i = 1, npvcm call hfill(407,tpvcm(i),epvcm(i),1.) enddo endif endif endif pvmy=1 C my = 1. * END