CompSpecAna.h

Go to the documentation of this file.

#ifndef _COMPSPEC_ANA_
#define _COMPSPEC_ANA_
//project 
#include "PostProcBase.h"
#include "CompSpecEvent.h"
#include "CutManager.h"
#include "CutBuilder.h"
#include "ConfigBuilder.h"
#include "HistHelpers.h"
#include "COMPASSWFAna.h"
//root
#include "TH1.h"
#include "TH2.h"
#include "TF1.h"
#include "TGraphErrors.h"
#include "TFitResult.h"
//c++
#include <string>
#include <map>
#include <random>
#include <fstream>
#include <TStyle.h>
#include <cmath>
#include <nlohmann/json.hpp>
 
class CompSpecAna : public PostProcBase<CompSpecEvent>
{
public:
  //default constructor
  CompSpecAna();
 
  CompSpecAna(std::string a_jsonConfigName);
 
  //used to read the config after the class has already been instantiated.
  void loadConfig(std::string a_jsonConfigName);
 
  // used to build and store histogram of target scintillator 
  // light yield vs electron energy deposition in target scintillator
  void buildTargetLightYield(bool a_applyCuts = true);
 
  void buildRelativeLightYield(std::string a_outfile="relativeLightYield.root",
                               bool a_applyCuts = true,
                               bool a_debug = true);
 
  void plotWFs(int a_numWFs = 10);
 
  std::map<string, bool> m_cutState;
 
private:
  ExperimentConfig m_targetConfig;
  ExperimentConfig m_scatterConfig;
  PointSourceInfo m_sourceConfig;
  nlohmann::json m_jsonConfig;
 
  CutManager m_cutManager;
 
  bool passesCuts(CompSpecEvent* a_event);
  bool m_isSim;
 
  //these are used to functionalize the templated base class
  void setTreeName();
  void setBranchName();
 
  //used to store 2D hists from relative light yield procedure
  std::map<int, TH2*> m_hists;
 
  // used to store relative light yield curve
  TGraphErrors m_relativeLightYield;
};
 
 
#endif