40 |
* (All the needed quantities for couples and singlets can be evaluated |
* (All the needed quantities for couples and singlets can be evaluated |
41 |
* by means of the routine xyz_PAM) |
* by means of the routine xyz_PAM) |
42 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
43 |
|
|
44 |
|
|
45 |
REAL*8 AL |
REAL*8 AL |
46 |
REAL*8 XM,YM,ZM |
REAL*8 XM,YM,ZM |
47 |
REAL*8 XM_A,YM_A!,ZM_A |
REAL*8 XM_A,YM_A!,ZM_A |
48 |
REAL*8 XM_B,YM_B!,ZM_B |
REAL*8 XM_B,YM_B!,ZM_B |
49 |
REAL*8 RESX,RESY |
REAL*8 RESX,RESY |
50 |
|
REAL*8 TAILX,TAILY |
51 |
REAL*8 XGOOD,YGOOD |
REAL*8 XGOOD,YGOOD |
52 |
c REAL*8 DEDXTRK !(1) |
REAL*8 DEDXTRK_X |
53 |
REAL*8 DEDXTRK_X !(1) |
REAL*8 DEDXTRK_Y |
|
REAL*8 DEDXTRK_Y !(1) |
|
54 |
REAL*8 ZINI ! ***PP*** |
REAL*8 ZINI ! ***PP*** |
55 |
REAL*8 PFIXED ! ***PP*** |
REAL*8 PFIXED ! ***PP*** |
56 |
REAL*8 CHI2 |
REAL*8 CHI2 |
59 |
REAL*8 COV |
REAL*8 COV |
60 |
REAL*8 FACT |
REAL*8 FACT |
61 |
INTEGER NSTEP |
INTEGER NSTEP |
62 |
INTEGER IDCAND !(2) |
INTEGER IDCAND |
63 |
INTEGER TRACKMODE |
INTEGER TRACKMODE |
64 |
INTEGER ISTEPMIN |
INTEGER ISTEPMIN |
65 |
COMMON/TRACK/ |
COMMON/TRACK/ |
68 |
$ ,XM_A(NPLANES),YM_A(NPLANES)!,ZM_A(NPLANES) !single cl => segment AB |
$ ,XM_A(NPLANES),YM_A(NPLANES)!,ZM_A(NPLANES) !single cl => segment AB |
69 |
$ ,XM_B(NPLANES),YM_B(NPLANES)!,ZM_B(NPLANES) ! |
$ ,XM_B(NPLANES),YM_B(NPLANES)!,ZM_B(NPLANES) ! |
70 |
$ ,RESX(NPLANES),RESY(NPLANES) !spatial resolution |
$ ,RESX(NPLANES),RESY(NPLANES) !spatial resolution |
71 |
|
$ ,TAILX(NPLANES),TAILY(NPLANES) !tails in the resolution |
72 |
$ ,XGOOD(NPLANES),YGOOD(NPLANES) !mask of included views |
$ ,XGOOD(NPLANES),YGOOD(NPLANES) !mask of included views |
73 |
c $ ,DEDXTRK(NPLANES) !energy release |
c $ ,DEDXTRK(NPLANES) !energy release |
74 |
$ ,DEDXTRK_X(NPLANES),DEDXTRK_Y(NPLANES) !(1) |
$ ,DEDXTRK_X(NPLANES),DEDXTRK_Y(NPLANES) |
75 |
$ ,ZINI |
$ ,ZINI |
76 |
$ ,PFIXED |
$ ,PFIXED |
77 |
$ ,CHI2 |
$ ,CHI2 |
88 |
cccccc 17/08/2006 modified by elena vannuccini ---> (2) |
cccccc 17/08/2006 modified by elena vannuccini ---> (2) |
89 |
c****************************************************** |
c****************************************************** |
90 |
|
|
91 |
|
|
92 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
93 |
* The FIT info includes: |
* The FIT info includes: |
94 |
* - REDUCED chi2 |
* - REDUCED chi2 |
113 |
|
|
114 |
|
|
115 |
|
|
|
|
|
|
|
|
116 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
117 |
* Other commons needed by the minimization procedures: |
* Other commons needed by the minimization procedures: |
118 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
119 |
* computation of chi^2 and its derivatives |
* computation of chi^2 and its derivatives |
120 |
REAL*8 CHI2D,CHI2DD !first and second derivatives of chi^2 |
REAL*8 DAL |
121 |
|
REAL*8 CHI2D,CHI2DOLD,CHI2DD !first and second derivatives of chi^2 |
122 |
REAL*8 DXDAL,DYDAL !first derivatives of x and y on AL-pha |
REAL*8 DXDAL,DYDAL !first derivatives of x and y on AL-pha |
123 |
REAL*8 STEPAL ! |
REAL*8 STEPAL ! |
124 |
***PP*** REAL*8 ZINI !z coordinate of reference plane |
***PP*** REAL*8 ZINI !z coordinate of reference plane |
125 |
COMMON/MINICOM/ |
COMMON/MINICOM/ |
126 |
$ CHI2D(5),CHI2DD(5,5) |
$ DAL(5) |
127 |
|
$ ,CHI2D(5),CHI2DOLD(5),CHI2DD(5,5) |
128 |
$ ,DXDAL(NPLANES,5),DYDAL(NPLANES,5) |
$ ,DXDAL(NPLANES,5),DYDAL(NPLANES,5) |
129 |
$ ,STEPAL(5) |
$ ,STEPAL(5) |
130 |
***PP*** $ ,ZINI |
***PP*** $ ,ZINI |
131 |
|
|
132 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
133 |
* convergence conditions |
* convergence conditions |
134 |
REAL*8 ALTOL !tolerances on AL-pha parameters |
REAL*8 ALTOL !tolerances on AL-pha parameters |
138 |
$ ALTOL(5) |
$ ALTOL(5) |
139 |
$ ,ALMAX(5),ALMIN(5) |
$ ,ALMAX(5),ALMIN(5) |
140 |
$ ,ISTEPMAX |
$ ,ISTEPMAX |
141 |
|
|
142 |
|
|
143 |
* ------------------------------------------------------------------------- |
* ------------------------------------------------------------------------- |
144 |
* variables for the tracking (routine posxy) |
* variables for the tracking (routine posxy) |
145 |
REAL*8 VECT,VECTINI,VOUT !vectors for Runge-Kutta method |
REAL*8 VECT,VECTINI,VOUT !vectors for Runge-Kutta method |