/[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.6 by pam-fi, Tue Nov 7 15:55:10 2006 UTC revision 1.12 by pam-fi, Thu Dec 4 14:33:48 2008 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
# Line 50  c      parameter (ncp_max        = ncoup Line 54  c      parameter (ncp_max        = ncoup
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  *     2 = n.couples  > ncouplemax
60  *     3 = n.doublets > ndblt_max  *     3 = n.doublets > ndblt_max
61  *     4 = n.triplets > ntrpt_max  *     4 = n.triplets > ntrpt_max
# Line 77  c      parameter (ncp_max        = ncoup Line 82  c      parameter (ncp_max        = ncoup
82  *      0 = not used  *      0 = not used
83        integer cl_used(nclstrmax_level2)        integer cl_used(nclstrmax_level2)
84        integer ncl_view(nviews)  !n.clusters per plane        integer ncl_view(nviews)  !n.clusters per plane
85          real dedx_x_min,dedx_y_min
86  *     -----------------------------------------------------------  *     -----------------------------------------------------------
87        common/clusters/cl_good,cl_used,ncl_view        common/clusters/cl_good,cl_used,ncl_view,dedx_x_min,dedx_y_min
88    
89  *     -----------------------------------------------------------  *     -----------------------------------------------------------
90  *     --- COUPLES -----------------------------------------------  *     --- COUPLES -----------------------------------------------
# Line 131  c      parameter (ncp_max        = ncoup Line 137  c      parameter (ncp_max        = ncoup
137        common/singlets/ncls,cls,cl_single        common/singlets/ncls,cls,cl_single
138    
139    
140          logical RECOVER_SINGLETS,RECOVER_NUCLEI,SECOND_SEARCH
141          common/recover/RECOVER_SINGLETS,RECOVER_NUCLEI,SECOND_SEARCH
142    
143  c+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +  c+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
144  c + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +  c + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
145  c------------------------------------------------------------------------  c------------------------------------------------------------------------
# Line 258  c     variable related to track selectio Line 267  c     variable related to track selectio
267  c------------------------------------------------------------------------  c------------------------------------------------------------------------
268    
269  *     track candidates  *     track candidates
270        PARAMETER (NTRACKSMAX=ncloyz_max*ncloxz_max)  cc      PARAMETER (NTRACKSMAX=ncloyz_max*ncloxz_max)
271          PARAMETER (NTRACKSMAX=10000)
272        INTEGER NTRACKS           !number of track candidates        INTEGER NTRACKS           !number of track candidates
273        REAL AL_STORE(5,NTRACKSMAX)        REAL AL_STORE(5,NTRACKSMAX)
274        REAL XM_STORE(NPLANES,NTRACKSMAX)        REAL XM_STORE(NPLANES,NTRACKSMAX)
# Line 275  c--------------------------------------- Line 285  c---------------------------------------
285        REAL YGOOD_STORE(NPLANES,NTRACKSMAX)        REAL YGOOD_STORE(NPLANES,NTRACKSMAX)
286        INTEGER CP_STORE(NPLANES,NTRACKSMAX)        INTEGER CP_STORE(NPLANES,NTRACKSMAX)
287        INTEGER CLS_STORE(NPLANES,NTRACKSMAX)        INTEGER CLS_STORE(NPLANES,NTRACKSMAX)
288          INTEGER SENSOR_STORE(NPLANES,NTRACKSMAX)
289          INTEGER LADDER_STORE(NPLANES,NTRACKSMAX)
290          REAL BX_STORE(NPLANES,NTRACKSMAX)
291          REAL BY_STORE(NPLANES,NTRACKSMAX)
292        REAL RCHI2_STORE(NTRACKSMAX)        REAL RCHI2_STORE(NTRACKSMAX)
293        common/track_candidates/NTRACKS,AL_STORE        common/trackcandidates/NTRACKS,AL_STORE
294       $     ,XV_STORE,YV_STORE,ZV_STORE       $     ,XV_STORE,YV_STORE,ZV_STORE
295       $     ,XM_STORE,YM_STORE,ZM_STORE       $     ,XM_STORE,YM_STORE,ZM_STORE
296       $     ,RESX_STORE,RESY_STORE       $     ,RESX_STORE,RESY_STORE
297       $     ,AXV_STORE,AYV_STORE       $     ,AXV_STORE,AYV_STORE
298       $     ,XGOOD_STORE,YGOOD_STORE       $     ,XGOOD_STORE,YGOOD_STORE
299       $     ,CP_STORE,CLS_STORE,RCHI2_STORE       $     ,CP_STORE,CLS_STORE
300         $     ,SENSOR_STORE,LADDER_STORE
301         $     ,BX_STORE,BY_STORE
302         $     ,RCHI2_STORE
303    
304    
305  *     best-candidate selection  *     best-candidate selection
# Line 295  c     CUTS and other parameters Line 312  c     CUTS and other parameters
312  c------------------------------------------------------------------------  c------------------------------------------------------------------------
313        PARAMETER (PIGR=3.14)     !159265359)!(1)        PARAMETER (PIGR=3.14)     !159265359)!(1)
314  *     -----------------------------------------------------  *     -----------------------------------------------------
315  *     cuts on cluster signal  *     cuts on cluster signal (MIP)
316  *     -----------------------------------------------------  *     -----------------------------------------------------
317        parameter (dedx_x_min=0.)        parameter (dedx_x_min_mip=0.)
318        parameter (dedx_y_min=0.)        parameter (dedx_y_min_mip=0.)
319          parameter (dedx_x_min_nuclei=5.)
320          parameter (dedx_y_min_nuclei=5.)
321          parameter (dedx_min_nuclei=9.1) !cut to apply nuclei recovery
322          parameter (ddedx_min_nuclei=0.35)!cut to identify bad tracked nuclei
323    
324  *     -----------------------------------------------------  *     -----------------------------------------------------
325  *     number of GOOD strips around MAXS  *     number of GOOD strips around MAXS
326  *     (NB this cut has been removed in track fitting but  *     (NB this cut has been removed in track fitting but
# Line 344  c$$$ Line 366  c$$$
366  *     Parameter normalization constants, needed to evaluate  *     Parameter normalization constants, needed to evaluate
367  *     distances in parameter space  *     distances in parameter space
368  *     -----------------------------------------------------  *     -----------------------------------------------------
369        parameter (Dalfaxz3=1.)  
370        parameter (Dalfayz1=0.8864e-1)  *     --- David ground
371        parameter (Dalfayz2=0.6204e-3)  c$$$      parameter (Dalfayz1=0.8864e-1)
372        parameter (Dalfaxz1=0.2909e-1)  c$$$      parameter (Dalfayz2=0.6204e-3)
373        parameter (Dalfaxz2=0.2759e-2)  c$$$      parameter (Dalfaxz1=0.2909e-1)
374    c$$$      parameter (Dalfaxz2=0.2759e-2)
375    *     --- David flight
376          parameter (Dalfayz1=0.64e-3)
377          parameter (Dalfayz2=0.54e-3)
378          parameter (Dalfaxz1=0.50e-3)
379          parameter (Dalfaxz2=0.12e-3)
380          parameter (Dalfaxz3=0.86e-5)
381    
382    
383  *     -----------------------------------------------------  *     -----------------------------------------------------
384  *     Cut on normalized distances in parameter space.  *     Cut on normalized distances in parameter space.
385  *     Doublets/triplets are recursively included in a cloud  *     Doublets/triplets are recursively included in a cloud
386  *     if the distance from any of the points already included  *     if the distance from any of the points already included
387  *     is less than this cut.  *     is less than this cut.
388  *     -----------------------------------------------------  *     -----------------------------------------------------
389        parameter(cutystart=0.3)    *     --- David ground
390        parameter(cutystep=0.3)      c$$$      parameter(cutystart=0.3)  
391        parameter(cutxstart=1.)  c$$$      parameter(cutystep=0.3)    
392        parameter(cutxstep=1.)    c$$$      parameter(cutxstart=1.)
393        parameter(maxcuty=100.)  c$$$      parameter(cutxstep=1.)  
394        parameter(maxcutx=150.)  c$$$      parameter(maxcuty=100.)
395    c$$$      parameter(maxcutx=150.)
396    c$$$      parameter(nstepx=50)     !inclusion-cut increasing steps
397    c$$$      parameter(nstepy=50)
398    *     --- David flight (preliminary)
399    c$$$      parameter(cutystart=7.)  
400    c$$$      parameter(cutystep=5.)              
401    c$$$      parameter(cutxstart=5.)
402    c$$$      parameter(cutxstep=2.)  
403    c$$$      parameter(maxcuty=200.)
404    c$$$      parameter(maxcutx=150.)  
405    c$$$      parameter(nstepx=5)       !inclusion-cut increasing steps
406    c$$$      parameter(nstepy=5)      
407    *     --- David flight
408          parameter(cutystart=30.)  
409          parameter(cutystep=10.)              !buoni???
410          parameter(cutxstart=5.)
411          parameter(cutxstep=2.)  
412          parameter(maxcuty=1000.)
413          parameter(maxcutx=1000.)  
414          parameter(nstepx=10)     !inclusion-cut increasing steps
415          parameter(nstepy=30)
416    
       parameter(nstepx=50)     !inclusion-cut increasing steps  
       parameter(nstepy=50)  
417    
418        real cutdistyz            !y0 / tg theta_yz space          real cutdistyz            !y0 / tg theta_yz space  
419        real cutdistxz            !x0 / tg theta_xz space          real cutdistxz            !x0 / tg theta_xz space  
# Line 406  c      parameter (defmax=10000.) !GV-1 Line 456  c      parameter (defmax=10000.) !GV-1
456  *     cut to include new couple or single clusters in the  *     cut to include new couple or single clusters in the
457  *     track fitting, after the first fit  *     track fitting, after the first fit
458  *     -----------------------------------------------------  *     -----------------------------------------------------
459        parameter (clinc=3)  c      parameter (clinc=3)
460          parameter (clinc=7)
461    
462    
463    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.23