| 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 |
| 149 |
COMMON/TRACKING/VECT(7),VECTINI(7),VOUT(7),TOLL,STEPMAX, |
COMMON/TRACKING/VECT(7),VECTINI(7),VOUT(7),TOLL,STEPMAX, |
| 150 |
$ CHARGE |
$ CHARGE |
| 151 |
|
|
| 152 |
|
* ------------------------------------------------------------------------- |
| 153 |
|
* variables for the DeltaB coorection |
| 154 |
|
REAL*8 DELTA0,DELTA1 ! DELTA-DEF = DELTA0 + DELTA1 x Y |
| 155 |
|
COMMON/DELTAB/DELTA0,DELTA1 |
| 156 |
|
|
| 157 |
|
|