COSMOS v7.655  COSMOSv7655
(AirShowerMC)
Functions/Subroutines
mnmnos.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine mnmnos (FCN, FUTIL)
 

Function/Subroutine Documentation

◆ mnmnos()

subroutine mnmnos ( external  FCN,
external  FUTIL 
)

Definition at line 10 of file mnmnos.f.

References a, b, c, d, e, f, g, h, softenpik::half, i, j, m, maxp, mnmatu(), mnmnot(), mnprin(), n, o, p, parameter(), to, up(), x, xs, and z.

Referenced by mnexcm().

10 *
11 * $Id: d506dp.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
12 *
13 * $Log: d506dp.inc,v $
14 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
15 * Minuit
16 *
17 *
18 *
19 *
20 * d506dp.inc
21 *
22 C ************ DOUBLE PRECISION VERSION *************
23  IMPLICIT DOUBLE PRECISION (a-h,o-z)
24 CC Performs a MINOS error analysis on those parameters for
25 CC which it is requested on the MINOS command by calling
26 CC MNMNOT for each parameter requested.
27 *
28 * $Id: d506cm.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
29 *
30 * $Log: d506cm.inc,v $
31 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
32 * Minuit
33 *
34 *
35 *
36 *
37 * d506cm.inc
38 *
39  parameter(mne=100 , mni=50)
40  parameter(mnihl=mni*(mni+1)/2)
41  CHARACTER*10 cpnam
42  COMMON
43  1/mn7nam/ cpnam(mne)
44  2/mn7ext/ u(mne) ,alim(mne) ,blim(mne)
45  3/mn7err/ erp(mni) ,ern(mni) ,werr(mni) ,globcc(mni)
46  4/mn7inx/ nvarl(mne) ,niofex(mne),nexofi(mni)
47  5/mn7int/ x(mni) ,xt(mni) ,dirin(mni)
48  6/mn7fx2/ xs(mni) ,xts(mni) ,dirins(mni)
49  7/mn7der/ grd(mni) ,g2(mni) ,gstep(mni) ,gin(mne) ,dgrd(mni)
50  8/mn7fx3/ grds(mni) ,g2s(mni) ,gsteps(mni)
51  9/mn7fx1/ ipfix(mni) ,npfix
52  a/mn7var/ vhmat(mnihl)
53  b/mn7vat/ vthmat(mnihl)
54  c/mn7sim/ p(mni,mni+1),pstar(mni),pstst(mni) ,pbar(mni),prho(mni)
55 C
56  parameter(maxdbg=10, maxstk=10, maxcwd=20, maxp=30, maxcpt=101)
57  parameter(zero=0.0, one=1.0, half=0.5)
58  COMMON
59  d/mn7npr/ maxint ,npar ,maxext ,nu
60  e/mn7iou/ isysrd ,isyswr ,isyssa ,npagwd ,npagln ,newpag
61  e/mn7io2/ istkrd(maxstk) ,nstkrd ,istkwr(maxstk) ,nstkwr
62  f/mn7tit/ cfrom ,cstatu ,ctitl ,cword ,cundef ,cvrsn ,covmes
63  g/mn7flg/ isw(7) ,idbg(0:maxdbg) ,nblock ,icomnd
64  h/mn7min/ amin ,up ,edm ,fval3 ,epsi ,apsi ,dcovar
65  i/mn7cnv/ nfcn ,nfcnmx ,nfcnlc ,nfcnfr ,itaur,istrat,nwrmes(2)
66  j/mn7arg/ word7(maxp)
67  k/mn7log/ lwarn ,lrepor ,limset ,lnolim ,lnewmn ,lphead
68  l/mn7cns/ epsmac ,epsma2 ,vlimlo ,vlimhi ,undefi ,bigedm,updflt
69  m/mn7rpt/ xpt(maxcpt) ,ypt(maxcpt)
70  n/mn7cpt/ chpt(maxcpt)
71  o/mn7xcr/ xmidcr ,ymidcr ,xdircr ,ydircr ,ke1cr ,ke2cr
72  CHARACTER ctitl*50, cword*(maxcwd), cundef*10, cfrom*8,
73  + cvrsn*6, covmes(0:3)*22, cstatu*10, chpt*1
74  LOGICAL lwarn, lrepor, limset, lnolim, lnewmn, lphead
75  EXTERNAL fcn,futil
76  IF (npar .LE. 0) GO TO 700
77  ngood = 0
78  nbad = 0
79  nfcnmi = nfcn
80 C . loop over parameters requested
81  DO 570 knt= 1, npar
82  IF (int(word7(2)) .EQ. 0) THEN
83  ilax = nexofi(knt)
84  ELSE
85  IF (knt .GE. 7) GO TO 580
86  ilax = int(word7(knt+1))
87  IF (ilax .EQ. 0) GO TO 580
88  IF (ilax .GT. 0 .AND. ilax .LE. nu) THEN
89  IF (niofex(ilax) .GT. 0) GO TO 565
90  ENDIF
91  WRITE (isyswr,564) ilax
92  564 FORMAT (' PARAMETER NUMBER ',i5,' NOT VARIABLE. IGNORED.')
93  GO TO 570
94  ENDIF
95  565 CONTINUE
96 C calculate one pair of M E's
97  ilax2 = 0
98  CALL mnmnot(fcn,ilax,ilax2,val2pl,val2mi,futil)
99  IF (lnewmn) GO TO 650
100 C update NGOOD and NBAD
101  iin = niofex(ilax)
102  IF (erp(iin) .GT. zero) THEN
103  ngood=ngood+1
104  ELSE
105  nbad=nbad+1
106  ENDIF
107  IF (ern(iin) .LT. zero) THEN
108  ngood=ngood+1
109  ELSE
110  nbad=nbad+1
111  ENDIF
112  570 CONTINUE
113 C end of loop . . . . . . .
114  580 CONTINUE
115 C . . . . printout final values .
116  cfrom = 'MINOS '
117  nfcnfr = nfcnmi
118  cstatu= 'UNCHANGED '
119  IF (ngood.EQ.0.AND.nbad.EQ.0) GO TO 700
120  IF (ngood.GT.0.AND.nbad.EQ.0) cstatu='SUCCESSFUL'
121  IF (ngood.EQ.0.AND.nbad.GT.0) cstatu='FAILURE '
122  IF (ngood.GT.0.AND.nbad.GT.0) cstatu='PROBLEMS '
123  IF (isw(5) .GE. 0) CALL mnprin(4,amin)
124  IF (isw(5) .GE. 2) CALL mnmatu(0)
125  GO TO 900
126 C . . . new minimum found . . . .
127  650 CONTINUE
128  cfrom = 'MINOS '
129  nfcnfr = nfcnmi
130  cstatu= 'NEW MINIMU'
131  IF (isw(5) .GE. 0) CALL mnprin(4,amin)
132  WRITE (isyswr,675)
133  675 FORMAT(/50h new minimum found. go back to minimization step./1h ,
134  +60(1h=)/60x,1hv/60x,1hv/60x,1hv/57x,7hvvvvvvv/58x,5hvvvvv/59x,
135  +3hvvv/60x,1hv//)
136  GO TO 900
137  700 WRITE (isyswr,'(A)') ' THERE ARE NO MINOS ERRORS TO CALCULATE.'
138  900 RETURN
parameter
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
z
nodes z
Definition: cstdatmos0.fNew2.h:93
up
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data h g *is for param c g data up(2, 1)/7.0d0/
e
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
i
nodes i
Definition: cstdatmos0.fNew2.h:93
nrl_atmos::h
real(4), dimension(:), allocatable, save h
Definition: cNRLAtmos.f:28
mnmnot
subroutine mnmnot(FCN, ILAX, ILAX2, VAL2PL, VAL2MI, FUTIL)
Definition: mnmnot.f:10
mnprin
subroutine mnprin(INKODE, FVAL)
Definition: mnprin.f:10
xs
! constants thru Cosmos real ! if multiplied to deg radian Torad ! light velocity m sec ! infinty ! kg m2 *Togpcm2 g cm2 ! g cm2 *Tokgpm2 kg m2 ! cm *Tom m ! m *Tocm cm ! g cm3 *Tokgpm3 kg m3 ! kg m3 *Togpcm3 g cm3 ! sec *Tonsec nsec ! Tesla m ! Avogadro *A2deninv ! mfp *n * xs
Definition: Zglobalc.h:18
p
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
Definition: cblkHeavy.h:7
j
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
to
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer to
Definition: Zfit.h:15
o
struct ob o[NpMax]
Definition: Zprivate.h:34
d
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
mnmatu
subroutine mnmatu(KODE)
Definition: mnmatu.f:10
m
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
Definition: ZlatfitBD.h:35
nrl_atmos::a
real(4), save a
Definition: cNRLAtmos.f:20
g
dE dx *! Nuc Int sampling table g
Definition: cblkMuInt.h:130
nrl_atmos::b
real(4), save b
Definition: cNRLAtmos.f:21
maxp
!onst int maxp
Definition: Zprivate.h:3
n
integer n
Definition: Zcinippxc.h:1
softenpik::half
integer, parameter half
Definition: csoftenPiK.f:108
x
! structure defining a particle at production ! Basic idea of what is to be contained in ! the particle structue is that dynamical ones should be included those derivable from the particle code ! is not included ******************************************************type fmom momentum sequence union map real e endmap map real * x
Definition: Zptcl.h:21
f
dE dx *! Nuc Int sampling table f
Definition: cblkMuInt.h:130
c
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130
Here is the call graph for this function:
Here is the caller graph for this function: