EIC code displayed by LXR master/​include/​Pythia8/​SigmaCompositeness.h	Source navigation Diff markup Identifier search General search
	Version: master test Revert   Change

File indexing completed on 2025-01-18 10:06:29

 // SigmaCompositeness.h is a part of the PYTHIA event generator.
 // Copyright (C) 2024 Torbjorn Sjostrand.
 // PYTHIA is licenced under the GNU GPL v2 or later, see COPYING for details.
 // Please respect the MCnet Guidelines, see GUIDELINES for details.
 
 // Header file for compositiness-process differential cross sections.
 // Contains classes derived from SigmaProcess via Sigma(1/2)Process.
 
 #ifndef Pythia8_SigmaCompositeness_H
 #define Pythia8_SigmaCompositeness_H
 
 #include "Pythia8/SigmaProcess.h"
 
 namespace Pythia8 {
 
 //==========================================================================
 
 // A derived class for q g -> q^* (excited quark state).
 
 class Sigma1qg2qStar : public Sigma1Process {
 
 public:
 
   // Constructor.
   Sigma1qg2qStar(int idqIn) : idq(idqIn), idRes(), codeSave(), mRes(),
     GammaRes(), m2Res(), GamMRat(), Lambda(), coupFcol(), widthIn(),
     sigBW(), qStarPtr() {}
 
   // Initialize process.
   virtual void initProc();
 
   // Calculate flavour-independent parts of cross section.
   virtual void sigmaKin();
 
   // Evaluate sigmaHat(sHat).
   virtual double sigmaHat();
 
   // Select flavour, colour and anticolour.
   virtual void setIdColAcol();
 
   // Evaluate weight for q* decay angles (else inactive).
   virtual double weightDecay( Event& process, int iResBeg, int iResEnd);
 
   // Info on the subprocess.
   virtual string name()       const {return nameSave;}
   virtual int    code()       const {return codeSave;}
   virtual string inFlux()     const {return "qg";}
   virtual int    resonanceA() const {return idRes;}
 
 private:
 
   // Parameters set at initialization or for current kinematics.
   int    idq, idRes, codeSave;
   string nameSave;
   double mRes, GammaRes, m2Res, GamMRat, Lambda, coupFcol, widthIn, sigBW;
 
   // Pointer to properties of the particle species, to access decay channels.
   ParticleDataEntryPtr qStarPtr;
 
 };
 
 //==========================================================================
 
 // A derived class for l gamma -> l^* (excited lepton state).
 
 class Sigma1lgm2lStar : public Sigma1Process {
 
 public:
 
   // Constructor.
   Sigma1lgm2lStar(int idlIn) : idl(idlIn), idRes(), codeSave(), mRes(),
     GammaRes(), m2Res(), GamMRat(), Lambda(), coupChg(), widthIn(), sigBW(),
     qStarPtr() {}
 
   // Initialize process.
   virtual void initProc();
 
   // Calculate flavour-independent parts of cross section.
   virtual void sigmaKin();
 
   // Evaluate sigmaHat(sHat).
   virtual double sigmaHat();
 
   // Select flavour, colour and anticolour.
   virtual void setIdColAcol();
 
   // Evaluate weight for l* decay angles (else inactive).
   virtual double weightDecay( Event& process, int iResBeg, int iResEnd);
 
   // Info on the subprocess.
   virtual string name()       const {return nameSave;}
   virtual int    code()       const {return codeSave;}
   virtual string inFlux()     const {return "fgm";}
   virtual int    resonanceA() const {return idRes;}
 
 private:
 
   // Parameters set at initialization or for current kinematics.
   int    idl, idRes, codeSave;
   string nameSave;
   double mRes, GammaRes, m2Res, GamMRat, Lambda, coupChg, widthIn, sigBW;
 
   // Pointer to properties of the particle species, to access decay channels.
   ParticleDataEntryPtr qStarPtr;
 
 };
 
 //==========================================================================
 
 // A derived class for q q' -> q^* q' (excited quark state).
 
 class Sigma2qq2qStarq : public Sigma2Process {
 
 public:
 
   // Constructor.
   Sigma2qq2qStarq(int idqIn) : idq(idqIn), idRes(), codeSave(), Lambda(),
     preFac(), openFracPos(), openFracNeg(), sigmaA(), sigmaB() {}
 
   // Initialize process.
   virtual void initProc();
 
   // Calculate flavour-independent parts of cross section.
   virtual void sigmaKin();
 
   // Evaluate sigmaHat(sHat).
   virtual double sigmaHat();
 
   // Select flavour, colour and anticolour.
   virtual void setIdColAcol();
 
   // Evaluate weight for q* decay angles (else inactive).
   virtual double weightDecay(Event& process, int iResBeg, int iResEnd);
 
   // Info on the subprocess.
   virtual string name()       const {return nameSave;}
   virtual int    code()       const {return codeSave;}
   virtual string inFlux()     const {return "qq";}
   virtual int    id3Mass()    const {return idRes;}
 
 private:
 
   // Parameters set at initialization or for current kinematics.
   int    idq, idRes, codeSave;
   string nameSave;
   double Lambda, preFac, openFracPos, openFracNeg, sigmaA, sigmaB;
 
 };
 
 //==========================================================================
 
 // A derived class for q qbar -> l^* lbar (excited lepton state).
 
 class Sigma2qqbar2lStarlbar : public Sigma2Process {
 
 public:
 
   // Constructor.
   Sigma2qqbar2lStarlbar(int idlIn) : idl(idlIn), idRes(), codeSave(), Lambda(),
     preFac(), openFracPos(), openFracNeg(), sigma() {}
 
   // Initialize process.
   virtual void initProc();
 
   // Calculate flavour-independent parts of cross section.
   virtual void sigmaKin();
 
   // Evaluate sigmaHat(sHat).
   virtual double sigmaHat() {return sigma;}
 
   // Select flavour, colour and anticolour.
   virtual void setIdColAcol();
 
   // Evaluate weight for l* decay angles (else inactive).
   virtual double weightDecay(Event& process, int iResBeg, int iResEnd);
 
   // Info on the subprocess.
   virtual string name()       const {return nameSave;}
   virtual int    code()       const {return codeSave;}
   virtual string inFlux()     const {return "qqbarSame";}
   virtual int    id3Mass()    const {return idRes;}
 
 private:
 
   // Parameters set at initialization or for current kinematics.
   int    idl, idRes, codeSave;
   string nameSave;
   double Lambda, preFac, openFracPos, openFracNeg, sigma;
 
 };
 
 //==========================================================================
 
 // A derived class for q qbar -> lStar lStarBar.
 // Code contributed by Olga Igonkina.
 
 class Sigma2qqbar2lStarlStarBar: public Sigma2Process {
 
 public:
 
   // Constructor.
   Sigma2qqbar2lStarlStarBar(int idlIn) : idl(idlIn), idRes(), codeSave(),
     Lambda(), preFac(), openFracPos(), openFracNeg(), sigma() {}
 
   // Initialize process.
   void initProc();
 
   // Calculate flavour-independent parts of cross section.
   void sigmaKin();
 
   // Evaluate sigmaHat(sHat).
   virtual double sigmaHat() {return sigma;}
 
   // Select flavour, colour and anticolour.
   void setIdColAcol();
 
   // Evaluate weight for l* decay angles (else inactive).
   virtual double weightDecay(Event& process, int iResBeg, int iResEnd);
 
   // Info on the subprocess.
   virtual string name()       const {return nameSave;}
   virtual int    code()       const {return codeSave;}
   virtual string inFlux()     const {return "qqbarSame";}
   virtual int    id3Mass()    const {return idRes;}
   virtual int    id4Mass()    const {return idRes;}
 
 private:
 
   // Parameters set at initialization or for current kinematics.
   int    idl, idRes, codeSave;
   string nameSave;
   double Lambda, preFac, openFracPos, openFracNeg, sigma;
 
 };
 
 //==========================================================================
 
 // A derived class for q q -> q q (quark contact interactions).
 // Based on, Sigma2qq2qq (QCD).
 
 class Sigma2QCqq2qq : public Sigma2Process {
 
 public:
 
   // Constructor.
   Sigma2QCqq2qq() : sigT(), sigU(), sigTU(), sigST(), sigSum(), sigQCSTU(),
     sigQCUTS(), qCLambda2(), qCetaLL(), qCetaRR(), qCetaLR(){}
 
   // Initialize process.
   virtual void initProc();
 
   // Calculate flavour-independent parts of cross section.
   virtual void sigmaKin();
 
   // Evaluate d(sigmaHat)/d(tHat).
   virtual double sigmaHat();
 
   // Select flavour, colour and anticolour.
   virtual void setIdColAcol();
 
   // Info on the subprocess.
   virtual string name()   const {return "q q(bar)' -> (QC) -> q q(bar)'";}
   virtual int    code()   const {return 4201;}
   virtual string inFlux() const {return "qq";}
 
  private:
 
   // Values stored for colour flow selection.
   double sigT, sigU, sigTU, sigST, sigSum, sigQCSTU, sigQCUTS;
 
   // Compositeness parameters.
   double qCLambda2;
   int    qCetaLL, qCetaRR, qCetaLR;
 
 };
 
 //==========================================================================
 
 // A derived class for q qbar -> q' qbar' (quark contact interactions).
 // Based on, Sigma2qqbar2qqbarNew(QCD).
 // Note: This process give the same contributions for q == q' and q != q'.
 
 class Sigma2QCqqbar2qqbar : public Sigma2Process {
 
 public:
 
   // Constructor.
   Sigma2QCqqbar2qqbar() : qCnQuarkNew(), idNew(), mNew(), m2New(), sigS(),
     sigma(), qCLambda2(), qCetaLL(), qCetaRR(), qCetaLR(){}
 
   // Initialize process.
   virtual void initProc();
 
   // Calculate flavour-independent parts of cross section.
   virtual void sigmaKin();
 
   // Evaluate d(sigmaHat)/d(tHat).
   virtual double sigmaHat() {return sigma;}
 
   // Select flavour, colour and anticolour.
   virtual void setIdColAcol();
 
   // Info on the subprocess.
   virtual string name()   const {return "q qbar -> (QC) -> q' qbar' (uds)";}
   virtual int    code()   const {return 4202;}
   virtual string inFlux() const {return "qqbarSame";}
 
  private:
 
   // Number of outgoing quark flavours to be considered, given that
   // matrix elements are calculated in the massless approximation.
   int    qCnQuarkNew;
 
   // Values stored for colour flow selection.
   int    idNew;
   double mNew, m2New, sigS, sigma;
 
   // Compositeness parameters.
   double qCLambda2;
   int    qCetaLL, qCetaRR, qCetaLR;
 
 };
 
 //==========================================================================
 
 // A derived class for f fbar -> l lbar
 // (contact interactions).
 // Does not include t-channel contributions relevant for e^+e^- to e^+e^-
 
 class Sigma2QCffbar2llbar : public Sigma2Process {
 
 public:
 
   // Constructor: bool Graviton  = true, to use LED graviton settings.
   Sigma2QCffbar2llbar (int idIn, int codeIn) : idNew(idIn), codeNew(codeIn),
     qCmNew(), qCmNew2(), qCmZ(), qCmZ2(), qCGZ(), qCGZ2(), sigma0(),
     qCLambda2(), qCetaLL(), qCetaRR(), qCetaLR(), qCetaRL(), qCPropGm(),
     qCrePropZ(), qCimPropZ() {}
 
   // Initialize process.
   virtual void initProc();
 
   // Calculate flavour-independent parts of cross section;
   // first step when inflavours unknown.
   virtual void sigmaKin();
 
   // Evaluate sigmaHat(sHat); second step for given inflavours.
   virtual double sigmaHat();
 
   // Select flavour, colour and anticolour.
   virtual void setIdColAcol();
 
   // Info on the subprocess.
   virtual string name()       const {return nameNew;}
   virtual int    code()       const {return codeNew;}
   virtual string inFlux()     const {return "ffbarSame";}
   virtual bool   isSChannel() const {return true;}
 
 private:
 
   // Process values.
   string nameNew;
   int    idNew, codeNew;
   double qCmNew, qCmNew2, qCmZ, qCmZ2, qCGZ, qCGZ2, sigma0;
 
   // Compositeness parameters.
   double qCLambda2;
   int    qCetaLL, qCetaRR, qCetaLR, qCetaRL;
   double qCPropGm, qCrePropZ, qCimPropZ;
 
 };
 
 //==========================================================================
 
 } // end namespace Pythia8
 
 #endif // Pythia8_SigmaCompositeness_H

This page was automatically generated by the 2.3.7 LXR engine. The LXR team