Calibrator.h

Go to the documentation of this file.

// $Id: Calibrator.h,v 1.6 2009/10/19 17:15:21 jpaley Exp $
//
// Singleton class that handles calibrations for a particular 
// validity context (eg: detector+timestamp).
//

#ifndef CALIBRATOR_H
#define CALIBRATOR_H

#include "TObject.h"
#include "TTimeStamp.h"
#include "TH1D.h"
#include <iostream>
#include <vector>
#include <map>
#include <string>

#include "Geometry/Geometry.h"
#include "RawData/DAQHeader.h"
#include "RawData/RawDigit.h"
#include "RecoBase/CellHit.h"
#include "RecoBase/RecoHit.h"

using namespace rawdata;
using namespace recobase;

namespace calib {

  struct attlen {
    float a1, l1;
    float a2, l2;
    float a1sig, l1sig;
    float a2sig, l2sig;
  };
  
  class Calibrator {

  public:
    Calibrator(); 

    // Get a Calibrator instance here (or use EDMUtils::GetCalibrator())
    static Calibrator& Instance(TTimeStamp timestamp,
                                rawdata::det_id_ detector,
                                const geo::Geometry &geometry);

    // Get the timestamp or detector ID currently assumed by the Calibrator
    TTimeStamp       GetTimeStamp() const { return fCurrTimeStamp; }
    rawdata::det_id_ GetDetector () const { return fCurrDetector;  }


  private:
    // private methods and data members for Calibrator class
    void SetContext(TTimeStamp timestamp,
                    rawdata::det_id_ detector,
                    const geo::Geometry &geometry);    
    static Calibrator*   fInstance;
    TTimeStamp           fCurrTimeStamp;
    rawdata::det_id_     fCurrDetector;
    const geo::Geometry* fCurrGeo;

  //---------------------------------------------
  // RawDigit to CellHit methods
  //---------------------------------------------
  public:
    // Main thing you'll want:
    void  MakeCellHit    (const RawDigit *rawdigit, CellHit *cellhit); // Do everything

    float GetPE            (const RawDigit *rawdigit);
    float GetPECorrFactor  (const RawDigit *rawdigit);
    float GetTNS           (const RawDigit *rawdigit);

  //---------------------------------------------
  // CellHit to RecoHit methods
  //---------------------------------------------
  // (NOTE: w != dist_to_readout.  Those are
  // different arguments.)
  //---------------------------------------------
  public:
    // Main thing you'll want:

    // Do everything
    void  MakeRecoHit(const CellHit *cellhit, RecoHit *recohit, float w); 

    float GetT       (const CellHit *cellhit, float dist_to_readout);
    float GetPEAtt   (const CellHit *cellhit, float dist_to_readout);
    float GetMIPScale(const CellHit *cellhit); // for PEAtt-->MIP conversion
    float GetGeVScale(const CellHit *cellhit); // for MIP-->GeV conversion
    
    float                fPEACperMIP;
    float                fPEACperMIPErr;
    std::vector<float>   fMIPtoGeV;
    std::vector<float>   fMIPtoGeVErr;
    std::map<int,float>  fADCtoPE;
    std::map<int,float>  fPEtoSigCor;
    std::map<int,attlen> fSCtoPEAC;
    std::vector<float>   fFTDistNDX;
    std::vector<float>   fFTDistNDY;
    std::vector<float>   fFTDistFDX;
    std::vector<float>   fFTDistFDY;

  private:
    bool fLoadedFineTimeDistros;

    // Internal utilites for above method(s).  If you need this
    // functionality, look to Geometry, not Calibrator.
    void GetXYZD(const CellHit *cellhit, float w, double *xyzd);

    bool LoadADCtoPE();
    bool LoadPEtoSigCor();
    bool LoadSCtoPEAC();
    bool LoadPEACtoMIP();
    bool LoadMIPtoGeV();
    bool LoadFineTimeDistros();

  };
}
#endif
  

---