1c1
< # $Id: blastrc,v 1.41 2004/05/05 16:01:43 blast Exp $ --*-Shell-script-*--
---
> # $Id: blastrc,v 1.1.2.28 2003/07/02 14:11:08 blast Exp $ --*-Shell-script-*--
16d15
< *.DstFile:              $ANALDIR/dst-#.root
18,19d16
< *.SQL:			off
< #*.Raw.EleMap:		sql
21d17
< *.EleMap.Archive:	/home/daq/blast/pro2003/elog/archive/Electronics/
22a19,21
> *.MagField.File:	bgrid2.blast
> *.MagField.Magnitude:	-1.0
> *.MagField.Epics:	off
25,26c24
< *.Time2Dis3rdOrd.File:	WC_pos_vs_t_parameters
< *.Time2Dis6thOrd.File:  WC_pos_vs_t_parameters_2
---
> *.Time2Distance.File:	WC_pos_vs_t_parameters
29,30c27
< *.EpicsRecon.EpiCal:    epics.cal
< *.EventList:		no
---
> *.EventList:		false
32,33c29
< *.Run:			0
< *.StartEvent:		1
---
> *.StartEvent:		0
36,37d31
< *.ShowFiles		no
< *.Chain:
49,57d42
< lrn.TitleBar:		true		# display status in terminal
< lrn.Charge:		true		# fill charge tree
< 
< dst.OutFile:		$ANALDIR/dst-#.root
< dst.lz:			on		# lazy fitter
< dst.ff:			on		# fast fitter
< dst.nw:			on		# newton fitter
< dst.pid			on		# pid reconstruction
< dst.cscale:		100		# number of events per '.'
64c49,52
< d2_flrn.OutFile:	$ANALDIR/d2_flr-#.root
---
> flrn.DataFile:		$ANALDIR/lr-#.root
> sngls.OutFile:		$ANALDIR/sngls-#.root
> d2eep_flrn.OutFile:	$ANALDIR/d2eep-#.root
> d2_tof_flrn.OutFile:	$ANALDIR/d2_tof-#.root
68,69d55
< charge.OutFile:	        $ANALDIR/charge-#.chg
< crap.OutFile:	        $ANALDIR/crap-#.crap
71,107c57
< nsed.DrawMC:		on
< nsed.DrawPMT:		on
< nsed.DrawHits:		on
< nsed.DrawStubs:		on
< nsed.DrawSegments:	on
< nsed.DrawLinks:		off
< nsed.DrawFast:		on
< nsed.DrawNewt:		on
< nsed.DrawLazy:		on
< 
< nsed.PrintRaw:		off
< nsed.PrintHits:		off
< nsed.PrintClusters:	off
< nsed.PrintStubs:	off
< nsed.PrintSegments:	off
< nsed.PrintWcTracks:	off
< nsed.PrintFastBack:	off
< nsed.PrintNewt:		off
< nsed.PrintMass:		off
< nsed.PrintLazy:		off
< nsed.PrintData:		off
< nsed.PrintMC:		off
< 
< nsed.OnlyWc:		off
< nsed.OnlyTracks:	off
< nsed.TraceBack:         on
< 
< dst.DataFile:		$ANALDIR/dst-#.root
< dst.Chain:		dst
< lr.DataFile:		$ANALDIR/lr-#.root
< lr.Chain:		lr
< flr.DataFile:		$ANALDIR/flr-#.root
< flr.Chain:		flr
< flrpm.DataFile:		$ANALDIR/flr-#.root
< flrpm.Chain:		flrpm
< elas.DataFile:		$ANALDIR/elas-#.root
< elas.Chain:		elas
---
> ##charge.OutFile:     	$ANALDIR/charge-#.chg
113,114c63
< *.Recon.MaxTrack:       40
< *.Recon.MaxFastChi:	1000000		# worst Chi^2 allowed for a fast track
---
> *.Recon.MaxFastChi:	3e3		# worst Chi^2 allowed for a fast track
117,126c66,71
< *.Recon.PrefitCut1:	1000000		# prefit and select goodtracks
< *.Recon.PrefitCut2:	700000		# prefit and select goodtracks
< *.Recon.PrefitCut3:	100000		# prefit and select goodtracks
< *.Recon.PrefitCut4:	50000		# prefit and select goodtracks
< *.Recon.Max_Chi2:       6000		# maximum chi square allowed in newton
< *.Recon.Max_Z:          50		# maximum |Z| allowed in newton
< *.Recon.PRes:           0.03		# p res within which two tracks are considered same, used in screen tracks
< *.Recon.TRes:           0.01		# th resolution 0.5deg
< *.Recon.FRes:           0.01		# ph resolution 0.5deg
< *.Recon.ZRes:           1		# z resolution (cm)
---
> *.Recon.PrefitCut1:	5		# prefit and select goodtracks
> *.Recon.PrefitCut2:	0.5		# prefit and select goodtracks
> *.Recon.PrefitCut3:	0.1		# prefit and select goodtracks
> *.Recon.PrefitCut4:	200		# prefit and select goodtracks
> *.Recon.Max_Chi2:       150		# maximum chi square allowedin newton,
> *.Recon.Max_Z:          50		# maximum chi square allowedin newton,
128,133c73,78
< *.Recon.Refit:		on ###		# refit with pid/beta information
< *.Recon.ScreenTrack:	on		# screen tracks after newton fit
< nsed.Recon.ScreenTrack: off
< *.Recon.PickZ:		on		# screen tracks based on common vertex
< nsed.Recon.PickZ:       off
< *.Recon.PickBest:	off		# pick best track in each sector
---
> *.Recon.Refit:		on		# refit with pid/beta information
> *.Recon.PickZ:		on
> *.Recon.Lazy:		on
> *.Recon.LazySwim:	off
> #nsed.Recon.LazySwim:	on
> *.Recon.LazyNtuple:	off
147,197d91
< # Geometry
< # --------------
< 
< *.GeomCC.DetColor:	4
< *.GeomCC.DetWidth:	1
< *.GeomCC.HitColor:	4
< *.GeomCC.HitWidth:	2
< 
< *.GeomLG.DetColor:	7
< *.GeomLG.DetWidth:	1
< *.GeomLG.HitColor:	7
< *.GeomLG.HitWidth:	2
< 
< *.GeomNC.DetColor:	6
< *.GeomNC.DetWidth:	1
< *.GeomNC.HitColor:	6
< *.GeomNC.HitWidth:	2
< 
< *.GeomRD.DetColor:	2
< *.GeomRD.DetWidth:	1
< *.GeomRD.HitColor:	4
< *.GeomRD.HitWidth:	2
< 
< *.GeomSC.DetColor:	5
< *.GeomSC.DetWidth:	1
< *.GeomSC.HitColor:	2
< *.GeomSC.HitWidth:	2
< 
< *.GeomLD.DetColor:	6
< *.GeomLD.DetWidth:	1
< *.GeomLD.HitColor:	6
< *.GeomLD.HitWidth:	3
< 
< *.GeomWC.DetColor:	3
< *.GeomWC.DetWidth:	1
< 
< UtilPad.Zoom:		2		# zoom factor (0 print coords instead)
< UtilPad.TopX1:		-240		# default top view
< UtilPad.TopY1:		-480
< UtilPad.TopX2:		720
< UtilPad.TopY2:		480
< UtilPad.FrontX1:	-480		# default front view
< UtilPad.FrontY1:	-480
< UtilPad.FrontX2:	480
< UtilPad.FrontY2:	480
< 
< *.Vector3.Color:	1
< *.Vector3.Style:	5
< *.Vector3.ArrowSize	0.02
< 
< # --------------
201,206c95,97
< *.MagField.Class:	TBLMagField
< *.MagField.File:	bgrid2.blast
< *.MagField.Magnitude:	1.		# scale fieldmap file by this factor
< *.MagField.Epics:	on		# read new fMagnitude from coda file
< *.FieldCoil.NSeg:	100		# analytical field
< *.FieldConst.Bx:	0.		# fixed field (same in both sectors)
---
> *.Field.Class:		TBLMagField
> *.FieldCoil.NSeg:	100
> *.FieldConst.Bx:	0.
209,210c100
< *.FieldRho.B:		1.5		# azimuthal field
< *.FieldTr.Dx:		0.		# TMagField w/translation
---
> *.FieldTr.Dx:		0.
213c103
< *.FieldRandom.Seed:	11		# TMagField w/random errors
---
> *.FieldRandom.Seed:	11
215c105
< *.FieldRandom.Dx:	0.		# sigma
---
> *.FieldRandom.Dx:	0.
224,227c114,115
< *.WireCal.Color		3		# all wires
< *.WireCal.Marker	6		#  dot
< *.Wire.Color:		1		# single hit
< *.Wire.Marker:		25		#  box
---
> *.WireCal.Color:	1
> *.WireCal.Marker:	25
232,233c120,121
< *.Hit.Vel:		0.00182		# 2.02/1100.
< *.Hit.MaxHits:		1000		# max hits / event
---
> *.Hit.Vel:		0.0182		# 20.2/1100.
> *.Hit.MaxHits:		500		# max hits / event
236,240d123
< *.Hit.SixthOrder:       on              # on = 6th order Garfield, off = 3rd
< *.Hit.TdcMin:		2000		# 2000=8us (normal range=2us)
< *.Hit.TdcMax:		6600		# full range (close to the wire)
< *.Hit.DistMin:		0		# closest to wire
< *.Hit.DistMax:		9999		# farthest from wire
250a134
> lrn.Stub.Cut:		180.		# keep all chisq's < Cut
252,264c136,138
< *.Stub.CUTTIGHT:	10.
< *.Stub.CUTFORWARD:      320.            # 320 for 797 ev 950
< *.Stub.CUTBACKWARD:     120.
< *.Stub.MaxStub:		500		# max stubs / event
< *.Stub.MaxSlope10:     -2.0		# max slope pointing downstream in ch0
< *.Stub.MaxSlope11:     -2.6		# max slope pointing downstream in ch1
< *.Stub.MaxSlope12:     -3.5		# max slope pointing downstream in ch2
< *.Stub.MaxSlope20:      0.5		# max slope pointing upstream in ch0
< *.Stub.MaxSlope21:      1.5
< *.Stub.MaxSlope22:      2.0
< *.Stub.MaxDx:           60		# when extended to next chamber, max
< *.Stub.Select_Stub:     on
< *.Stub.Screen_Stub:     off
---
> *.Stub.NHits:		1		# minimal # of hits to form a stub
> *.Stub.MaxStub:		1000		# max stubs / event
> *.Stub.MaxSlope:        3
270,273c144
< *.Segment.MaxDPhi:	0.4		# max dPhi between segments(pi/12)
< *.Segment.MaxDX:	50		# max dx between segments(pi/12)
< *.Segment.MaxFake::	3		# max # of fake segments in one chamber
< *.Segment.MaxSegChi2:   50              # max chi squared of segment toloerable
---
> *.Segment.MaxDPhi:	0.22;		# max dPhi between segments(pi/12)
275,283c146,150
< *.Segment.Select_Seg:	on		# call select_seg() to clean up
< *.Segment.Slope1:       0.47		# phi of an segment center, tan(25 deg)=0.47
< *.Segment.Slope2:       2.5		# min tilt of seg projected onto XY plane. >=6
< *.Segment.SlopeX11:      0.0
< *.Segment.SlopeX12:      1.5
< *.Segment.SlopeX13:      2.5
< *.Segment.SlopeX21:      -2.0
< *.Segment.SlopeX22:      -2.7
< *.Segment.SlopeX23:      -3.5
---
> *.Segment.Slope1:       0.47		# phi of an segment center,25 deg
> *.Segment.Slope2:       2.0		# phi of an segment 60 deg
> *.Segment.Slope3:       0.8		# theta toward upstream
> *.Segment.Slope4:      -2.2             # theta toward downstream
> *.Segment.Slope5:       1.4             # phi between track and segment
285,287c152
< *.Segment.Color0:	7
< *.Segment.Color1:	8
< *.Segment.Color2:	9
---
> *.Segment.Color:	7
291,292c156,157
< *.WcTrack.MaxTrks:	200		# max fast tracks / event
< *.WcTrack.MaxXPhi:	0.5		# max chi square in phi
---
> *.WcTrack.MaxTrks:	500		# max fast tracks / event
> *.WcTrack.MaxXPhi:	0.05		# max chi square in phi
295c160
< *.WcTrack.Color:	41		# plain circle track
---
> *.WcTrack.Color:	6
297,300d161
< *.WcTrack.Extend:	20		# fast fit track past end segs (cm)
< *.WcTrack.SimColor:	3		# simulated track
< *.WcTrack.SimWidth:	1
< *.WcTrack.Style:	4		# segments
307c168
<  nsed.SimTrack.MaxRadius: 600.
---
> nsed.SimTrack.MaxRadius: 500.
309c170
<  nsed.SimTrack.MaxSteps: 200
---
> nsed.SimTrack.MaxSteps:	 100
311c172
<  nsed.SimTrack.MaxPathLen: 600.
---
> nsed.SimTrack.MaxPathLen: 500.
313c174
<  nsed.SimTrack.MaxX:	 500.
---
> nsed.SimTrack.MaxX:	 400.
315c176
<  nsed.SimTrack.MaxZ:	 600.
---
> nsed.SimTrack.MaxZ:	 500.
319,323d179
< *.SimTrack.Color:	1
< *.SimTrack.Width:	1
< 
< *.Raw.MCColor:		4
< *.Raw.MCWidth:		1
325,327c181
< *.FitTrack.Spline:	true		# interpolate simulations using spline
< *.FitTrack.ChamberFit:	true		# simulate tracks only in chamber
< *.FitTrack.PhiZero:	0.		# scale ff phi towards center plane
---
> *.FitTrack.Spline:	true
329,335d182
< *.FitTrack.OnVertex:	true
< *.FitTrack.CorrectBeta:		true
< *.FitTrack.CorrectDistance:	false
< *.FitTrack.CorrectImpact:	false
< *.FitTrack.CorrectY:		true
< *.FitTrack.CorrectB:		true
< 
345,348d191
< *.FitTrack.Color:	2
< *.FitTrack.Width:	1
< *.FitTrack.VertexColor	1
< *.FitTrack.VertexStyle	5
353c196
< *.Newt.Refit:		true		# refit with negative charge
---
> *.Newt.Refit:		true
369,382d211
< *.PID.MAX_NC_X2:	160000		# when a WC track cross NC, and is less 30 cm away from
< 					# a NC hit, the hit is considered by the WC track, not a neutron
< *.Track.CorrectBeta:    on
< 
< *.Lazy:			on		# use lazy fitter
< *.Lazy.Swim:		off		# swim track back to vertex
<  nsed.Lazy.Swim: 	on
< *.Lazy.Ntuple:		off		# output ntuple
< *.Lazy.PhiTOF:		true		# use tof info for phi
< *.Lazy.MaxStub:		500		# maximum stubs used per sector
< *.Lazy.Color:		1
< *.Lazy.Width:		1
< *.Lazy.Markers:		true		# draw markers for segs and tofs
< 
389d217
< *.Stub.PrintFlag:       0
394,395c222,226
<  nsed.Recon.PrintFlag:	 1
< *.Lazy.Print:		0		# 1=summary 2=details
---
> nsed.Recon.PrintFlag:    0
> *.Recon.LazyPrint:	0		# 1=summary 2=details
> 
> *.Newt.PrintFlag:	0
> nsed.Newt.PrintFlag:	 0
397,399d227
< *.Newt.PrintFlag:	0		# 1:flow information, results
<  nsed.Newt.PrintFlag:	 0		# 2:iterations
< 					# 3:derivatives,data
400a229,230
> 
> #wc cutoff at 2000