subroutine NtupleInit(filename) implicit none save include 'hbook.inc' include 'simulate.inc' character*80 filename,directory character*16 NtupleTag(80),name,title integer*4 m,io,recl,bank,id,status parameter(recl = 1024) parameter(bank = 8000) parameter(io = 29) parameter(name = 'SimcNtuple') parameter(title = 'SIMTUPLE') NtupleID = defaultID id = NtupleID NtupleIO = io NtupleName = name call HCDIR(directory,'R') !CERNLIB read current directory call HROPEN(io,name,filename,'N',recl,status) !CERNLIB !directory set to "//TUPLE" if (status.ne.0)then write(6,*) 'HROPEN error: istat=',status stop endif m = 0 m = m+1 NtupleTag(m) = 'hsdelta' ! 1 m = m+1 NtupleTag(m) = 'hsyptar' ! 2 m = m+1 NtupleTag(m) = 'hsxptar' ! 3 m = m+1 NtupleTag(m) = 'hsytar' ! 4 m = m+1 NtupleTag(m) = 'hsxfp' ! 5 m = m+1 NtupleTag(m) = 'hsxpfp' ! 6 m = m+1 NtupleTag(m) = 'hsyfp' ! 7 m = m+1 NtupleTag(m) = 'hsypfp' ! 8 m = m+1 NtupleTag(m) = 'hsdeltai' ! 9 m = m+1 NtupleTag(m) = 'hsyptari' ! 10 m = m+1 NtupleTag(m) = 'hsxptari' ! 11 m = m+1 NtupleTag(m) = 'hsytari' ! 12 m = m+1 NtupleTag(m) = 'ssdelta' ! 13 m = m+1 NtupleTag(m) = 'ssyptar' ! 14 m = m+1 NtupleTag(m) = 'ssxptar' ! 15 m = m+1 NtupleTag(m) = 'ssytar' ! 16 m = m+1 NtupleTag(m) = 'ssxfp' ! 17 m = m+1 NtupleTag(m) = 'ssxpfp' ! 18 m = m+1 NtupleTag(m) = 'ssyfp' ! 19 m = m+1 NtupleTag(m) = 'ssypfp' ! 20 m = m+1 NtupleTag(m) = 'ssdeltai' ! 21 m = m+1 NtupleTag(m) = 'ssyptari' ! 22 m = m+1 NtupleTag(m) = 'ssxptari' ! 23 m = m+1 NtupleTag(m) = 'ssytari' ! 24 m = m+1 NtupleTag(m) = 'q' ! 25 m = m+1 NtupleTag(m) = 'nu' ! 26 m = m+1 NtupleTag(m) = 'Q2' ! 27 m = m+1 NtupleTag(m) = 'W' ! 28 m = m+1 NtupleTag(m) = 'epsilon' ! 29 m = m+1 NtupleTag(m) = 'Em' ! 30 m = m+1 NtupleTag(m) = 'Pm' ! 31 m = m+1 NtupleTag(m) = 'thetapq' ! 32 m = m+1 NtupleTag(m) = 'phipq' ! 33 if (doing_pion .or. doing_kaon .or. doing_delta) then m = m+1 NtupleTag(m) = 'missmass' ! 34 m = m+1 NtupleTag(m) = 'mmnuc' ! 35 m = m+1 NtupleTag(m) = 'phad' ! 36 m = m+1 NtupleTag(m) = 't' ! 37 m = m+1 NtupleTag(m) = 'pmpar' ! 38 m = m+1 NtupleTag(m) = 'pmper' ! 39 m = m+1 NtupleTag(m) = 'pmoop' ! 40 m = m+1 NtupleTag(m) = 'fry' ! 41 !+y is up. m = m+1 NtupleTag(m) = 'radphot' ! 42 m = m+1 NtupleTag(m) = 'pfermi' ! 43 m = m+1 NtupleTag(m) = 'siglab' ! 44 m = m+1 NtupleTag(m) = 'sigcm' ! 45 m = m+1 NtupleTag(m) = 'Weight' ! 46 m = m+1 NtupleTag(m) = 'jacobian' ! 47 m = m+1 NtupleTag(m) = 'decdist' ! 48 m = m+1 NtupleTag(m) = 'Mhadron' ! 49 m = m+1 NtupleTag(m) = 'pdotqhat' ! 50 m = m+1 NtupleTag(m) = 'Q2i' ! 51 m = m+1 NtupleTag(m) = 'Wi' ! 52 m = m+1 NtupleTag(m) = 'ti' ! 53 m = m+1 NtupleTag(m) = 'thetapqi' ! 54 m = m+1 NtupleTag(m) = 'phipqi' ! 55 m = m+1 NtupleTag(m) = 'thetacm' ! 56 m = m+1 NtupleTag(m) = 'phicm' ! 57 m = m+1 NtupleTag(m) = 'phadcm' ! 58 if (doing_kaon) then m = m+1 NtupleTag(m) = 'saghai' ! 59 m = m+1 NtupleTag(m) = 'factor' ! 60 endif else if (doing_semi.or.doing_rho) then m = m+1 NtupleTag(m) = 'missmass' ! 34 <- Wprime for semi-inclusive folks m = m+1 NtupleTag(m) = 'ppi' ! 35 m = m+1 NtupleTag(m) = 't' ! 36 m = m+1 NtupleTag(m) = 'fry' ! 37 !+y is up. m = m+1 NtupleTag(m) = 'radphot' ! 38 m = m+1 NtupleTag(m) = 'siglab' ! 39 m = m+1 NtupleTag(m) = 'sigcent' ! 40 m = m+1 NtupleTag(m) = 'Weight' ! 41 m = m+1 NtupleTag(m) = 'decdist' ! 42 m = m+1 NtupleTag(m) = 'Mhadron' ! 43 m = m+1 NtupleTag(m) = 'z' ! 44 m = m+1 NtupleTag(m) = 'zi' ! 45 m = m+1 NtupleTag(m) = 'pt2' ! 46 m = m+1 NtupleTag(m) = 'pt2i' ! 47 m = m+1 NtupleTag(m) = 'xbj' ! 48 m = m+1 NtupleTag(m) = 'xbji' ! 49 m = m+1 NtupleTag(m) = 'thqi' ! 50 m = m+1 NtupleTag(m) = 'sighad' ! 51 m = m+1 NtupleTag(m) = 'jacobian' ! 52 m = m+1 NtupleTag(m) = 'centjac' ! 53 m = m+1 NtupleTag(m) = 'pfermi' ! 54 m = m+1 NtupleTag(m) = 'xfermi' ! 55 if(doing_rho) then m = m+1 NtupleTag(m) = 'Mrho' ! 56 m = m+1 NtupleTag(m) = 'Thrho' ! 57 endif else if (doing_hyd_elast .or. doing_deuterium .or. doing_heavy) then m = m+1 NtupleTag(m) = 'corrsing' ! 34 m = m+1 NtupleTag(m) = 'Pmx' ! 35 !for Heepcheck m = m+1 NtupleTag(m) = 'Pmy' ! 36 !for Heepcheck m = m+1 NtupleTag(m) = 'Pmz' ! 37 !for Heepcheck m = m+1 NtupleTag(m) = 'PmPar' ! 38 m = m+1 NtupleTag(m) = 'PmPer' ! 39 m = m+1 NtupleTag(m) = 'PmOop' ! 40 m = m+1 NtupleTag(m) = 'fry' ! 41 !+y is up. m = m+1 NtupleTag(m) = 'radphot' ! 42 m = m+1 NtupleTag(m) = 'sigcc' ! 43 m = m+1 NtupleTag(m) = 'Weight' ! 44 endif ! else !used to be the if (doing_phsp) option. ! m=m+1 ! NtupleTag(m)='gd' ! m=m+1 ! NtupleTag(m)='gt' ! m=m+1 ! NtupleTag(m)='gp' ! m=m+1 ! NtupleTag(m)='gy' ! m=m+1 ! NtupleTag(m)='rd' ! m=m+1 ! NtupleTag(m)='rt' ! m=m+1 ! NtupleTag(m)='rp' ! m=m+1 ! NtupleTag(m)='ry' ! m=m+1 ! NtupleTag(m)='w' ! endif NtupleSize = m call HBOOKN(id,title,NtupleSize,name,bank,NtupleTag) !create Ntuple call HCDIR(NtupleDirectory,'R') !record Ntuple directory call HCDIR(directory,' ') !reset CERNLIB directory return end •r=•r=•r=íA5s¯?^Y?^Y?€žr=žr=P ÀA@T¯?•r=•r=    ¤¹9o?¹9o?¹9o?˜3Èÿ¡½ˆ®?¤ùˆ=¤ùˆ=/usr/local/cern/š=§€&£ˆ=&£ˆ=&£ˆ=/"Ü$R¤2ÏH? ÕH?2ÏH? ÕH?dT¦gÿ¡2d?ñ‰=ñ‰=/w/cache±|ºÿ¡q®®?ñ‰=ñ‰=/u/apps²|ºÿ¡q®®?ñ‰=ñ‰=/u/site³|ºÿ¡ Lª?ñ‰=ñ‰=/u/scratch´|ºÿ¡ÛD­?ñ‰=ñ‰=/w/workµ|ºÿ¡ 2d?ñ‰=ñ‰=/u/mirror¶|ºÿ¡Kª?ñ‰=ñ‰=/w/mss·|º¤‡‹=I‡‹=‡‹=I‡‹=°|º