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 */