TEvtProb Class Reference

#include <TEvtProb.hh>

Collaboration diagram for TEvtProb:

Public Member Functions
	TEvtProb ()
	TEvtProb (const char *pathtoXSW, double ebeam, const char *pathtoPDFSet, int PDFMember=0, TVar::VerbosityLevel verbosity_=TVar::ERROR)
	TEvtProb (const TEvtProb &other)
	~TEvtProb ()
void	Set_LHAgrid (const char *path, int pdfmember=0)
void	SetProcess (TVar::Process proc, TVar::MatrixElement me, TVar::Production prod)
void	SetVerbosity (TVar::VerbosityLevel tmp)
void	SetLeptonInterf (TVar::LeptonInterference tmp)
void	SetCandidateDecayMode (TVar::CandidateDecayMode mode)
void	SetCurrentCandidateFromIndex (unsigned int icand)
void	SetCurrentCandidate (MELACandidate *cand)
void	AllowSeparateWWCouplings (bool doAllow=false)
void	ResetMass (double inmass, int ipart)
void	ResetWidth (double inwidth, int ipart)
void	SetZprimeMassWidth (double inmass, double inwidth)
void	SetWprimeMassWidth (double inmass, double inwidth)
void	ResetQuarkMasses ()
void	ResetMCFM_EWKParameters (double ext_Gf, double ext_aemmz, double ext_mW, double ext_mZ, double ext_xW, int ext_ewscheme=3)
void	ResetCouplings ()
void	SetPrimaryHiggsMass (double mass)
void	SetHiggsMass (double mass, double wHiggs=-1., int whichResonance=-1)
void	SetRenFacScaleMode (TVar::EventScaleScheme renormalizationSch, TVar::EventScaleScheme factorizationSch, double ren_sf, double fac_sf)
void	ResetRenFacScaleMode ()
void	SetInputEvent (SimpleParticleCollection_t *pDaughters, SimpleParticleCollection_t *pAssociated=0, SimpleParticleCollection_t *pMothers=0, bool isGen=false)
void	AppendTopCandidate (SimpleParticleCollection_t *TopDaughters)
void	ResetInputEvent ()
void	ResetIORecord ()
double	XsecCalc_XVV ()
double	XsecCalc_VVXVV ()
double	XsecCalcXJJ ()
double	XsecCalcXJ ()
double	XsecCalc_VX (bool includeHiggsDecay)
double	XsecCalc_TTX (int topProcess, int topDecay)
double	GetXPropagator (TVar::ResonancePropagatorScheme scheme)
MELAHXSWidth const *	GetHXSWidthEstimator () const
SpinZeroCouplings *	GetSelfDSpinZeroCouplings ()
SpinOneCouplings *	GetSelfDSpinOneCouplings ()
SpinTwoCouplings *	GetSelfDSpinTwoCouplings ()
VprimeCouplings *	GetSelfDVprimeCouplings ()
aTQGCCouplings *	GetSelfDaTQGCCouplings ()
AZffCouplings *	GetSelfDAZffCouplings ()
double	GetPrimaryHiggsMass ()
double	GetPrimaryMass (int ipart)
double	GetPrimaryWidth (int ipart)
double	GetHiggsWidthAtPoleMass (double mass)
MelaIO *	GetIORecord ()
MELACandidate *	GetCurrentCandidate ()
int	GetCurrentCandidateIndex ()
int	GetNCandidates ()
std::vector< MELATopCandidate_t * > *	GetTopCandidates ()

Protected Member Functions
void	InitializeMCFM ()
void	InitializeJHUGen (const char *pathtoPDFSet, int PDFMember)
void	CrossInitialize ()
bool	CheckInputPresent ()
void	SetRcdCandPtr ()
bool	CheckSelfDCouplings_Hgg ()
bool	CheckSelfDCouplings_Hqq ()
bool	CheckSelfDCouplings_Htt ()
bool	CheckSelfDCouplings_Hbb ()
bool	CheckSelfDCouplings_HVV ()
bool	CheckSelfDCouplings_aTQGC ()
bool	CheckSelfDCouplings_AZff ()
void	Build ()

Protected Attributes
const char *	pathtoPDFSet_
int	PDFMember_
TVar::Process	process
TVar::MatrixElement	matrixElement
TVar::Production	production
TVar::VerbosityLevel	verbosity
TVar::LeptonInterference	leptonInterf
double	PrimaryHMass
double	_hmass
double	_hwidth
double	_h2mass
double	_h2width
double	EBEAM
MELAHXSWidth	myCSW_
TVar::event_scales_type	event_scales
SpinZeroCouplings	selfDSpinZeroCoupl
SpinOneCouplings	selfDSpinOneCoupl
SpinTwoCouplings	selfDSpinTwoCoupl
VprimeCouplings	selfDVprimeCoupl
aTQGCCouplings	selfDaTQGCCoupl
AZffCouplings	selfDAZffCoupl
MelaIO	RcdME
MELACandidate *	melaCand
std::vector< MELAParticle * >	particleList
std::vector< MELACandidate * >	candList
std::vector< MELATopCandidate_t * >	topCandList

Detailed Description

Definition at line 48 of file TEvtProb.hh.

Constructor & Destructor Documentation

TEvtProb() [1/3]

TEvtProb::TEvtProb ( )

inline

Definition at line 53 of file TEvtProb.hh.

{};

TEvtProb() [2/3]

TEvtProb::TEvtProb	(	const char *	pathtoXSW,
		double	ebeam,
		const char *	pathtoPDFSet,
		int	PDFMember = 0,
		TVar::VerbosityLevel	verbosity_ = TVar::ERROR
	)

Definition at line 26 of file TEvtProb.cc.

    :
  pathtoPDFSet_(pathtoPDFSet),
  PDFMember_(PDFMember),
  verbosity(verbosity_),
  EBEAM(ebeam),
  myCSW_(pathtoXSW)
{
  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb constructor" << endl;
 
  /***** Build everything except MELAHXSWidth *****/
  Build();
 
  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb constructor" << endl;
}

TEvtProb() [3/3]

TEvtProb::TEvtProb

(

const TEvtProb &

other

)

Definition at line 42 of file TEvtProb.cc.

                                        :
pathtoPDFSet_(other.pathtoPDFSet_),
PDFMember_(other.PDFMember_),
verbosity(other.verbosity),
EBEAM(other.EBEAM),
myCSW_(other.myCSW_)
{
 
}

~TEvtProb()

TEvtProb::~TEvtProb

(

)

Definition at line 74 of file TEvtProb.cc.

                   {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb destructor" << endl;
 
  ResetInputEvent();
 
  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb destructor" << endl;
}

Member Function Documentation

AllowSeparateWWCouplings()

void TEvtProb::AllowSeparateWWCouplings

(

bool

doAllow = false

)

Definition at line 215 of file TEvtProb.cc.

{ SetJHUGenDistinguishWWCouplings(doAllow); selfDSpinZeroCoupl.allow_WWZZSeparation(doAllow); }

AppendTopCandidate()

void TEvtProb::AppendTopCandidate

(

SimpleParticleCollection_t *

TopDaughters

)

Definition at line 287 of file TEvtProb.cc.

                                                                         {
  if (!CheckInputPresent()){
    if (verbosity>=TVar::INFO) MELAout << "TEvtProb::AppendTopCandidate: No MELACandidates are present to append this top!" << endl;
    return;
  }
  MELATopCandidate_t* cand = ConvertThreeBodyDecayCandidate(
    TopDaughters,
    &particleList, &topCandList // push_back is done automatically
    );
  if (cand) melaCand->addAssociatedTop(cand);
}

Build()

void TEvtProb::Build ( )

protected

Definition at line 51 of file TEvtProb.cc.

                    {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb::Build" << endl;
 
  /***** Initialize lepton interference scheme *****/
  SetLeptonInterf(TVar::DefaultLeptonInterf);
 
  /***** Initialize MCFM *****/
  InitializeMCFM();
 
  /***** Initialize JHUGen *****/
  InitializeJHUGen(pathtoPDFSet_, PDFMember_);
 
  /***** Cross-initializations *****/
  CrossInitialize();
 
  /***** Initialize input event properties *****/
  ResetInputEvent();
  SetCandidateDecayMode(TVar::CandidateDecay_ZZ);
 
  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb::Build" << endl;
}

CheckInputPresent()

bool TEvtProb::CheckInputPresent ( )

protected

Definition at line 389 of file TEvtProb.cc.

                                {
  if (!melaCand){
    if (verbosity>=TVar::INFO) MELAout
      << "TEvtProb::CheckInputPresent: melaCand==" << melaCand << " is nullPtr!"
      << endl;
    if (candList.empty()) return false;
    else{
      SetCurrentCandidateFromIndex(candList.size()-1);
      if (verbosity>=TVar::INFO) MELAout << "TEvtProb::CheckInputPresent: melaCand now points to the latest candidate (cand[" << (candList.size()-1) << "])" << endl;
    }
  }
  SetRcdCandPtr(); // If input event is present, set the RcdME pointer to melaCand
  return true;
}

CheckSelfDCouplings_aTQGC()

bool TEvtProb::CheckSelfDCouplings_aTQGC ( )

protected

Definition at line 1140 of file TEvtProb.cc.

                                        {
  for (int vv = 0; vv < SIZE_ATQGC; vv++){
    if (
      (selfDaTQGCCoupl.aTQGCcoupl)[vv][1] != 0 || (selfDaTQGCCoupl.aTQGCcoupl)[vv][0] != 0
      ){
      return true;
    }
  }
  return false;
}

CheckSelfDCouplings_AZff()

bool TEvtProb::CheckSelfDCouplings_AZff ( )

protected

Definition at line 1150 of file TEvtProb.cc.

                                       {
  for (int vv = 0; vv < SIZE_AZff; vv++){
    if (
      (selfDAZffCoupl.AZffcoupl)[vv][1] != 0 || (selfDAZffCoupl.AZffcoupl)[vv][0] != 0
      ){
      return true;
    }
  }
  return false;
}

CheckSelfDCouplings_Hbb()

bool TEvtProb::CheckSelfDCouplings_Hbb ( )

protected

Definition at line 1090 of file TEvtProb.cc.

                                      {
  if (_hmass>=0. && _hwidth>0.){
    for (int vv = 0; vv < SIZE_HQQ; vv++){
      if (
        (selfDSpinZeroCoupl.Hbbcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hbbcoupl)[vv][0] != 0
        ||
        (selfDSpinZeroCoupl.Hb4b4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hb4b4coupl)[vv][0] != 0
        ){
        return true;
      }
    }
  }
  if (_h2mass>=0. && _h2width>0.){
    for (int vv = 0; vv < SIZE_HQQ; vv++){
      if (
        (selfDSpinZeroCoupl.H2bbcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2bbcoupl)[vv][0] != 0
        ||
        (selfDSpinZeroCoupl.H2b4b4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2b4b4coupl)[vv][0] != 0
        ){
        return true;
      }
    }
  }
  return false;
}

CheckSelfDCouplings_Hgg()

bool TEvtProb::CheckSelfDCouplings_Hgg ( )

protected

Definition at line 1019 of file TEvtProb.cc.

                                      {
  if (_hmass>=0. && _hwidth>0.){
    for (int vv = 0; vv < SIZE_HGG; vv++){
      if (
        (selfDSpinZeroCoupl.Hggcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hggcoupl)[vv][0] != 0
        ||
        (selfDSpinZeroCoupl.Hg4g4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hg4g4coupl)[vv][0] != 0
        ){
        return true;
      }
    }
  }
  if (_h2mass>=0. && _h2width>0.){
    for (int vv = 0; vv < SIZE_HGG; vv++){
      if (
        (selfDSpinZeroCoupl.H2ggcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2ggcoupl)[vv][0] != 0
        ||
        (selfDSpinZeroCoupl.H2g4g4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2g4g4coupl)[vv][0] != 0
        ){
        return true;
      }
    }
  }
  return false;
}

CheckSelfDCouplings_Hqq()

bool TEvtProb::CheckSelfDCouplings_Hqq ( )

protected

Definition at line 1044 of file TEvtProb.cc.

                                      {
  if (_hmass>=0. && _hwidth>0.){
    for (int vv = 0; vv < SIZE_HQQ; vv++){
      if (
        (selfDSpinZeroCoupl.Hqqcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hqqcoupl)[vv][0] != 0
        ){
        return true;
      }
    }
  }
  if (_h2mass>=0. && _h2width>0.){
    for (int vv = 0; vv < SIZE_HQQ; vv++){
      if (
        (selfDSpinZeroCoupl.H2qqcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2qqcoupl)[vv][0] != 0
        ){
        return true;
      }
    }
  }
  return false;
}

CheckSelfDCouplings_Htt()

bool TEvtProb::CheckSelfDCouplings_Htt ( )

protected

Definition at line 1065 of file TEvtProb.cc.

                                      {
  if (_hmass>=0. && _hwidth>0.){
    for (int vv = 0; vv < SIZE_HQQ; vv++){
      if (
        (selfDSpinZeroCoupl.Httcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Httcoupl)[vv][0] != 0
        ||
        (selfDSpinZeroCoupl.Ht4t4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Ht4t4coupl)[vv][0] != 0
        ){
        return true;
      }
    }
  }
  if (_h2mass>=0. && _h2width>0.){
    for (int vv = 0; vv < SIZE_HQQ; vv++){
      if (
        (selfDSpinZeroCoupl.H2ttcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2ttcoupl)[vv][0] != 0
        ||
        (selfDSpinZeroCoupl.H2t4t4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2t4t4coupl)[vv][0] != 0
        ){
        return true;
      }
    }
  }
  return false;
}

CheckSelfDCouplings_HVV()

bool TEvtProb::CheckSelfDCouplings_HVV ( )

protected

Definition at line 1115 of file TEvtProb.cc.

                                      {
  if (_hmass>=0. && _hwidth>0.){
    for (int vv = 0; vv < SIZE_HVV; vv++){
      if (
        (selfDSpinZeroCoupl.Hzzcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hzzcoupl)[vv][0] != 0
        ||
        (selfDSpinZeroCoupl.Hwwcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hwwcoupl)[vv][0] != 0
        ){
        return true;
      }
    } // No need to check c_q**2. If these are 0, z_q**2 do not have any effect.
  }
  if (_h2mass>=0. && _h2width>0.){
    for (int vv = 0; vv < SIZE_HVV; vv++){
      if (
        (selfDSpinZeroCoupl.H2zzcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2zzcoupl)[vv][0] != 0
        ||
        (selfDSpinZeroCoupl.H2wwcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2wwcoupl)[vv][0] != 0
        ){
        return true;
      }
    } // No need to check c_q**2. If these are 0, z_q**2 do not have any effect.
  }
  return false;
}

CrossInitialize()

void TEvtProb::CrossInitialize ( )

protected

Definition at line 149 of file TEvtProb.cc.

                              {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb::CrossInitialize" << endl;
 
  /**** Initialize alphaS *****/
  const double GeV=1./100.;
  double init_Q;
  int init_nl, init_nf;
  __modjhugenmela_MOD_getpdfconstants(&init_Q, &init_nl, &init_nf);
  init_Q /= GeV;
  if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::TEvtProb: Initializing the PDF with initial Q=" << init_Q << ", nloops=" << init_nl << ", nf=" << init_nf << endl;
  SetAlphaS(init_Q, init_Q, 1, 1, init_nl, init_nf, "cteq6_l"); // "cteq6_l" is just some dummy variable
 
  /**** Initialize MCFM CKM initializers in cabib.f using JHUGen defaults *****/
  int i, j;
  i=2; j=1;
  cabib_.Vud = __modparameters_MOD_ckmbare(&i, &j);
  i=2; j=3;
  cabib_.Vus = __modparameters_MOD_ckmbare(&i, &j);
  i=2; j=5;
  cabib_.Vub = __modparameters_MOD_ckmbare(&i, &j);
  i=4; j=1;
  cabib_.Vcd = __modparameters_MOD_ckmbare(&i, &j);
  i=4; j=3;
  cabib_.Vcs = __modparameters_MOD_ckmbare(&i, &j);
  i=4; j=5;
  cabib_.Vcb = __modparameters_MOD_ckmbare(&i, &j);
 
  /***** Initialize schemes *****/
  ResetCouplings();
  ResetRenFacScaleMode();
  SetPrimaryHiggsMass(125.); // Should come after InitializeMCFM and InitializeJHUGen
 
  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb::CrossInitialize" << endl;
}

GetCurrentCandidate()

MELACandidate * TEvtProb::GetCurrentCandidate

(

)

Definition at line 376 of file TEvtProb.cc.

{ return melaCand; }

GetCurrentCandidateIndex()

int TEvtProb::GetCurrentCandidateIndex

(

)

Definition at line 377 of file TEvtProb.cc.

                                      {
  if (!melaCand) return -1;
  for (unsigned int icand=0; icand<candList.size(); icand++){
    if (candList.at(icand)==melaCand) return (int)icand;
  }
  return -1;
}

GetHiggsWidthAtPoleMass()

double TEvtProb::GetHiggsWidthAtPoleMass

(

double

mass

)

Definition at line 371 of file TEvtProb.cc.

                                                   {
  if (mass>0.) return myCSW_.HiggsWidth(mass);
  else return -1.;
}

GetHXSWidthEstimator()

MELAHXSWidth const * TEvtProb::GetHXSWidthEstimator

(

)

const

Definition at line 355 of file TEvtProb.cc.

{ return &myCSW_; }

GetIORecord()

MelaIO * TEvtProb::GetIORecord

(

)

Definition at line 375 of file TEvtProb.cc.

{ return RcdME.getRef(); }

GetNCandidates()

int TEvtProb::GetNCandidates

(

)

Definition at line 384 of file TEvtProb.cc.

{ return (static_cast<int>(candList.size())); }

GetPrimaryHiggsMass()

double TEvtProb::GetPrimaryHiggsMass

(

)

Definition at line 362 of file TEvtProb.cc.

{ return PrimaryHMass; }

GetPrimaryMass()

double TEvtProb::GetPrimaryMass

(

int

ipart

)

Definition at line 363 of file TEvtProb.cc.

                                        {
  if (PDGHelpers::isAHiggs(ipart)) return GetPrimaryHiggsMass();
  else return TUtil::GetMass(ipart);
}

GetPrimaryWidth()

double TEvtProb::GetPrimaryWidth

(

int

ipart

)

Definition at line 367 of file TEvtProb.cc.

                                         {
  if (PDGHelpers::isAHiggs(ipart)) return myCSW_.HiggsWidth(GetPrimaryHiggsMass());
  else return TUtil::GetDecayWidth(ipart);
}

GetSelfDaTQGCCouplings()

aTQGCCouplings * TEvtProb::GetSelfDaTQGCCouplings

(

)

Definition at line 360 of file TEvtProb.cc.

{ return selfDaTQGCCoupl.getRef(); }

GetSelfDAZffCouplings()

AZffCouplings * TEvtProb::GetSelfDAZffCouplings

(

)

Definition at line 361 of file TEvtProb.cc.

{ return selfDAZffCoupl.getRef(); }

GetSelfDSpinOneCouplings()

SpinOneCouplings * TEvtProb::GetSelfDSpinOneCouplings

(

)

Definition at line 357 of file TEvtProb.cc.

{ return selfDSpinOneCoupl.getRef(); }

GetSelfDSpinTwoCouplings()

SpinTwoCouplings * TEvtProb::GetSelfDSpinTwoCouplings

(

)

Definition at line 358 of file TEvtProb.cc.

{ return selfDSpinTwoCoupl.getRef(); }

GetSelfDSpinZeroCouplings()

SpinZeroCouplings * TEvtProb::GetSelfDSpinZeroCouplings

(

)

Definition at line 356 of file TEvtProb.cc.

{ return selfDSpinZeroCoupl.getRef(); }

GetSelfDVprimeCouplings()

VprimeCouplings * TEvtProb::GetSelfDVprimeCouplings

(

)

Definition at line 359 of file TEvtProb.cc.

{ return selfDVprimeCoupl.getRef(); }

GetTopCandidates()

std::vector< MELATopCandidate_t * > * TEvtProb::GetTopCandidates

(

)

Definition at line 385 of file TEvtProb.cc.

{ return &topCandList; }

GetXPropagator()

double TEvtProb::GetXPropagator

(

TVar::ResonancePropagatorScheme

scheme

)

Definition at line 1011 of file TEvtProb.cc.

   {
  if (!CheckInputPresent()) return 0.;
  return TUtil::ResonancePropagator(melaCand->m(), scheme);
}

InitializeJHUGen()

void TEvtProb::InitializeJHUGen ( const char *  pathtoPDFSet, int  PDFMember  )

protected

Definition at line 142 of file TEvtProb.cc.

                                                                      {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb::InitializeJHUGen" << endl;
 
  InitJHUGenMELA(pathtoPDFSet, PDFMember, 2.*EBEAM);
 
  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb::InitializeJHUGen" << endl;
}

InitializeMCFM()

void TEvtProb::InitializeMCFM ( )

protected

Definition at line 82 of file TEvtProb.cc.

                             {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb::InitializeMCFM" << endl;
 
  mcfm_init_((char *)"input.DAT", (char *)"./");
  if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::TEvtProb: mcfm_init successful" << endl;
  energy_.sqrts = 2.*EBEAM;
  ResetQuarkMasses();
  ResetMCFM_EWKParameters(
    1.16639E-05, // Gf
    1./128., // alpha_EW
    80.399, // mW
    91.1876, // mZ
    0.23119, // xw=sin**2(thetaW)
    3 // MCFM EW scheme
    );
  qlinit_();
  // First resonance constant parameters
  spinzerohiggs_anomcoupl_.LambdaBSM=1000;
  spinzerohiggs_anomcoupl_.Lambda_Q=10000;
  spinzerohiggs_anomcoupl_.Lambda_z1=10000;
  spinzerohiggs_anomcoupl_.Lambda_z2=10000;
  spinzerohiggs_anomcoupl_.Lambda_z3=10000;
  spinzerohiggs_anomcoupl_.Lambda_z4=10000;
  spinzerohiggs_anomcoupl_.Lambda_z11=100; spinzerohiggs_anomcoupl_.Lambda_z21=100; spinzerohiggs_anomcoupl_.Lambda_z31=100; spinzerohiggs_anomcoupl_.Lambda_z41=100;
  spinzerohiggs_anomcoupl_.Lambda_z12=100; spinzerohiggs_anomcoupl_.Lambda_z22=100; spinzerohiggs_anomcoupl_.Lambda_z32=100; spinzerohiggs_anomcoupl_.Lambda_z42=100;
  spinzerohiggs_anomcoupl_.Lambda_z10=100; spinzerohiggs_anomcoupl_.Lambda_z20=100; spinzerohiggs_anomcoupl_.Lambda_z30=100; spinzerohiggs_anomcoupl_.Lambda_z40=100;
  spinzerohiggs_anomcoupl_.Lambda_zgs1=10000;
  spinzerohiggs_anomcoupl_.Lambda_w1=10000;
  spinzerohiggs_anomcoupl_.Lambda_w2=10000;
  spinzerohiggs_anomcoupl_.Lambda_w3=10000;
  spinzerohiggs_anomcoupl_.Lambda_w4=10000;
  spinzerohiggs_anomcoupl_.Lambda_w11=100; spinzerohiggs_anomcoupl_.Lambda_w21=100; spinzerohiggs_anomcoupl_.Lambda_w31=100; spinzerohiggs_anomcoupl_.Lambda_w41=100;
  spinzerohiggs_anomcoupl_.Lambda_w12=100; spinzerohiggs_anomcoupl_.Lambda_w22=100; spinzerohiggs_anomcoupl_.Lambda_w32=100; spinzerohiggs_anomcoupl_.Lambda_w42=100;
  spinzerohiggs_anomcoupl_.Lambda_w10=100; spinzerohiggs_anomcoupl_.Lambda_w20=100; spinzerohiggs_anomcoupl_.Lambda_w30=100; spinzerohiggs_anomcoupl_.Lambda_w40=100;
  // Second resonance constant parameters
  spinzerohiggs_anomcoupl_.Lambda2BSM=1000;
  spinzerohiggs_anomcoupl_.Lambda2_Q=10000;
  spinzerohiggs_anomcoupl_.Lambda2_z1=10000;
  spinzerohiggs_anomcoupl_.Lambda2_z2=10000;
  spinzerohiggs_anomcoupl_.Lambda2_z3=10000;
  spinzerohiggs_anomcoupl_.Lambda2_z4=10000;
  spinzerohiggs_anomcoupl_.Lambda2_z11=100; spinzerohiggs_anomcoupl_.Lambda2_z21=100; spinzerohiggs_anomcoupl_.Lambda2_z31=100; spinzerohiggs_anomcoupl_.Lambda2_z41=100;
  spinzerohiggs_anomcoupl_.Lambda2_z12=100; spinzerohiggs_anomcoupl_.Lambda2_z22=100; spinzerohiggs_anomcoupl_.Lambda2_z32=100; spinzerohiggs_anomcoupl_.Lambda2_z42=100;
  spinzerohiggs_anomcoupl_.Lambda2_z10=100; spinzerohiggs_anomcoupl_.Lambda2_z20=100; spinzerohiggs_anomcoupl_.Lambda2_z30=100; spinzerohiggs_anomcoupl_.Lambda2_z40=100;
  spinzerohiggs_anomcoupl_.Lambda2_zgs1=10000;
  spinzerohiggs_anomcoupl_.Lambda2_w1=10000;
  spinzerohiggs_anomcoupl_.Lambda2_w2=10000;
  spinzerohiggs_anomcoupl_.Lambda2_w3=10000;
  spinzerohiggs_anomcoupl_.Lambda2_w4=10000;
  spinzerohiggs_anomcoupl_.Lambda2_w11=100; spinzerohiggs_anomcoupl_.Lambda2_w21=100; spinzerohiggs_anomcoupl_.Lambda2_w31=100; spinzerohiggs_anomcoupl_.Lambda2_w41=100;
  spinzerohiggs_anomcoupl_.Lambda2_w12=100; spinzerohiggs_anomcoupl_.Lambda2_w22=100; spinzerohiggs_anomcoupl_.Lambda2_w32=100; spinzerohiggs_anomcoupl_.Lambda2_w42=100;
  spinzerohiggs_anomcoupl_.Lambda2_w10=100; spinzerohiggs_anomcoupl_.Lambda2_w20=100; spinzerohiggs_anomcoupl_.Lambda2_w30=100; spinzerohiggs_anomcoupl_.Lambda2_w40=100;
  // Switches for spin-0 Higgs couplings
  spinzerohiggs_anomcoupl_.channeltoggle_stu=2; spinzerohiggs_anomcoupl_.vvhvvtoggle_vbfvh=2;
  spinzerohiggs_anomcoupl_.AnomalCouplPR=1; spinzerohiggs_anomcoupl_.AnomalCouplDK=1;
  // Constant parameters for all processes
  qlfirst_.qlfirst=false;
 
  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb::InitializeMCFM" << endl;
}

ResetCouplings()

void TEvtProb::ResetCouplings

(

)

Definition at line 328 of file TEvtProb.cc.

                             {
  selfDSpinZeroCoupl.reset();
  selfDSpinOneCoupl.reset();
  selfDSpinTwoCoupl.reset();
  selfDVprimeCoupl.reset();
  selfDaTQGCCoupl.reset();
  selfDAZffCoupl.reset();
  AllowSeparateWWCouplings(false);
  ResetAmplitudeIncludes();
}

ResetInputEvent()

void TEvtProb::ResetInputEvent

(

)

Definition at line 338 of file TEvtProb.cc.

                              {
  RcdME.melaCand = 0;
  melaCand = 0;
 
  // Clear bookkeeping objects
  // Clean MELACandidates first since they contain all other objects
  for (MELACandidate* tmpPart:candList) delete tmpPart;
  candList.clear();
  // Clean MELATopCandidates next since they contain MELAParticles
  for (MELATopCandidate_t* tmpPart:topCandList) delete tmpPart;
  topCandList.clear();
  // Clean all remaining MELAParticles
  for (MELAParticle* tmpPart:particleList) delete tmpPart;
  particleList.clear();
}

ResetIORecord()

void TEvtProb::ResetIORecord

(

)

Definition at line 310 of file TEvtProb.cc.

{ RcdME.reset(); }

ResetMass()

void TEvtProb::ResetMass	(	double	inmass,
		int	ipart
	)

Definition at line 312 of file TEvtProb.cc.

{ TUtil::SetMass(inmass, ipart); }

ResetMCFM_EWKParameters()

void TEvtProb::ResetMCFM_EWKParameters	(	double	ext_Gf,
		double	ext_aemmz,
		double	ext_mW,
		double	ext_mZ,
		double	ext_xW,
		int	ext_ewscheme = 3
	)

Definition at line 327 of file TEvtProb.cc.

{ SetEwkCouplingParameters(ext_Gf, ext_aemmz, ext_mW, ext_mZ, ext_xW, ext_ewscheme); }

ResetQuarkMasses()

void TEvtProb::ResetQuarkMasses

(

)

Definition at line 316 of file TEvtProb.cc.

                               {
  ResetMass(1e-3, 1); // d
  ResetMass(5e-3, 2); // u
  ResetMass(1e-1, 3); // s
  ResetMass(1.275, 4); // c
  ResetMass(4.75, 5); // b
  ResetMass(173.2, 6); // t
  ResetMass(1e5, 7); // bprime
  ResetMass(1e5, 8); // tprime
  //coupling_(); // Already called at b and t twice!
}

ResetRenFacScaleMode()

void TEvtProb::ResetRenFacScaleMode

(

)

Definition at line 311 of file TEvtProb.cc.

{ SetRenFacScaleMode(TVar::DefaultScaleScheme, TVar::DefaultScaleScheme, 0.5, 0.5); }

ResetWidth()

void TEvtProb::ResetWidth	(	double	inwidth,
		int	ipart
	)

Definition at line 315 of file TEvtProb.cc.

{ TUtil::SetDecayWidth(inwidth, ipart); }

Set_LHAgrid()

void TEvtProb::Set_LHAgrid	(	const char *	path,
		int	pdfmember = 0
	)

Definition at line 185 of file TEvtProb.cc.

                                                         {
  char path_nnpdf_c[200];
  sprintf(path_nnpdf_c, "%s", path);
  int pathLength = strlen(path_nnpdf_c);
  nnpdfdriver_(path_nnpdf_c, &pathLength);
  nninitpdf_(&pdfmember);
}

SetCandidateDecayMode()

void TEvtProb::SetCandidateDecayMode

(

TVar::CandidateDecayMode

mode

)

Definition at line 208 of file TEvtProb.cc.

{ PDGHelpers::setCandidateDecayMode(mode); }

SetCurrentCandidate()

void TEvtProb::SetCurrentCandidate

(

MELACandidate *

cand

)

Definition at line 303 of file TEvtProb.cc.

                                                     {
  melaCand = cand;
  if (verbosity>=TVar::INFO && !melaCand) MELAout << "TEvtProb::SetCurrentCandidate: BE CAREFUL! !melaCand!" << endl;
  if (verbosity>=TVar::INFO && GetCurrentCandidateIndex()<0) MELAout << "TEvtProb::SetCurrentCandidate: The current candidate is not in the list of candidates. It is the users' responsibility to delete this candidate and all of its associated particles." << endl;
}

SetCurrentCandidateFromIndex()

void TEvtProb::SetCurrentCandidateFromIndex

(

unsigned int

icand

)

Definition at line 299 of file TEvtProb.cc.

                                                             {
  if (candList.size()>icand) melaCand = candList.at(icand);
  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::SetCurrentCandidateFromIndex: icand=" << icand << ">=candList.size()=" << candList.size() << endl;
}

SetHiggsMass()

void TEvtProb::SetHiggsMass	(	double	mass,
		double	wHiggs = -1.,
		int	whichResonance = -1
	)

Definition at line 217 of file TEvtProb.cc.

                                                                         {
  // Regular, first resonance
  if (whichResonance==1 || whichResonance==-1){
    if (mass<0.){
      _hmass = -1;
      _hwidth = 0;
    }
    else if (wHiggs<0.){
      _hmass = mass;
      _hwidth = myCSW_.HiggsWidth(_hmass);
    }
    else{
      _hmass = mass;
      _hwidth = wHiggs;
    }
    masses_mcfm_.hmass = _hmass;
    masses_mcfm_.hwidth = _hwidth;
 
    if (_hmass<0.) SetJHUGenHiggsMassWidth(0., _hwidth);
    else SetJHUGenHiggsMassWidth(_hmass, _hwidth);
  }
 
  // Second resonance
  if (whichResonance==2){
    if (mass<0.){
      _h2mass = -1;
      _h2width = 0;
    }
    else if (wHiggs<0.){
      _h2mass = mass;
      _h2width = myCSW_.HiggsWidth(_h2mass);
    }
    else{
      _h2mass = mass;
      _h2width = wHiggs;
    }
    spinzerohiggs_anomcoupl_.h2mass = _h2mass;
    spinzerohiggs_anomcoupl_.h2width = _h2width;
    //SetJHUGenHiggsMassWidth(_h2mass, _h2width); // Second resonance is not implemented in JHUGen yet.
  }
  else if (whichResonance==-1){
    _h2mass = -1;
    _h2width = 0;
    spinzerohiggs_anomcoupl_.h2mass = _h2mass;
    spinzerohiggs_anomcoupl_.h2width = _h2width;
    //SetJHUGenHiggsMassWidth(_h2mass, _h2width); // Second resonance is not implemented in JHUGen yet.
  }
 
  if (verbosity>=TVar::DEBUG) MELAout
    << "TEvtProb::SetHiggsMass(" << mass << ", " << wHiggs << ", " << whichResonance << "):\n"
    << '\t' << "hmass: " << _hmass << ", " << _hwidth << '\n'
    << '\t' << "h2mass: " << _h2mass << ", " << _h2width << '\n'
    << '\t' << "MCFM hmass: " << masses_mcfm_.hmass << ", " << masses_mcfm_.hwidth << '\n'
    << '\t' << "MCFM h2mass: " << spinzerohiggs_anomcoupl_.h2mass << ", " << spinzerohiggs_anomcoupl_.h2width << endl;
}

SetInputEvent()

void TEvtProb::SetInputEvent	(	SimpleParticleCollection_t *	pDaughters,
		SimpleParticleCollection_t *	pAssociated = 0,
		SimpleParticleCollection_t *	pMothers = 0,
		bool	isGen = false
	)

Definition at line 272 of file TEvtProb.cc.

   {
  MELACandidate* cand = ConvertVectorFormat(
    pDaughters,
    pAssociated,
    pMothers,
    isGen,
    &particleList, &candList // push_back is done automatically
    );
  if (cand) melaCand=cand;
}

SetLeptonInterf()

void TEvtProb::SetLeptonInterf

(

TVar::LeptonInterference

tmp

)

Definition at line 207 of file TEvtProb.cc.

{ leptonInterf = tmp; }

SetPrimaryHiggsMass()

void TEvtProb::SetPrimaryHiggsMass

(

double

mass

)

Definition at line 216 of file TEvtProb.cc.

{ PrimaryHMass=mass; SetHiggsMass(PrimaryHMass, -1., -1); }

SetProcess()

void TEvtProb::SetProcess	(	TVar::Process	proc,
		TVar::MatrixElement	me,
		TVar::Production	prod
	)

Definition at line 192 of file TEvtProb.cc.

                                                                                    {
  matrixElement = me;
  production = prod;
  // In case s-channel processes are passed for JHUGen ME, flip them back to JHUGen-specific productions.
  if (matrixElement == TVar::JHUGen){
    if (production == TVar::Had_ZH_S) production = TVar::Had_ZH;
    else if (production == TVar::Had_WH_S) production = TVar::Had_WH;
    else if (production == TVar::Lep_ZH_S) production = TVar::Lep_ZH;
    else if (production == TVar::Lep_WH_S) production = TVar::Lep_WH;
    else if (production == TVar::JJVBF_S) production = TVar::JJVBF;
    else if (production == TVar::JJQCD_S) production = TVar::JJQCD;
  }
  process = proc;
}

SetRcdCandPtr()

void TEvtProb::SetRcdCandPtr ( )

protected

Definition at line 298 of file TEvtProb.cc.

{ RcdME.melaCand = melaCand; }

SetRenFacScaleMode()

void TEvtProb::SetRenFacScaleMode	(	TVar::EventScaleScheme	renormalizationSch,
		TVar::EventScaleScheme	factorizationSch,
		double	ren_sf,
		double	fac_sf
	)

Definition at line 209 of file TEvtProb.cc.

                                                                                                                                             {
  event_scales.renomalizationScheme = renormalizationSch;
  event_scales.factorizationScheme = factorizationSch;
  event_scales.ren_scale_factor = ren_sf;
  event_scales.fac_scale_factor = fac_sf;
}

SetVerbosity()

void TEvtProb::SetVerbosity

(

TVar::VerbosityLevel

tmp

)

Definition at line 206 of file TEvtProb.cc.

{ verbosity = tmp; }

SetWprimeMassWidth()

void TEvtProb::SetWprimeMassWidth	(	double	inmass,
		double	inwidth
	)

Definition at line 314 of file TEvtProb.cc.

{ this->ResetMass(inmass, 34); this->ResetWidth(inwidth, 34); }

SetZprimeMassWidth()

void TEvtProb::SetZprimeMassWidth	(	double	inmass,
		double	inwidth
	)

Definition at line 313 of file TEvtProb.cc.

{ this->ResetMass(inmass, 32); this->ResetWidth(inwidth, 32); }

XsecCalc_TTX()

double TEvtProb::XsecCalc_TTX	(	int	topProcess,
		int	topDecay
	)

Definition at line 975 of file TEvtProb.cc.

   {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_TTX" << endl;
  double dXsec = 0;
  ResetIORecord();
  if (!CheckInputPresent()) return dXsec;
 
  if (matrixElement == TVar::JHUGen){
    double Hqqcoupl[SIZE_HQQ][2]={ { 0 } };
    if (process == TVar::HSMHiggs) Hqqcoupl[gHIGGS_KAPPA][0] = 1.;
    else if (process == TVar::H0minus) Hqqcoupl[gHIGGS_KAPPA_TILDE][0] = 1.;
    else if (process == TVar::SelfDefine_spin0){
      for (int i=0; i<SIZE_HQQ; i++){
        for (int j=0; j<2; j++){
          if ((selfDSpinZeroCoupl.Httcoupl)[i][j]!=0. && production == TVar::ttH) Hqqcoupl[i][j] = (selfDSpinZeroCoupl.Httcoupl)[i][j];
          else if ((selfDSpinZeroCoupl.Hbbcoupl)[i][j]!=0. && production == TVar::bbH) Hqqcoupl[i][j] = (selfDSpinZeroCoupl.Hbbcoupl)[i][j];
          else Hqqcoupl[i][j] = (selfDSpinZeroCoupl.Hqqcoupl)[i][j];
        }
      }
    }
    SetJHUGenSpinZeroQQCouplings(Hqqcoupl);
 
    if (production == TVar::ttH) dXsec = TTHiggsMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, topDecay, topProcess, verbosity);
    else if (production == TVar::bbH) dXsec = BBHiggsMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, topProcess, verbosity);
    else if (verbosity >= TVar::ERROR) MELAout << "TEvtProb::XsecCalc_TTX only supports ttH and bbH productions for the moment." << endl;
    if (verbosity >= TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_TTX: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
  }
  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::XsecCalc_TTX: Non-JHUGen ttH MEs are not supported!" << endl;
 
  ResetCouplings();
  ResetRenFacScaleMode();
  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_TTX" << endl;
  return dXsec;
}

XsecCalc_VVXVV()

double TEvtProb::XsecCalc_VVXVV

(

)

Definition at line 729 of file TEvtProb.cc.

                               {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_VVXVV" << endl;
  double dXsec=0;
  ResetIORecord();
  if (!CheckInputPresent()) return dXsec;
 
  bool useMCFM = matrixElement == TVar::MCFM;
  bool needBSMHiggs=false;
  bool needATQGC=false;
  bool calculateME=false;
  if (useMCFM){
    if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_VVXVV: Try MCFM" << endl;
    needBSMHiggs = CheckSelfDCouplings_HVV();
    needATQGC = CheckSelfDCouplings_aTQGC();
    if (needBSMHiggs || needATQGC) SetLeptonInterf(TVar::InterfOn); // All anomalous coupling computations have to use lepton interference
 
    calculateME = (
      production == TVar::Had_WH || production == TVar::Had_ZH || production == TVar::Lep_WH || production == TVar::Lep_ZH || production == TVar::JJVBF || production == TVar::JJEW || production == TVar::JJQCD || production == TVar::JJEWQCD
      || production == TVar::Had_WH_S || production == TVar::Had_ZH_S || production == TVar::Lep_WH_S || production == TVar::Lep_ZH_S || production == TVar::JJVBF_S || production == TVar::JJEW_S/* || production == TVar::JJQCD_S*/ || production == TVar::JJEWQCD_S
      || production == TVar::Had_WH_TU || production == TVar::Had_ZH_TU || production == TVar::Lep_WH_TU || production == TVar::Lep_ZH_TU || production == TVar::JJVBF_TU || production == TVar::JJEW_TU/* || production == TVar::JJQCD_TU*/ || production == TVar::JJEWQCD_TU
      );
    if (calculateME){
      SetMCFMSpinZeroCouplings(needBSMHiggs, &selfDSpinZeroCoupl, false);
      SetMCFMaTQGCCouplings(needATQGC, &selfDaTQGCCoupl);
      dXsec = SumMatrixElementPDF(process, production, matrixElement, leptonInterf, &event_scales, &RcdME, EBEAM, verbosity);
    }
    else if (verbosity>=TVar::INFO) MELAout << "TEvtProb::XsecCalc_VVXVV: MCFM_chooser failed to determine the process configuration." << endl;
  }
  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::XsecCalc_VVXVV: Non-MCFM Mes are not supported." << endl;
 
  if (verbosity >= TVar::DEBUG) MELAout
    << "Process " << TVar::ProcessName(process)
    << " TEvtProb::XsecCalc_VVXVV: dXsec=" << dXsec
    << endl;
 
  ResetCouplings(); // Should come first
  if (useMCFM){ // Set defaults. Should come next...
    if (needBSMHiggs || needATQGC) SetLeptonInterf(TVar::DefaultLeptonInterf);
    if (calculateME){
      // ... because of these!
      SetMCFMSpinZeroCouplings(false, &selfDSpinZeroCoupl, true);
      SetMCFMaTQGCCouplings(false, &selfDaTQGCCoupl);
    }
  }
  ResetRenFacScaleMode();
  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_VVXVV" << endl;
  return dXsec;
}

XsecCalc_VX()

double TEvtProb::XsecCalc_VX

(

bool

includeHiggsDecay

)

Definition at line 899 of file TEvtProb.cc.

   {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_VX" << endl;
  if (matrixElement == TVar::MCFM) return XsecCalc_VVXVV();
  double dXsec = 0;
  ResetIORecord();
  if (!CheckInputPresent()) return dXsec;
 
  AllowSeparateWWCouplings(false); // HZZ couplings are used for both in spin-0
 
  if (matrixElement == TVar::JHUGen){
    double Hvvcoupl[SIZE_HVV][2] ={ { 0 } };
    double HvvLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ] ={ { 0 } };
    int HvvCLambda_qsq[SIZE_HVV_CQSQ] ={ 0 };
    double Hvvpcoupl[SIZE_HVV][2] = { { 0 } };
    double Hvpvpcoupl[SIZE_HVV][2] = { { 0 } };
    double Zpffcoupl[SIZE_Vpff][2] = { { 0 } };
    double Wpffcoupl[SIZE_Vpff][2] = { { 0 } };
    double M_Zprime = -1;
    double Ga_Zprime = 0;
    double M_Wprime = -1;
    double Ga_Wprime = 0;
 
    for (int ic=0; ic<SIZE_HVV_LAMBDAQSQ; ic++){ for (int ik=0; ik<SIZE_HVV_CQSQ; ik++) HvvLambda_qsq[ic][ik]=100.; }
 
    if (process == TVar::HSMHiggs) Hvvcoupl[gHIGGS_VV_1][0] = 1.;
    else if (process == TVar::H0_g1prime2) Hvvcoupl[gHIGGS_VV_1_PRIME2][0] = 1.;
    else if (process == TVar::H0hplus) Hvvcoupl[gHIGGS_VV_2][0] = 1.;
    else if (process == TVar::H0minus) Hvvcoupl[gHIGGS_VV_4][0] = 1.;
    else if (process == TVar::H0_Zgsg1prime2) Hvvcoupl[gHIGGS_ZA_1_PRIME2][0] = 1.;
    else if (process == TVar::H0_Zgs) Hvvcoupl[gHIGGS_ZA_2][0] = 1.;
    else if (process == TVar::H0_Zgs_PS) Hvvcoupl[gHIGGS_ZA_4][0] = 1.;
    else if (process == TVar::H0_gsgs) Hvvcoupl[gHIGGS_AA_2][0] = 1.;
    else if (process == TVar::H0_gsgs_PS) Hvvcoupl[gHIGGS_AA_4][0] = 1.;
    else if (process == TVar::SelfDefine_spin0){
      for (int i=0; i<SIZE_HVV; i++){
        for (int j=0; j<2; j++){
          Hvvcoupl[i][j] = (selfDSpinZeroCoupl.Hzzcoupl)[i][j];
 
          Hvvpcoupl[i][j] = (selfDSpinZeroCoupl.Hzzpcoupl)[i][j];
          Hvpvpcoupl[i][j] = (selfDSpinZeroCoupl.Hzpzpcoupl)[i][j];
        }
      }
      for (int i=0; i<SIZE_Vpff; i++){
        for (int j=0; j<2; j++){
          Zpffcoupl[i][j] = (selfDVprimeCoupl.Zpffcoupl)[i][j];
          Wpffcoupl[i][j] = (selfDVprimeCoupl.Wpffcoupl)[i][j];
        }
      }
      M_Zprime = selfDVprimeCoupl.M_Zprime;
      Ga_Zprime = selfDVprimeCoupl.Ga_Zprime;
      M_Wprime = selfDVprimeCoupl.M_Wprime;
      Ga_Wprime = selfDVprimeCoupl.Ga_Wprime;
      for (int j=0; j<SIZE_HVV_CQSQ; j++){
        for (int i=0; i<SIZE_HVV_LAMBDAQSQ; i++) HvvLambda_qsq[i][j] = (selfDSpinZeroCoupl.HzzLambda_qsq)[i][j];
        HvvCLambda_qsq[j] = (selfDSpinZeroCoupl.HzzCLambda_qsq)[j];
      }
    }
    SetJHUGenSpinZeroVVCouplings(Hvvcoupl, Hvvpcoupl, Hvpvpcoupl, HvvCLambda_qsq, HvvLambda_qsq, false);
    SetJHUGenVprimeContactCouplings(Zpffcoupl, Wpffcoupl);
    SetZprimeMassWidth(M_Zprime, Ga_Zprime);
    SetWprimeMassWidth(M_Wprime, Ga_Wprime);
 
    dXsec = VHiggsMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, includeHiggsDecay, verbosity);
    if (verbosity >= TVar::DEBUG) MELAout << "TEVtProb::XsecCalc_VX: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
  } // end of JHUGen matrix element calculations
  else if (verbosity>=TVar::ERROR) MELAerr << "TEVtProb::XsecCalc_VX: Non-JHUGen VH MEs are not supported!" << endl;
 
  ResetCouplings();
  ResetRenFacScaleMode();
  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_VX" << endl;
  return dXsec;
}

XsecCalc_XVV()

double TEvtProb::XsecCalc_XVV

(

)

Definition at line 407 of file TEvtProb.cc.

                             {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_XVV" << endl;
  double dXsec=0;
  ResetIORecord();
  if (!CheckInputPresent()) return dXsec;
 
  bool useMCFM = matrixElement == TVar::MCFM;
  bool calculateME=false;
  bool needBSMHiggs=false;
  bool needAZff=CheckSelfDCouplings_AZff();
  if (useMCFM){
    if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XVV: Try MCFM" << endl;
    needBSMHiggs = CheckSelfDCouplings_Hgg() || CheckSelfDCouplings_Htt() || CheckSelfDCouplings_Hbb() || CheckSelfDCouplings_HVV();
    if (needBSMHiggs || needAZff) SetLeptonInterf(TVar::InterfOn); // All anomalous coupling computations have to use lepton interference
 
    calculateME = (
      production == TVar::ZZGG ||
      production == TVar::ZZQQB ||
      production == TVar::ZZQQB_STU ||
      production == TVar::ZZQQB_S ||
      production == TVar::ZZQQB_TU ||
      production == TVar::ZZINDEPENDENT ||
      (production == TVar::JJQCD && process == TVar::bkgZJets) // Allow Z+jets to be computed here.
      );
    if (calculateME){
      SetMCFMSpinZeroCouplings(needBSMHiggs, &selfDSpinZeroCoupl, false);
      SetMCFMAZffCouplings(needAZff, &selfDAZffCoupl);
      dXsec = SumMatrixElementPDF(process, production, matrixElement, leptonInterf, &event_scales, &RcdME, EBEAM, verbosity);
    }
    else if (verbosity>=TVar::INFO) MELAout << "TEvtProb::XsecCalc_XVV: MCFM_chooser failed to determine the process configuration." << endl;
  }
  else if (matrixElement == TVar::JHUGen){
    if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XVV::Try JHUGen" << endl;
    AllowSeparateWWCouplings(false); // HZZ couplings are used for both in spin-0
    // all the possible couplings
    double Hggcoupl[SIZE_HGG][2] ={ { 0 } };
    double Hvvcoupl[SIZE_HVV][2] ={ { 0 } };
    double HvvLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ] ={ { 0 } };
    int HvvCLambda_qsq[SIZE_HVV_CQSQ] ={ 0 };
 
    double Hvvpcoupl[SIZE_HVV][2] = { { 0 } };
    double Hvpvpcoupl[SIZE_HVV][2] = { { 0 } };
 
    double Zqqcoupl[SIZE_ZQQ][2] ={ { 0 } };
    double Zvvcoupl[SIZE_ZVV][2] ={ { 0 } };
 
    double Gqqcoupl[SIZE_GQQ][2] ={ { 0 } };
    double Gggcoupl[SIZE_GGG][2] ={ { 0 } };
    double Gvvcoupl[SIZE_GVV][2] ={ { 0 } };
    double Gvvpcoupl[SIZE_GVV][2] ={ { 0 } };
    double Gvpvpcoupl[SIZE_GVV][2] ={ { 0 } };
 
    double Zpffcoupl[SIZE_Vpff][2] ={ { 0 } };
    double Wpffcoupl[SIZE_Vpff][2] ={ { 0 } };
    double M_Zprime = -1;
    double Ga_Zprime = 0;
    double M_Wprime = -1;
    double Ga_Wprime = 0;
 
    double AZffcoupl[SIZE_AZff][2] ={ { 0 } };
    for (int j=0; j<2; j++){
      for (int i=0; i<SIZE_AZff; i++){
        AZffcoupl[i][j] = (selfDAZffCoupl.AZffcoupl)[i][j];
      }
    }
    SetJHUGenAZffCouplings(needAZff, AZffcoupl);
 
    //
    // set spin 0 default numbers
    //
    // By default set the Spin 0 couplings for SM case (0+m)
    Hggcoupl[gHIGGS_GG_2][0]=1.0;
    for (int ic=0; ic<SIZE_HVV_LAMBDAQSQ; ic++){ for (int ik=0; ik<SIZE_HVV_CQSQ; ik++) HvvLambda_qsq[ic][ik]=100.; }
    //
    // set spin 1 default numbers
    //
    Zqqcoupl[gZPRIME_QQ_LEFT][0]=1.0;
    Zqqcoupl[gZPRIME_QQ_RIGHT][0]=1.0;
    //
    // set spin 2 default numbers
    //
    Gqqcoupl[gGRAVITON_QQ_LEFT][0]=1.0;
    Gqqcoupl[gGRAVITON_QQ_RIGHT][0]=1.0;
 
    bool isSpinZero = false;
    bool isSpinOne = false;
    bool isSpinTwo = false;
 
    // 0m+
    if (process == TVar::HSMHiggs){
      Hvvcoupl[gHIGGS_VV_1][0]=1.;
      isSpinZero = true;
    }
    // 0+L1
    else if (process == TVar::H0_g1prime2){
      Hvvcoupl[gHIGGS_VV_1_PRIME2][0] = 1.;
      isSpinZero = true;
    }
    // 0h+
    else if (process == TVar::H0hplus){
      Hvvcoupl[gHIGGS_VV_2][0] = 1.;
      isSpinZero = true;
    }
    // 0-
    else if (process == TVar::H0minus){
      Hvvcoupl[gHIGGS_VV_4][0] = 1.;
      isSpinZero = true;
    }
    else if (process == TVar::H0_Zgsg1prime2){
      Hvvcoupl[gHIGGS_ZA_1_PRIME2][0] = 1.;
      isSpinZero = true;
    }
    else if (process == TVar::H0_Zgs){
      Hvvcoupl[gHIGGS_ZA_2][0] = 1.;
      isSpinZero = true;
    }
    else if (process == TVar::H0_Zgs_PS){
      Hvvcoupl[gHIGGS_ZA_4][0] = 1.;
      isSpinZero = true;
    }
    else if (process == TVar::H0_gsgs){
      Hvvcoupl[gHIGGS_AA_2][0] = 1.;
      isSpinZero = true;
    }
    else if (process == TVar::H0_gsgs_PS){
      Hvvcoupl[gHIGGS_AA_4][0] = 1.;
      isSpinZero = true;
    }
    else if (process == TVar::SelfDefine_spin0){
      for (int j=0; j<2; j++){
        for (int i=0; i<SIZE_HGG; i++) Hggcoupl[i][j] = (selfDSpinZeroCoupl.Hggcoupl)[i][j];
        for (int i=0; i<SIZE_HVV; i++){
          Hvvcoupl[i][j] = (selfDSpinZeroCoupl.Hzzcoupl)[i][j];
 
          Hvvpcoupl[i][j] = (selfDSpinZeroCoupl.Hzzpcoupl)[i][j];
          Hvpvpcoupl[i][j] = (selfDSpinZeroCoupl.Hzpzpcoupl)[i][j];
        }
      }
      for (int j=0; j<SIZE_HVV_CQSQ; j++){
        for (int i=0; i<SIZE_HVV_LAMBDAQSQ; i++) HvvLambda_qsq[i][j] = (selfDSpinZeroCoupl.HzzLambda_qsq)[i][j];
        HvvCLambda_qsq[j] = (selfDSpinZeroCoupl.HzzCLambda_qsq)[j];
      }
 
      if (verbosity>=TVar::DEBUG_VERBOSE){
        for (int j=0; j<2; j++){
          for (int i=0; i<SIZE_HGG; i++){ if ((selfDSpinZeroCoupl.Hggcoupl)[i][j]!=0.) MELAout << "Hggcoupl[" << i << "][" << j << "] = " << (selfDSpinZeroCoupl.Hggcoupl)[i][j] << endl; }
          for (int i=0; i<SIZE_HVV; i++){
            if ((selfDSpinZeroCoupl.Hzzcoupl)[i][j]!=0.) MELAout << "Hvvcoupl[" << i << "][" << j << "] = " << (selfDSpinZeroCoupl.Hzzcoupl)[i][j] << endl;
 
            if ((selfDSpinZeroCoupl.Hzzpcoupl)[i][j]!=0.) MELAout << "Hzzpcoupl[" << i << "][" << j << "] = " << (selfDSpinZeroCoupl.Hzzpcoupl)[i][j] << endl;
            if ((selfDSpinZeroCoupl.Hzpzpcoupl)[i][j]!=0.) MELAout << "Hzpzpcoupl[" << i << "][" << j << "] = " << (selfDSpinZeroCoupl.Hzpzpcoupl)[i][j] << endl;
          }
        }
        for (int j=0; j<SIZE_HVV_CQSQ; j++){
          for (int i=0; i<SIZE_HVV_LAMBDAQSQ; i++){ if ((selfDSpinZeroCoupl.HzzLambda_qsq)[i][j]!=0.) MELAout << "HvvLambda_qsq[" << i << "][" << j << "] = " << (selfDSpinZeroCoupl.HzzLambda_qsq)[i][j] << endl; }
          if ((selfDSpinZeroCoupl.HzzCLambda_qsq)[j]!=0.) MELAout << "HvvCLambda_qsq[" << j << "] = " << (selfDSpinZeroCoupl.HzzCLambda_qsq)[j] << endl;
        }
      }
 
      isSpinZero = true;
    }
 
    // 1-
    else if (process == TVar::H1minus){
      Zvvcoupl[gZPRIME_VV_1][0]=1.;
      isSpinOne = true;
    }
    // 1+
    else if (process == TVar::H1plus){
      Zvvcoupl[gZPRIME_VV_2][0]=1.;
      isSpinOne = true;
    }
    else if (process == TVar::SelfDefine_spin1){
      for (int j=0; j<2; j++){
        for (int i=0; i<SIZE_ZQQ; i++) Zqqcoupl[i][j] = (selfDSpinOneCoupl.Zqqcoupl)[i][j];
        for (int i=0; i<SIZE_ZVV; i++) Zvvcoupl[i][j] = (selfDSpinOneCoupl.Zvvcoupl)[i][j];
      }
      isSpinOne = true;
    }
 
    // 2m+
    else if (process == TVar::H2_g1g5){
      Gggcoupl[gGRAVITON_GG_1][0]=1.;
      Gvvcoupl[gGRAVITON_VV_1][0]=1.;
      Gvvcoupl[gGRAVITON_VV_5][0]=1.;
      isSpinTwo = true;
    }
    // 2b+
    else if (process == TVar::H2_g1){
      Gggcoupl[gGRAVITON_GG_1][0]=1.;
      Gvvcoupl[gGRAVITON_VV_1][0]=1.;
      isSpinTwo = true;
    }
    else if (process == TVar::H2_g2){
      Gggcoupl[gGRAVITON_GG_2][0]=1.;
      Gvvcoupl[gGRAVITON_VV_2][0]=1.;
      isSpinTwo = true;
    }
    // 2h3plus
    else if (process == TVar::H2_g3){
      Gggcoupl[gGRAVITON_GG_3][0]=1.;
      Gvvcoupl[gGRAVITON_VV_3][0]=1.;
      isSpinTwo = true;
    }
    // 2h+
    else if (process == TVar::H2_g4){
      Gggcoupl[gGRAVITON_GG_4][0]=1.;
      Gvvcoupl[gGRAVITON_VV_4][0]=1.;
      isSpinTwo = true;
    }
    // 2b+
    else if (process == TVar::H2_g5){
      Gggcoupl[gGRAVITON_GG_1][0]=1.;
      Gvvcoupl[gGRAVITON_VV_5][0]=1.;
      isSpinTwo = true;
    }
    // 2h6+
    else if (process == TVar::H2_g6){
      Gggcoupl[gGRAVITON_GG_1][0]=1.;
      Gvvcoupl[gGRAVITON_VV_6][0]=1.;
      isSpinTwo = true;
    }
    // 2h7plus
    else if (process == TVar::H2_g7){
      Gggcoupl[gGRAVITON_GG_1][0]=1.;
      Gvvcoupl[gGRAVITON_VV_7][0]=1.;
      isSpinTwo = true;
    }
    // 2h-
    else if (process == TVar::H2_g8){
      Gggcoupl[gGRAVITON_GG_5][0]=1.;
      Gvvcoupl[gGRAVITON_VV_8][0]=1.;
      isSpinTwo = true;
    }
    // 2h9minus
    else if (process == TVar::H2_g9){
      Gggcoupl[gGRAVITON_GG_5][0]=1.;
      Gvvcoupl[gGRAVITON_VV_9][0]=1.;
      isSpinTwo = true;
    }
    // 2h10minus
    else if (process == TVar::H2_g10){
      Gggcoupl[gGRAVITON_GG_5][0]=1.;
      Gvvcoupl[gGRAVITON_VV_10][0]=1.;
      isSpinTwo = true;
    }
    else if (process == TVar::SelfDefine_spin2){
      for (int j=0; j<2; j++){
        for (int i=0; i<SIZE_GGG; i++) Gggcoupl[i][j] = (selfDSpinTwoCoupl.Gggcoupl)[i][j];
        for (int i=0; i<SIZE_GQQ; i++) Gqqcoupl[i][j] = (selfDSpinTwoCoupl.Gqqcoupl)[i][j];
        for (int i=0; i<SIZE_GVV; i++){
          Gvvcoupl[i][j] = (selfDSpinTwoCoupl.Gvvcoupl)[i][j];
          Gvvpcoupl[i][j] = (selfDSpinTwoCoupl.Gvvpcoupl)[i][j];
          Gvpvpcoupl[i][j] = (selfDSpinTwoCoupl.Gvpvpcoupl)[i][j];
        }
      }
      isSpinTwo = true;
    }
 
    // Vprime / contact couplings
    if (process == TVar::SelfDefine_spin0 || process == TVar::SelfDefine_spin2){
      for (int j=0; j<2; j++){
        for (int i=0; i<SIZE_Vpff; i++){
          Zpffcoupl[i][j] = (selfDVprimeCoupl.Zpffcoupl)[i][j];
          Wpffcoupl[i][j] = (selfDVprimeCoupl.Wpffcoupl)[i][j];
        }
      }
      M_Zprime = selfDVprimeCoupl.M_Zprime;
      Ga_Zprime = selfDVprimeCoupl.Ga_Zprime;
      M_Wprime = selfDVprimeCoupl.M_Wprime;
      Ga_Wprime = selfDVprimeCoupl.Ga_Wprime;
 
      if (verbosity>=TVar::DEBUG_VERBOSE){
        for (int j=0; j<2; j++){
          for (int i=0; i<SIZE_Vpff; i++){
            if ((selfDVprimeCoupl.Zpffcoupl)[i][j]!=0.) MELAout << "Zpffcoupl[" << i << "][" << j << "] = " << (selfDVprimeCoupl.Zpffcoupl)[i][j] << endl;
            if ((selfDVprimeCoupl.Wpffcoupl)[i][j]!=0.) MELAout << "Wpffcoupl[" << i << "][" << j << "] = " << (selfDVprimeCoupl.Wpffcoupl)[i][j] << endl;
          }
        }
        MELAout << "M_Zprime = " << selfDVprimeCoupl.M_Zprime << endl;
        MELAout << "Ga_Zprime = " << selfDVprimeCoupl.Ga_Zprime << endl;
        MELAout << "M_Wprime = " << selfDVprimeCoupl.M_Wprime << endl;
        MELAout << "Ga_Wprime = " << selfDVprimeCoupl.Ga_Wprime << endl;
      }
    }
 
    if (isSpinZero){
      SetJHUGenSpinZeroGGCouplings(Hggcoupl);
      SetJHUGenSpinZeroVVCouplings(Hvvcoupl, Hvvpcoupl, Hvpvpcoupl, HvvCLambda_qsq, HvvLambda_qsq, false);
      SetJHUGenVprimeContactCouplings(Zpffcoupl, Wpffcoupl);
      SetZprimeMassWidth(M_Zprime, Ga_Zprime);
      SetWprimeMassWidth(M_Wprime, Ga_Wprime);
    }
    else if (isSpinOne) SetJHUGenSpinOneCouplings(Zqqcoupl, Zvvcoupl);
    else if (isSpinTwo){
      SetJHUGenSpinTwoCouplings(Gggcoupl, Gvvcoupl, Gvvpcoupl, Gvpvpcoupl, Gqqcoupl);
      SetJHUGenVprimeContactCouplings(Zpffcoupl, Wpffcoupl);
      SetZprimeMassWidth(M_Zprime, Ga_Zprime);
      SetWprimeMassWidth(M_Wprime, Ga_Wprime);
    }
 
    if (isSpinZero || isSpinOne || isSpinTwo) dXsec = JHUGenMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, verbosity);
    else if (verbosity>=TVar::ERROR) MELAerr
      << "TEvtProb::XsecCalc_XVV: JHUGen ME is not spin zero, one or two! The process is described by "
      << "Process: " << TVar::ProcessName(process) << ", Production: " << TVar::ProductionName(production) << ", and ME: " << TVar::MatrixElementName(matrixElement)
      << endl;
  } // end of JHUGen matrix element calculations
 
  if (verbosity >= TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XVV: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
 
  if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XVV::Reset couplings" << endl;
  ResetCouplings(); // Should come first
  if (useMCFM){ // Set defaults. Should come next...
    if (needBSMHiggs) SetLeptonInterf(TVar::DefaultLeptonInterf);
    if (calculateME) SetMCFMSpinZeroCouplings(false, &selfDSpinZeroCoupl, true); // ... because of this!
  }
  ResetRenFacScaleMode();
  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_XVV" << endl;
  return dXsec;
}

XsecCalcXJ()

double TEvtProb::XsecCalcXJ

(

)

Definition at line 868 of file TEvtProb.cc.

                           {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_XJ" << endl;
  double dXsec = 0;
  ResetIORecord();
  if (!CheckInputPresent()) return dXsec;
 
  // Calculate the ME
  if (matrixElement == TVar::JHUGen){
    // Not currently supported, but implement anyway
    double Hggcoupl[SIZE_HGG][2] ={ { 0 } };
    if (process == TVar::HSMHiggs) Hggcoupl[gHIGGS_GG_2][0] = 1.;
    else if (process == TVar::H0minus) Hggcoupl[gHIGGS_GG_4][0] = 1.;
    else if (process == TVar::SelfDefine_spin0){
      for (int j=0; j<2; j++){
        for (int i=0; i<SIZE_HGG; i++) Hggcoupl[i][j] = (selfDSpinZeroCoupl.Hggcoupl)[i][j];
      }
    }
    SetJHUGenSpinZeroGGCouplings(Hggcoupl);
 
    dXsec = HJJMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, verbosity);
    if (verbosity >= TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XJ: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
  }
  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::XsecCalc_XJ: Non-JHUGen HJ MEs are not supported!" << endl;
 
  ResetCouplings();
  ResetRenFacScaleMode();
  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_XJ" << endl;
  return dXsec;
}

XsecCalcXJJ()

double TEvtProb::XsecCalcXJJ

(

)

Definition at line 779 of file TEvtProb.cc.

                            {
  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_XJJ" << endl;
  if (matrixElement == TVar::MCFM) return XsecCalc_VVXVV();
  double dXsec = 0;
  ResetIORecord();
  if (!CheckInputPresent()) return dXsec;
 
  if (matrixElement == TVar::JHUGen){
    if (production == TVar::JJQCD){
      double Hggcoupl[SIZE_HGG][2] ={ { 0 } };
      if (process == TVar::HSMHiggs) Hggcoupl[gHIGGS_GG_2][0] = 1.;
      else if (process == TVar::H0minus) Hggcoupl[gHIGGS_GG_4][0] = 1.;
      else if (process == TVar::SelfDefine_spin0){
        for (int j=0; j<2; j++){
          for (int i=0; i<SIZE_HGG; i++) Hggcoupl[i][j] = (selfDSpinZeroCoupl.Hggcoupl)[i][j];
        }
      }
      SetJHUGenSpinZeroGGCouplings(Hggcoupl);
    }
    else if (production == TVar::JJVBF){
      double Hzzcoupl[SIZE_HVV][2] ={ { 0 } };
      double Hwwcoupl[SIZE_HVV][2] ={ { 0 } };
      double HzzLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ] ={ { 0 } };
      double HwwLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ] ={ { 0 } };
      int HzzCLambda_qsq[SIZE_HVV_CQSQ] ={ 0 };
      int HwwCLambda_qsq[SIZE_HVV_CQSQ] ={ 0 };
      double Hzzpcoupl[SIZE_HVV][2] = { { 0 } };
      double Hzpzpcoupl[SIZE_HVV][2] = { { 0 } };
      double Zpffcoupl[SIZE_Vpff][2] = { { 0 } };
      double Hwwpcoupl[SIZE_HVV][2] = { { 0 } };
      double Hwpwpcoupl[SIZE_HVV][2] = { { 0 } };
      double Wpffcoupl[SIZE_Vpff][2] = { { 0 } };
      double M_Zprime = -1;
      double Ga_Zprime = 0;
      double M_Wprime = -1;
      double Ga_Wprime = 0;
 
      for (int ic=0; ic<SIZE_HVV_LAMBDAQSQ; ic++){ for (int ik=0; ik<SIZE_HVV_CQSQ; ik++){ HzzLambda_qsq[ic][ik]=100.; HwwLambda_qsq[ic][ik]=100.; } }
      SetJHUGenDistinguishWWCouplings(false);
 
      if (process == TVar::HSMHiggs){ Hzzcoupl[gHIGGS_VV_1][0] = 1.; Hwwcoupl[gHIGGS_VV_1][0] = 1.; }
      else if (process == TVar::H0_g1prime2){ Hzzcoupl[gHIGGS_VV_1_PRIME2][0] = 1.; Hwwcoupl[gHIGGS_VV_1_PRIME2][0] = 1.; }
      else if (process == TVar::H0hplus){ Hzzcoupl[gHIGGS_VV_2][0] = 1.; Hwwcoupl[gHIGGS_VV_2][0] = 1.; }
      else if (process == TVar::H0minus){ Hzzcoupl[gHIGGS_VV_4][0] = 1.; Hwwcoupl[gHIGGS_VV_4][0] = 1.; }
      else if (process == TVar::H0_Zgsg1prime2){ Hzzcoupl[gHIGGS_ZA_1_PRIME2][0] = 1.; }
      else if (process == TVar::H0_Zgs){ Hzzcoupl[gHIGGS_ZA_2][0] = 1.; }
      else if (process == TVar::H0_Zgs_PS){ Hzzcoupl[gHIGGS_ZA_4][0] = 1.; }
      else if (process == TVar::H0_gsgs){ Hzzcoupl[gHIGGS_AA_2][0] = 1.; }
      else if (process == TVar::H0_gsgs_PS){ Hzzcoupl[gHIGGS_AA_4][0] = 1.; }
      else if (process == TVar::SelfDefine_spin0){
        for (int j=0; j<2; j++){
          for (int i=0; i<SIZE_HVV; i++){
            Hzzcoupl[i][j] = (selfDSpinZeroCoupl.Hzzcoupl)[i][j]; Hwwcoupl[i][j] = (selfDSpinZeroCoupl.Hwwcoupl)[i][j];
            Hzzpcoupl[i][j] = (selfDSpinZeroCoupl.Hzzpcoupl)[i][j]; Hwwpcoupl[i][j] = (selfDSpinZeroCoupl.Hwwpcoupl)[i][j];
            Hzpzpcoupl[i][j] = (selfDSpinZeroCoupl.Hzpzpcoupl)[i][j]; Hwpwpcoupl[i][j] = (selfDSpinZeroCoupl.Hwpwpcoupl)[i][j];
          }
          for (int i=0; i<SIZE_Vpff; i++){
            Zpffcoupl[i][j] = (selfDVprimeCoupl.Zpffcoupl)[i][j]; Wpffcoupl[i][j] = (selfDVprimeCoupl.Wpffcoupl)[i][j];
          }
        }
        M_Zprime = selfDVprimeCoupl.M_Zprime;
        Ga_Zprime = selfDVprimeCoupl.Ga_Zprime;
        M_Wprime = selfDVprimeCoupl.M_Wprime;
        Ga_Wprime = selfDVprimeCoupl.Ga_Wprime;
        for (int j=0; j<SIZE_HVV_CQSQ; j++){
          for (int i=0; i<SIZE_HVV_LAMBDAQSQ; i++){ HzzLambda_qsq[i][j] = (selfDSpinZeroCoupl.HzzLambda_qsq)[i][j]; HwwLambda_qsq[i][j] = (selfDSpinZeroCoupl.HwwLambda_qsq)[i][j]; }
          HzzCLambda_qsq[j] = (selfDSpinZeroCoupl.HzzCLambda_qsq)[j]; HwwCLambda_qsq[j] = (selfDSpinZeroCoupl.HwwCLambda_qsq)[j];
        }
        SetJHUGenDistinguishWWCouplings(selfDSpinZeroCoupl.separateWWZZcouplings);
      }
      SetJHUGenSpinZeroVVCouplings(Hzzcoupl, Hzzpcoupl, Hzpzpcoupl, HzzCLambda_qsq, HzzLambda_qsq, false);
      SetJHUGenSpinZeroVVCouplings(Hwwcoupl, Hwwpcoupl, Hwpwpcoupl, HwwCLambda_qsq, HwwLambda_qsq, true); // Set the WW couplings regardless of SetJHUGenDistinguishWWCouplings(false/true) because of how JHUGen handles this true flag.
      SetJHUGenVprimeContactCouplings(Zpffcoupl, Wpffcoupl);
      SetZprimeMassWidth(M_Zprime, Ga_Zprime);
      SetWprimeMassWidth(M_Wprime, Ga_Wprime);
    }
 
    dXsec = HJJMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, verbosity);
    if (verbosity >= TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XJJ: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
  }
  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::XsecCalc_XJJ: Non-JHUGen vbfMELA MEs are not supported!" << endl;
 
  ResetCouplings();
  ResetRenFacScaleMode();
  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_XJJ" << endl;
  return dXsec;
}

Member Data Documentation

_h2mass

double TEvtProb::_h2mass

protected

Definition at line 149 of file TEvtProb.hh.

_h2width

double TEvtProb::_h2width

protected

Definition at line 150 of file TEvtProb.hh.

_hmass

double TEvtProb::_hmass

protected

Definition at line 147 of file TEvtProb.hh.

_hwidth

double TEvtProb::_hwidth

protected

Definition at line 148 of file TEvtProb.hh.

candList

std::vector<MELACandidate*> TEvtProb::candList

protected

Definition at line 165 of file TEvtProb.hh.

EBEAM

double TEvtProb::EBEAM

protected

Definition at line 151 of file TEvtProb.hh.

event_scales

TVar::event_scales_type TEvtProb::event_scales

protected

Definition at line 153 of file TEvtProb.hh.

leptonInterf

TVar::LeptonInterference TEvtProb::leptonInterf

protected

Definition at line 145 of file TEvtProb.hh.

matrixElement

TVar::MatrixElement TEvtProb::matrixElement

protected

Definition at line 142 of file TEvtProb.hh.

melaCand

MELACandidate* TEvtProb::melaCand

protected

Definition at line 163 of file TEvtProb.hh.

myCSW_

MELAHXSWidth TEvtProb::myCSW_

protected

Definition at line 152 of file TEvtProb.hh.

particleList

std::vector<MELAParticle*> TEvtProb::particleList

protected

Definition at line 164 of file TEvtProb.hh.

pathtoPDFSet_

const char* TEvtProb::pathtoPDFSet_

protected

Definition at line 139 of file TEvtProb.hh.

PDFMember_

int TEvtProb::PDFMember_

protected

Definition at line 140 of file TEvtProb.hh.

PrimaryHMass

double TEvtProb::PrimaryHMass

protected

Definition at line 146 of file TEvtProb.hh.

process

TVar::Process TEvtProb::process

protected

Definition at line 141 of file TEvtProb.hh.

production

TVar::Production TEvtProb::production

protected

Definition at line 143 of file TEvtProb.hh.

RcdME

MelaIO TEvtProb::RcdME

protected

Definition at line 161 of file TEvtProb.hh.

selfDaTQGCCoupl

aTQGCCouplings TEvtProb::selfDaTQGCCoupl

protected

Definition at line 159 of file TEvtProb.hh.

selfDAZffCoupl

AZffCouplings TEvtProb::selfDAZffCoupl

protected

Definition at line 160 of file TEvtProb.hh.

selfDSpinOneCoupl

SpinOneCouplings TEvtProb::selfDSpinOneCoupl

protected

Definition at line 156 of file TEvtProb.hh.

selfDSpinTwoCoupl

SpinTwoCouplings TEvtProb::selfDSpinTwoCoupl

protected

Definition at line 157 of file TEvtProb.hh.

selfDSpinZeroCoupl

SpinZeroCouplings TEvtProb::selfDSpinZeroCoupl

protected

Definition at line 155 of file TEvtProb.hh.

selfDVprimeCoupl

VprimeCouplings TEvtProb::selfDVprimeCoupl

protected

Definition at line 158 of file TEvtProb.hh.

topCandList

std::vector<MELATopCandidate_t*> TEvtProb::topCandList

protected

Definition at line 166 of file TEvtProb.hh.

verbosity

TVar::VerbosityLevel TEvtProb::verbosity

protected

Definition at line 144 of file TEvtProb.hh.

---

The documentation for this class was generated from the following files:

• MELA/interface/TEvtProb.hh
• MELA/src/TEvtProb.cc