PhotonStrucFunc.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!

\class    genie::PhotonStrucFunc

\brief    Structure function using photon PDFs of nucleons

\author   Alfonso Garcia <aagarciasoto \at km3net.de>
          IFIC & Harvard University

\created  Dec 8, 2021

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

#ifndef _PHOTON_STRUC_FUNC_H_
#define _PHOTON_STRUC_FUNC_H_

#include "Framework/Numerical/Spline.h"

#include <map>

using std::map;

namespace genie {

  struct SF_x {
    double Fp1,Fm1;
    double Fp2,Fm2;
    double Fp3,Fm3;
  };

  class PhotonStrucFunc
  {
    public:

      // ................................................................
      // GLRES form factor type
      //

      class PhotonStrucFuncTable 
      {
        public:
          PhotonStrucFuncTable() { }
          ~PhotonStrucFuncTable() { /* note: should delete the grids! */ }
          map< int, genie::Spline * > Table;
      };

      // ................................................................

      static PhotonStrucFunc * Instance(void);

      double EvalSF  ( int hitnuc, int hitlep, double x ) { return fSFTables[hitnuc].Table[hitlep]->Evaluate(x); }

    private:

      // Ctors & dtor
      PhotonStrucFunc();
      PhotonStrucFunc(const PhotonStrucFunc &);
     ~PhotonStrucFunc();

      // Self
      static PhotonStrucFunc * fgInstance;

      // These map holds all SF tables (interaction channel is the key)
      map<int, PhotonStrucFuncTable> fSFTables;


      // singleton cleaner
      struct Cleaner {
        void DummyMethodAndSilentCompiler(){}
          ~Cleaner(){
          if (PhotonStrucFunc::fgInstance !=0){
            delete PhotonStrucFunc::fgInstance;
            PhotonStrucFunc::fgInstance = 0;
          }
        }
      };
      friend struct Cleaner;
  };

} // genie namespace

#endif // _PHOTON_STRUC_FUNC_H_