ZZMatrixElement Class Reference

#include <ZZMatrixElement.h>

Collaboration diagram for ZZMatrixElement:

Public Member Functions
	ZZMatrixElement (const char *pathtoPDFSet, int PDFMember, const char *pathtoHiggsCSandWidth, double ebeam, TVar::VerbosityLevel verbosity)
	ZZMatrixElement (const ZZMatrixElement &other)
	~ZZMatrixElement ()
void	computeXS (float &mevalue)
void	computeProdXS_VVHVV (float &mevalue)
void	computeProdXS_JJH (float &mevalue)
void	computeProdXS_JH (float &mevalue)
void	computeProdXS_VH (float &mevalue, bool includeHiggsDecay=false)
void	computeProdXS_ttH (float &mevalue, int topProcess, int topDecay=0)
void	get_XPropagator (TVar::ResonancePropagatorScheme scheme, float &prop)
void	set_Process (TVar::Process process_, TVar::MatrixElement me_, TVar::Production production_)
void	set_Verbosity (TVar::VerbosityLevel verbosity_)
void	set_LeptonInterference (TVar::LeptonInterference myLepInterf)
void	set_TempCandidate (SimpleParticleCollection_t *pDaughters, SimpleParticleCollection_t *pAssociated=0, SimpleParticleCollection_t *pMothers=0, bool isGen=false)
void	set_RenFacScaleMode (TVar::EventScaleScheme renormalizationSch, TVar::EventScaleScheme factorizationSch, double ren_sf, double fac_sf)
void	set_LHAgrid (const char *path, int pdfmember=0)
void	set_PrimaryHiggsMass (double mh)
void	set_CandidateDecayMode (TVar::CandidateDecayMode mode)
void	set_CurrentCandidateFromIndex (unsigned int icand)
void	set_CurrentCandidate (MELACandidate *cand)
void	set_InputEvent (SimpleParticleCollection_t *pDaughters, SimpleParticleCollection_t *pAssociated=0, SimpleParticleCollection_t *pMothers=0, bool isGen=false)
void	append_TopCandidate (SimpleParticleCollection_t *TopDaughters)
void	set_mHiggs (double mh_, int index)
void	set_wHiggs (double gah_, int index)
void	set_mHiggs_wHiggs (double mh_, double gah_, int index)
void	reset_Mass (double inmass, int ipart)
void	reset_Width (double inmass, int ipart)
void	reset_QuarkMasses ()
void	reset_MCFM_EWKParameters (double ext_Gf, double ext_aemmz, double ext_mW, double ext_mZ, double ext_xW, int ext_ewscheme=3)
void	resetPerEvent ()
void	reset_InputEvent ()
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)
void	set_SpinZeroContact (double selfDHzzpcoupl[SIZE_HVV][2], double selfDHzpzpcoupl[SIZE_HVV][2], double selfDHwwpcoupl[SIZE_HVV][2], double selfDHwpwpcoupl[SIZE_HVV][2])
void	set_SpinOneCouplings (double selfDZqqcoupl[SIZE_ZQQ][2], double selfDZvvcoupl[SIZE_ZVV][2])
void	set_SpinTwoCouplings (double selfDGqqcoupl[SIZE_GQQ][2], double selfDGggcoupl[SIZE_GGG][2], double selfDGvvcoupl[SIZE_GVV][2])
void	set_SpinTwoContact (double selfDGvvpcoupl[SIZE_GVV][2], double selfDGvpvpcoupl[SIZE_GVV][2])
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)
void	set_aTQGCCouplings (double selfDaTQGCcoupl[SIZE_ATQGC][2])
void	set_AZffCouplings (double selfDAZffcoupl[SIZE_AZff][2])
std::vector< TLorentzVector >	Calculate4Momentum (double Mx, double M1, double M2, double theta, double theta1, double theta2, double Phi1, double Phi)
MelaIO *	get_IORecord ()
double	get_PrimaryMass (int ipart)
double	get_PrimaryHiggsMass ()
double	get_PrimaryWidth (int ipart)
double	get_HiggsWidthAtPoleMass (double mass)
MELACandidate *	get_CurrentCandidate ()
int	get_CurrentCandidateIndex ()
int	get_NCandidates ()
std::vector< MELATopCandidate_t * > *	get_TopCandidateCollection ()

Protected Member Functions
void	build ()

Protected Attributes
TVar::Process	processModel
TVar::MatrixElement	processME
TVar::Production	processProduction
TVar::VerbosityLevel	processVerbosity
TVar::LeptonInterference	processLeptonInterference
double	EBEAM
double	mHiggs [nSupportedHiggses]
double	wHiggs [nSupportedHiggses]
TEvtProb	Xcal2
SpinZeroCouplings *	selfD_SpinZeroCouplings
SpinOneCouplings *	selfD_SpinOneCouplings
SpinTwoCouplings *	selfD_SpinTwoCouplings
VprimeCouplings *	selfD_VprimeCouplings
aTQGCCouplings *	selfD_aTQGCCouplings
AZffCouplings *	selfD_AZffCouplings
MELACandidate *	melaCand
std::vector< MELAParticle * >	tmpPartList
std::vector< MELACandidate * >	tmpCandList

Detailed Description

Definition at line 9 of file ZZMatrixElement.h.

Constructor & Destructor Documentation

ZZMatrixElement() [1/2]

ZZMatrixElement::ZZMatrixElement	(	const char *	pathtoPDFSet,
		int	PDFMember,
		const char *	pathtoHiggsCSandWidth,
		double	ebeam,
		TVar::VerbosityLevel	verbosity
	)

Definition at line 12 of file ZZMatrixElement.cc.

    :
  processVerbosity(verbosity),
  processLeptonInterference(TVar::DefaultLeptonInterf),
  EBEAM(ebeam),
  Xcal2(pathtoHiggsCSandWidth, EBEAM, pathtoPDFSet, PDFMember, verbosity),
  melaCand(0)
{
  if (processVerbosity>=TVar::DEBUG) MELAout << "Begin ZZMatrixElement constructor" << endl;
  build();
  if (processVerbosity>=TVar::DEBUG) MELAout << "End ZZMatrixElement constructor" << endl;
}

ZZMatrixElement() [2/2]

ZZMatrixElement::ZZMatrixElement

(

const ZZMatrixElement &

other

)

Definition at line 29 of file ZZMatrixElement.cc.

                                                             :
processVerbosity(other.processVerbosity),
processLeptonInterference(other.processLeptonInterference),
EBEAM(other.EBEAM),
Xcal2(other.Xcal2),
melaCand(0) // 0 is correct in the copy constructor
{
  if (processVerbosity>=TVar::DEBUG) MELAout << "Begin ZZMatrixElement copy constructor" << endl;
  build();
  if (processVerbosity>=TVar::DEBUG) MELAout << "End ZZMatrixElement copy constructor" << endl;
}

~ZZMatrixElement()

ZZMatrixElement::~ZZMatrixElement

(

)

Definition at line 59 of file ZZMatrixElement.cc.

                                 {
  if (processVerbosity>=TVar::DEBUG) MELAout << "Begin ZZMatrixElement destructor" << endl;
  resetPerEvent();
  reset_InputEvent();
  if (processVerbosity>=TVar::DEBUG) MELAout << "End ZZMatrixElement destructor" << endl;
}

Member Function Documentation

append_TopCandidate()

void ZZMatrixElement::append_TopCandidate

(

SimpleParticleCollection_t *

TopDaughters

)

Definition at line 178 of file ZZMatrixElement.cc.

{ Xcal2.AppendTopCandidate(TopDaughters); }

build()

void ZZMatrixElement::build ( )

protected

Definition at line 41 of file ZZMatrixElement.cc.

                           {
  if (processVerbosity>=TVar::DEBUG) MELAout << "Begin ZZMatrixElement::build" << endl;
 
  // Set default parameters explicitly
  set_PrimaryHiggsMass(125.);
  set_mHiggs(125., 0); set_wHiggs(-1., 0);
  set_mHiggs(-1., 1); set_wHiggs(-1, 1);
 
  selfD_SpinZeroCouplings = Xcal2.GetSelfDSpinZeroCouplings();
  selfD_SpinOneCouplings = Xcal2.GetSelfDSpinOneCouplings();
  selfD_SpinTwoCouplings = Xcal2.GetSelfDSpinTwoCouplings();
  selfD_VprimeCouplings = Xcal2.GetSelfDVprimeCouplings();
  selfD_aTQGCCouplings = Xcal2.GetSelfDaTQGCCouplings();
  selfD_AZffCouplings = Xcal2.GetSelfDAZffCouplings();
 
  if (processVerbosity>=TVar::DEBUG) MELAout << "End ZZMatrixElement::build" << endl;
}

Calculate4Momentum()

std::vector< TLorentzVector > ZZMatrixElement::Calculate4Momentum	(	double	Mx,
		double	M1,
		double	M2,
		double	theta,
		double	theta1,
		double	theta2,
		double	Phi1,
		double	Phi
	)

Definition at line 67 of file ZZMatrixElement.cc.

                                                                                                                                                                 {
  double phi1, phi2;
  phi1=TMath::Pi()-Phi1;
  phi2=Phi1+Phi;
 
  double gamma1=1, gamma2=1, beta1=0, beta2=0;
 
  if (M1>0. && Mx>0.){
    gamma1=(Mx*Mx+M1*M1-M2*M2)/(2*Mx*M1);
    beta1=sqrt(1.-1./(gamma1*gamma1));
  }
  if (M2>0. && Mx>0.){
    gamma2=(Mx*Mx-M1*M1+M2*M2)/(2*Mx*M2);
    beta2=sqrt(1.-1./(gamma2*gamma2));
  }
 
  //gluon 4 vectors
  TLorentzVector p1CM(0, 0, Mx/2, Mx/2);
  TLorentzVector p2CM(0, 0, -Mx/2, Mx/2);
 
  //vector boson 4 vectors
  TLorentzVector kZ1(gamma1*M1*sin(theta)*beta1, 0, gamma1*M1*cos(theta)*beta1, gamma1*M1);
  TLorentzVector kZ2(-gamma2*M2*sin(theta)*beta2, 0, -gamma2*M2*cos(theta)*beta2, gamma2*M2);
 
  //Rotation and Boost matrices. Note gamma1*beta1*M1=gamma2*beta2*M2.
 
  TLorentzRotation Z1ToZ, Z2ToZ;
 
  Z1ToZ.Boost(0, 0, beta1);
  Z2ToZ.Boost(0, 0, beta2);
  Z1ToZ.RotateY(theta);
  Z2ToZ.RotateY(TMath::Pi()+theta);
 
 
  //fermion 4 vectors in vector boson rest frame
 
  TLorentzVector p3Z1((M1/2)*sin(theta1)*cos(phi1), (M1/2)*sin(theta1)*sin(phi1), (M1/2)*cos(theta1), (M1/2)*1);
  TLorentzVector p4Z1(-(M1/2)*sin(theta1)*cos(phi1), -(M1/2)*sin(theta1)*sin(phi1), -(M1/2)*cos(theta1), (M1/2)*1);
  TLorentzVector p5Z2((M2/2)*sin(theta2)*cos(phi2), (M2/2)*sin(theta2)*sin(phi2), (M2/2)*cos(theta2), (M2/2)*1);
  TLorentzVector p6Z2(-(M2/2)*sin(theta2)*cos(phi2), -(M2/2)*sin(theta2)*sin(phi2), -(M2/2)*cos(theta2), (M2/2)*1);
 
  // fermions 4 vectors in CM frame
 
  TLorentzVector p3CM, p4CM, p5CM, p6CM;
 
  p3CM=Z1ToZ*p3Z1;
  p4CM=Z1ToZ*p4Z1;
  p5CM=Z2ToZ*p5Z2;
  p6CM=Z2ToZ*p6Z2;
 
  vector<TLorentzVector> p;
 
  p.push_back(p3CM);
  p.push_back(p4CM);
  p.push_back(p5CM);
  p.push_back(p6CM);
 
  return p;
}

computeProdXS_JH()

void ZZMatrixElement::computeProdXS_JH

(

float &

mevalue

)

Definition at line 398 of file ZZMatrixElement.cc.

   {
  melaCand = get_CurrentCandidate();
 
  if (melaCand){
    mevalue  = Xcal2.XsecCalcXJ();
  }
 
  resetPerEvent();
  return;
}

computeProdXS_JJH()

void ZZMatrixElement::computeProdXS_JJH

(

float &

mevalue

)

Definition at line 384 of file ZZMatrixElement.cc.

   {
  melaCand = get_CurrentCandidate();
 
  if (melaCand){
    mevalue  = Xcal2.XsecCalcXJJ();
  }
 
  resetPerEvent();
  return;
}

computeProdXS_ttH()

void ZZMatrixElement::computeProdXS_ttH	(	float &	mevalue,
		int	topProcess,
		int	topDecay = 0
	)

Definition at line 435 of file ZZMatrixElement.cc.

   {
  melaCand = get_CurrentCandidate();
 
  if (melaCand){
    double zzmass = melaCand->m();
    if (processME == TVar::MCFM){ for (int jh=0; jh<(int)nSupportedHiggses; jh++) Xcal2.SetHiggsMass(mHiggs[jh], wHiggs[jh], jh+1); }
    else Xcal2.SetHiggsMass(zzmass, wHiggs[0], -1);
 
    if (processProduction == TVar::ttH || processProduction == TVar::bbH){
      mevalue  = Xcal2.XsecCalc_TTX(
        topProcess, topDecay
        );
    }
  }
 
  resetPerEvent();
  return;
}

computeProdXS_VH()

void ZZMatrixElement::computeProdXS_VH	(	float &	mevalue,
		bool	includeHiggsDecay = false
	)

Definition at line 412 of file ZZMatrixElement.cc.

   {
  melaCand = get_CurrentCandidate();
 
  if (melaCand){
    double zzmass = melaCand->m();
    if (processME == TVar::MCFM){ for (int jh=0; jh<(int)nSupportedHiggses; jh++) Xcal2.SetHiggsMass(mHiggs[jh], wHiggs[jh], jh+1); }
    else Xcal2.SetHiggsMass(zzmass, wHiggs[0], -1);
 
    mevalue  = Xcal2.XsecCalc_VX(
      includeHiggsDecay
      );
  }
 
  resetPerEvent();
  return;
}

computeProdXS_VVHVV()

void ZZMatrixElement::computeProdXS_VVHVV

(

float &

mevalue

)

Definition at line 364 of file ZZMatrixElement.cc.

   {
  melaCand = get_CurrentCandidate();
 
  if (melaCand){
    double zzmass = melaCand->m();
    if (processME == TVar::MCFM){
      for (int jh=0; jh<(int)nSupportedHiggses; jh++) Xcal2.SetHiggsMass(mHiggs[jh], wHiggs[jh], jh+1);
    }
    else Xcal2.SetHiggsMass(zzmass, wHiggs[0], -1);
 
    mevalue = Xcal2.XsecCalc_VVXVV();
  }
 
  resetPerEvent();
  return;
}

computeXS()

void ZZMatrixElement::computeXS

(

float &

mevalue

)

Definition at line 344 of file ZZMatrixElement.cc.

   {
  melaCand = get_CurrentCandidate();
 
  if (melaCand){
    double zzmass = melaCand->m();
    if (processME == TVar::MCFM){
      for (int jh=0; jh<(int)nSupportedHiggses; jh++) Xcal2.SetHiggsMass(mHiggs[jh], wHiggs[jh], jh+1);
    }
    else Xcal2.SetHiggsMass(zzmass, wHiggs[0], -1);
 
    mevalue = Xcal2.XsecCalc_XVV();
  }
 
  resetPerEvent();
  return;
}

get_CurrentCandidate()

MELACandidate * ZZMatrixElement::get_CurrentCandidate

(

)

Definition at line 232 of file ZZMatrixElement.cc.

{ return Xcal2.GetCurrentCandidate(); }

get_CurrentCandidateIndex()

int ZZMatrixElement::get_CurrentCandidateIndex

(

)

Definition at line 233 of file ZZMatrixElement.cc.

{ return Xcal2.GetCurrentCandidateIndex(); }

get_HiggsWidthAtPoleMass()

double ZZMatrixElement::get_HiggsWidthAtPoleMass

(

double

mass

)

Definition at line 231 of file ZZMatrixElement.cc.

{ return Xcal2.GetHiggsWidthAtPoleMass(mass); }

get_IORecord()

MelaIO * ZZMatrixElement::get_IORecord

(

)

Definition at line 228 of file ZZMatrixElement.cc.

{ return Xcal2.GetIORecord(); }

get_NCandidates()

int ZZMatrixElement::get_NCandidates

(

)

Definition at line 234 of file ZZMatrixElement.cc.

{ return Xcal2.GetNCandidates(); }

get_PrimaryHiggsMass()

double ZZMatrixElement::get_PrimaryHiggsMass ( )

inline

Definition at line 141 of file ZZMatrixElement.h.

{ return get_PrimaryMass(25); }

get_PrimaryMass()

double ZZMatrixElement::get_PrimaryMass

(

int

ipart

)

Definition at line 229 of file ZZMatrixElement.cc.

{ return Xcal2.GetPrimaryMass(ipart); }

get_PrimaryWidth()

double ZZMatrixElement::get_PrimaryWidth

(

int

ipart

)

Definition at line 230 of file ZZMatrixElement.cc.

{ return Xcal2.GetPrimaryWidth(ipart); }

get_TopCandidateCollection()

std::vector< MELATopCandidate_t * > * ZZMatrixElement::get_TopCandidateCollection

(

)

Definition at line 235 of file ZZMatrixElement.cc.

{ return Xcal2.GetTopCandidates(); }

get_XPropagator()

void ZZMatrixElement::get_XPropagator	(	TVar::ResonancePropagatorScheme	scheme,
		float &	prop
	)

Definition at line 459 of file ZZMatrixElement.cc.

                                                                                      {
  prop=0.;
  melaCand = get_CurrentCandidate();
 
  if (melaCand){
    Xcal2.SetHiggsMass(mHiggs[0], wHiggs[0], -1);
    prop=Xcal2.GetXPropagator(scheme);
  }
 
  resetPerEvent();
  return;
}

reset_InputEvent()

void ZZMatrixElement::reset_InputEvent

(

)

Definition at line 225 of file ZZMatrixElement.cc.

{ Xcal2.ResetInputEvent(); }

reset_Mass()

void ZZMatrixElement::reset_Mass	(	double	inmass,
		int	ipart
	)

Definition at line 197 of file ZZMatrixElement.cc.

{ Xcal2.ResetMass(inmass, ipart); }

reset_MCFM_EWKParameters()

void ZZMatrixElement::reset_MCFM_EWKParameters	(	double	ext_Gf,
		double	ext_aemmz,
		double	ext_mW,
		double	ext_mZ,
		double	ext_xW,
		int	ext_ewscheme = 3
	)

Definition at line 200 of file ZZMatrixElement.cc.

                                                                                                                                            {
  Xcal2.ResetMCFM_EWKParameters(ext_Gf, ext_aemmz, ext_mW, ext_mZ, ext_xW, ext_ewscheme);
}

reset_QuarkMasses()

void ZZMatrixElement::reset_QuarkMasses

(

)

Definition at line 199 of file ZZMatrixElement.cc.

{ Xcal2.ResetQuarkMasses(); }

reset_Width()

void ZZMatrixElement::reset_Width	(	double	inmass,
		int	ipart
	)

Definition at line 198 of file ZZMatrixElement.cc.

{ Xcal2.ResetWidth(inwidth, ipart); }

resetPerEvent()

void ZZMatrixElement::resetPerEvent

(

)

Definition at line 205 of file ZZMatrixElement.cc.

                                   {
  // Protection against forgetfulness; custom width has to be set per-computation
  set_mHiggs(Xcal2.GetPrimaryHiggsMass(), 0); // Sets mHiggs[0]
  if (wHiggs[0]>=0.) set_wHiggs(-1., 0);
 
  if (mHiggs[1]>=0.) set_mHiggs(-1., 1);
  if (wHiggs[1]>=0.) set_wHiggs(-1., 1);
 
  Xcal2.SetHiggsMass(mHiggs[0], -1, -1);
 
  // Return back to default lepton interference settings after each calculation
  if (processLeptonInterference != TVar::DefaultLeptonInterf) set_LeptonInterference(TVar::DefaultLeptonInterf);
 
  // Delete the temporary input objects owned
  for (unsigned int ic=0; ic<tmpCandList.size(); ic++){ if (tmpCandList.at(ic)) delete tmpCandList.at(ic); }
  //for (unsigned int itc=0; itc<tmpTopCandList.size(); itc++){ if (tmpTopCandList.at(itc)) delete tmpTopCandList.at(itc); }
  for (unsigned int ip=0; ip<tmpPartList.size(); ip++){ if (tmpPartList.at(ip)) delete tmpPartList.at(ip); }
  melaCand=nullptr;
}

set_aTQGCCouplings()

void ZZMatrixElement::set_aTQGCCouplings

(

double

selfDaTQGCcoupl[SIZE_ATQGC][2]

)

Definition at line 330 of file ZZMatrixElement.cc.

 {
  for (int ic=0; ic<SIZE_ATQGC; ic++) selfD_aTQGCCouplings->SetATQGCCouplings(ic, selfDaTQGCcoupl[ic][0], selfDaTQGCcoupl[ic][1]);
}

set_AZffCouplings()

void ZZMatrixElement::set_AZffCouplings

(

double

selfDAZffcoupl[SIZE_AZff][2]

)

Definition at line 335 of file ZZMatrixElement.cc.

 {
  for (int ic=0; ic<SIZE_AZff; ic++) selfD_AZffCouplings->SetAZffCouplings(ic, selfDAZffcoupl[ic][0], selfDAZffcoupl[ic][1]);
}

set_CandidateDecayMode()

void ZZMatrixElement::set_CandidateDecayMode

(

TVar::CandidateDecayMode

mode

)

Definition at line 141 of file ZZMatrixElement.cc.

{ Xcal2.SetCandidateDecayMode(mode); }

set_CurrentCandidate()

void ZZMatrixElement::set_CurrentCandidate

(

MELACandidate *

cand

)

Definition at line 144 of file ZZMatrixElement.cc.

{ Xcal2.SetCurrentCandidate(cand); }

set_CurrentCandidateFromIndex()

void ZZMatrixElement::set_CurrentCandidateFromIndex

(

unsigned int

icand

)

Definition at line 143 of file ZZMatrixElement.cc.

{ Xcal2.SetCurrentCandidateFromIndex(icand); }

set_InputEvent()

void ZZMatrixElement::set_InputEvent	(	SimpleParticleCollection_t *	pDaughters,
		SimpleParticleCollection_t *	pAssociated = 0,
		SimpleParticleCollection_t *	pMothers = 0,
		bool	isGen = false
	)

Definition at line 145 of file ZZMatrixElement.cc.

   { // Adds a new candidate to Xcal2
  Xcal2.SetInputEvent(
    pDaughters,
    pAssociated,
    pMothers,
    isGen
    );
}

set_LeptonInterference()

void ZZMatrixElement::set_LeptonInterference

(

TVar::LeptonInterference

myLepInterf

)

Definition at line 135 of file ZZMatrixElement.cc.

{ processLeptonInterference = lepInterf_; Xcal2.SetLeptonInterf(processLeptonInterference); }

set_LHAgrid()

void ZZMatrixElement::set_LHAgrid	(	const char *	path,
		int	pdfmember = 0
	)

Definition at line 137 of file ZZMatrixElement.cc.

{ Xcal2.Set_LHAgrid(path, pdfmember); }

set_mHiggs()

void ZZMatrixElement::set_mHiggs	(	double	mh_,
		int	index
	)

Definition at line 180 of file ZZMatrixElement.cc.

                                                     {
  if (index<nSupportedHiggses && index>=0) mHiggs[index] = mh_;
  else MELAerr << "ZZMatrixElement::set_mHiggs: Only resonances 0 (regular) and 1 (additional, possibly high-mass) are supported" << endl;
}

set_mHiggs_wHiggs()

void ZZMatrixElement::set_mHiggs_wHiggs	(	double	mh_,
		double	gah_,
		int	index
	)

Definition at line 188 of file ZZMatrixElement.cc.

                                                                         {
  if (index<nSupportedHiggses && index>=0){
    mHiggs[index] = mh_;
    wHiggs[index] = gah_;
  }
  else MELAerr << "ZZMatrixElement::set_mHiggs_wHiggs: Only resonances 0 (regular) and 1 (additional, possibly high-mass) are supported" << endl;
}

set_PrimaryHiggsMass()

void ZZMatrixElement::set_PrimaryHiggsMass

(

double

mh

)

Definition at line 142 of file ZZMatrixElement.cc.

{ Xcal2.SetPrimaryHiggsMass(mh); }

set_Process()

void ZZMatrixElement::set_Process	(	TVar::Process	process_,
		TVar::MatrixElement	me_,
		TVar::Production	production_
	)

Definition at line 128 of file ZZMatrixElement.cc.

                                                                                                        {
  processModel = process_;
  processME = me_;
  processProduction = production_;
  Xcal2.SetProcess(processModel, processME, processProduction);
}

set_RenFacScaleMode()

void ZZMatrixElement::set_RenFacScaleMode	(	TVar::EventScaleScheme	renormalizationSch,
		TVar::EventScaleScheme	factorizationSch,
		double	ren_sf,
		double	fac_sf
	)

Definition at line 138 of file ZZMatrixElement.cc.

                                                                                                                                                     {
  Xcal2.SetRenFacScaleMode(renormalizationSch, factorizationSch, ren_sf, fac_sf);
}

set_SpinOneCouplings()

void ZZMatrixElement::set_SpinOneCouplings	(	double	selfDZqqcoupl[SIZE_ZQQ][2],
		double	selfDZvvcoupl[SIZE_ZVV][2]
	)

Definition at line 290 of file ZZMatrixElement.cc.

   {
  for (int ic=0; ic<SIZE_ZQQ; ic++) selfD_SpinOneCouplings->SetZQQCouplings(ic, selfDZqqcoupl[ic][0], selfDZqqcoupl[ic][1]);
  for (int ic=0; ic<SIZE_ZVV; ic++) selfD_SpinOneCouplings->SetZVVCouplings(ic, selfDZvvcoupl[ic][0], selfDZvvcoupl[ic][1]);
}

set_SpinTwoContact()

void ZZMatrixElement::set_SpinTwoContact	(	double	selfDGvvpcoupl[SIZE_GVV][2],
		double	selfDGvpvpcoupl[SIZE_GVV][2]
	)

Definition at line 306 of file ZZMatrixElement.cc.

 {
  for (int ic=0; ic<SIZE_GVV; ic++){
    selfD_SpinTwoCouplings->SetGVVpCouplings(ic, selfDGvvpcoupl[ic][0], selfDGvvpcoupl[ic][1]);
    selfD_SpinTwoCouplings->SetGVpVpCouplings(ic, selfDGvpvpcoupl[ic][0], selfDGvpvpcoupl[ic][1]);
  }
}

set_SpinTwoCouplings()

void ZZMatrixElement::set_SpinTwoCouplings	(	double	selfDGqqcoupl[SIZE_GQQ][2],
		double	selfDGggcoupl[SIZE_GGG][2],
		double	selfDGvvcoupl[SIZE_GVV][2]
	)

Definition at line 297 of file ZZMatrixElement.cc.

   {
  for (int ic=0; ic<SIZE_GQQ; ic++) selfD_SpinTwoCouplings->SetGQQCouplings(ic, selfDGqqcoupl[ic][0], selfDGqqcoupl[ic][1]);
  for (int ic=0; ic<SIZE_GGG; ic++) selfD_SpinTwoCouplings->SetGGGCouplings(ic, selfDGggcoupl[ic][0], selfDGggcoupl[ic][1]);
  for (int ic=0; ic<SIZE_GVV; ic++) selfD_SpinTwoCouplings->SetGVVCouplings(ic, selfDGvvcoupl[ic][0], selfDGvvcoupl[ic][1]);
}

set_SpinZeroContact()

void ZZMatrixElement::set_SpinZeroContact	(	double	selfDHzzpcoupl[SIZE_HVV][2],
		double	selfDHzpzpcoupl[SIZE_HVV][2],
		double	selfDHwwpcoupl[SIZE_HVV][2],
		double	selfDHwpwpcoupl[SIZE_HVV][2]
	)

Definition at line 277 of file ZZMatrixElement.cc.

   {
  for (int ic=0; ic<SIZE_HVV; ic++){
    selfD_SpinZeroCouplings->SetHVVpCouplings(ic, selfDHzzpcoupl[ic][0], selfDHzzpcoupl[ic][1]);
    selfD_SpinZeroCouplings->SetHVpVpCouplings(ic, selfDHzpzpcoupl[ic][0], selfDHzpzpcoupl[ic][1]);
    selfD_SpinZeroCouplings->SetHVVpCouplings(ic, selfDHwwpcoupl[ic][0], selfDHwwpcoupl[ic][1], true);
    selfD_SpinZeroCouplings->SetHVpVpCouplings(ic, selfDHwpwpcoupl[ic][0], selfDHwpwpcoupl[ic][1], true);
  }
}

set_SpinZeroCouplings()

void ZZMatrixElement::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 at line 238 of file ZZMatrixElement.cc.

   {
  Xcal2.AllowSeparateWWCouplings(diffHWW);
  for (int jh=1; jh<=(int)nSupportedHiggses; jh++){
    for (int ic=0; ic<SIZE_HGG; ic++) selfD_SpinZeroCouplings->SetHGGCouplings(ic, selfDHggcoupl[jh-1][ic][0], selfDHggcoupl[jh-1][ic][1], 1, jh);
    for (int ic=0; ic<SIZE_HGG; ic++) selfD_SpinZeroCouplings->SetHGGCouplings(ic, selfDHg4g4coupl[jh-1][ic][0], selfDHg4g4coupl[jh-1][ic][1], 2, jh);
 
    for (int ic=0; ic<SIZE_HQQ; ic++) selfD_SpinZeroCouplings->SetHQQCouplings(ic, selfDHqqcoupl[jh-1][ic][0], selfDHqqcoupl[jh-1][ic][1], 0, jh);
    for (int ic=0; ic<SIZE_HQQ; ic++) selfD_SpinZeroCouplings->SetHQQCouplings(ic, selfDHbbcoupl[jh-1][ic][0], selfDHbbcoupl[jh-1][ic][1], 5, jh);
    for (int ic=0; ic<SIZE_HQQ; ic++) selfD_SpinZeroCouplings->SetHQQCouplings(ic, selfDHttcoupl[jh-1][ic][0], selfDHttcoupl[jh-1][ic][1], 6, jh);
    for (int ic=0; ic<SIZE_HQQ; ic++) selfD_SpinZeroCouplings->SetHQQCouplings(ic, selfDHb4b4coupl[jh-1][ic][0], selfDHb4b4coupl[jh-1][ic][1], 7, jh);
    for (int ic=0; ic<SIZE_HQQ; ic++) selfD_SpinZeroCouplings->SetHQQCouplings(ic, selfDHt4t4coupl[jh-1][ic][0], selfDHt4t4coupl[jh-1][ic][1], 8, jh);
 
    for (int ic=0; ic<SIZE_HVV; ic++) selfD_SpinZeroCouplings->SetHVVCouplings(ic, selfDHzzcoupl[jh-1][ic][0], selfDHzzcoupl[jh-1][ic][1], false, jh);
    for (int ic=0; ic<SIZE_HVV; ic++) selfD_SpinZeroCouplings->SetHVVCouplings(ic, selfDHwwcoupl[jh-1][ic][0], selfDHwwcoupl[jh-1][ic][1], true, jh);
    for (int ik=0; ik<SIZE_HVV_CQSQ; ik++){
      for (int ig=0; ig<SIZE_HVV_LAMBDAQSQ; ig++){
        selfD_SpinZeroCouplings->SetHVVLambdaQ2(ig, ik, selfDHzzLambda_qsq[jh-1][ig][ik], false, jh);
        selfD_SpinZeroCouplings->SetHVVLambdaQ2(ig, ik, selfDHwwLambda_qsq[jh-1][ig][ik], true, jh);
      }
      selfD_SpinZeroCouplings->SetHVVSignCQ2(ik, selfDHzzCLambda_qsq[jh-1][ik], false, jh);
      selfD_SpinZeroCouplings->SetHVVSignCQ2(ik, selfDHwwCLambda_qsq[jh-1][ik], true, jh);
    }
  }
}

set_TempCandidate()

void ZZMatrixElement::set_TempCandidate	(	SimpleParticleCollection_t *	pDaughters,
		SimpleParticleCollection_t *	pAssociated = 0,
		SimpleParticleCollection_t *	pMothers = 0,
		bool	isGen = false
	)

Definition at line 159 of file ZZMatrixElement.cc.

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

set_Verbosity()

void ZZMatrixElement::set_Verbosity

(

TVar::VerbosityLevel

verbosity_

)

Definition at line 134 of file ZZMatrixElement.cc.

{ processVerbosity = verbosity_; Xcal2.SetVerbosity(verbosity_); }

set_VprimeContactCouplings()

void ZZMatrixElement::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 at line 315 of file ZZMatrixElement.cc.

 {
  for (int ic=0; ic<SIZE_Vpff; ic++){
    selfD_VprimeCouplings->SetVpffCouplings(ic, selfDZpffcoupl[ic][0], selfDZpffcoupl[ic][1]);
    selfD_VprimeCouplings->SetVpffCouplings(ic, selfDWpffcoupl[ic][0], selfDWpffcoupl[ic][1], true);
  }
  selfD_VprimeCouplings->SetZPrimeMassWidth(M_Zprime, Ga_Zprime);
  selfD_VprimeCouplings->SetWPrimeMassWidth(M_Wprime, Ga_Wprime);
}

set_wHiggs()

void ZZMatrixElement::set_wHiggs	(	double	gah_,
		int	index
	)

Definition at line 184 of file ZZMatrixElement.cc.

                                                      {
  if (index<nSupportedHiggses && index>=0) wHiggs[index] = (double)gah_;
  else MELAerr << "ZZMatrixElement::set_wHiggs: Only resonances 0 (regular) and 1 (additional, possibly high-mass) are supported" << endl;
}

Member Data Documentation

EBEAM

double ZZMatrixElement::EBEAM

protected

Definition at line 157 of file ZZMatrixElement.h.

melaCand

MELACandidate* ZZMatrixElement::melaCand

protected

Definition at line 169 of file ZZMatrixElement.h.

mHiggs

double ZZMatrixElement::mHiggs[nSupportedHiggses]

protected

Definition at line 158 of file ZZMatrixElement.h.

processLeptonInterference

TVar::LeptonInterference ZZMatrixElement::processLeptonInterference

protected

Definition at line 155 of file ZZMatrixElement.h.

processME

TVar::MatrixElement ZZMatrixElement::processME

protected

Definition at line 152 of file ZZMatrixElement.h.

processModel

TVar::Process ZZMatrixElement::processModel

protected

Definition at line 151 of file ZZMatrixElement.h.

processProduction

TVar::Production ZZMatrixElement::processProduction

protected

Definition at line 153 of file ZZMatrixElement.h.

processVerbosity

TVar::VerbosityLevel ZZMatrixElement::processVerbosity

protected

Definition at line 154 of file ZZMatrixElement.h.

selfD_aTQGCCouplings

aTQGCCouplings* ZZMatrixElement::selfD_aTQGCCouplings

protected

Definition at line 166 of file ZZMatrixElement.h.

selfD_AZffCouplings

AZffCouplings* ZZMatrixElement::selfD_AZffCouplings

protected

Definition at line 167 of file ZZMatrixElement.h.

selfD_SpinOneCouplings

SpinOneCouplings* ZZMatrixElement::selfD_SpinOneCouplings

protected

Definition at line 163 of file ZZMatrixElement.h.

selfD_SpinTwoCouplings

SpinTwoCouplings* ZZMatrixElement::selfD_SpinTwoCouplings

protected

Definition at line 164 of file ZZMatrixElement.h.

selfD_SpinZeroCouplings

SpinZeroCouplings* ZZMatrixElement::selfD_SpinZeroCouplings

protected

Definition at line 162 of file ZZMatrixElement.h.

selfD_VprimeCouplings

VprimeCouplings* ZZMatrixElement::selfD_VprimeCouplings

protected

Definition at line 165 of file ZZMatrixElement.h.

tmpCandList

std::vector<MELACandidate*> ZZMatrixElement::tmpCandList

protected

Definition at line 173 of file ZZMatrixElement.h.

tmpPartList

std::vector<MELAParticle*> ZZMatrixElement::tmpPartList

protected

Definition at line 170 of file ZZMatrixElement.h.

wHiggs

double ZZMatrixElement::wHiggs[nSupportedHiggses]

protected

Definition at line 159 of file ZZMatrixElement.h.

Xcal2

TEvtProb ZZMatrixElement::Xcal2

protected

Definition at line 160 of file ZZMatrixElement.h.

---

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

• MELA/interface/ZZMatrixElement.h
• MELA/src/ZZMatrixElement.cc