lmpint.F

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

C **********************************************************************

      REAL FUNCTION lmpint(PEXTI,I)

C...This is the MINUIT routine PINTF.
CC        CALCULATES THE INTERNAL PARAMETER VALUE LMPINT CORRESPONDING
CC        TO THE EXTERNAL VALUE PEXTI FOR PARAMETER I.
CC
      COMMON
     1/lmmine/ erp(30)  ,ern(30)
     2/lmpari/ x(15)    ,xt(15)   ,dirin(15) ,maxint     ,npar
     3/lmpare/ u(30)              ,werr(30)  ,maxext     ,nu
     4/lmlimi/ alim(30) ,blim(30) ,lcode(30) ,lcorsp(30) ,limset
     5/lmvari/ v(15,15)
     7/lmfix / ipfix(15),xs(15)   ,xts(15)   ,dirins(15) ,npfix
     7/lmfix2/ grds(15) ,g2s(15)  ,gsteps(15),aberfs(15)
     c/lmcasc/ y(16)    ,jh       ,jl
     f/lmderi/ gin(30)  ,grd(15)  ,g2(15)    ,gstep(15)  ,aberf(15)
     g/lmsimv/ p(15,16) ,pstar(15),pstst(15) ,pbar(15)   ,prho(15)
     j/lmvart/ vt(15,15)
      COMMON
     6/lmunit/ isysrd   ,isyswr   ,isyspu
     8/lmtitl/ title(13),date(2)  ,isw(7)    ,nblock
     9/lmconv/ epsi     ,apsi     ,vtest     ,nstepq     ,nfcn ,nfcnmx
     a/lmcard/ cword    ,cword2   ,cword3    ,word7(7)
     b/lmmini/ amin     ,up       ,newmin    ,itaur      ,sigma,epsmac
      COMMON /lpflag/ lst3
      INTEGER lst3
      SAVE /lpflag/

      DATA big, small / 1.570796326795 , -1.570796326795 /
      igo = lcode(i)
      go to(100,200,300,400),igo
C--       IGO = 1  MEANS NO LIMITS
  100 lmpint = pexti
      go to 800
  200 CONTINUE
  300 CONTINUE
C--       IGO = 4  MEANS THERE ARE TWO LIMITS
  400 alimi = alim(i)
      blimi = blim(i)
      IF (pexti-alimi)  440,500,460
  440 a = small
  450 lmpint = a
      pexti = alimi + 0.5* (blimi-alimi) *(sin(a) +1.0)
      limset=1
      IF(lst3.GE.1) WRITE (isyswr,241)  i
      go to 800
  460 IF (blimi-pexti)  470,520,480
  470 a = big
      go to 450
  480 yy=2.0*(pexti-alimi)/(blimi-alimi) - 1.0
      lmpint = atan(yy/sqrt(1.0- yy**2) )
      go to 800
  500 lmpint = small
      go to 800
  520 lmpint = big
  800 RETURN
  241 FORMAT(' WARNING - VARIABLE',i3,' HAS BEEN BROUGHT BACK IN',
     +'SIDE LIMITS BY LMPINT.')
      END