DMELUtils.h

Go to the documentation of this file.

#ifndef _DMEL_UTILS_H_
#define _DMEL_UTILS_H_
 
#include "Framework/Interaction/Interaction.h"
#include "Physics/NuclearState/NuclearModelI.h"
#include "Framework/EventGen/XSecAlgorithmI.h"
 
#include "TLorentzVector.h"
#include "TVector3.h"
 
#include <string>
 
namespace genie {
 
  // Enumerated type used to specify the method for determining the off-shell energy
  // of the hit nucleon for quasielastic events
  typedef enum EDMELEvGenBindingMode {
 
    // Use removal energy from the nuclear model
    kUseNuclearModel,
 
    // Calculate binding energy assuming that the remnant nucleus is left in its
    // ground state
    kUseGroundStateRemnant,
 
    // Leave the struck nucleon on shell, effectively ignoring its binding energy
    kOnShell
  } DMELEvGen_BindingMode_t;
 
  namespace utils {
 
    double EnergyDeltaFunctionSolutionDMEL(const Interaction& inter);
 
    DMELEvGen_BindingMode_t StringToDMELBindingMode( const std::string& mode_str );
 
    double ComputeFullDMELPXSec(Interaction* interaction,
      const NuclearModelI* nucl_model, const XSecAlgorithmI* xsec_model,
      double cos_theta_0, double phi_0, double& Eb,
      DMELEvGen_BindingMode_t hitNucleonBindingMode, double min_angle_EM = 0.,
      bool bind_nucleon = true);
 
    double CosTheta0Max(const genie::Interaction& interaction);
 
    void BindHitNucleon(Interaction& interaction, const NuclearModelI& nucl_model,
      double& Eb, DMELEvGen_BindingMode_t hitNucleonBindingMode);
  }
}
 
#endif