HJetsAnalysis.h

// -*- C++ -*-
#ifndef Herwig_HJetsAnalysis_H
#define Herwig_HJetsAnalysis_H
//
// This is the declaration of the HJetsAnalysis class.
//
 
#include "Herwig/Analysis/JetsPlusAnalysis.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class HJetsAnalysis: public Herwig::JetsPlusAnalysis {
 
public:
 
  HJetsAnalysis();
 
public:
 
  virtual void reconstructHardObjects(ParticleVector&);
 
protected:
 
  virtual void analyzeSpecial(long id, double weight);
 
  virtual void finalize(XML::Element&);
 
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:
 
  Statistics::Histogram theHiggsYStar;
 
  Statistics::Histogram theThirdJetYStar;
 
  Statistics::Histogram theFourthJetYStar;
 
  Statistics::Histogram theJet12HiggsDeltaPhi;
 
  Statistics::Histogram theJeppeDeltaPhi;
 
protected:
 
  double yStar(const LorentzMomentum& jet1, const LorentzMomentum& jet2,
               const LorentzMomentum& obj) const {
    double y1 = jet1.rapidity();
    double y2 = jet2.rapidity();
    double res =
      obj.rapidity() - 0.5*(y1+y2);
    return res/(y1-y2);
  }
 
  Statistics::Histogram& higgsYStar() {
    if ( !theHiggsYStar.bins().empty() )
      return theHiggsYStar;
    return theHiggsYStar =
      Statistics::Histogram("HiggsYStar",Statistics::Histogram::regularBinEdges(-6,6,120),false,false);
  }
 
  Statistics::Histogram& thirdJetYStar() {
    if ( !theThirdJetYStar.bins().empty() )
      return theThirdJetYStar;
    return theThirdJetYStar =
      Statistics::Histogram("Jet3YStar",Statistics::Histogram::regularBinEdges(-6,6,120),false,false);
  }
 
  Statistics::Histogram& fourthJetYStar() {
    if ( !theFourthJetYStar.bins().empty() )
      return theFourthJetYStar;
    return theFourthJetYStar =
      Statistics::Histogram("Jet4YStar",Statistics::Histogram::regularBinEdges(-6,6,120),false,false);
  }
 
  Statistics::Histogram& jet12HiggsDeltaPhi() {
    if ( !theJet12HiggsDeltaPhi.bins().empty() )
      return theJet12HiggsDeltaPhi;
    return theJet12HiggsDeltaPhi =
      Statistics::Histogram("Jet12hDeltaPhi",
                            Statistics::Histogram::regularBinEdges(-Constants::pi,Constants::pi,32),
                            make_pair(-Constants::pi,Constants::pi));
  }
 
  Statistics::Histogram& jeppeDeltaPhi() {
    if ( !theJeppeDeltaPhi.bins().empty() )
      return theJeppeDeltaPhi;
    return theJeppeDeltaPhi =
      Statistics::Histogram("JeppeDeltaPhi",
                            Statistics::Histogram::regularBinEdges(-Constants::pi,Constants::pi,32),
                            make_pair(-Constants::pi,Constants::pi));
  }
 
private:
 
  HJetsAnalysis & operator=(const HJetsAnalysis &) = delete;
 
};
 
}
 
#endif /* Herwig_HJetsAnalysis_H */