+++ /dev/null
-// SigmaSUSY.h is a part of the PYTHIA event generator.
-// Copyright (C) 2012 Torbjorn Sjostrand.
-// Main authors of this file: N. Desai, P. Skands
-// PYTHIA is licenced under the GNU GPL version 2, see COPYING for details.
-// Please respect the MCnet Guidelines, see GUIDELINES for details.
-
-// Header file for Supersymmetric process differential cross sections.
-// Contains classes derived from SigmaProcess via Sigma2Process.
-
-#ifndef Pythia8_SigmaSUSY_H
-#define Pythia8_SigmaSUSY_H
-
-#include "PhaseSpace.h"
-#include "PythiaComplex.h"
-#include "SigmaProcess.h"
-#include "SusyCouplings.h"
-#include "SusyLesHouches.h"
-
-namespace Pythia8 {
-
-//==========================================================================
-
-// A derived class for q qbar -> neutralino_i neutralino_j.
-
-class Sigma2qqbar2chi0chi0 : public Sigma2Process {
-
-public:
-
- // Constructor.
- Sigma2qqbar2chi0chi0() {};
-
- // Constructor.
- Sigma2qqbar2chi0chi0(int id3chiIn, int id4chiIn, int codeIn) {
-
- // Save ordering indices and process code
- id3chi = id3chiIn;
- id4chi = id4chiIn;
- codeSave = codeIn;
-
-
- // Construct id codes from ordering indices.
- id3 = 1000022;
- if (id3chi == 2) id3 = 1000023;
- if (id3chi == 3) id3 = 1000025;
- if (id3chi == 4) id3 = 1000035;
- if (id3chi == 5) id3 = 1000045;
- id4 = 1000022;
- if (id4chi == 2) id4 = 1000023;
- if (id4chi == 3) id4 = 1000025;
- if (id4chi == 4) id4 = 1000035;
- if (id4chi == 5) id4 = 1000045;
-
- }
-
- // 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 nameSave;}
- virtual int code() const {return codeSave;}
- virtual string inFlux() const {return "qq";}
- virtual int id3Mass() const {return abs(id3);}
- virtual int id4Mass() const {return abs(id4);}
- virtual int resonanceA() const {return 23;}
- virtual bool isSUSY() const {return true;}
-
- protected:
-
- // Basic process information
- int id3chi, id4chi, codeSave;
- string nameSave;
-
- // Values stored for later use
- double sigma0, ui, uj, ti, tj, openFracPair;
- complex propZ;
-
- CoupSUSY* coupSUSYPtr;
-
-};
-
-//==========================================================================
-
-// A derived class for q qbar -> neutralino_i chargino_j.
-
-class Sigma2qqbar2charchi0 : public Sigma2qqbar2chi0chi0 {
-
-public:
-
- // Constructor.
- Sigma2qqbar2charchi0(int id3chiIn, int id4chiIn, int codeIn) {
-
- // Save ordering indices and process code
- id3chi = id3chiIn;
- id4chi = id4chiIn;
- codeSave = codeIn;
-
- // Construct id codes from ordering indices.
- id3 = (abs(id3chi) == 2) ? 1000037 : 1000024;
- if (id3chi < 0) id3 = -id3;
-
- id4 = 1000022;
- if (id4chi == 2) id4 = 1000023;
- if (id4chi == 3) id4 = 1000025;
- if (id4chi == 4) id4 = 1000035;
- if (id4chi == 5) id4 = 1000045;
-
- }
-
- // Calculate flavour-independent parts of cross section.
- virtual void sigmaKin();
-
- // Evaluate d(sigmaHat)/d(tHat).
- virtual double sigmaHat();
-
- virtual int resonanceA() const {return 24;}
-
-protected :
-
- complex propW;
-
-};
-
-//==========================================================================
-
-// A derived class for q qbar -> chargino+_i chargino-_j.
-
-class Sigma2qqbar2charchar : public Sigma2qqbar2chi0chi0 {
-
-public:
-
- // Constructor.
- Sigma2qqbar2charchar(int id3chiIn, int id4chiIn, int codeIn) {
-
- // Save ordering indices and process code
- id3chi = id3chiIn;
- id4chi = id4chiIn;
- codeSave = codeIn;
-
- // Construct id codes from ordering indices.
- id3 = (abs(id3chi) == 2) ? 1000037 : 1000024;
- id4 = (abs(id4chi) == 2) ? -1000037 : -1000024;
-
- }
-
- // Calculate flavour-independent parts of cross section.
- virtual void sigmaKin();
-
- // Evaluate d(sigmaHat)/d(tHat).
- virtual double sigmaHat();
-
-};
-
-//==========================================================================
-
-// A derived class for q g -> neutralino_i squark_j (and cc)
-
-class Sigma2qg2chi0squark : public Sigma2Process {
-
-public:
-
- // Constructor.
- Sigma2qg2chi0squark() { };
-
- // Constructor.
- Sigma2qg2chi0squark(int id3chiIn, int id4sqIn, bool isUp, int codeIn) {
-
- // Save ordering indices and process code
- id3chi = id3chiIn;
- id4sq = id4sqIn;
- codeSave = codeIn;
-
- // Construct id codes from ordering indices.
- id3 = 1000022;
- if (id3chi == 2) id3 = 1000023;
- if (id3chi == 3) id3 = 1000025;
- if (id3chi == 4) id3 = 1000035;
- if (id3chi == 5) id3 = 1000045;
- id4 = 1000001 + (isUp ? 1 : 0);
- if (id4sq == 2) id4 = 1000003 + (isUp ? 1 : 0);
- if (id4sq == 3) id4 = 1000005 + (isUp ? 1 : 0);
- if (id4sq == 4) id4 = 2000001 + (isUp ? 1 : 0);
- if (id4sq == 5) id4 = 2000003 + (isUp ? 1 : 0);
- if (id4sq == 6) id4 = 2000005 + (isUp ? 1 : 0);
-
- }
-
- // 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 nameSave;}
- virtual int code() const {return codeSave;}
- virtual string inFlux() const {return "qg";}
- virtual int id3Mass() const {return abs(id3);}
- virtual int id4Mass() const {return abs(id4);}
- virtual bool isSUSY() const {return true;}
-
- protected:
-
- // Basic process information
- int id3chi, id4sq, codeSave;
- string nameSave;
-
- // Values stored for later use
- double sigma0, ui, uj, ti, tj, openFracPair;
-
- //SUSY couplings
- CoupSUSY* coupSUSYPtr;
-
-};
-
-//==========================================================================
-
-// A derived class for q g -> chargino_i squark_j (incl cc)
-
-class Sigma2qg2charsquark : public Sigma2qg2chi0squark {
-
-public:
-
- // Constructor.
- Sigma2qg2charsquark(int id3chiIn, int id4sqIn, bool isUp, int codeIn) {
-
- // Save ordering indices and process code
- id3chi = id3chiIn;
- id4sq = id4sqIn;
- codeSave = codeIn;
-
- // Construct id codes from ordering indices.
- id3Sav = 1000024;
- if (abs(id3chi) == 2) id3Sav = 1000037;
- if (isUp) id3Sav = -id3Sav;
- id4Sav = 1000001 + (isUp ? 1 : 0);
- if (id4sq == 2) id4Sav = 1000003 + (isUp ? 1 : 0);
- if (id4sq == 3) id4Sav = 1000005 + (isUp ? 1 : 0);
- if (id4sq == 4) id4Sav = 2000001 + (isUp ? 1 : 0);
- if (id4sq == 5) id4Sav = 2000003 + (isUp ? 1 : 0);
- if (id4sq == 6) id4Sav = 2000005 + (isUp ? 1 : 0);
-
- // Initial values, can be swapped to charge conjugates event by event.
- id3 = id3Sav;
- id4 = id4Sav;
-
- }
-
- // Initialize process.
- virtual void initProc();
-
- // Evaluate d(sigmaHat)/d(tHat).
- virtual double sigmaHat();
-
- // Select flavour, colour and anticolour.
- virtual void setIdColAcol();
-
- private:
-
- // Basic process information
- int id3Sav, id4Sav;
-
-};
-
-//==========================================================================
-
-// A derived class for q q' -> ~q_i ~q_j
-
-class Sigma2qq2squarksquark : public Sigma2Process {
-
-public:
-
- // Constructor.
- Sigma2qq2squarksquark() {}
-
- // Constructor.
- Sigma2qq2squarksquark(int id3In, int id4In, int codeIn) {
-
- // Save ordering indices and process code
- id3Sav = id3In;
- id4Sav = id4In;
- codeSave = codeIn;
- // Initial values (flipped for c.c.)
- id3 = id3Sav;
- id4 = id4Sav;
-
- }
-
- // 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 nameSave;}
- virtual int code() const {return codeSave;}
- virtual string inFlux() const {return "qq";}
- virtual int id3Mass() const {return abs(id3Sav);}
- virtual int id4Mass() const {return abs(id4Sav);}
- virtual bool isSUSY() const {return true;}
-
-private:
-
- // Basic process information
- int id3Sav, id4Sav, codeSave, iGen3, iGen4, nNeut;
- string nameSave;
- bool isUD;
-
- // Storage of mass squares
- double m2Glu;
- vector<double> m2Neut, m2Char;
-
- // Flavor-independent prefactors.
- double sigmaChar, sigmaNeut, sigmaGlu;
- double sigmaCharNeut, sigmaCharGlu, sigmaNeutGlu;
- double openFracPair;
-
- // Point-by-point info
- double tGlu, uGlu;
- vector<double> tNeut, uNeut, tChar, uChar;
- double sumCt, sumCu, sumNt, sumNu, sumGt, sumGu, sumInterference;
-
- //SUSY couplings
- CoupSUSY* coupSUSYPtr;
-};
-
-//==========================================================================
-
-// A derived class for q qbar' -> ~q_i ~q*_j
-
-class Sigma2qqbar2squarkantisquark : public Sigma2Process {
-
-public:
-
- // Constructor.
- Sigma2qqbar2squarkantisquark() {}
-
- // Constructor.
- Sigma2qqbar2squarkantisquark(int id3In, int id4In, int codeIn) {
-
- // Save ordering indices and process code
- // (always store squark first, antisquark second)
- id3Sav = abs(id3In);
- id4Sav = -abs(id4In);
- codeSave = codeIn;
- // Initial values
- id3 = id3Sav;
- id4 = id4Sav;
-
- }
-
- // 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 nameSave;}
- virtual int code() const {return codeSave;}
- virtual string inFlux() const {return "qq";}
- virtual int id3Mass() const {return abs(id3Sav);}
- virtual int id4Mass() const {return abs(id4Sav);}
- virtual bool isSUSY() const {return true;}
-
-private:
-
- // Basic process information
- int id3Sav, id4Sav, codeSave, iGen3, iGen4, nNeut;
- string nameSave;
- bool isUD, isCC;
-
- // Storage of mass squares
- double m2Glu;
- vector<double> m2Neut;
-
- // Flavor-independent prefactors: EW, strong, and interference
- double xW;
- double openFracPair;
- double sigmaEW, sigmaGlu, sigmaEWG;
-
- // Point-by-point info
- double tGlu, uGlu;
- vector<double> tNeut, uNeut;
- complex propZW;
- double sumColS, sumColT, sumColSCC, sumColTCC, sumInterference;
-
- //SUSY couplings
- CoupSUSY* coupSUSYPtr;
-
-};
-
-//==========================================================================
-
-// A derived class for g g -> ~q ~q*
-
-class Sigma2gg2squarkantisquark : public Sigma2Process {
-
-public:
-
- // Constructor.
- Sigma2gg2squarkantisquark() {
- }
-
- // Constructor.
- Sigma2gg2squarkantisquark(int id34In, int codeIn) {
-
- // Save ordering indices and process code
- // (always store squark first, antisquark second)
- id3Sav = abs(id34In);
- id4Sav = -abs(id34In);
- codeSave = codeIn;
- // Initial values
- id3 = id3Sav;
- id4 = id4Sav;
-
- }
-
- // 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 nameSave;}
- virtual int code() const {return codeSave;}
- virtual string inFlux() const {return "gg";}
- virtual int id3Mass() const {return abs(id3Sav);}
- virtual int id4Mass() const {return abs(id4Sav);}
- virtual bool isSUSY() const {return true;}
-
-private:
-
- // Basic process information
- int id3Sav, id4Sav, codeSave;
- string nameSave;
- double sigma, m2Sq, openFracPair;
-
- // Color flow info
- double sumColT, sumColU, sumInterference;
-
- //SUSY couplings
- CoupSUSY* coupSUSYPtr;
-
-};
-
-//==========================================================================
-
-// A derived class for q g -> ~q ~g
-
-class Sigma2qg2squarkgluino : public Sigma2Process {
-
-public:
-
- // Constructor.
- Sigma2qg2squarkgluino() {}
-
- // Constructor.
- Sigma2qg2squarkgluino(int id3In, int codeIn) {
-
- // Save ordering indices and process code
- id3Sav = abs(id3In);
- codeSave = codeIn;
- // Initial values
- id3 = id3Sav;
- id4 = 1000021;
-
- }
-
- // 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 nameSave;}
- virtual int code() const {return codeSave;}
- virtual string inFlux() const {return "qg";}
- virtual int id3Mass() const {return abs(id3Sav);}
- virtual int id4Mass() const {return 1000021;}
- virtual bool isSUSY() const {return true;}
-
-private:
-
- // Basic process information
- int id3Sav, codeSave;
- string nameSave;
- double sigmaA, sigmaB, comFacHat, m2Glu, m2Sq, openFracPair;
-
- //SUSY couplings
- CoupSUSY* coupSUSYPtr;
-
-};
-
-//==========================================================================
-
-// A derived class for g g -> gluino gluino.
-
-class Sigma2gg2gluinogluino : public Sigma2Process {
-
-public:
-
- // Constructor.
- Sigma2gg2gluinogluino() {
- }
-
- // 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 "g g -> gluino gluino";}
- virtual int code() const {return 1201;}
- virtual string inFlux() const {return "gg";}
- virtual int id3Mass() const {return 1000021;}
- virtual int id4Mass() const {return 1000021;}
- virtual bool isSUSY() const {return true;}
-
-private:
-
- // Values stored for process type and colour flow selection.
- double sigTS, sigUS, sigTU, sigSum, sigma, openFracPair;
-
- //SUSY couplings
- CoupSUSY* coupSUSYPtr;
-
-};
-
-//==========================================================================
-
-// A derived class for q qbar -> gluino gluino.
-
-class Sigma2qqbar2gluinogluino : public Sigma2Process {
-
-public:
-
- // Constructor.
- Sigma2qqbar2gluinogluino() {
-
-}
-
- // 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 qbar -> gluino gluino";}
- virtual int code() const {return 1202;}
- virtual string inFlux() const {return "qqbarSame";}
- virtual int id3Mass() const {return 1000021;}
- virtual int id4Mass() const {return 1000021;}
- virtual bool isSUSY() const {return true;}
-
-private:
-
- // Values stored for process type and colour flow selection.
- double openFracPair, s34Avg, sigS, tHG, uHG, tHG2, uHG2;
-
- //SUSY couplings
- CoupSUSY* coupSUSYPtr;
-
-};
-
-//==========================================================================
-
-class Sigma1qq2antisquark : public Sigma1Process {
-public:
-
- // Constructor.
- Sigma1qq2antisquark() {}
-
-
- Sigma1qq2antisquark(int id3In) {
-
- idRes = id3In;
-
- }
-
- // 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 nameSave;}
- virtual int code() const {return codeSave;}
- virtual string inFlux() const {return "qq";}
- virtual bool isSUSY() const {return true;}
- virtual bool isRPV() const {return true;}
- virtual int resonanceA() const {return idRes;}
-
-private:
-
- // Values stored for process type and colour flow selection.
- double mRes, GammaRes, m2Res, GamMRat, sigBW, widthOut;
- int squarkType, codeSave, idRes;
- string nameSave;
-
- //SUSY couplings
- CoupSUSY* coupSUSYPtr;
-
-};
-
-
-//==========================================================================
-} // end namespace Pythia8
-
-
-
-
-
-
-
-#endif // Pythia8_SigmaSUSY_H
-