JetRegion.h

// -*- C++ -*-
//
// JetRegion.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
// Copyright (C) 2009-2019 Simon Platzer
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ThePEG_JetRegion_H
#define ThePEG_JetRegion_H
//
// This is the declaration of the JetRegion class.
//
 
#include "ThePEG/Handlers/HandlerBase.h"
#include "ThePEG/EventRecord/Particle.h"
#include "ThePEG/Cuts/Cuts.h"
 
namespace ThePEG {
 
class JetRegion: public HandlerBase {
 
public:
 
  JetRegion();
 
public:
 
  Energy ptMin() const { return thePtMin; }
 
  Energy ptMax() const { return thePtMax; }
 
  const vector<pair<double,double> >& yRanges() const { return theYRanges; }
 
  const vector<int>& accepts() const { return theAccepts; }
 
  bool fuzzy() const { return theFuzzy; }
 
  double cutWeight() const { return theCutWeight; }
 
  bool lessThanEnergy(Energy a, Energy b, double& weight) const;
 
  bool lessThanRapidity(double a, double b, double& weight) const;
 
public:
 
  virtual void describe() const;
 
  virtual bool matches(tcCutsPtr parent, int n, const LorentzMomentum& p,
                       double yHat = 0.0);
 
  bool didMatch() { return theDidMatch; }
 
  virtual void reset() { theDidMatch = false; }
 
  int lastNumber() const { return theLastNumber; }
 
  const LorentzMomentum& lastMomentum() const { return theLastMomentum; }
 
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:
 
  string doYRange(string);
 
  Energy thePtMin;
 
  Energy thePtMax;
 
  vector<pair<double,double> > theYRanges;
 
  vector<int> theAccepts;
 
  bool theDidMatch;
 
  int theLastNumber;
 
  LorentzMomentum theLastMomentum;
 
  bool theFuzzy;
 
  double theCutWeight;
 
  Energy theEnergyCutWidth;
 
  double theRapidityCutWidth;
 
  JetRegion & operator=(const JetRegion &) = delete;
 
};
 
}
 
#endif /* ThePEG_JetRegion_H */