Cluster.h

// -*- C++ -*-
//
// Cluster.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_Cluster_H
#define HERWIG_Cluster_H
 
#include <ThePEG/EventRecord/Particle.h>
#include "Herwig/Utilities/EnumParticles.h"
#include "CluHadConfig.h"
#include "ClusterHadronizationHandler.fh"
#include "Cluster.fh"
#include "ThePEG/Utilities/ClassDescription.h"
 
namespace Herwig {
using namespace ThePEG;
 
class Cluster : public Particle {
  
protected:
 
  Cluster();
 
  friend struct ClassTraits<Cluster>;
 
public:
  Cluster(tcEventPDPtr);
  
  Cluster(tPPtr part1, tPPtr part2, tPPtr part3 = tPPtr());    
  
  Cluster(const Particle &);
  
  unsigned int numComponents() const
  { return _numComp; }
  
  Energy sumConstituentMasses() const;
    
  tPPtr particle(int i) const;
 
  tPPtr particleB(int i) const;
  
  tPPtr colParticle(bool anti = false) const;
 
  tPPtr antiColParticle() const;
  
  bool isPerturbative(int) const;
  
  bool isBeamRemnant(int) const;
  
  void setBeamRemnant(int,bool);
  
  int clusterId() const { return _id; }
 
public:
 
  bool isBeamCluster() const;
  
  void flagAsReshuffled()
  { _hasReshuffled = true; }
 
  void isBeamCluster(tPPtr part);
  
  bool isAvailable() const
  { return _isAvailable; }
  
  void isAvailable(bool inputAvailable)
  { _isAvailable = inputAvailable; }
  
  bool isStatusInitial() const
  { return parents().empty(); }
  
   bool isReadyToDecay() const
  { return children().empty(); }
  
  bool isRedefined() const { 
    return ( children().size() == 1 
             && children()[0]->id() == ParticleID::Cluster );
  }
  
  bool hasBeenReshuffled() const
  { return _hasReshuffled; }
  
  bool isStatusFinal() const;
 
public:
 
  static LorentzPoint calculateX(tPPtr q1, tPPtr q2);
  
protected:
  
  virtual PPtr clone() const;
 
  virtual PPtr fullclone() const;
 
public:
 
  void persistentOutput(PersistentOStream &) const;
 
  void persistentInput(PersistentIStream &, int);
 
private:   
  Cluster & operator=(const Cluster &) = delete;
  
  void calculateP();
  
  bool initPerturbative(tPPtr p);
 
  static ClassDescription<Cluster> initCluster; 
 
 
 
  bool        _isAvailable;        
  bool        _hasReshuffled;      
  ParticleVector _component;       
  tParticleVector _original;        
  vector<bool> _isBeamRemnant;     
  vector<bool> _isPerturbative;    
  unsigned int _numComp;                    
  long _id;                        
};
  
} // end namespace Herwig  
 
#include "ThePEG/Utilities/ClassTraits.h"
 
namespace ThePEG {
 
template <>
struct BaseClassTrait<Herwig::Cluster,1> {
  typedef Particle NthBase;
};
 
template <>
struct ClassTraits<Herwig::Cluster>:
  public ClassTraitsBase<Herwig::Cluster> {
  static string className() { return "Herwig::Cluster"; }
  static TPtr create() { return TPtr::Create(Herwig::Cluster()); }
};
 
}
 
 
 
 
#endif // HERWIG_Cluster_H