JHUGen MELA
v2.4.1
Matrix element calculations as used in JHUGen. MELA is an important tool that was used for the Higgs boson discovery and for precise measurements of its structure and interactions. Please see the website https://spin.pha.jhu.edu/ and papers cited there for more details, and kindly cite those papers when using this code.
|
Go to the documentation of this file. 1 #ifndef ZZMatrixElement_ZZMatrixElement_h
2 #define ZZMatrixElement_ZZMatrixElement_h
13 const char* pathtoPDFSet,
15 const char* pathtoHiggsCSandWidth,
36 bool includeHiggsDecay=
false
59 void set_LHAgrid(
const char* path,
int pdfmember=0);
81 void reset_MCFM_EWKParameters(
double ext_Gf,
double ext_aemmz,
double ext_mW,
double ext_mZ,
double ext_xW,
int ext_ewscheme=3);
136 std::vector<TLorentzVector>
Calculate4Momentum(
double Mx,
double M1,
double M2,
double theta,
double theta1,
double theta2,
double Phi1,
double Phi);
void reset_Width(double inmass, int ipart)
double selfDHg4g4coupl[nSupportedHiggses][SIZE_HGG][2]
double selfDHzzpcoupl[SIZE_HVV][2]
double selfDaTQGCcoupl[SIZE_ATQGC][2]
std::vector< MELAParticle * > tmpPartList
double selfDHzzLambda_qsq[nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]
double selfDGggcoupl[SIZE_GGG][2]
void get_XPropagator(TVar::ResonancePropagatorScheme scheme, float &prop)
double selfDZpffcoupl[SIZE_Vpff][2]
void set_mHiggs_wHiggs(double mh_, double gah_, int index)
double selfDHb4b4coupl[nSupportedHiggses][SIZE_HQQ][2]
double selfDZqqcoupl[SIZE_ZQQ][2]
std::vector< MELATopCandidate_t * > * get_TopCandidateCollection()
void reset_Mass(double inmass, int ipart)
TVar::LeptonInterference processLeptonInterference
double selfDGqqcoupl[SIZE_GQQ][2]
double selfDHt4t4coupl[nSupportedHiggses][SIZE_HQQ][2]
void computeProdXS_JH(float &mevalue)
double selfDHbbcoupl[nSupportedHiggses][SIZE_HQQ][2]
double wHiggs[nSupportedHiggses]
void set_Verbosity(TVar::VerbosityLevel verbosity_)
void set_SpinZeroContact(double selfDHzzpcoupl[SIZE_HVV][2], double selfDHzpzpcoupl[SIZE_HVV][2], double selfDHwwpcoupl[SIZE_HVV][2], double selfDHwpwpcoupl[SIZE_HVV][2])
void set_SpinOneCouplings(double selfDZqqcoupl[SIZE_ZQQ][2], double selfDZvvcoupl[SIZE_ZVV][2])
TVar::Process processModel
double get_HiggsWidthAtPoleMass(double mass)
double selfDHzpzpcoupl[SIZE_HVV][2]
std::vector< TLorentzVector > Calculate4Momentum(double Mx, double M1, double M2, double theta, double theta1, double theta2, double Phi1, double Phi)
double mHiggs[nSupportedHiggses]
void set_SpinZeroCouplings(double selfDHggcoupl[nSupportedHiggses][SIZE_HGG][2], double selfDHg4g4coupl[nSupportedHiggses][SIZE_HGG][2], double selfDHqqcoupl[nSupportedHiggses][SIZE_HQQ][2], double selfDHbbcoupl[nSupportedHiggses][SIZE_HQQ][2], double selfDHttcoupl[nSupportedHiggses][SIZE_HQQ][2], double selfDHb4b4coupl[nSupportedHiggses][SIZE_HQQ][2], double selfDHt4t4coupl[nSupportedHiggses][SIZE_HQQ][2], double selfDHzzcoupl[nSupportedHiggses][SIZE_HVV][2], double selfDHwwcoupl[nSupportedHiggses][SIZE_HVV][2], double selfDHzzLambda_qsq[nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ], double selfDHwwLambda_qsq[nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ], int selfDHzzCLambda_qsq[nSupportedHiggses][SIZE_HVV_CQSQ], int selfDHwwCLambda_qsq[nSupportedHiggses][SIZE_HVV_CQSQ], bool diffHWW=false)
int selfDHwwCLambda_qsq[nSupportedHiggses][SIZE_HVV_CQSQ]
TVar::MatrixElement processME
std::vector< MELACandidate * > tmpCandList
void append_TopCandidate(SimpleParticleCollection_t *TopDaughters)
void set_CandidateDecayMode(TVar::CandidateDecayMode mode)
void set_TempCandidate(SimpleParticleCollection_t *pDaughters, SimpleParticleCollection_t *pAssociated=0, SimpleParticleCollection_t *pMothers=0, bool isGen=false)
double get_PrimaryMass(int ipart)
void set_AZffCouplings(double selfDAZffcoupl[SIZE_AZff][2])
void set_SpinTwoCouplings(double selfDGqqcoupl[SIZE_GQQ][2], double selfDGggcoupl[SIZE_GGG][2], double selfDGvvcoupl[SIZE_GVV][2])
int get_CurrentCandidateIndex()
void set_Process(TVar::Process process_, TVar::MatrixElement me_, TVar::Production production_)
MELACandidate * get_CurrentCandidate()
double selfDHzzcoupl[nSupportedHiggses][SIZE_HVV][2]
void set_mHiggs(double mh_, int index)
void set_VprimeContactCouplings(double selfDZpffcoupl[SIZE_Vpff][2], double selfDWpffcoupl[SIZE_Vpff][2], double M_Zprime, double Ga_Zprime, double M_Wprime, double Ga_Wprime)
double selfDHqqcoupl[nSupportedHiggses][SIZE_HQQ][2]
AZffCouplings * selfD_AZffCouplings
TVar::VerbosityLevel processVerbosity
VprimeCouplings * selfD_VprimeCouplings
void set_InputEvent(SimpleParticleCollection_t *pDaughters, SimpleParticleCollection_t *pAssociated=0, SimpleParticleCollection_t *pMothers=0, bool isGen=false)
void set_CurrentCandidateFromIndex(unsigned int icand)
int selfDHzzCLambda_qsq[nSupportedHiggses][SIZE_HVV_CQSQ]
void computeProdXS_VH(float &mevalue, bool includeHiggsDecay=false)
void computeProdXS_JJH(float &mevalue)
ZZMatrixElement(const char *pathtoPDFSet, int PDFMember, const char *pathtoHiggsCSandWidth, double ebeam, TVar::VerbosityLevel verbosity)
double selfDHggcoupl[nSupportedHiggses][SIZE_HGG][2]
void computeProdXS_VVHVV(float &mevalue)
void set_RenFacScaleMode(TVar::EventScaleScheme renormalizationSch, TVar::EventScaleScheme factorizationSch, double ren_sf, double fac_sf)
aTQGCCouplings * selfD_aTQGCCouplings
double get_PrimaryHiggsMass()
double selfDHwwpcoupl[SIZE_HVV][2]
void computeXS(float &mevalue)
void set_LHAgrid(const char *path, int pdfmember=0)
double selfDHttcoupl[nSupportedHiggses][SIZE_HQQ][2]
ResonancePropagatorScheme
SpinZeroCouplings * selfD_SpinZeroCouplings
double selfDWpffcoupl[SIZE_Vpff][2]
void set_CurrentCandidate(MELACandidate *cand)
double selfDGvpvpcoupl[SIZE_GVV][2]
double selfDAZffcoupl[SIZE_AZff][2]
void set_LeptonInterference(TVar::LeptonInterference myLepInterf)
void set_aTQGCCouplings(double selfDaTQGCcoupl[SIZE_ATQGC][2])
SpinOneCouplings * selfD_SpinOneCouplings
std::vector< SimpleParticle_t > SimpleParticleCollection_t
double selfDHwpwpcoupl[SIZE_HVV][2]
void reset_MCFM_EWKParameters(double ext_Gf, double ext_aemmz, double ext_mW, double ext_mZ, double ext_xW, int ext_ewscheme=3)
void computeProdXS_ttH(float &mevalue, int topProcess, int topDecay=0)
void set_SpinTwoContact(double selfDGvvpcoupl[SIZE_GVV][2], double selfDGvpvpcoupl[SIZE_GVV][2])
double get_PrimaryWidth(int ipart)
double selfDHwwcoupl[nSupportedHiggses][SIZE_HVV][2]
TVar::Production processProduction
double selfDGvvpcoupl[SIZE_GVV][2]
double selfDHwwLambda_qsq[nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]
double selfDZvvcoupl[SIZE_ZVV][2]
SpinTwoCouplings * selfD_SpinTwoCouplings
double selfDGvvcoupl[SIZE_GVV][2]
void set_PrimaryHiggsMass(double mh)
void set_wHiggs(double gah_, int index)