Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://xmm.vilspa.esa.es/sas/6.5.0/doc/rgssources/node30.html
Дата изменения: Tue Aug 16 05:31:58 2005 Дата индексирования: Sat Dec 22 12:15:46 2007 Кодировка: Поисковые слова: hst |
read parameters; ################# start format upgrade section. if (filemode eq 'modify') { # convert old-style source list if necessary: if (! instrument_supplied) { if (hasAttribute(set, 'INSTRUME')) { instrument = stringAttribute(oldListSet, 'INSTRUME'); } else { # attempt to get it from file name: if (! instrument_in_filename) {call error;} } } call CAL_setState(instrument=instrumentId); rename table to 'SRCLIST'; get attitude from *_REF keywords; write_generic_keywords; if (writeobskwds) {write_observation_keywords;} if (writeexpkwds) { if (! exposure_supplied) { if (hasAttribute(set, 'EXP_ID')) { exposure = stringAttribute(oldListSet, 'EXP_ID'); } else { # attempt to get it from file name: if (! exposure_in_filename) {call error;} } } } # Convert columns: n_offaxis_sources = 0; foreach (source) { fixed_on_sky = true; # default process = false; # default bkg_exclude = false; # default if (old_label eq 'proposal') { new_label = 'PROPOSAL'; new_ra = old_ra; new_dec = old_dec; } elsif (label eq 'onaxis') { new_label = 'ON_AXIS'; new_ra = attitude_ra; new_dec = attitude_dec; delta_disp = 0; delta_xdsp = 0; new_fov_phi = 0; new_fov_r = 0; fixed_on_sky = false; } elsif (label eq 'epic') { new_label = 'OLDEPIC'.id_band.ml_id_src; new_ra = old_ra; new_dec = old_dec; epic_file = number of bands found; } elsif (label eq 'offaxis') { n_offaxis_sources = n_offaxis_sources + 1; new_label = 'OFFAXIS'.n_offaxis_sources; delta_disp = rgs_disp; delta_xdsp = rgs_xdsp; calculate new_fov_phi, new_fov_r; fixed_on_sky = false; } else { # assume is user-added source. new_label = upper_case(old_label); new_ra = old_ra; new_dec = old_dec; } new_rate = old_rate; if (src_select > 0) {process = true;} if (back_select > 0) {bkg_exclude = true;} if (fixed_on_sky) { calculate delta_disp; calculate delta_xdsp; calculate new_fov_phi; calculate new_fov_r; } else { calculate new_ra; calculate new_dec; } } calculate_confusion(src_data, prime_index); delete old columns; write column data to new columns; write epic keywords; } ################# end format upgrade section. # get the instrument: if (filemode eq 'modify') { instrument = stringAttribute(oldListSet, 'INSTRUME'); } else { # get instrument from parameter; if (instrument parameter at default value) {call error;} } # get the exposure number: expid_is_valid = true; if (filemode eq 'modify') { exposure = int32Attribute(oldListSet, 'EXP_ID'); } else { # get exposure from parameter; if (exposure parameter at default value) { expid_is_valid = false; } } call CAL_setState(instrument=instrumentId); # get the pointing: if (filemode .eq. 'modify' && ! changeAttitude) { # read attitude from old source list; } elsif (attitudeStyle .eq. 'user') { # read attitude from attra, attdec and attapos parameters; } else { if (!expid_is_valid) {call error;} call OAL_setState(instrument=instrument, exposureNr=exposure); if (attitudeStyle .eq. 'start') { # get attitude from attitude at start of observation; } elsif (attitudeStyle .eq. 'mean') { # get attitude from mean of attitude history file values; } elsif (attitudeStyle .eq. 'median') { # get attitude from median of attitude history file values; } } allocate(src_data(max_number_sources)); if (filemode eq 'create') { proposal position -> src_data(1); attitude position -> src_data(2); } else { # filemode eq 'modify' open(old rgs set, 'READ'); call check(old rgs list); # This checks that all label names are unique, # and that there is a proposal and attitude source. old rgs sources -> src_data; release(old rgs set); if (changeattitude) {attitude position -> src_data(index of onaxis src);} } if (withepicset) { clone(old epic set); if (enablefilter) {spatial filtering of cloned epic set;} foreach(epic source) { label = 'epiclabelprefix'.column('*_ID_SRC'); foreach(source already in src_data) { if (label eq src_data(LABEL)) { no_label_clash = F; if (clobberonlabel) { last; } else { call error; } } } if (no_label_clash) { epic source -> src_data; } } release(old epic set); } if (addusersource) { foreach(source already in src_data) { if (label eq src_data(LABEL)) { no_label_clash = F; if (clobberonlabel) { last; } else { call error; } } } if (no_label_clash) { user source -> src_data; } } # Possible routine in this place that filters out (or at least identifies) # all pairs of too-close sources. if (changeattitude or null in column) { foreach(source in src_data) { # Calculate ra/dec or delta_disp/delta_xdsp/fov_phi/fov_r, as appropriate; } } # Write the data to file: if (filemode eq 'create') { open(new rgs set, 'CREATE'); # write keywords to header: write_generic_keywords; if (writeobskwds) {write_observation_keywords;} if (writeexpkwds && expid_is_valid) {write_exposure_keywords;} } else { # filemode eq 'modify' open(old rgs set, 'MODIFY'); if (changeattitude) { rewrite attitude attributes (#RA_PNT etc); delete all region extensions; set PROCESS and BKG_EXCLUDE to false; } delete all table entries; } src_data -> new rgs set; if (withepicset) { write_epic_specific_keywords; } if (enablefilter) {check that all members of the src list are within the fov, warn if not;} # Now calculate prime source index value: if (filemode eq 'modify' && !changeprime) { prime_index = prime_index in old file; } else { if (addusersource && userasprime) { prime_index = index of (last) user-added source; } else { if (primestyle eq 'label') { prime_index = index of source with label=primesourcelabel; } elsif(primestyle eq 'index') { prime_index = primesourceindex; } elsif(primestyle eq 'brightest') { prime_index = index of brightest source; } elsif(primestyle eq 'auto') { if (proposal source is within FOV) { prime_index = index of proposal source; } elsif (there are some epic sources within FOV) { prime_index = index of brightest within-FOV epic source; } else { prime_index = index of pointing (ONAXIS) source; } } } } else { prime_index = previous prime_index; } if (filemode eq 'create' || changeprime) { prime_index -> #PRIMESRC; } call calculate_confusion(src_data, prime_index); # Write the confusion data to file: # release the data set.