GReWeightINukeParams.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!

\class   genie::rew::GReWeightINukeParams

\brief   Helper class for cross section model reweighting

\author  Costas Andreopoulos <costas.andreopoulos \at stfc.ac.uk>
         University of Liverpool & STFC Rutherford Appleton Lab

         Jim Dobson <J.Dobson07 \at imperial.ac.uk>
         Imperial College London

\created Sep 10, 2009

\cpright  Copyright (c) 2003-2017, GENIE Neutrino MC Generator Collaboration
          For the full text of the license visit http://copyright.genie-mc.org
          or see $GENIE/LICENSE
*/
//____________________________________________________________________________

#ifndef _G_REWEIGHT_INTRANUKE_PARAMS_H_
#define _G_REWEIGHT_INTRANUKE_PARAMS_H_

#include <map>

#include "PDG/PDGUtils.h"
#include "ReWeight/GSyst.h"

using std::map;

class TLorentzVector;

namespace genie {
namespace rew   {

 class GReWeightINukeParams {

 public:

   typedef enum EHadronType {
      kRwINukeUndefined = 0,
      kRwINukePion,
      kRwINukeNucl,
   } HadronType_t;

   static HadronType_t HadronTypeFromPdg(int pdgc) {
     if(pdg::IsPion   (pdgc)) return kRwINukePion;
     if(pdg::IsNucleon(pdgc)) return kRwINukeNucl;
     return kRwINukeUndefined;
   }

   GReWeightINukeParams();
  ~GReWeightINukeParams();

   class Fates;
   class MFP;

   Fates * FateParams         (int pdgc) const;         ///<
   MFP *   MeanFreePathParams (int pdgc) const;         ///<
   void    Reset              (void);                   ///<
   void    Reconfigure        (void);                   ///<
   double  ChisqPenalty       (void) const;             ///<
   void    SetTwkDial         (GSyst_t s, double val);  ///<

   //.........................................................................
   //
   // nested class: Fates
   //
   //.........................................................................

   class Fates {
   public :
     Fates(HadronType_t hadtype = kRwINukeUndefined);
    ~Fates();

     double ScaleFactor   (GSyst_t s, const TLorentzVector & p4) const; ///< see next
     double ScaleFactor   (GSyst_t s, double KE=-1.) const;             ///< fate fraction scale factor = 1 + twk_dial * fractional_err
     bool   IsIncluded    (GSyst_t s) const;                            ///< is included?
     bool   IsCushionTerm (GSyst_t s) const;                            ///< is it a cushion term?
     bool   IsTweaked     (GSyst_t s) const;                            ///< is included & tweaked to non-def value?
     bool   IsTweaked     (void) const;                                 ///< is any param tweaked
     void   Reset         (void);                                       ///<
     void   Reconfigure   (void);                                       ///<
     double ChisqPenalty  (void) const;                                 ///<
     void   SetTwkDial    (GSyst_t s, double val);                      ///< 

   private:    

     bool   IsHandled       (GSyst_t s) const;
     void   AddCushionTerms (void);
     double ActualTwkDial   (GSyst_t s, double KE=-1.) const;  ///< actual tweaking dial for input systematic at input kinetic energy

     HadronType_t         fHadType;           ///<
     map<GSyst_t, double> fSystValuesUser;    ///< List of systematics included & values set by the user
     mutable 
     map<GSyst_t, double> fSystValuesActual;  ///< List of systematics included & values actually used (user values limited to physical range)
     map<GSyst_t, bool>   fIsCushion;         ///< cushion term flag

   }; // Fates nested class


   //.........................................................................
   //
   // nested class: MFP
   //
   //.........................................................................

   class MFP {
   public :
     MFP(HadronType_t hadtype = kRwINukeUndefined);
    ~MFP();

     double ScaleFactor   (void) const;  ///< mean free path scale factor = 1 + twk_dial * fractional_err
     double TwkDial       (void) const;  ///< current value of mfp tweak dial
     bool   IsIncluded    (void) const;  ///<
     bool   IsTweaked     (void) const;  ///<
     double ChisqPenalty  (void) const;  ///<
     void   Reset         (void);        ///<
     void   SetTwkDial    (double val);  ///<

   private:    
     HadronType_t fHadType;     ///<
     GSyst_t      fSyst;        ///<
     double       fTwkDial;     ///<
     bool         fIsIncluded;  ///<

   }; // MFP nested class


 private:

    void Init(void);

    Fates * fParmPionFates;
    Fates * fParmNuclFates; 
    MFP *   fParmPionMFP;   
    MFP *   fParmNuclMFP;

 }; //GReWeightINukeParams

}      // rew   namespace
}      // genie namespace

#endif // _G_REWEIGHT_INTRANUKE_PARAMS_H_