ScalarMesonFactorizedDecayer.h

// -*- C++ -*-
//
// ScalarMesonFactorizedDecayer.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_ScalarMesonFactorizedDecayer_H
#define HERWIG_ScalarMesonFactorizedDecayer_H
//
// This is the declaration of the ScalarMesonFactorizedDecayer class.
//
 
#include "Herwig/Decay/DecayIntegrator.h"
#include "Herwig/Decay/WeakCurrents/WeakCurrent.h"
#include "Herwig/Decay/FormFactors/ScalarFormFactor.h"
#include "ThePEG/StandardModel/StandardModelBase.h"
#include "ThePEG/Helicity/LorentzPolarizationVector.h"
#include "Herwig/Models/StandardModel/StandardCKM.h"
 
namespace Herwig {
using namespace ThePEG;
 
class ScalarMesonFactorizedDecayer: public DecayIntegrator {
 
public:
 
  ScalarMesonFactorizedDecayer();
 
public:
 
  virtual int modeNumber(bool & cc, tcPDPtr parent, 
                         const tPDVector & children) const;
  
  virtual bool accept(tcPDPtr parent, const tPDVector & children) const;
 
  double me2(const int ichan,const Particle & part,
             const tPDVector & outgoing,
             const vector<Lorentz5Momentum> & momenta,
             MEOption meopt) const;
 
  virtual void constructSpinInfo(const Particle & part,
                                 ParticleVector outgoing) const;
 
  virtual void dataBaseOutput(ofstream & os,bool header) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
protected:
 
  virtual void doinit();
 
  virtual void doinitrun();
 
  virtual void rebind(const TranslationMap & trans)
   ;
 
  virtual IVector getReferences();
 
 
private:
 
  void findModes(unsigned int imode, tPDVector & incoming,
                 vector<tPDVector> & outgoing,
                 vector<unsigned int> & loc,vector<bool> & cc);
 
private:
 
  ScalarMesonFactorizedDecayer & operator=(const ScalarMesonFactorizedDecayer &) = delete;
 
private:
 
  vector<WeakCurrentPtr> _current;
 
  vector<ScalarFormFactorPtr> _form;
 
  double _a1b;
 
  double _a2b;
 
  double _a1c;
 
  double _a2c;
 
  vector<vector<unsigned int> > _currentmapA;
 
  vector<vector<unsigned int> > _currentmapB;
 
  vector<vector<unsigned int> > _formmapA;
 
  vector<vector<unsigned int> > _formmapB;
 
  vector<vector<unsigned int> > _formpart;
 
  vector<vector<Complex> > _CKMfact;
 
  vector<int> _wgtloc;
 
  vector<double> _wgtmax;
 
  vector<double> _weights;
 
  Ptr<StandardCKM>::pointer _ckm;
 
  mutable RhoDMatrix _rho;
 
  mutable vector<vector<Helicity::LorentzPolarizationVector> > _vectors;
 
  mutable vector<vector<Helicity::LorentzTensor<double>    > > _tensors;
 
};
 
}
 
#endif /* HERWIG_ScalarMesonFactorizedDecayer_H */