[evla-sw-discuss] OSRO python scripting: RestFreqX and VelocityX Intents
Bryan Butler
bbutler at nrao.edu
Mon Mar 8 11:24:34 EST 2010
i've been looking at the SDM, and how they store frequencies, so we can
figure out where these end up (what rich needs to do with them in MCAF).
as usual i'm confused by what's in there, but i'll put it down here and
maybe rich or michael will respond and clarify.
rest frequencies and velocities go in the Source table. the following
optional elements are the ones of interest:
Name Type Comment
numLines int Nlines, the number of line transitions
transition string[Nlines] the names of the transitions
restFrequency Frequency[Nlines] the rest frequencies
sysVel Speed[Nlines] the systemic velocity
rangeVel Speed[2] the velocity range
numLines : Number of line transitions for which information is given.
transition : The line of each spectral line transition.
restFrequency : The rest frequency for each transition considered. LABREST
sysVel : The source radial velocity in each of the line transitions. LSR
rangeVel : Range of radial velocities present in the source (for
spectroscopy).
so, as far as i can tell, MCAF should just collect up the unique values
of the RestFreq.... quantities and stuff them into Source.restFrequency,
and do the same with the Velocity... quantities into sysVel. i don't
know what happens if there is not the same number of unique RestFreq....
and Velocity.... quantities.
-bryan
David Harland wrote, On 3/5/10 08:50:
> David Harland wrote:
>> David Harland wrote:
>>> I'm modifying the code that turns an OPT-created observation into a
>>> python script. Below are some lines i want to address:
>>>
>>> intents = Intention()
>>>
>>> intents.addIntent('RestFreqA=26500')
>>> intents.addIntent('RestFreqC=26500')
>>>
>>> intents.addIntent('VelocityA=500')
>>> intents.addIntent('VelocityC=500')
>>>
>>> I want to know if we will be changing the RestFreqX and VelocityX
>>> strings.
>>>
>>> When we start doing the four 3-bit outputs, will we have RestFreqA1,
>>> RestFreqA2, ..., RestFreqD2?
>>>
>>> For the 8-bit outputs, will we continue to use RestFreqA, or will it
>>> be renamed RestFreqA0?
>>>
>>>
>> I received no feedback on the above, so below is how its been coded.
>> Showing RestFreq only below; Velocity uses same pattern.
>>
>> Case 1: Two 8-bit pairs
>>
>> intents.addIntent('RestFreqA0C0=#####')
>> intents.addIntent('RestFreqB0D0=#####')
>>
>> Case 2: Four 3-bit pairs
>>
>> intents.addIntent('RestFreqA1C1=#####')
>> intents.addIntent('RestFreqA2C2=#####')
>> intents.addIntent('RestFreqB1D1=#####')
>> intents.addIntent('RestFreqB2D2=#####')
>>
>> Case 3: One 8-bit pair & Two 3-bit pairs
>>
>> intents.addIntent('RestFreqA0C0=#####')
>> intents.addIntent('RestFreqB1D1=#####')
>> intents.addIntent('RestFreqB2D2=#####')
>>
>> or
>>
>> intents.addIntent('RestFreqA1C1=#####')
>> intents.addIntent('RestFreqA2C2=#####')
>> intents.addIntent('RestFreqB0D0=#####')
>>
>>
>>
> No replied in October, and no one objected to the above in November, so
> upcoming Doppler observations will have the lines as shown above.
> Furthermore, it was unclear as to what velocity to put in the velocity
> intent -- that of the source against its rest frame or the total
> velocity of that source relative to the VLA at the time of observation
> (which for now is the SB start time). I am providing the latter.
> _______________________________________________
> evla-sw-discuss mailing list
> evla-sw-discuss at listmgr.cv.nrao.edu
> http://listmgr.cv.nrao.edu/mailman/listinfo/evla-sw-discuss
More information about the evla-sw-discuss
mailing list