Changes between Version 25 and Version 26 of Notes/ExampleReduction


Ignore:
Timestamp:
04/01/09 13:30:49 (15 years ago)
Author:
JamieStevens
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Notes/ExampleReduction

    v25 v26  
    348348need to make the beam 3 times the size of the area being cleaned, instead of using options=double which simply doubles
    349349the beam size compared to the image. To do this we need to make the image 3 times bigger than the primary beam and then
    350 cut the central part of the image out using imframe. Since the parameters we need to do this will be different for each
     350specify the normal imaging region during mfclean. Since the parameters we need to do this will be different for each
    351351central frequency, we will control dac_repeat.pl with a parameter file.
    352352
     
    381381
    382382{{{
    383 4540 cell=26.672,22.525 imsize=255,302
    384 4668 cell=25.265,21.494 imsize=262,308
    385 4796 cell=24.635,20.907 imsize=261,308
    386 4924 cell=23.992,20.378 imsize=261,308
    387 5052 cell=23.463,19.901 imsize=260,307
    388 5180 cell=22.910,19.483 imsize=260,306
    389 5308 cell=22.394,19.020 imsize=260,306
    390 5436 cell=21.931,18.553 imsize=259,306
    391 5564 cell=21.531,18.208 imsize=258,305
    392 5692 cell=21.018,17.827 imsize=258,304
    393 5820 cell=20.568,17.428 imsize=258,304
    394 5948 cell=20.152,17.082 imsize=257,304
    395 6076 cell=19.800,16.781 imsize=257,303
    396 6204 cell=19.379,16.469 imsize=257,302
    397 6332 cell=19.053,16.110 imsize=256,303
    398 6460 cell=19.330,16.151 imsize=247,296
     3834540 cell=8.891,7.508 imsize=255,302
     3844668 cell=8.422,7.165 imsize=262,308
     3854796 cell=8.212,6.969 imsize=261,308
     3864924 cell=7.997,6.793 imsize=261,308
     3875052 cell=7.821,6.634 imsize=260,307
     3885180 cell=7.637,6.494 imsize=260,306
     3895308 cell=7.465,6.340 imsize=260,306
     3905436 cell=7.310,6.184 imsize=259,306
     3915564 cell=7.177,6.069 imsize=258,305
     3925692 cell=7.006,5.942 imsize=258,304
     3935820 cell=6.856,5.809 imsize=258,304
     3945948 cell=6.717,5.694 imsize=257,304
     3956076 cell=6.600,5.594 imsize=257,303
     3966204 cell=6.460,5.490 imsize=257,302
     3976332 cell=6.351,5.370 imsize=256,303
     3986460 cell=6.443,5.384 imsize=247,296
    399399}}}
    400400
     
    403403dac_repeat.pl --task invert --repeat vis --repeat-freq rnovalmc --repeat map --repeat-freq rnovalmc_m2_%f.map --repeat beam --repeat-freq rnovalmc_m2_%f.beam "select=-auto,-ant(6)" stokes=ii options=mfs,sdb slop=1 --args-file dac_invert_args.txt
    404404
    405 This makes the triple-sized maps and beams, and now we have to chop out the normal sized map
    406 using imframe. We make the file dac_imframe_args.txt that looks like this:
     405This makes the triple-sized maps and beams, and now we have to specify the normal sized map
     406while using mfclean. We make the file dac_mfclean_args.txt that looks like this:
    407407
    408408{{{
     
    425425}}}
    426426
    427 And we run imframe using dac_repeat.pl like this:
    428 
    429 dac_repeat.pl --task imframe --repeat in --repeat-freq rnovalmc_m2_%f.map --repeat out --repeat-freq rnovalmc_m2_small_%f.map --args-file dac_imframe_args.txt
    430 
     427This is handy since mfclean also requires that the "region being cleaned should be reasonably
     428centred in the map, and should have an appreciable guard band around it to the map edge (of
     429size comparable to the width of the region being cleaned)".
     430
     431We run mfclean using dac_repeat.pl like this:
     432
     433dac_repeat.pl --task mfclean --repeat map --repeat-freq rnovalmc_m2_%f.map --repeat beam --repeat-freq rnovalmc_m2_%f.beam --repeat out --repeat-freq rnovalmc_m2_%f.clean --args-file dac_mfclean_args.txt
     434
     435We can still use restor with MFS datasets, but it will create a triple-sized image, although
     436of course only the central region will be cleaned.
     437
     438dac_repeat.pl --task restor --repeat model --repeat-freq rnovalmc_m2_%f.clean --repeat beam --repeat-freq rnovalmc_m2_%f.beam --repeat map --repeat-freq rnovalmc_m2_%f.map --repeat out --repeat-freq rnovalmc_m2_%f.restor
     439
     440If you only want to see the primary beam in your restor'd images, then you can use imframe
     441using the same dac_mfclean_args.txt:
     442
     443dac_repeat.pl --task imframe --repeat in --repeat-freq rnovalmc_m2_%f.restor --repeat out --repeat-freq rnovalmc_m2_small_%f.restor --args-file dac_mfclean_args.txt
     444