TimePixCalibrationProcessor.h

#ifndef TIMEPIXCALIBRATIONPROCESSOR_H
#define TIMEPIXCALIBRATIONPROCESSOR_H

// C++
#include <string>
#include <fstream>

// LCIO
#include "lcio.h"

// Marlin
#include "marlin/Processor.h"
#include "TimePixPixelMode.h"
#include "TimePixPixelInterpolationProcessor.h"
#include "TimePixCalibParam.h"
#include "lccd/ConditionsMap.hh"

//AIDA
#include <AIDA/AIDA.h>

namespace AIDA
{
class IHistogram2D;
}

namespace marlintpc
{

/* the following comment is tranfered to doxygen documentation
 * which begins with a second asterix, or in case of on line comment an extra third slash
 */

class TimePixCalibrationProcessor : public marlin::Processor
{
    public:

        /* This method will be called by the marlin package
         * It returns a processor of the currend type (here TimePixCalibrationProcessor)
         */
        Processor* newProcessor()
        {
            return new TimePixCalibrationProcessor;
        }

        /* the default constructor
         * here the processor parameters are registered to the marlin package
         * other initialisation should be placed in the init method
         */
        TimePixCalibrationProcessor();

        /* Called at the begin of the job before anything is read.
         * Use to initialize the processor, e.g. book histograms
         */
        void init();

        /* Called for every run.
         * in this processor it is not used
         */
        void processRunHeader(lcio::LCRunHeader* run);

        /* Called for every event - the working horse.
         * Here the real conversion take place
         */
        void processEvent(lcio::LCEvent* evt);

        /* Called after data processing for clean up.
         * e.g. saving contol histograms, delete created objects
         * (which are not stored in a collection/event), etc.
         */
        void end();


    protected:

        /* normally the doxygen documentation is placed before the method/variable
         * but for one line comments in case of variables is may be more convenient to you
         * to place it after the variable. In this case you use an arrow to tell doxygen
         * to which variable the comment belongs to.
         */

        /* For input and output maps: 4 histograms are created in the processor for each of the chips,
         * to compare the pixel map before and after the calibration.
         * (for single chips only one of the histograms will be full, while the others are empty)
         */
        AIDA::IHistogram2D * _inputMapHisto0;
        AIDA::IHistogram2D * _outputMapHisto0;
        AIDA::IHistogram2D * _inputMapHisto1;
        AIDA::IHistogram2D * _outputMapHisto1;
        AIDA::IHistogram2D * _inputMapHisto2;
        AIDA::IHistogram2D * _outputMapHisto2;
        AIDA::IHistogram2D * _inputMapHisto3;
        AIDA::IHistogram2D * _outputMapHisto3;

        std::string _inputColName;  
        std::string _inputCalibrationData;  
        std::string _outputColName; 

        std::string _dataFileList;  
        std::ifstream _infile;    

        int _nColumn;  
        int _nRow;  
        int _nParameters;  
        int _minNoPixelToInterpolate;  
        double _vTHL;  

        bool _outputIsTransient;

        int _nRun;  
        int _nEvt;  

        void createInputPixelMap(LCCollection* inCol);

        float calibrateTimePixel(int chipID, int pixelID, float pixelValue);

        float calibrateChargePixel(int chipID, float pixelValue);

        float interpolatePixelCharge(int chipID, int pixelID);

        tpcconddata::TimePixMap* map;    
        std::map< tpcconddata::TimePixPixelMode::PixelKey, float> inMap;


        /* Called to read the calibration parameters by their index (param)
         *
         */
        float getParameter(int chipID, int param);

        bool _isFirstEvent;
        LCCollection* _vCalib;

};
} // class declaration
#endif