Tags:
create new tag
view all tags

boloSource

The boloSource program takes PACS level 1 frames and subtracts embedded- and foreground galactic sources. BoloSource preserves the statistical properties of the diffuse background while 'disturbing' sources are being subtracted.

The boloSource algorithm works on the flux calibrated detector timeline, it can separate sky signal from detector 1/f noise on a statistical basis and applies an adaptive interpolation method which simulates the power spectrum of instrument noise. These interpolated source-free timelines are re-projected onto a map using inversion algorithms standard for Herschel data processing. Besides the good quality of point-source subtracted maps, the boloSource algorithm provides an alternative way of source photometry in crowded regions, and can mitigate the impact of reconstruction artifacts of mapping algorithms.

For more information please refer to:

First results with the boloSource() algorithm: Photometry of faint standard stars observed by Herschel/PACS by Marton, G. Vavrek, R. Kiss, Cs. & Muller, Th. G.

and

boloSource(): Motivation, concept, & implementation, first results, and perspectives by Roland Vavrek and Gabor Marton. Presented at the Map-Maker Workshop at ESAC, on January 31, 2013.

The boloSource program uses the Unimap map maker, for more information on Unimap please see:

Unimap documentation homepage

and

Unimap user manual

Using boloSource

After installing the plugin in HIPE and unpacking the folders in the jar file go to the scripts/unimap/ folder. This folder includes; MCRInstaller.bin - the Matlab Compiler Runtime installer, run_unimap.sh - the shell for launching unimap, unimap - the unimap executable, and unimap_par.txt - the unimap parameter file (this folder also includes the Unimap user manual where the Unimap installation is also discussed on page 20).

Then install Matlab Compiler Runtime. First make MCRInstaller.bin executable and then launch it. It will tell you the directory of the library, for example "/opt/MATLAB/MATLAB_Compiler_Runtime/v78/". In the BOLOSOURCE.py script change the MCR variable to the directory of the library.

Next make the run_unimap.sh and unimap files executable. The directory where these files and the unimap_par.txt files are located is the "unimap_path" and should be the same as the unimap_path variable in the BOLOSOURCE.py script.

After this is done change the parameters in BOLOSOURCE.py. The parameters are explained below:

  • dirHome: is the location of the un-jarred boloSource folder
  • unimap_path: is the location of run_unimap.sh, unimap, and unimap_par.txt as explained above
  • MCR: is the location of the Matlab Compiler Runtime library
  • LocalFrames: is the name of the folder where the local level 1 frames are located, the local level 1 frames should be named "obsid.fits", for example "1342238774.fits"
  • LocalFramesPath: is the directory of the LocalFrames folder
  • obsids: are the obsids that you want to be run through the program
  • BoloSourceDataDir: is the directory in which the results of the boloSource program will be saved in
  • SourceMaskfilename: is the directory and name of the text or fits source mask file

Finally in the last line of the script is a function with two additional parameter that can be defined aside from the ones mentioned above, source and bolo.

  • source: if source is equal to HSA, boloSource will retrieve the level 1 frames for the obsids stated above from the Herschel Science Archive, if source is equal to Local, boloSource will use the local level 1 frames for the obsids stated above that are placed the in the LocalFrames folder

  • bolo: if bolo is equal to 1, boloSource will be run on the local Level 1 frames placed in the LocalFrames folder, if bolo is equal to 2, boloSource will be run on the local Level 1 frames placed in the LocalFrames folder, and creates maps of the interpolated product, the original product, the diffmap, the noise map, artifact and no artifact maps, and coverage map
Edit | Attach | Watch | Print version | History: r3 < r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions
Topic revision: r3 - 2013-11-20 - ThomasWilson
 
This site is powered by the TWiki collaboration platform Powered by Perl