ডিপথট:Turning half of the tiles of all stations off

ডিপথট থেকে

Prediction and BBS calibration

Sky model: 3C196.sky

sim.parset

Strategy.InputColumn = DATA
Strategy.ChunkSize = 300
Strategy.Baseline = [CR]S*&
Strategy.Steps = [predict]

Step.predict.Operation = PREDICT
Step.predict.Model.Sources = []
Step.predict.Model.Beam.Enable = T
Step.predict.Output.Column = DATA

cal.parset

Strategy.InputColumn = DATA
Strategy.ChunkSize = 300
Strategy.Baselines = [CR]S*&
Strategy.UseSolver = F
Strategy.Correlations = []
Strategy.Steps = [solve,correct]

Step.solve.Operation = SOLVE
Step.solve.Model.Sources = []
Step.solve.Model.Cache.Enable = T
Step.solve.Model.Gain.Enable = T
Step.solve.Model.Beam.Enable = T
Step.solve.Solve.Parms = ["Gain*"]
Step.solve.Solve.ExclParms = []
Step.solve.Solve.CalibrationGroups = []
Step.solve.Solve.CellSize.Freq = 0
Step.solve.Solve.CellSize.Time = 1
Step.solve.Solve.CellChunkSize = 160
Step.solve.Solve.UVRange = []
Step.solve.Solve.PropagateSolutions = F
Step.solve.Solve.Options.MaxIter = 40
Step.solve.Solve.Options.EpsValue = 1e-9
Step.solve.Solve.Options.EpsDerivative = 1e-9
Step.solve.Solve.Options.ColFactor = 1e-9
Step.solve.Solve.Options.LMFactor = 1.0
Step.solve.Solve.Options.BalancedEqs = F
Step.solve.Solve.Options.UseSVD = T

Step.correct.Operation = CORRECT
Step.correct.Model.Sources = []
Step.correct.Model.Cache.Enable = T
Step.correct.Model.Gain.Enable = T
Step.correct.Model.Beam.Enable = F
Step.correct.Output.Column = CORRECTED_DATA

Bash script to run the experiment:

 0 # 3C196.tilex03.MS, sim.parset, cal.parset, 3C196.sky
 1 # Turn half of the tiles of all the stations off before calibrating
 2 ms=ms/3C196.tilex03.MS
 3 rm -rf $ms
 4 echo ... creating $ms
 5 cp -r ms/3C196.root.MS $ms
 6 echo ... predicting in DATA of $ms
 7 calibrate-stand-alone -v -f $ms pset/sim.parset sky/3C196.sky &> log/3C196.tilex03.sim.log
 8 echo ... turning half of the tiles of all the stations off in $ms
 9 python py/elementflag.py $ms w all half
10 echo ... calibrating and saving in CORRECTED_DATA of $ms
11 calibrate-stand-alone -v -f $ms pset/cal.parset sky/3C196.sky &> log/3C196.tilex03.cal.log

Inspect visibilities

uvplot.py -i ms/3C196.tilex03.MS/ -c DATA-CORRECTED_DATA -d /png -e 0..3 -m

3C196.tilex03.visibilities.png

python py/vis.py ms/3C196.tilex03.MS

Residual RMS in Jy and %, I: 1 (3.8%), Q: 0.058 (2.6%), U: 0.065 (2.6%), V: 0.0037 (19%).

Imaging

SAGECal for 10 directions