NeutronAngularDistAna.h

Go to the documentation of this file.

//ROOT includes
#include "TH1.h"
#include "TH2.h"
#include "TMatrixD.h"
#include "TFitResult.h"
//c++includes
#include <cstdint>
#include <vector>
#include <utility>
#include <deque>
#include <map>
using std::vector;
using std::map;
using std::pair;
// externals
// #include "../../BasicSupport/external/eigen/Eigen/Dense"
#include <Eigen/Dense>
 
struct EffData
{
    EffData();
    vector<double> m_lowBinEdges;
    vector<double> m_effs;
    vector<double> m_effsErr;
    double m_angle;
    vector<double> getEfficiency(double a_lowE, double a_highE);
    void drawEff();
};
class NeutronAngularDistAna
{
public:
    NeutronAngularDistAna();
    void readEfficiencyData(std::string a_effFile);
    void uniBinEfficiency(int a_detID, int a_nBins, double a_lowE=1);
    TH2F* drawEffCorrAngDist(std::string a_expFile,
                                                    std::string a_bkgFile);
    void addExperimentalTOFData(std::string a_fileName);
    // stores scint/scint in m_tofHistBKGOrgOrg
    void addExperimentalTOFBKGData(std::string a_fileName);
    void setBackgroundScaler(double a_bkgScaler,
                                                        bool a_scintScint=false);
    std::vector<double> getEnerDistTOF(int a_detID,
                                                                        double a_tof,
                                                                        bool a_draw=false,
                                                                        bool a_useEff=true,
                                                                        double a_tofWidth=8);
    std::vector<double> getEnerDistTOF(double a_angle,
                                                                        double a_angleWidth,
                                                                        double a_tof,
                                                                        bool a_draw=false,
                                                                        bool a_useEff=true,
                                                                        double a_tofWidth = 8);
    void drawRFExpBkg(int a_detID, bool a_drawSub=false);
 
    void plotEnerDistTOF(int a_detID);
    void plotEnerDistTOF(double a_angle,
                                            double a_angleWidth);
    map<double, std::pair<double,double>> plotEnerDistTOF(int a_detID,
                                                Eigen::VectorXf a_calcYield,
                                                double a_normalization = -1);
    void plot2DEnerDistTOF(int a_detID);
    void getAngDistTOF(int a_tof, 
                                            bool a_normalize=false,
                                            int a_tofWidth=8);
    EffData getEffData(int a_detID);
    TH2F* getScintScintHist(int a_detID,
                                                    bool a_normEff=false);
    private:
    map<int, EffData> m_detEff;
    map<int, TH2*> m_tofHists;
    map<int, TH2*> m_tofHistsBKG;
    map<int, TH2*> m_tofHistsOrgOrg;
    map<int, TH2*> m_tofHistsBKGOrgOrg;
 
    map<int,map<double,vector<double>>> m_maxwellFitRes;
    double m_bkgScaler;
    double m_bkgScalerOrgOrg;
};
 
// [0]*2*sqrt(x/3.14)*pow(1./[1],1.5)*exp(-x/[1])+[2]