ThinTargetMC.h

Go to the documentation of this file.

#ifndef THINTARGETMC_H
#define THINTARGETMC_H

#include "TH2D.h"
#include "TFile.h"
#include <vector>
#include <string>

namespace NeutrinoFluxReweight{
  
  /*! \class ThinTargetMC
   *  \brief A class to manage the MC value for thin target 
   */
  class ThinTargetMC{

  private:
    ThinTargetMC();

  public:
    static ThinTargetMC* getInstance();
    
    //! MC value for this HP production
    double getMCval_pC_X(double incP, double xf,double pt, int pdgcode);

    //! Get the MC roduction cross-section pC->pi, K:
    double getMCxs_pC_piK(int genid, double inc_mom);
    //! Get the MC roduction cross-section pC->n, p:
    double getMCxs_pC_nucleon(int genid, int pdg, double inc_mom);

    //! Vector of the scaling histograms:
    std::vector< std::vector<TH2F*> > hTTScl; 
    //! Vector of the scaling histograms for neutrons:
    std::vector<TH1F*> hTTScl_n; 
    
  private:

    TFile* fpC_x[8];     //convention: 0: pip. 1: pim.  2: kap. 3: kam. 4: prt 5: neu. 6: K0L. 7:K0S
    std::vector<TH2D*> vpC_x[8];
    std::vector<TH1D*> vpC_n;
    
    std::vector<std::string> spart_prod;
    std::vector<int> mom_inc;
    
    //convention: 0: prt. 1: neu. 
    std::vector<std::string> spart_qe_corr;
    TFile* fqe_corr[2];
    std::vector<TH2D*> vqe_corr_p;
    std::vector<TH1D*> vqe_corr_n;

    std::vector<TFile*> fTTscale;
    
    static ThinTargetMC* instance;
    
  };

  
  
}
#endif