HardProcessAnalysis.h

// -*- C++ -*-
#ifndef Herwig_HardProcessAnalysis_H
#define Herwig_HardProcessAnalysis_H
//
// This is the declaration of the HardProcessAnalysis class.
//
 
#include "ThePEG/Handlers/AnalysisHandler.h"
#include "Herwig/Utilities/Histogram.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class HardProcessAnalysis: public AnalysisHandler {
 
public:
 
  HardProcessAnalysis();
 
public:
 
  virtual void analyze(tEventPtr event, long ieve, int loop, int state);
 
 
protected:
 
  virtual void doinitrun();
 
  virtual void dofinish();
 
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
 
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
 
 
private:
 
  HardProcessAnalysis & operator=(const HardProcessAnalysis &) = delete;
 
  struct Histograms {
 
    Histograms() {}
 
    explicit Histograms(Energy ECM, unsigned int theNBins);
 
    void fill(const Lorentz5Momentum& p, double weight);
 
    void finalize(ostream& dat,
                  ostream& plot,
                  const string& subpro,
                  size_t legid,
                  double norm,
                  bool theUnitWeights);
 
    HistogramPtr transverse;
 
    HistogramPtr rapidity;
 
    HistogramPtr phi;
 
  };
 
  struct AllHistograms {
 
    vector<Histograms> outgoing;
 
    HistogramPtr x1;
 
    HistogramPtr x2;
 
    HistogramPtr sshat;
 
    HistogramPtr rapidity;
 
    double sumWeights;
 
  };
 
  map<vector<string>,AllHistograms> histogramData;
 
  double sumWeights;
 
  void fill(PPair, ParticleVector, double);
 
  unsigned int theNBins;
 
  bool theUnitWeights;
 
  bool theSplitInitialStates;
 
  bool thePartonsAreJets;
 
};
 
}
 
#endif /* Herwig_HardProcessAnalysis_H */