nbw::skzpReweight Class Reference

#include "skzpReweight.h"

Classes
struct	LessThan
struct	mapkey

Public Member Functions
	skzpReweight (std::string fpath="/nova/data/flux/SKZPdata/fluka05ptxf.root", std::string bpath="/nova/data/flux/SKZPdata/IPNDhists.root", int flag=2)
	~skzpReweight ()
void	SetParams (std::vector< double > fpar, std::vector< double > bpar)
void	SetFlukParams (std::vector< double > fpar)
void	SetBeamParams (std::vector< double > bpar)
double	GetWeight (const simb::MCFlux *mcf, double Enu, int det, int beam)
double	GetFlukWeight (const simb::MCFlux *mcf)
double	GetFlukWeight (int ptype, double pT, double pz)
double	GetBeamWeight (int ntype, double Enu, int det=1, int beam=2)

Private Types
typedef std::map< double, double >	WeightMap_t

Private Member Functions
void	FlukConfig ()
Conventions::ParticleType_t	GeantToEnum (int ptype)
std::string	PartEnumToString (Conventions::ParticleType_t ptype)
void	BeamConfig ()
std::string	GetHname (int inu, int eff, int beam, int det)
void	FillVector (TH1D *hist, int ntype, int eff, int beam, int det)
void	FillVector (TH1F *hist, int ntype, int eff, int beam, int det)
std::string	BeamSysToString (Conventions::BeamSys_t bstype)
std::string	BeamTypeToString (Conventions::BeamType_t btype)
std::string	DetTypeToString (Conventions::DetType_t dtype)

Private Attributes
std::vector< double >	fFPar
std::string	fFpath
std::vector< Conventions::ParticleType_t >	fPlist
std::map< Conventions::ParticleType_t, TH2F * >	fPTPZ
std::map< Conventions::ParticleType_t, TH2F * >	fWeightedPTPZ
std::map< Conventions::ParticleType_t, TH2F * >	fWeightHist
std::map< Conventions::ParticleType_t, double >	fMeanPT
std::map< Conventions::ParticleType_t, double >	fN
std::map< Conventions::ParticleType_t, double >	fNWeighted
std::map< Conventions::ParticleType_t, double >	fWeightedMeanPT
std::map< Conventions::ParticleType_t, int >	fNBinsY
std::map< Conventions::ParticleType_t, int >	fNBinsX
std::vector< double >	fBPar
std::string	fBpath
TFile *	fBeamSysFile
std::map< mapkey, WeightMap_t, LessThan >	fBeamSysMap
int	fBflag

Detailed Description

Definition at line 23 of file skzpReweight.h.

Member Typedef Documentation

typedef std::map<double, double> nbw::skzpReweight::WeightMap_t

private

Definition at line 119 of file skzpReweight.h.

Constructor & Destructor Documentation

nbw::skzpReweight::skzpReweight	(	std::string	fpath = "/nova/data/flux/SKZPdata/fluka05ptxf.root",
		std::string	bpath = "/nova/data/flux/SKZPdata/IPNDhists.root",
		int	flag = 2
	)

Definition at line 18 of file skzpReweight.cxx.

References BeamConfig(), fBflag, fBPar, fBpath, fFPar, fFpath, and FlukConfig().

  {
    //default parameters specified by minos-doc-7146
    fFPar.push_back(1.56);
    fFPar.push_back(-6.42);
    fFPar.push_back(1.11);
    fFPar.push_back(0.13);
    fFPar.push_back(1.00);
    fFPar.push_back(1.25);
    fFPar.push_back(3.50);
    fFPar.push_back(4.83);
    fFPar.push_back(1.51);
    fFPar.push_back(0.29);
    fFPar.push_back(0.97);
    fFPar.push_back(2.16);
    fFPar.push_back(1.04);
    fFPar.push_back(-0.89);
    fFPar.push_back(0.88);
    fFPar.push_back(0.05);

    fBPar.push_back(-3.85);
    fBPar.push_back(1.39);

    fFpath = fpath;
    fBpath = bpath;
    fBflag = flag;
    FlukConfig();
    if (fBflag>0)
      BeamConfig();
  }

nbw::skzpReweight::~skzpReweight

(

)

Definition at line 50 of file skzpReweight.cxx.

References fFPar, fMeanPT, fN, fNBinsX, fNBinsY, fNWeighted, fPlist, fPTPZ, fWeightedMeanPT, fWeightedPTPZ, and fWeightHist.

  {
    //Deconstructs members for Fluk
    for (std::vector<Conventions::ParticleType_t>::iterator itPlist=fPlist.begin();itPlist!=fPlist.end(); itPlist++)
    {
      delete fWeightHist[*itPlist];
      delete fPTPZ[*itPlist];
      delete fWeightedPTPZ[*itPlist];
    }
    fWeightHist.clear();
    fPTPZ.clear();
    fWeightedPTPZ.clear();
    fMeanPT.clear();
    fWeightedMeanPT.clear();
    fN.clear();
    fNWeighted.clear();
    fPlist.clear();
    fFPar.clear();
    fNBinsX.clear();
    fNBinsY.clear();
  }

Member Function Documentation

void nbw::skzpReweight::BeamConfig ( )

private

Definition at line 254 of file skzpReweight.cxx.

References fBeamSysFile, fBflag, fBpath, FillVector(), GetHname(), hist, hist2, nbw::Conventions::kBeamEnd, nbw::Conventions::kBeamSysEnd, nbw::Conventions::kDetEnd, nbw::Conventions::kMINOSrat, nbw::Conventions::kNOvArat, and nbw::Conventions::kUnknownDet.

Referenced by SetBeamParams(), SetParams(), and skzpReweight().

  {
    std::cout <<"BeamConfig reading data from: "<<fBpath<<std::endl;
  
    bool FoundHist = false;
    TDirectory *save = gDirectory;
    fBeamSysFile = new TFile(fBpath.c_str());
    if (fBeamSysFile->IsZombie()){
      std::cout << "Don't recognize path: " << fBpath << std::endl;
      return;
    }

    int ntype[]    ={    56,       55  , 53  , 52};
    for (int inu = 0; inu < 4; ++inu) {
      //'End' enums are there so one can just change conventions when needed.
      for (int eff=1;eff<Conventions::kBeamSysEnd;eff++) {
        for (int beam=1;beam<Conventions::kBeamEnd;beam++) {
          for (int det=1;det<Conventions::kDetEnd;det++) {
            std::string hname = GetHname(inu,eff,beam,det);
            //I know this is ugly, but it works and I don't know why the dynamic_cast does not work like it is supposed to. I'll fix it later.
            TH1D* hist=dynamic_cast<TH1D*> (fBeamSysFile->Get(hname.c_str()));
            if (hist)
            {
              FoundHist = true;
              if (det!=Conventions::kMINOSrat || det!=Conventions::kNOvArat)
                FillVector(hist,ntype[inu],eff,beam,det);
              else //for far/near error bar internally use kUnknownDet
                FillVector(hist,ntype[inu],eff,beam,Conventions::kUnknownDet);
            }
            TH1F* hist2=dynamic_cast<TH1F*> (fBeamSysFile->Get(hname.c_str()));
            if (hist2)
            {
              FoundHist = true;
              if (det!=Conventions::kMINOSrat || det!=Conventions::kNOvArat)
                FillVector(hist2,ntype[inu],eff,beam,det);
              else //for far/near error bar internally use kUnknownDet
                FillVector(hist2,ntype[inu],eff,beam,Conventions::kUnknownDet);
            }
          }
        }
      }
    }
    if (!FoundHist)
      std::cout<<"No histograms found in " << fBpath <<  " with name formatted as expected by flag: " << fBflag << std::endl;
    fBeamSysFile->Close();
    delete fBeamSysFile;
    fBeamSysFile = 0;
    save->cd();
    return;
  }

std::string nbw::skzpReweight::BeamSysToString ( Conventions::BeamSys_t  bstype )

private

Definition at line 413 of file skzpReweight.cxx.

References fBflag, nbw::Conventions::kHornIDist, nbw::Conventions::kHornIMiscal, and nbw::Conventions::kUnknownEff.

Referenced by GetHname().

  {
    if(fBflag>=0 || fBflag<=2)
    {
      switch (bstype)
        {
        case Conventions::kHornIMiscal: return "HornIMiscal"; break;
        case Conventions::kHornIDist  : return "HornIDist";   break;
        case Conventions::kUnknownEff : return "Unknown";     break;
        default                       : return "Unknown";     break;
        }
    }
    return "Unknown";
  }

std::string nbw::skzpReweight::BeamTypeToString ( Conventions::BeamType_t  btype )

private

Definition at line 428 of file skzpReweight.cxx.

References fBflag, nbw::Conventions::kLE, nbw::Conventions::kLE010z000i, nbw::Conventions::kLE010z170i, nbw::Conventions::kLE010z185i, nbw::Conventions::kLE010z185iL, nbw::Conventions::kLE010z200i, nbw::Conventions::kLE100z200i, nbw::Conventions::kLE150z200i, nbw::Conventions::kLE250z200i, and nbw::Conventions::kUnknownBeam.

Referenced by GetHname().

  {
    if(fBflag == 0)
    {
      switch (btype)
        {
        case Conventions::kLE         : return "LE";          break;
        case Conventions::kLE010z185i : return "LE010z185i";  break;
        case Conventions::kLE100z200i : return "LE100z200i";  break;
        case Conventions::kLE250z200i : return "LE250z200i";  break;
        case Conventions::kLE010z185iL: return "LE010z185iL"; break;
        case Conventions::kLE010z170i : return "LE010z170i";  break;
        case Conventions::kLE010z200i : return "LE010z200i";  break;
        case Conventions::kLE010z000i : return "LE010z000i";  break;
        case Conventions::kLE150z200i : return "LE150z200i";  break;
        case Conventions::kUnknownBeam: return "Unknown";     break;
        default          : return "Unknown";     break;
        }
    }
    else if (fBflag == 1 || fBflag == 2)
    {
      switch (btype)
        {
        case Conventions::kLE         : return "L";                break;
        case Conventions::kLE010z185i : return "L010z185i";        break;
        case Conventions::kLE100z200i : return "L100z200i";        break;
        case Conventions::kLE250z200i : return "L250z200i";        break;
        case Conventions::kLE010z185iL: return "L010z185i_lowint"; break;
        case Conventions::kLE010z170i : return "L010z170i";        break;
        case Conventions::kLE010z200i : return "L010z200i";        break;
        case Conventions::kLE010z000i : return "L010z000i";        break;
        case Conventions::kLE150z200i : return "L150z200i";        break;
        case Conventions::kUnknownBeam: return "Unknown";          break;
        default                       : return "Unknown";          break;
        }
    }
    return "Unknown";
  }

std::string nbw::skzpReweight::DetTypeToString ( Conventions::DetType_t  dtype )

private

Definition at line 467 of file skzpReweight.cxx.

References fBflag, nbw::Conventions::kIPND, nbw::Conventions::kMINOSfd, nbw::Conventions::kMINOSnd, nbw::Conventions::kMINOSrat, nbw::Conventions::kNOvAfd, nbw::Conventions::kNOvAnd, nbw::Conventions::kNOvArat, and nbw::Conventions::kUnknownDet.

Referenced by GetHname().

  {
    if(fBflag == 1)
    {
      switch (dtype)
        {
        case Conventions::kMINOSnd   : return "Near";    break;
        case Conventions::kMINOSfd   : return "Far";     break;
        case Conventions::kMINOSrat  : return "FN";      break;
        case Conventions::kUnknownDet: return "Unknown"; break;
        default                      : return "Unknown"; break;
        }
    }
    else if (fBflag == 0 || fBflag == 2)
    {
      switch (dtype)
        {
        case Conventions::kNOvAnd    : return "NOvAnd";   break;
        case Conventions::kNOvAfd    : return "NOvAfd";   break;
        case Conventions::kIPND      : return "IPND";     break;
        case Conventions::kMINOSnd   : return "MINOSnd";  break;
        case Conventions::kMINOSfd   : return "MINOSfd";  break;
        case Conventions::kNOvArat   : return "NOvArat";  break;
        case Conventions::kMINOSrat  : return "MINOSrat"; break;
        case Conventions::kUnknownDet: return "Unknown";  break;
        default                      : return "Unknown";  break;
        }
    }
    return "Unknown";
  }

void nbw::skzpReweight::FillVector ( TH1D *  hist, int  ntype, int  eff, int  beam, int  det  )

private

Definition at line 323 of file skzpReweight.cxx.

References nbw::skzpReweight::mapkey::BeamDex, nbw::skzpReweight::mapkey::DetDex, nbw::skzpReweight::mapkey::EffDex, fBeamSysMap, and nbw::skzpReweight::mapkey::NuDex.

Referenced by BeamConfig().

  { 
    mapkey dex;
    dex.NuDex=ntype;
    dex.EffDex=eff;
    dex.BeamDex=beam;
    dex.DetDex=det;

    std::map<mapkey, WeightMap_t, LessThan>::iterator dexit = fBeamSysMap.find(dex);
    if (dexit == fBeamSysMap.end())
    {
      WeightMap_t wmap;
      //Endex is the upper-edge for the energy of each bin
      double EnDex=0;
      for (int i=1; i<=hist->GetNbinsX();i++)
      {
        EnDex+=hist->GetBinWidth(i);
        wmap.insert(WeightMap_t::value_type(EnDex,hist->GetBinContent(i)));
      }
      fBeamSysMap.insert(std::map<mapkey, WeightMap_t, LessThan >::value_type(dex, wmap));
    }
    else
      std::cout << "Already loaded this histogram!" << std::endl;

    return;
  }

void nbw::skzpReweight::FillVector ( TH1F *  hist, int  ntype, int  eff, int  beam, int  det  )

private

Definition at line 351 of file skzpReweight.cxx.

References nbw::skzpReweight::mapkey::BeamDex, nbw::skzpReweight::mapkey::DetDex, nbw::skzpReweight::mapkey::EffDex, fBeamSysMap, and nbw::skzpReweight::mapkey::NuDex.

  {
    mapkey dex;
    dex.NuDex=ntype;
    dex.EffDex=eff;
    dex.BeamDex=beam;
    dex.DetDex=det;

    std::map<mapkey, WeightMap_t, LessThan>::iterator dexit = fBeamSysMap.find(dex);
    if (dexit == fBeamSysMap.end())
    {
      WeightMap_t wmap;
      //Endex is the upper-edge for the energy of each bin
      double EnDex=0;
      for (int i=1; i<hist->GetNbinsX()+1;i++)
      {
        EnDex+=hist->GetBinWidth(i);
        wmap.insert(WeightMap_t::value_type(EnDex,hist->GetBinContent(i)));
      }
      fBeamSysMap.insert(std::map<mapkey, WeightMap_t, LessThan >::value_type(dex, wmap));
    }
    else
      std::cout << "Already loaded this histogram!" << std::endl;

    return;
  }

void nbw::skzpReweight::FlukConfig ( )

private

Definition at line 146 of file skzpReweight.cxx.

References fFpath, fMeanPT, fN, fNWeighted, fPlist, fPTPZ, fWeightedMeanPT, fWeightedPTPZ, fWeightHist, GetFlukWeight(), hist, hist2, nbw::Conventions::kK0L, nbw::Conventions::kKMinus, nbw::Conventions::kKPlus, nbw::Conventions::kPiMinus, nbw::Conventions::kPiPlus, PartEnumToString(), and pt.

Referenced by GetFlukWeight(), SetFlukParams(), SetParams(), and skzpReweight().

  {
    std::cout <<"FlukConfig reading data from: "<<fFpath<<std::endl;

    //Builds maps and histograms for Fluk from file
    //In fluka05ptxf.root, file and histogram names refer to xf, but they actually mean pz, which is proportional to xF:The conversion takes place in this code and not the histogram.

    TFile* hFile=new TFile(fFpath.c_str());
    
    fPlist.push_back(Conventions::kPiPlus);
    fPlist.push_back(Conventions::kPiMinus);
    fPlist.push_back(Conventions::kKPlus);
    fPlist.push_back(Conventions::kKMinus);
    fPlist.push_back(Conventions::kK0L);
  
    for (std::vector<Conventions::ParticleType_t>::iterator itPlist=fPlist.begin();itPlist!=fPlist.end(); itPlist++)
    {
      TH2F* hist=dynamic_cast <TH2F*> (hFile->Get(Form("hF05ptxf%s",PartEnumToString(*itPlist).c_str()))->Clone());
      hist->SetDirectory(0);
      TH2F* hist2=dynamic_cast <TH2F*> (hist->Clone(Form("hWeightedPTXF%s",PartEnumToString(*itPlist).c_str())));
      hist2->SetDirectory(0);
      hist2->SetTitle(Form("%s weighted pt-pz",PartEnumToString(*itPlist).c_str()));
      fWeightHist[*itPlist]=dynamic_cast<TH2F*> (hist->Clone(Form("hWeight%s",PartEnumToString(*itPlist).c_str())));

      fWeightHist[*itPlist]->Divide(hist2);
      fWeightHist[*itPlist]->SetDirectory(0);

      std::pair<Conventions::ParticleType_t, TH2F* > p(*itPlist, hist);
      std::pair<Conventions::ParticleType_t, TH2F* > p2(*itPlist, hist2);
      fPTPZ.insert(p);
      fWeightedPTPZ.insert(p2);

      fPTPZ[*itPlist]->SetDirectory(0);
      fWeightedPTPZ[*itPlist]->SetDirectory(0);

      fMeanPT[*itPlist]=fPTPZ[*itPlist]->ProjectionY()->GetMean()*1000.;
      fWeightedMeanPT[*itPlist]=fWeightedPTPZ[*itPlist]->ProjectionY()->GetMean()*1000.;
            
      double N=fPTPZ[*itPlist]->ProjectionY()->GetSumOfWeights();
      double wN=fWeightedPTPZ[*itPlist]->ProjectionY()->GetSumOfWeights();
    
      fN[*itPlist]=N;
      fNWeighted[*itPlist]=wN;
    }
    hFile->Close();
    delete hFile;
    hFile=0;
    //Updates maps and histograms for Fluk with GetFlukWeight (the weight for KOL depends on fN and fNWeighted for KPlus and KMinus).
    for (std::vector<Conventions::ParticleType_t>::iterator itPlist=fPlist.begin(); itPlist!=fPlist.end(); itPlist++)
    {
      double pt,pz;
      double meanpt(0.);
      double sum(0.);
      for (int i=0;i<fPTPZ[*itPlist]->GetNbinsX()+1;i++)
      {
        for (int j=0;j<fPTPZ[*itPlist]->GetNbinsY()+1;j++)
        {
          pz=fPTPZ[*itPlist]->GetXaxis()->GetBinCenter(i);
          pt=fPTPZ[*itPlist]->GetYaxis()->GetBinCenter(j);
          fWeightedPTPZ[*itPlist]->SetBinContent(i,j,(fPTPZ[*itPlist]->GetBinContent(i,j)*GetFlukWeight(*itPlist,pt,pz)));
          fWeightHist[*itPlist]->SetBinContent(i,j,GetFlukWeight(*itPlist,pt,pz));
          meanpt+=fWeightedPTPZ[*itPlist]->GetBinContent(i,j)*pt;
          sum+=fWeightedPTPZ[*itPlist]->GetBinContent(i,j);
        }
      }
      meanpt/=sum;
      fWeightedMeanPT[*itPlist]=meanpt*1000.; //GeV to MeV
      fNWeighted[*itPlist]=sum;
      meanpt=0.;
      sum=0.;
    }
    return;
  }

Conventions::ParticleType_t nbw::skzpReweight::GeantToEnum ( int  ptype )

private

Definition at line 379 of file skzpReweight.cxx.

References nbw::Conventions::kK0L, nbw::Conventions::kKMinus, nbw::Conventions::kKPlus, nbw::Conventions::kPiMinus, nbw::Conventions::kPiPlus, and nbw::Conventions::kUnknown.

Referenced by GetFlukWeight().

  {
    switch (ptype)
      {
      case 8   : return Conventions::kPiPlus;  break;
      case 211 : return Conventions::kPiPlus;  break;    
      case 9   : return Conventions::kPiMinus; break;
      case -211: return Conventions::kPiMinus; break;
      case 11  : return Conventions::kKPlus;   break;
      case 321 : return Conventions::kKPlus;   break;
      case 12  : return Conventions::kKMinus;  break;
      case -321: return Conventions::kKMinus;  break;
      case 10  : return Conventions::kK0L;     break;
      case 130 : return Conventions::kK0L;     break;
      default  : return Conventions::kUnknown; break;
      }
    return Conventions::kUnknown;
  }

double nbw::skzpReweight::GetBeamWeight	(	int	ntype,
		double	Enu,
		int	det = 1,
		int	beam = 2
	)

Definition at line 221 of file skzpReweight.cxx.

References nbw::skzpReweight::mapkey::BeamDex, nbw::skzpReweight::mapkey::DetDex, nbw::skzpReweight::mapkey::EffDex, evd::details::end(), fBeamSysMap, fBPar, nbw::skzpReweight::mapkey::NuDex, w, and weight.

Referenced by GetFlukWeight(), and GetWeight().

  {
    double weight=1.;
    
    if (ntype == 14) ntype = 56;
    else if (ntype == -14) ntype = 55;
    else if (ntype == 12) ntype = 53;
    else if (ntype == -12) ntype = 52;

    struct mapkey dex;
    dex.NuDex=ntype;
    dex.BeamDex=beam;
    dex.DetDex=det;

    for(int sys = 1; sys <= 2; ++sys){//Loops through each beam-focusing correction applied
      double w = 0.;
      dex.EffDex=sys;
      
      std::map<mapkey, WeightMap_t, LessThan>::iterator dexit = fBeamSysMap.find(dex);
      if (dexit != fBeamSysMap.end()){
        for (WeightMap_t::iterator EnDex = (dexit->second).begin();  EnDex!=(dexit->second).end(); EnDex++){
          if(EnDex->first > Enu){//the first EnDex that goes over is mapped to the weight value
            w = EnDex->second;
            break;
          }
        }
      }
      weight *= std::abs(w)*fBPar[sys-1]+1.;
    }    
    return weight;
  }

double nbw::skzpReweight::GetFlukWeight ( const simb::MCFlux *  mcf )

inline

Definition at line 45 of file skzpReweight.h.

References FlukConfig(), simb::MCFlux::ftptype, simb::MCFlux::ftpx, simb::MCFlux::ftpy, simb::MCFlux::ftpz, GeantToEnum(), GetBeamWeight(), PartEnumToString(), and pt.

Referenced by FlukConfig(), and GetWeight().

                                               {
      double pt = sqrt(mcf->ftpx*mcf->ftpx + mcf->ftpy*mcf->ftpy);
      return GetFlukWeight(mcf->ftptype,pt,mcf->ftpz); };

double nbw::skzpReweight::GetFlukWeight	(	int	ptype,
		double	pT,
		double	pz
	)

Definition at line 73 of file skzpReweight.cxx.

References B, fFPar, fN, fNWeighted, GeantToEnum(), nbw::Conventions::kK0L, nbw::Conventions::kKMinus, nbw::Conventions::kKPlus, nbw::Conventions::kPiMinus, nbw::Conventions::kPiPlus, and weight.

  {
    double weight=1.;
    double xF=pz/120.;
    double A,B,C;
    double Ap,Bp,Cp;
    Conventions::ParticleType_t eptype=GeantToEnum(ptype);

    // This is SKZP parameterization
    if (xF>1.||xF<0.) return weight;
    if (pT>1.||pT<0.) return weight;

    if (eptype==Conventions::kPiPlus || eptype==Conventions::kPiMinus)
    {
      //Calculate weight for pions
      if (pT<0.03) // for low pT
        pT=0.03;
      // calculate the A, B and C in SKZP parameterization
      // A, B and C are best fit to Fluka 05
      
      A=-0.00761*pow((1.-xF),4.045)*(1.+9620.*xF)*pow(xF,-2.975);
      B=0.05465*pow((1.-xF),2.675)*(1.+69590.*xF)*pow(xF,-3.144);
      
      if (xF<0.22)
        C=-7.058/pow(xF,-0.1419)+9.188;
      else
        C = (3.008/exp((xF-0.1984)*3.577)) + 2.616*xF+0.1225;
      
      // scale/skew A, B and C
      Ap=(fFPar[0]+fFPar[1]*xF)*A;
      Bp=(fFPar[2]+fFPar[3]*xF)*B;
      Cp=(fFPar[4]+fFPar[5]*xF)*C;
      
      weight=(Ap+Bp*pT)/(A+B*pT)*exp(-(Cp-C)*pow(pT,3./2.));
      if(eptype == Conventions::kPiMinus && pz > 4)
        weight *= ( fFPar[12] + fFPar[13] * xF );
    }
    else if (eptype==Conventions::kKPlus || eptype==Conventions::kKMinus)
    {
      //Calculate weight for kaons
      if (pT<0.05) // for low pT
        pT=0.05;
      // calculate the A, B and C in SKZP parameterization
      // A, B and C are best fit to Fluka 05
          
      A=-0.005187*pow((1.-xF),4.119)*(1.+2170.*xF)*pow(xF,-2.767);
      B=0.4918*pow((1.-xF),2.672)*(1.+1373.*xF)*pow(xF,-2.927);
      
      if (xF<0.22)
        C=-16.10/pow(xF,-0.04582)+17.92;
      else
        C = (6.905/exp((xF+0.163)*6.718)) - 0.4257*xF + 2.486;
          
      // scale/skew A, B and C
      Ap=(fFPar[6]+fFPar[7]*xF)*A;
      Bp=(fFPar[8]+fFPar[9]*xF)*B;
      Cp=(fFPar[10]+fFPar[11]*xF)*C;
          
      weight=(Ap+Bp*pT)/(A+B*pT)*exp(-(Cp-C)*pow(pT,3./2.));
      if(eptype==Conventions::kKMinus)
        weight*=(fFPar[14]+fFPar[15]*xF);
    }
    else if (eptype==Conventions::kK0L)
    {
      //N(K0L) approximately given with (N(K+)+3*N(K-))/4
      weight= ((fNWeighted[Conventions::kKPlus]+3.*fNWeighted[Conventions::kKMinus])/(fN[Conventions::kKPlus]+3.*fN[Conventions::kKMinus]));
    }
    
    if (weight>10.) weight=10.;
    return weight;
  }

std::string nbw::skzpReweight::GetHname ( int  inu, int  eff, int  beam, int  det  )

private

Definition at line 306 of file skzpReweight.cxx.

References BeamSysToString(), BeamTypeToString(), DetTypeToString(), and fBflag.

Referenced by BeamConfig().

  {
    std::string hname;
    if (fBflag == 1 || fBflag == 2)
    {
      std::string nus[]={"NuMu","NuMuBar","NuE","NuEBar"};
      hname=nus[inu]+"_";
    }
    hname=hname + BeamSysToString(Conventions::BeamSys_t(eff));
    if (fBflag == 1 || fBflag == 2) hname = hname + "_";
    hname=hname + BeamTypeToString(Conventions::BeamType_t(beam));
    if (fBflag == 1 || fBflag == 2) hname = hname + "_";
    hname=hname + DetTypeToString(Conventions::DetType_t(det));
    return hname;
  }

double nbw::skzpReweight::GetWeight ( const simb::MCFlux *  mcf, double  Enu, int  det, int  beam  )

inline

Definition at line 41 of file skzpReweight.h.

References simb::MCFlux::fntype, simb::MCFlux::ftptype, simb::MCFlux::ftpx, simb::MCFlux::ftpy, simb::MCFlux::ftpz, GetBeamWeight(), GetFlukWeight(), and pt.

                                                                          {      
      double pt = sqrt(mcf->ftpx*mcf->ftpx + mcf->ftpy*mcf->ftpy);
      return GetFlukWeight(mcf->ftptype,pt,mcf->ftpz)*GetBeamWeight(mcf->fntype,Enu,beam,det); };

std::string nbw::skzpReweight::PartEnumToString ( Conventions::ParticleType_t  ptype )

private

Definition at line 398 of file skzpReweight.cxx.

References nbw::Conventions::kK0L, nbw::Conventions::kKMinus, nbw::Conventions::kKPlus, nbw::Conventions::kPiMinus, nbw::Conventions::kPiPlus, and nbw::Conventions::kUnknown.

Referenced by FlukConfig(), and GetFlukWeight().

  {
    switch (ptype)
      {
      case Conventions::kPiPlus : return "PiPlus";  break;
      case Conventions::kPiMinus: return "PiMinus"; break;
      case Conventions::kKPlus  : return "KPlus";   break;
      case Conventions::kKMinus : return "KMinus";  break;
      case Conventions::kK0L    : return "K0L";     break;
      case Conventions::kUnknown: return "Unknown"; break;
      default                   : return "Unknown"; break;
      }
    return "Unknown";
  }

void nbw::skzpReweight::SetBeamParams ( std::vector< double >  bpar )

inline

Definition at line 37 of file skzpReweight.h.

References BeamConfig(), and fBPar.

                                               {
      fBPar=bpar; BeamConfig();
      return; };

void nbw::skzpReweight::SetFlukParams ( std::vector< double >  fpar )

inline

Definition at line 34 of file skzpReweight.h.

References fFPar, and FlukConfig().

                                               {
      fFPar=fpar; FlukConfig();
      return; };

void nbw::skzpReweight::SetParams ( std::vector< double >  fpar, std::vector< double >  bpar  )

inline

Definition at line 30 of file skzpReweight.h.

References BeamConfig(), fBPar, fFPar, and FlukConfig().

                                                                   {
      fFPar=fpar; FlukConfig();
      fBPar=bpar; BeamConfig();
      return; };

Member Data Documentation

TFile* nbw::skzpReweight::fBeamSysFile

private

Definition at line 117 of file skzpReweight.h.

Referenced by BeamConfig().

std::map<mapkey, WeightMap_t, LessThan > nbw::skzpReweight::fBeamSysMap

private

Definition at line 121 of file skzpReweight.h.

Referenced by FillVector(), and GetBeamWeight().

int nbw::skzpReweight::fBflag

private

Definition at line 126 of file skzpReweight.h.

Referenced by BeamConfig(), BeamSysToString(), BeamTypeToString(), DetTypeToString(), GetHname(), and skzpReweight().

std::vector<double> nbw::skzpReweight::fBPar

private

Definition at line 115 of file skzpReweight.h.

Referenced by GetBeamWeight(), SetBeamParams(), SetParams(), and skzpReweight().

std::string nbw::skzpReweight::fBpath

private

Definition at line 116 of file skzpReweight.h.

Referenced by BeamConfig(), and skzpReweight().

std::vector<double> nbw::skzpReweight::fFPar

private

Definition at line 59 of file skzpReweight.h.

Referenced by GetFlukWeight(), SetFlukParams(), SetParams(), skzpReweight(), and ~skzpReweight().

std::string nbw::skzpReweight::fFpath

private

Definition at line 60 of file skzpReweight.h.

Referenced by FlukConfig(), and skzpReweight().

std::map<Conventions::ParticleType_t, double > nbw::skzpReweight::fMeanPT

private

Definition at line 65 of file skzpReweight.h.

Referenced by FlukConfig(), and ~skzpReweight().

std::map<Conventions::ParticleType_t, double > nbw::skzpReweight::fN

private

Definition at line 66 of file skzpReweight.h.

Referenced by FlukConfig(), GetFlukWeight(), and ~skzpReweight().

std::map<Conventions::ParticleType_t, int > nbw::skzpReweight::fNBinsX

private

Definition at line 69 of file skzpReweight.h.

Referenced by ~skzpReweight().

std::map<Conventions::ParticleType_t, int > nbw::skzpReweight::fNBinsY

private

Definition at line 69 of file skzpReweight.h.

Referenced by ~skzpReweight().

std::map<Conventions::ParticleType_t, double > nbw::skzpReweight::fNWeighted

private

Definition at line 67 of file skzpReweight.h.

Referenced by FlukConfig(), GetFlukWeight(), and ~skzpReweight().

std::vector<Conventions::ParticleType_t> nbw::skzpReweight::fPlist

private

Definition at line 61 of file skzpReweight.h.

Referenced by FlukConfig(), and ~skzpReweight().

std::map<Conventions::ParticleType_t, TH2F* > nbw::skzpReweight::fPTPZ

private

Definition at line 62 of file skzpReweight.h.

Referenced by FlukConfig(), and ~skzpReweight().

std::map<Conventions::ParticleType_t, double > nbw::skzpReweight::fWeightedMeanPT

private

Definition at line 68 of file skzpReweight.h.

Referenced by FlukConfig(), and ~skzpReweight().

std::map<Conventions::ParticleType_t, TH2F* > nbw::skzpReweight::fWeightedPTPZ

private

Definition at line 63 of file skzpReweight.h.

Referenced by FlukConfig(), and ~skzpReweight().

std::map<Conventions::ParticleType_t, TH2F* > nbw::skzpReweight::fWeightHist

private

Definition at line 64 of file skzpReweight.h.

Referenced by FlukConfig(), and ~skzpReweight().

---

The documentation for this class was generated from the following files:

• nutools/NuBeamWeights/skzpReweight.h
• nutools/NuBeamWeights/skzpReweight.cxx