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