Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:29:54

0001 
0002 c*******************************************************

0003 c... open files for recording grid generated by vegas and 

0004 c... for recording the intermediate running information

0005 c... e.g. (.grid)---> for grid; (.cs)---> for cross-section.

0006 c*******************************************************

0007 
0008       subroutine upopenfile(imix,imixtype,ibcstate,ioctet)
0009       implicit none
0010         integer imix,imixtype,ibcstate,ioctet
0011 
0012 c**************************************************

0013 c...output files to store the running information for the integrated

0014 c...cross-section. note here the p-wave states are only for gluon-

0015 c...gluon fusion subprocess.

0016 
0017         if(imix.eq.1) then
0018 c        open(unit=3,file='data/mix.cs',status='unknown')

0019          open(unit=3,file='mix.cs',status='unknown')
0020         end if
0021 
0022 c--------------------------------------------

0023 c...file about the generated grade by vegas.

0024         if(imix.eq.1) then
0025          if(imixtype.eq.1) then
0026 c         open(unit=36,file='data/grade1s0.grid',status='unknown')

0027 c         open(unit=37,file='data/grade3s1.grid',status='unknown')

0028 c         open(unit=38,file='data/grade1p1.grid',status='unknown')

0029 c         open(unit=39,file='data/grade3p0.grid',status='unknown')

0030 c         open(unit=46,file='data/grade3p1.grid',status='unknown')

0031 c         open(unit=47,file='data/grade3p2.grid',status='unknown')

0032 c         open(unit=48,file='data/grade81s.grid',status='unknown')

0033 c         open(unit=49,file='data/grade83s.grid',status='unknown')

0034           open(unit=36,file='grade1s0.grid',status='unknown')
0035           open(unit=37,file='grade3s1.grid',status='unknown')
0036           open(unit=38,file='grade1p1.grid',status='unknown')
0037           open(unit=39,file='grade3p0.grid',status='unknown')
0038           open(unit=46,file='grade3p1.grid',status='unknown')
0039           open(unit=47,file='grade3p2.grid',status='unknown')
0040           open(unit=48,file='grade81s.grid',status='unknown')
0041           open(unit=49,file='grade83s.grid',status='unknown')
0042          end if
0043          if(imixtype.eq.2) then
0044 c         open(unit=36,file='data/grade1s0.grid',status='unknown')

0045 c         open(unit=37,file='data/grade3s1.grid',status='unknown')

0046           open(unit=36,file='grade1s0.grid',status='unknown')
0047           open(unit=37,file='grade3s1.grid',status='unknown')
0048          end if
0049          if(imixtype.eq.3) then
0050 c         open(unit=38,file='data/grade1p1.grid',status='unknown')

0051 c         open(unit=39,file='data/grade3p0.grid',status='unknown')

0052 c         open(unit=46,file='data/grade3p1.grid',status='unknown')

0053 c         open(unit=47,file='data/grade3p2.grid',status='unknown')

0054 c         open(unit=48,file='data/grade81s.grid',status='unknown')

0055 c         open(unit=49,file='data/grade83s.grid',status='unknown')

0056           open(unit=38,file='grade1p1.grid',status='unknown')
0057           open(unit=39,file='grade3p0.grid',status='unknown')
0058           open(unit=46,file='grade3p1.grid',status='unknown')
0059           open(unit=47,file='grade3p2.grid',status='unknown')
0060           open(unit=48,file='grade81s.grid',status='unknown')
0061           open(unit=49,file='grade83s.grid',status='unknown')
0062          end if
0063         end if
0064 
0065 c******************************************************

0066 c...output files to store the running information for the integrated 

0067 c...cross-section. note here the p-wave states are only for gluon-

0068 c...gluon fusion subprocess.

0069 
0070       if(imix.eq.0) then
0071          if(ibcstate.eq.1.and.ioctet.eq.0) then
0072 c           open(unit=3,file='data/1s0.cs',status='unknown')

0073 c            open(unit=11,file='data/grade1s0.grid',status='unknown')

0074             open(unit=3,file='1s0.cs',status='unknown')
0075             open(unit=11,file='grade1s0.grid',status='unknown')
0076          end if
0077          if(ibcstate.eq.2.and.ioctet.eq.0) then
0078 c           open(unit=3,file='data/3s1.cs',status='unknown')

0079 c            open(unit=11,file='data/grade3s1.grid',status='unknown')

0080             open(unit=3,file='3s1.cs',status='unknown')
0081             open(unit=11,file='grade3s1.grid',status='unknown')
0082          end if
0083          if(ibcstate.eq.1.and.ioctet.eq.1) then
0084 c           open(unit=3,file='data/oct1s0.cs',status='unknown')

0085 c           open(unit=11,file='data/grade81s.grid',status='unknown')

0086            open(unit=3,file='oct1s0.cs',status='unknown')
0087            open(unit=11,file='grade81s.grid',status='unknown')
0088          end if
0089          if(ibcstate.eq.2.and.ioctet.eq.1) then
0090 c           open(unit=3,file='data/oct3s1.cs',status='unknown')

0091 c           open(unit=11,file='data/grade83s.grid',status='unknown')

0092            open(unit=3,file='oct3s1.cs',status='unknown')
0093            open(unit=11,file='grade83s.grid',status='unknown')
0094          end if
0095          if(ibcstate.eq.3) then
0096 c           open(unit=3,file='data/1p1.cs',status='unknown')

0097 c           open(unit=11,file='data/grade1p1.grid',status='unknown')

0098            open(unit=3,file='1p1.cs',status='unknown')
0099            open(unit=11,file='grade1p1.grid',status='unknown')
0100          end if
0101          if(ibcstate.eq.4) then
0102 c          open(unit=3,file='data/3p0.cs',status='unknown')

0103 c           open(unit=11,file='data/grade3p0.grid',status='unknown')

0104            open(unit=3,file='3p0.cs',status='unknown')
0105            open(unit=11,file='grade3p0.grid',status='unknown')
0106          end if
0107          if(ibcstate.eq.5) then
0108 c          open(unit=3,file='data/3p1.cs',status='unknown')

0109 c           open(unit=11,file='data/grade3p1.grid',status='unknown')

0110            open(unit=3,file='3p1.cs',status='unknown')
0111            open(unit=11,file='grade3p1.grid',status='unknown')
0112          end if
0113          if(ibcstate.eq.6) then
0114 c          open(unit=3,file='data/3p2.cs',status='unknown')

0115 c           open(unit=11,file='data/grade3p2.grid',status='unknown')

0116            open(unit=3,file='3p2.cs',status='unknown')
0117            open(unit=11,file='grade3p2.grid',status='unknown')
0118          end if
0119       end if
0120 
0121         end