/[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.4 by pam-fi, Fri Oct 27 16:08:18 2006 UTC revision 1.9 by pam-fi, Fri Apr 27 10:39:57 2007 UTC
# Line 19  c     * charge correlation Line 19  c     * charge correlation
19  c      c    
20  c------------------------------------------------------------------------  c------------------------------------------------------------------------
21    
22  *     maximum (total) number of stored clusters  *     maximum (total) number of STORED clusters
23        parameter (nclstrmax_level2 = nclstrmax)        parameter (nclstrmax_level2 = nclstrmax) ! same as level1
 c      parameter (nclstrmax_level2=10*nplanes)  
 c      parameter (nclstrmax_level2=5*nplanes)  
24    
25  *     -----------------------------------------------------------  *     -----------------------------------------------------------
26  *     maximum number of cluster (per view) required to perform  *     maximum number of cluster (per view) required to perform
27  *     track serching with full Hough transform  *     track serching with full Hough transform
28        parameter (nclustermax=6)  *     (singlets recovered during track refinement)
29          parameter (nclusterlimit=6)
30    *     -----------------------------------------------------------
31    *     maximum number of couples (per plane) required to perform
32    *     track serching with full Hough transform
33    *     (couples recovered during track refinement)
34          parameter (ncouplelimit=8)
35    
36  *     -----------------------------------------------------------  *     -----------------------------------------------------------
37  *     maximum number of couples per plane and total  *     maximum number of STORED couples, per plane and total
38  *     (only to dimension the vectors)  *     (to dimension the vectors)
39  c      parameter (ncouplemax=2*(nclstrmax_level2/nplanes)**2)  c      parameter (ncouplemax=2*(nclstrmax_level2/nplanes)**2)
40  c      parameter (ncouplemaxtot=ncouplemax*nplanes)  c      parameter (ncouplemaxtot=ncouplemax*nplanes)
41         parameter (ncouplemax=nclustermax*nclustermax)         parameter (ncouplemax    = nclusterlimit*nclusterlimit)
42         parameter (ncouplemaxtot=ncouplemax*nplanes)         parameter (ncouplemaxtot = ncouplemax*nplanes)
43  *     -----------------------------------------------------------  *     -----------------------------------------------------------
44  *     maximum number of STORED couples  *     maximum number of STORED couples
45  *     (if the identified couples exceeds this number the  *     (if the identified couples exceeds this number the
# Line 43  c      parameter (ncouplemaxtot=ncouplem Line 47  c      parameter (ncouplemaxtot=ncouplem
47  *     tagged as not good)  *     tagged as not good)
48  c      parameter (ncp_max=10*nplanes)  c      parameter (ncp_max=10*nplanes)
49  c      parameter (ncp_max=10*nplanes)  c      parameter (ncp_max=10*nplanes)
50        parameter (ncp_max=ncouplemaxtot)  c      parameter (ncp_max        = ncouplemaxtot)
51    
52    
53    
54  *     -----------------------------------------------------------  *     -----------------------------------------------------------
55  *     mask of views  *     mask of views
56  *     0 = ok  *     0 = ok
57  *     1 = n.clusters > nclustermax  *     1 = n.clusters > nclusterlimit    
58    *     8 = n.couples  > ncouplelimit
59    *     2 = n.couples  > ncouplemax
60    *     3 = n.doublets > ndblt_max
61    *     4 = n.triplets > ntrpt_max
62    *     5 = n.clouds yz > ncloyz_max
63    *     6 = n.clouds xz > ncloxz_max
64    *     7 = n.track candidates > NTRACKSMAX
65        integer mask_view(nviews)        integer mask_view(nviews)
66        common/maskview/mask_view        common/maskview/mask_view
67    
# Line 269  c--------------------------------------- Line 280  c---------------------------------------
280        REAL YGOOD_STORE(NPLANES,NTRACKSMAX)        REAL YGOOD_STORE(NPLANES,NTRACKSMAX)
281        INTEGER CP_STORE(NPLANES,NTRACKSMAX)        INTEGER CP_STORE(NPLANES,NTRACKSMAX)
282        INTEGER CLS_STORE(NPLANES,NTRACKSMAX)        INTEGER CLS_STORE(NPLANES,NTRACKSMAX)
283          INTEGER SENSOR_STORE(NPLANES,NTRACKSMAX)
284          INTEGER LADDER_STORE(NPLANES,NTRACKSMAX)
285          REAL BX_STORE(NPLANES,NTRACKSMAX)
286          REAL BY_STORE(NPLANES,NTRACKSMAX)
287        REAL RCHI2_STORE(NTRACKSMAX)        REAL RCHI2_STORE(NTRACKSMAX)
288        common/track_candidates/NTRACKS,AL_STORE        common/track_candidates/NTRACKS,AL_STORE
289       $     ,XV_STORE,YV_STORE,ZV_STORE       $     ,XV_STORE,YV_STORE,ZV_STORE
# Line 276  c--------------------------------------- Line 291  c---------------------------------------
291       $     ,RESX_STORE,RESY_STORE       $     ,RESX_STORE,RESY_STORE
292       $     ,AXV_STORE,AYV_STORE       $     ,AXV_STORE,AYV_STORE
293       $     ,XGOOD_STORE,YGOOD_STORE       $     ,XGOOD_STORE,YGOOD_STORE
294       $     ,CP_STORE,CLS_STORE,RCHI2_STORE       $     ,CP_STORE,CLS_STORE,SENSOR_STORE,LADDER_STORE
295         $     ,BX_STORE,BY_STORE
296         $     ,RCHI2_STORE
297    
298    
299  *     best-candidate selection  *     best-candidate selection
# Line 338  c$$$ Line 355  c$$$
355  *     Parameter normalization constants, needed to evaluate  *     Parameter normalization constants, needed to evaluate
356  *     distances in parameter space  *     distances in parameter space
357  *     -----------------------------------------------------  *     -----------------------------------------------------
358        parameter (Dalfaxz3=1.)  
359        parameter (Dalfayz1=0.8864e-1)  *     --- David ground
360        parameter (Dalfayz2=0.6204e-3)  c$$$      parameter (Dalfayz1=0.8864e-1)
361        parameter (Dalfaxz1=0.2909e-1)  c$$$      parameter (Dalfayz2=0.6204e-3)
362        parameter (Dalfaxz2=0.2759e-2)  c$$$      parameter (Dalfaxz1=0.2909e-1)
363    c$$$      parameter (Dalfaxz2=0.2759e-2)
364    *     --- David flight
365          parameter (Dalfayz1=0.64e-3)
366          parameter (Dalfayz2=0.54e-3)
367          parameter (Dalfaxz1=0.50e-3)
368          parameter (Dalfaxz2=0.12e-3)
369          parameter (Dalfaxz3=0.86e-5)
370    
371    
372  *     -----------------------------------------------------  *     -----------------------------------------------------
373  *     Cut on normalized distances in parameter space.  *     Cut on normalized distances in parameter space.
374  *     Doublets/triplets are recursively included in a cloud  *     Doublets/triplets are recursively included in a cloud
375  *     if the distance from any of the points already included  *     if the distance from any of the points already included
376  *     is less than this cut.  *     is less than this cut.
377  *     -----------------------------------------------------  *     -----------------------------------------------------
378        parameter(cutystart=0.3)    *     --- David ground
379        parameter(cutystep=0.3)      c$$$      parameter(cutystart=0.3)  
380        parameter(cutxstart=1.)  c$$$      parameter(cutystep=0.3)    
381        parameter(cutxstep=1.)    c$$$      parameter(cutxstart=1.)
382        parameter(maxcuty=100.)  c$$$      parameter(cutxstep=1.)  
383        parameter(maxcutx=150.)  c$$$      parameter(maxcuty=100.)
384    c$$$      parameter(maxcutx=150.)
385    c$$$      parameter(nstepx=50)     !inclusion-cut increasing steps
386    c$$$      parameter(nstepy=50)
387    *     --- David flight (preliminary)
388    c$$$      parameter(cutystart=7.)  
389    c$$$      parameter(cutystep=5.)              
390    c$$$      parameter(cutxstart=5.)
391    c$$$      parameter(cutxstep=2.)  
392    c$$$      parameter(maxcuty=200.)
393    c$$$      parameter(maxcutx=150.)  
394    c$$$      parameter(nstepx=5)       !inclusion-cut increasing steps
395    c$$$      parameter(nstepy=5)      
396    *     --- David flight
397          parameter(cutystart=30.)  
398          parameter(cutystep=10.)              !buoni???
399          parameter(cutxstart=5.)
400          parameter(cutxstep=2.)  
401          parameter(maxcuty=1000.)
402          parameter(maxcutx=1000.)  
403          parameter(nstepx=10)     !inclusion-cut increasing steps
404          parameter(nstepy=30)
405    
       parameter(nstepx=50)     !inclusion-cut increasing steps  
       parameter(nstepy=50)  
406    
407        real cutdistyz            !y0 / tg theta_yz space          real cutdistyz            !y0 / tg theta_yz space  
408        real cutdistxz            !x0 / tg theta_xz space          real cutdistxz            !x0 / tg theta_xz space  
# Line 393  c$$$ Line 438  c$$$
438  *     maximum value of deflection from Hough transform  *     maximum value of deflection from Hough transform
439  *     accepted in order to perform the track fit  *     accepted in order to perform the track fit
440  *     -----------------------------------------------------  *     -----------------------------------------------------
441        parameter (defmax=10000.) !GV-1  c      parameter (defmax=10000.) !GV-1
442          parameter (defmax=100.) !GV-1
443    
444  *     -----------------------------------------------------  *     -----------------------------------------------------
445  *     cut to include new couple or single clusters in the  *     cut to include new couple or single clusters in the
446  *     track fitting, after the first fit  *     track fitting, after the first fit
447  *     -----------------------------------------------------  *     -----------------------------------------------------
448        parameter (clinc=10)  c      parameter (clinc=3)
449          parameter (clinc=7)
450    
451    
452    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.23