WeakCurrent.h

// -*- C++ -*-
//
// WeakCurrent.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_WeakCurrent_H
#define HERWIG_WeakCurrent_H
//
// This is the declaration of the WeakCurrent class.
//
#include "ThePEG/Interface/Interfaced.h"
#include "WeakCurrent.fh"
#include "Herwig/Decay/DecayIntegrator.h"
#include "Herwig/Decay/PhaseSpaceMode.h"
#include "Herwig/Decay/PhaseSpaceChannel.h"
#include "ThePEG/Helicity/LorentzPolarizationVector.h"
#include "Herwig/Decay/IsoSpin.h"
 
namespace Herwig {
using namespace ThePEG;
using ThePEG::Helicity::LorentzPolarizationVector;
using ThePEG::Helicity::LorentzPolarizationVectorE;
  
class WeakCurrent: public Interfaced {
 
public:
 
  WeakCurrent() : numberModes_(0) {}
 
public:
 
  virtual bool createMode(int icharge, tcPDPtr resonance,
                          FlavourInfo flavour,
                          unsigned int imode,PhaseSpaceModePtr mode,
                          unsigned int iloc,int ires,
                          PhaseSpaceChannel phase, Energy upp )=0;
 
  virtual tPDVector particles(int icharge, unsigned int imode, int iq, int ia)=0;
 
  unsigned int numberOfModes() const {return quark_.size();}
 
  virtual vector<LorentzPolarizationVectorE> 
  current(tcPDPtr resonance,
          FlavourInfo flavour,
          const int imode, const int ichan,Energy & scale,
          const tPDVector & outgoing,
          const vector<Lorentz5Momentum> & momenta,
          DecayIntegrator::MEOption meopt) const=0;
 
  virtual void constructSpinInfo(ParticleVector decay) const;
 
  virtual bool accept(vector<int> id)=0;
 
  virtual unsigned int decayMode(vector<int> id)=0;
 
  void decayModeInfo(unsigned int imode, int& iq, int& ia) const {
    if(imode<quark_.size()) {
      iq=quark_[imode];
      ia=antiquark_[imode];
    }
    else {
      iq=0;
      ia=0;
    }
  }
 
  virtual void dataBaseOutput(ofstream & os,bool header,bool create) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  void addDecayMode(int iq,int ia) {
    quark_.push_back(iq);
    antiquark_.push_back(ia);
  }
 
  void setInitialModes(unsigned int nmodes) {numberModes_=nmodes;}
 
private:
 
  WeakCurrent & operator=(const WeakCurrent &) = delete;
 
private:
 
  vector<int> quark_;
 
  vector<int> antiquark_;
 
  unsigned int numberModes_;
 
};
 
}
 
 
#endif /* HERWIG_WeakCurrent_H */