/[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.1 by mocchiut, Fri May 19 13:15:53 2006 UTC revision 1.9 by pam-fi, Fri Apr 27 10:39:57 2007 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  c      parameter (nclstrmax_level2=10*nplanes)  *     maximum (total) number of STORED clusters
23        parameter (nclstrmax_level2=5*nplanes)        parameter (nclstrmax_level2 = nclstrmax) ! same as level1
24  *     -----------------------------------------------------------  
25  *     maximum number of couples per plane and total  *     -----------------------------------------------------------
26  *     (only to dimension the vectors)  *     maximum number of cluster (per view) required to perform
27        parameter (ncouplemax=(nclstrmax_level2/nplanes)**2)  *     track serching with full Hough transform
28        parameter (ncouplemaxtot=ncouplemax*nplanes)  *     (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 STORED couples, per plane and total
38    *     (to dimension the vectors)
39    c      parameter (ncouplemax=2*(nclstrmax_level2/nplanes)**2)
40    c      parameter (ncouplemaxtot=ncouplemax*nplanes)
41           parameter (ncouplemax    = nclusterlimit*nclusterlimit)
42           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
46  *     track identification is not performed and the event is  *     track identification is not performed and the event is
47  *     tagged as not good)  *     tagged as not good)
48        parameter (ncp_max=10*nplanes)  c      parameter (ncp_max=10*nplanes)
49  c      parameter (ncp_max=ncouplemaxtot)  c      parameter (ncp_max=10*nplanes)
50    c      parameter (ncp_max        = ncouplemaxtot)
51    
52    
53    
54    *     -----------------------------------------------------------
55    *     mask of views
56    *     0 = ok
57    *     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)
66          common/maskview/mask_view
67    
68    
69  *     -----------------------------------------------------------  *     -----------------------------------------------------------
70  *     --- CLUSTERS ----------------------------------------------  *     --- CLUSTERS ----------------------------------------------
# Line 46  c      parameter (ncp_max=ncouplemaxtot) Line 80  c      parameter (ncp_max=ncouplemaxtot)
80  *     mask for clusters associated to a track  *     mask for clusters associated to a track
81  *      1 = used  *      1 = used
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
85  *     -----------------------------------------------------------  *     -----------------------------------------------------------
86        common/clusters/cl_good,cl_used        common/clusters/cl_good,cl_used,ncl_view
87    
88  *     -----------------------------------------------------------  *     -----------------------------------------------------------
89  *     --- COUPLES -----------------------------------------------  *     --- COUPLES -----------------------------------------------
# Line 216  c     are referred to the mini notation Line 251  c     are referred to the mini notation
251  c  c
252  c------------------------------------------------------------------------        c------------------------------------------------------------------------      
253    
254    *     -----------------------------------------------------
255    *     default p.f.a.
256    *     used to fit track candidates
257    *     -----------------------------------------------------
258        character*5 PFAdef        character*5 PFAdef
259  c      common/PFAdefault/PFAdef        parameter (PFAdef='COG1')
260        parameter (PFAdef='COG')  
         
261  c------------------------------------------------------------------------  c------------------------------------------------------------------------
262  c     variable related to track selection and fitting  c     variable related to track selection and fitting
263  c------------------------------------------------------------------------  c------------------------------------------------------------------------
# Line 243  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 250  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 258  c--------------------------------------- Line 301  c---------------------------------------
301    
302    
303    
304    c------------------------------------------------------------------------
305    c     CUTS and other parameters
306    c------------------------------------------------------------------------
307          PARAMETER (PIGR=3.14)     !159265359)!(1)
308    *     -----------------------------------------------------
309    *     cuts on cluster signal
310    *     -----------------------------------------------------
311          parameter (dedx_x_min=0.)
312          parameter (dedx_y_min=0.)
313    *     -----------------------------------------------------
314    *     number of GOOD strips around MAXS
315    *     (NB this cut has been removed in track fitting but
316    *     not in calibration)
317    *     -----------------------------------------------------
318          parameter (ngoodstr=2)  
319    *     -----------------------------------------------------
320    *     first cuts on the parameters determined with
321    *     the hough transform, in order to exclude from the
322    *     analysis non physical dublets and triplets
323    *     (out of the geometrical acceptance of the apparatus)
324    *     -----------------------------------------------------
325    *     geometrical acceptance whit 6 planes              
326    c$$$      parameter (alfyz1_max=8.7)  !Y0  (real acc+1cm)    
327    c$$$      parameter (alfyz2_max=0.47) !tg theta-yz (18+7deg)
328    c$$$      parameter (alfxz1_max=9.5)  !X0  (real acc+1cm)    
329    c$$$      parameter (alfxz2_max=0.47) !tg theta-xz (18+7deg)
330    
331    c$$$*     geometrical acceptance whit the last 5 planes      
332    c$$$      parameter (alfyz1_max=12.4) !Y0  (real acc+1cm)  
333    c$$$      parameter (alfyz2_max=0.55) !tg theta-yz (22+7deg)
334    c$$$      parameter (alfxz1_max=13.7) !X0  (real acc+1cm)    
335    c$$$      parameter (alfxz2_max=0.55) !tg theta-xz (22+7deg)
336    c$$$
337    c$$$*     geometrical acceptance whit the last 4 planes      
338    c$$$      parameter (alfyz1_max=18.6) !Y0  (real acc+1cm)  
339    c$$$      parameter (alfyz2_max=0.70) !tg theta-yz (28+7deg)
340    c$$$      parameter (alfxz1_max=20.6) !X0  (real acc+1cm)  
341    c$$$      parameter (alfxz2_max=0.70) !tg theta-xz (28+7deg)
342    c$$$
343    *     geometrical acceptance whit the last 3 planes    
344          parameter (alfyz1_max=31.0) !Y0  (real acc+1cm)  
345          parameter (alfyz2_max=1.04) !tg theta-yz (39+7deg)
346          parameter (alfxz1_max=34.4) !X0  (real acc+1cm)  
347          parameter (alfxz2_max=1.04) !tg theta-xz (39+7deg)
348    *     -----------------------------------------------------
349    *     cut on the position of the circle center:
350    *     center inside the spectrometer volume is not physical
351    *     ==> THE TRIPLET IS REJECTED
352    *     -----------------------------------------------------
353          parameter (xclimit=8.)    !cm
354    *     -----------------------------------------------------
355    *     Parameter normalization constants, needed to evaluate
356    *     distances in parameter space
357    *     -----------------------------------------------------
358    
359    *     --- David ground
360    c$$$      parameter (Dalfayz1=0.8864e-1)
361    c$$$      parameter (Dalfayz2=0.6204e-3)
362    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.
374    *     Doublets/triplets are recursively included in a cloud
375    *     if the distance from any of the points already included
376    *     is less than this cut.
377    *     -----------------------------------------------------
378    *     --- David ground
379    c$$$      parameter(cutystart=0.3)  
380    c$$$      parameter(cutystep=0.3)    
381    c$$$      parameter(cutxstart=1.)
382    c$$$      parameter(cutxstep=1.)  
383    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    
406    
407          real cutdistyz            !y0 / tg theta_yz space  
408          real cutdistxz            !x0 / tg theta_xz space  
409          common/cutdi/cutdistyz,cutdistxz
410    *     (NB deflection is not considered in the selection criteria)
411    *     --------------------------------------------------
412    *     cloud selection
413    *     --------------------------------------------------
414          parameter(x_min_start=4)
415          parameter(x_min_step=1)  
416          
417          integer ncpxz_min
418          integer nptxz_min         !for a circle        
419          integer nplxz_min         !for a circle      
420    
421          common/cutxclouds/ncpxz_min,nptxz_min,nplxz_min
422    
423          parameter(y_min_start=3)
424    
425          integer ncpyz_min      
426          integer nptyz_min         !for a straight line
427          integer nplyz_min         !for a circle        
428          common/cutyclouds/ncpyz_min,nptyz_min,nplyz_min
429    *     --------------------------------------------------
430    
431    *     -----------------------------------------------------
432    *     minimum number of matching couples required to combine
433    *     XZ-YZ clouds and perform the fit
434    *     -----------------------------------------------------
435          parameter (ncpok=3)      
436    
437    *     -----------------------------------------------------
438    *     maximum value of deflection from Hough transform
439    *     accepted in order to perform the track fit
440    *     -----------------------------------------------------
441    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
446    *     track fitting, after the first fit
447    *     -----------------------------------------------------
448    c      parameter (clinc=3)
449          parameter (clinc=7)
450    
451    
452    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.23