/[PAMELA software]/DarthVader/TrackerLevel2/inc/F77/common_mini_2.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/inc/F77/common_mini_2.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by mocchiut, Fri May 19 13:15:53 2006 UTC revision 1.10 by pam-fi, Tue Jan 15 14:27:24 2008 UTC
# Line 40  c      PARAMETER (PIGR=3.141592654)!???? Line 40  c      PARAMETER (PIGR=3.141592654)!????
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
       REAL*8 XGOOD,YGOOD  
49        REAL*8 RESX,RESY        REAL*8 RESX,RESY
50  c     REAL*8 DEDXTRK            !(1)        REAL*8 TAILX,TAILY
51        REAL*8 DEDXTRK_X          !(1)        REAL*8 XGOOD,YGOOD
52        REAL*8 DEDXTRK_Y          !(1)        REAL*8 DEDXTRK_X          
53          REAL*8 DEDXTRK_Y        
54          REAL*8 ZINI     ! ***PP***
55          REAL*8 PFIXED   ! ***PP***
56          REAL*8 CHI2
57          REAL*8 XV,YV,ZV
58          REAL*8 AXV,AYV
59          REAL*8 COV
60          REAL*8 FACT
61          INTEGER NSTEP
62          INTEGER IDCAND            
63          INTEGER TRACKMODE
64          INTEGER ISTEPMIN
65        COMMON/TRACK/        COMMON/TRACK/
66       $      AL(5)                                      !track parameters       $      AL(5)                                      !track parameters
67       $     ,XM(NPLANES),YM(NPLANES),ZM(NPLANES)        !couple    => xyz point       $     ,XM(NPLANES),YM(NPLANES),ZM(NPLANES)        !couple    => xyz point
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
76         $     ,PFIXED
77         $     ,CHI2
78         $     ,XV(NPLANES),YV(NPLANES),ZV(NPLANES)
79         $     ,AXV(NPLANES),AYV(NPLANES)
80         $     ,COV(5,5)
81         $     ,FACT                ! factor scale for tracking precision
82         $     ,NSTEP
83         $     ,IDCAND              !(2) track-candidate ID (to relate some level1 info)
84         $     ,TRACKMODE           ! tracking mode, 0=normal, 1=with bremsstrahlung
85         $     ,ISTEPMIN            ! minimum number of step during tracking
86  c******************************************************  c******************************************************
87  cccccc 06/10/2005 modified by elena vannuccini ---> (1)  cccccc 06/10/2005 modified by elena vannuccini ---> (1)
88    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
# Line 70  c*************************************** Line 97  c***************************************
97  *     - projected angles  *     - projected angles
98  *     - parameter covariance matrix  *     - parameter covariance matrix
99  *     -------------------------------------------------------------------------  *     -------------------------------------------------------------------------
100        REAL*8 CHI2  c$$$      REAL*8 ZINI     ! ***PP***
101        REAL*8 XV,YV,ZV  c$$$      REAL*8 PFIXED   ! ***PP***
102        REAL*8 AXV,AYV  c$$$      REAL*8 CHI2
103        REAL*8 COV  c$$$      REAL*8 XV,YV,ZV
104        COMMON/FIT_INFO/  c$$$      REAL*8 AXV,AYV
105       $      CHI2  c$$$      REAL*8 COV
106       $     ,XV(NPLANES),YV(NPLANES),ZV(NPLANES)  c$$$      COMMON/FIT_INFO/
107       $     ,AXV(NPLANES),AYV(NPLANES)  c$$$     $      ZINI
108       $     ,COV(5,5)              c$$$     $     ,PFIXED
109    c$$$     $     ,CHI2
110    c$$$     $     ,XV(NPLANES),YV(NPLANES),ZV(NPLANES)
111    c$$$     $     ,AXV(NPLANES),AYV(NPLANES)
112    c$$$     $     ,COV(5,5)            
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        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       $     ,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
# Line 106  c*************************************** Line 138  c***************************************
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
# Line 115  c*************************************** Line 149  c***************************************
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,DLT        ! DELTA-DEF = DELTA0 + DELTA1 x Y
155          COMMON/DELTAB/DELTA0,DELTA1,DLT
156        
157            
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.23