BaseCalib.h

Go to the documentation of this file.

/***************************************************************************
 *   Copyright (C) 2005-2006 by Christophe Theisen & Olivier Stezowski     *
 *   christophe.theisen(AT)cea.fr, stezow(AT)ipnl.in2p3.fr                       *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/
 
#ifndef GW_BASECALIB_H
#define GW_BASECALIB_H

//#include<stdio>
#include<fstream>
#include<iostream>
#include<iomanip>
using namespace std;

#ifndef ROOT_TH1_H
#include "TH1.h"
#endif
#ifndef ROOT_TGRAPHERRORS_H
#include "TGraphErrors.h"
#endif
#ifndef ROOT_TF1_H
#include "TF1.h"
#endif
#ifndef GW_H1CALIBRATOR_H
#include "H1Calibrator.h"
#endif

namespace Gw {

class BaseCalib {

public:
        BaseCalib();
        virtual ~BaseCalib();

        void AlignMax(TH1 *histo, TF1 *function, const Double_t AlignValue, 
                                  const Bool_t invert=false, const Double_t xmin=0, const Double_t xmax=0); 
  
        void AlignMax(TH1 *histo, TF1 *function, const Double_t AlignValue,
                                  const Double_t ProportionalCoef, const Bool_t invert=false, const Double_t xmin=0, const Double_t xmax=0); 

        void SetPoint(const Int_t index, const Double_t channel, const Double_t energy) { fTabulated[index] = energy; fRaw[index] = channel; }
        void AddPoint(const Double_t channel, const Double_t energy);
  
    void SetEnergy(const Int_t index, const Double_t energy)   { fTabulated[index] = energy; }
    void SetChannel(const Int_t index, const Double_t channel) { fRaw[index] = channel; }
  
        Double_t GetEnergy(const Int_t index) const  { return fTabulated[index]; }
        Double_t GetChannel(const Int_t index) const { return fRaw[index]; }

        void Calibrate(TF1 *function) const;
        void Calibrate(const char* filename, TF1 *function);
        void Calibrate(TF1 *function, const Int_t nbpoints, const Double_t *raw, const Double_t *tabulated);
        void Calibrate(Int_t dimension);

        void DrawResults(const TF1 *function) const;

        static void CheckCalibration(TH1 *histo, const Double_t value, Double_t xmin=-1111, Double_t xmax=-1111, const Int_t color=2);
        static void CheckCalibration(TH1 *histo, const Int_t nbvalue, const Double_t *values, Double_t xmin=-1111, Double_t xmax=-1111, const Int_t color=2);
        static void CheckCalibration(TH1 *histo, const char *SourceName, Double_t xmin=-1111, Double_t xmax=-1111, const Int_t color=2);
protected:

        TArrayD      fTabulated;
        TArrayD      fRaw;

        H1Calibrator fCalibrator;

        ClassDef(BaseCalib,0); // a BaseCalib
        
};

}

#endif