LinearCalibrator.h

Go to the documentation of this file.

#ifndef LINEARCALIBRATOR_HH_
#define LINEARCALIBRATOR_HH_

#include "RecoParticleFlow/PFClusterTools/interface/Calibrator.h"


#include "TMatrixD.h"
#include "TVectorD.h"

namespace pftools {
class LinearCalibrator : public Calibrator {
public:
    LinearCalibrator();
    virtual ~LinearCalibrator();


    /*
     * Note: covariant return type w.r.t. Calibrator class: the overloading has changed 
     * the return type but this IS allowed with modern compilers.
     * See documentation in Calibrator.h
     */
    LinearCalibrator* clone() const;
    LinearCalibrator* create() const;

protected:
    
    virtual std::map<DetectorElementPtr, double>
            getCalibrationCoefficientsCore() throw(PFToolsException&);


    LinearCalibrator(const LinearCalibrator& lc);
    /*
     * Converts the particle deposits into a useful matrix formulation.
     */
    virtual void initEijMatrix(TMatrixD& eij, TVectorD& truthE);

    /*
     * Utility method to extract the unique number of detected elements.
     */
    virtual void populateDetElIndex();

    virtual TVectorD& getProjections(const TMatrixD& eij, TVectorD& proj,
            const TVectorD& truthE) const;

    virtual TMatrixD& getHessian(const TMatrixD& eij, TMatrixD& hess,
            const TVectorD& truthE) const;

    /*
     * Map to convert detector element to array row/column index.
     */
    std::map<DetectorElementPtr, unsigned> myDetElIndex;

};
}

#endif /*LINEARCALIBRATOR_HH_*/