BertuzzoDNuCOHPXSec.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
  \class    genie::BertuzzoDNuCOHPXSec

  \brief    Differential cross section for v+As coherent elastic scattering.Coherent DNu. \n
            Is a concrete implementation of the XSecAlgorithmI interface. \n

  \ref      E.Bertuzzo, S.Jana, P.A.N.Machado, R. Zukanovich Funcal
            PhysRevLett.121.241801 (2018)

  \author   Author: Iker de Icaza <i.de-icaza-astiz \at sussex.ac.uk>
            University of Sussex

            Costas Andreopoulos <constantinos.andreopoulos \at cern.ch>
            University of Liverpool & STFC Rutherford Appleton Laboratory

  \created  June 12, 2020

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

*/
//____________________________________________________________________________

#ifndef _BERTUZZO_DNu_COH_CROSS_SECTION_H_
#define _BERTUZZO_DNu_COH_CROSS_SECTION_H_

#include "Framework/EventGen/XSecAlgorithmI.h"
#include "Physics/DarkNeutrino/XSection/EngelFormFactor.h"

namespace genie {

class XSecIntegratorI;

class BertuzzoDNuCOHPXSec : public XSecAlgorithmI {

public:
  BertuzzoDNuCOHPXSec();
  BertuzzoDNuCOHPXSec(string config);
  virtual ~BertuzzoDNuCOHPXSec();

  // XSecAlgorithmI interface implementation
  double XSec            (const Interaction * i, KinePhaseSpace_t k) const override;
  double Integral        (const Interaction * i) const override;
  bool   ValidProcess    (const Interaction * i) const override;
  bool   ValidKinematics (const Interaction * i) const override;

  // Override the Algorithm::Configure methods to load configuration
  // data to private data members
  void Configure (const Registry & config) override;
  void Configure (string param_set) override;

private:

  void LoadConfig(void);

  const XSecIntegratorI * fXSecIntegrator;  ///< cross section integrator
  const EngelFormFactor * fFF; ///< Engel Form Factor algorithm

  double fEps2;

  std::array<double, 4> fMixing2s;
  // mixing angles between the neutrinos (including a 4rth sterile one)
  // to the dark neutrino

  double fAlpha_D;

  double fDNuMass, fDNuMass2;
  double fDMediatorMass, fDMediatorMass2;

};

} // genie namespace


#endif // _BERTUZZO_DNu_COH_CROSS_SECTION_H_