Using SupremePhoto - Tutorial
The SupremePhoto plugin consists in 2 tasks:
SupremePhotoGetDataFromObs
The supremePhotoGetDataFromObs task will retrieve the needed informations from the observations.
Input Parameters
| obsids: | (List of long) Determine a list of obsids. |
| bandName: | (String='PSW' or 'PMW' or 'PLW') Band to process. |
| poolName: | (String) Name of the pool from which to get the observations locally. |
| useHsa: | (Boolean) If you prefer to get the observations from the Herschel Science Archive and not locally. |
Output
This task will return a ListContext containing all the need informations for the supremePhotoMapMaking task:
For each obsid, a product will contain Signal, Ra, Dec, Mask, Time, Channels and Turnovers tables, as well as obsid, Ra_nom, Dec_nom and posAngle metadata.
SupremePhotoGetDataFromObs Task in HIPE.
SupremePhotoMapMaking
supremePhotoMapMaking is the main task of the plugin. This task will make the map.
Input Parameters
| AlphaStep: | (Integer) Alpha step angle in arcsec. Advice: Choose 3 for PSW, 5 for PMW and 7 for PLW. |
| RhosOp: |
(String='GetRhos' or 'SetRhos' )
'GetRhos' will activate rhos estimation jointly with the map. In this case the RhosVal is ignored. 'SetRhos' will fix Rhos to RhosVal |
| RhosVal: | (real positive) See above. |
| RhonOp: | (String='GetRhon' or 'SetRhon') see RhosOp |
| RhonVal: | (real positive) See above. |
| OffsetOp: | (String='GetOffsetLeg', 'GetOffsetAll' or 'ZeroOffset')
'GetOffsetXXX' or will activate joint estimation of detectors offsets either per scan leg (Leg) or in a global way for whole observation (All). 'ZeroOffset' will deactivate offset estimation and supposes that the offset was aleardy substracted from the data. |
| BeamOp: | (String='Simulated' or 'Empiric') Determine which beam to use.
Two beam models are possible: Simulated or Emperic. |
| PointingOp: | (String= 'Nearest' or 'Bilinear') Determine the pointing model.
Two pointing model are available: Nearest and Bilinear |
| StopMode: | (String= 'Automatic' or 'Manual') Determine stoping mode.
Two modes are available: Automatic (based on free energy evolution) or manual. |
| StopValue: | (real positive) Determine the stoping value according to stopping mode.
In automatic mode, StopValue correspond to the tolerance in Free energy variation. In manual mode, StopValue correspond to the number of iterations. |
| dirBeams: | (String) Directory that contains all the P[S,M,L]W beam files. You can download them here. |
| display: | (Integer='0' or '1') Display the evolution of the Map-Making. |
Output
This task will return some maps and additionnal informations about the computation:
| Map: | Output map & a projected map of error. |
| StdMap: | Standard deviation map from the estimator. |
| CoAddMap: | Coadded map using Bilinear projecter. |
| Extra: |
Rhon: Estimated noise Rhos: Estimated source Offset: Estimated Offsets VarRhon: the variance of Rhon estimator VarRhos: the variance of Rhos estimator VarOffset: the variance of Offsets estimator FreeEnergy: The log of free energy evolution during iterations EqBeam: Equivalent beam for the map |
SupremePhotoMapMaking Task in HIPE and its default values.
Using SupremePhoto - Tutorial Section
Here are some script examples to use in HIPE.
Example 1 - Map Making - Using HSA and Manual StopMode
clear(all=True) #USING HSA listObsids = [1342192100] # list of obsids useHsa = True # Use Herschel Science Archive Band = "PLW" # band to process # Retrieve the informations from the observations data = supremePhotoGetDataFromObs(obsids=listObsids,bandName=Band,useHsa=useHsa) #MAP MAKING # Main SupremePhoto task that return the map and additionnal informations about the computation dirBeams = "~/SupremePhotoPlugin/beams" # directory containing beams files AlphaStep = 7 # alpha step angle in arcsec. We recommend these values: PSW:3 PMW:5 PLW:7 # Exemple 1 NbrItr = 100 # number of iteration (for the default stopMode="manual") Map, StdMap, CoAddMap, Extra = supremePhotoMapMaking(data=data,AlphaStep=AlphaStep,StopValue=NbrItr,dirBeams=dirBeams) # Save the map as a fits file dirSupreme="~/SupremePhotoPlugin/" simpleFitsWriter(Map,dirSupreme+poolName+"_SupremePhoto_"+Band+"_"+str(AlphaStep)+"_"+str(NbrItr)+".fits")
Example 2 - Map Making - Using a Pool and Automatic StopMode
clear(all=True) #USING A POOL poolName = "HH_dense_core_SPIRE" # poolname listObsids = [1342192100] # list of obsids Band = "PLW" # band to process # Retrieve the informations from the observations data = supremePhotoGetDataFromObs(obsids=listObsids,bandName=Band,poolName=poolName) #MAP MAKING # Main SupremePhoto task that return the map and additionnal informations about the computation dirBeams = "~/SupremePhotoPlugin/beams" # directory containing beams files AlphaStep = 7 # alpha step angle in arcsec. We recommend these values: PSW:3 PMW:5 PLW:7 StopMode = "Automatic" StopValue = 10**(-5) Map, StdMap, CoAddMap, Extra = supremePhotoMapMaking(data=data,AlphaStep=AlphaStep,StopMode=StopMode,StopValue=StopValue,dirBeams=dirBeams) # Save the map as a fits file dirSupreme="~/SupremePhotoPlugin/" simpleFitsWriter(Map,dirSupreme+poolName+"_SupremePhoto_"+Band+"_"+str(AlphaStep)+"_"+str(NbrItr)+".fits")
