WeakDecayCurrent.h

// -*- C++ -*-
//
// WeakDecayCurrent.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_WeakDecayCurrent_H
#define HERWIG_WeakDecayCurrent_H
//
// This is the declaration of the WeakDecayCurrent class.
//
#include "ThePEG/Interface/Interfaced.h"
#include "WeakDecayCurrent.fh"
#include "Herwig++/Decay/DecayPhaseSpaceMode.h"
#include "Herwig++/Decay/DecayPhaseSpaceChannel.h"
#include "ThePEG/Helicity/LorentzPolarizationVector.h"

namespace Herwig {
using namespace ThePEG;
using ThePEG::Helicity::LorentzPolarizationVector;
using ThePEG::Helicity::LorentzPolarizationVectorE;

class WeakDecayCurrent: public Interfaced {

public:

  WeakDecayCurrent() : _numbermodes(0) {}

public:

  virtual bool createMode(int icharge,unsigned int imode,DecayPhaseSpaceModePtr mode,
                          unsigned int iloc,unsigned int ires,
                          DecayPhaseSpaceChannelPtr 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(const int imode, const int ichan,Energy & scale, 
          const ParticleVector & decay, DecayIntegrator::MEOption meopt) const=0;

  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:

  static AbstractClassDescription<WeakDecayCurrent> initWeakDecayCurrent;

  WeakDecayCurrent & operator=(const WeakDecayCurrent &);

private:

  vector<int> _quark;

  vector<int> _antiquark;

  unsigned int _numbermodes;

};

}


#include "ThePEG/Utilities/ClassTraits.h"

namespace ThePEG {

template <>
struct BaseClassTrait<Herwig::WeakDecayCurrent,1> {
  typedef Interfaced NthBase;
};

template <>
struct ClassTraits<Herwig::WeakDecayCurrent>
  : public ClassTraitsBase<Herwig::WeakDecayCurrent> {
  static string className() { return "Herwig::WeakDecayCurrent"; }
};

}

#endif /* HERWIG_WeakDecayCurrent_H */