LightYield.hh

Go to the documentation of this file.

#ifndef LIGHT_YIELD_h
#define LIGHT_YIELD_h

#include "Rtypes.h"
#include <map>
#include <vector>
#include "TH1.h"
class TGraphErrors;
class TTree;


class LightYieldGraph{
  //these stick around
  std::vector<double> _x, _y, _ex, _ey; 
  std::map<int, std::vector<double> > _chan_spe; 
  std::map<int, std::vector<double> > _chan_spe_err; 
  TH1* histo; 
public:
  LightYieldGraph() : histo(0) {} 
  virtual ~LightYieldGraph(){ if(histo) delete histo; }
  int AddRun(int run, double epeak, int nbins=200, 
             double emin=1000, double emax=2500);
  int AddRuns(const std::vector<int>& runs, double epeak, int nbins=200, 
              double emin=1000, double emax=2500);

  int AddFile(const char* filename, double epeak, int nbins=200, 
             double emin=1000, double emax=2500);
  int SetAliasesFromLocalData(TTree* Events, bool draw = false);
  
  TH1* DrawSpectrum(const char* filename, int nbins = 200, double xmin=0, 
                    double xmax = 4500);
  TGraphErrors* DrawGraph();
  
  TGraphErrors* DrawChannelSpeGraph(int channel);
  
  static int FitPhotopeak(TH1* h, double epeak, 
                          double& yield, double& yield_err);

  ClassDef(LightYieldGraph,1);
};

TGraphErrors* PlotNa22LightYield();
#endif