McPairEnergyAverage.h

#ifndef SIMP_NOPAIR
#ifndef MCMD_MC_PAIR_ENERGY_AVERAGE_H
#define MCMD_MC_PAIR_ENERGY_AVERAGE_H

/*
* Simpatico - Simulation Package for Polymeric and Molecular Liquids
*
* Copyright 2010 - 2017, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/

#include <mcMd/analyzers/SystemAnalyzer.h>   // base class template
#include <mcMd/mcSimulation/McSystem.h>          // base template parameter
#include <util/accumulators/Average.h>           // member
#include <mcMd/neighbor/CellList.h>              // member
#include <mcMd/analyzers/util/PairSelector.h>  // member

#include <cstdio> 

namespace McMd
{

   using namespace Util;

   class McPairEnergyAverage : public SystemAnalyzer<McSystem>
   {
   
   public:

      McPairEnergyAverage(McSystem& system);

      virtual void readParameters(std::istream& in);

      virtual void loadParameters(Serializable::IArchive& ar);

      virtual void save(Serializable::OArchive& ar);
  
      template <class Archive>
      void serialize(Archive& ar, const unsigned int version);

      virtual void setup();
   
      /* 
      * Evaluate energy per particle, and add to ensemble. 
      */
      virtual void sample(long iStep);
   
      virtual void output();

   private:

      std::ofstream outputFile_;

      Average  accumulator_;

      int nSamplePerBlock_;

      mutable CellList::NeighborArray neighbors_;

      PairSelector  selector_;

      // Has readParam been called?
      bool isInitialized_;

   };

   /*
   * Serialize to/from an archive. 
   */
   template <class Archive>
   void McPairEnergyAverage::serialize(Archive& ar, const unsigned int version)
   {
      Analyzer::serialize(ar, version);
      ar & nSamplePerBlock_;
      ar & selector_;
      ar & accumulator_;
   }

}
#endif 
#endif