MatchboxAmplitudelnuqqbar.h

// -*- C++ -*-
//
// MatchboxAmplitudelnuqqbar.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_MatchboxAmplitudelnuqqbar_H
#define Herwig_MatchboxAmplitudelnuqqbar_H
//
// This is the declaration of the MatchboxAmplitudelnuqqbar class.
//
 
#include "Herwig/MatrixElement/Matchbox/Base/MatchboxAmplitude.h"
#include "Herwig/MatrixElement/Matchbox/Builtin/Amplitudes/MatchboxCurrents.h"
#include "ThePEG/StandardModel/StandardModelBase.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/EventRecord/Particle.h"
#include "ThePEG/Repository/UseRandom.h"
#include "ThePEG/Repository/EventGenerator.h"
#include "ThePEG/Utilities/DescribeClass.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/PDT/EnumParticles.h"
#include "Herwig/MatrixElement/Matchbox/Utility/SU2Helper.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MatchboxAmplitudelnuqqbar: public MatchboxAmplitude, public MatchboxCurrents {
 
public:
 
  MatchboxAmplitudelnuqqbar();
 
  virtual unsigned int orderInGs() const { return 0; }
 
  virtual unsigned int orderInGem() const { return 2; }
 
  virtual bool haveOneLoop() const { return true; }
  
  virtual bool canHandle(const PDVector&) const;
  
  virtual void prepareAmplitudes(Ptr<MatchboxMEBase>::tcptr);
  
  virtual Complex evaluate(size_t, const vector<int>&, Complex&);
  
  virtual Complex evaluateOneLoop(size_t, const vector<int>&);
  
  virtual bool isBDK() const { return true; }
  
  virtual Energy2 mu2() const { return lastSHat(); }
 
  virtual void flushCaches() {
    MatchboxCurrents::reset();
    MatchboxAmplitude::flushCaches();
  }
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
  
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
  
  virtual void doinit();
 
  virtual void doinitrun();
  
private:
  
  bool theDiagonal;
  
  vector< vector<Complex> > theCKM;
 
  MatchboxAmplitudelnuqqbar & operator=(const MatchboxAmplitudelnuqqbar &) = delete;
 
};
 
}
 
#endif /* Herwig_MatchboxAmplitudelnuqqbar_H */