QPMDMDISPXSec.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
 
\class    genie::QPMDMDISPXSec
 
\brief    Computes DMDIS differential cross sections.
          Is a concrete implementation of the XSecAlgorithmI interface.
 
\ref      E.A.Paschos and J.Y.Yu, Phys.Rev.D 65.03300
 
\author   Joshua Berger <jberger \at physics.wisc.edu
          University of Wisconsin-Madison
 
          Costas Andreopoulos <c.andreopoulos \at cern.ch>
          University of Liverpool
 
\created  September 4, 2017
 
\cpright  Copyright (c) 2003-2025, The GENIE Collaboration
          For the full text of the license visit http://copyright.genie-mc.org
          
*/
//____________________________________________________________________________
 
#ifndef _DMDIS_PARTON_MODEL_PARTIAL_XSEC_H_
#define _DMDIS_PARTON_MODEL_PARTIAL_XSEC_H_
 
#include "Framework/EventGen/XSecAlgorithmI.h"
#include "Physics/DeepInelastic/XSection/DISStructureFunc.h"
 
namespace genie {
 
class DISStructureFuncModelI;
class HadronizationModelI;
class XSecIntegratorI;
 
class QPMDMDISPXSec : public XSecAlgorithmI {
 
public:
  QPMDMDISPXSec();
  QPMDMDISPXSec(string config);
  virtual ~QPMDMDISPXSec();
 
  // XSecAlgorithmI interface implementation
  double XSec            (const Interaction * i, KinePhaseSpace_t k) const;
  double Integral        (const Interaction * i) const;
  bool   ValidProcess    (const Interaction * i) const;
 
  // overload the Algorithm::Configure() methods to load private data
  // members from configuration options
  void Configure(const Registry & config);
  void Configure(string config);
 
private:
  void   LoadConfig                  (void);
  double DMDISRESJoinSuppressionFactor (const Interaction * in) const;
 
  mutable DISStructureFunc fDISSF;
  bool                     fInInitPhase;
 
  const DISStructureFuncModelI * fDISSFModel;         ///< SF model
  const HadronizationModelI *    fHadronizationModel; ///< hadronic multip. model
  const XSecIntegratorI *        fXSecIntegrator;     ///< diff. xsec integrator
 
  const XSecAlgorithmI * fCharmProdModel;
 
  bool   fUsingDisResJoin;  ///< use a DMDIS/RES joining scheme?
  bool   fUseCache;         ///< cache reduction factors used in joining scheme
  double fWcut;             ///< apply DMDIS/RES joining scheme < Wcut
  double fCCScale;          ///< cross section scaling factor for CC processes
  double fNCScale;          ///< cross section scaling factor for NC processes
  double fEMScale;          ///< cross section scaling factor for EM processes
  double fSin48w;           ///< sin^4(Weingberg angle)
  int    fVelMode;          ///< velcoity dependence for xsec
  double fMedMass;          ///< Mediator mass
  double fgzp;              ///< Coupling to the mediator Zprime
  double fQchiL;            ///< Left-handed DM charge
  double fQchiR;            ///< Right-handed DM charge
  double fQchiS;            ///< Scalar DM charge
};
 
}       // genie namespace
#endif  // _DMDIS_PARTON_MODEL_PARTIAL_XSEC_H_