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

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

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

revision 1.2 by pam-fi, Mon Oct 2 13:12:47 2006 UTC revision 1.4 by pam-fi, Fri Oct 27 16:08:18 2006 UTC
# Line 6  C Line 6  C
6  C     Author: Vannuccini Elena  C     Author: Vannuccini Elena
7  C     Date: 2004/2005  C     Date: 2004/2005
8  C  C
9    
10    c      include 'level1.f'
11  c+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-  c+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-
12  c     SPEP 1                          CLUSTERS ----> COUPLES and SINGLETS  c     SPEP 1                          CLUSTERS ----> COUPLES and SINGLETS
13  c+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-  c+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-
# Line 16  c     * geometrical criteria Line 18  c     * geometrical criteria
18  c     * charge correlation  c     * charge correlation
19  c      c    
20  c------------------------------------------------------------------------  c------------------------------------------------------------------------
21  *     cut on maximum number of cluster allowed for track finding  
22        parameter (nclstrmax_level2=10*nplanes)  *     maximum (total) number of stored clusters
23          parameter (nclstrmax_level2 = nclstrmax)
24    c      parameter (nclstrmax_level2=10*nplanes)
25  c      parameter (nclstrmax_level2=5*nplanes)  c      parameter (nclstrmax_level2=5*nplanes)
26    
27    *     -----------------------------------------------------------
28    *     maximum number of cluster (per view) required to perform
29    *     track serching with full Hough transform
30          parameter (nclustermax=6)
31    
32  *     -----------------------------------------------------------  *     -----------------------------------------------------------
33  *     maximum number of couples per plane and total  *     maximum number of couples per plane and total
34  *     (only to dimension the vectors)  *     (only to dimension the vectors)
35        parameter (ncouplemax=(nclstrmax_level2/nplanes)**2)  c      parameter (ncouplemax=2*(nclstrmax_level2/nplanes)**2)
36        parameter (ncouplemaxtot=ncouplemax*nplanes)  c      parameter (ncouplemaxtot=ncouplemax*nplanes)
37           parameter (ncouplemax=nclustermax*nclustermax)
38           parameter (ncouplemaxtot=ncouplemax*nplanes)
39  *     -----------------------------------------------------------  *     -----------------------------------------------------------
40  *     maximum number of STORED couples  *     maximum number of STORED couples
41  *     (if the identified couples exceeds this number the  *     (if the identified couples exceeds this number the
42  *     track identification is not performed and the event is  *     track identification is not performed and the event is
43  *     tagged as not good)  *     tagged as not good)
44  c      parameter (ncp_max=10*nplanes)  c      parameter (ncp_max=10*nplanes)
45        parameter (ncp_max=20*nplanes)  c      parameter (ncp_max=10*nplanes)
46  c      parameter (ncp_max=ncouplemaxtot)        parameter (ncp_max=ncouplemaxtot)
47    
48    
49    
50    *     -----------------------------------------------------------
51    *     mask of views
52    *     0 = ok
53    *     1 = n.clusters > nclustermax
54          integer mask_view(nviews)
55          common/maskview/mask_view
56    
57    
58  *     -----------------------------------------------------------  *     -----------------------------------------------------------
59  *     --- CLUSTERS ----------------------------------------------  *     --- CLUSTERS ----------------------------------------------
# Line 47  c      parameter (ncp_max=ncouplemaxtot) Line 69  c      parameter (ncp_max=ncouplemaxtot)
69  *     mask for clusters associated to a track  *     mask for clusters associated to a track
70  *      1 = used  *      1 = used
71  *      0 = not used  *      0 = not used
72        integer cl_used(nclstrmax_level2)        integer cl_used(nclstrmax_level2)
73          integer ncl_view(nviews)  !n.clusters per plane
74  *     -----------------------------------------------------------  *     -----------------------------------------------------------
75        common/clusters/cl_good,cl_used        common/clusters/cl_good,cl_used,ncl_view
76    
77  *     -----------------------------------------------------------  *     -----------------------------------------------------------
78  *     --- COUPLES -----------------------------------------------  *     --- COUPLES -----------------------------------------------
# Line 217  c     are referred to the mini notation Line 240  c     are referred to the mini notation
240  c  c
241  c------------------------------------------------------------------------        c------------------------------------------------------------------------      
242    
243    *     -----------------------------------------------------
244    *     default p.f.a.
245    *     used to fit track candidates
246    *     -----------------------------------------------------
247        character*5 PFAdef        character*5 PFAdef
248  c      common/PFAdefault/PFAdef        parameter (PFAdef='COG1')
249        parameter (PFAdef='COG')  
         
250  c------------------------------------------------------------------------  c------------------------------------------------------------------------
251  c     variable related to track selection and fitting  c     variable related to track selection and fitting
252  c------------------------------------------------------------------------  c------------------------------------------------------------------------
# Line 259  c--------------------------------------- Line 284  c---------------------------------------
284    
285    
286    
287    c------------------------------------------------------------------------
288    c     CUTS and other parameters
289    c------------------------------------------------------------------------
290          PARAMETER (PIGR=3.14)     !159265359)!(1)
291    *     -----------------------------------------------------
292    *     cuts on cluster signal
293    *     -----------------------------------------------------
294          parameter (dedx_x_min=0.)
295          parameter (dedx_y_min=0.)
296    *     -----------------------------------------------------
297    *     number of GOOD strips around MAXS
298    *     (NB this cut has been removed in track fitting but
299    *     not in calibration)
300    *     -----------------------------------------------------
301          parameter (ngoodstr=2)  
302    *     -----------------------------------------------------
303    *     first cuts on the parameters determined with
304    *     the hough transform, in order to exclude from the
305    *     analysis non physical dublets and triplets
306    *     (out of the geometrical acceptance of the apparatus)
307    *     -----------------------------------------------------
308    *     geometrical acceptance whit 6 planes              
309    c$$$      parameter (alfyz1_max=8.7)  !Y0  (real acc+1cm)    
310    c$$$      parameter (alfyz2_max=0.47) !tg theta-yz (18+7deg)
311    c$$$      parameter (alfxz1_max=9.5)  !X0  (real acc+1cm)    
312    c$$$      parameter (alfxz2_max=0.47) !tg theta-xz (18+7deg)
313    
314    c$$$*     geometrical acceptance whit the last 5 planes      
315    c$$$      parameter (alfyz1_max=12.4) !Y0  (real acc+1cm)  
316    c$$$      parameter (alfyz2_max=0.55) !tg theta-yz (22+7deg)
317    c$$$      parameter (alfxz1_max=13.7) !X0  (real acc+1cm)    
318    c$$$      parameter (alfxz2_max=0.55) !tg theta-xz (22+7deg)
319    c$$$
320    c$$$*     geometrical acceptance whit the last 4 planes      
321    c$$$      parameter (alfyz1_max=18.6) !Y0  (real acc+1cm)  
322    c$$$      parameter (alfyz2_max=0.70) !tg theta-yz (28+7deg)
323    c$$$      parameter (alfxz1_max=20.6) !X0  (real acc+1cm)  
324    c$$$      parameter (alfxz2_max=0.70) !tg theta-xz (28+7deg)
325    c$$$
326    *     geometrical acceptance whit the last 3 planes    
327          parameter (alfyz1_max=31.0) !Y0  (real acc+1cm)  
328          parameter (alfyz2_max=1.04) !tg theta-yz (39+7deg)
329          parameter (alfxz1_max=34.4) !X0  (real acc+1cm)  
330          parameter (alfxz2_max=1.04) !tg theta-xz (39+7deg)
331    *     -----------------------------------------------------
332    *     cut on the position of the circle center:
333    *     center inside the spectrometer volume is not physical
334    *     ==> THE TRIPLET IS REJECTED
335    *     -----------------------------------------------------
336          parameter (xclimit=8.)    !cm
337    *     -----------------------------------------------------
338    *     Parameter normalization constants, needed to evaluate
339    *     distances in parameter space
340    *     -----------------------------------------------------
341          parameter (Dalfaxz3=1.)
342          parameter (Dalfayz1=0.8864e-1)
343          parameter (Dalfayz2=0.6204e-3)
344          parameter (Dalfaxz1=0.2909e-1)
345          parameter (Dalfaxz2=0.2759e-2)
346    *     -----------------------------------------------------
347    *     Cut on normalized distances in parameter space.
348    *     Doublets/triplets are recursively included in a cloud
349    *     if the distance from any of the points already included
350    *     is less than this cut.
351    *     -----------------------------------------------------
352          parameter(cutystart=0.3)  
353          parameter(cutystep=0.3)    
354          parameter(cutxstart=1.)
355          parameter(cutxstep=1.)  
356          parameter(maxcuty=100.)
357          parameter(maxcutx=150.)
358    
359          parameter(nstepx=50)     !inclusion-cut increasing steps
360          parameter(nstepy=50)
361    
362          real cutdistyz            !y0 / tg theta_yz space  
363          real cutdistxz            !x0 / tg theta_xz space  
364          common/cutdi/cutdistyz,cutdistxz
365    *     (NB deflection is not considered in the selection criteria)
366    *     --------------------------------------------------
367    *     cloud selection
368    *     --------------------------------------------------
369          parameter(x_min_start=4)
370          parameter(x_min_step=1)  
371          
372          integer ncpxz_min
373          integer nptxz_min         !for a circle        
374          integer nplxz_min         !for a circle      
375    
376          common/cutxclouds/ncpxz_min,nptxz_min,nplxz_min
377    
378          parameter(y_min_start=3)
379    
380          integer ncpyz_min      
381          integer nptyz_min         !for a straight line
382          integer nplyz_min         !for a circle        
383          common/cutyclouds/ncpyz_min,nptyz_min,nplyz_min
384    *     --------------------------------------------------
385    
386    *     -----------------------------------------------------
387    *     minimum number of matching couples required to combine
388    *     XZ-YZ clouds and perform the fit
389    *     -----------------------------------------------------
390          parameter (ncpok=3)      
391    
392    *     -----------------------------------------------------
393    *     maximum value of deflection from Hough transform
394    *     accepted in order to perform the track fit
395    *     -----------------------------------------------------
396          parameter (defmax=10000.) !GV-1
397    
398    *     -----------------------------------------------------
399    *     cut to include new couple or single clusters in the
400    *     track fitting, after the first fit
401    *     -----------------------------------------------------
402          parameter (clinc=10)
403    
404    
405    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.23