SpectrumAnalysis.h

Go to the documentation of this file.

#ifndef _SPECTRUM_ANALYSIS_H_
#define _SPECTRUM_ANALYSIS_H_
 
#include "TH1.h"
#include "TFile.h"
 
#include <vector>
#include <map>
#include <string>
using std::map;
using std::vector;
using std::string;
 
struct Peak
{
  Peak();
  double m_amplitude  = 0.;
  double m_mean       = 0.;
  double m_stdDev     = 0.;
  double m_chi2ndf    = 0.;
  double m_binWidth   = 0.;
  double m_stdDevErr  = 0.;
  double m_ampErr     = 0.;
  void printPeak();
  bool checkResolution();
  double m_maxRes     = 2.; // percent deltaE/E
};
 
class SpectrumAnalysis
{
public:
  SpectrumAnalysis();
 
  void loadHists(string a_fname, int a_modID, double a_countTime=0);
  void loadSimHists(string a_fname);
  void loadBackgroundHists(string a_fname, int a_modID, double a_countTime);
  void subtractBackgroundSpectra();
  vector<Peak> FitPeaksMan(int a_channel, 
                          double a_lower, 
                          double a_higher,// int a_numPeaks,
                          vector<double> a_peakMeans,
                          bool a_fullClover=false, 
                          bool a_addBack=false, 
                          bool a_checkPeakRes=true);
 
  vector<Peak> FitPeaksManPol1(int a_channel, 
                              int a_lower, 
                              int a_higher, 
                              // int a_numPeaks, 
                              vector<double> a_peakMeans,
                              bool a_checkPeakRes=true,
                              bool a_fullClover=false, 
                              bool a_addBack=false);
  void printPeakFits();
  vector<Peak> FitPeaksErf(int a_channel, 
                          double a_lower, 
                          double a_higher, 
                          int a_numPeaks,
                          bool a_fullClover=false, 
                          bool a_addBack=false);
 
  map<double, std::pair<double,double>> Eu152Efficiency(int a_channel, 
                                                        double a_activity, 
                                                        double a_countTime,
                                                        double a_simBrTot=1,
                                                        bool a_res=false,
                                                        bool a_fullClover=false,
                                                        bool a_addBack=false,
                                                        bool a_draw =false);
 
  void calcAddBackFactor();
  map<double, std::pair<double,double>> U235Efficiency(int a_channel, 
                                                      double a_activity,
                                                      double a_activityErr, 
                                                      double a_countTime);
 
  map<int, TH1*> getHists();
  map<int,vector<TH1*>> getCloverHists();
  map<int, vector<Peak>> getPeakFits();
 
  void ChannelResolution(int a_channel);
/* 
  void findPeaks(int a_channel, double a_minHeight, int a_numPeaks);
 
  void FitPeakBasic(int a_channel, double a_peakCentroid, double a_low= 5 , double a_high=5);;
  void FitDoublePeak(int a_channel, double a_peakCentroid, double a_peakSeparation, 
                                      double a_relativeHeight, double a_fitWindow);
 
  map<string, vector<vector<double>>> m_zeroCrossings;
  void writeToConfig();
*/
 
 
private:
  double m_countTime=0;
  double m_bkgCountTime=0;
  vector<double> m_energy;
 
  map<int, TH1*> m_hists;
  map<int, TH1*> m_bkgHists;
  map<int, vector<TH1*>> m_cloverHists;
  map<int, vector<double>> m_calibration;
  
  map<int, vector<Peak>> m_peakFits;
};
#endif