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.
Functions | Variables
Coupling arrays

Self-define arrays for MELA. More...

Functions

static void Mela::cleanLinkedFiles ()
 
void Mela::printLogo () const
 
void Mela::setSpinZeroCouplings ()
 
void Mela::setSpinOneCouplings ()
 
void Mela::setSpinTwoCouplings ()
 
void Mela::setATQGCCouplings ()
 
void Mela::setAZffCouplings ()
 
bool Mela::configureAnalyticalPDFs ()
 
void Mela::reset_SelfDCouplings ()
 
void Mela::reset_PAux ()
 
void Mela::reset_CandRef ()
 
void Mela::constructDggr (float bkg_VAMCFM_noscale, float ggzz_VAMCFM_noscale, float ggHZZ_prob_pure_noscale, float ggHZZ_prob_int_noscale, float widthScale, float &myDggr)
 
void Mela::getPConstantHandles ()
 
void Mela::deletePConstantHandles ()
 
MelaPConstant * Mela::getPConstantHandle (TVar::MatrixElement me_, TVar::Production prod_, TVar::Process proc_, TString relpath, TString spname, const bool useSqrts=false)
 
void Mela::deletePConstantHandle (MelaPConstant *&handle)
 
void Mela::computeConstant (float &prob)
 
void Mela::setConstant ()
 
float Mela::getConstant_JHUGenUndecayed ()
 
float Mela::getConstant_4l ()
 
float Mela::getConstant_2l2q ()
 
float Mela::getConstant_4q ()
 
float Mela::getConstant_FourFermionDecay (const int &decid)
 

Variables

double Mela::selfDHggcoupl [nSupportedHiggses][SIZE_HGG][2]
 
double Mela::selfDHg4g4coupl [nSupportedHiggses][SIZE_HGG][2]
 
double Mela::selfDHqqcoupl [nSupportedHiggses][SIZE_HQQ][2]
 
double Mela::selfDHbbcoupl [nSupportedHiggses][SIZE_HQQ][2]
 
double Mela::selfDHttcoupl [nSupportedHiggses][SIZE_HQQ][2]
 
double Mela::selfDHb4b4coupl [nSupportedHiggses][SIZE_HQQ][2]
 
double Mela::selfDHt4t4coupl [nSupportedHiggses][SIZE_HQQ][2]
 
double Mela::selfDHzzcoupl [nSupportedHiggses][SIZE_HVV][2]
 
double Mela::selfDHwwcoupl [nSupportedHiggses][SIZE_HVV][2]
 
double Mela::selfDHzzLambda_qsq [nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]
 
double Mela::selfDHwwLambda_qsq [nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]
 
int Mela::selfDHzzCLambda_qsq [nSupportedHiggses][SIZE_HVV_CQSQ]
 
int Mela::selfDHwwCLambda_qsq [nSupportedHiggses][SIZE_HVV_CQSQ]
 
bool Mela::differentiate_HWW_HZZ
 
double Mela::selfDHzzpcoupl [SIZE_HVV][2]
 
double Mela::selfDHzpzpcoupl [SIZE_HVV][2]
 
double Mela::selfDZpffcoupl [SIZE_Vpff][2]
 
double Mela::selfDHwwpcoupl [SIZE_HVV][2]
 
double Mela::selfDHwpwpcoupl [SIZE_HVV][2]
 
double Mela::selfDWpffcoupl [SIZE_Vpff][2]
 
double Mela::selfDM_Zprime
 
double Mela::selfDGa_Zprime
 
double Mela::selfDM_Wprime
 
double Mela::selfDGa_Wprime
 
double Mela::selfDZqqcoupl [SIZE_ZQQ][2]
 
double Mela::selfDZvvcoupl [SIZE_ZVV][2]
 
double Mela::selfDGqqcoupl [SIZE_GQQ][2]
 
double Mela::selfDGggcoupl [SIZE_GGG][2]
 
double Mela::selfDGvvcoupl [SIZE_GVV][2]
 
double Mela::selfDGvvpcoupl [SIZE_GVV][2]
 
double Mela::selfDGvpvpcoupl [SIZE_GVV][2]
 
double Mela::selfDaTQGCcoupl [SIZE_ATQGC][2]
 
double Mela::selfDAZffcoupl [SIZE_AZff][2]
 
double Mela::LHCsqrts
 
TVar::Process Mela::myModel_
 
TVar::MatrixElement Mela::myME_
 
TVar::Production Mela::myProduction_
 
TVar::LeptonInterference Mela::myLepInterf_
 
TVar::VerbosityLevel Mela::myVerbosity_
 
ZZMatrixElement * Mela::ZZME
 
SuperDijetMela * Mela::superDijet
 
float Mela::auxiliaryProb
 
MELACandidate * Mela::melaCand
 
MelaPConstant * Mela::pAvgSmooth_JHUGen_JQCD_HSMHiggs [TVar::nFermionMassRemovalSchemes-1]
 
MelaPConstant * Mela::pAvgSmooth_JHUGen_JJQCD_HSMHiggs [TVar::nFermionMassRemovalSchemes-1]
 
MelaPConstant * Mela::pAvgSmooth_JHUGen_JJVBF_HSMHiggs [TVar::nFermionMassRemovalSchemes-1]
 
MelaPConstant * Mela::pAvgSmooth_JHUGen_Had_ZH_HSMHiggs [TVar::nFermionMassRemovalSchemes-1]
 
MelaPConstant * Mela::pAvgSmooth_JHUGen_Had_WH_HSMHiggs [TVar::nFermionMassRemovalSchemes-1]
 
MelaPConstant * Mela::pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4mu
 
MelaPConstant * Mela::pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4e
 
MelaPConstant * Mela::pAvgSmooth_JHUGen_ZZGG_HSMHiggs_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_ZZGG_bkgZZ_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJVBF_bkgZZ_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_WH_bkgZZ_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJQCD_bkgZZ_4e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e
 
MelaPConstant * Mela::pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q
 

Detailed Description

Self-define arrays for MELA.

Attention
The first dimension (of size [TMCFM::nSupportedHiggses=2]) supports a second resonance present in MCFM

Function Documentation

◆ cleanLinkedFiles()

void Mela::cleanLinkedFiles ( )
static

Definition at line 162 of file Mela.cc.

162  {
163  remove("ffwarn.dat");
164  remove("br.sm1");
165  remove("br.sm2");
166  remove("input.DAT");
167  remove("process.DAT");
168  remove("Pdfdata/cteq6l1.tbl");
169  remove("Pdfdata/cteq6l.tbl");
170  remove("Pdfdata/NNPDF30_lo_as_0130.LHgrid");
171  rmdir("Pdfdata");
172 }

◆ computeConstant()

void Mela::computeConstant ( float &  prob)
protected

Definition at line 2169 of file Mela.cc.

2169  {
2170  float pConst=1.;
2171  setConstant();
2172  getConstant(pConst);
2173  prob *= pConst;
2174 }

◆ configureAnalyticalPDFs()

bool Mela::configureAnalyticalPDFs ( )
protected

Definition at line 1969 of file Mela.cc.

1969  {
1970  //
1971  // Configure the analytical calculations
1972  //
1973  bool noPass=false;
1974  pdf=0;
1975 
1976  if (myModel_==TVar::bkgZZ) pdf = qqZZmodel;
1977  else if (myProduction_ == TVar::JJQCD || myProduction_ == TVar::JJVBF);
1978  else if (myProduction_ == TVar::Lep_ZH || myProduction_ == TVar::Lep_WH || myProduction_ == TVar::Had_ZH || myProduction_ == TVar::Had_WH || myProduction_ == TVar::GammaH); // To be implemented
1979  else if (
1980  myModel_ == TVar::HSMHiggs
1981  || myModel_ == TVar::H0minus || myModel_ == TVar::D_g1g4 || myModel_ == TVar::D_g1g4_pi_2
1982  || myModel_ == TVar::H0hplus || myModel_ == TVar::D_g1g2 || myModel_ == TVar::D_g1g2_pi_2
1983  || myModel_ == TVar::H0_g1prime2 || myModel_ == TVar::D_g1g1prime2
1984  || myModel_ == TVar::SelfDefine_spin0
1985  ){
1986  pdf = (RooAbsPdf*)ggSpin0Model->getPDF();
1987  ggSpin0Model->makeParamsConst(false);
1988  ggSpin0Model->resetHypotheses();
1989 
1990  // Add the hypotheses with best-guess coefficients
1991  // ZZ/WW
1992  if (
1993  myModel_ == TVar::HSMHiggs
1994  || myModel_ == TVar::D_g1g1prime2
1995  || myModel_ == TVar::D_g1g2 || myModel_ == TVar::D_g1g2_pi_2
1996  || myModel_ == TVar::D_g1g4 || myModel_ == TVar::D_g1g4_pi_2
1997  || myModel_ == TVar::D_zzzg_g1prime2 || myModel_ == TVar::D_zzzg_g1prime2_pi_2 || myModel_ == TVar::D_zzzg || myModel_ == TVar::D_zzzg_PS
1998  || myModel_ == TVar::D_zzgg || myModel_ == TVar::D_zzgg_PS
1999  ) ggSpin0Model->addHypothesis(0, 0);
2000  if (myModel_ == TVar::H0_g1prime2 || myModel_ == TVar::D_g1g1prime2) ggSpin0Model->addHypothesis(0, 2);
2001  if (myModel_ == TVar::H0hplus || myModel_ == TVar::D_g1g2 || myModel_ == TVar::D_g1g2_pi_2) ggSpin0Model->addHypothesis(1, 0, (myModel_ == TVar::D_g1g2_pi_2 ? TMath::Pi() : 0.));
2002  if (myModel_ == TVar::H0minus || myModel_ == TVar::D_g1g4 || myModel_ == TVar::D_g1g4_pi_2) ggSpin0Model->addHypothesis(3, 0, (myModel_ == TVar::D_g1g4_pi_2 ? TMath::Pi() : 0.));
2003  // ZG/ZGs
2004  if (myModel_ == TVar::H0_Zgsg1prime2 || myModel_ == TVar::D_zzzg_g1prime2 || myModel_ == TVar::D_zzzg_g1prime2_pi_2) ggSpin0Model->addHypothesis(4, 2, (myModel_ == TVar::D_zzzg_g1prime2_pi_2 ? TMath::Pi() : 0.));
2005  if (myModel_ == TVar::H0_Zgs || myModel_ == TVar::D_zzzg) ggSpin0Model->addHypothesis(5, 0);
2006  if (myModel_ == TVar::H0_Zgs_PS || myModel_ == TVar::D_zzzg_PS) ggSpin0Model->addHypothesis(7, 0);
2007  // GG/GGs/GsGs
2008  if (myModel_ == TVar::H0_gsgs || myModel_ == TVar::D_zzgg) ggSpin0Model->addHypothesis(8, 0);
2009  if (myModel_ == TVar::H0_gsgs_PS || myModel_ == TVar::D_zzgg_PS) ggSpin0Model->addHypothesis(10, 0);
2010  // Self-defined spin-0
2011  if (myModel_ == TVar::SelfDefine_spin0){
2012  for (int im=0; im<2; im++){
2013  ((RooRealVar*)ggSpin0Model->couplings.g1List[0][im])->setVal(selfDHzzcoupl[0][0][im]);
2014  ((RooRealVar*)ggSpin0Model->couplings.g2List[0][im])->setVal(selfDHzzcoupl[0][1][im]);
2015  ((RooRealVar*)ggSpin0Model->couplings.g3List[0][im])->setVal(selfDHzzcoupl[0][2][im]);
2016  ((RooRealVar*)ggSpin0Model->couplings.g4List[0][im])->setVal(selfDHzzcoupl[0][3][im]);
2017 
2018  ((RooRealVar*)ggSpin0Model->couplings.gzgs2List[0][im])->setVal(selfDHzzcoupl[0][4][im]);
2019  ((RooRealVar*)ggSpin0Model->couplings.gzgs3List[0][im])->setVal(selfDHzzcoupl[0][5][im]);
2020  ((RooRealVar*)ggSpin0Model->couplings.gzgs4List[0][im])->setVal(selfDHzzcoupl[0][6][im]);
2021  ((RooRealVar*)ggSpin0Model->couplings.ggsgs2List[0][im])->setVal(selfDHzzcoupl[0][7][im]);
2022  ((RooRealVar*)ggSpin0Model->couplings.ggsgs3List[0][im])->setVal(selfDHzzcoupl[0][8][im]);
2023  ((RooRealVar*)ggSpin0Model->couplings.ggsgs4List[0][im])->setVal(selfDHzzcoupl[0][9][im]);
2024 
2025  ((RooRealVar*)ggSpin0Model->couplings.g1List[1][im])->setVal(selfDHzzcoupl[0][10][im]);
2026  ((RooRealVar*)ggSpin0Model->couplings.g1List[2][im])->setVal(selfDHzzcoupl[0][11][im]);
2027  ((RooRealVar*)ggSpin0Model->couplings.g1List[3][im])->setVal(selfDHzzcoupl[0][12][im]);
2028  ((RooRealVar*)ggSpin0Model->couplings.g1List[4][im])->setVal(selfDHzzcoupl[0][13][im]);
2029  ((RooRealVar*)ggSpin0Model->couplings.g1List[5][im])->setVal(selfDHzzcoupl[0][14][im]);
2030 
2031  ((RooRealVar*)ggSpin0Model->couplings.g2List[1][im])->setVal(selfDHzzcoupl[0][15][im]);
2032  ((RooRealVar*)ggSpin0Model->couplings.g2List[2][im])->setVal(selfDHzzcoupl[0][16][im]);
2033  ((RooRealVar*)ggSpin0Model->couplings.g2List[3][im])->setVal(selfDHzzcoupl[0][17][im]);
2034  ((RooRealVar*)ggSpin0Model->couplings.g2List[4][im])->setVal(selfDHzzcoupl[0][18][im]);
2035  ((RooRealVar*)ggSpin0Model->couplings.g2List[5][im])->setVal(selfDHzzcoupl[0][19][im]);
2036 
2037  ((RooRealVar*)ggSpin0Model->couplings.g3List[1][im])->setVal(selfDHzzcoupl[0][20][im]);
2038  ((RooRealVar*)ggSpin0Model->couplings.g3List[2][im])->setVal(selfDHzzcoupl[0][21][im]);
2039  ((RooRealVar*)ggSpin0Model->couplings.g3List[3][im])->setVal(selfDHzzcoupl[0][22][im]);
2040  ((RooRealVar*)ggSpin0Model->couplings.g3List[4][im])->setVal(selfDHzzcoupl[0][23][im]);
2041  ((RooRealVar*)ggSpin0Model->couplings.g3List[5][im])->setVal(selfDHzzcoupl[0][24][im]);
2042 
2043  ((RooRealVar*)ggSpin0Model->couplings.g4List[1][im])->setVal(selfDHzzcoupl[0][25][im]);
2044  ((RooRealVar*)ggSpin0Model->couplings.g4List[2][im])->setVal(selfDHzzcoupl[0][26][im]);
2045  ((RooRealVar*)ggSpin0Model->couplings.g4List[3][im])->setVal(selfDHzzcoupl[0][27][im]);
2046  ((RooRealVar*)ggSpin0Model->couplings.g4List[4][im])->setVal(selfDHzzcoupl[0][28][im]);
2047  ((RooRealVar*)ggSpin0Model->couplings.g4List[5][im])->setVal(selfDHzzcoupl[0][29][im]);
2048 
2049  ((RooRealVar*)ggSpin0Model->couplings.gzgs1List[0][im])->setVal(selfDHzzcoupl[0][30][im]); // Zgs1_prime2
2050 
2051  ((RooRealVar*)ggSpin0Model->couplings.g1List[6][im])->setVal(selfDHzzcoupl[0][31][im]);
2052  ((RooRealVar*)ggSpin0Model->couplings.g1List[7][im])->setVal(selfDHzzcoupl[0][32][im]);
2053  ((RooRealVar*)ggSpin0Model->couplings.g2List[6][im])->setVal(selfDHzzcoupl[0][33][im]);
2054  ((RooRealVar*)ggSpin0Model->couplings.g2List[7][im])->setVal(selfDHzzcoupl[0][34][im]);
2055  ((RooRealVar*)ggSpin0Model->couplings.g3List[6][im])->setVal(selfDHzzcoupl[0][35][im]);
2056  ((RooRealVar*)ggSpin0Model->couplings.g3List[7][im])->setVal(selfDHzzcoupl[0][36][im]);
2057  ((RooRealVar*)ggSpin0Model->couplings.g4List[6][im])->setVal(selfDHzzcoupl[0][37][im]);
2058  ((RooRealVar*)ggSpin0Model->couplings.g4List[7][im])->setVal(selfDHzzcoupl[0][38][im]);
2059  }
2060  for (int qoqtqz=0; qoqtqz<SIZE_HVV_CQSQ; qoqtqz++){ // 0==q1, 1==q2, 2==q12
2061  ((RooRealVar*)ggSpin0Model->couplings.Lambda_z1qsq[qoqtqz])->setVal(selfDHzzLambda_qsq[0][0][qoqtqz]);
2062  ((RooRealVar*)ggSpin0Model->couplings.Lambda_z2qsq[qoqtqz])->setVal(selfDHzzLambda_qsq[0][1][qoqtqz]);
2063  ((RooRealVar*)ggSpin0Model->couplings.Lambda_z3qsq[qoqtqz])->setVal(selfDHzzLambda_qsq[0][2][qoqtqz]);
2064  ((RooRealVar*)ggSpin0Model->couplings.Lambda_z4qsq[qoqtqz])->setVal(selfDHzzLambda_qsq[0][3][qoqtqz]);
2065  ((RooRealVar*)ggSpin0Model->couplings.cLambda_qsq[qoqtqz])->setVal(selfDHzzCLambda_qsq[0][qoqtqz]);
2066  }
2067  }
2068  ggSpin0Model->makeParamsConst(true);
2069  }
2070  else if (!spin1Model->configure(myModel_)){
2071  pdf = spin1Model->PDF;
2072  // Self-defined spin-1
2073  if (myModel_ == TVar::SelfDefine_spin1){
2074  for (int i=0; i<SIZE_ZVV; i++){ if (selfDZvvcoupl[i][1]!=0){ if (myVerbosity_>=TVar::ERROR) MELAerr << "Mela::configureAnalyticalPDFs: MELA does not support complex couplings for spin-1 at the moment! " << endl; noPass=true; break; } }
2075  if (!noPass){
2076  spin1Model->g1Val->setVal(selfDZvvcoupl[0][0]);
2077  spin1Model->g2Val->setVal(selfDZvvcoupl[1][0]);
2078  }
2079  }
2080  }
2081  else if (
2082  myModel_ == TVar::H2_g1
2083  || myModel_ == TVar::H2_g1g5
2084  || myModel_ == TVar::H2_g2
2085  || myModel_ == TVar::H2_g3
2086  || myModel_ == TVar::H2_g4
2087  || myModel_ == TVar::H2_g5
2088  || myModel_ == TVar::H2_g6
2089  || myModel_ == TVar::H2_g7
2090  || myModel_ == TVar::H2_g8
2091  || myModel_ == TVar::H2_g9
2092  || myModel_ == TVar::H2_g10
2093  || myModel_ == TVar::SelfDefine_spin2
2094  ){
2095  pdf = (RooAbsPdf*)spin2Model->getPDF();
2096  spin2Model->makeParamsConst(false);
2097  spin2Model->resetHypotheses();
2098  // Add the hypotheses with best-guess coefficients
2099  // ZZ/WW
2100  if (
2101  myModel_ == TVar::H2_g1
2102  || myModel_ == TVar::H2_g1g5
2103  ) spin2Model->addHypothesis(0, 1.);
2104  if (
2105  myModel_ == TVar::H2_g1g5
2106  || myModel_ == TVar::H2_g5
2107  ) spin2Model->addHypothesis(4, 1.);
2108  if (myModel_ == TVar::H2_g2) spin2Model->addHypothesis(1, 1.);
2109  if (myModel_ == TVar::H2_g3) spin2Model->addHypothesis(2, 1.);
2110  if (myModel_ == TVar::H2_g4) spin2Model->addHypothesis(3, 1.);
2111  if (myModel_ == TVar::H2_g5) spin2Model->addHypothesis(4, 1.);
2112  if (myModel_ == TVar::H2_g6) spin2Model->addHypothesis(5, 1.);
2113  if (myModel_ == TVar::H2_g7) spin2Model->addHypothesis(6, 1.);
2114  if (myModel_ == TVar::H2_g8) spin2Model->addHypothesis(7, 1.);
2115  if (myModel_ == TVar::H2_g9) spin2Model->addHypothesis(8, 1.);
2116  if (myModel_ == TVar::H2_g10) spin2Model->addHypothesis(10, 1.);
2117  // Self-defined spin-2
2118  if (myModel_ == TVar::SelfDefine_spin2){
2119  for (int i=0; i<SIZE_GVV; i++){ if (selfDGvvcoupl[i][1]!=0.){ if (myVerbosity_>=TVar::ERROR) MELAerr << "Mela::configureAnalyticalPDFs: MELA does not support complex couplings for spin-2 at the moment! " << endl; noPass=true; break; } }
2120  if (!noPass){
2121  for (int ig=0; ig<SIZE_GVV; ig++){
2122  for (int im=0; im<2; im++) ((RooRealVar*)spin2Model->couplings.bList[ig][im])->setVal(selfDGvvcoupl[ig][im]);
2123  }
2124  }
2125  }
2126  if (!noPass){
2127  if (myProduction_ == TVar::ZZQQB){
2128  spin2Model->setTensorPolarization(1, 1.);
2129  spin2Model->setTensorPolarization(2, 0.);
2130  }
2131  else{
2132  if (myModel_ == TVar::SelfDefine_spin2){
2133  double c1 = 2*selfDGggcoupl[0][0] + 2.*selfDGggcoupl[1][0];
2134  double c2 = -0.5*selfDGggcoupl[0][0] + selfDGggcoupl[2][0] + 2.*selfDGggcoupl[3][0];
2135  double c5 = 0./*4*selfDGggcoupl[7][0]*/;
2136  Double_t fppReal = 1./sqrt(6.) * (c1/4.*2. + 2.*c2);
2137  Double_t fppImag = 1./sqrt(6.) * c5;
2138  Double_t fmmReal = 1./sqrt(6.) * (c1/4.*2. + 2.*c2);
2139  Double_t fmmImag = 1./sqrt(6.)* c5;
2140  Double_t fmpReal = 1./4.*c1*2.;
2141  Double_t fmpImag = 0;
2142  Double_t fpp = fppImag*fppImag + fppReal*fppReal;
2143  Double_t fmm = fmmImag*fmmImag + fmmReal*fmmReal;
2144  Double_t fmp = fmpImag*fmpImag + fmpReal*fmpReal;
2145  spin2Model->setTensorPolarization(1, 0.); // This is wrong in the strict sense of what "SelfDefine_spin2" is.
2146  spin2Model->setTensorPolarization(2, 2.*fmp/(fmm+fpp+2.*fmp));
2147  }
2148  else{
2149  spin2Model->setTensorPolarization(1, 0.);
2150  spin2Model->setTensorPolarization(2, 1.);
2151  }
2152  }
2153  spin2Model->makeParamsConst(true);
2154  }
2155  }
2156  else{
2157  MELAout << "Mela::configureAnalyticalPDFs -> ERROR TVar::Process not applicable!!! ME: " << myME_ << ", model: " << myModel_ << endl;
2158  noPass=true;
2159  }
2160 
2161  if (pdf==0) noPass=true;
2162  return (!noPass);
2163 }

◆ constructDggr()

void Mela::constructDggr ( float  bkg_VAMCFM_noscale,
float  ggzz_VAMCFM_noscale,
float  ggHZZ_prob_pure_noscale,
float  ggHZZ_prob_int_noscale,
float  widthScale,
float &  myDggr 
)
protected

Definition at line 1920 of file Mela.cc.

1927  {
1928  float total_sig_ME = (widthScale * ggHZZ_prob_pure_noscale + sqrt(widthScale) * ggHZZ_prob_int_noscale + ggzz_VAMCFM_noscale);
1929  float total_bkg_ME = bkg_VAMCFM_noscale;
1930  float kd_denominator = (total_sig_ME+total_bkg_ME);
1931  if (kd_denominator>0.) myDggr = total_sig_ME/(total_sig_ME+total_bkg_ME);
1932  else myDggr=-99.;
1933 }

◆ deletePConstantHandle()

void Mela::deletePConstantHandle ( MelaPConstant *&  handle)
protected

Definition at line 2785 of file Mela.cc.

2785  {
2786  if (myVerbosity_>=TVar::DEBUG) MELAout << "Mela::deletePConstantHandle: Deleting PConstant handle " << handle->GetSplineName() << " at " << handle->GetFileName() << endl;
2787  delete handle; handle=0;
2788  if (myVerbosity_>=TVar::DEBUG) MELAout << "End Mela::deletePConstantHandle." << endl;
2789 }

◆ deletePConstantHandles()

void Mela::deletePConstantHandles ( )
protected

Definition at line 2729 of file Mela.cc.

2729  {
2730  for (int isch=(unsigned int)(TVar::nFermionMassRemovalSchemes-2); isch>=0; isch--){
2731  if (isch==0 || pAvgSmooth_JHUGen_JJQCD_HSMHiggs[isch]!=pAvgSmooth_JHUGen_JJQCD_HSMHiggs[0]) deletePConstantHandle(pAvgSmooth_JHUGen_JJQCD_HSMHiggs[isch]);
2732  if (isch==0 || pAvgSmooth_JHUGen_JQCD_HSMHiggs[isch]!=pAvgSmooth_JHUGen_JQCD_HSMHiggs[0]) deletePConstantHandle(pAvgSmooth_JHUGen_JQCD_HSMHiggs[isch]);
2733  if (isch==0 || pAvgSmooth_JHUGen_JJVBF_HSMHiggs[isch]!=pAvgSmooth_JHUGen_JJVBF_HSMHiggs[0]) deletePConstantHandle(pAvgSmooth_JHUGen_JJVBF_HSMHiggs[isch]);
2734  if (isch==0 || pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[isch]!=pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[0]) deletePConstantHandle(pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[isch]);
2735  if (isch==0 || pAvgSmooth_JHUGen_Had_WH_HSMHiggs[isch]!=pAvgSmooth_JHUGen_Had_WH_HSMHiggs[0]) deletePConstantHandle(pAvgSmooth_JHUGen_Had_WH_HSMHiggs[isch]);
2736  }
2737  //
2738  deletePConstantHandle(pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q);
2739  //
2740  deletePConstantHandle(pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4mu);
2741  deletePConstantHandle(pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4e);
2742  deletePConstantHandle(pAvgSmooth_JHUGen_ZZGG_HSMHiggs_2mu2e);
2743  //
2744  deletePConstantHandle(pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu);
2745  deletePConstantHandle(pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e);
2746  deletePConstantHandle(pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e);
2747  //
2748  deletePConstantHandle(pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4mu);
2749  deletePConstantHandle(pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4e);
2750  deletePConstantHandle(pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_2mu2e);
2751  //
2752  deletePConstantHandle(pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4mu);
2753  deletePConstantHandle(pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4e);
2754  deletePConstantHandle(pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_2mu2e);
2755  //
2756  deletePConstantHandle(pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4mu);
2757  deletePConstantHandle(pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4e);
2758  deletePConstantHandle(pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_2mu2e);
2759  //
2760  deletePConstantHandle(pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu);
2761  deletePConstantHandle(pAvgSmooth_MCFM_ZZGG_bkgZZ_4e);
2762  deletePConstantHandle(pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e);
2763  //
2764  deletePConstantHandle(pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu);
2765  deletePConstantHandle(pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e);
2766  deletePConstantHandle(pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e);
2767  //
2768  deletePConstantHandle(pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu);
2769  deletePConstantHandle(pAvgSmooth_MCFM_JJVBF_bkgZZ_4e);
2770  deletePConstantHandle(pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e);
2771  //
2772  deletePConstantHandle(pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu);
2773  deletePConstantHandle(pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e);
2774  deletePConstantHandle(pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e);
2775  //
2776  deletePConstantHandle(pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu);
2777  deletePConstantHandle(pAvgSmooth_MCFM_Had_WH_bkgZZ_4e);
2778  deletePConstantHandle(pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e);
2779  //
2780  deletePConstantHandle(pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu);
2781  deletePConstantHandle(pAvgSmooth_MCFM_JJQCD_bkgZZ_4e);
2782  deletePConstantHandle(pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e);
2783  //
2784 }

◆ getConstant_2l2q()

float Mela::getConstant_2l2q ( )
protected

Definition at line 2297 of file Mela.cc.

2297  {
2298  float constant = 1;
2299  if (!melaCand) return constant;
2300  int decid = 1;
2301  if (PDGHelpers::isALepton(melaCand->getSortedV(0)->getDaughter(0)->id)) decid *= melaCand->getSortedV(0)->getDaughter(0)->id*melaCand->getSortedV(0)->getDaughter(1)->id;
2302  if (PDGHelpers::isALepton(melaCand->getSortedV(1)->getDaughter(0)->id)) decid *= melaCand->getSortedV(1)->getDaughter(0)->id*melaCand->getSortedV(1)->getDaughter(1)->id;
2303  decid = abs(decid);
2304  return getConstant_FourFermionDecay(decid);
2305 }

◆ getConstant_4l()

float Mela::getConstant_4l ( )
protected

Definition at line 2286 of file Mela.cc.

2286  {
2287  float constant = 1;
2288  if (!melaCand) return constant;
2289  int decid = abs(
2290  melaCand->getSortedV(0)->getDaughter(0)->id*
2291  melaCand->getSortedV(0)->getDaughter(1)->id*
2292  melaCand->getSortedV(1)->getDaughter(0)->id*
2293  melaCand->getSortedV(1)->getDaughter(1)->id
2294  );
2295  return getConstant_FourFermionDecay(decid);
2296 }

◆ getConstant_4q()

float Mela::getConstant_4q ( )
protected

Definition at line 2306 of file Mela.cc.

2306  {
2307  float constant = 1;
2308  if (!melaCand) return constant;
2309  const int decid = 121;
2310  return getConstant_FourFermionDecay(decid);
2311 }

◆ getConstant_FourFermionDecay()

float Mela::getConstant_FourFermionDecay ( const int &  decid)
protected

Definition at line 2313 of file Mela.cc.

2313  {
2314  float constant = 1;
2315 
2316  const bool is4mu = (decid==28561);
2317  const bool is4e = (decid==14641 || decid==50625); // Use 4e for 4tau as well (I don't know why you would do this, but anyway)
2318  const bool is2mu2e = (decid==20449 || decid==27225 || decid==38025 || decid==169 || decid==121); // Use 2e2mu for 2e2tau and 2mu2tau as well
2319 
2320  const unsigned int nPossibleHandles=6;
2321  MelaPConstant* pchandle[nPossibleHandles]={ 0 };
2322 
2323  float constant_tmp=0;
2324  if (myME_ == TVar::JHUGen){
2325  if (myProduction_ == TVar::ZZGG){
2326  if (
2327  myModel_==TVar::HSMHiggs
2328  ||
2329  myModel_==TVar::H0_g1prime2
2330  ||
2331  myModel_==TVar::H0hplus
2332  ||
2333  myModel_==TVar::H0minus
2334  ||
2335  myModel_==TVar::H0_Zgsg1prime2
2336  ||
2337  myModel_==TVar::H0_Zgs
2338  ||
2339  myModel_==TVar::H0_Zgs_PS
2340  ||
2341  myModel_==TVar::H0_gsgs
2342  ||
2343  myModel_==TVar::H0_gsgs_PS
2344  ||
2345  myModel_==TVar::SelfDefine_spin0
2346  ){
2347  if (is2mu2e) pchandle[0] = pAvgSmooth_JHUGen_ZZGG_HSMHiggs_2mu2e;
2348  else if (is4mu) pchandle[0] = pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4mu;
2349  else if (is4e) pchandle[0] = pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4e;
2350  }
2351  }
2352  }
2353  else if (myME_ == TVar::MCFM){
2354  // ZZQQB
2355  if (myProduction_ == TVar::ZZQQB){
2356  if (myModel_ == TVar::bkgZZ){
2357  if (is2mu2e) pchandle[0] = pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e;
2358  else if (is4mu) pchandle[0] = pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu;
2359  else if (is4e) pchandle[0] = pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e;
2360  }
2361  }
2362  // ZZGG
2363  else if (myProduction_ == TVar::ZZGG){
2364  if (myModel_ == TVar::bkgZZ){
2365  if (is2mu2e) pchandle[0] = pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e;
2366  else if (is4mu) pchandle[0] = pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu;
2367  else if (is4e) pchandle[0] = pAvgSmooth_MCFM_ZZGG_bkgZZ_4e;
2368  }
2369  else if (myModel_ == TVar::HSMHiggs){
2370  if (is2mu2e) pchandle[0] = pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e;
2371  else if (is4mu) pchandle[0] = pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu;
2372  else if (is4e) pchandle[0] = pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e;
2373  }
2374  // Sum signal and bkg
2375  else if (myModel_ == TVar::bkgZZ_SMHiggs){
2376  if (is2mu2e){
2377  pchandle[0] = pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e;
2378  pchandle[1] = pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e;
2379  }
2380  else if (is4mu){
2381  pchandle[0] = pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu;
2382  pchandle[1] = pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu;
2383  }
2384  else if (is4e){
2385  pchandle[0] = pAvgSmooth_MCFM_ZZGG_bkgZZ_4e;
2386  pchandle[1] = pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e;
2387  }
2388  }
2389  }
2390  // JJEW and components
2391  else if (
2392  myProduction_ == TVar::JJVBF || myProduction_ == TVar::Had_WH || myProduction_ == TVar::Had_ZH || myProduction_ == TVar::JJEW
2393  ||
2394  myProduction_ == TVar::JJVBF_S || myProduction_ == TVar::Had_WH_S || myProduction_ == TVar::Had_ZH_S || myProduction_ == TVar::JJEW_S
2395  ||
2396  myProduction_ == TVar::JJVBF_TU || myProduction_ == TVar::Had_WH_TU || myProduction_ == TVar::Had_ZH_TU || myProduction_ == TVar::JJEW_TU
2397  ){
2398  MelaPConstant* hvbf=0;
2399  MelaPConstant* hwh=0;
2400  MelaPConstant* hzh=0;
2401  MelaPConstant* hvbs=0;
2402  MelaPConstant* hwzz=0;
2403  MelaPConstant* hzzz=0;
2404  bool isEW = (myProduction_ == TVar::JJEW || myProduction_ == TVar::JJEW_S || myProduction_ == TVar::JJEW_TU);
2405  bool hasVBF = isEW || (myProduction_ == TVar::JJVBF || myProduction_ == TVar::JJVBF_S || myProduction_ == TVar::JJVBF_TU);
2406  bool hasZH = isEW || (myProduction_ == TVar::Had_ZH || myProduction_ == TVar::Had_ZH_S || myProduction_ == TVar::Had_ZH_TU);
2407  bool hasWH = isEW || (myProduction_ == TVar::Had_WH || myProduction_ == TVar::Had_WH_S || myProduction_ == TVar::Had_WH_TU);
2408  if (is2mu2e){
2409  hvbf = pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_2mu2e;
2410  hwh = pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_2mu2e;
2411  hzh = pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_2mu2e;
2412  hvbs = pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e;
2413  hwzz = pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e;
2414  hzzz = pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e;
2415  }
2416  else if (is4mu){
2417  hvbf = pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4mu;
2418  hwh = pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4mu;
2419  hzh = pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4mu;
2420  hvbs = pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu;
2421  hwzz = pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu;
2422  hzzz = pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu;
2423  }
2424  else if (is4e){
2425  hvbf = pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4e;
2426  hwh = pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4e;
2427  hzh = pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4e;
2428  hvbs = pAvgSmooth_MCFM_JJVBF_bkgZZ_4e;
2429  hwzz = pAvgSmooth_MCFM_Had_WH_bkgZZ_4e;
2430  hzzz = pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e;
2431  }
2432  if (myModel_ == TVar::HSMHiggs || myModel_ == TVar::bkgZZ_SMHiggs){
2433  if (hasVBF) pchandle[0]=hvbf;
2434  if (hasZH) pchandle[1]=hzh;
2435  if (hasWH) pchandle[2]=hwh;
2436  }
2437  if (myModel_ == TVar::bkgZZ || myModel_ == TVar::bkgZZ_SMHiggs){
2438  if (hasVBF) pchandle[3]=hvbs;
2439  if (hasZH) pchandle[4]=hzzz;
2440  if (hasWH) pchandle[5]=hwzz;
2441  }
2442  }
2443  else if (myProduction_ == TVar::JJQCD){
2444  if (myModel_ == TVar::bkgZJets){
2445  pchandle[0] = pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q; // Only option at the moment
2446  }
2447  else if (myModel_ == TVar::bkgZZ){
2448  if (is2mu2e){
2449  pchandle[0] = pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e;
2450  }
2451  else if (is4mu){
2452  pchandle[0] = pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu;
2453  }
2454  else if (is4e){
2455  pchandle[0] = pAvgSmooth_MCFM_JJQCD_bkgZZ_4e;
2456  }
2457  }
2458  }
2459  }
2460 
2461  bool hasNullHandle=true;
2462  for (unsigned int ihandle=0; ihandle<nPossibleHandles; ihandle++){ if (pchandle[ihandle]!=0){ constant_tmp += pchandle[ihandle]->Eval(getIORecord(), myVerbosity_); hasNullHandle=false; } }
2463  if (hasNullHandle) return constant;
2464 
2465  constant = constant_tmp;
2466  return constant;
2467 }

◆ getConstant_JHUGenUndecayed()

float Mela::getConstant_JHUGenUndecayed ( )
protected

Definition at line 2251 of file Mela.cc.

2251  {
2252  float constant = 1;
2253  if (!melaCand) return constant;
2254 
2255  MelaPConstant* pchandle=0;
2256  unsigned int iarray=0;
2257  if (TUtil::JetMassScheme == TVar::ConserveDifermionMass) iarray=0; // First element points to the case when the difermion invariant mass is conserved in mass removal scheme
2258  else if (TUtil::JetMassScheme == TVar::MomentumToEnergy) iarray=1; // Second element points to the case when the 3-momentum vector magnitude is scaled to energy in mass removal scheme
2259  if (myProduction_ == TVar::JQCD){
2260  pchandle = pAvgSmooth_JHUGen_JQCD_HSMHiggs[iarray];
2261  }
2262  else if (myProduction_ == TVar::JJQCD){
2263  pchandle = pAvgSmooth_JHUGen_JJQCD_HSMHiggs[iarray];
2264  }
2265  else if (myProduction_ == TVar::JJVBF){
2266  pchandle = pAvgSmooth_JHUGen_JJVBF_HSMHiggs[iarray];
2267  }
2268  else if (myProduction_ == TVar::Had_ZH){
2269  pchandle = pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[iarray];
2270  }
2271  else if (myProduction_ == TVar::Had_WH){
2272  pchandle = pAvgSmooth_JHUGen_Had_WH_HSMHiggs[iarray];
2273  }
2274  /*
2275  else if (myProduction_ == TVar::Lep_ZH)
2276  else if (myProduction_ == TVar::Lep_WH)
2277  else if (myProduction_ == TVar::GammaH)
2278  else if (myProduction_ == TVar::ttH)
2279  else if (myProduction_ == TVar::bbH)
2280  */
2281  else return constant;
2282 
2283  constant = pchandle->Eval(getIORecord(), myVerbosity_);
2284  return constant;
2285 }

◆ getPConstantHandle()

MelaPConstant * Mela::getPConstantHandle ( TVar::MatrixElement  me_,
TVar::Production  prod_,
TVar::Process  proc_,
TString  relpath,
TString  spname,
const bool  useSqrts = false 
)
protected

Definition at line 2665 of file Mela.cc.

2672  {
2673  if (myVerbosity_>=TVar::DEBUG) MELAout << "Begin Mela::getPConstantHandle" << endl;
2674 
2675  MelaPConstant* pchandle=0;
2676  string cfile_fullpath;
2677 
2678  // Get data/ path
2679  if (myVerbosity_>=TVar::DEBUG) MELAout << "Mela::getPConstantHandle: relpath and spline name: " << relpath << ", " << spname << endl;
2680  const string MELAPKGPATH = TVar::GetMELAPath();
2681  const string path = MELAPKGPATH + "data/";
2682  if (myVerbosity_>=TVar::DEBUG) MELAout << "Mela::getPConstantHandle: path and spline name: " << path << ", " << spname << endl;
2683 
2684  if (useSqrts){ // Loop over possible sqrts values to get the closest one
2685  const unsigned int npossiblesqrts=3;
2686  const double possible_sqrts[npossiblesqrts]={ 7, 8, 13 };
2687  vector<double> trysqrts;
2688  for (unsigned isq=0; isq<npossiblesqrts; isq++){
2689  double val = possible_sqrts[isq];
2690  double diff = fabs(LHCsqrts-val);
2691  bool inserted=false;
2692  for (std::vector<double>::iterator it = trysqrts.begin(); it<trysqrts.end(); it++){
2693  if (fabs((*it)-LHCsqrts)>diff){
2694  inserted=true;
2695  trysqrts.insert(it, val);
2696  break;
2697  }
2698  }
2699  if (!inserted) trysqrts.push_back(val);
2700  }
2701  for (auto& dsqrts:trysqrts){
2702  TString strsqrts = Form("%s_%.0f%s", relpath.Data(), dsqrts, "TeV");
2703  cfile_fullpath = path;
2704  cfile_fullpath.append(strsqrts.Data());
2705  cfile_fullpath.append(".root");
2706  pchandle = new MelaPConstant(me_, prod_, proc_, cfile_fullpath.c_str(), spname.Data());
2707  if (pchandle->IsValid()){
2708  if (myVerbosity_>=TVar::DEBUG) MELAout << "Mela::getPConstantHandle: Full path and spline name: " << cfile_fullpath << ", " << spname << " is valid." << endl;
2709  break;
2710  }
2711  else{
2712  if (myVerbosity_>=TVar::DEBUG) MELAout << "Mela::getPConstantHandle: Full path and spline name: " << cfile_fullpath << ", " << spname << " is invalid." << endl;
2713  deletePConstantHandle(pchandle);
2714  }
2715  }
2716  }
2717  else{
2718  cfile_fullpath = path;
2719  cfile_fullpath.append(relpath.Data());
2720  cfile_fullpath.append(".root");
2721  pchandle = new MelaPConstant(me_, prod_, proc_, cfile_fullpath.c_str(), spname.Data());
2722  if (!pchandle->IsValid()) deletePConstantHandle(pchandle);
2723  }
2724  if (myVerbosity_>=TVar::DEBUG) MELAout << "Mela::getPConstantHandle: Full path and spline name: " << cfile_fullpath << ", " << spname << endl;
2725  if (myVerbosity_>=TVar::DEBUG && pchandle==0) MELAerr << "Mela::getPConstantHandle: Handle of " << spname << " from " << cfile_fullpath << " is invalid!" << endl;
2726  if (myVerbosity_>=TVar::DEBUG) MELAout << "End Mela::getPConstantHandle" << endl;
2727  return pchandle;
2728 }

◆ getPConstantHandles()

void Mela::getPConstantHandles ( )
protected

Definition at line 2470 of file Mela.cc.

2470  {
2471  if (myVerbosity_>=TVar::DEBUG) MELAout << "Begin Mela::getPConstantHandles" << endl;
2472 
2473  // Initialize the handles to 0
2474  for (unsigned int isch=0; isch<(unsigned int)(TVar::nFermionMassRemovalSchemes-1); isch++){
2475  pAvgSmooth_JHUGen_JJQCD_HSMHiggs[isch]=0;
2476  pAvgSmooth_JHUGen_JQCD_HSMHiggs[isch]=0;
2477  pAvgSmooth_JHUGen_JJVBF_HSMHiggs[isch]=0;
2478  pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[isch]=0;
2479  pAvgSmooth_JHUGen_Had_WH_HSMHiggs[isch]=0;
2480  }
2481  //
2482  pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q=0;
2483  //
2484  pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4mu=0;
2485  pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4e=0;
2486  pAvgSmooth_JHUGen_ZZGG_HSMHiggs_2mu2e=0;
2487  //
2488  pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu=0;
2489  pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e=0;
2490  pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e=0;
2491  //
2492  pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4mu=0;
2493  pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4e=0;
2494  pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_2mu2e=0;
2495  //
2496  pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4mu=0;
2497  pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4e=0;
2498  pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_2mu2e=0;
2499  //
2500  pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4mu=0;
2501  pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4e=0;
2502  pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_2mu2e=0;
2503  //
2504  pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu=0;
2505  pAvgSmooth_MCFM_ZZGG_bkgZZ_4e=0;
2506  pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e=0;
2507  //
2508  pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu=0;
2509  pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e=0;
2510  pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e=0;
2511  //
2512  pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu=0;
2513  pAvgSmooth_MCFM_JJVBF_bkgZZ_4e=0;
2514  pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e=0;
2515  //
2516  pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu=0;
2517  pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e=0;
2518  pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e=0;
2519  //
2520  pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu=0;
2521  pAvgSmooth_MCFM_Had_WH_bkgZZ_4e=0;
2522  pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e=0;
2523  //
2524  pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu=0;
2525  pAvgSmooth_MCFM_JJQCD_bkgZZ_4e=0;
2526  pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e=0;
2527  //
2528 
2529  TString filename, spname;
2530 
2531  // Fill versions with difermion correction, set to ConserveDifermionMass if others don't exist.
2532  filename = "pAvgSmooth_JHUGen_JJQCD_HSMHiggs";
2533  spname = "P_ConserveDifermionMass";
2534  pAvgSmooth_JHUGen_JJQCD_HSMHiggs[0] = getPConstantHandle(TVar::JHUGen, TVar::JJQCD, TVar::HSMHiggs, filename, spname, true);
2535  spname = "P_MomentumToEnergy";
2536  pAvgSmooth_JHUGen_JJQCD_HSMHiggs[1] = getPConstantHandle(TVar::JHUGen, TVar::JJQCD, TVar::HSMHiggs, filename, spname, true);
2537  if (pAvgSmooth_JHUGen_JJQCD_HSMHiggs[1]==0) pAvgSmooth_JHUGen_JJQCD_HSMHiggs[1]=pAvgSmooth_JHUGen_JJQCD_HSMHiggs[0];
2538  //
2539  filename = "pAvgSmooth_JHUGen_JQCD_HSMHiggs";
2540  spname = "P_ConserveDifermionMass";
2541  pAvgSmooth_JHUGen_JQCD_HSMHiggs[0] = getPConstantHandle(TVar::JHUGen, TVar::JQCD, TVar::HSMHiggs, filename, spname, true);
2542  spname = "P_MomentumToEnergy";
2543  pAvgSmooth_JHUGen_JQCD_HSMHiggs[1] = getPConstantHandle(TVar::JHUGen, TVar::JQCD, TVar::HSMHiggs, filename, spname, true);
2544  if (pAvgSmooth_JHUGen_JQCD_HSMHiggs[1]==0) pAvgSmooth_JHUGen_JQCD_HSMHiggs[1]=pAvgSmooth_JHUGen_JQCD_HSMHiggs[0];
2545  //
2546  filename = "pAvgSmooth_JHUGen_JJVBF_HSMHiggs";
2547  spname = "P_ConserveDifermionMass";
2548  pAvgSmooth_JHUGen_JJVBF_HSMHiggs[0] = getPConstantHandle(TVar::JHUGen, TVar::JJVBF, TVar::HSMHiggs, filename, spname, true);
2549  spname = "P_MomentumToEnergy";
2550  pAvgSmooth_JHUGen_JJVBF_HSMHiggs[1] = getPConstantHandle(TVar::JHUGen, TVar::JJVBF, TVar::HSMHiggs, filename, spname, true);
2551  if (pAvgSmooth_JHUGen_JJVBF_HSMHiggs[1]==0) pAvgSmooth_JHUGen_JJVBF_HSMHiggs[1]=pAvgSmooth_JHUGen_JJVBF_HSMHiggs[0];
2552  //
2553  filename = "pAvgSmooth_JHUGen_Had_ZH_HSMHiggs";
2554  spname = "P_ConserveDifermionMass";
2555  pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[0] = getPConstantHandle(TVar::JHUGen, TVar::Had_ZH, TVar::HSMHiggs, filename, spname, true);
2556  spname = "P_MomentumToEnergy";
2557  //pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[1] = getPConstantHandle(TVar::JHUGen, TVar::Had_ZH, TVar::HSMHiggs, filename, spname, true);
2558  if (pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[1]==0) pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[1]=pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[0];
2559  //
2560  filename = "pAvgSmooth_JHUGen_Had_WH_HSMHiggs";
2561  spname = "P_ConserveDifermionMass";
2562  pAvgSmooth_JHUGen_Had_WH_HSMHiggs[0] = getPConstantHandle(TVar::JHUGen, TVar::Had_WH, TVar::HSMHiggs, filename, spname, true);
2563  spname = "P_MomentumToEnergy";
2564  //pAvgSmooth_JHUGen_Had_WH_HSMHiggs[1] = getPConstantHandle(TVar::JHUGen, TVar::Had_WH, TVar::HSMHiggs, filename, spname, true);
2565  if (pAvgSmooth_JHUGen_Had_WH_HSMHiggs[1]==0) pAvgSmooth_JHUGen_Had_WH_HSMHiggs[1]=pAvgSmooth_JHUGen_Had_WH_HSMHiggs[0];
2566  //
2567  //
2568  filename = "pAvgSmooth_MCFM_JJQCD_bkgZJets_13TeV_2l2q"; // 13 TeV is a placeholder for all energies.
2569  spname = "P_ConserveDifermionMass";
2570  pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q = getPConstantHandle(TVar::MCFM, TVar::JJQCD, TVar::bkgZJets, filename, spname);
2571  //
2572  //
2573  filename = "pAvgSmooth_JHUGen_ZZGG_HSMHiggs";
2574  spname = "P_ConserveDifermionMass_4mu";
2575  pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4mu = getPConstantHandle(TVar::JHUGen, TVar::ZZGG, TVar::HSMHiggs, filename, spname);
2576  spname = "P_ConserveDifermionMass_4e";
2577  pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4e = getPConstantHandle(TVar::JHUGen, TVar::ZZGG, TVar::HSMHiggs, filename, spname);
2578  spname = "P_ConserveDifermionMass_2mu2e";
2579  pAvgSmooth_JHUGen_ZZGG_HSMHiggs_2mu2e = getPConstantHandle(TVar::JHUGen, TVar::ZZGG, TVar::HSMHiggs, filename, spname);
2580  //
2581  filename = "pAvgSmooth_MCFM_ZZGG_HSMHiggs";
2582  spname = "P_ConserveDifermionMass_4mu";
2583  pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu = getPConstantHandle(TVar::MCFM, TVar::ZZGG, TVar::HSMHiggs, filename, spname);
2584  spname = "P_ConserveDifermionMass_4e";
2585  pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e = getPConstantHandle(TVar::MCFM, TVar::ZZGG, TVar::HSMHiggs, filename, spname);
2586  spname = "P_ConserveDifermionMass_2mu2e";
2587  pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e = getPConstantHandle(TVar::MCFM, TVar::ZZGG, TVar::HSMHiggs, filename, spname);
2588  //
2589  filename = "pAvgSmooth_MCFM_JJVBF_S_HSMHiggs";
2590  spname = "P_ConserveDifermionMass_4mu";
2591  pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4mu = getPConstantHandle(TVar::MCFM, TVar::JJVBF_S, TVar::HSMHiggs, filename, spname, true);
2592  spname = "P_ConserveDifermionMass_4e";
2593  pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4e = getPConstantHandle(TVar::MCFM, TVar::JJVBF_S, TVar::HSMHiggs, filename, spname, true);
2594  spname = "P_ConserveDifermionMass_2mu2e";
2595  pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_2mu2e = getPConstantHandle(TVar::MCFM, TVar::JJVBF_S, TVar::HSMHiggs, filename, spname, true);
2596  //
2597  filename = "pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs";
2598  spname = "P_ConserveDifermionMass_4mu";
2599  pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4mu = getPConstantHandle(TVar::MCFM, TVar::Had_ZH_S, TVar::HSMHiggs, filename, spname, true);
2600  spname = "P_ConserveDifermionMass_4e";
2601  pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4e = getPConstantHandle(TVar::MCFM, TVar::Had_ZH_S, TVar::HSMHiggs, filename, spname, true);
2602  spname = "P_ConserveDifermionMass_2mu2e";
2603  pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_2mu2e = getPConstantHandle(TVar::MCFM, TVar::Had_ZH_S, TVar::HSMHiggs, filename, spname, true);
2604  //
2605  filename = "pAvgSmooth_MCFM_Had_WH_S_HSMHiggs";
2606  spname = "P_ConserveDifermionMass_4mu";
2607  pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4mu = getPConstantHandle(TVar::MCFM, TVar::Had_WH_S, TVar::HSMHiggs, filename, spname, true);
2608  spname = "P_ConserveDifermionMass_4e";
2609  pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4e = getPConstantHandle(TVar::MCFM, TVar::Had_WH_S, TVar::HSMHiggs, filename, spname, true);
2610  spname = "P_ConserveDifermionMass_2mu2e";
2611  pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_2mu2e = getPConstantHandle(TVar::MCFM, TVar::Had_WH_S, TVar::HSMHiggs, filename, spname, true);
2612  //
2613  //
2614  filename = "pAvgSmooth_MCFM_ZZGG_bkgZZ";
2615  spname = "P_ConserveDifermionMass_4mu";
2616  pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu = getPConstantHandle(TVar::MCFM, TVar::ZZGG, TVar::bkgZZ, filename, spname);
2617  spname = "P_ConserveDifermionMass_4e";
2618  pAvgSmooth_MCFM_ZZGG_bkgZZ_4e = getPConstantHandle(TVar::MCFM, TVar::ZZGG, TVar::bkgZZ, filename, spname);
2619  spname = "P_ConserveDifermionMass_2mu2e";
2620  pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e = getPConstantHandle(TVar::MCFM, TVar::ZZGG, TVar::bkgZZ, filename, spname);
2621  //
2622  filename = "pAvgSmooth_MCFM_ZZQQB_bkgZZ";
2623  spname = "P_ConserveDifermionMass_4mu";
2624  pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu = getPConstantHandle(TVar::MCFM, TVar::ZZQQB, TVar::bkgZZ, filename, spname);
2625  spname = "P_ConserveDifermionMass_4e";
2626  pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e = getPConstantHandle(TVar::MCFM, TVar::ZZQQB, TVar::bkgZZ, filename, spname);
2627  spname = "P_ConserveDifermionMass_2mu2e";
2628  pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e = getPConstantHandle(TVar::MCFM, TVar::ZZQQB, TVar::bkgZZ, filename, spname);
2629  //
2630  filename = "pAvgSmooth_MCFM_JJVBF_bkgZZ";
2631  spname = "P_ConserveDifermionMass_4mu";
2632  pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu = getPConstantHandle(TVar::MCFM, TVar::JJVBF, TVar::bkgZZ, filename, spname, true);
2633  spname = "P_ConserveDifermionMass_4e";
2634  pAvgSmooth_MCFM_JJVBF_bkgZZ_4e = getPConstantHandle(TVar::MCFM, TVar::JJVBF, TVar::bkgZZ, filename, spname, true);
2635  spname = "P_ConserveDifermionMass_2mu2e";
2636  pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e = getPConstantHandle(TVar::MCFM, TVar::JJVBF, TVar::bkgZZ, filename, spname, true);
2637  //
2638  filename = "pAvgSmooth_MCFM_Had_ZH_bkgZZ";
2639  spname = "P_ConserveDifermionMass_4mu";
2640  pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu = getPConstantHandle(TVar::MCFM, TVar::Had_ZH, TVar::bkgZZ, filename, spname, true);
2641  spname = "P_ConserveDifermionMass_4e";
2642  pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e = getPConstantHandle(TVar::MCFM, TVar::Had_ZH, TVar::bkgZZ, filename, spname, true);
2643  spname = "P_ConserveDifermionMass_2mu2e";
2644  pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e = getPConstantHandle(TVar::MCFM, TVar::Had_ZH, TVar::bkgZZ, filename, spname, true);
2645  //
2646  filename = "pAvgSmooth_MCFM_Had_WH_bkgZZ";
2647  spname = "P_ConserveDifermionMass_4mu";
2648  pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu = getPConstantHandle(TVar::MCFM, TVar::Had_WH, TVar::bkgZZ, filename, spname, true);
2649  spname = "P_ConserveDifermionMass_4e";
2650  pAvgSmooth_MCFM_Had_WH_bkgZZ_4e = getPConstantHandle(TVar::MCFM, TVar::Had_WH, TVar::bkgZZ, filename, spname, true);
2651  spname = "P_ConserveDifermionMass_2mu2e";
2652  pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e = getPConstantHandle(TVar::MCFM, TVar::Had_WH, TVar::bkgZZ, filename, spname, true);
2653  //
2654  filename = "pAvgSmooth_MCFM_JJQCD_bkgZZ";
2655  spname = "P_ConserveDifermionMass_4mu";
2656  pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu = getPConstantHandle(TVar::MCFM, TVar::JJQCD, TVar::bkgZZ, filename, spname, true);
2657  spname = "P_ConserveDifermionMass_4e";
2658  pAvgSmooth_MCFM_JJQCD_bkgZZ_4e = getPConstantHandle(TVar::MCFM, TVar::JJQCD, TVar::bkgZZ, filename, spname, true);
2659  spname = "P_ConserveDifermionMass_2mu2e";
2660  pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e = getPConstantHandle(TVar::MCFM, TVar::JJQCD, TVar::bkgZZ, filename, spname, true);
2661  //
2662 
2663  if (myVerbosity_>=TVar::DEBUG) MELAout << "End Mela::getPConstantHandles" << endl;
2664 }

◆ printLogo()

void Mela::printLogo ( ) const
protected

Definition at line 278 of file Mela.cc.

278  {
279  vector<string> logolines;
280  logolines.push_back("MELA (Matrix Element Likelihood Approach)");
281  logolines.push_back("");
282  logolines.push_back("Data analysis and Monte Carlo weights package");
283  logolines.push_back("for analyses of resonances produced at pp, ppbar, and e+e- colliders, featuring:");
284  logolines.push_back("");
285  logolines.push_back("* JHUGenMELA *");
286  logolines.push_back("Signal calculations based on analytical pdf.s, and JHU Generator (JHUGen) matrix elements");
287  logolines.push_back("(See JHUGen credits below)");
288  logolines.push_back("");
289  logolines.push_back("* MCFM *");
290  logolines.push_back("Signal, background, and interference calculations, modified based on JHUGen matrix elements");
291  logolines.push_back("(See MCFM credits below)");
292  logolines.push_back("");
293  logolines.push_back("For more details: http://spin.pha.jhu.edu");
294  logolines.push_back("");
295  size_t maxlinesize = 0;
296  for (auto const& l:logolines) maxlinesize = std::max(maxlinesize, l.length());
297  MELAout.writeCentered("", '*', maxlinesize+10); MELAout << endl;
298  unsigned int iline=0;
299  for (auto const& l:logolines){
300  MELAout << '*';
301  MELAout.writeCentered(l, ' ', maxlinesize+8);
302  MELAout << '*' << endl;
303  if (iline==0){ MELAout.writeCentered("", '*', maxlinesize+10); MELAout << endl; }
304  iline++;
305  }
306  MELAout.writeCentered("", '*', maxlinesize+10); MELAout << endl;
307 }

◆ reset_CandRef()

void Mela::reset_CandRef ( )
protected

Definition at line 550 of file Mela.cc.

550 { melaCand=0; }

◆ reset_PAux()

void Mela::reset_PAux ( )
protected

Definition at line 555 of file Mela.cc.

555 { auxiliaryProb=1.; } // SuperProb reset

◆ reset_SelfDCouplings()

void Mela::reset_SelfDCouplings ( )
protected

Definition at line 426 of file Mela.cc.

426  {
427  // We have a lot of them, now even more!
428 
429  //****Spin-0****//
430  differentiate_HWW_HZZ=false;
431  // Loop over the number of supported resonances
432  for (int jh=0; jh<(int)nSupportedHiggses; jh++){
433  for (int im=0; im<2; im++){
434  for (int ic=0; ic<SIZE_HGG; ic++) selfDHggcoupl[jh][ic][im]=0;
435  for (int ic=0; ic<SIZE_HGG; ic++) selfDHg4g4coupl[jh][ic][im]=0;
436 
437  for (int ic=0; ic<SIZE_HQQ; ic++) selfDHqqcoupl[jh][ic][im]=0;
438  for (int ic=0; ic<SIZE_HQQ; ic++) selfDHbbcoupl[jh][ic][im]=0;
439  for (int ic=0; ic<SIZE_HQQ; ic++) selfDHttcoupl[jh][ic][im]=0;
440  for (int ic=0; ic<SIZE_HQQ; ic++) selfDHb4b4coupl[jh][ic][im]=0;
441  for (int ic=0; ic<SIZE_HQQ; ic++) selfDHt4t4coupl[jh][ic][im]=0;
442 
443  for (int ic=0; ic<SIZE_HVV; ic++){
444  selfDHzzcoupl[jh][ic][im] = 0;
445  selfDHwwcoupl[jh][ic][im] = 0;
446  }
447  }
448  for (int ik=0; ik<SIZE_HVV_CQSQ; ik++){
449  selfDHzzCLambda_qsq[jh][ik]=0;
450  selfDHwwCLambda_qsq[jh][ik]=0;
451  for (int ic=0; ic<SIZE_HVV_LAMBDAQSQ; ic++){ // These default values do not matter as long as the c's are 0.
452  selfDHzzLambda_qsq[jh][ic][ik] = 100.;
453  selfDHwwLambda_qsq[jh][ic][ik] = 100.;
454  }
455  }
456  }
457  // Spin-0 contact terms
458  for (int im=0; im<2; im++){
459  for (int ic=0; ic<SIZE_HVV; ic++){
460  selfDHzzpcoupl[ic][im] = 0;
461  selfDHzpzpcoupl[ic][im] = 0;
462  selfDHwwpcoupl[ic][im] = 0;
463  selfDHwpwpcoupl[ic][im] = 0;
464  }
465  }
466 
467  //****Spin-1****//
468  for (int im=0; im<2; im++){
469  for (int ic=0; ic<SIZE_ZVV; ic++) selfDZvvcoupl[ic][im] = 0;
470  for (int ic=0; ic<SIZE_ZQQ; ic++) selfDZqqcoupl[ic][im] = 0;
471  }
472 
473  //****Spin-2****//
474  for (int im=0; im<2; im++){
475  for (int ic=0; ic<SIZE_GVV; ic++){
476  selfDGvvcoupl[ic][im] = 0;
477  selfDGvvpcoupl[ic][im] = 0;
478  selfDGvpvpcoupl[ic][im] = 0;
479  }
480  for (int ic=0; ic<SIZE_GGG; ic++) selfDGggcoupl[ic][im] = 0;
481  for (int ic=0; ic<SIZE_GQQ; ic++) selfDGqqcoupl[ic][im] = 0;
482  }
483 
484  // Vprime / contact couplings
485  for (int im=0; im<2; im++){
486  for (int ic=0; ic<SIZE_Vpff; ic++){
487  selfDZpffcoupl[ic][im] = 0;
488  selfDWpffcoupl[ic][im] = 0;
489  }
490  }
491  selfDM_Zprime = -1;
492  selfDGa_Zprime = 0;
493  selfDM_Wprime = -1;
494  selfDGa_Wprime = 0;
495 
496  // aTQGC couplings
497  for (int im=0; im<2; im++){
498  for (int ic=0; ic<SIZE_ATQGC; ic++) selfDaTQGCcoupl[ic][im] = 0;
499  }
500 
501  // AZff couplings
502  for (int im=0; im<2; im++){
503  for (int ic=0; ic<SIZE_AZff; ic++) selfDAZffcoupl[ic][im] = 0;
504  }
505 
506  // Did I tell you that we have a lot of them?
507 }

◆ setATQGCCouplings()

void Mela::setATQGCCouplings ( )
protected

Definition at line 413 of file Mela.cc.

413  {
414  ZZME->set_aTQGCCouplings(
415  selfDaTQGCcoupl
416  );
417 }

◆ setAZffCouplings()

void Mela::setAZffCouplings ( )
protected

Definition at line 418 of file Mela.cc.

418  {
419  ZZME->set_AZffCouplings(
420  selfDAZffcoupl
421  );
422 }

◆ setConstant()

void Mela::setConstant ( )
protected

Definition at line 2175 of file Mela.cc.

2175  {
2176  float constant = 1;
2177  if (!melaCand){ if (myVerbosity_>=TVar::DEBUG) MELAout << "Mela::getConstant: !melaCand" << endl; }
2178  else{
2179  if ( // Undecayed Higgs MEs from JHUGen
2180  myME_ == TVar::JHUGen
2181  &&
2182  (
2183  myProduction_ == TVar::JQCD
2184  ||
2185  myProduction_ == TVar::JJQCD || myProduction_ == TVar::JJVBF
2186  ||
2187  myProduction_ == TVar::Lep_ZH || myProduction_ == TVar::Had_ZH
2188  ||
2189  myProduction_ == TVar::Lep_WH || myProduction_ == TVar::Had_WH
2190  ||
2191  myProduction_ == TVar::GammaH
2192  ||
2193  myProduction_ == TVar::ttH || myProduction_ == TVar::bbH
2194  )
2195  ) constant = getConstant_JHUGenUndecayed();
2196  else if ( // H->4l/2l2l
2197  melaCand->getSortedV(0)->getNDaughters()==2
2198  &&
2199  melaCand->getSortedV(1)->getNDaughters()==2
2200  &&
2201  PDGHelpers::isALepton(melaCand->getSortedV(0)->getDaughter(0)->id) && PDGHelpers::isALepton(melaCand->getSortedV(0)->getDaughter(1)->id)
2202  &&
2203  PDGHelpers::isALepton(melaCand->getSortedV(1)->getDaughter(0)->id) && PDGHelpers::isALepton(melaCand->getSortedV(1)->getDaughter(1)->id)
2204  ) constant = getConstant_4l();
2205  else if ( // H->2l2q
2206  melaCand->getSortedV(0)->getNDaughters()==2
2207  &&
2208  melaCand->getSortedV(1)->getNDaughters()==2
2209  &&
2210  (
2211  (
2212  PDGHelpers::isALepton(melaCand->getSortedV(0)->getDaughter(0)->id) && PDGHelpers::isALepton(melaCand->getSortedV(0)->getDaughter(1)->id)
2213  &&
2214  PDGHelpers::isAJet(melaCand->getSortedV(1)->getDaughter(0)->id) && PDGHelpers::isAJet(melaCand->getSortedV(1)->getDaughter(1)->id)
2215  &&
2216  !PDGHelpers::isAGluon(melaCand->getSortedV(1)->getDaughter(0)->id) && !PDGHelpers::isAGluon(melaCand->getSortedV(1)->getDaughter(1)->id)
2217  )
2218  ||
2219  (
2220  PDGHelpers::isALepton(melaCand->getSortedV(1)->getDaughter(0)->id) && PDGHelpers::isALepton(melaCand->getSortedV(1)->getDaughter(1)->id)
2221  &&
2222  PDGHelpers::isAJet(melaCand->getSortedV(0)->getDaughter(0)->id) && PDGHelpers::isAJet(melaCand->getSortedV(0)->getDaughter(1)->id)
2223  &&
2224  !PDGHelpers::isAGluon(melaCand->getSortedV(0)->getDaughter(0)->id) && !PDGHelpers::isAGluon(melaCand->getSortedV(0)->getDaughter(1)->id)
2225  )
2226  )
2227  ) constant = getConstant_2l2q();
2228  else if ( // H->4q
2229  melaCand->getSortedV(0)->getNDaughters()==2
2230  &&
2231  melaCand->getSortedV(1)->getNDaughters()==2
2232  &&
2233  (
2234  (
2235  PDGHelpers::isAJet(melaCand->getSortedV(0)->getDaughter(0)->id) && PDGHelpers::isAJet(melaCand->getSortedV(0)->getDaughter(1)->id)
2236  &&
2237  PDGHelpers::isAJet(melaCand->getSortedV(1)->getDaughter(0)->id) && PDGHelpers::isAJet(melaCand->getSortedV(1)->getDaughter(1)->id)
2238  &&
2239  !PDGHelpers::isAGluon(melaCand->getSortedV(0)->getDaughter(0)->id) && !PDGHelpers::isAGluon(melaCand->getSortedV(0)->getDaughter(1)->id)
2240  &&
2241  !PDGHelpers::isAGluon(melaCand->getSortedV(1)->getDaughter(0)->id) && !PDGHelpers::isAGluon(melaCand->getSortedV(1)->getDaughter(1)->id)
2242  )
2243  )
2244  ) constant = getConstant_4q();
2245  }
2246  if (std::isnan(constant) || std::isinf(constant) || constant<=0.) constant=0;
2247  else constant=1./constant;
2248  if (myVerbosity_>=TVar::DEBUG) MELAout << "Mela::getConstant: Constant is " << constant << endl;
2249  getIORecord()->setMEConst(constant);
2250 }

◆ setSpinOneCouplings()

void Mela::setSpinOneCouplings ( )
protected

Definition at line 398 of file Mela.cc.

398  {
399  ZZME->set_SpinOneCouplings(selfDZqqcoupl, selfDZvvcoupl);
400 }

◆ setSpinTwoCouplings()

void Mela::setSpinTwoCouplings ( )
protected

Definition at line 401 of file Mela.cc.

401  {
402  ZZME->set_SpinTwoCouplings(selfDGqqcoupl, selfDGggcoupl, selfDGvvcoupl);
403  ZZME->set_SpinTwoContact(selfDGvvpcoupl, selfDGvpvpcoupl);
404  ZZME->set_VprimeContactCouplings(
405  selfDZpffcoupl,
406  selfDWpffcoupl,
407  selfDM_Zprime,
408  selfDGa_Zprime,
409  selfDM_Wprime,
410  selfDGa_Wprime
411  );
412 }

◆ setSpinZeroCouplings()

void Mela::setSpinZeroCouplings ( )
protected

Definition at line 366 of file Mela.cc.

366  {
367  ZZME->set_SpinZeroCouplings(
368  selfDHggcoupl,
369  selfDHg4g4coupl,
370  selfDHqqcoupl,
371  selfDHbbcoupl,
372  selfDHttcoupl,
373  selfDHb4b4coupl,
374  selfDHt4t4coupl,
375  selfDHzzcoupl,
376  selfDHwwcoupl,
377  selfDHzzLambda_qsq,
378  selfDHwwLambda_qsq,
379  selfDHzzCLambda_qsq,
380  selfDHwwCLambda_qsq,
381  differentiate_HWW_HZZ
382  );
383  ZZME->set_SpinZeroContact(
384  selfDHzzpcoupl,
385  selfDHzpzpcoupl,
386  selfDHwwpcoupl,
387  selfDHwpwpcoupl
388  );
389  ZZME->set_VprimeContactCouplings(
390  selfDZpffcoupl,
391  selfDWpffcoupl,
392  selfDM_Zprime,
393  selfDGa_Zprime,
394  selfDM_Wprime,
395  selfDGa_Wprime
396  );
397 }

Variable Documentation

◆ auxiliaryProb

float Mela::auxiliaryProb
protected

Definition at line 1023 of file Mela.h.

◆ differentiate_HWW_HZZ

bool Mela::differentiate_HWW_HZZ

Definition at line 980 of file Mela.h.

◆ LHCsqrts

double Mela::LHCsqrts
protected

Definition at line 1012 of file Mela.h.

◆ melaCand

MELACandidate* Mela::melaCand
protected

Definition at line 1025 of file Mela.h.

◆ myLepInterf_

TVar::LeptonInterference Mela::myLepInterf_
protected

Definition at line 1016 of file Mela.h.

◆ myME_

TVar::MatrixElement Mela::myME_
protected

Definition at line 1014 of file Mela.h.

◆ myModel_

TVar::Process Mela::myModel_
protected

Definition at line 1013 of file Mela.h.

◆ myProduction_

TVar::Production Mela::myProduction_
protected

Definition at line 1015 of file Mela.h.

◆ myVerbosity_

TVar::VerbosityLevel Mela::myVerbosity_
protected

Definition at line 1017 of file Mela.h.

◆ pAvgSmooth_JHUGen_Had_WH_HSMHiggs

MelaPConstant* Mela::pAvgSmooth_JHUGen_Had_WH_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
protected

Definition at line 1038 of file Mela.h.

◆ pAvgSmooth_JHUGen_Had_ZH_HSMHiggs

MelaPConstant* Mela::pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
protected

Definition at line 1036 of file Mela.h.

◆ pAvgSmooth_JHUGen_JJQCD_HSMHiggs

MelaPConstant* Mela::pAvgSmooth_JHUGen_JJQCD_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
protected

Definition at line 1032 of file Mela.h.

◆ pAvgSmooth_JHUGen_JJVBF_HSMHiggs

MelaPConstant* Mela::pAvgSmooth_JHUGen_JJVBF_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
protected

Definition at line 1034 of file Mela.h.

◆ pAvgSmooth_JHUGen_JQCD_HSMHiggs

MelaPConstant* Mela::pAvgSmooth_JHUGen_JQCD_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
protected

Definition at line 1030 of file Mela.h.

◆ pAvgSmooth_JHUGen_ZZGG_HSMHiggs_2mu2e

MelaPConstant* Mela::pAvgSmooth_JHUGen_ZZGG_HSMHiggs_2mu2e
protected

Definition at line 1043 of file Mela.h.

◆ pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4e

MelaPConstant* Mela::pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4e
protected

Definition at line 1042 of file Mela.h.

◆ pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4mu

MelaPConstant* Mela::pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4mu
protected

Definition at line 1041 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e
protected

Definition at line 1079 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_WH_bkgZZ_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_WH_bkgZZ_4e
protected

Definition at line 1078 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu
protected

Definition at line 1077 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_2mu2e
protected

Definition at line 1059 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4e
protected

Definition at line 1058 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4mu
protected

Definition at line 1057 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e
protected

Definition at line 1075 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e
protected

Definition at line 1074 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu
protected

Definition at line 1073 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_2mu2e
protected

Definition at line 1055 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4e
protected

Definition at line 1054 of file Mela.h.

◆ pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4mu
protected

Definition at line 1053 of file Mela.h.

◆ pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q

MelaPConstant* Mela::pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q
protected

Definition at line 1085 of file Mela.h.

◆ pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e
protected

Definition at line 1083 of file Mela.h.

◆ pAvgSmooth_MCFM_JJQCD_bkgZZ_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_JJQCD_bkgZZ_4e
protected

Definition at line 1082 of file Mela.h.

◆ pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu
protected

Definition at line 1081 of file Mela.h.

◆ pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e
protected

Definition at line 1071 of file Mela.h.

◆ pAvgSmooth_MCFM_JJVBF_bkgZZ_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_JJVBF_bkgZZ_4e
protected

Definition at line 1070 of file Mela.h.

◆ pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu
protected

Definition at line 1069 of file Mela.h.

◆ pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_2mu2e
protected

Definition at line 1051 of file Mela.h.

◆ pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4e
protected

Definition at line 1050 of file Mela.h.

◆ pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4mu
protected

Definition at line 1049 of file Mela.h.

◆ pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e
protected

Definition at line 1063 of file Mela.h.

◆ pAvgSmooth_MCFM_ZZGG_bkgZZ_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_ZZGG_bkgZZ_4e
protected

Definition at line 1062 of file Mela.h.

◆ pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu
protected

Definition at line 1061 of file Mela.h.

◆ pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e
protected

Definition at line 1047 of file Mela.h.

◆ pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e
protected

Definition at line 1046 of file Mela.h.

◆ pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu
protected

Definition at line 1045 of file Mela.h.

◆ pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e

MelaPConstant* Mela::pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e
protected

Definition at line 1067 of file Mela.h.

◆ pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e

MelaPConstant* Mela::pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e
protected

Definition at line 1066 of file Mela.h.

◆ pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu

MelaPConstant* Mela::pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu
protected

Definition at line 1065 of file Mela.h.

◆ selfDaTQGCcoupl

double Mela::selfDaTQGCcoupl[SIZE_ATQGC][2]

Definition at line 1001 of file Mela.h.

◆ selfDAZffcoupl

double Mela::selfDAZffcoupl[SIZE_AZff][2]

Definition at line 1003 of file Mela.h.

◆ selfDGa_Wprime

double Mela::selfDGa_Wprime

Definition at line 990 of file Mela.h.

◆ selfDGa_Zprime

double Mela::selfDGa_Zprime

Definition at line 988 of file Mela.h.

◆ selfDGggcoupl

double Mela::selfDGggcoupl[SIZE_GGG][2]

Definition at line 996 of file Mela.h.

◆ selfDGqqcoupl

double Mela::selfDGqqcoupl[SIZE_GQQ][2]

Definition at line 995 of file Mela.h.

◆ selfDGvpvpcoupl

double Mela::selfDGvpvpcoupl[SIZE_GVV][2]

Definition at line 999 of file Mela.h.

◆ selfDGvvcoupl

double Mela::selfDGvvcoupl[SIZE_GVV][2]

Definition at line 997 of file Mela.h.

◆ selfDGvvpcoupl

double Mela::selfDGvvpcoupl[SIZE_GVV][2]

Definition at line 998 of file Mela.h.

◆ selfDHb4b4coupl

double Mela::selfDHb4b4coupl[nSupportedHiggses][SIZE_HQQ][2]

Definition at line 972 of file Mela.h.

◆ selfDHbbcoupl

double Mela::selfDHbbcoupl[nSupportedHiggses][SIZE_HQQ][2]

Definition at line 970 of file Mela.h.

◆ selfDHg4g4coupl

double Mela::selfDHg4g4coupl[nSupportedHiggses][SIZE_HGG][2]

Couplings between the Higgs and PLACEHOLDER

Definition at line 967 of file Mela.h.

◆ selfDHggcoupl

double Mela::selfDHggcoupl[nSupportedHiggses][SIZE_HGG][2]

Couplings between the higgs and gluons.

Definition at line 965 of file Mela.h.

◆ selfDHqqcoupl

double Mela::selfDHqqcoupl[nSupportedHiggses][SIZE_HQQ][2]

Couplings between the Higgs and quarks

Definition at line 969 of file Mela.h.

◆ selfDHt4t4coupl

double Mela::selfDHt4t4coupl[nSupportedHiggses][SIZE_HQQ][2]

Definition at line 973 of file Mela.h.

◆ selfDHttcoupl

double Mela::selfDHttcoupl[nSupportedHiggses][SIZE_HQQ][2]

Definition at line 971 of file Mela.h.

◆ selfDHwpwpcoupl

double Mela::selfDHwpwpcoupl[SIZE_HVV][2]

Definition at line 985 of file Mela.h.

◆ selfDHwwCLambda_qsq

int Mela::selfDHwwCLambda_qsq[nSupportedHiggses][SIZE_HVV_CQSQ]

Definition at line 979 of file Mela.h.

◆ selfDHwwcoupl

double Mela::selfDHwwcoupl[nSupportedHiggses][SIZE_HVV][2]

Definition at line 975 of file Mela.h.

◆ selfDHwwLambda_qsq

double Mela::selfDHwwLambda_qsq[nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]

Definition at line 977 of file Mela.h.

◆ selfDHwwpcoupl

double Mela::selfDHwwpcoupl[SIZE_HVV][2]

Definition at line 984 of file Mela.h.

◆ selfDHzpzpcoupl

double Mela::selfDHzpzpcoupl[SIZE_HVV][2]

Definition at line 982 of file Mela.h.

◆ selfDHzzCLambda_qsq

int Mela::selfDHzzCLambda_qsq[nSupportedHiggses][SIZE_HVV_CQSQ]

Definition at line 978 of file Mela.h.

◆ selfDHzzcoupl

double Mela::selfDHzzcoupl[nSupportedHiggses][SIZE_HVV][2]

Definition at line 974 of file Mela.h.

◆ selfDHzzLambda_qsq

double Mela::selfDHzzLambda_qsq[nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]

Definition at line 976 of file Mela.h.

◆ selfDHzzpcoupl

double Mela::selfDHzzpcoupl[SIZE_HVV][2]

Definition at line 981 of file Mela.h.

◆ selfDM_Wprime

double Mela::selfDM_Wprime

Definition at line 989 of file Mela.h.

◆ selfDM_Zprime

double Mela::selfDM_Zprime

Definition at line 987 of file Mela.h.

◆ selfDWpffcoupl

double Mela::selfDWpffcoupl[SIZE_Vpff][2]

Definition at line 986 of file Mela.h.

◆ selfDZpffcoupl

double Mela::selfDZpffcoupl[SIZE_Vpff][2]

Definition at line 983 of file Mela.h.

◆ selfDZqqcoupl

double Mela::selfDZqqcoupl[SIZE_ZQQ][2]

Definition at line 992 of file Mela.h.

◆ selfDZvvcoupl

double Mela::selfDZvvcoupl[SIZE_ZVV][2]

Definition at line 993 of file Mela.h.

◆ superDijet

SuperDijetMela* Mela::superDijet
protected

Definition at line 1020 of file Mela.h.

◆ ZZME

ZZMatrixElement* Mela::ZZME
protected

Definition at line 1019 of file Mela.h.

TVar::ERROR
@ ERROR
Definition: TVar.hh:49
Mela::pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4e
MelaPConstant * pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4e
Definition: Mela.h:1054
Mela::pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e
MelaPConstant * pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e
Definition: Mela.h:1074
RooSpinZero::modelCouplings::gzgs4List
RooAbsReal * gzgs4List[1][2]
Definition: RooSpinZero.h:19
Mela::getConstant
void getConstant(float &prob)
This function returns a multiplicative constant to the matrix element calculation.
Definition: Mela.cc:2168
TensorPdfFactory::couplings
RooSpinTwo::modelCouplings couplings
Definition: TensorPdfFactory.h:12
TVar::JJVBF
@ JJVBF
Definition: TVar.hh:72
PDGHelpers::isALepton
bool isALepton(const int id)
Definition: PDGHelpers.cc:62
TVar::D_zzzg_g1prime2
@ D_zzzg_g1prime2
Definition: TVar.hh:145
Mela::getConstant_FourFermionDecay
float getConstant_FourFermionDecay(const int &decid)
Definition: Mela.cc:2313
Mela::pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu
MelaPConstant * pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu
Definition: Mela.h:1065
TVar::Had_ZH_S
@ Had_ZH_S
Definition: TVar.hh:86
TVar::D_g1g4
@ D_g1g4
Definition: TVar.hh:139
MELAParticle::getDaughter
MELAParticle * getDaughter(int index) const
Definition: MELAParticle.cc:68
Mela::pAvgSmooth_JHUGen_JQCD_HSMHiggs
MelaPConstant * pAvgSmooth_JHUGen_JQCD_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
Definition: Mela.h:1030
SIZE_GVV
@ SIZE_GVV
Definition: raw_couplings.txt:158
Mela::pAvgSmooth_MCFM_JJQCD_bkgZZ_4e
MelaPConstant * pAvgSmooth_MCFM_JJQCD_bkgZZ_4e
Definition: Mela.h:1082
TVar::H0_Zgs_PS
@ H0_Zgs_PS
Definition: TVar.hh:132
TensorPdfFactory::resetHypotheses
virtual void resetHypotheses()
Definition: TensorPdfFactory.cc:70
Mela::selfDGqqcoupl
double selfDGqqcoupl[SIZE_GQQ][2]
Definition: Mela.h:995
Mela::selfDGvpvpcoupl
double selfDGvpvpcoupl[SIZE_GVV][2]
Definition: Mela.h:999
TensorPdfFactory::setTensorPolarization
virtual void setTensorPolarization(int ig, double initval)
Definition: TensorPdfFactory.cc:63
SIZE_HVV_CQSQ
@ SIZE_HVV_CQSQ
Definition: raw_couplings.txt:74
Mela::selfDM_Zprime
double selfDM_Zprime
Definition: Mela.h:987
Mela::getPConstantHandle
MelaPConstant * getPConstantHandle(TVar::MatrixElement me_, TVar::Production prod_, TVar::Process proc_, TString relpath, TString spname, const bool useSqrts=false)
Definition: Mela.cc:2665
TVar::H0minus
@ H0minus
Definition: TVar.hh:129
TensorPdfFactory::addHypothesis
virtual void addHypothesis(int ig, double initval, double iphase=0)
Definition: TensorPdfFactory.cc:56
Mela::pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e
MelaPConstant * pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e
Definition: Mela.h:1075
MELAParticle::getNDaughters
int getNDaughters() const
Definition: MELAParticle.h:50
Mela::pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e
MelaPConstant * pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e
Definition: Mela.h:1046
Mela::selfDHwwCLambda_qsq
int selfDHwwCLambda_qsq[nSupportedHiggses][SIZE_HVV_CQSQ]
Definition: Mela.h:979
Mela::pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4mu
MelaPConstant * pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4mu
Definition: Mela.h:1049
TVar::nFermionMassRemovalSchemes
@ nFermionMassRemovalSchemes
Definition: TVar.hh:115
Mela::auxiliaryProb
float auxiliaryProb
Definition: Mela.h:1023
SIZE_GGG
@ SIZE_GGG
Definition: raw_couplings.txt:131
Mela::getIORecord
MelaIO * getIORecord()
Returns the MELAIO object, and by consequence, the entire parton-by-parton matrix element record.
Definition: Mela.cc:544
Mela::pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu
MelaPConstant * pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu
Definition: Mela.h:1077
Mela::pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_2mu2e
MelaPConstant * pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_2mu2e
Definition: Mela.h:1059
MELACandidate::getSortedV
MELAParticle * getSortedV(int index) const
Definition: MELACandidate.cc:152
Mela::pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu
MelaPConstant * pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu
Definition: Mela.h:1061
Mela::selfDHttcoupl
double selfDHttcoupl[nSupportedHiggses][SIZE_HQQ][2]
Definition: Mela.h:971
Mela::pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e
MelaPConstant * pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e
Definition: Mela.h:1083
TVar::GammaH
@ GammaH
Definition: TVar.hh:102
RooSpinZero::modelCouplings::g2List
RooAbsReal * g2List[8][2]
Definition: RooSpinZero.h:12
RooSpinZero::modelCouplings::gzgs3List
RooAbsReal * gzgs3List[1][2]
Definition: RooSpinZero.h:18
TVar::H2_g10
@ H2_g10
Definition: TVar.hh:161
Mela::pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4mu
MelaPConstant * pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4mu
Definition: Mela.h:1057
TVar::SelfDefine_spin2
@ SelfDefine_spin2
Definition: TVar.hh:182
Mela::pAvgSmooth_JHUGen_Had_WH_HSMHiggs
MelaPConstant * pAvgSmooth_JHUGen_Had_WH_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
Definition: Mela.h:1038
Mela::pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e
MelaPConstant * pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e
Definition: Mela.h:1066
Mela::selfDHzpzpcoupl
double selfDHzpzpcoupl[SIZE_HVV][2]
Definition: Mela.h:982
SIZE_ZQQ
@ SIZE_ZQQ
Definition: raw_couplings.txt:107
TVar::bkgZJets
@ bkgZJets
Definition: TVar.hh:165
VectorPdfFactory::g1Val
RooRealVar * g1Val
Definition: VectorPdfFactory.h:24
Mela::selfDGa_Zprime
double selfDGa_Zprime
Definition: Mela.h:988
TVar::JJQCD
@ JJQCD
Definition: TVar.hh:71
Mela::selfDWpffcoupl
double selfDWpffcoupl[SIZE_Vpff][2]
Definition: Mela.h:986
Mela::getConstant_4l
float getConstant_4l()
Definition: Mela.cc:2286
Mela::selfDZpffcoupl
double selfDZpffcoupl[SIZE_Vpff][2]
Definition: Mela.h:983
Mela::selfDaTQGCcoupl
double selfDaTQGCcoupl[SIZE_ATQGC][2]
Definition: Mela.h:1001
Mela::getConstant_JHUGenUndecayed
float getConstant_JHUGenUndecayed()
Definition: Mela.cc:2251
TVar::ZZQQB
@ ZZQQB
Definition: TVar.hh:62
ZZMatrixElement::set_SpinZeroContact
void set_SpinZeroContact(double selfDHzzpcoupl[SIZE_HVV][2], double selfDHzpzpcoupl[SIZE_HVV][2], double selfDHwwpcoupl[SIZE_HVV][2], double selfDHwpwpcoupl[SIZE_HVV][2])
Definition: ZZMatrixElement.cc:277
Mela::pAvgSmooth_JHUGen_ZZGG_HSMHiggs_2mu2e
MelaPConstant * pAvgSmooth_JHUGen_ZZGG_HSMHiggs_2mu2e
Definition: Mela.h:1043
ZZMatrixElement::set_SpinOneCouplings
void set_SpinOneCouplings(double selfDZqqcoupl[SIZE_ZQQ][2], double selfDZvvcoupl[SIZE_ZVV][2])
Definition: ZZMatrixElement.cc:290
TVar::D_g1g4_pi_2
@ D_g1g4_pi_2
Definition: TVar.hh:140
TVar::D_zzzg_PS
@ D_zzzg_PS
Definition: TVar.hh:143
modmisc::isnan
logical function isnan(x)
Definition: mod_Misc.F90:380
TVar::D_zzzg_g1prime2_pi_2
@ D_zzzg_g1prime2_pi_2
Definition: TVar.hh:146
Mela::qqZZmodel
RooqqZZ_JHU_ZgammaZZ_fast * qqZZmodel
Definition: Mela.h:950
Mela::pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4mu
MelaPConstant * pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_4mu
Definition: Mela.h:1053
Mela::selfDHqqcoupl
double selfDHqqcoupl[nSupportedHiggses][SIZE_HQQ][2]
Definition: Mela.h:969
testME_all.iline
int iline
Definition: testME_all.py:166
Mela::setConstant
void setConstant()
Definition: Mela.cc:2175
ZZMatrixElement::set_SpinZeroCouplings
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)
Definition: ZZMatrixElement.cc:238
TVar::H2_g7
@ H2_g7
Definition: TVar.hh:158
SIZE_ZVV
@ SIZE_ZVV
Definition: raw_couplings.txt:114
Mela::deletePConstantHandle
void deletePConstantHandle(MelaPConstant *&handle)
Definition: Mela.cc:2785
Mela::ggSpin0Model
ScalarPdfFactory_HVV * ggSpin0Model
Definition: Mela.h:947
Mela::pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e
MelaPConstant * pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e
Definition: Mela.h:1063
MelaPConstant::IsValid
bool IsValid()
Definition: MelaPConstant.h:43
TVar::ZZGG
@ ZZGG
Definition: TVar.hh:61
Mela::selfDGvvcoupl
double selfDGvvcoupl[SIZE_GVV][2]
Definition: Mela.h:997
TVar::MomentumToEnergy
@ MomentumToEnergy
Definition: TVar.hh:114
ZZMatrixElement::set_AZffCouplings
void set_AZffCouplings(double selfDAZffcoupl[SIZE_AZff][2])
Definition: ZZMatrixElement.cc:335
ZZMatrixElement::set_SpinTwoCouplings
void set_SpinTwoCouplings(double selfDGqqcoupl[SIZE_GQQ][2], double selfDGggcoupl[SIZE_GGG][2], double selfDGvvcoupl[SIZE_GVV][2])
Definition: ZZMatrixElement.cc:297
TVar::H0_Zgsg1prime2
@ H0_Zgsg1prime2
Definition: TVar.hh:130
nSupportedHiggses
@ nSupportedHiggses
Definition: TMCFM.hh:15
MELAStreamHelpers::MELAout
MELAOutputStreamer MELAout
Mela::pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e
MelaPConstant * pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e
Definition: Mela.h:1067
TVar::D_zzgg
@ D_zzgg
Definition: TVar.hh:142
TVar::D_zzgg_PS
@ D_zzgg_PS
Definition: TVar.hh:144
TVar::H2_g5
@ H2_g5
Definition: TVar.hh:155
Mela::pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e
MelaPConstant * pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e
Definition: Mela.h:1079
Mela::differentiate_HWW_HZZ
bool differentiate_HWW_HZZ
Definition: Mela.h:980
Mela::melaCand
MELACandidate * melaCand
Definition: Mela.h:1025
ZZMatrixElement::set_VprimeContactCouplings
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)
Definition: ZZMatrixElement.cc:315
SIZE_HQQ
@ SIZE_HQQ
Definition: raw_couplings.txt:5
ScalarPdfFactory::addHypothesis
virtual void addHypothesis(int ig, int ilam, double iphase=0, double altparam_fracval=0)
Definition: ScalarPdfFactory.cc:736
TVar::H0_g1prime2
@ H0_g1prime2
Definition: TVar.hh:127
Mela::getConstant_2l2q
float getConstant_2l2q()
Definition: Mela.cc:2297
anonymous_namespace{TCouplingsBase.hh}::SIZE_HGG
@ SIZE_HGG
Definition: TCouplingsBase.hh:40
Mela::getConstant_4q
float getConstant_4q()
Definition: Mela.cc:2306
Mela::pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4e
MelaPConstant * pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4e
Definition: Mela.h:1042
Mela::myModel_
TVar::Process myModel_
Definition: Mela.h:1013
TVar::H2_g3
@ H2_g3
Definition: TVar.hh:153
Mela::myME_
TVar::MatrixElement myME_
Definition: Mela.h:1014
MelaPConstant::GetSplineName
TString GetSplineName()
Definition: MelaPConstant.h:46
modparameters::c5
complex(8), parameter, public c5
Definition: mod_Parameters.F90:993
testME_all.int
int
Definition: testME_all.py:13
TVar::H0_gsgs_PS
@ H0_gsgs_PS
Definition: TVar.hh:134
MelaPConstant::Eval
double Eval(const MelaIO *RcdME, TVar::VerbosityLevel verbosity) const
Definition: MELAPConstant.cc:59
RooSpinZero::modelCouplings::Lambda_z4qsq
RooAbsReal * Lambda_z4qsq[SIZE_HVV_CQSQ]
Definition: RooSpinZero.h:35
Mela::selfDHg4g4coupl
double selfDHg4g4coupl[nSupportedHiggses][SIZE_HGG][2]
Definition: Mela.h:967
VectorPdfFactory::configure
int configure(TVar::Process model_)
Definition: VectorPdfFactory.h:70
RooSpinZero::modelCouplings::Lambda_z1qsq
RooAbsReal * Lambda_z1qsq[SIZE_HVV_CQSQ]
Definition: RooSpinZero.h:32
TVar::HSMHiggs
@ HSMHiggs
Definition: TVar.hh:126
TVar::H2_g2
@ H2_g2
Definition: TVar.hh:152
Mela::pAvgSmooth_MCFM_Had_WH_bkgZZ_4e
MelaPConstant * pAvgSmooth_MCFM_Had_WH_bkgZZ_4e
Definition: Mela.h:1078
ScalarPdfFactory_HVV::getPDF
RooSpinZero * getPDF()
Definition: ScalarPdfFactory_HVV.h:22
Mela::ZZME
ZZMatrixElement * ZZME
Definition: Mela.h:1019
TVar::GetMELAPath
std::string GetMELAPath()
Definition: TVar.cc:121
TVar::H2_g9
@ H2_g9
Definition: TVar.hh:160
Mela::pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_2mu2e
MelaPConstant * pAvgSmooth_MCFM_Had_ZH_S_HSMHiggs_2mu2e
Definition: Mela.h:1055
Mela::LHCsqrts
double LHCsqrts
Definition: Mela.h:1012
TVar::H0_gsgs
@ H0_gsgs
Definition: TVar.hh:133
Mela::spin2Model
TensorPdfFactory_ppHVV * spin2Model
Definition: Mela.h:949
Mela::pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4mu
MelaPConstant * pAvgSmooth_JHUGen_ZZGG_HSMHiggs_4mu
Definition: Mela.h:1041
ScalarPdfFactory_HVV::makeParamsConst
void makeParamsConst(bool yesNo=true)
Definition: ScalarPdfFactory_HVV.cc:149
Mela::selfDZqqcoupl
double selfDZqqcoupl[SIZE_ZQQ][2]
Definition: Mela.h:992
RooSpinZero::modelCouplings::ggsgs4List
RooAbsReal * ggsgs4List[1][2]
Definition: RooSpinZero.h:22
Mela::selfDHb4b4coupl
double selfDHb4b4coupl[nSupportedHiggses][SIZE_HQQ][2]
Definition: Mela.h:972
Mela::selfDHwpwpcoupl
double selfDHwpwpcoupl[SIZE_HVV][2]
Definition: Mela.h:985
l
double l[nf]
Definition: TMCFM.hh:274
Mela::pAvgSmooth_JHUGen_JJQCD_HSMHiggs
MelaPConstant * pAvgSmooth_JHUGen_JJQCD_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
Definition: Mela.h:1032
TVar::Had_ZH_TU
@ Had_ZH_TU
Definition: TVar.hh:97
RooSpinZero::modelCouplings::gzgs2List
RooAbsReal * gzgs2List[1][2]
Definition: RooSpinZero.h:17
TVar::MCFM
@ MCFM
Definition: TVar.hh:56
TVar::D_g1g2
@ D_g1g2
Definition: TVar.hh:137
Mela::pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4e
MelaPConstant * pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_4e
Definition: Mela.h:1050
SIZE_Vpff
@ SIZE_Vpff
Definition: raw_couplings.txt:100
Mela::selfDHzzLambda_qsq
double selfDHzzLambda_qsq[nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]
Definition: Mela.h:976
VectorPdfFactory::PDF
RooSpinOne_7D * PDF
Definition: VectorPdfFactory.h:22
TVar::JJEW_S
@ JJEW_S
Definition: TVar.hh:84
ScalarPdfFactory::resetHypotheses
virtual void resetHypotheses()
Definition: ScalarPdfFactory.cc:890
Mela::selfDHwwcoupl
double selfDHwwcoupl[nSupportedHiggses][SIZE_HVV][2]
Definition: Mela.h:975
Mela::pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu
MelaPConstant * pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu
Definition: Mela.h:1069
TVar::ConserveDifermionMass
@ ConserveDifermionMass
Definition: TVar.hh:113
Mela::pAvgSmooth_MCFM_ZZGG_bkgZZ_4e
MelaPConstant * pAvgSmooth_MCFM_ZZGG_bkgZZ_4e
Definition: Mela.h:1062
Mela::selfDGvvpcoupl
double selfDGvvpcoupl[SIZE_GVV][2]
Definition: Mela.h:998
Mela::pdf
RooAbsPdf * pdf
Definition: Mela.h:946
SIZE_ATQGC
@ SIZE_ATQGC
Definition: raw_couplings.txt:176
TensorPdfFactory_ppHVV::getPDF
RooSpinTwo * getPDF()
Definition: TensorPdfFactory_ppHVV.h:15
Mela::selfDAZffcoupl
double selfDAZffcoupl[SIZE_AZff][2]
Definition: Mela.h:1003
TVar::JJVBF_TU
@ JJVBF_TU
Definition: TVar.hh:94
Mela::selfDHzzCLambda_qsq
int selfDHzzCLambda_qsq[nSupportedHiggses][SIZE_HVV_CQSQ]
Definition: Mela.h:978
TVar::H2_g1
@ H2_g1
Definition: TVar.hh:151
Mela::pAvgSmooth_MCFM_JJVBF_bkgZZ_4e
MelaPConstant * pAvgSmooth_MCFM_JJVBF_bkgZZ_4e
Definition: Mela.h:1070
Mela::myProduction_
TVar::Production myProduction_
Definition: Mela.h:1015
TVar::JHUGen
@ JHUGen
Definition: TVar.hh:57
SIZE_HVV_LAMBDAQSQ
@ SIZE_HVV_LAMBDAQSQ
Definition: raw_couplings.txt:66
TVar::Had_ZH
@ Had_ZH
Definition: TVar.hh:75
Mela::selfDZvvcoupl
double selfDZvvcoupl[SIZE_ZVV][2]
Definition: Mela.h:993
TVar::H2_g1g5
@ H2_g1g5
Definition: TVar.hh:156
TVar::H2_g8
@ H2_g8
Definition: TVar.hh:159
VectorPdfFactory::g2Val
RooRealVar * g2Val
Definition: VectorPdfFactory.h:25
Mela::pAvgSmooth_JHUGen_JJVBF_HSMHiggs
MelaPConstant * pAvgSmooth_JHUGen_JJVBF_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
Definition: Mela.h:1034
SpinPdfFactory::makeParamsConst
virtual void makeParamsConst(bool yesNo)
Definition: SpinPdfFactory.cc:154
MELAStreamHelpers::MELAerr
MELAOutputStreamer MELAerr
Mela::selfDM_Wprime
double selfDM_Wprime
Definition: Mela.h:989
TVar::Had_WH_S
@ Had_WH_S
Definition: TVar.hh:87
TVar::JJEW
@ JJEW
Definition: TVar.hh:73
TVar::D_g1g1prime2
@ D_g1g1prime2
Definition: TVar.hh:136
Mela::selfDGa_Wprime
double selfDGa_Wprime
Definition: Mela.h:990
MELAOutputStreamer::writeCentered
void writeCentered(const T &val, char fillch=' ', std::streamsize gapsize=0)
Definition: MELAOutputStreamer.h:131
Mela::selfDHt4t4coupl
double selfDHt4t4coupl[nSupportedHiggses][SIZE_HQQ][2]
Definition: Mela.h:973
Mela::selfDHwwLambda_qsq
double selfDHwwLambda_qsq[nSupportedHiggses][SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]
Definition: Mela.h:977
TVar::bkgZZ_SMHiggs
@ bkgZZ_SMHiggs
Definition: TVar.hh:170
Mela::pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e
MelaPConstant * pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e
Definition: Mela.h:1047
TVar::bkgZZ
@ bkgZZ
Definition: TVar.hh:166
MelaIO::setMEConst
void setMEConst(const double &val)
Definition: MelaIO.h:108
SIZE_AZff
@ SIZE_AZff
Definition: raw_couplings.txt:194
modparameters::c1
complex(8), parameter, public c1
Definition: mod_Parameters.F90:988
TVar::Had_WH_TU
@ Had_WH_TU
Definition: TVar.hh:98
TVar::H0_Zgs
@ H0_Zgs
Definition: TVar.hh:131
Mela::pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu
MelaPConstant * pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu
Definition: Mela.h:1081
TVar::ttH
@ ttH
Definition: TVar.hh:66
Mela::selfDHwwpcoupl
double selfDHwwpcoupl[SIZE_HVV][2]
Definition: Mela.h:984
Mela::myVerbosity_
TVar::VerbosityLevel myVerbosity_
Definition: Mela.h:1017
RooSpinZero::modelCouplings::Lambda_z3qsq
RooAbsReal * Lambda_z3qsq[SIZE_HVV_CQSQ]
Definition: RooSpinZero.h:34
TVar::SelfDefine_spin0
@ SelfDefine_spin0
Definition: TVar.hh:180
Mela::pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu
MelaPConstant * pAvgSmooth_MCFM_ZZGG_HSMHiggs_4mu
Definition: Mela.h:1045
RooSpinZero::modelCouplings::gzgs1List
RooAbsReal * gzgs1List[1][2]
Definition: RooSpinZero.h:16
PDGHelpers::isAGluon
bool isAGluon(const int id)
Definition: PDGHelpers.cc:58
TVar::JQCD
@ JQCD
Definition: TVar.hh:69
RooSpinZero::modelCouplings::g3List
RooAbsReal * g3List[8][2]
Definition: RooSpinZero.h:13
Mela::selfDHzzcoupl
double selfDHzzcoupl[nSupportedHiggses][SIZE_HVV][2]
Definition: Mela.h:974
RooSpinZero::modelCouplings::g1List
RooAbsReal * g1List[8][2]
Definition: RooSpinZero.h:11
TUtil::JetMassScheme
TVar::FermionMassRemoval JetMassScheme
Definition: TUtil.cc:28
TVar::H2_g6
@ H2_g6
Definition: TVar.hh:157
Mela::pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q
MelaPConstant * pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q
Definition: Mela.h:1085
TVar::Lep_ZH
@ Lep_ZH
Definition: TVar.hh:77
MelaPConstant::GetFileName
TString GetFileName()
Definition: MelaPConstant.h:45
ScalarPdfFactory::couplings
RooSpinZero::modelCouplings couplings
Definition: ScalarPdfFactory.h:12
Mela::pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu
MelaPConstant * pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu
Definition: Mela.h:1073
RooSpinZero::modelCouplings::ggsgs3List
RooAbsReal * ggsgs3List[1][2]
Definition: RooSpinZero.h:21
Mela::spin1Model
VectorPdfFactory * spin1Model
Definition: Mela.h:948
Mela::selfDHzzpcoupl
double selfDHzzpcoupl[SIZE_HVV][2]
Definition: Mela.h:981
SIZE_HVV
@ SIZE_HVV
Definition: raw_couplings.txt:57
Mela::pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4e
MelaPConstant * pAvgSmooth_MCFM_Had_WH_S_HSMHiggs_4e
Definition: Mela.h:1058
TVar::DEBUG
@ DEBUG
Definition: TVar.hh:51
ZZMatrixElement::set_aTQGCCouplings
void set_aTQGCCouplings(double selfDaTQGCcoupl[SIZE_ATQGC][2])
Definition: ZZMatrixElement.cc:330
TVar::D_zzzg
@ D_zzzg
Definition: TVar.hh:141
TVar::H2_g4
@ H2_g4
Definition: TVar.hh:154
TVar::D_g1g2_pi_2
@ D_g1g2_pi_2
Definition: TVar.hh:138
MelaPConstant
Definition: MelaPConstant.h:9
Mela::selfDHggcoupl
double selfDHggcoupl[nSupportedHiggses][SIZE_HGG][2]
Definition: Mela.h:965
RooSpinZero::modelCouplings::Lambda_z2qsq
RooAbsReal * Lambda_z2qsq[SIZE_HVV_CQSQ]
Definition: RooSpinZero.h:33
RooSpinZero::modelCouplings::g4List
RooAbsReal * g4List[8][2]
Definition: RooSpinZero.h:14
Mela::pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_2mu2e
MelaPConstant * pAvgSmooth_MCFM_JJVBF_S_HSMHiggs_2mu2e
Definition: Mela.h:1051
ZZMatrixElement::set_SpinTwoContact
void set_SpinTwoContact(double selfDGvvpcoupl[SIZE_GVV][2], double selfDGvpvpcoupl[SIZE_GVV][2])
Definition: ZZMatrixElement.cc:306
TVar::Had_WH
@ Had_WH
Definition: TVar.hh:76
RooSpinZero::modelCouplings::cLambda_qsq
RooAbsReal * cLambda_qsq[SIZE_HVV_CQSQ]
Definition: RooSpinZero.h:36
TVar::JJEW_TU
@ JJEW_TU
Definition: TVar.hh:95
Mela::selfDHbbcoupl
double selfDHbbcoupl[nSupportedHiggses][SIZE_HQQ][2]
Definition: Mela.h:970
MELAParticle::id
int id
Definition: MELAParticle.h:17
TVar::JJVBF_S
@ JJVBF_S
Definition: TVar.hh:83
TVar::Lep_WH
@ Lep_WH
Definition: TVar.hh:78
TVar::SelfDefine_spin1
@ SelfDefine_spin1
Definition: TVar.hh:181
TVar::H0hplus
@ H0hplus
Definition: TVar.hh:128
RooSpinTwo::modelCouplings::bList
RooAbsReal * bList[SIZE_GVV][2]
Definition: RooSpinTwo.h:11
Mela::selfDGggcoupl
double selfDGggcoupl[SIZE_GGG][2]
Definition: Mela.h:996
TVar::bbH
@ bbH
Definition: TVar.hh:67
RooSpinZero::modelCouplings::ggsgs2List
RooAbsReal * ggsgs2List[1][2]
Definition: RooSpinZero.h:20
Mela::pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e
MelaPConstant * pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e
Definition: Mela.h:1071
SIZE_GQQ
@ SIZE_GQQ
Definition: raw_couplings.txt:121
PDGHelpers::isAJet
bool isAJet(const int id)
Definition: PDGHelpers.cc:18
modparameters::c2
complex(8), parameter, public c2
Definition: mod_Parameters.F90:989
Mela::pAvgSmooth_JHUGen_Had_ZH_HSMHiggs
MelaPConstant * pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[TVar::nFermionMassRemovalSchemes-1]
Definition: Mela.h:1036
Generated by   doxygen 1.8.17