examples/JCalibrate/JHVGainGraph.cc

Go to the documentation of this file.

#include <string>
#include <vector>
#include <istream>

#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

#include "JROOT/JManager.hh"
#include "JGizmo/JRootObjectID.hh"

#include "JDetector/JCalibration.hh"

#include "JCalibrate/JFitToT.hh"
#include "JCalibrate/JHVGainGraph.hh"


/**
 * \file
 * 
 * Example program to interpolate input high-voltage-gain data.
 * \author bjjung
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;

  using JCALIBRATE::JHVGainGraph;

  typedef map<double, pair<double, double>>      map_t;
  typedef JRange<double>                      JRange_t;
  
  map_t             input;
  JRootObjectID     output;

  double            minDistanceHV = 2 * 3.14;
  JRange_t          rangeHV       = JRange_t(-1500 + 1e-2,
                                               -80 - 1e-2);
  JRange_t          rangeGain     = JRange_t(FITTOT_GAIN_MIN + 1e-2,
                                             FITTOT_GAIN_MAX - 1e-2);
  double            targetGain    = NOMINAL_GAIN;
  

  bool              loglog        = true;
  double            base          = 10.0;
  
  int               debug;

  try {

    JParser<> zap("Example program to interpolate high-voltage versus gain ASCII data.");

    JProperties properties;

    properties.insert(gmake_property(minDistanceHV));
    properties.insert(gmake_property(rangeHV));    
    properties.insert(gmake_property(rangeGain));
    properties.insert(gmake_property(targetGain));
    properties.insert(gmake_property(loglog));
    properties.insert(gmake_property(base));
    
    zap['f'] = make_field(input);
    zap['o'] = make_field(output)     = JPARSER::initialised();
    zap['@'] = make_field(properties) = JPARSER::initialised();
    zap['d'] = make_field(debug)      = 2;

    zap(argc, argv);
  }
  catch (const exception& error) {
    FATAL(error.what() << endl);
  }

  JHVGainGraph::setMinHVDistance(minDistanceHV);
  JHVGainGraph::setHVRange      (rangeHV);
  JHVGainGraph::setGainRange    (rangeGain);
  
  JHVGainGraph graph(output.getObjectName(), loglog, base);

  NOTICE(RIGHT(35) << "High-Voltage [V]" << RIGHT(31) << "gain [-]" << endl);

  for (map_t::const_iterator i = input.cbegin(); i != input.cend(); ++i) {

    const double HV        = i->first;
    const double gain      = i->second.first;
    const double gainError = i->second.second;
    
    NOTICE(FIXED(35,2) << HV << FIXED(22,2) << gain << " +/- " << FIXED(3,2) << gainError << endl);
    
    graph.AddPoint(HV, gain, gainError);
  }
      
  NOTICE(FILL(100,'-') << " Interpolated high-voltage" << setfill(' ') << endl);

  try {
    
    const double targetHV      = graph.getTargetHV     (targetGain);
    const double targetHVError = graph.getTargetHVError(targetGain);
  
    NOTICE(FIXED(35,2) << right << targetHV << " +/- " <<
           FIXED(8, 2) << left  << targetHVError <<
           FIXED(9, 2) << right << targetGain << endl);
  }
  catch (const exception& error) {
    FATAL(error.what() << endl);
  }

  TFile file(output.getFilename().c_str(), "RECREATE");
  
  graph.Write();
  file.Close();
}