COSMOS v7.655  COSMOSv7655
(AirShowerMC)
mnwerr.f
Go to the documentation of this file.
1 *
2 * $Id: mnwerr.F,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
3 *
4 * $Log: mnwerr.F,v $
5 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
6 * Minuit
7 *
8 *
9  SUBROUTINE mnwerr
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 Calculates the WERR, external parameter errors,
25 CC and the global correlation coefficients, to be called
26 CC whenever a new covariance matrix is available.
27 CC
28 *
29 * $Id: d506cm.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
30 *
31 * $Log: d506cm.inc,v $
32 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
33 * Minuit
34 *
35 *
36 *
37 *
38 * d506cm.inc
39 *
40  parameter(mne=100 , mni=50)
41  parameter(mnihl=mni*(mni+1)/2)
42  CHARACTER*10 CPNAM
43  COMMON
44  1/mn7nam/ cpnam(mne)
45  2/mn7ext/ u(mne) ,alim(mne) ,blim(mne)
46  3/mn7err/ erp(mni) ,ern(mni) ,werr(mni) ,globcc(mni)
47  4/mn7inx/ nvarl(mne) ,niofex(mne),nexofi(mni)
48  5/mn7int/ x(mni) ,xt(mni) ,dirin(mni)
49  6/mn7fx2/ xs(mni) ,xts(mni) ,dirins(mni)
50  7/mn7der/ grd(mni) ,g2(mni) ,gstep(mni) ,gin(mne) ,dgrd(mni)
51  8/mn7fx3/ grds(mni) ,g2s(mni) ,gsteps(mni)
52  9/mn7fx1/ ipfix(mni) ,npfix
53  a/mn7var/ vhmat(mnihl)
54  b/mn7vat/ vthmat(mnihl)
55  c/mn7sim/ p(mni,mni+1),pstar(mni),pstst(mni) ,pbar(mni),prho(mni)
56 C
57  parameter(maxdbg=10, maxstk=10, maxcwd=20, maxp=30, maxcpt=101)
58  parameter(zero=0.0, one=1.0, half=0.5)
59  COMMON
60  d/mn7npr/ maxint ,npar ,maxext ,nu
61  e/mn7iou/ isysrd ,isyswr ,isyssa ,npagwd ,npagln ,newpag
62  e/mn7io2/ istkrd(maxstk) ,nstkrd ,istkwr(maxstk) ,nstkwr
63  f/mn7tit/ cfrom ,cstatu ,ctitl ,cword ,cundef ,cvrsn ,covmes
64  g/mn7flg/ isw(7) ,idbg(0:maxdbg) ,nblock ,icomnd
65  h/mn7min/ amin ,up ,edm ,fval3 ,epsi ,apsi ,dcovar
66  i/mn7cnv/ nfcn ,nfcnmx ,nfcnlc ,nfcnfr ,itaur,istrat,nwrmes(2)
67  j/mn7arg/ word7(maxp)
68  k/mn7log/ lwarn ,lrepor ,limset ,lnolim ,lnewmn ,lphead
69  l/mn7cns/ epsmac ,epsma2 ,vlimlo ,vlimhi ,undefi ,bigedm,updflt
70  m/mn7rpt/ xpt(maxcpt) ,ypt(maxcpt)
71  n/mn7cpt/ chpt(maxcpt)
72  o/mn7xcr/ xmidcr ,ymidcr ,xdircr ,ydircr ,ke1cr ,ke2cr
73  CHARACTER CTITL*50, CWORD*(maxcwd), CUNDEF*10, CFROM*8,
74  + cvrsn*6, covmes(0:3)*22, cstatu*10, chpt*1
75  LOGICAL LWARN, LREPOR, LIMSET, LNOLIM, LNEWMN, LPHEAD
76 C calculate external error if v exists
77  IF (isw(2) .GE. 1) THEN
78  DO 100 l= 1, npar
79  ndex = l*(l+1)/2
80  dx = sqrt(abs(vhmat(ndex)*up))
81  i = nexofi(l)
82  IF (nvarl(i) .GT. 1) THEN
83  al = alim(i)
84  ba = blim(i) - al
85  du1 = al + 0.5 *(sin(x(l)+dx) +1.0) * ba - u(i)
86  du2 = al + 0.5 *(sin(x(l)-dx) +1.0) * ba - u(i)
87  IF (dx .GT. 1.0) du1 = ba
88  dx = 0.5 * (abs(du1) + abs(du2))
89  ENDIF
90  werr(l) = dx
91  100 CONTINUE
92  ENDIF
93 C global correlation coefficients
94  IF (isw(2) .GE. 1) THEN
95  DO 130 i= 1, npar
96  globcc(i) = 0.
97  k1 = i*(i-1)/2
98  DO 130 j= 1, i
99  k = k1 + j
100  p(i,j) = vhmat(k)
101  130 p(j,i) = p(i,j)
102  CALL mnvert(p,maxint,maxint,npar,ierr)
103  IF (ierr .EQ. 0) THEN
104  DO 150 iin= 1, npar
105  ndiag = iin*(iin+1)/2
106  denom = p(iin,iin)*vhmat(ndiag)
107  IF (denom.LE.one .AND. denom.GE.zero) THEN
108  globcc(iin) = 0.
109  ELSE
110  globcc(iin) = sqrt(1.0-1.0/denom)
111  ENDIF
112  150 CONTINUE
113  ENDIF
114  ENDIF
115  RETURN
116  END
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
dx
integer npitbl real *nx dx real dx
Definition: Zcinippxc.h:10
mnwerr
subroutine mnwerr
Definition: mnwerr.f:10
mnvert
subroutine mnvert(A, L, M, N, IFAIL)
Definition: mnvert.f:25
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
o
struct ob o[NpMax]
Definition: Zprivate.h:34
d
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
a
nodes a
Definition: cstdatmos0.fNew2.h:93
b
dE dx *! Nuc Int sampling table b
Definition: cblkMuInt.h:130
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
h
dE dx *! Nuc Int sampling table h
Definition: cblkMuInt.h:130
g
dE dx *! Nuc Int sampling table g
Definition: cblkMuInt.h:130
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