/data/refman/pasoursint/CMSSW_5_3_10/src/DataFormats/ParticleFlowReco/interface/CalibrationResultWrapper.h

Go to the documentation of this file.

#ifndef CALIBRATIONRESULTWRAPPER_H_
#define CALIBRATIONRESULTWRAPPER_H_

//#include <boost/shared_ptr.hpp>

#include "DataFormats/ParticleFlowReco/interface/CalibrationProvenance.h"

namespace pftools {

class CalibrationResultWrapper {
public:

//      typedef boost::shared_ptr<CalibrationResultWrapper>
//                      CalibrationResultWrapperPtr;

        CalibrationResultWrapper() {
                reset();
        }

        virtual ~CalibrationResultWrapper() {
        }


        void reset() {
                resetCore();
        }


        void compute() {
                computeCore();
        }

        double bias() const {
                return (particleEnergy_ -  truthEnergy_) / truthEnergy_;
        }

        double ratio() const {
                return(particleEnergy_/truthEnergy_);
        }

        /*
         * Which calibrator made this?
         */
        CalibrationProvenance provenance_;
        /*
         * What energy was this particle optimised to?
         */
        double truthEnergy_;

        /*
         * Calibrated ecal deposition
         */
        double ecalEnergy_;

        /*
         * Calibrated hcal deposition
         */
        double hcalEnergy_;

        /*
         * Calibrated particle energy (not necessarily ecal + hcal!)
         */
        double particleEnergy_;

        /*
         * What objects did this optimise on?
         */
        CalibrationTarget target_;

        /*
         * (reco - truth)/truth
         */
        double bias_;

        /*
         * reco/truth
         */
        double ratio_;

        /*
        * Target function contribution
        */
        double targetFuncContrib_;

        double a_;
        double b_;
        double c_;

private:

        virtual void computeCore() {
                bias_ = bias();
                ratio_ = ratio();
        }

        virtual void resetCore() {
                truthEnergy_ = 0;
                ecalEnergy_ = 0;
                hcalEnergy_ = 0;
                particleEnergy_ = 0;
                provenance_ = UNCALIBRATED;
                target_ = UNDEFINED;
                bias_ = 0;
                ratio_ = 1.0;
                targetFuncContrib_ = 0;
                a_ = 0.0;
                b_ = 1.0;
                c_ = 1.0;
        }

};

}

#endif /*CALIBRATIONRESULTWRAPPER_H_*/