lsst::meas::extensions::psfex::detail::PsfexPsfFactory Class Reference

Inheritance diagram for lsst::meas::extensions::psfex::detail::PsfexPsfFactory:

Public Member Functions
virtual std::shared_ptr< table::io::Persistable >	read (InputArchive const &archive, CatalogVector const &catalogs) const
	PsfexPsfFactory (std::string const &name)

Detailed Description

Definition at line 377 of file PsfexPsf.cc.

Constructor & Destructor Documentation

PsfexPsfFactory()

lsst::meas::extensions::psfex::detail::PsfexPsfFactory::PsfexPsfFactory ( std::string const & name )

inlineexplicit

Definition at line 464 of file PsfexPsf.cc.

: table::io::PersistableFactory(name) {}

Member Function Documentation

read()

virtual std::shared_ptr< table::io::Persistable > lsst::meas::extensions::psfex::detail::PsfexPsfFactory::read ( InputArchive const & archive, CatalogVector const & catalogs ) const

inlinevirtual

Definition at line 381 of file PsfexPsf.cc.

                                                                         {
    LSST_ARCHIVE_ASSERT(catalogs.size() == 2u);
 
    std::shared_ptr<PsfexPsf> result(new PsfexPsf());
 
    int ndim, ngroup, ncoeff;
    int size_size, comp_size, context_size;
    {
        PsfexPsfSchema1 const & keys = PsfexPsfSchema1();
        LSST_ARCHIVE_ASSERT(catalogs[0].size() == 1u);
        LSST_ARCHIVE_ASSERT(catalogs[0].getSchema() == keys.schema);
        table::BaseRecord const & record = catalogs[0].front();
 
        // fields in _poly
        ndim = record.get(keys.ndim);
        ngroup = record.get(keys.ngroup);
        ncoeff = record.get(keys.ncoeff);
        // Other scalars
        result->_averagePosition = record.get(keys.averagePosition);
        result->_pixstep = record.get(keys._pixstep);
        // sizes of vectors
        size_size = record.get(keys._size_size);
        comp_size = record.get(keys._comp_size);
        context_size = record.get(keys._context_size);
    }
    // Now we can read the data
    {
        PsfexPsfSchema2 const keys(ndim, ngroup, ncoeff,
                                   size_size, comp_size, context_size);
 
        LSST_ARCHIVE_ASSERT(catalogs[1].size() == 1u);
        LSST_ARCHIVE_ASSERT(catalogs[1].getSchema() == keys.schema);
        table::BaseRecord const & record = catalogs[1].front();
 
        // _poly
        std::vector<int> group;
        {
            int const *begin = record.getElement(keys.group);
            group.assign(begin, begin + ndim);
 
            for (int i = 0; i != ndim; ++i) {
                ++group[i];             // poly_init subtracts 1 from each element.  Sigh.
            }
        }
        std::vector<int> degree;
        {
            int const *begin = record.getElement(keys.degree);
            degree.assign(begin, begin + ngroup);
        }
        result->_poly = poly_init(&group[0], group.size(), &degree[0], degree.size());
        LSST_ARCHIVE_ASSERT(result->_poly->ncoeff == ncoeff);
 
        {
            double const *begin = record.getElement(keys.basis);
            std::copy(begin, begin + ncoeff, result->_poly->basis);
        }
        {
            double const *begin = record.getElement(keys.coeff);
            std::copy(begin, begin + ncoeff, result->_poly->coeff);
        }
        // vectors
        {
            int const *begin = record.getElement(keys._size);
            result->_size.assign(begin, begin + size_size);
        }
        {
            float const *begin = record.getElement(keys._comp);
            result->_comp.assign(begin, begin + comp_size);
        }
        {
            double const *begin1 = record.getElement(keys._context_first);
            double const *begin2 = record.getElement(keys._context_second);
            result->_context.resize(context_size);
            for (int i = 0; i != context_size; ++i) {
                result->_context[i].first  = begin1[i];
                result->_context[i].second = begin2[i];
            }
        }
    }
 
    return result;
}

---

The documentation for this class was generated from the following file:

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-8.0.0/Linux64/meas_extensions_psfex/g2305ad1205+a8e3196225/src/PsfexPsf.cc