SupremeSpectro

The SupremeSpectro Plugin v1.4 is now available :

SupremePhoto

The SupremePhoto Plugin v1.4 is now available :

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")