BeamMatrix.h

Go to the documentation of this file.

// $Id: BeamMatrix.h,v 1.1 2009/12/07 18:47:46 pittam Exp $
//
// A class t
//
// R. Pittam 06/2009
#ifndef BEAMMATRIX_H
#define BEAMMATRIX_H

//#include <map>
#include <string>
#include <vector>

#include "TH1.h"
#include "TH2.h"
//#include "NCUtils/NCType.h"

//#include "TH3.h"
//#include "TCanvas.h"
//#include "TFile.h"

//#include "AnalysisNtuples/ANtpTruthInfoBeam.h"
//#include "AnalysisNtuples/ANtpRecoInfo.h"
//#include "AnalysisNtuples/ANtpAnalysisInfo.h"

#include "NCUtils/Extrapolation/NCBeam.h"
//#include "NCUtils/NCOscProb.h"

class BeamMatrix : public TNamed
{
 public:
  BeamMatrix() 
    : fBMHistsMC(0), fBMHistsData(0), fBMHists2D(0), fBMHists2DOsc(0),  fFD_BeamMatrixfit(0), fFD_data(0), mScaleNCbm(0) { };
    //,fND_XsecByParams(0), fND_XSectionFitRew(0), fND_XSectionFit(0) { };
  
  BeamMatrix(std::string bmname, std::string title, std::string name);
  //BeamMatrix(TString bmname, TString title, TString name);
  
  std::vector<TH1D*> fBMHistsMC; 
  std::vector<TH1D*> fBMHistsData ;
  std::vector<TH2D*> fBMHists2D ;
  std::vector<TH2D*> fBMHists2DOsc ;
  std::vector<TH1D*> fFD_BeamMatrixfit ;
  std::vector<TH1D*> fFD_data;

  std::vector<Double_t> mScaleNCbm;
  //std::vector<Double_t> mScaleNCErr;
 
  //double mScaleNeg[4];// positive and negative errors 
  // from MINOS
  //double mScalePos[4];
  //double corr[4];


  /*
  std::vector<std::vector<TH1D* > > fND_XsecByParams ;
  std::vector< TH1D* > fND_XSectionFitRew; 
  std::vector< TH1D* > fND_XSectionFit;   
  */
  std::vector<TH1D* > GetMCHists(){ return fBMHistsMC ;}
  std::vector<TH1D* > GetDataPredHists(){ return fBMHistsData ;}
  std::vector<TH2D* > Get2DHists(){ return fBMHists2D ;}
  std::vector<TH2D* > Get2DOscHists(){ return fBMHists2DOsc ;}
  std::vector<TH1D* > GetFitHists(){ return fFD_BeamMatrixfit ;}
  std::vector<TH1D* > GetFDDataHists(){ return fFD_data ;}

  std::vector<Double_t> GetmScaleNC(){ return mScaleNCbm ;}
  void SetmScaleNCparam(int param , double val);//{ 
    // MSG("BeamMatrix",Msg::kInfo) << "param " << param << " val " << val << " mScaleNCbm.at("<<param<<") " << mScaleNCbm.at(param) <<endl; 
    //mScaleNCbm.at(param) = val ; 
    //MSG("BeamMatrix",Msg::kInfo) << "param " << param << " val " << val << " mScaleNCbm.at("<<param<<") " << mScaleNCbm.at(param) <<endl; 
    //return;}
  //std::vector<Double_t> GetmScaleNC(){ return mScaleNC ;}
  //std::vector<std::vector<TH1D* > > GetNDXSecByParams {return fND_XsecByParams; } 
  //std::vector<TH1D* > GetNDXSecFitRewHists(){ return fND_XSectionFitRew ;}
  //std::vector<TH1D* > GetNDXsecFitHists(){ return fND_XSectionFit ;}
  
 private:

  double fTrue_bin_width;

  ClassDef(BeamMatrix, 1);


};
#endif 

---