--- DarthVader/TrackerLevel2/inc/F77/common_momanhough.f 2006/05/19 13:15:53 1.1.1.1 +++ DarthVader/TrackerLevel2/inc/F77/common_momanhough.f 2006/11/21 17:13:52 1.8 @@ -6,6 +6,8 @@ C Author: Vannuccini Elena C Date: 2004/2005 C + +c include 'level1.f' c+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+- c SPEP 1 CLUSTERS ----> COUPLES and SINGLETS c+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+- @@ -16,21 +18,53 @@ c * charge correlation c c------------------------------------------------------------------------ -* cut on maximum number of cluster allowed for track finding -c parameter (nclstrmax_level2=10*nplanes) - parameter (nclstrmax_level2=5*nplanes) -* ----------------------------------------------------------- -* maximum number of couples per plane and total -* (only to dimension the vectors) - parameter (ncouplemax=(nclstrmax_level2/nplanes)**2) - parameter (ncouplemaxtot=ncouplemax*nplanes) + +* maximum (total) number of STORED clusters + parameter (nclstrmax_level2 = nclstrmax) ! same as level1 + +* ----------------------------------------------------------- +* maximum number of cluster (per view) required to perform +* track serching with full Hough transform +* (singlets recovered during track refinement) + parameter (nclusterlimit=6) +* ----------------------------------------------------------- +* maximum number of couples (per plane) required to perform +* track serching with full Hough transform +* (couples recovered during track refinement) + parameter (ncouplelimit=8) + +* ----------------------------------------------------------- +* maximum number of STORED couples, per plane and total +* (to dimension the vectors) +c parameter (ncouplemax=2*(nclstrmax_level2/nplanes)**2) +c parameter (ncouplemaxtot=ncouplemax*nplanes) + parameter (ncouplemax = nclusterlimit*nclusterlimit) + parameter (ncouplemaxtot = ncouplemax*nplanes) * ----------------------------------------------------------- * maximum number of STORED couples * (if the identified couples exceeds this number the * track identification is not performed and the event is * tagged as not good) - parameter (ncp_max=10*nplanes) -c parameter (ncp_max=ncouplemaxtot) +c parameter (ncp_max=10*nplanes) +c parameter (ncp_max=10*nplanes) +c parameter (ncp_max = ncouplemaxtot) + + + +* ----------------------------------------------------------- +* mask of views +* 0 = ok +* 1 = n.clusters > nclusterlimit +* 8 = n.couples > ncouplelimit +* 2 = n.couples > ncouplemax +* 3 = n.doublets > ndblt_max +* 4 = n.triplets > ntrpt_max +* 5 = n.clouds yz > ncloyz_max +* 6 = n.clouds xz > ncloxz_max +* 7 = n.track candidates > NTRACKSMAX + integer mask_view(nviews) + common/maskview/mask_view + * ----------------------------------------------------------- * --- CLUSTERS ---------------------------------------------- @@ -46,9 +80,10 @@ * mask for clusters associated to a track * 1 = used * 0 = not used - integer cl_used(nclstrmax_level2) + integer cl_used(nclstrmax_level2) + integer ncl_view(nviews) !n.clusters per plane * ----------------------------------------------------------- - common/clusters/cl_good,cl_used + common/clusters/cl_good,cl_used,ncl_view * ----------------------------------------------------------- * --- COUPLES ----------------------------------------------- @@ -216,11 +251,13 @@ c c------------------------------------------------------------------------ - +* ----------------------------------------------------- +* default p.f.a. +* used to fit track candidates +* ----------------------------------------------------- character*5 PFAdef -c common/PFAdefault/PFAdef - parameter (PFAdef='COG') - + parameter (PFAdef='COG1') + c------------------------------------------------------------------------ c variable related to track selection and fitting c------------------------------------------------------------------------ @@ -258,5 +295,143 @@ +c------------------------------------------------------------------------ +c CUTS and other parameters +c------------------------------------------------------------------------ + PARAMETER (PIGR=3.14) !159265359)!(1) +* ----------------------------------------------------- +* cuts on cluster signal +* ----------------------------------------------------- + parameter (dedx_x_min=0.) + parameter (dedx_y_min=0.) +* ----------------------------------------------------- +* number of GOOD strips around MAXS +* (NB this cut has been removed in track fitting but +* not in calibration) +* ----------------------------------------------------- + parameter (ngoodstr=2) +* ----------------------------------------------------- +* first cuts on the parameters determined with +* the hough transform, in order to exclude from the +* analysis non physical dublets and triplets +* (out of the geometrical acceptance of the apparatus) +* ----------------------------------------------------- +* geometrical acceptance whit 6 planes +c$$$ parameter (alfyz1_max=8.7) !Y0 (real acc+1cm) +c$$$ parameter (alfyz2_max=0.47) !tg theta-yz (18+7deg) +c$$$ parameter (alfxz1_max=9.5) !X0 (real acc+1cm) +c$$$ parameter (alfxz2_max=0.47) !tg theta-xz (18+7deg) + +c$$$* geometrical acceptance whit the last 5 planes +c$$$ parameter (alfyz1_max=12.4) !Y0 (real acc+1cm) +c$$$ parameter (alfyz2_max=0.55) !tg theta-yz (22+7deg) +c$$$ parameter (alfxz1_max=13.7) !X0 (real acc+1cm) +c$$$ parameter (alfxz2_max=0.55) !tg theta-xz (22+7deg) +c$$$ +c$$$* geometrical acceptance whit the last 4 planes +c$$$ parameter (alfyz1_max=18.6) !Y0 (real acc+1cm) +c$$$ parameter (alfyz2_max=0.70) !tg theta-yz (28+7deg) +c$$$ parameter (alfxz1_max=20.6) !X0 (real acc+1cm) +c$$$ parameter (alfxz2_max=0.70) !tg theta-xz (28+7deg) +c$$$ +* geometrical acceptance whit the last 3 planes + parameter (alfyz1_max=31.0) !Y0 (real acc+1cm) + parameter (alfyz2_max=1.04) !tg theta-yz (39+7deg) + parameter (alfxz1_max=34.4) !X0 (real acc+1cm) + parameter (alfxz2_max=1.04) !tg theta-xz (39+7deg) +* ----------------------------------------------------- +* cut on the position of the circle center: +* center inside the spectrometer volume is not physical +* ==> THE TRIPLET IS REJECTED +* ----------------------------------------------------- + parameter (xclimit=8.) !cm +* ----------------------------------------------------- +* Parameter normalization constants, needed to evaluate +* distances in parameter space +* ----------------------------------------------------- + parameter (Dalfaxz3=1.) + +* --- David ground +c$$$ parameter (Dalfayz1=0.8864e-1) +c$$$ parameter (Dalfayz2=0.6204e-3) +c$$$ parameter (Dalfaxz1=0.2909e-1) +c$$$ parameter (Dalfaxz2=0.2759e-2) +* --- David flight (preliminary) + parameter (Dalfayz1=0.64e-3) + parameter (Dalfayz2=0.54e-3) + parameter (Dalfaxz1=0.50e-3) + parameter (Dalfaxz2=0.12e-3) + + +* ----------------------------------------------------- +* Cut on normalized distances in parameter space. +* Doublets/triplets are recursively included in a cloud +* if the distance from any of the points already included +* is less than this cut. +* ----------------------------------------------------- +* --- David ground +c$$$ parameter(cutystart=0.3) +c$$$ parameter(cutystep=0.3) +c$$$ parameter(cutxstart=1.) +c$$$ parameter(cutxstep=1.) +c$$$ parameter(maxcuty=100.) +c$$$ parameter(maxcutx=150.) +c$$$ parameter(nstepx=50) !inclusion-cut increasing steps +c$$$ parameter(nstepy=50) +* --- David flight (preliminary) + parameter(cutystart=7.) + parameter(cutystep=5.) + parameter(cutxstart=5.) + parameter(cutxstep=2.) + parameter(maxcuty=200.) + parameter(maxcutx=150.) + parameter(nstepx=5) !inclusion-cut increasing steps + parameter(nstepy=5) + + + + real cutdistyz !y0 / tg theta_yz space + real cutdistxz !x0 / tg theta_xz space + common/cutdi/cutdistyz,cutdistxz +* (NB deflection is not considered in the selection criteria) +* -------------------------------------------------- +* cloud selection +* -------------------------------------------------- + parameter(x_min_start=4) + parameter(x_min_step=1) + + integer ncpxz_min + integer nptxz_min !for a circle + integer nplxz_min !for a circle + + common/cutxclouds/ncpxz_min,nptxz_min,nplxz_min + + parameter(y_min_start=3) + + integer ncpyz_min + integer nptyz_min !for a straight line + integer nplyz_min !for a circle + common/cutyclouds/ncpyz_min,nptyz_min,nplyz_min +* -------------------------------------------------- + +* ----------------------------------------------------- +* minimum number of matching couples required to combine +* XZ-YZ clouds and perform the fit +* ----------------------------------------------------- + parameter (ncpok=3) + +* ----------------------------------------------------- +* maximum value of deflection from Hough transform +* accepted in order to perform the track fit +* ----------------------------------------------------- +c parameter (defmax=10000.) !GV-1 + parameter (defmax=100.) !GV-1 + +* ----------------------------------------------------- +* cut to include new couple or single clusters in the +* track fitting, after the first fit +* ----------------------------------------------------- + parameter (clinc=3) +