pyfact.f

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

 
C*********************************************************************
 
C...PYFACT
C...Multiplies histogram contents by factor.
 
      SUBROUTINE pyfact(ID,F)
 
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. Multiply all contents bins.
      IF(id.LE.0.OR.id.GT.ihist(1)) CALL pyerrm(28,
     &'(PYFACT:) not allowed histogram number')
      is=indx(id)
      IF(is.EQ.0) CALL pyerrm(28,
     &'(PYFACT:) scaling unbooked histogram')
      DO 100 ix=is+6,is+8+nint(bin(is+1))
        bin(ix)=f*bin(ix)
  100 CONTINUE
 
      RETURN
      END