PhysicsTools::VarProcessor Class Reference

Common base class for variable processors. More...

#include <VarProcessor.h>

Inheritance diagram for PhysicsTools::VarProcessor:

Classes
class	Config
	Helper class for discriminator computer set-up procedure. More...
class	ConfigCtx
class	ConfIterator
	Iterator to loop over the input/output variable configuration. More...
struct	Dummy
class	LoopCtx
	Hold context information for looping processors. More...
struct	ValueIterator

Public Types
enum	LoopStatus { kStop, kNext, kReset, kSkip }
typedef Dummy *	PluginFunctionPrototype()

Public Member Functions
void	configure (ConfigCtx &config)
	called from the discriminator computer to configure processor More...
void	deriv (double *input, int *conf, double *output, int *outConf, int *loop, LoopCtx &ctx, unsigned int offset, unsigned int in, unsigned int out, std::vector< double > &deriv) const
	run the processor evaluation pass on this processor and compute derivatives More...
void	eval (double *input, int *conf, double *output, int *outConf, int *loop, LoopCtx &loopCtx, unsigned int offset) const
	run the processor evaluation pass on this processor More...
virtual LoopStatus	loop (double *output, int *outConf, unsigned int nOutput, LoopCtx &ctx, unsigned int &nOffset) const
virtual	~VarProcessor ()
Public Member Functions inherited from PhysicsTools::ProcessRegistry< Base_t, CalibBase_t, Parent_t >::Factory
template<>
VarProcessor *	create (const char *name, const Calibration::VarProcessor *calib, const MVAComputer *parent)
template<>
VarProcessor *	create (const char *, const Calibration::VarProcessor *, const MVAComputer *)

Protected Types
typedef ProcessRegistry< VarProcessor, Calibration::VarProcessor, const MVAComputer >	Registry

Protected Member Functions
virtual void	configure (ConfIterator iter, unsigned int n)=0
	virtual configure method, implemented in actual processor More...
virtual ConfigCtx::Context *	configureLoop (ConfigCtx::Context *ctx, ConfigCtx::iterator begin, ConfigCtx::iterator cur, ConfigCtx::iterator end)
	virtual loop configure method More...
virtual std::vector< double >	deriv (ValueIterator iter, unsigned int n) const
	virtual derivative evaluation method, implemented in actual processor More...
virtual void	eval (ValueIterator iter, unsigned int n) const =0
	virtual evaluation method, implemented in actual processor More...
	VarProcessor (const char *name, const Calibration::VarProcessor *calib, const MVAComputer *computer)

Protected Attributes
const MVAComputer *	computer

Private Attributes
BitSet	inputVars
	bit set to select the input variables to be passed to this processor More...
unsigned int	nInputVars

Additional Inherited Members
Static Public Member Functions inherited from PhysicsTools::ProcessRegistry< Base_t, CalibBase_t, Parent_t >::Factory
static Base_t *	create (const char *name, const CalibBase_t *calib, Parent_t *parent=0)

Detailed Description

Common base class for variable processors.

This base class for variable processors manages the common interface to the global discriminator class and how the individual variable processors are interconnected.

Definition at line 36 of file VarProcessor.h.

Member Typedef Documentation

typedef Dummy* PhysicsTools::VarProcessor::PluginFunctionPrototype()

Definition at line 141 of file VarProcessor.h.

typedef ProcessRegistry<VarProcessor, Calibration::VarProcessor, const MVAComputer> PhysicsTools::VarProcessor::Registry

protected

Definition at line 313 of file VarProcessor.h.

Member Enumeration Documentation

enum PhysicsTools::VarProcessor::LoopStatus

Enumerator
kStop
kNext
kReset
kSkip

Definition at line 132 of file VarProcessor.h.

{ kStop, kNext, kReset, kSkip };

Constructor & Destructor Documentation

PhysicsTools::VarProcessor::~VarProcessor ( )

virtual

Definition at line 40 of file VarProcessor.cc.

References inputVars, and nInputVars.

                              {
    inputVars = BitSet(0);
    nInputVars = 0;
  }

PhysicsTools::VarProcessor::VarProcessor ( const char *  name, const Calibration::VarProcessor *  calib, const MVAComputer *  computer  )

protected

Definition at line 37 of file VarProcessor.cc.

      : computer(computer), inputVars(Calibration::convert(calib->inputVars)), nInputVars(inputVars.bits()) {}

Member Function Documentation

void PhysicsTools::VarProcessor::configure

(

ConfigCtx &

config

)

called from the discriminator computer to configure processor

Definition at line 45 of file VarProcessor.cc.

References PhysicsTools::VarProcessor::ConfigCtx::begin(), looper::config, configureLoop(), PhysicsTools::VarProcessor::ConfigCtx::ctx, PhysicsTools::VarProcessor::ConfigCtx::end(), inputVars, PhysicsTools::BitSet::iter(), PhysicsTools::VarProcessor::ConfigCtx::loop, loop(), nInputVars, PhysicsTools::VarProcessor::ConfigCtx::size(), and PhysicsTools::BitSet::size().

Referenced by PhysicsTools::MVAComputer::setup().

                                                {
    ConfigCtx::size_type pos = config.size();
    if (pos != inputVars.size())
      return;

    ConfIterator iter(inputVars.iter(), config);
    configure(iter, nInputVars);

    VarProcessor *loop = config.loop ? config.loop : this;
    ConfigCtx::Context *ctx = loop->configureLoop(config.ctx, config.begin(), config.begin() + pos, config.end());

    if (ctx != config.ctx) {
      delete config.ctx;
      config.ctx = ctx;
    }

    if (config.loop && !ctx)
      config.loop = nullptr;
    else if (!config.loop && ctx)
      config.loop = this;
  }

virtual void PhysicsTools::VarProcessor::configure ( ConfIterator  iter, unsigned int  n  )

protectedpure virtual

virtual configure method, implemented in actual processor

VarProcessor::ConfigCtx::Context * PhysicsTools::VarProcessor::configureLoop ( ConfigCtx::Context *  ctx, ConfigCtx::iterator  begin, ConfigCtx::iterator  cur, ConfigCtx::iterator  end  )

protectedvirtual

virtual loop configure method

Definition at line 72 of file VarProcessor.cc.

Referenced by configure().

                                                                                       {
    return nullptr;
  }

void PhysicsTools::VarProcessor::deriv	(	double *	input,
		int *	conf,
		double *	output,
		int *	outConf,
		int *	loop,
		LoopCtx &	ctx,
		unsigned int	offset,
		unsigned int	in,
		unsigned int	out,
		std::vector< double > &	deriv
	)		const

run the processor evaluation pass on this processor and compute derivatives

Definition at line 103 of file VarProcessor.cc.

References PhysicsTools::VarProcessor::ConfigCtx::begin(), gather_cfg::cout, PhysicsTools::VarProcessor::ConfigCtx::ctx, relativeConstraints::empty, PhysicsTools::VarProcessor::ConfigCtx::end(), eval(), ntuplemaker::fill, mps_fire::i, recoMuon::in, input, inputVars, createfilelist::int, PhysicsTools::BitSet::iter(), dqmiolumiharvest::j, dqmdumpme::k, PhysicsTools::VarProcessor::ConfigCtx::loop, visualization-live-secondInstance_cfg::m, makeMuonMisalignmentScenario::matrix, dqmiodumpmetadata::n, nInputVars, hltrates_dqm_sourceclient-live_cfg::offset, MillePedeFileConverter_cfg::out, convertSQLitetoXML_cfg::output, AlCaHLTBitMon_ParallelJobs::p, data-class-funcs::q, PhysicsTools::VarProcessor::ConfigCtx::size(), AlCaHLTBitMon_QueryRunRegistry::string, and edm::typeDemangle().

Referenced by eval(), and PhysicsTools::MVAComputer::DerivContext::eval().

                                                           {
    ValueIterator iter(inputVars.iter(), input, conf, output, outConf, loop, ctx, offset);

    eval(iter, nInputVars);

    std::vector<double> matrix = this->deriv(iter, nInputVars);

    unsigned int size = 0;
    while (iter)
      size += (iter++).size();
    bool empty = matrix.empty();
    assert(size != 0 || empty);
    unsigned int out = empty ? 0 : (matrix.size() / size);

    if (matrix.size() != out * size || (out > 1 && (int)out != outConf[out_] - outConf[0]))
      throw cms::Exception("VarProcessor")
          << "Derivative matrix implausible size in " << typeid(*this).name() << "." << std::endl;

#ifdef DEBUG_DERIV
    if (!matrix.empty()) {
      std::string demangledName;
      edm::typeDemangle(typeid(*this).name(), demangledName);
      std::cout << demangledName << std::endl;
      for (unsigned int i = 0; i < out; i++) {
        for (unsigned int j = 0; j < size; j++)
          std::cout << matrix.at(i * size + j) << "\t";
        std::cout << std::endl;
      }
      std::cout << "----------------" << std::endl;
    }

    std::cout << "======= in = " << in << ", size = " << size << ", out = " << out << ", matrix = " << matrix.size()
              << std::endl;
#endif

    unsigned int sz = (outConf[out_] - in) * in;
    unsigned int oldSz = deriv.size();
    if (oldSz < sz)
      deriv.resize(sz);

    double *begin = &deriv.front() + (outConf[0] - in + offset) * in;
    double *end = begin + out * in;
    if (begin < &deriv.front() + oldSz)
      std::fill(begin, end, 0.0);

    if (matrix.empty())
      return;

    double *m0 = &matrix.front();
    BitSet::Iterator cur = inputVars.iter();
    for (unsigned int i = 0; i < nInputVars; i++, ++cur) {
#ifdef DEBUG_DERIV
      std::cout << " inputvar " << i << std::endl;
#endif
      int *curConf = conf + cur();
      unsigned int pos = *curConf;
#ifdef DEBUG_DERIV
      std::cout << " -> cur = " << cur() << ", pos = " << pos << std::endl;
#endif
      if (loop && curConf >= loop) {
        pos += offset;
        loop = nullptr;
      }

      unsigned int n = loop ? (curConf[1] - curConf[0]) : 1;
      for (unsigned int j = 0; j < n; m0++, j++, pos++) {
#ifdef DEBUG_DERIV
        std::cout << "  multip " << j << std::endl;
#endif
        double *p = begin;
        if (pos >= in) {
#ifdef DEBUG_DERIV
          std::cout << "   deriv " << (pos - in) << std::endl;
#endif
          const double *q = &deriv.front() + (pos - in) * in;
          for (const double *m = m0; p < end; m += size, q -= in)
            for (unsigned int k = 0; k < in; k++)
              *p++ += *m * *q++;
        } else {
#ifdef DEBUG_DERIV
          std::cout << "   in " << pos << std::endl;
#endif
          for (const double *m = m0; p < end; m += size, p += in)
            p[pos] += *m;
        }
      }
    }

#ifdef DEBUG_DERIV
    std::cout << "================" << std::endl;
    for (const double *p = &deriv.front(); p != &deriv.front() + deriv.size();) {
      for (unsigned int j = 0; j < in; j++)
        std::cout << *p++ << "\t";
      std::cout << std::endl;
    }
    std::cout << "================" << std::endl;
#endif
  }

virtual std::vector<double> PhysicsTools::VarProcessor::deriv ( ValueIterator  iter, unsigned int  n  ) const

inlineprotectedvirtual

virtual derivative evaluation method, implemented in actual processor

Definition at line 330 of file VarProcessor.h.

{ return std::vector<double>(); }

void PhysicsTools::VarProcessor::eval ( double *  input, int *  conf, double *  output, int *  outConf, int *  loop, LoopCtx &  loopCtx, unsigned int  offset  ) const

inline

run the processor evaluation pass on this processor

Definition at line 114 of file VarProcessor.h.

References deriv(), recoMuon::in, input, inputVars, PhysicsTools::BitSet::iter(), loop(), nInputVars, hltrates_dqm_sourceclient-live_cfg::offset, MillePedeFileConverter_cfg::out, and convertSQLitetoXML_cfg::output.

Referenced by deriv(), and PhysicsTools::MVAComputer::EvalContext::eval().

                                                                                                                        {
      ValueIterator iter(inputVars.iter(), input, conf, output, outConf, loop, loopCtx, offset);
      eval(iter, nInputVars);
    }

virtual void PhysicsTools::VarProcessor::eval ( ValueIterator  iter, unsigned int  n  ) const

protectedpure virtual

virtual evaluation method, implemented in actual processor

virtual LoopStatus PhysicsTools::VarProcessor::loop ( double *  output, int *  outConf, unsigned int  nOutput, LoopCtx &  ctx, unsigned int &  nOffset  ) const

inlinevirtual

Definition at line 134 of file VarProcessor.h.

References kStop.

Referenced by configure(), and eval().

                                                                                                       {
      return kStop;
    }

Member Data Documentation

const MVAComputer* PhysicsTools::VarProcessor::computer

protected

Definition at line 333 of file VarProcessor.h.

BitSet PhysicsTools::VarProcessor::inputVars

private

bit set to select the input variables to be passed to this processor

Definition at line 337 of file VarProcessor.h.

Referenced by configure(), deriv(), eval(), and ~VarProcessor().

unsigned int PhysicsTools::VarProcessor::nInputVars

private

Definition at line 338 of file VarProcessor.h.

Referenced by configure(), deriv(), eval(), and ~VarProcessor().