sh_zN_grid_source_const.h

Go to the documentation of this file.

// -*- C++ -*-
/*! \file
 *  \brief Shell source construction
 * this source implements a ZN (were N = Nc) for baryonic states 
 * It also works on mesoning 2 quark states.
 * The idea is to increase statistics with a grid source construction
 * This trick has been used by the Kentucky group. 
 * I think all I am doing is implementing what they have already done.
 * I know that there are many grid sources in Chroma but
 * I what one that works as I want...
 */
 
#ifndef __sh_zN_grid_source_const_h__
#define __sh_zN_grid_source_const_h__
 
#include "meas/sources/source_construction.h"
#include "io/xml_group_reader.h"
 
namespace Chroma
{
 
  //! Name and registration
  /*! @ingroup sources */
  namespace ShellZnGridQuarkSourceConstEnv
  {
    bool registerAll();
 
    //! Return the name
    std::string getName();
 
    //! Point source parameters
    /*! @ingroup sources */
    struct Params
    {
      Params();
      Params(XMLReader& in, const std::string& path);
      void writeXML(XMLWriter& out, const std::string& path) const;
    
      bool             quark_smear_lastP;    /*!< Flag controlling order of smearing */
 
      GroupXML_t       quark_smearing;       /*!< xml std::string holding smearing params */
      GroupXML_t       quark_displacement;   /*!< displacement xml */
      GroupXML_t       link_smearing;        /*!< link smearing xml */
 
      int              j_decay;              /*!< Decay direction */
      
      multi1d<int>     t_srce;               /*!< source location */
      multi1d<int>     grd   ;               /*!< grid spacing negative means no grid in this direction. */
      multi1d<multi1d<int> > moms ; /*!<  If not empty add a linear combination of momenta to the source */
      Seed ran_seed; /*!< Set the seed to this value */
    };
 
 
    //! Shell source construction
    /*! @ingroup sources
     *
     * Create a shell quark source
     */
    template<typename T>
    class SourceConst : public QuarkSourceConstruction<T>
    {
    public:
      //! Full constructor
      SourceConst(const Params& p) : params(p) {}
 
      //! Construct the source
      T operator()(const multi1d<LatticeColorMatrix>& u) const;
 
    private:
      //! Hide partial constructor
      SourceConst() {}
 
    private:
      Params  params;   /*!< source params */
    };
 
  }  // end namespace ShellQuarkSourceConstEnv
 
 
  //! Reader
  /*! @ingroup sources */
  void read(XMLReader& xml, const std::string& path, ShellZnGridQuarkSourceConstEnv::Params& param);
 
  //! Writer
  /*! @ingroup sources */
  void write(XMLWriter& xml, const std::string& path, const ShellZnGridQuarkSourceConstEnv::Params& param);
 
 
}  // end namespace Chroma
 
 
#endif