LSSTApplications  20.0.0
LSSTDataManagementBasePackage
Namespaces | Functions
psfexAdaptors.cc File Reference
#include "boost/format.hpp"
#include "lsst/pex/exceptions.h"
#include "lsst/meas/extensions/psfex/Field.hh"
#include "globals.h"
#include "context.h"
#include "prefs.h"
#include "sample.h"

Go to the source code of this file.

Namespaces

 lsst
 A base class for image defects.
 
 lsst::meas
 
 lsst::meas::extensions
 
 lsst::meas::extensions::psfex
 

Functions

setstruct * load_samples (char **filenames, int catindex, int ncat, int ext, int next, contextstruct *context)
 
void lsst::meas::extensions::psfex::makeit (std::vector< std::shared_ptr< Field > > &fields_, std::vector< std::shared_ptr< Set > > const &sets)
 

Function Documentation

◆ load_samples()

setstruct* load_samples ( char **  filenames,
int  catindex,
int  ncat,
int  ext,
int  next,
contextstruct *  context 
)

Definition at line 14 of file psfexAdaptors.cc.

16 {
17  /*
18  * The C version of this is called two ways:
19  * catindex == 0, ncat == ncat Read all catalogues
20  * catindex == c, ncat == 1 Read only catalogue c
21  */
22  setstruct *completeSet = reinterpret_cast<setstruct *>(filenames[catindex + 0]);
23  /*
24  * Make a new set, which may be a subset of the completeSet
25  */
26  setstruct *set = init_set(context);
27  set->fwhm = completeSet->fwhm;
28  for (int i = 0; i != completeSet->vigdim; ++i) {
29  set->vigsize[i] = completeSet->vigsize[i];
30  }
31  for (int i = 0; i != completeSet->ncontext; ++i) {
32  strcpy(set->contextname[i], completeSet->contextname[i]);
33  set->contextoffset[i] = completeSet->contextoffset[i];
34  set->contextscale[i] = completeSet->contextscale[i];
35  }
36  /*
37  * Count how many samples we'll be including
38  */
39  int nsample_keep = 0;
40  for (int i = 0; i != ncat; ++i) {
41  setstruct *s = reinterpret_cast<setstruct *>(filenames[catindex + i]);
42  for (int j = 0; j != completeSet->nsample; ++j) {
43  samplestruct const *samp = s->sample[j];
44  if (ext == ALL_EXTENSIONS || ext == samp->extindex) {
45  ++nsample_keep;
46  }
47  }
48  }
49 
50  set->samples_owner = 0;
51  malloc_samples(set, nsample_keep);
52  for (int i = 0; i != ncat; ++i) {
53  setstruct *s = reinterpret_cast<setstruct *>(filenames[catindex + i]);
54  for (int j = 0; j != completeSet->nsample; ++j) {
55  samplestruct *samp = s->sample[j];
56  if (ext == ALL_EXTENSIONS || ext == samp->extindex) {
57  set->sample[set->nsample++] = samp;
58  }
59  }
60  }
61 
62  return set;
63 }
std::strcpy
T strcpy(T... args)
set
daf::base::PropertySet * set
Definition: fits.cc:912