# Define the path for CFT-only track finding.
#
# D.Adams.
# 18jun99 1240
#
# It is assumed that the following layers are already defined:
# CFTXi i = 1, 2, ..., 8
# CFTUi where i = 1, 2, 3, 4
# CFTVi where i = 1, 2, 3, 4

# Uncomment the following line to enable debugging messages.
# debug = true

# Local definitions.
bfield = 2.0
ptmin = 0.45
trkvec = double( 0.0 0.0 0.0 0.0 0.0 )
trkerr = double(
   20.0
    0.0 100.0
    0.0   0.0  20.1
    0.0   0.0   0.0  10.0
    0.0   0.0   0.0   0.0  20.0
  )

# Starting track surface.
[ srft SurfCylinder radius=60.0 ]

# Propagator.
[ propc PropCyl bfield=$bfield ]

# Simple fitters
[ fitk AddFitKalman ]
[ fitp AddFitCyl_Phi ]
[ fitpp AddFitCyl_Phi_Phi
  bfield=$bfield
  ptmin=$ptmin
  fitter=*fitk
  propagator=@0
]
[ fitppp AddFitCyl_Phi_Phi_Phi nsigma=5.0 fitter=*fitk ]
[ fits AddFitCyl_PhiZ zmin=-1000. zmax=+1000. ]
[ fitss AddFitCyl_PhiZ_PhiZ
  bfield = $bfield
  fitter = *fitk
  propagator = *propc
  surface = @srft
  error_matrix = $trkerr
  flag = "DR_ALPHA"
  z0max = 200.0
]
  
# Combined fitter.
[ fit_start AddFitStarter max_hits=5
  default_fitter=@fitk
  hit_types = string(
    "HitCylPhi" ":"
    "HitCylPhi" "HitCylPhi" ":"
    "HitCylPhi" "HitCylPhi" "HitCylPhi" ":"
    "HitCylPhi" "HitCylPhiZ" ":"
    "HitCylPhi" "HitCylPhi" "HitCylPhiZ" ":"
    "HitCylPhi" "HitCylPhi" "HitCylPhi" "HitCylPhiZ" ":"
    "HitCylPhi" "HitCylPhi" "HitCylPhi" "HitCylPhi" "HitCylPhiZ" ":"

    "HitCylPhi" "HitCylPhiZ" "HitCylPhiZ" ":"
    "HitCylPhiZ" "HitCylPhi" "HitCylPhiZ" ":"
    "HitCylPhi" "HitCylPhiZ" "HitCylPhi" "HitCylPhiZ" ":"
    "HitCylPhi" "HitCylPhi" "HitCylPhi" "HitCylPhiZ" "HitCylPhiZ" ":"
  )
  fitters = @( fitp fitpp fitppp fits fits fits fits
               fitss fitss fitss fitss )
]

# Checkers.
[ chkcs06 CheckChiSquare max_chisq=10.0 ]
[ chkcs07 CheckChiSquare max_chisq=15.0 ]
[ chkcs08 CheckChiSquare max_chisq=18.0 ]
[ chkcs09 CheckChiSquare max_chisq=20.0 ]
[ chkcs10 CheckChiSquare max_chisq=22.0 ]
[ chkcs11 CheckChiSquare max_chisq=24.0 ]
[ chkcs12 CheckChiSquare max_chisq=26.0 ]
[ chkcs13 CheckChiSquare max_chisq=28.0 ]
[ chkcs14 CheckChiSquare max_chisq=30.0 ]
[ chkcs15 CheckChiSquare max_chisq=32.0 ]
[ chkcs16 CheckChiSquare max_chisq=34.0 ]

# Filters.
[ filter_shares FilterShares min_share_count=5 ]

# Stops.
[ stop_end PathStop filters=@( filter_shares ) ]

# Paths.

[ head Path parent=@0 checkers=@( ) stop=@0 end=false ]

[ pcftx8 Path parent=@head layer=*CFTX8 propagator=*propc fitter=*fit_start
  checkers=@( ) stop=@0 end=false ]

[ pcftx7 Path parent=@pcftx8 layer=*CFTX7 propagator=*propc fitter=*fit_start
  checkers=@( ) stop=@0 end=false ]

[ pcftx6 Path parent=@pcftx7 layer=*CFTX6 propagator=*propc fitter=*fit_start
  checkers=@( ) stop=@0 end=false ]

[ pcftv4 Path parent=@pcftx6 layer=*CFTV4 propagator=*propc fitter=*fit_start
  checkers=@( ) stop=@0 end=false ]

[ pcftu4 Path parent=@pcftv4 layer=*CFTU4 propagator=*propc fitter=*fit_start
  checkers=@( ) stop=@0 end=false ]

[ pcftv3 Path parent=@pcftu4 layer=*CFTV3 propagator=*propc fitter=*fitk
  checkers=@( chkcs06 ) stop=@0 end=false ]

[ pcftu3 Path parent=@pcftv3 layer=*CFTU3 propagator=*propc fitter=*fitk
  checkers=@( chkcs07 ) stop=@0 end=false ]

[ pcftx5 Path parent=@pcftu3 layer=*CFTX5 propagator=*propc fitter=*fitk
  checkers=@( chkcs08 ) stop=@0 end=false ]

[ pcftv2 Path parent=@pcftx5 layer=*CFTV2 propagator=*propc fitter=*fitk
  checkers=@( chkcs09 ) stop=@0 end=false ]

[ pcftx4 Path parent=@pcftv2 layer=*CFTX4 propagator=*propc fitter=*fitk
  checkers=@( chkcs10 ) stop=@0 end=false ]

[ pcftu2 Path parent=@pcftx4 layer=*CFTU2 propagator=*propc fitter=*fitk
  checkers=@( chkcs11 ) stop=@0 end=false ]

[ pcftx3 Path parent=@pcftu2 layer=*CFTX3 propagator=*propc fitter=*fitk
  checkers=@( chkcs12 ) stop=@0 end=false ]

[ pcftv1 Path parent=@pcftx3 layer=*CFTV1 propagator=*propc fitter=*fitk
  checkers=@( chkcs13 ) stop=@0 end=false ]

[ pcftx2 Path parent=@pcftv1 layer=*CFTX2 propagator=*propc fitter=*fitk
  checkers=@( chkcs14 ) stop=@0 end=false ]

[ pcftu1 Path parent=@pcftx2 layer=*CFTU1 propagator=*propc fitter=*fitk
  checkers=@( chkcs15 ) stop=@0 end=false ]

[ pcftx1 Path parent=@pcftu1 layer=*CFTX1 propagator=*propc fitter=*fitk
  checkers=@( chkcs16 ) stop=@stop_end end=false ]