pyfill.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pyfill.f

 
C*********************************************************************
 
C...PYFILL
C...Fills entry in histogram.
 
      SUBROUTINE pyfill(ID,X,W)
 
C...Double precision declaration.
      IMPLICIT DOUBLE PRECISION(a-h, o-z)
      IMPLICIT INTEGER(i-n)
C...Commonblock.
      common/pybins/ihist(4),indx(1000),bin(20000)
      SAVE /pybins/
 
C...Find initial address in memory. Increase number of entries.
      IF(id.LE.0.OR.id.GT.ihist(1)) CALL pyerrm(28,
     &'(PYFILL:) not allowed histogram number')
      is=indx(id)
      IF(is.EQ.0) CALL pyerrm(28,
     &'(PYFILL:) filling unbooked histogram')
      bin(is+5)=bin(is+5)+1d0
 
C...Find bin in x, including under/overflow, and fill.
      IF(x.LT.bin(is+2)) THEN
        bin(is+6)=bin(is+6)+w
      ELSEIF(x.GE.bin(is+3)) THEN
        bin(is+8)=bin(is+8)+w
      ELSE
        bin(is+7)=bin(is+7)+w
        ix=(x-bin(is+2))/bin(is+4)
        ix=max(0,min(nint(bin(is+1))-1,ix))
        bin(is+9+ix)=bin(is+9+ix)+w
      ENDIF
 
      RETURN
      END