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

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

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

revision 1.2 by pam-fi, Thu Sep 28 14:04:39 2006 UTC revision 1.11 by pam-fi, Wed Aug 22 07:03:45 2007 UTC
# Line 7  Line 7 
7  *      *    
8  *************************************************************************  *************************************************************************
9    
10          parameter (clcutx=7.)     !cluster seed cut
11          parameter (clcuty=6.)
12          parameter (incutx=4.)     !cut to include strips in cluster definition
13          parameter (incuty=4.)
14    
15  c------------------------------------------------------------------------  c------------------------------------------------------------------------
16  c      c    
17  c     pedestal, sigma and badstrip variables  c     some parameters to configure data reduction
18  c      c    
19  c------------------------------------------------------------------------  c------------------------------------------------------------------------
20    c$$$      real clcutx               !cluster seed cut
21    c$$$      real clcuty
22    c$$$      real incutx               !cut to include strips in cluster definition
23    c$$$      real incuty
24          integer pfaid             !if of PFA
25          common/sw/
26         $     pfaid
27    c$$$     $     ,clcutx,clcuty
28    c$$$     $     ,incutx,incuty
29    
30          save/sw/
31    
32        parameter (id_hi_bad=100) !badstrip histo identifier  c------------------------------------------------------------------------
33        parameter (id_hi_ped=200) !pedestal histo identifier  c    
34        parameter (id_hi_sig=300) !sigma histo identifier  c     pedestal, sigma and badstrip variables
35    c    
36    c------------------------------------------------------------------------
37    
38        real pedestal(nviews,nva1_view,nstrips_va1) !pedestal value        real pedestal(nviews,nva1_view,nstrips_va1) !pedestal value
39        real pedestal_t(nviews,nva1_view,nstrips_va1) !pedestal truncated value        real pedestal_t(nviews,nva1_view,nstrips_va1) !pedestal truncated value
# Line 32  c--------------------------------------- Line 50  c---------------------------------------
50        save/pedsigbad/        save/pedsigbad/
51    
52    
       parameter (clcutx=7.)     !cluster seed cut  
       parameter (clcuty=7.)  
       parameter (incutx=4.)     !cut to include strips in cluster definition  
       parameter (incuty=4.)  
53    
54    c------------------------------------------------------------------------
55    c    
56    c     VA1-mask
57    c    
58    c------------------------------------------------------------------------
59    c     VA1 mask, from the DB
60    c      0 = force mask
61    c      1 = ---> run-by-run mask ---> event-by-event mask
62    c     -1 = ---> event-by-event mask
63    c     -------------------------------------
64          integer mask_vk(nviews,nva1_view)    
65    c     -------------------------------------
66    c     VA1 mask run-by-run (set to 0 if <SIG> < SIG_min)
67    c      0 = force mask
68    c      1 = ---> event-by-event mask
69    c     -------------------------------------
70          integer mask_vk_run(nviews,nva1_view)
71    c     -------------------------------------
72    c     VA1 mask event-by-event (set to 0 if CN computation fails)
73    c      0 = force mask
74    c      1 = search clusters
75    c     -------------------------------------
76          integer mask_vk_ev(nviews,nva1_view)
77    
78        parameter (id_hi_mask=400) !mask histo identifier        integer mask(nviews,nva1_view,nstrips_va1) !effective mask
   
       integer mask(nviews,nva1_view,nstrips_va1)  
       integer mask_vk(nviews,nva1_view)  
79    
80        common/mask/mask,mask_vk  *     set from outside F77 routines
81          common/mask/mask_vk,mask_vk_run
82        save/mask/        save/mask/
83    *     set event-by-event
84          common/maskev/mask_vk_ev,mask
85    
86  c------------------------------------------------------------------------  c------------------------------------------------------------------------
87  c      c    
# Line 52  c     eta2, eta3, eta4 pfa correction pa Line 89  c     eta2, eta3, eta4 pfa correction pa
89  c      c    
90  c------------------------------------------------------------------------  c------------------------------------------------------------------------
91  *     angular binning  *     angular binning
 *     NB the angular binning is defined in common_preanalysis.f  
 *     here are defined variables to retrieve calibration info  
 *     (for processing levels next to 1)  
92        parameter (nangmax=21)    !maximum number of angular bins        parameter (nangmax=21)    !maximum number of angular bins
93        integer nangbin        integer nangbin
94        real angL(nangmax),angR(nangmax) !Left and Right bin limits        real angL(nangmax),angR(nangmax) !Left and Right bin limits
95    
96        parameter (netavalmax=150)   !eta2 points        parameter (netavalmax=500)   !eta2 points
97        integer netaval        integer netaval
98        real eta2(netavalmax,nangmax)        real eta2(netavalmax,nangmax)
99        real feta2(netavalmax,nviews,nladders_view,nangmax)        real feta2(netavalmax,nviews,nladders_view,nangmax)
# Line 67  c--------------------------------------- Line 101  c---------------------------------------
101        real feta3(netavalmax,nviews,nladders_view,nangmax)        real feta3(netavalmax,nviews,nladders_view,nangmax)
102        real eta4(netavalmax,nangmax)        real eta4(netavalmax,nangmax)
103        real feta4(netavalmax,nviews,nladders_view,nangmax)        real feta4(netavalmax,nviews,nladders_view,nangmax)
104          real fcorr(nviews,nladders_view,nangmax)
105        common/angbinning/nangbin,angL,angR  *     limit of application of pfa
106        save/angbinning/        real e2fax,e2tax,e3fax,e3tax,e4fax,e4tax
107        common/pfa/netaval,eta2,feta2,eta3,feta3,eta4,feta4        real e2fay,e2tay,e3fay,e3tay,e4fay,e4tay
108          
109    
110    c$$$      common/angbinning/nangbin,angL,angR
111    c$$$      save/angbinning/
112    c$$$      common/pfa/netaval,eta2,feta2,eta3,feta3,eta4,feta4
113    c$$$      save/pfa/
114        
115          common/pfa/nangbin,angL,angR
116         $     ,netaval,eta2,feta2,eta3,feta3,eta4,feta4
117         $     ,fcorr
118         $     ,e2fax,e2tax,e3fax,e3tax,e4fax,e4tax
119         $     ,e2fay,e2tay,e3fay,e3tay,e4fay,e4tay
120        save/pfa/        save/pfa/
121    
122    c$$$*     limit of application of pfa
123    c$$$      parameter (e2fax=0.)
124    c$$$      parameter (e2tax=10.)
125    c$$$      parameter (e3fax=10.)
126    c$$$      parameter (e3tax=15.)
127    c$$$      parameter (e4fax=90.)
128    c$$$      parameter (e4tax=90.)
129    c$$$
130    c$$$      parameter (e2fay=0.)
131    c$$$      parameter (e2tay=20.)
132    c$$$      parameter (e3fay=90.)
133    c$$$      parameter (e3tay=90.)
134    c$$$      parameter (e4fay=90.)
135    c$$$      parameter (e4tay=90.)
136    
137  c------------------------------------------------------------------------  c------------------------------------------------------------------------
138  c      c    
139  c     CHARGE CORRELATION  c     CHARGE CORRELATION
# Line 80  c--------------------------------------- Line 142  c---------------------------------------
142    
143  c      parameter (chcut=3.)      !cut to associate two clusters  c      parameter (chcut=3.)      !cut to associate two clusters
144  c      parameter (chcut=4.)      !cut to associate two clusters  c      parameter (chcut=4.)      !cut to associate two clusters
145        parameter (chcut=10.)      !cut to associate two clusters        parameter (chcut=50.)      !cut to associate two clusters
146        parameter (chsatx=2500.)   !saturation limit        parameter (chsatx=1700.)   !saturation limit
147        parameter (chsaty=2500.)   !saturation limit        parameter (chsaty=2500.)   !saturation limit
148          parameter (chmipx=200.)   !mip limit
149          parameter (chmipy=200.)   !mip limit
150  *     above saturation limit correlation conditions are not applied  *     above saturation limit correlation conditions are not applied
151        real kch(nplanes,nladders_view) !angular coeff        real kch(nplanes,nladders_view) !angular coeff
152        real cch(nplanes,nladders_view) !const        real cch(nplanes,nladders_view) !const
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.23