PartonBinInstance.h

// -*- C++ -*-
//
// PartonBinInstance.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// 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_PartonBinInstance_H
#define THEPEG_PartonBinInstance_H
// This is the declaration of the PartonBinInstance class.
 
#include "ThePEG/Config/ThePEG.h"
#include "ThePEG/PDF/PartonBin.h"
 
namespace ThePEG {
 
ThePEG_DECLARE_CLASS_POINTERS(PartonBinInstance,PBIPtr);
typedef pair<PBIPtr,PBIPtr> PBIPair;
 
ThePEG_DECLARE_CLASS_POINTERS(RemInfoBase,RemIPtr);
 
class PartonBinInstance: public PersistentBase {
 
public:
 
  PartonBinInstance();
 
  PartonBinInstance(const PartonBinInstance &);
 
  virtual ~PartonBinInstance();
 
  PartonBinInstance(tcPBPtr, tPBIPtr = tPBIPtr());
 
  PartonBinInstance(tPPtr parton, tcPBPtr pb, Energy2 scale = ZERO);
 
  PartonBinInstance(tPPtr particle, tPPtr parton, tcPBPtr pb, 
                    Energy2 scale = ZERO);
 
 
public:
 
  tcPBPtr bin() const { return theBin; }
 
  const PartonVector & bins() const { return theBins; }
 
  tcPDPtr particleData() const { return bin()->particle(); }
 
  tcPDPtr partonData() const { return bin()->parton(); }
 
  tPBIPtr incoming() const { return theIncoming; }
 
  tPBIPtr getFirst();
 
  tcPDFPtr pdf() const { return bin()->pdf(); }
 
  tcRemHPtr remnantHandler() const { return bin()->remnantHandler(); }
 
  bool hasPoleIn1() const;
 
  void reset(double lx = 0, Energy2 Q2 = ZERO);
 
  void prepare();
 
  void generate(const double * r);
 
  double jacobian() const { return theJacobian; }
 
  void jacobian(double j) { theJacobian = j; }
 
  tPPtr particle() const { return theParticle; }
 
  void particle(tPPtr p) { theParticle = p; }
 
  tPPtr parton() const { return theParton; }
 
  void parton(tPPtr p) { theParton = p; }
 
  const PVector & partons() const { return thePartons; }
 
  double xi() const {
    if ( theXi < 0.0 ) theXi = exp(-li());
    return theXi;
  }
 
 
  double eps() const {
    if ( theEps < 0.0 ) theEps =  Math::exp1m(-li());
    return theEps;
  }
 
  double li() const { return theLi; }
 
  void li(double lx) {
    theLi = lx;
    theXi = theEps = -1.0;
  }
 
 
  double x() const {
    if ( theX < 0.0 ) theX = exp(-l());
    return theX;
  }
 
 
  double l() const { return theL; }
 
  void l(double lx) {
    theL = lx;
    theX = -1.0;
  }
 
 
  Energy2 scale() const { return theScale; }
  
 
  void scale(Energy2 s) { theScale = s; }
 
  const TransverseMomentum & kT() const { return theKT; }
 
  double remnantWeight() const { return theRemnantWeight; }
 
  void remnantWeight(double w) { theRemnantWeight = w; }
 
  const PVector & remnants() const { return theRemnants; }
 
  void remnants(const PVector & rems) { theRemnants = rems; }
 
  tRemIPtr remnantInfo() const { return theRemInfo; }
 
  void remnantInfo(tRemIPtr ri) { theRemInfo = ri; }
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
private:
 
  cPBPtr theBin;
 
  PartonVector theBins;
 
  PBIPtr theIncoming;
 
  double theJacobian;
 
  PPtr theParticle;
 
  PPtr theParton;
 
  PVector thePartons;
 
  mutable double theXi;
  mutable double theEps;
  double theLi;
 
  mutable double theX;
  double theL;
 
  Energy2 theScale;
 
  TransverseMomentum theKT;
 
  double theRemnantWeight;
 
  PVector theRemnants;
 
  RemIPtr theRemInfo;
 
private:
 
  static ClassDescription<PartonBinInstance> initPartonBinInstance;
 
  PartonBinInstance & operator=(const PartonBinInstance &) = delete;
 
};
 
class RemInfoBase: public Base {
public:
  virtual ~RemInfoBase() {}
};
 
}
 
 
namespace ThePEG {
 
template <>
struct BaseClassTrait<PartonBinInstance,1>: public ClassTraitsType {
  typedef Base NthBase;
};
 
template <>
struct ClassTraits<PartonBinInstance>:
    public ClassTraitsBase<PartonBinInstance> {
  static string className() { return "ThePEG::PartonBinInstance"; }
};
 
}
 
#endif /* THEPEG_PartonBinInstance_H */