libthepeg-dev 1.8.0-1.1 / usr / include / ThePEG / Analysis / XSecCheck.h

// -*- C++ -*-
#ifndef THEPEG_XSecCheck_H
#define THEPEG_XSecCheck_H
//
// This is the declaration of the XSecCheck class.
//

#include "ThePEG/Handlers/AnalysisHandler.h"

namespace ThePEG {

/**
 * The XSecCheck class is a simple analysis class used for testing
 * purposes. If the total cross section does not match the one
 * specified by TargetXSec, an exception will
 * be thrown.
 *
 * @see \ref XSecCheckInterfaces "The interfaces"
 * defined for XSecCheck.
 */
class XSecCheck: public AnalysisHandler {

public:

  /** @name Standard constructors and destructors. */
  //@{
  /**
   * The default constructor.
   */
  XSecCheck() : target(ZERO), tol(0.01), sumw(0.0) {}

  /**
   * The destructor.
   */
  virtual ~XSecCheck();
  //@}

public:

  /**
   * The exception class used if the target cross section was not met.
   */
  struct UnexpectedXSec: public Exception {};

public:

  /** @name Virtual functions required by the AnalysisHandler class. */
  //@{
  /**
   * Analyze a given Event. Note that a fully generated event
   * may be presented several times, if it has been manipulated in
   * between. The default version of this function will call transform
   * to make a lorentz transformation of the whole event, then extract
   * all final state particles and call analyze(tPVector) of this
   * analysis object and those of all associated analysis objects. The
   * default version will not, however, do anything on events which
   * have not been fully generated, or have been manipulated in any
   * way.
   * @param event pointer to the Event to be analyzed.
   * @param ieve the event number.
   * @param loop the number of times this event has been presented.
   * If negative the event is now fully generated.
   * @param state a number different from zero if the event has been
   * manipulated in some way since it was last presented.
   */
  virtual void analyze(tEventPtr event, long ieve, int loop, int state);
  //@}

public:

  /** @name Functions used by the persistent I/O system. */
  //@{
  /**
   * Function used to write out object persistently.
   * @param os the persistent output stream written to.
   */
  void persistentOutput(PersistentOStream & os) const;

  /**
   * Function used to read in object persistently.
   * @param is the persistent input stream read from.
   * @param version the version number of the object when written.
   */
  void persistentInput(PersistentIStream & is, int version);
  //@}

  /**
   * The standard Init function used to initialize the interfaces.
   * Called exactly once for each class by the class description system
   * before the main function starts or
   * when this class is dynamically loaded.
   */
  static void Init();

protected:

  /** @name Clone Methods. */
  //@{
  /**
   * Make a simple clone of this object.
   * @return a pointer to the new object.
   */
  virtual IBPtr clone() const;

  /** Make a clone of this object, possibly modifying the cloned object
   * to make it sane.
   * @return a pointer to the new object.
   */
  virtual IBPtr fullclone() const;
  //@}



protected:

  /** @name Standard Interfaced functions. */
  //@{
  /**
   * Initialize this object. Called in the run phase just before
   * a run begins.
   */
  virtual void doinitrun() {
    AnalysisHandler::doinitrun();
    sumw = 0.0;
  }

  /**
   * Finalize this object. Called in the run phase just after a
   * run has ended. Used eg. to write out statistics.
   */
  virtual void dofinish();
  //@}

private:

  /**
   * The expected total cross section.
   */
  CrossSection target;

  /**
   * The relative tolerance accepted when comparing the total cross
   * section with the target.
   */
  double tol;

  /**
   * The sum of the weights of the events analyzed.
   */
  double sumw;

private:

  /**
   * The static object used to initialize the description of this class.
   * Indicates that this is a concrete class with persistent data.
   */
  static ClassDescription<XSecCheck> initXSecCheck;

  /**
   * The assignment operator is private and must never be called.
   * In fact, it should not even be implemented.
   */
  XSecCheck & operator=(const XSecCheck &);

};

}

#include "ThePEG/Utilities/ClassTraits.h"

namespace ThePEG {

/** @cond TRAITSPECIALIZATIONS */

/** This template specialization informs ThePEG about the
 *  base classes of XSecCheck. */
template <>
struct BaseClassTrait<XSecCheck,1> {
  /** Typedef of the first base class of XSecCheck. */
  typedef AnalysisHandler NthBase;
};

/** This template specialization informs ThePEG about the name of
 *  the XSecCheck class and the shared object where it is defined. */
template <>
struct ClassTraits<XSecCheck>
  : public ClassTraitsBase<XSecCheck> {
  /** Return a platform-independent class name */
  static string className() { return "ThePEG::XSecCheck"; }
  /**
   * The name of a file containing the dynamic library where the class
   * XSecCheck is implemented. It may also include several, space-separated,
   * libraries if the class XSecCheck depends on other classes (base classes
   * excepted). In this case the listed libraries will be dynamically
   * linked in the order they are specified.
   */
  static string library() { return "XSecCheck.so"; }
};

/** @endcond */

}

#endif /* THEPEG_XSecCheck_H */