MECScaleVsW.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!

\class    genie::MECScaleVsW

\brief    This class is responsible to compute the MEC scaling factor given 
          Q0, Q3. The scaling is done as a function of the hadronic invariant
          mass.

\author   Code contributed by J.Tena Vidal and M.Roda

\created  June, 2020

\cpright  Copyright (c) 2003-2020, The GENIE Collaboration
          For the full text of the license visit http://copyright.genie-mc.org
*/
//____________________________________________________________________________

#ifndef _MEC_SCALE_VS_W_H_
#define _MEC_SCALE_VS_W_H_

#include "Physics/Common/XSecScaleI.h"
#include <TSpline.h>

namespace genie {
  
  class MECScaleVsW: public XSecScaleI {

    using weight_type_map = std::map<double,double> ;
    using weight_type_pair = std::pair<double,double> ;

  public:
    MECScaleVsW();
    MECScaleVsW(string config);
    virtual ~MECScaleVsW();
    
    // This function returns the scaling value at a given Q0 Q3:
    virtual double GetScaling( const Interaction & ) const override ; 
    
  protected:
    
    // Load algorithm configuration
    virtual void LoadConfig (void) override ;

    // This function returns the scaling value at a given Q0 Q3:
    virtual double GetScaling( const double Q0, const double Q3 ) const ; 

    // This function adds the limits of the phase space if they are not set by the user
    weight_type_map GetMapWithLimits( const double Q0, const double Q3 ) const ;

    // Thist function calculates the scale factor value at W as a linear interpolation
    // between two W values (Wmin,Wmax) with weights (scale_min,scale_max).
    virtual double ScaleFunction( const double W, const weight_type_pair min, const weight_type_pair max ) const ;

 private: 
    double fDefaultWeight ; 
    weight_type_map fWeightsMap ;
    // Adding Spline to handle the limits of W1:
    TSpline3 fW1_Q0Q3_limits ; 

    double fLowLimitWeight ; 
    double fUpperLimitWeight ; 

  };
  
}       // genie namespace
#endif  // _MEC_SCALE_VS_W_H_