UA5Handler.h

// -*- C++ -*-
//
// UA5Handler.h is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_UA5_H_
#define HERWIG_UA5_H_

#include <ThePEG/Handlers/HadronizationHandler.h>
#include "Herwig++/Hadronization/CluHadConfig.h"
#include <ThePEG/Vectors/LorentzRotation.h>

namespace Herwig {

using namespace ThePEG;

class UA5Handler : public HadronizationHandler {

public:

   UA5Handler();

  static void Init();

public:

  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

public:

  virtual void handle(EventHandler &eh, const tPVector &tagged,
                      const Hint &hint) 
   ;

protected:

  virtual IBPtr clone() const;

  virtual IBPtr fullclone() const;

private:

  void performDecay(PPtr parent,int & totalcharge,int & numbercharge) const;

  void decayCluster(ClusterPtr cluster, bool single) const;

  void insertParticle(PPtr particle,StepPtr step,bool all) const;

  inline double negativeBinomial(int N, double mean, double ek) const;
  
  inline double meanMultiplicity(Energy E) const;
  
  unsigned int multiplicity(Energy E) const;
  
  void generateMomentum(tClusterPtr clu1,tClusterPtr clu2,
                        const ClusterVector &clusters, Energy CME,
                        const Lorentz5Momentum & cm) const;
  
  void generateCylindricalPS(const ClusterVector &clusters, Energy CME) const;
  
  LorentzRotation rotate(const LorentzMomentum &p) const;
  
  template <typename T>
  inline T gaussDistribution(T mean, T stdev) const;
  
  inline double randUng(double A, double B) const;
  
  template <typename T>
  inline void randAzm(T pt, T &px, T &py) const;
  
  inline Energy randExt(Energy AM0,InvEnergy B) const;
  
private:
  
  static ClassDescription<UA5Handler> initUA5Handler;
  
  UA5Handler& operator=(const UA5Handler &);
  
private:

  ClusterFissionerPtr clusterFissioner;
  
  ClusterDecayerPtr   clusterDecayer;
  
  double _n1;
  
  double _n2;
  
  double _n3;
  
  double _k1;
  
  double _k2;
  
  Energy _m1;
  
  InvEnergy _m2;
  
  InvEnergy _p1;
  
  InvEnergy _p2;
  
  InvEnergy _p3;
  
  double _probSoft;
  
  double _enhanceCM; 
  
  unsigned int _maxtries;
  

  bool _needWarning;
};

}

namespace ThePEG {

template<>
struct BaseClassTrait<Herwig::UA5Handler,1> { 
  typedef HadronizationHandler NthBase;
};

template<>
struct ClassTraits<Herwig::UA5Handler> :
  public ClassTraitsBase<Herwig::UA5Handler> {
    static string className() { return "Herwig::UA5Handler"; }
    static string library() { return "HwUA5.so"; }
};

}

#include "UA5Handler.icc"

#endif