Back to home page

Project CMSSW displayed by LXR

 
 

    


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

0001 
0002 c...pyfill all the data.

0003       subroutine uppyfill(idwtup,generate,xwgtup,pt,eta,st,y,pseta)
0004       implicit double precision(a-h, o-z)
0005         implicit integer(i-n)
0006 
0007         common/ptpass/ptmin,ptmax,crossmax,etamin,etamax,
0008      &  smin,smax,ymin,ymax,psetamin,psetamax
0009         common/histcol/inx
0010         common/counter/ibcstate,nev
0011       common/extraz/zfactor,zmin,zmax
0012       logical generate
0013 
0014 c      write (91,'(d20.10)') pseta

0015         
0016                  if(idwtup.eq.1.and.generate) then
0017                             call pyfill(1,pt,1.0d0)
0018                     call pyfill(2,eta,1.0d0)
0019                     call pyfill(3,st,1.0d0)
0020                     call pyfill(4,y,1.0d0)
0021                             call pyfill(5,pseta,1.0d0)
0022                                 call pyfill(6,zfactor,1.0d0)
0023 c...added output files. record differential distribution of rapidity

0024 c...with certain pt-cut.

0025                   if(pt .gt. 5.0d+0) then
0026                                     call pyfill(11,eta,1.0d0)
0027                                         call pyfill(41,st,1.0d0)
0028                     end if
0029                                 if(pt .gt. 2.0d+1) then
0030                                     call pyfill(12,eta,1.0d0)
0031                                         call pyfill(42,st,1.0d0)
0032                     end if
0033                                 if(pt .gt. 3.5d+1) then
0034                                     call pyfill(13,eta,1.0d0)
0035                                         call pyfill(43,st,1.0d0)
0036                                 end if
0037                                 if(pt .gt. 5.0d+1) then
0038                                     call pyfill(14,eta,1.0d0)
0039                                         call pyfill(44,st,1.0d0)
0040                                 end if
0041                                 if(pt .gt. 1.0d+2) then
0042                                     call pyfill(15,eta,1.0d0)
0043                                         call pyfill(45,st,1.0d0)
0044                     end if
0045 c...added output files. record differential distribution of pt

0046 c...and shat with certain rapidity-cut.

0047                   if(dabs(eta) .le. 5.0d-1) then
0048                                     call pyfill(21,pt,1.0d0)
0049                         call pyfill(31,st,1.0d0)
0050                     end if
0051                     if(dabs(eta) .le. 1.0d+0) then
0052                                     call pyfill(22,pt,1.0d0)
0053                         call pyfill(32,st,1.0d0)
0054                     end if
0055                     if(dabs(eta) .le. 1.5d+0) then
0056                                     call pyfill(23,pt,1.0d0)
0057                         call pyfill(33,st,1.0d0)
0058                     end if
0059                                 if(dabs(eta) .le. 2.0d+0) then
0060                                     call pyfill(24,pt,1.0d0)
0061                         call pyfill(34,st,1.0d0)
0062                     end if
0063                                 if(dabs(eta) .le. 2.5d+0) then
0064                                     call pyfill(25,pt,1.0d0)
0065                         call pyfill(35,st,1.0d0)
0066                     end if
0067 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

0068                                 if(dabs(eta).ge.5.0d-1 .and. dabs(eta).le.1.0d+0) then
0069                                     call pyfill(52,pt,1.0d0)
0070                     end if
0071                                 if(dabs(eta).ge.1.0d+0 .and. dabs(eta).le.1.5d+0) then
0072                                     call pyfill(53,pt,1.0d0)
0073                     end if
0074                                 if(dabs(eta).ge.1.5d+0 .and. dabs(eta).le.2.0d+0) then
0075                                     call pyfill(54,pt,1.0d0)
0076                     end if
0077                                 if(dabs(eta).ge.2.0d+0 .and. dabs(eta).le.2.5d+0) then
0078                                     call pyfill(55,pt,1.0d0)
0079                     end if
0080               else
0081                                 call pyfill(1,pt,xwgtup/nev*inx/(ptmax-ptmin))
0082                     call pyfill(2,eta,xwgtup/nev*inx/(etamax-etamin))
0083                     call pyfill(3,st,xwgtup/nev*inx/(smax-smin))
0084                     call pyfill(4,y,xwgtup/nev*inx/(ymax-ymin))
0085                             call pyfill(5,pseta,xwgtup/nev*inx/(psetamax-psetamin))
0086                                 call pyfill(6,zfactor,xwgtup/nev*inx/(zmax-zmin))
0087 c...added output files. record differential distribution of rapidity

0088 c...with certain pt-cut.

0089                   if(pt .ge. 5.0d+0) then
0090                             call pyfill(11,eta,xwgtup/nev*inx/(etamax-etamin))
0091                                         call pyfill(41,st,xwgtup/nev*inx/(smax-smin))
0092                     end if
0093                                 if(pt .ge. 2.0d+1) then
0094                             call pyfill(12,eta,xwgtup/nev*inx/(etamax-etamin))
0095                                         call pyfill(42,st,xwgtup/nev*inx/(smax-smin))
0096                     end if
0097                                 if(pt .ge. 3.5d+1) then
0098                             call pyfill(13,eta,xwgtup/nev*inx/(etamax-etamin))
0099                                     call pyfill(43,st,xwgtup/nev*inx/(smax-smin))
0100                                 end if
0101                                 if(pt .ge. 5.0d+1) then
0102                             call pyfill(14,eta,xwgtup/nev*inx/(etamax-etamin))
0103                                     call pyfill(44,st,xwgtup/nev*inx/(smax-smin))
0104                                 end if
0105                                 if(pt .ge. 1.0d+2) then
0106                             call pyfill(15,eta,xwgtup/nev*inx/(etamax-etamin))
0107                         call pyfill(45,st,xwgtup/nev*inx/(smax-smin))
0108                     end if
0109 c...added output files. record differential distribution of pt

0110 c...and shat with certain rapidity-cut.

0111                   if(dabs(eta) .le. 5.0d-1) then
0112                                     call pyfill(21,pt,xwgtup/nev*inx/(ptmax-ptmin))
0113                         call pyfill(31,st,xwgtup/nev*inx/(smax-smin))
0114                     end if
0115                     if(dabs(eta) .le. 1.0d+0) then
0116                                     call pyfill(22,pt,xwgtup/nev*inx/(ptmax-ptmin))
0117                         call pyfill(32,st,xwgtup/nev*inx/(smax-smin))
0118                     end if
0119                     if(dabs(eta) .le. 1.5d+0) then
0120                                     call pyfill(23,pt,xwgtup/nev*inx/(ptmax-ptmin))
0121                         call pyfill(33,st,xwgtup/nev*inx/(smax-smin))
0122                     end if
0123                                 if(dabs(eta) .le. 2.0d+0) then
0124                                     call pyfill(24,pt,xwgtup/nev*inx/(ptmax-ptmin))
0125                         call pyfill(34,st,xwgtup/nev*inx/(smax-smin))
0126                     end if
0127                                 if(dabs(eta) .le. 2.5d+0) then
0128                                     call pyfill(25,pt,xwgtup/nev*inx/(ptmax-ptmin))
0129                         call pyfill(35,st,xwgtup/nev*inx/(smax-smin))
0130                     end if
0131 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

0132                                 if(dabs(eta).ge.5.0d-1 .and. dabs(eta).le.1.0d+0) then
0133                                     call pyfill(52,pt,xwgtup/nev*inx/(ptmax-ptmin))
0134                     end if
0135                                 if(dabs(eta).ge.1.0d+0 .and. dabs(eta).le.1.5d+0) then
0136                                     call pyfill(53,pt,xwgtup/nev*inx/(ptmax-ptmin))
0137                     end if
0138                                 if(dabs(eta).ge.1.5d+0 .and. dabs(eta).le.2.0d+0) then
0139                                     call pyfill(54,pt,xwgtup/nev*inx/(ptmax-ptmin))
0140                     end if
0141                                 if(dabs(eta).ge.2.0d+0 .and. dabs(eta).le.2.5d+0) then
0142                                     call pyfill(55,pt,xwgtup/nev*inx/(ptmax-ptmin))
0143                     end if
0144 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

0145               end if
0146 
0147       end