JHUGen MELA  JHUGen v7.5.6, MELA v2.4.2
Matrix element calculations as used in JHUGen.
TEvtProb.cc
Go to the documentation of this file.
1 //-----------------------------------------------------------------------------
2 //
3 // Class EventProb Module
4 //
5 // EventProb Module
6 //
7 // March 21 2011
8 // S. Jindariani ([email protected])
9 // Y. Gao ([email protected])
10 // K. Burkett ([email protected])
11 //-----------------------------------------------------------------------------
12 
13 #include "TEvtProb.hh"
14 #include "MELAStreamHelpers.hh"
15 #include "MadMela.h"
16 
17 
18 using namespace std;
19 using namespace TUtil;
20 using MELAStreamHelpers::MELAout;
21 using MELAStreamHelpers::MELAerr;
22 
23 
24 //-----------------------------------------------------------------------------
25 // Constructors and Destructor
26 //-----------------------------------------------------------------------------
27 TEvtProb::TEvtProb(
28  const char* pathtoXSW, double ebeam, const char* pathtoPDFSet, int PDFMember, TVar::VerbosityLevel verbosity_
29  ) :
30  pathtoPDFSet_(pathtoPDFSet),
31  PDFMember_(PDFMember),
32  verbosity(verbosity_),
33  EBEAM(ebeam),
34  myCSW_(pathtoXSW)
35 {
36  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb constructor" << endl;
37 
38  /***** Build everything except MELAHXSWidth *****/
39  Build();
40 
41  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb constructor" << endl;
42 }
43 TEvtProb::TEvtProb(const TEvtProb& other) :
44 pathtoPDFSet_(other.pathtoPDFSet_),
45 PDFMember_(other.PDFMember_),
46 verbosity(other.verbosity),
47 EBEAM(other.EBEAM),
48 myCSW_(other.myCSW_)
49 {
50 
51 }
52 void TEvtProb::Build(){
53  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb::Build" << endl;
54 
55  /***** Initialize lepton interference scheme *****/
56  SetLeptonInterf(TVar::DefaultLeptonInterf);
57 
58  /***** Initialize MCFM *****/
59  InitializeMCFM();
60 
61  /***** Initialize JHUGen *****/
62  InitializeJHUGen(pathtoPDFSet_, PDFMember_);
63 
64  /***** Cross-initializations *****/
65  CrossInitialize();
66 
67  /***** Initialize input event properties *****/
68  ResetInputEvent();
69  SetCandidateDecayMode(TVar::CandidateDecay_ZZ);
70 
71  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb::Build" << endl;
72 }
73 
74 
75 TEvtProb::~TEvtProb(){
76  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb destructor" << endl;
77 
78  ResetInputEvent();
79 
80  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb destructor" << endl;
81 }
82 
83 void TEvtProb::InitializeMCFM(){
84  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb::InitializeMCFM" << endl;
85 
86  mcfm_init_((char *)"input.DAT", (char *)"./");
87  if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::TEvtProb: mcfm_init successful" << endl;
88  energy_.sqrts = 2.*EBEAM;
89  ResetQuarkMasses();
90  ResetMCFM_EWKParameters(
91  1.16639E-05, // Gf
92  1./128., // alpha_EW
93  80.399, // mW
94  91.1876, // mZ
95  0.23119, // xw=sin**2(thetaW)
96  3 // MCFM EW scheme
97  );
98  qlinit_();
99  // First resonance constant parameters
100  spinzerohiggs_anomcoupl_.LambdaBSM=1000;
101  spinzerohiggs_anomcoupl_.Lambda_Q=10000;
102  spinzerohiggs_anomcoupl_.Lambda_z1=10000;
103  spinzerohiggs_anomcoupl_.Lambda_z2=10000;
104  spinzerohiggs_anomcoupl_.Lambda_z3=10000;
105  spinzerohiggs_anomcoupl_.Lambda_z4=10000;
106  spinzerohiggs_anomcoupl_.Lambda_z11=100; spinzerohiggs_anomcoupl_.Lambda_z21=100; spinzerohiggs_anomcoupl_.Lambda_z31=100; spinzerohiggs_anomcoupl_.Lambda_z41=100;
107  spinzerohiggs_anomcoupl_.Lambda_z12=100; spinzerohiggs_anomcoupl_.Lambda_z22=100; spinzerohiggs_anomcoupl_.Lambda_z32=100; spinzerohiggs_anomcoupl_.Lambda_z42=100;
108  spinzerohiggs_anomcoupl_.Lambda_z10=100; spinzerohiggs_anomcoupl_.Lambda_z20=100; spinzerohiggs_anomcoupl_.Lambda_z30=100; spinzerohiggs_anomcoupl_.Lambda_z40=100;
109  spinzerohiggs_anomcoupl_.Lambda_zgs1=10000;
110  spinzerohiggs_anomcoupl_.Lambda_w1=10000;
111  spinzerohiggs_anomcoupl_.Lambda_w2=10000;
112  spinzerohiggs_anomcoupl_.Lambda_w3=10000;
113  spinzerohiggs_anomcoupl_.Lambda_w4=10000;
114  spinzerohiggs_anomcoupl_.Lambda_w11=100; spinzerohiggs_anomcoupl_.Lambda_w21=100; spinzerohiggs_anomcoupl_.Lambda_w31=100; spinzerohiggs_anomcoupl_.Lambda_w41=100;
115  spinzerohiggs_anomcoupl_.Lambda_w12=100; spinzerohiggs_anomcoupl_.Lambda_w22=100; spinzerohiggs_anomcoupl_.Lambda_w32=100; spinzerohiggs_anomcoupl_.Lambda_w42=100;
116  spinzerohiggs_anomcoupl_.Lambda_w10=100; spinzerohiggs_anomcoupl_.Lambda_w20=100; spinzerohiggs_anomcoupl_.Lambda_w30=100; spinzerohiggs_anomcoupl_.Lambda_w40=100;
117  // Second resonance constant parameters
118  spinzerohiggs_anomcoupl_.Lambda2BSM=1000;
119  spinzerohiggs_anomcoupl_.Lambda2_Q=10000;
120  spinzerohiggs_anomcoupl_.Lambda2_z1=10000;
121  spinzerohiggs_anomcoupl_.Lambda2_z2=10000;
122  spinzerohiggs_anomcoupl_.Lambda2_z3=10000;
123  spinzerohiggs_anomcoupl_.Lambda2_z4=10000;
124  spinzerohiggs_anomcoupl_.Lambda2_z11=100; spinzerohiggs_anomcoupl_.Lambda2_z21=100; spinzerohiggs_anomcoupl_.Lambda2_z31=100; spinzerohiggs_anomcoupl_.Lambda2_z41=100;
125  spinzerohiggs_anomcoupl_.Lambda2_z12=100; spinzerohiggs_anomcoupl_.Lambda2_z22=100; spinzerohiggs_anomcoupl_.Lambda2_z32=100; spinzerohiggs_anomcoupl_.Lambda2_z42=100;
126  spinzerohiggs_anomcoupl_.Lambda2_z10=100; spinzerohiggs_anomcoupl_.Lambda2_z20=100; spinzerohiggs_anomcoupl_.Lambda2_z30=100; spinzerohiggs_anomcoupl_.Lambda2_z40=100;
127  spinzerohiggs_anomcoupl_.Lambda2_zgs1=10000;
128  spinzerohiggs_anomcoupl_.Lambda2_w1=10000;
129  spinzerohiggs_anomcoupl_.Lambda2_w2=10000;
130  spinzerohiggs_anomcoupl_.Lambda2_w3=10000;
131  spinzerohiggs_anomcoupl_.Lambda2_w4=10000;
132  spinzerohiggs_anomcoupl_.Lambda2_w11=100; spinzerohiggs_anomcoupl_.Lambda2_w21=100; spinzerohiggs_anomcoupl_.Lambda2_w31=100; spinzerohiggs_anomcoupl_.Lambda2_w41=100;
133  spinzerohiggs_anomcoupl_.Lambda2_w12=100; spinzerohiggs_anomcoupl_.Lambda2_w22=100; spinzerohiggs_anomcoupl_.Lambda2_w32=100; spinzerohiggs_anomcoupl_.Lambda2_w42=100;
134  spinzerohiggs_anomcoupl_.Lambda2_w10=100; spinzerohiggs_anomcoupl_.Lambda2_w20=100; spinzerohiggs_anomcoupl_.Lambda2_w30=100; spinzerohiggs_anomcoupl_.Lambda2_w40=100;
135  // Switches for spin-0 Higgs couplings
136  spinzerohiggs_anomcoupl_.channeltoggle_stu=2; spinzerohiggs_anomcoupl_.vvhvvtoggle_vbfvh=2;
137  spinzerohiggs_anomcoupl_.AnomalCouplPR=1; spinzerohiggs_anomcoupl_.AnomalCouplDK=1;
138  // Constant parameters for all processes
139  qlfirst_.qlfirst=false;
140 
141  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb::InitializeMCFM" << endl;
142 }
143 void TEvtProb::InitializeJHUGen(const char* pathtoPDFSet, int PDFMember){
144  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb::InitializeJHUGen" << endl;
145 
146  InitJHUGenMELA(pathtoPDFSet, PDFMember, 2.*EBEAM);
147 
148  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb::InitializeJHUGen" << endl;
149 }
150 void TEvtProb::CrossInitialize(){
151  if (verbosity>=TVar::DEBUG) MELAout << "Begin TEvtProb::CrossInitialize" << endl;
152 
153  /**** Initialize alphaS *****/
154  const double GeV=1./100.;
155  double init_Q;
156  int init_nl, init_nf;
157  __modjhugenmela_MOD_getpdfconstants(&init_Q, &init_nl, &init_nf);
158  init_Q /= GeV;
159  if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::TEvtProb: Initializing the PDF with initial Q=" << init_Q << ", nloops=" << init_nl << ", nf=" << init_nf << endl;
160  SetAlphaS(init_Q, init_Q, 1, 1, init_nl, init_nf, "cteq6_l"); // "cteq6_l" is just some dummy variable
161 
162  /**** Initialize MCFM CKM initializers in cabib.f using JHUGen defaults *****/
163  int i, j;
164  i=2; j=1;
165  cabib_.Vud = __modparameters_MOD_ckmbare(&i, &j);
166  i=2; j=3;
167  cabib_.Vus = __modparameters_MOD_ckmbare(&i, &j);
168  i=2; j=5;
169  cabib_.Vub = __modparameters_MOD_ckmbare(&i, &j);
170  i=4; j=1;
171  cabib_.Vcd = __modparameters_MOD_ckmbare(&i, &j);
172  i=4; j=3;
173  cabib_.Vcs = __modparameters_MOD_ckmbare(&i, &j);
174  i=4; j=5;
175  cabib_.Vcb = __modparameters_MOD_ckmbare(&i, &j);
176 
177  /***** Initialize schemes *****/
178  ResetCouplings();
179  ResetRenFacScaleMode();
180  SetPrimaryHiggsMass(125.); // Should come after InitializeMCFM and InitializeJHUGen
181 
182  if (verbosity>=TVar::DEBUG) MELAout << "End TEvtProb::CrossInitialize" << endl;
183 }
184 
185 // Set NNPDF driver path
186 void TEvtProb::Set_LHAgrid(const char* path, int pdfmember){
187  char path_nnpdf_c[200];
188  sprintf(path_nnpdf_c, "%s", path);
189  int pathLength = strlen(path_nnpdf_c);
190  nnpdfdriver_(path_nnpdf_c, &pathLength);
191  nninitpdf_(&pdfmember);
192 }
193 void TEvtProb::SetProcess(TVar::Process proc, TVar::MatrixElement me, TVar::Production prod){
194  matrixElement = me;
195  production = prod;
196  // In case s-channel processes are passed for JHUGen ME, flip them back to JHUGen-specific productions.
197  if (matrixElement == TVar::JHUGen){
198  if (production == TVar::Had_ZH_S) production = TVar::Had_ZH;
199  else if (production == TVar::Had_WH_S) production = TVar::Had_WH;
200  else if (production == TVar::Lep_ZH_S) production = TVar::Lep_ZH;
201  else if (production == TVar::Lep_WH_S) production = TVar::Lep_WH;
202  else if (production == TVar::JJVBF_S) production = TVar::JJVBF;
203  else if (production == TVar::JJQCD_S) production = TVar::JJQCD;
204  }
205  process = proc;
206 }
207 void TEvtProb::SetVerbosity(TVar::VerbosityLevel tmp){ verbosity = tmp; }
208 void TEvtProb::SetLeptonInterf(TVar::LeptonInterference tmp){ leptonInterf = tmp; }
209 void TEvtProb::SetCandidateDecayMode(TVar::CandidateDecayMode mode){ PDGHelpers::setCandidateDecayMode(mode); }
210 void TEvtProb::SetRenFacScaleMode(TVar::EventScaleScheme renormalizationSch, TVar::EventScaleScheme factorizationSch, double ren_sf, double fac_sf){
211  event_scales.renomalizationScheme = renormalizationSch;
212  event_scales.factorizationScheme = factorizationSch;
213  event_scales.ren_scale_factor = ren_sf;
214  event_scales.fac_scale_factor = fac_sf;
215 }
216 const TVar::event_scales_type& TEvtProb::GetRenFacScaleMode() const{
217  return this->event_scales;
218 }
219 void TEvtProb::AllowSeparateWWCouplings(bool doAllow){ SetJHUGenDistinguishWWCouplings(doAllow); selfDSpinZeroCoupl.allow_WWZZSeparation(doAllow); }
220 void TEvtProb::SetPrimaryHiggsMass(double mass){ PrimaryHMass=mass; SetHiggsMass(PrimaryHMass, -1., -1); }
221 void TEvtProb::SetHiggsMass(double mass, double wHiggs, int whichResonance){
222  // Regular, first resonance
223  if (whichResonance==1 || whichResonance==-1){
224  if (mass<0.){
225  _hmass = -1;
226  _hwidth = 0;
227  }
228  else if (wHiggs<0.){
229  _hmass = mass;
230  _hwidth = myCSW_.HiggsWidth(_hmass);
231  }
232  else{
233  _hmass = mass;
234  _hwidth = wHiggs;
235  }
236  masses_mcfm_.hmass = _hmass;
237  masses_mcfm_.hwidth = _hwidth;
238 
239  madMela::mad_masses_.mdl_mh = _hmass;
240  madMela::widths_.mdl_wh = _hwidth;
241 
242  if (_hmass<0.) SetJHUGenHiggsMassWidth(0., _hwidth);
243  else SetJHUGenHiggsMassWidth(_hmass, _hwidth);
244  }
245 
246  // Second resonance
247  if (whichResonance==2){
248  if (mass<0.){
249  _h2mass = -1;
250  _h2width = 0;
251  }
252  else if (wHiggs<0.){
253  _h2mass = mass;
254  _h2width = myCSW_.HiggsWidth(_h2mass);
255  }
256  else{
257  _h2mass = mass;
258  _h2width = wHiggs;
259  }
260  spinzerohiggs_anomcoupl_.h2mass = _h2mass;
261  spinzerohiggs_anomcoupl_.h2width = _h2width;
262  //SetJHUGenHiggsMassWidth(_h2mass, _h2width); // Second resonance is not implemented in JHUGen yet.
263  }
264  else if (whichResonance==-1){
265  _h2mass = -1;
266  _h2width = 0;
267  spinzerohiggs_anomcoupl_.h2mass = _h2mass;
268  spinzerohiggs_anomcoupl_.h2width = _h2width;
269  //SetJHUGenHiggsMassWidth(_h2mass, _h2width); // Second resonance is not implemented in JHUGen yet.
270  }
271 
272  if (verbosity>=TVar::DEBUG) MELAout
273  << "TEvtProb::SetHiggsMass(" << mass << ", " << wHiggs << ", " << whichResonance << "):\n"
274  << '\t' << "hmass: " << _hmass << ", " << _hwidth << '\n'
275  << '\t' << "h2mass: " << _h2mass << ", " << _h2width << '\n'
276  << '\t' << "MCFM hmass: " << masses_mcfm_.hmass << ", " << masses_mcfm_.hwidth << '\n'
277  << '\t' << "MCFM h2mass: " << spinzerohiggs_anomcoupl_.h2mass << ", " << spinzerohiggs_anomcoupl_.h2width << endl;
278 }
279 void TEvtProb::SetInputEvent(
280  SimpleParticleCollection_t* pDaughters,
281  SimpleParticleCollection_t* pAssociated,
282  SimpleParticleCollection_t* pMothers,
283  bool isGen
284  ){
285  MELACandidate* cand = ConvertVectorFormat(
286  pDaughters,
287  pAssociated,
288  pMothers,
289  isGen,
290  &particleList, &candList // push_back is done automatically
291  );
292  if (cand) melaCand=cand;
293 }
294 void TEvtProb::AppendTopCandidate(SimpleParticleCollection_t* TopDaughters){
295  if (!CheckInputPresent()){
296  if (verbosity>=TVar::INFO) MELAout << "TEvtProb::AppendTopCandidate: No MELACandidates are present to append this top!" << endl;
297  return;
298  }
299  MELATopCandidate_t* cand = ConvertThreeBodyDecayCandidate(
300  TopDaughters,
301  &particleList, &topCandList // push_back is done automatically
302  );
303  if (cand) melaCand->addAssociatedTop(cand);
304 }
305 void TEvtProb::SetRcdCandPtr(){ RcdME.melaCand = melaCand; }
306 void TEvtProb::SetCurrentCandidateFromIndex(unsigned int icand){
307  if (candList.size()>icand) melaCand = candList.at(icand);
308  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::SetCurrentCandidateFromIndex: icand=" << icand << ">=candList.size()=" << candList.size() << endl;
309 }
310 void TEvtProb::SetCurrentCandidate(MELACandidate* cand){
311  melaCand = cand;
312  if (verbosity>=TVar::INFO && !melaCand) MELAout << "TEvtProb::SetCurrentCandidate: BE CAREFUL! !melaCand!" << endl;
313  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;
314 }
315 
316 // Reset functions
317 void TEvtProb::ResetIORecord(){ RcdME.reset(); }
318 void TEvtProb::ResetRenFacScaleMode(){ SetRenFacScaleMode(TVar::DefaultScaleScheme, TVar::DefaultScaleScheme, 0.5, 0.5); }
319 void TEvtProb::ResetMass(double inmass, int ipart){ TUtil::SetMass(inmass, ipart); }
320 void TEvtProb::SetZprimeMassWidth(double inmass, double inwidth){ this->ResetMass(inmass, 32); this->ResetWidth(inwidth, 32); }
321 void TEvtProb::SetWprimeMassWidth(double inmass, double inwidth){ this->ResetMass(inmass, 34); this->ResetWidth(inwidth, 34); }
322 void TEvtProb::ResetWidth(double inwidth, int ipart){ TUtil::SetDecayWidth(inwidth, ipart); }
323 void TEvtProb::ResetQuarkMasses(){
324  ResetMass(1e-3, 1); // d
325  ResetMass(5e-3, 2); // u
326  ResetMass(1e-1, 3); // s
327  ResetMass(1.275, 4); // c
328  ResetMass(4.75, 5); // b
329  ResetMass(173.2, 6); // t
330  ResetMass(1e5, 7); // bprime
331  ResetMass(1e5, 8); // tprime
332  //coupling_(); // Already called at b and t twice!
333 }
334 void TEvtProb::ResetMCFM_EWKParameters(double ext_Gf, double ext_aemmz, double ext_mW, double ext_mZ, double ext_xW, int ext_ewscheme){ SetEwkCouplingParameters(ext_Gf, ext_aemmz, ext_mW, ext_mZ, ext_xW, ext_ewscheme); }
335 void TEvtProb::ResetCouplings(){
336  selfDSpinZeroCoupl.reset();
337  selfDSpinOneCoupl.reset();
338  selfDSpinTwoCoupl.reset();
339  selfDVprimeCoupl.reset();
340  selfDaTQGCCoupl.reset();
341  selfDAZffCoupl.reset();
342  AllowSeparateWWCouplings(false);
343  ResetAmplitudeIncludes();
344 }
345 void TEvtProb::ResetInputEvent(){
346  RcdME.melaCand = 0;
347  melaCand = 0;
348 
349  // Clear bookkeeping objects
350  // Clean MELACandidates first since they contain all other objects
351  for (MELACandidate* tmpPart:candList) delete tmpPart;
352  candList.clear();
353  // Clean MELATopCandidates next since they contain MELAParticles
354  for (MELATopCandidate_t* tmpPart:topCandList) delete tmpPart;
355  topCandList.clear();
356  // Clean all remaining MELAParticles
357  for (MELAParticle* tmpPart:particleList) delete tmpPart;
358  particleList.clear();
359 }
360 
361 // Get-functions
362 MELAHXSWidth const* TEvtProb::GetHXSWidthEstimator() const{ return &myCSW_; }
363 SpinZeroCouplings* TEvtProb::GetSelfDSpinZeroCouplings(){ return selfDSpinZeroCoupl.getRef(); }
364 SpinOneCouplings* TEvtProb::GetSelfDSpinOneCouplings(){ return selfDSpinOneCoupl.getRef(); }
365 SpinTwoCouplings* TEvtProb::GetSelfDSpinTwoCouplings(){ return selfDSpinTwoCoupl.getRef(); }
366 VprimeCouplings* TEvtProb::GetSelfDVprimeCouplings(){ return selfDVprimeCoupl.getRef(); }
367 aTQGCCouplings* TEvtProb::GetSelfDaTQGCCouplings(){ return selfDaTQGCCoupl.getRef(); }
368 AZffCouplings* TEvtProb::GetSelfDAZffCouplings() { return selfDAZffCoupl.getRef(); }
369 double TEvtProb::GetPrimaryHiggsMass(){ return PrimaryHMass; }
370 double TEvtProb::GetPrimaryMass(int ipart){
371  if (PDGHelpers::isAHiggs(ipart)) return GetPrimaryHiggsMass();
372  else return TUtil::GetMass(ipart);
373 }
374 double TEvtProb::GetPrimaryWidth(int ipart){
375  if (PDGHelpers::isAHiggs(ipart)) return myCSW_.HiggsWidth(GetPrimaryHiggsMass());
376  else return TUtil::GetDecayWidth(ipart);
377 }
378 double TEvtProb::GetHiggsWidthAtPoleMass(double mass){
379  if (mass>0.) return myCSW_.HiggsWidth(mass);
380  else return -1.;
381 }
382 MelaIO* TEvtProb::GetIORecord(){ return RcdME.getRef(); }
383 MELACandidate* TEvtProb::GetCurrentCandidate(){ return melaCand; }
384 int TEvtProb::GetCurrentCandidateIndex(){
385  if (!melaCand) return -1;
386  for (unsigned int icand=0; icand<candList.size(); icand++){
387  if (candList.at(icand)==melaCand) return (int)icand;
388  }
389  return -1;
390 }
391 int TEvtProb::GetNCandidates(){ return (static_cast<int>(candList.size())); }
392 std::vector<MELATopCandidate_t*>* TEvtProb::GetTopCandidates(){ return &topCandList; }
393 
394 
395 // Check/test functions
396 bool TEvtProb::CheckInputPresent(){
397  if (!melaCand){
398  if (verbosity>=TVar::INFO) MELAout
399  << "TEvtProb::CheckInputPresent: melaCand==" << melaCand << " is nullPtr!"
400  << endl;
401  if (candList.empty()) return false;
402  else{
403  SetCurrentCandidateFromIndex(candList.size()-1);
404  if (verbosity>=TVar::INFO) MELAout << "TEvtProb::CheckInputPresent: melaCand now points to the latest candidate (cand[" << (candList.size()-1) << "])" << endl;
405  }
406  }
407  SetRcdCandPtr(); // If input event is present, set the RcdME pointer to melaCand
408  return true;
409 }
410 
411 
412 // ME functions
413 // Cross-section calculations for H + 0 jet
414 double TEvtProb::XsecCalc_XVV(){
415  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_XVV" << endl;
416  double dXsec=0;
417  ResetIORecord();
418  if (!CheckInputPresent()) return dXsec;
419  bool useMCFM = matrixElement == TVar::MCFM;
420  bool calculateME=false;
421  bool needBSMHiggs=false;
422  bool needAZff=CheckSelfDCouplings_AZff();
423  if (useMCFM){
424  if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XVV: Try MCFM" << endl;
425  needBSMHiggs = CheckSelfDCouplings_Hgg() || CheckSelfDCouplings_Htt() || CheckSelfDCouplings_Hbb() || CheckSelfDCouplings_HVV();
426  if (needBSMHiggs || needAZff) SetLeptonInterf(TVar::InterfOn); // All anomalous coupling computations have to use lepton interference
427 
428  calculateME = (
429  production == TVar::ZZGG ||
430  production == TVar::ZZQQB ||
431  production == TVar::ZZQQB_STU ||
432  production == TVar::ZZQQB_S ||
433  production == TVar::ZZQQB_TU ||
434  production == TVar::ZZINDEPENDENT ||
435  (production == TVar::JJQCD && process == TVar::bkgZJets) // Allow Z+jets to be computed here.
436  );
437  if (calculateME){
438  SetMCFMSpinZeroCouplings(needBSMHiggs, &selfDSpinZeroCoupl, false);
439  SetMCFMAZffCouplings(needAZff, &selfDAZffCoupl);
440  dXsec = SumMatrixElementPDF(process, production, matrixElement, leptonInterf, &event_scales, &RcdME, EBEAM, verbosity);
441  }
442  else if (verbosity>=TVar::INFO) MELAout << "TEvtProb::XsecCalc_XVV: MCFM_chooser failed to determine the process configuration." << endl;
443  }
444  else if (matrixElement == TVar::JHUGen){
445  if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XVV::Try JHUGen" << endl;
446  AllowSeparateWWCouplings(false); // HZZ couplings are used for both in spin-0
447  // all the possible couplings
448  double Hggcoupl[SIZE_HGG][2] ={ { 0 } };
449  double Hvvcoupl[SIZE_HVV][2] ={ { 0 } };
450  double HvvLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ] ={ { 0 } };
451  int HvvCLambda_qsq[SIZE_HVV_CQSQ] ={ 0 };
452 
453  double Hvvpcoupl[SIZE_HVV][2] = { { 0 } };
454  double Hvpvpcoupl[SIZE_HVV][2] = { { 0 } };
455 
456  double Zqqcoupl[SIZE_ZQQ][2] ={ { 0 } };
457  double Zvvcoupl[SIZE_ZVV][2] ={ { 0 } };
458 
459  double Gqqcoupl[SIZE_GQQ][2] ={ { 0 } };
460  double Gggcoupl[SIZE_GGG][2] ={ { 0 } };
461  double Gvvcoupl[SIZE_GVV][2] ={ { 0 } };
462  double Gvvpcoupl[SIZE_GVV][2] ={ { 0 } };
463  double Gvpvpcoupl[SIZE_GVV][2] ={ { 0 } };
464 
465  double Zpffcoupl[SIZE_Vpff][2] ={ { 0 } };
466  double Wpffcoupl[SIZE_Vpff][2] ={ { 0 } };
467  double M_Zprime = -1;
468  double Ga_Zprime = 0;
469  double M_Wprime = -1;
470  double Ga_Wprime = 0;
471 
472  double AZffcoupl[SIZE_AZff][2] ={ { 0 } };
473  for (int j=0; j<2; j++){
474  for (int i=0; i<SIZE_AZff; i++){
475  AZffcoupl[i][j] = (selfDAZffCoupl.AZffcoupl)[i][j];
476  }
477  }
478  SetJHUGenAZffCouplings(needAZff, AZffcoupl);
479 
480  //
481  // set spin 0 default numbers
482  //
483  // By default set the Spin 0 couplings for SM case (0+m)
484  Hggcoupl[gHIGGS_GG_2][0]=1.0;
485  for (int ic=0; ic<SIZE_HVV_LAMBDAQSQ; ic++){ for (int ik=0; ik<SIZE_HVV_CQSQ; ik++) HvvLambda_qsq[ic][ik]=100.; }
486  //
487  // set spin 1 default numbers
488  //
489  Zqqcoupl[gZPRIME_QQ_LEFT][0]=1.0;
490  Zqqcoupl[gZPRIME_QQ_RIGHT][0]=1.0;
491  //
492  // set spin 2 default numbers
493  //
494  Gqqcoupl[gGRAVITON_QQ_LEFT][0]=1.0;
495  Gqqcoupl[gGRAVITON_QQ_RIGHT][0]=1.0;
496 
497  bool isSpinZero = false;
498  bool isSpinOne = false;
499  bool isSpinTwo = false;
500 
501  // 0m+
502  if (process == TVar::HSMHiggs){
503  Hvvcoupl[gHIGGS_VV_1][0]=1.;
504  isSpinZero = true;
505  }
506  // 0+L1
507  else if (process == TVar::H0_g1prime2){
508  Hvvcoupl[gHIGGS_VV_1_PRIME2][0] = 1.;
509  isSpinZero = true;
510  }
511  // 0h+
512  else if (process == TVar::H0hplus){
513  Hvvcoupl[gHIGGS_VV_2][0] = 1.;
514  isSpinZero = true;
515  }
516  // 0-
517  else if (process == TVar::H0minus){
518  Hvvcoupl[gHIGGS_VV_4][0] = 1.;
519  isSpinZero = true;
520  }
521  else if (process == TVar::H0_Zgsg1prime2){
522  Hvvcoupl[gHIGGS_ZA_1_PRIME2][0] = 1.;
523  isSpinZero = true;
524  }
525  else if (process == TVar::H0_Zgs){
526  Hvvcoupl[gHIGGS_ZA_2][0] = 1.;
527  isSpinZero = true;
528  }
529  else if (process == TVar::H0_Zgs_PS){
530  Hvvcoupl[gHIGGS_ZA_4][0] = 1.;
531  isSpinZero = true;
532  }
533  else if (process == TVar::H0_gsgs){
534  Hvvcoupl[gHIGGS_AA_2][0] = 1.;
535  isSpinZero = true;
536  }
537  else if (process == TVar::H0_gsgs_PS){
538  Hvvcoupl[gHIGGS_AA_4][0] = 1.;
539  isSpinZero = true;
540  }
541  else if (process == TVar::SelfDefine_spin0){
542  for (int j=0; j<2; j++){
543  for (int i=0; i<SIZE_HGG; i++) Hggcoupl[i][j] = (selfDSpinZeroCoupl.Hggcoupl)[i][j];
544  for (int i=0; i<SIZE_HVV; i++){
545  Hvvcoupl[i][j] = (selfDSpinZeroCoupl.Hzzcoupl)[i][j];
546 
547  Hvvpcoupl[i][j] = (selfDSpinZeroCoupl.Hzzpcoupl)[i][j];
548  Hvpvpcoupl[i][j] = (selfDSpinZeroCoupl.Hzpzpcoupl)[i][j];
549  }
550  }
551  for (int j=0; j<SIZE_HVV_CQSQ; j++){
552  for (int i=0; i<SIZE_HVV_LAMBDAQSQ; i++) HvvLambda_qsq[i][j] = (selfDSpinZeroCoupl.HzzLambda_qsq)[i][j];
553  HvvCLambda_qsq[j] = (selfDSpinZeroCoupl.HzzCLambda_qsq)[j];
554  }
555 
556  if (verbosity>=TVar::DEBUG_VERBOSE){
557  for (int j=0; j<2; j++){
558  for (int i=0; i<SIZE_HGG; i++){ if ((selfDSpinZeroCoupl.Hggcoupl)[i][j]!=0.) MELAout << "Hggcoupl[" << i << "][" << j << "] = " << (selfDSpinZeroCoupl.Hggcoupl)[i][j] << endl; }
559  for (int i=0; i<SIZE_HVV; i++){
560  if ((selfDSpinZeroCoupl.Hzzcoupl)[i][j]!=0.) MELAout << "Hvvcoupl[" << i << "][" << j << "] = " << (selfDSpinZeroCoupl.Hzzcoupl)[i][j] << endl;
561 
562  if ((selfDSpinZeroCoupl.Hzzpcoupl)[i][j]!=0.) MELAout << "Hzzpcoupl[" << i << "][" << j << "] = " << (selfDSpinZeroCoupl.Hzzpcoupl)[i][j] << endl;
563  if ((selfDSpinZeroCoupl.Hzpzpcoupl)[i][j]!=0.) MELAout << "Hzpzpcoupl[" << i << "][" << j << "] = " << (selfDSpinZeroCoupl.Hzpzpcoupl)[i][j] << endl;
564  }
565  }
566  for (int j=0; j<SIZE_HVV_CQSQ; j++){
567  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; }
568  if ((selfDSpinZeroCoupl.HzzCLambda_qsq)[j]!=0.) MELAout << "HvvCLambda_qsq[" << j << "] = " << (selfDSpinZeroCoupl.HzzCLambda_qsq)[j] << endl;
569  }
570  }
571 
572  isSpinZero = true;
573  }
574 
575  // 1-
576  else if (process == TVar::H1minus){
577  Zvvcoupl[gZPRIME_VV_1][0]=1.;
578  isSpinOne = true;
579  }
580  // 1+
581  else if (process == TVar::H1plus){
582  Zvvcoupl[gZPRIME_VV_2][0]=1.;
583  isSpinOne = true;
584  }
585  else if (process == TVar::SelfDefine_spin1){
586  for (int j=0; j<2; j++){
587  for (int i=0; i<SIZE_ZQQ; i++) Zqqcoupl[i][j] = (selfDSpinOneCoupl.Zqqcoupl)[i][j];
588  for (int i=0; i<SIZE_ZVV; i++) Zvvcoupl[i][j] = (selfDSpinOneCoupl.Zvvcoupl)[i][j];
589  }
590  isSpinOne = true;
591  }
592 
593  // 2m+
594  else if (process == TVar::H2_g1g5){
595  Gggcoupl[gGRAVITON_GG_1][0]=1.;
596  Gvvcoupl[gGRAVITON_VV_1][0]=1.;
597  Gvvcoupl[gGRAVITON_VV_5][0]=1.;
598  isSpinTwo = true;
599  }
600  // 2b+
601  else if (process == TVar::H2_g1){
602  Gggcoupl[gGRAVITON_GG_1][0]=1.;
603  Gvvcoupl[gGRAVITON_VV_1][0]=1.;
604  isSpinTwo = true;
605  }
606  else if (process == TVar::H2_g2){
607  Gggcoupl[gGRAVITON_GG_2][0]=1.;
608  Gvvcoupl[gGRAVITON_VV_2][0]=1.;
609  isSpinTwo = true;
610  }
611  // 2h3plus
612  else if (process == TVar::H2_g3){
613  Gggcoupl[gGRAVITON_GG_3][0]=1.;
614  Gvvcoupl[gGRAVITON_VV_3][0]=1.;
615  isSpinTwo = true;
616  }
617  // 2h+
618  else if (process == TVar::H2_g4){
619  Gggcoupl[gGRAVITON_GG_4][0]=1.;
620  Gvvcoupl[gGRAVITON_VV_4][0]=1.;
621  isSpinTwo = true;
622  }
623  // 2b+
624  else if (process == TVar::H2_g5){
625  Gggcoupl[gGRAVITON_GG_1][0]=1.;
626  Gvvcoupl[gGRAVITON_VV_5][0]=1.;
627  isSpinTwo = true;
628  }
629  // 2h6+
630  else if (process == TVar::H2_g6){
631  Gggcoupl[gGRAVITON_GG_1][0]=1.;
632  Gvvcoupl[gGRAVITON_VV_6][0]=1.;
633  isSpinTwo = true;
634  }
635  // 2h7plus
636  else if (process == TVar::H2_g7){
637  Gggcoupl[gGRAVITON_GG_1][0]=1.;
638  Gvvcoupl[gGRAVITON_VV_7][0]=1.;
639  isSpinTwo = true;
640  }
641  // 2h-
642  else if (process == TVar::H2_g8){
643  Gggcoupl[gGRAVITON_GG_5][0]=1.;
644  Gvvcoupl[gGRAVITON_VV_8][0]=1.;
645  isSpinTwo = true;
646  }
647  // 2h9minus
648  else if (process == TVar::H2_g9){
649  Gggcoupl[gGRAVITON_GG_5][0]=1.;
650  Gvvcoupl[gGRAVITON_VV_9][0]=1.;
651  isSpinTwo = true;
652  }
653  // 2h10minus
654  else if (process == TVar::H2_g10){
655  Gggcoupl[gGRAVITON_GG_5][0]=1.;
656  Gvvcoupl[gGRAVITON_VV_10][0]=1.;
657  isSpinTwo = true;
658  }
659  else if (process == TVar::SelfDefine_spin2){
660  for (int j=0; j<2; j++){
661  for (int i=0; i<SIZE_GGG; i++) Gggcoupl[i][j] = (selfDSpinTwoCoupl.Gggcoupl)[i][j];
662  for (int i=0; i<SIZE_GQQ; i++) Gqqcoupl[i][j] = (selfDSpinTwoCoupl.Gqqcoupl)[i][j];
663  for (int i=0; i<SIZE_GVV; i++){
664  Gvvcoupl[i][j] = (selfDSpinTwoCoupl.Gvvcoupl)[i][j];
665  Gvvpcoupl[i][j] = (selfDSpinTwoCoupl.Gvvpcoupl)[i][j];
666  Gvpvpcoupl[i][j] = (selfDSpinTwoCoupl.Gvpvpcoupl)[i][j];
667  }
668  }
669  isSpinTwo = true;
670  }
671 
672  // Vprime / contact couplings
673  if (process == TVar::SelfDefine_spin0 || process == TVar::SelfDefine_spin2){
674  for (int j=0; j<2; j++){
675  for (int i=0; i<SIZE_Vpff; i++){
676  Zpffcoupl[i][j] = (selfDVprimeCoupl.Zpffcoupl)[i][j];
677  Wpffcoupl[i][j] = (selfDVprimeCoupl.Wpffcoupl)[i][j];
678  }
679  }
680  M_Zprime = selfDVprimeCoupl.M_Zprime;
681  Ga_Zprime = selfDVprimeCoupl.Ga_Zprime;
682  M_Wprime = selfDVprimeCoupl.M_Wprime;
683  Ga_Wprime = selfDVprimeCoupl.Ga_Wprime;
684 
685  if (verbosity>=TVar::DEBUG_VERBOSE){
686  for (int j=0; j<2; j++){
687  for (int i=0; i<SIZE_Vpff; i++){
688  if ((selfDVprimeCoupl.Zpffcoupl)[i][j]!=0.) MELAout << "Zpffcoupl[" << i << "][" << j << "] = " << (selfDVprimeCoupl.Zpffcoupl)[i][j] << endl;
689  if ((selfDVprimeCoupl.Wpffcoupl)[i][j]!=0.) MELAout << "Wpffcoupl[" << i << "][" << j << "] = " << (selfDVprimeCoupl.Wpffcoupl)[i][j] << endl;
690  }
691  }
692  MELAout << "M_Zprime = " << selfDVprimeCoupl.M_Zprime << endl;
693  MELAout << "Ga_Zprime = " << selfDVprimeCoupl.Ga_Zprime << endl;
694  MELAout << "M_Wprime = " << selfDVprimeCoupl.M_Wprime << endl;
695  MELAout << "Ga_Wprime = " << selfDVprimeCoupl.Ga_Wprime << endl;
696  }
697  }
698 
699  if (isSpinZero){
700  SetJHUGenSpinZeroGGCouplings(Hggcoupl);
701  SetJHUGenSpinZeroVVCouplings(Hvvcoupl, Hvvpcoupl, Hvpvpcoupl, HvvCLambda_qsq, HvvLambda_qsq, false);
702  SetJHUGenVprimeContactCouplings(Zpffcoupl, Wpffcoupl);
703  SetZprimeMassWidth(M_Zprime, Ga_Zprime);
704  SetWprimeMassWidth(M_Wprime, Ga_Wprime);
705  }
706  else if (isSpinOne) SetJHUGenSpinOneCouplings(Zqqcoupl, Zvvcoupl);
707  else if (isSpinTwo){
708  SetJHUGenSpinTwoCouplings(Gggcoupl, Gvvcoupl, Gvvpcoupl, Gvpvpcoupl, Gqqcoupl);
709  SetJHUGenVprimeContactCouplings(Zpffcoupl, Wpffcoupl);
710  SetZprimeMassWidth(M_Zprime, Ga_Zprime);
711  SetWprimeMassWidth(M_Wprime, Ga_Wprime);
712  }
713 
714  if (isSpinZero || isSpinOne || isSpinTwo) dXsec = JHUGenMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, verbosity);
715  else if (verbosity>=TVar::ERROR) MELAerr
716  << "TEvtProb::XsecCalc_XVV: JHUGen ME is not spin zero, one or two! The process is described by "
717  << "Process: " << TVar::ProcessName(process) << ", Production: " << TVar::ProductionName(production) << ", and ME: " << TVar::MatrixElementName(matrixElement)
718  << endl;
719  } // end of JHUGen matrix element calculations
720  else if (matrixElement == TVar::MADGRAPH){
721  calculateME = (
722  production == TVar::ZZGG ||
723  production == TVar::ZZQQB
724  );
725  if (calculateME){
726  SetMadgraphSpinZeroCouplings(&selfDSpinZeroCoupl);
727  // madMela::ggFSIG_update_all_coup_(); //calculates all couplings
728  dXsec = TUtil::MadgraphMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, verbosity);
729  }
730  else if (verbosity>=TVar::INFO) MELAout << "TEvtProb::XsecCalc_XVV: Madgraph_chooser failed to determine the process configuration." << endl;
731  }
732  if (verbosity >= TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XVV: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
733 
734  if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XVV::Reset couplings" << endl;
735  ResetCouplings(); // Should come first
736  if (useMCFM){ // Set defaults. Should come next...
737  if (needBSMHiggs) SetLeptonInterf(TVar::DefaultLeptonInterf);
738  if (calculateME) SetMCFMSpinZeroCouplings(false, &selfDSpinZeroCoupl, true); // ... because of this!
739  }
740  ResetRenFacScaleMode();
741  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_XVV" << endl;
742  return dXsec;
743 }
744 
745 // Cross-section calculations for H(->VV) + 2 jets
746 double TEvtProb::XsecCalc_VVXVV(){
747  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_VVXVV" << endl;
748  double dXsec=0;
749  ResetIORecord();
750  if (!CheckInputPresent()) return dXsec;
751 
752  bool useMCFM = matrixElement == TVar::MCFM;
753  bool needBSMHiggs=false;
754  bool needATQGC=false;
755  bool calculateME=false;
756  if (useMCFM){
757  if (verbosity>=TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_VVXVV: Try MCFM" << endl;
758  needBSMHiggs = CheckSelfDCouplings_HVV();
759  needATQGC = CheckSelfDCouplings_aTQGC();
760  if (needBSMHiggs || needATQGC) SetLeptonInterf(TVar::InterfOn); // All anomalous coupling computations have to use lepton interference
761 
762  calculateME = (
763  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
764  || 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
765  || 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
766  );
767  if (calculateME){
768  SetMCFMSpinZeroCouplings(needBSMHiggs, &selfDSpinZeroCoupl, false);
769  SetMCFMaTQGCCouplings(needATQGC, &selfDaTQGCCoupl);
770  dXsec = SumMatrixElementPDF(process, production, matrixElement, leptonInterf, &event_scales, &RcdME, EBEAM, verbosity);
771  }
772  else if (verbosity>=TVar::INFO) MELAout << "TEvtProb::XsecCalc_VVXVV: MCFM_chooser failed to determine the process configuration." << endl;
773  }
774  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::XsecCalc_VVXVV: Non-MCFM Mes are not supported." << endl;
775 
776  if (verbosity >= TVar::DEBUG) MELAout
777  << "Process " << TVar::ProcessName(process)
778  << " TEvtProb::XsecCalc_VVXVV: dXsec=" << dXsec
779  << endl;
780 
781  ResetCouplings(); // Should come first
782  if (useMCFM){ // Set defaults. Should come next...
783  if (needBSMHiggs || needATQGC) SetLeptonInterf(TVar::DefaultLeptonInterf);
784  if (calculateME){
785  // ... because of these!
786  SetMCFMSpinZeroCouplings(false, &selfDSpinZeroCoupl, true);
787  SetMCFMaTQGCCouplings(false, &selfDaTQGCCoupl);
788  }
789  }
790  ResetRenFacScaleMode();
791  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_VVXVV" << endl;
792  return dXsec;
793 }
794 
795 // Cross-section calculations for H + 2 jets
796 double TEvtProb::XsecCalcXJJ(){
797  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_XJJ" << endl;
798  if (matrixElement == TVar::MCFM) return XsecCalc_VVXVV();
799  double dXsec = 0;
800  ResetIORecord();
801  if (!CheckInputPresent()) return dXsec;
802 
803  if (matrixElement == TVar::JHUGen){
804  if (production == TVar::JJQCD){
805  double Hggcoupl[SIZE_HGG][2] ={ { 0 } };
806  if (process == TVar::HSMHiggs) Hggcoupl[gHIGGS_GG_2][0] = 1.;
807  else if (process == TVar::H0minus) Hggcoupl[gHIGGS_GG_4][0] = 1.;
808  else if (process == TVar::SelfDefine_spin0){
809  for (int j=0; j<2; j++){
810  for (int i=0; i<SIZE_HGG; i++) Hggcoupl[i][j] = (selfDSpinZeroCoupl.Hggcoupl)[i][j];
811  }
812  }
813  SetJHUGenSpinZeroGGCouplings(Hggcoupl);
814  }
815  else if (production == TVar::JJVBF){
816  double Hzzcoupl[SIZE_HVV][2] ={ { 0 } };
817  double Hwwcoupl[SIZE_HVV][2] ={ { 0 } };
818  double HzzLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ] ={ { 0 } };
819  double HwwLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ] ={ { 0 } };
820  int HzzCLambda_qsq[SIZE_HVV_CQSQ] ={ 0 };
821  int HwwCLambda_qsq[SIZE_HVV_CQSQ] ={ 0 };
822  double Hzzpcoupl[SIZE_HVV][2] = { { 0 } };
823  double Hzpzpcoupl[SIZE_HVV][2] = { { 0 } };
824  double Zpffcoupl[SIZE_Vpff][2] = { { 0 } };
825  double Hwwpcoupl[SIZE_HVV][2] = { { 0 } };
826  double Hwpwpcoupl[SIZE_HVV][2] = { { 0 } };
827  double Wpffcoupl[SIZE_Vpff][2] = { { 0 } };
828  double M_Zprime = -1;
829  double Ga_Zprime = 0;
830  double M_Wprime = -1;
831  double Ga_Wprime = 0;
832 
833  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.; } }
834  SetJHUGenDistinguishWWCouplings(false);
835 
836  if (process == TVar::HSMHiggs){ Hzzcoupl[gHIGGS_VV_1][0] = 1.; Hwwcoupl[gHIGGS_VV_1][0] = 1.; }
837  else if (process == TVar::H0_g1prime2){ Hzzcoupl[gHIGGS_VV_1_PRIME2][0] = 1.; Hwwcoupl[gHIGGS_VV_1_PRIME2][0] = 1.; }
838  else if (process == TVar::H0hplus){ Hzzcoupl[gHIGGS_VV_2][0] = 1.; Hwwcoupl[gHIGGS_VV_2][0] = 1.; }
839  else if (process == TVar::H0minus){ Hzzcoupl[gHIGGS_VV_4][0] = 1.; Hwwcoupl[gHIGGS_VV_4][0] = 1.; }
840  else if (process == TVar::H0_Zgsg1prime2){ Hzzcoupl[gHIGGS_ZA_1_PRIME2][0] = 1.; }
841  else if (process == TVar::H0_Zgs){ Hzzcoupl[gHIGGS_ZA_2][0] = 1.; }
842  else if (process == TVar::H0_Zgs_PS){ Hzzcoupl[gHIGGS_ZA_4][0] = 1.; }
843  else if (process == TVar::H0_gsgs){ Hzzcoupl[gHIGGS_AA_2][0] = 1.; }
844  else if (process == TVar::H0_gsgs_PS){ Hzzcoupl[gHIGGS_AA_4][0] = 1.; }
845  else if (process == TVar::SelfDefine_spin0){
846  for (int j=0; j<2; j++){
847  for (int i=0; i<SIZE_HVV; i++){
848  Hzzcoupl[i][j] = (selfDSpinZeroCoupl.Hzzcoupl)[i][j]; Hwwcoupl[i][j] = (selfDSpinZeroCoupl.Hwwcoupl)[i][j];
849  Hzzpcoupl[i][j] = (selfDSpinZeroCoupl.Hzzpcoupl)[i][j]; Hwwpcoupl[i][j] = (selfDSpinZeroCoupl.Hwwpcoupl)[i][j];
850  Hzpzpcoupl[i][j] = (selfDSpinZeroCoupl.Hzpzpcoupl)[i][j]; Hwpwpcoupl[i][j] = (selfDSpinZeroCoupl.Hwpwpcoupl)[i][j];
851  }
852  for (int i=0; i<SIZE_Vpff; i++){
853  Zpffcoupl[i][j] = (selfDVprimeCoupl.Zpffcoupl)[i][j]; Wpffcoupl[i][j] = (selfDVprimeCoupl.Wpffcoupl)[i][j];
854  }
855  }
856  M_Zprime = selfDVprimeCoupl.M_Zprime;
857  Ga_Zprime = selfDVprimeCoupl.Ga_Zprime;
858  M_Wprime = selfDVprimeCoupl.M_Wprime;
859  Ga_Wprime = selfDVprimeCoupl.Ga_Wprime;
860  for (int j=0; j<SIZE_HVV_CQSQ; j++){
861  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]; }
862  HzzCLambda_qsq[j] = (selfDSpinZeroCoupl.HzzCLambda_qsq)[j]; HwwCLambda_qsq[j] = (selfDSpinZeroCoupl.HwwCLambda_qsq)[j];
863  }
864  SetJHUGenDistinguishWWCouplings(selfDSpinZeroCoupl.separateWWZZcouplings);
865  }
866  SetJHUGenSpinZeroVVCouplings(Hzzcoupl, Hzzpcoupl, Hzpzpcoupl, HzzCLambda_qsq, HzzLambda_qsq, false);
867  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.
868  SetJHUGenVprimeContactCouplings(Zpffcoupl, Wpffcoupl);
869  SetZprimeMassWidth(M_Zprime, Ga_Zprime);
870  SetWprimeMassWidth(M_Wprime, Ga_Wprime);
871  }
872 
873  dXsec = HJJMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, verbosity);
874  if (verbosity >= TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XJJ: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
875  }
876  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::XsecCalc_XJJ: Non-JHUGen vbfMELA MEs are not supported!" << endl;
877 
878  ResetCouplings();
879  ResetRenFacScaleMode();
880  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_XJJ" << endl;
881  return dXsec;
882 }
883 
884 // Cross-section calculations for H + 1 jet (only SM for the moment)
885 double TEvtProb::XsecCalcXJ(){
886  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_XJ" << endl;
887  double dXsec = 0;
888  ResetIORecord();
889  if (!CheckInputPresent()) return dXsec;
890 
891  // Calculate the ME
892  if (matrixElement == TVar::JHUGen){
893  // Not currently supported, but implement anyway
894  double Hggcoupl[SIZE_HGG][2] ={ { 0 } };
895  if (process == TVar::HSMHiggs) Hggcoupl[gHIGGS_GG_2][0] = 1.;
896  else if (process == TVar::H0minus) Hggcoupl[gHIGGS_GG_4][0] = 1.;
897  else if (process == TVar::SelfDefine_spin0){
898  for (int j=0; j<2; j++){
899  for (int i=0; i<SIZE_HGG; i++) Hggcoupl[i][j] = (selfDSpinZeroCoupl.Hggcoupl)[i][j];
900  }
901  }
902  SetJHUGenSpinZeroGGCouplings(Hggcoupl);
903 
904  dXsec = HJJMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, verbosity);
905  if (verbosity >= TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_XJ: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
906  }
907  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::XsecCalc_XJ: Non-JHUGen HJ MEs are not supported!" << endl;
908 
909  ResetCouplings();
910  ResetRenFacScaleMode();
911  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_XJ" << endl;
912  return dXsec;
913 }
914 
915 // Cross-section calculations for VH
916 double TEvtProb::XsecCalc_VX(
917  bool includeHiggsDecay
918  ){
919  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_VX" << endl;
920  if (matrixElement == TVar::MCFM) return XsecCalc_VVXVV();
921  double dXsec = 0;
922  ResetIORecord();
923  if (!CheckInputPresent()) return dXsec;
924 
925  AllowSeparateWWCouplings(false); // HZZ couplings are used for both in spin-0
926 
927  if (matrixElement == TVar::JHUGen){
928  double Hvvcoupl[SIZE_HVV][2] ={ { 0 } };
929  double HvvLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ] ={ { 0 } };
930  int HvvCLambda_qsq[SIZE_HVV_CQSQ] ={ 0 };
931  double Hvvpcoupl[SIZE_HVV][2] = { { 0 } };
932  double Hvpvpcoupl[SIZE_HVV][2] = { { 0 } };
933  double Zpffcoupl[SIZE_Vpff][2] = { { 0 } };
934  double Wpffcoupl[SIZE_Vpff][2] = { { 0 } };
935  double M_Zprime = -1;
936  double Ga_Zprime = 0;
937  double M_Wprime = -1;
938  double Ga_Wprime = 0;
939 
940  for (int ic=0; ic<SIZE_HVV_LAMBDAQSQ; ic++){ for (int ik=0; ik<SIZE_HVV_CQSQ; ik++) HvvLambda_qsq[ic][ik]=100.; }
941 
942  if (process == TVar::HSMHiggs) Hvvcoupl[gHIGGS_VV_1][0] = 1.;
943  else if (process == TVar::H0_g1prime2) Hvvcoupl[gHIGGS_VV_1_PRIME2][0] = 1.;
944  else if (process == TVar::H0hplus) Hvvcoupl[gHIGGS_VV_2][0] = 1.;
945  else if (process == TVar::H0minus) Hvvcoupl[gHIGGS_VV_4][0] = 1.;
946  else if (process == TVar::H0_Zgsg1prime2) Hvvcoupl[gHIGGS_ZA_1_PRIME2][0] = 1.;
947  else if (process == TVar::H0_Zgs) Hvvcoupl[gHIGGS_ZA_2][0] = 1.;
948  else if (process == TVar::H0_Zgs_PS) Hvvcoupl[gHIGGS_ZA_4][0] = 1.;
949  else if (process == TVar::H0_gsgs) Hvvcoupl[gHIGGS_AA_2][0] = 1.;
950  else if (process == TVar::H0_gsgs_PS) Hvvcoupl[gHIGGS_AA_4][0] = 1.;
951  else if (process == TVar::SelfDefine_spin0){
952  for (int i=0; i<SIZE_HVV; i++){
953  for (int j=0; j<2; j++){
954  Hvvcoupl[i][j] = (selfDSpinZeroCoupl.Hzzcoupl)[i][j];
955 
956  Hvvpcoupl[i][j] = (selfDSpinZeroCoupl.Hzzpcoupl)[i][j];
957  Hvpvpcoupl[i][j] = (selfDSpinZeroCoupl.Hzpzpcoupl)[i][j];
958  }
959  }
960  for (int i=0; i<SIZE_Vpff; i++){
961  for (int j=0; j<2; j++){
962  Zpffcoupl[i][j] = (selfDVprimeCoupl.Zpffcoupl)[i][j];
963  Wpffcoupl[i][j] = (selfDVprimeCoupl.Wpffcoupl)[i][j];
964  }
965  }
966  M_Zprime = selfDVprimeCoupl.M_Zprime;
967  Ga_Zprime = selfDVprimeCoupl.Ga_Zprime;
968  M_Wprime = selfDVprimeCoupl.M_Wprime;
969  Ga_Wprime = selfDVprimeCoupl.Ga_Wprime;
970  for (int j=0; j<SIZE_HVV_CQSQ; j++){
971  for (int i=0; i<SIZE_HVV_LAMBDAQSQ; i++) HvvLambda_qsq[i][j] = (selfDSpinZeroCoupl.HzzLambda_qsq)[i][j];
972  HvvCLambda_qsq[j] = (selfDSpinZeroCoupl.HzzCLambda_qsq)[j];
973  }
974  }
975  SetJHUGenSpinZeroVVCouplings(Hvvcoupl, Hvvpcoupl, Hvpvpcoupl, HvvCLambda_qsq, HvvLambda_qsq, false);
976  SetJHUGenVprimeContactCouplings(Zpffcoupl, Wpffcoupl);
977  SetZprimeMassWidth(M_Zprime, Ga_Zprime);
978  SetWprimeMassWidth(M_Wprime, Ga_Wprime);
979 
980  dXsec = VHiggsMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, includeHiggsDecay, verbosity);
981  if (verbosity >= TVar::DEBUG) MELAout << "TEVtProb::XsecCalc_VX: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
982  } // end of JHUGen matrix element calculations
983  else if (verbosity>=TVar::ERROR) MELAerr << "TEVtProb::XsecCalc_VX: Non-JHUGen VH MEs are not supported!" << endl;
984 
985  ResetCouplings();
986  ResetRenFacScaleMode();
987  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_VX" << endl;
988  return dXsec;
989 }
990 
991 // Cross-section calculations for ttH or bbH
992 double TEvtProb::XsecCalc_TTX(
993  int topProcess, int topDecay
994  ){
995  if (verbosity>=TVar::DEBUG) MELAout << "Begin XsecCalc_TTX" << endl;
996  double dXsec = 0;
997  ResetIORecord();
998  if (!CheckInputPresent()) return dXsec;
999 
1000  if (matrixElement == TVar::JHUGen){
1001  double Hqqcoupl[SIZE_HQQ][2]={ { 0 } };
1002  if (process == TVar::HSMHiggs) Hqqcoupl[gHIGGS_KAPPA][0] = 1.;
1003  else if (process == TVar::H0minus) Hqqcoupl[gHIGGS_KAPPA_TILDE][0] = 1.;
1004  else if (process == TVar::SelfDefine_spin0){
1005  for (int i=0; i<SIZE_HQQ; i++){
1006  for (int j=0; j<2; j++){
1007  if ((selfDSpinZeroCoupl.Httcoupl)[i][j]!=0. && production == TVar::ttH) Hqqcoupl[i][j] = (selfDSpinZeroCoupl.Httcoupl)[i][j];
1008  else if ((selfDSpinZeroCoupl.Hbbcoupl)[i][j]!=0. && production == TVar::bbH) Hqqcoupl[i][j] = (selfDSpinZeroCoupl.Hbbcoupl)[i][j];
1009  else Hqqcoupl[i][j] = (selfDSpinZeroCoupl.Hqqcoupl)[i][j];
1010  }
1011  }
1012  }
1013  SetJHUGenSpinZeroQQCouplings(Hqqcoupl);
1014 
1015  if (production == TVar::ttH) dXsec = TTHiggsMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, topDecay, topProcess, verbosity);
1016  else if (production == TVar::bbH) dXsec = BBHiggsMatEl(process, production, matrixElement, &event_scales, &RcdME, EBEAM, topProcess, verbosity);
1017  else if (verbosity >= TVar::ERROR) MELAout << "TEvtProb::XsecCalc_TTX only supports ttH and bbH productions for the moment." << endl;
1018  if (verbosity >= TVar::DEBUG) MELAout << "TEvtProb::XsecCalc_TTX: Process " << TVar::ProcessName(process) << " dXsec=" << dXsec << endl;
1019  }
1020  else if (verbosity>=TVar::ERROR) MELAerr << "TEvtProb::XsecCalc_TTX: Non-JHUGen ttH MEs are not supported!" << endl;
1021 
1022  ResetCouplings();
1023  ResetRenFacScaleMode();
1024  if (verbosity>=TVar::DEBUG) MELAout << "End XsecCalc_TTX" << endl;
1025  return dXsec;
1026 }
1027 
1028 double TEvtProb::GetXPropagator(
1029  TVar::ResonancePropagatorScheme scheme
1030  ){
1031  if (!CheckInputPresent()) return 0.;
1032  return TUtil::ResonancePropagator(melaCand->m(), scheme);
1033 }
1034 
1035 
1036 bool TEvtProb::CheckSelfDCouplings_Hgg(){
1037  if (_hmass>=0. && _hwidth>0.){
1038  for (int vv = 0; vv < SIZE_HGG; vv++){
1039  if (
1040  (selfDSpinZeroCoupl.Hggcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hggcoupl)[vv][0] != 0
1041  ||
1042  (selfDSpinZeroCoupl.Hg4g4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hg4g4coupl)[vv][0] != 0
1043  ){
1044  return true;
1045  }
1046  }
1047  }
1048  if (_h2mass>=0. && _h2width>0.){
1049  for (int vv = 0; vv < SIZE_HGG; vv++){
1050  if (
1051  (selfDSpinZeroCoupl.H2ggcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2ggcoupl)[vv][0] != 0
1052  ||
1053  (selfDSpinZeroCoupl.H2g4g4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2g4g4coupl)[vv][0] != 0
1054  ){
1055  return true;
1056  }
1057  }
1058  }
1059  return false;
1060 }
1061 bool TEvtProb::CheckSelfDCouplings_Hqq(){
1062  if (_hmass>=0. && _hwidth>0.){
1063  for (int vv = 0; vv < SIZE_HQQ; vv++){
1064  if (
1065  (selfDSpinZeroCoupl.Hqqcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hqqcoupl)[vv][0] != 0
1066  ){
1067  return true;
1068  }
1069  }
1070  }
1071  if (_h2mass>=0. && _h2width>0.){
1072  for (int vv = 0; vv < SIZE_HQQ; vv++){
1073  if (
1074  (selfDSpinZeroCoupl.H2qqcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2qqcoupl)[vv][0] != 0
1075  ){
1076  return true;
1077  }
1078  }
1079  }
1080  return false;
1081 }
1082 bool TEvtProb::CheckSelfDCouplings_Htt(){
1083  if (_hmass>=0. && _hwidth>0.){
1084  for (int vv = 0; vv < SIZE_HQQ; vv++){
1085  if (
1086  (selfDSpinZeroCoupl.Httcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Httcoupl)[vv][0] != 0
1087  ||
1088  (selfDSpinZeroCoupl.Ht4t4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Ht4t4coupl)[vv][0] != 0
1089  ){
1090  return true;
1091  }
1092  }
1093  }
1094  if (_h2mass>=0. && _h2width>0.){
1095  for (int vv = 0; vv < SIZE_HQQ; vv++){
1096  if (
1097  (selfDSpinZeroCoupl.H2ttcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2ttcoupl)[vv][0] != 0
1098  ||
1099  (selfDSpinZeroCoupl.H2t4t4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2t4t4coupl)[vv][0] != 0
1100  ){
1101  return true;
1102  }
1103  }
1104  }
1105  return false;
1106 }
1107 bool TEvtProb::CheckSelfDCouplings_Hbb(){
1108  if (_hmass>=0. && _hwidth>0.){
1109  for (int vv = 0; vv < SIZE_HQQ; vv++){
1110  if (
1111  (selfDSpinZeroCoupl.Hbbcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hbbcoupl)[vv][0] != 0
1112  ||
1113  (selfDSpinZeroCoupl.Hb4b4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hb4b4coupl)[vv][0] != 0
1114  ){
1115  return true;
1116  }
1117  }
1118  }
1119  if (_h2mass>=0. && _h2width>0.){
1120  for (int vv = 0; vv < SIZE_HQQ; vv++){
1121  if (
1122  (selfDSpinZeroCoupl.H2bbcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2bbcoupl)[vv][0] != 0
1123  ||
1124  (selfDSpinZeroCoupl.H2b4b4coupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2b4b4coupl)[vv][0] != 0
1125  ){
1126  return true;
1127  }
1128  }
1129  }
1130  return false;
1131 }
1132 bool TEvtProb::CheckSelfDCouplings_HVV(){
1133  if (_hmass>=0. && _hwidth>0.){
1134  for (int vv = 0; vv < SIZE_HVV; vv++){
1135  if (
1136  (selfDSpinZeroCoupl.Hzzcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hzzcoupl)[vv][0] != 0
1137  ||
1138  (selfDSpinZeroCoupl.Hwwcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.Hwwcoupl)[vv][0] != 0
1139  ){
1140  return true;
1141  }
1142  } // No need to check c_q**2. If these are 0, z_q**2 do not have any effect.
1143  }
1144  if (_h2mass>=0. && _h2width>0.){
1145  for (int vv = 0; vv < SIZE_HVV; vv++){
1146  if (
1147  (selfDSpinZeroCoupl.H2zzcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2zzcoupl)[vv][0] != 0
1148  ||
1149  (selfDSpinZeroCoupl.H2wwcoupl)[vv][1] != 0 || (selfDSpinZeroCoupl.H2wwcoupl)[vv][0] != 0
1150  ){
1151  return true;
1152  }
1153  } // No need to check c_q**2. If these are 0, z_q**2 do not have any effect.
1154  }
1155  return false;
1156 }
1157 bool TEvtProb::CheckSelfDCouplings_aTQGC(){
1158  for (int vv = 0; vv < SIZE_ATQGC; vv++){
1159  if (
1160  (selfDaTQGCCoupl.aTQGCcoupl)[vv][1] != 0 || (selfDaTQGCCoupl.aTQGCcoupl)[vv][0] != 0
1161  ){
1162  return true;
1163  }
1164  }
1165  return false;
1166 }
1167 bool TEvtProb::CheckSelfDCouplings_AZff(){
1168  for (int vv = 0; vv < SIZE_AZff; vv++){
1169  if (
1170  (selfDAZffCoupl.AZffcoupl)[vv][1] != 0 || (selfDAZffCoupl.AZffcoupl)[vv][0] != 0
1171  ){
1172  return true;
1173  }
1174  }
1175  return false;
1176 }
TVar::H1plus
@ H1plus
Definition: TVar.hh:148
SpinZeroCouplings::HzzCLambda_qsq
int HzzCLambda_qsq[SIZE_HVV_CQSQ]
Definition: TCouplings.hh:40
TVar::ERROR
@ ERROR
Definition: TVar.hh:49
MelaIO::getRef
MelaIO * getRef()
Definition: MelaIO.h:51
TEvtProb::XsecCalc_TTX
double XsecCalc_TTX(int topProcess, int topDecay)
Definition: TEvtProb.cc:992
TEvtProb::particleList
std::vector< MELAParticle * > particleList
Definition: TEvtProb.hh:165
SpinZeroCouplings::Hwpwpcoupl
double Hwpwpcoupl[SIZE_HVV][2]
Definition: TCouplings.hh:61
TEvtProb::leptonInterf
TVar::LeptonInterference leptonInterf
Definition: TEvtProb.hh:146
TVar::JJVBF
@ JJVBF
Definition: TVar.hh:72
TUtil::ResonancePropagator
double ResonancePropagator(double const &sqrts, TVar::ResonancePropagatorScheme scheme)
Definition: TUtil.cc:8194
TEvtProb::CheckInputPresent
bool CheckInputPresent()
Definition: TEvtProb.cc:396
TEvtProb::RcdME
MelaIO RcdME
Definition: TEvtProb.hh:162
aTQGCCouplings
Definition: TCouplings.hh:131
TUtil::SetJHUGenSpinZeroGGCouplings
void SetJHUGenSpinZeroGGCouplings(double Hggcoupl[SIZE_HGG][2])
Definition: TUtil.cc:4225
TVar::Had_ZH_S
@ Had_ZH_S
Definition: TVar.hh:86
TEvtProb::SetCandidateDecayMode
void SetCandidateDecayMode(TVar::CandidateDecayMode mode)
Definition: TEvtProb.cc:209
TEvtProb::GetHiggsWidthAtPoleMass
double GetHiggsWidthAtPoleMass(double mass)
Definition: TEvtProb.cc:378
MelaIO
Definition: MelaIO.h:8
TVar::LeptonInterference
LeptonInterference
Definition: TVar.hh:105
SIZE_GVV
@ SIZE_GVV
Definition: raw_couplings.txt:158
TEvtProb::ResetRenFacScaleMode
void ResetRenFacScaleMode()
Definition: TEvtProb.cc:318
TVar::H0_Zgs_PS
@ H0_Zgs_PS
Definition: TVar.hh:131
TEvtProb::GetNCandidates
int GetNCandidates()
Definition: TEvtProb.cc:391
energy_
struct @42 energy_
TVar::VerbosityLevel
VerbosityLevel
Definition: TVar.hh:47
TEvtProb::GetSelfDaTQGCCouplings
aTQGCCouplings * GetSelfDaTQGCCouplings()
Definition: TEvtProb.cc:367
TEvtProb::GetSelfDVprimeCouplings
VprimeCouplings * GetSelfDVprimeCouplings()
Definition: TEvtProb.cc:366
AZffCouplings::getRef
AZffCouplings * getRef()
Definition: TCouplings.cc:486
SpinZeroCouplings::H2zzcoupl
double H2zzcoupl[SIZE_HVV][2]
Definition: TCouplings.hh:51
masses_mcfm_
#define masses_mcfm_
Definition: TMCFM.hh:130
gHIGGS_KAPPA_TILDE
@ gHIGGS_KAPPA_TILDE
Definition: raw_couplings.txt:3
gGRAVITON_VV_10
@ gGRAVITON_VV_10
Definition: raw_couplings.txt:144
TEvtProb::verbosity
TVar::VerbosityLevel verbosity
Definition: TEvtProb.hh:145
TUtil::BBHiggsMatEl
double BBHiggsMatEl(const TVar::Process &process, const TVar::Production &production, const TVar::MatrixElement &matrixElement, TVar::event_scales_type *event_scales, MelaIO *RcdME, const double &EBEAM, int botProcess, TVar::VerbosityLevel verbosity)
Definition: TUtil.cc:7827
TEvtProb::candList
std::vector< MELACandidate * > candList
Definition: TEvtProb.hh:166
TEvtProb::XsecCalcXJ
double XsecCalcXJ()
Definition: TEvtProb.cc:885
SIZE_HVV_CQSQ
@ SIZE_HVV_CQSQ
Definition: raw_couplings.txt:74
TVar::JJEWQCD
@ JJEWQCD
Definition: TVar.hh:74
gGRAVITON_GG_3
@ gGRAVITON_GG_3
Definition: raw_couplings.txt:127
TVar::H0minus
@ H0minus
Definition: TVar.hh:128
TEvtProb::CheckSelfDCouplings_Hgg
bool CheckSelfDCouplings_Hgg()
Definition: TEvtProb.cc:1036
TVar::CandidateDecay_ZZ
@ CandidateDecay_ZZ
Definition: TVar.hh:41
MelaIO::reset
void reset()
Definition: MelaIO.h:31
TEvtProb::SetHiggsMass
void SetHiggsMass(double mass, double wHiggs=-1., int whichResonance=-1)
Definition: TEvtProb.cc:221
TUtil::SetJHUGenHiggsMassWidth
void SetJHUGenHiggsMassWidth(double MReso, double GaReso)
Definition: TUtil.cc:3402
VprimeCouplings::M_Wprime
double M_Wprime
Definition: TCouplings.hh:127
TUtil::HJJMatEl
double HJJMatEl(const TVar::Process &process, const TVar::Production &production, const TVar::MatrixElement &matrixElement, TVar::event_scales_type *event_scales, MelaIO *RcdME, const double &EBEAM, TVar::VerbosityLevel verbosity)
Definition: TUtil.cc:5383
TUtil::SetMCFMAZffCouplings
void SetMCFMAZffCouplings(bool useBSM, AZffCouplings const *Zcouplings)
Definition: TUtil.cc:4185
SpinOneCouplings::Zqqcoupl
double Zqqcoupl[SIZE_ZQQ][2]
Definition: TCouplings.hh:82
TVar::event_scales_type
Definition: TVar.hh:259
TVar::INFO
@ INFO
Definition: TVar.hh:50
TVar::Lep_WH_S
@ Lep_WH_S
Definition: TVar.hh:89
SIZE_GGG
@ SIZE_GGG
Definition: raw_couplings.txt:131
TEvtProb::SetRcdCandPtr
void SetRcdCandPtr()
Definition: TEvtProb.cc:305
aTQGCCouplings::aTQGCcoupl
double aTQGCcoupl[SIZE_ATQGC][2]
Definition: TCouplings.hh:143
VprimeCouplings::Ga_Wprime
double Ga_Wprime
Definition: TCouplings.hh:128
anonymous_namespace{TCouplingsBase.hh}::gHIGGS_GG_2
@ gHIGGS_GG_2
Definition: TCouplingsBase.hh:36
VprimeCouplings::M_Zprime
double M_Zprime
Definition: TCouplings.hh:125
aTQGCCouplings::reset
void reset()
Definition: TCouplings.cc:451
AZffCouplings
Definition: TCouplings.hh:146
SpinZeroCouplings::Ht4t4coupl
double Ht4t4coupl[SIZE_HQQ][2]
Definition: TCouplings.hh:33
TEvtProb::GetTopCandidates
std::vector< MELATopCandidate_t * > * GetTopCandidates()
Definition: TEvtProb.cc:392
TVar::ProductionName
TString ProductionName(TVar::Production temp)
Definition: TVar.cc:64
TEvtProb::ResetCouplings
void ResetCouplings()
Definition: TEvtProb.cc:335
MELAHXSWidth
Definition: MELAHXSWidth.h:10
TEvtProb::GetSelfDSpinOneCouplings
SpinOneCouplings * GetSelfDSpinOneCouplings()
Definition: TEvtProb.cc:364
TEvtProb::GetPrimaryWidth
double GetPrimaryWidth(int ipart)
Definition: TEvtProb.cc:374
VprimeCouplings::Ga_Zprime
double Ga_Zprime
Definition: TCouplings.hh:126
TUtil::SetAlphaS
void SetAlphaS(double &Q_ren, double &Q_fac, double multiplier_ren, double multiplier_fac, int mynloop, int mynflav, std::string mypartons)
TEvtProb::EBEAM
double EBEAM
Definition: TEvtProb.hh:152
TEvtProb::XsecCalc_XVV
double XsecCalc_XVV()
Definition: TEvtProb.cc:414
TVar::CandidateDecayMode
CandidateDecayMode
Definition: TVar.hh:37
TVar::H2_g10
@ H2_g10
Definition: TVar.hh:160
TEvtProb::_h2width
double _h2width
Definition: TEvtProb.hh:151
TEvtProb::CheckSelfDCouplings_Htt
bool CheckSelfDCouplings_Htt()
Definition: TEvtProb.cc:1082
TVar::SelfDefine_spin2
@ SelfDefine_spin2
Definition: TVar.hh:181
TEvtProb::GetPrimaryHiggsMass
double GetPrimaryHiggsMass()
Definition: TEvtProb.cc:369
SIZE_ZQQ
@ SIZE_ZQQ
Definition: raw_couplings.txt:107
TVar::bkgZJets
@ bkgZJets
Definition: TVar.hh:164
TEvtProb::ResetMCFM_EWKParameters
void ResetMCFM_EWKParameters(double ext_Gf, double ext_aemmz, double ext_mW, double ext_mZ, double ext_xW, int ext_ewscheme=3)
Definition: TEvtProb.cc:334
TEvtProb::ResetInputEvent
void ResetInputEvent()
Definition: TEvtProb.cc:345
TEvtProb::SetWprimeMassWidth
void SetWprimeMassWidth(double inmass, double inwidth)
Definition: TEvtProb.cc:321
TVar::ZZINDEPENDENT
@ ZZINDEPENDENT
Definition: TVar.hh:65
TUtil::SetMass
void SetMass(double inmass, int ipart)
Definition: TUtil.cc:1447
PDGHelpers::isAHiggs
bool isAHiggs(const int id)
Definition: PDGHelpers.cc:84
TVar::Process
Process
Definition: TVar.hh:124
TVar::EventScaleScheme
EventScaleScheme
Definition: TVar.hh:195
TVar::JJQCD
@ JJQCD
Definition: TVar.hh:71
TUtil::SetJHUGenVprimeContactCouplings
void SetJHUGenVprimeContactCouplings(double Zpffcoupl[SIZE_Vpff][2], double Wpffcoupl[SIZE_Vpff][2])
Definition: TUtil.cc:4231
TEvtProb::_h2mass
double _h2mass
Definition: TEvtProb.hh:150
SpinZeroCouplings::H2b4b4coupl
double H2b4b4coupl[SIZE_HQQ][2]
Definition: TCouplings.hh:49
TEvtProb::PDFMember_
int PDFMember_
Definition: TEvtProb.hh:141
TEvtProb::TEvtProb
TEvtProb()
Definition: TEvtProb.hh:53
TVar::ZZQQB
@ ZZQQB
Definition: TVar.hh:63
SpinZeroCouplings::HwwLambda_qsq
double HwwLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]
Definition: TCouplings.hh:39
TVar::H1minus
@ H1minus
Definition: TVar.hh:147
VprimeCouplings::Zpffcoupl
double Zpffcoupl[SIZE_Vpff][2]
Definition: TCouplings.hh:122
gHIGGS_AA_2
@ gHIGGS_AA_2
Definition: raw_couplings.txt:18
TEvtProb::InitializeMCFM
void InitializeMCFM()
Definition: TEvtProb.cc:83
SpinTwoCouplings::Gvpvpcoupl
double Gvpvpcoupl[SIZE_GVV][2]
Definition: TCouplings.hh:103
TUtil::SetMCFMSpinZeroCouplings
void SetMCFMSpinZeroCouplings(bool useBSM, SpinZeroCouplings const *Hcouplings, bool forceZZ)
Definition: TUtil.cc:3501
spinzerohiggs_anomcoupl_
struct @34 spinzerohiggs_anomcoupl_
gGRAVITON_VV_1
@ gGRAVITON_VV_1
Definition: raw_couplings.txt:135
TVar::DEBUG_VERBOSE
@ DEBUG_VERBOSE
Definition: TVar.hh:52
nnpdfdriver_
void nnpdfdriver_(char *gridfilename, int *lenfilename)
TUtil::SetJHUGenSpinOneCouplings
void SetJHUGenSpinOneCouplings(double Zqqcoupl[SIZE_ZQQ][2], double Zvvcoupl[SIZE_ZVV][2])
Definition: TUtil.cc:4227
TEvtProb::InitializeJHUGen
void InitializeJHUGen(const char *pathtoPDFSet, int PDFMember)
Definition: TEvtProb.cc:143
TVar::MADGRAPH
@ MADGRAPH
Definition: TVar.hh:59
TEvtProb::AppendTopCandidate
void AppendTopCandidate(SimpleParticleCollection_t *TopDaughters)
Definition: TEvtProb.cc:294
SpinOneCouplings::Zvvcoupl
double Zvvcoupl[SIZE_ZVV][2]
Definition: TCouplings.hh:81
TVar::H2_g7
@ H2_g7
Definition: TVar.hh:157
TEvtProb::CheckSelfDCouplings_Hbb
bool CheckSelfDCouplings_Hbb()
Definition: TEvtProb.cc:1107
SIZE_ZVV
@ SIZE_ZVV
Definition: raw_couplings.txt:114
TVar::ZZQQB_TU
@ ZZQQB_TU
Definition: TVar.hh:92
SpinZeroCouplings::separateWWZZcouplings
bool separateWWZZcouplings
Definition: TCouplings.hh:65
SpinZeroCouplings::Httcoupl
double Httcoupl[SIZE_HQQ][2]
Definition: TCouplings.hh:30
TEvtProb.hh
gHIGGS_ZA_2
@ gHIGGS_ZA_2
Definition: raw_couplings.txt:14
TUtil::SetDecayWidth
void SetDecayWidth(double inwidth, int ipart)
Definition: TUtil.cc:1497
TEvtProb::selfDSpinZeroCoupl
SpinZeroCouplings selfDSpinZeroCoupl
Definition: TEvtProb.hh:156
SpinZeroCouplings::Hzzpcoupl
double Hzzpcoupl[SIZE_HVV][2]
Definition: TCouplings.hh:58
TEvtProb::selfDVprimeCoupl
VprimeCouplings selfDVprimeCoupl
Definition: TEvtProb.hh:159
TVar::ZZGG
@ ZZGG
Definition: TVar.hh:62
TVar::ProcessName
TString ProcessName(TVar::Process temp)
Definition: TVar.cc:6
gGRAVITON_VV_9
@ gGRAVITON_VV_9
Definition: raw_couplings.txt:143
TVar::H0_Zgsg1prime2
@ H0_Zgsg1prime2
Definition: TVar.hh:129
TEvtProb::GetRenFacScaleMode
const TVar::event_scales_type & GetRenFacScaleMode() const
Definition: TEvtProb.cc:216
gGRAVITON_QQ_LEFT
@ gGRAVITON_QQ_LEFT
Definition: raw_couplings.txt:118
TEvtProb::event_scales
TVar::event_scales_type event_scales
Definition: TEvtProb.hh:154
TEvtProb::SetRenFacScaleMode
void SetRenFacScaleMode(TVar::EventScaleScheme renormalizationSch, TVar::EventScaleScheme factorizationSch, double ren_sf, double fac_sf)
Definition: TEvtProb.cc:210
madMela::widths_
struct madMela::@1 widths_
MELAStreamHelpers::MELAout
MELAOutputStreamer MELAout
SpinZeroCouplings::Hg4g4coupl
double Hg4g4coupl[SIZE_HGG][2]
Definition: TCouplings.hh:32
qlinit_
void qlinit_()
TEvtProb::myCSW_
MELAHXSWidth myCSW_
Definition: TEvtProb.hh:153
TVar::DefaultLeptonInterf
@ DefaultLeptonInterf
Definition: TVar.hh:106
TEvtProb::SetLeptonInterf
void SetLeptonInterf(TVar::LeptonInterference tmp)
Definition: TEvtProb.cc:208
MELAParticle::m
double m() const
Definition: MELAParticle.h:66
TVar::H2_g5
@ H2_g5
Definition: TVar.hh:154
TEvtProb::GetHXSWidthEstimator
MELAHXSWidth const * GetHXSWidthEstimator() const
Definition: TEvtProb.cc:362
gGRAVITON_GG_5
@ gGRAVITON_GG_5
Definition: raw_couplings.txt:129
TUtil::SetJHUGenSpinTwoCouplings
void SetJHUGenSpinTwoCouplings(double Gacoupl[SIZE_GGG][2], double Gvvcoupl[SIZE_GVV][2], double Gvvpcoupl[SIZE_GVV][2], double Gvpvpcoupl[SIZE_GVV][2], double qLeftRightcoupl[SIZE_GQQ][2])
Definition: TUtil.cc:4228
SIZE_HQQ
@ SIZE_HQQ
Definition: raw_couplings.txt:5
TEvtProb::production
TVar::Production production
Definition: TEvtProb.hh:144
TVar::H0_g1prime2
@ H0_g1prime2
Definition: TVar.hh:126
TUtil::ConvertVectorFormat
MELACandidate * ConvertVectorFormat(SimpleParticleCollection_t *pDaughters, SimpleParticleCollection_t *pAssociated, SimpleParticleCollection_t *pMothers, bool isGen, std::vector< MELAParticle * > *particleList, std::vector< MELACandidate * > *candList)
Definition: TUtil.cc:8675
TEvtProb::SetCurrentCandidateFromIndex
void SetCurrentCandidateFromIndex(unsigned int icand)
Definition: TEvtProb.cc:306
TEvtProb::melaCand
MELACandidate * melaCand
Definition: TEvtProb.hh:164
TVar::H2_g3
@ H2_g3
Definition: TVar.hh:152
gGRAVITON_VV_2
@ gGRAVITON_VV_2
Definition: raw_couplings.txt:136
qlfirst_
struct @11 qlfirst_
SpinZeroCouplings::Hzpzpcoupl
double Hzpzpcoupl[SIZE_HVV][2]
Definition: TCouplings.hh:59
anonymous_namespace{TCouplingsBase.hh}::SIZE_HGG
@ SIZE_HGG
Definition: TCouplingsBase.hh:40
SpinTwoCouplings::getRef
SpinTwoCouplings * getRef()
Definition: TCouplings.cc:349
TUtil::SetJHUGenSpinZeroVVCouplings
void SetJHUGenSpinZeroVVCouplings(double Hvvcoupl[SIZE_HVV][2], double Hvvpcoupl[SIZE_HVV][2], double Hvpvpcoupl[SIZE_HVV][2], int Hvvcoupl_cqsq[SIZE_HVV_CQSQ], double HvvLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ], bool useWWcoupl)
Definition: TUtil.cc:4219
MELAHXSWidth::HiggsWidth
double HiggsWidth(double mH) const
Definition: MELAHXSWidth.cc:72
TVar::H0_gsgs_PS
@ H0_gsgs_PS
Definition: TVar.hh:133
TEvtProb
Definition: TEvtProb.hh:48
TUtil::SetJHUGenDistinguishWWCouplings
void SetJHUGenDistinguishWWCouplings(bool doAllow)
Definition: TUtil.cc:3408
MadMela.h
This is the code that interfaces with the generated Madgraph matrix elements in fortran for specified...
TEvtProb::_hwidth
double _hwidth
Definition: TEvtProb.hh:149
SpinOneCouplings
Definition: TCouplings.hh:68
TEvtProb::_hmass
double _hmass
Definition: TEvtProb.hh:148
TVar::HSMHiggs
@ HSMHiggs
Definition: TVar.hh:125
VprimeCouplings::getRef
VprimeCouplings * getRef()
Definition: TCouplings.cc:417
TEvtProb::SetInputEvent
void SetInputEvent(SimpleParticleCollection_t *pDaughters, SimpleParticleCollection_t *pAssociated=0, SimpleParticleCollection_t *pMothers=0, bool isGen=false)
Definition: TEvtProb.cc:279
TVar::H2_g2
@ H2_g2
Definition: TVar.hh:151
TUtil::SetJHUGenAZffCouplings
void SetJHUGenAZffCouplings(bool needAZff, double AZffcoupl[SIZE_AZff][2])
Definition: TUtil.cc:4234
TEvtProb::matrixElement
TVar::MatrixElement matrixElement
Definition: TEvtProb.hh:143
TEvtProb::selfDAZffCoupl
AZffCouplings selfDAZffCoupl
Definition: TEvtProb.hh:161
TVar::Lep_WH_TU
@ Lep_WH_TU
Definition: TVar.hh:100
mcfm_init_
void mcfm_init_(char *inputfile, char *workdir)
SpinZeroCouplings::Hbbcoupl
double Hbbcoupl[SIZE_HQQ][2]
Definition: TCouplings.hh:31
AZffCouplings::AZffcoupl
double AZffcoupl[SIZE_AZff][2]
Definition: TCouplings.hh:158
TVar::H2_g9
@ H2_g9
Definition: TVar.hh:159
SpinZeroCouplings::reset
void reset()
Definition: TCouplings.cc:17
TEvtProb::XsecCalc_VVXVV
double XsecCalc_VVXVV()
Definition: TEvtProb.cc:746
TEvtProb::CheckSelfDCouplings_AZff
bool CheckSelfDCouplings_AZff()
Definition: TEvtProb.cc:1167
TEvtProb::GetSelfDSpinTwoCouplings
SpinTwoCouplings * GetSelfDSpinTwoCouplings()
Definition: TEvtProb.cc:365
SpinZeroCouplings::Hqqcoupl
double Hqqcoupl[SIZE_HQQ][2]
Definition: TCouplings.hh:29
TEvtProb::Build
void Build()
Definition: TEvtProb.cc:52
TEvtProb::CheckSelfDCouplings_HVV
bool CheckSelfDCouplings_HVV()
Definition: TEvtProb.cc:1132
TVar::H0_gsgs
@ H0_gsgs
Definition: TVar.hh:132
gHIGGS_VV_1_PRIME2
@ gHIGGS_VV_1_PRIME2
Definition: raw_couplings.txt:23
gHIGGS_VV_4
@ gHIGGS_VV_4
Definition: raw_couplings.txt:12
TUtil::GetDecayWidth
double GetDecayWidth(int ipart)
Definition: TUtil.cc:1619
SpinTwoCouplings::reset
void reset()
Definition: TCouplings.cc:325
aTQGCCouplings::getRef
aTQGCCouplings * getRef()
Definition: TCouplings.cc:461
TEvtProb::selfDSpinOneCoupl
SpinOneCouplings selfDSpinOneCoupl
Definition: TEvtProb.hh:157
gGRAVITON_GG_4
@ gGRAVITON_GG_4
Definition: raw_couplings.txt:128
TEvtProb::CrossInitialize
void CrossInitialize()
Definition: TEvtProb.cc:150
TEvtProb::XsecCalcXJJ
double XsecCalcXJJ()
Definition: TEvtProb.cc:796
gGRAVITON_VV_4
@ gGRAVITON_VV_4
Definition: raw_couplings.txt:138
TVar::Had_ZH_TU
@ Had_ZH_TU
Definition: TVar.hh:97
gGRAVITON_VV_5
@ gGRAVITON_VV_5
Definition: raw_couplings.txt:139
TVar::MCFM
@ MCFM
Definition: TVar.hh:56
TVar::Lep_ZH_S
@ Lep_ZH_S
Definition: TVar.hh:88
SIZE_Vpff
@ SIZE_Vpff
Definition: raw_couplings.txt:100
MELAParticle
Definition: MELAParticle.h:13
SpinZeroCouplings::Hzzcoupl
double Hzzcoupl[SIZE_HVV][2]
Definition: TCouplings.hh:36
TEvtProb::GetIORecord
MelaIO * GetIORecord()
Definition: TEvtProb.cc:382
TVar::JJEW_S
@ JJEW_S
Definition: TVar.hh:84
SpinZeroCouplings::H2wwcoupl
double H2wwcoupl[SIZE_HVV][2]
Definition: TCouplings.hh:52
TUtil::SetEwkCouplingParameters
void SetEwkCouplingParameters(double ext_Gf, double ext_aemmz, double ext_mW, double ext_mZ, double ext_xW, int ext_ewscheme)
Definition: TUtil.cc:1388
VprimeCouplings::Wpffcoupl
double Wpffcoupl[SIZE_Vpff][2]
Definition: TCouplings.hh:123
TEvtProb::ResetMass
void ResetMass(double inmass, int ipart)
Definition: TEvtProb.cc:319
TEvtProb::SetProcess
void SetProcess(TVar::Process proc, TVar::MatrixElement me, TVar::Production prod)
Definition: TEvtProb.cc:193
SpinZeroCouplings::Hwwpcoupl
double Hwwpcoupl[SIZE_HVV][2]
Definition: TCouplings.hh:60
TEvtProb::selfDSpinTwoCoupl
SpinTwoCouplings selfDSpinTwoCoupl
Definition: TEvtProb.hh:158
PDGHelpers::setCandidateDecayMode
void setCandidateDecayMode(TVar::CandidateDecayMode mode)
Definition: PDGHelpers.cc:232
MELACandidate::addAssociatedTop
void addAssociatedTop(MELATopCandidate_t *myParticle)
Definition: MELACandidate.cc:619
TVar::MatrixElement
MatrixElement
Definition: TVar.hh:55
TVar::event_scales_type::factorizationScheme
TVar::EventScaleScheme factorizationScheme
Definition: TVar.hh:261
SIZE_ATQGC
@ SIZE_ATQGC
Definition: raw_couplings.txt:176
TUtil::SetMadgraphSpinZeroCouplings
void SetMadgraphSpinZeroCouplings(SpinZeroCouplings const *Hcouplings)
Definition: TUtil.cc:3415
SpinZeroCouplings::allow_WWZZSeparation
void allow_WWZZSeparation(bool doAllow=true)
Definition: TCouplings.cc:16
SpinZeroCouplings::HzzLambda_qsq
double HzzLambda_qsq[SIZE_HVV_LAMBDAQSQ][SIZE_HVV_CQSQ]
Definition: TCouplings.hh:38
SpinZeroCouplings::getRef
SpinZeroCouplings * getRef()
Definition: TCouplings.cc:106
TVar::ZZQQB_STU
@ ZZQQB_STU
Definition: TVar.hh:64
TEvtProb::pathtoPDFSet_
const char * pathtoPDFSet_
Definition: TEvtProb.hh:140
TEvtProb::GetPrimaryMass
double GetPrimaryMass(int ipart)
Definition: TEvtProb.cc:370
TVar::JJVBF_TU
@ JJVBF_TU
Definition: TVar.hh:94
TVar::event_scales_type::fac_scale_factor
double fac_scale_factor
Definition: TVar.hh:263
TVar::H2_g1
@ H2_g1
Definition: TVar.hh:150
TEvtProb::GetCurrentCandidateIndex
int GetCurrentCandidateIndex()
Definition: TEvtProb.cc:384
SpinZeroCouplings::Hwwcoupl
double Hwwcoupl[SIZE_HVV][2]
Definition: TCouplings.hh:37
TEvtProb::SetVerbosity
void SetVerbosity(TVar::VerbosityLevel tmp)
Definition: TEvtProb.cc:207
gZPRIME_QQ_RIGHT
@ gZPRIME_QQ_RIGHT
Definition: raw_couplings.txt:105
TVar::JHUGen
@ JHUGen
Definition: TVar.hh:57
gHIGGS_AA_4
@ gHIGGS_AA_4
Definition: raw_couplings.txt:20
gGRAVITON_VV_3
@ gGRAVITON_VV_3
Definition: raw_couplings.txt:137
TEvtProb::Set_LHAgrid
void Set_LHAgrid(const char *path, int pdfmember=0)
Definition: TEvtProb.cc:186
SIZE_HVV_LAMBDAQSQ
@ SIZE_HVV_LAMBDAQSQ
Definition: raw_couplings.txt:66
TVar::Had_ZH
@ Had_ZH
Definition: TVar.hh:75
TVar::H2_g1g5
@ H2_g1g5
Definition: TVar.hh:155
TVar::H2_g8
@ H2_g8
Definition: TVar.hh:158
TUtil::ConvertThreeBodyDecayCandidate
MELAThreeBodyDecayCandidate * ConvertThreeBodyDecayCandidate(SimpleParticleCollection_t *tbdDaughters, std::vector< MELAParticle * > *particleList, std::vector< MELAThreeBodyDecayCandidate * > *tbdCandList)
Definition: TUtil.cc:8819
SpinZeroCouplings::H2g4g4coupl
double H2g4g4coupl[SIZE_HGG][2]
Definition: TCouplings.hh:47
TUtil::GetMass
double GetMass(int ipart)
Definition: TUtil.cc:1582
TEvtProb::selfDaTQGCCoupl
aTQGCCouplings selfDaTQGCCoupl
Definition: TEvtProb.hh:160
MELAStreamHelpers::MELAerr
MELAOutputStreamer MELAerr
gZPRIME_QQ_LEFT
@ gZPRIME_QQ_LEFT
Definition: raw_couplings.txt:104
TUtil::MadgraphMatEl
double MadgraphMatEl(const TVar::Process &process, const TVar::Production &production, const TVar::MatrixElement &matrixElement, TVar::event_scales_type *event_scales, MelaIO *RcdME, const double &EBEAM, TVar::VerbosityLevel verbosity)
Definition: TUtil.cc:5258
gZPRIME_VV_2
@ gZPRIME_VV_2
Definition: raw_couplings.txt:112
SpinZeroCouplings::H2bbcoupl
double H2bbcoupl[SIZE_HQQ][2]
Definition: TCouplings.hh:46
TEvtProb::SetZprimeMassWidth
void SetZprimeMassWidth(double inmass, double inwidth)
Definition: TEvtProb.cc:320
gGRAVITON_GG_2
@ gGRAVITON_GG_2
Definition: raw_couplings.txt:126
madMela::mad_masses_
struct madMela::@0 mad_masses_
TVar::ResonancePropagatorScheme
ResonancePropagatorScheme
Definition: TVar.hh:116
TVar::Had_WH_S
@ Had_WH_S
Definition: TVar.hh:87
SpinTwoCouplings::Gvvcoupl
double Gvvcoupl[SIZE_GVV][2]
Definition: TCouplings.hh:101
TVar::JJEW
@ JJEW
Definition: TVar.hh:73
TUtil::SetMCFMaTQGCCouplings
void SetMCFMaTQGCCouplings(bool useBSM, aTQGCCouplings const *couplings)
Definition: TUtil.cc:4134
VprimeCouplings
Definition: TCouplings.hh:108
TEvtProb::GetSelfDSpinZeroCouplings
SpinZeroCouplings * GetSelfDSpinZeroCouplings()
Definition: TEvtProb.cc:363
TEvtProb::GetCurrentCandidate
MELACandidate * GetCurrentCandidate()
Definition: TEvtProb.cc:383
SpinTwoCouplings::Gqqcoupl
double Gqqcoupl[SIZE_GQQ][2]
Definition: TCouplings.hh:104
SIZE_AZff
@ SIZE_AZff
Definition: raw_couplings.txt:194
gGRAVITON_VV_8
@ gGRAVITON_VV_8
Definition: raw_couplings.txt:142
TUtil::TTHiggsMatEl
double TTHiggsMatEl(const TVar::Process &process, const TVar::Production &production, const TVar::MatrixElement &matrixElement, TVar::event_scales_type *event_scales, MelaIO *RcdME, const double &EBEAM, int topDecay, int topProcess, TVar::VerbosityLevel verbosity)
Definition: TUtil.cc:7452
TEvtProb::ResetWidth
void ResetWidth(double inwidth, int ipart)
Definition: TEvtProb.cc:322
TVar::InterfOn
@ InterfOn
Definition: TVar.hh:107
gGRAVITON_VV_6
@ gGRAVITON_VV_6
Definition: raw_couplings.txt:140
__modparameters_MOD_ckmbare
double __modparameters_MOD_ckmbare(int *id1in, int *id2in)
VprimeCouplings::reset
void reset()
Definition: TCouplings.cc:393
TEvtProb::AllowSeparateWWCouplings
void AllowSeparateWWCouplings(bool doAllow=false)
Definition: TEvtProb.cc:219
gGRAVITON_GG_1
@ gGRAVITON_GG_1
Definition: raw_couplings.txt:125
TVar::Had_WH_TU
@ Had_WH_TU
Definition: TVar.hh:98
SpinTwoCouplings::Gvvpcoupl
double Gvvpcoupl[SIZE_GVV][2]
Definition: TCouplings.hh:102
TVar::event_scales_type::ren_scale_factor
double ren_scale_factor
Definition: TVar.hh:262
TVar::H0_Zgs
@ H0_Zgs
Definition: TVar.hh:130
TVar::ttH
@ ttH
Definition: TVar.hh:67
SpinZeroCouplings
Definition: TCouplings.hh:7
TVar::SelfDefine_spin0
@ SelfDefine_spin0
Definition: TVar.hh:179
TUtil::InitJHUGenMELA
void InitJHUGenMELA(const char *pathtoPDFSet, int PDFMember, double collider_sqrts)
Definition: TUtil.cc:3392
SpinOneCouplings::getRef
SpinOneCouplings * getRef()
Definition: TCouplings.cc:302
TUtil::SumMatrixElementPDF
double SumMatrixElementPDF(const TVar::Process &process, const TVar::Production &production, const TVar::MatrixElement &matrixElement, const TVar::LeptonInterference &leptonInterf, TVar::event_scales_type *event_scales, MelaIO *RcdME, const double &EBEAM, TVar::VerbosityLevel verbosity)
Definition: TUtil.cc:4244
gHIGGS_ZA_4
@ gHIGGS_ZA_4
Definition: raw_couplings.txt:16
TEvtProb::SetPrimaryHiggsMass
void SetPrimaryHiggsMass(double mass)
Definition: TEvtProb.cc:220
TEvtProb::XsecCalc_VX
double XsecCalc_VX(bool includeHiggsDecay)
Definition: TEvtProb.cc:916
TEvtProb::GetSelfDAZffCouplings
AZffCouplings * GetSelfDAZffCouplings()
Definition: TEvtProb.cc:368
anonymous_namespace{TCouplingsBase.hh}::gHIGGS_GG_4
@ gHIGGS_GG_4
Definition: TCouplingsBase.hh:38
SpinZeroCouplings::HwwCLambda_qsq
int HwwCLambda_qsq[SIZE_HVV_CQSQ]
Definition: TCouplings.hh:41
TEvtProb::CheckSelfDCouplings_Hqq
bool CheckSelfDCouplings_Hqq()
Definition: TEvtProb.cc:1061
TVar::H2_g6
@ H2_g6
Definition: TVar.hh:156
TVar::Lep_ZH
@ Lep_ZH
Definition: TVar.hh:77
TEvtProb::CheckSelfDCouplings_aTQGC
bool CheckSelfDCouplings_aTQGC()
Definition: TEvtProb.cc:1157
TVar::event_scales_type::renomalizationScheme
TVar::EventScaleScheme renomalizationScheme
Definition: TVar.hh:260
TEvtProb::~TEvtProb
~TEvtProb()
Definition: TEvtProb.cc:75
gHIGGS_VV_1
@ gHIGGS_VV_1
Definition: raw_couplings.txt:9
TVar::ZZQQB_S
@ ZZQQB_S
Definition: TVar.hh:81
nninitpdf_
void nninitpdf_(int *irep)
cabib_
struct @40 cabib_
MELACandidate
Definition: MELACandidate.h:7
TUtil::SetJHUGenSpinZeroQQCouplings
void SetJHUGenSpinZeroQQCouplings(double Hqqcoupl[SIZE_HQQ][2])
Definition: TUtil.cc:4226
AZffCouplings::reset
void reset()
Definition: TCouplings.cc:476
SIZE_HVV
@ SIZE_HVV
Definition: raw_couplings.txt:57
TVar::DEBUG
@ DEBUG
Definition: TVar.hh:51
TVar::JJEWQCD_S
@ JJEWQCD_S
Definition: TVar.hh:85
MELAStreamHelpers.hh
MELAThreeBodyDecayCandidate
Definition: MELAThreeBodyDecayCandidate.h:7
TVar::H2_g4
@ H2_g4
Definition: TVar.hh:153
TVar::JJEWQCD_TU
@ JJEWQCD_TU
Definition: TVar.hh:96
SimpleParticleCollection_t
std::vector< SimpleParticle_t > SimpleParticleCollection_t
Definition: TVar.hh:25
TUtil
Definition: TUtil.hh:49
TEvtProb::PrimaryHMass
double PrimaryHMass
Definition: TEvtProb.hh:147
gGRAVITON_QQ_RIGHT
@ gGRAVITON_QQ_RIGHT
Definition: raw_couplings.txt:119
SpinTwoCouplings
Definition: TCouplings.hh:85
SpinZeroCouplings::H2ttcoupl
double H2ttcoupl[SIZE_HQQ][2]
Definition: TCouplings.hh:45
TVar::DefaultScaleScheme
@ DefaultScaleScheme
Definition: TVar.hh:196
TEvtProb::topCandList
std::vector< MELATopCandidate_t * > topCandList
Definition: TEvtProb.hh:167
TEvtProb::SetCurrentCandidate
void SetCurrentCandidate(MELACandidate *cand)
Definition: TEvtProb.cc:310
TVar::Had_WH
@ Had_WH
Definition: TVar.hh:76
TVar::JJEW_TU
@ JJEW_TU
Definition: TVar.hh:95
TEvtProb::GetXPropagator
double GetXPropagator(TVar::ResonancePropagatorScheme scheme)
Definition: TEvtProb.cc:1028
MelaIO::melaCand
MELACandidate * melaCand
Definition: MelaIO.h:29
TVar::Production
Production
Definition: TVar.hh:61
hto_masses::me
real *8, parameter me
Definition: CALLING_cpHTO.f:77
SpinZeroCouplings::H2t4t4coupl
double H2t4t4coupl[SIZE_HQQ][2]
Definition: TCouplings.hh:48
TVar::JJVBF_S
@ JJVBF_S
Definition: TVar.hh:83
TVar::Lep_WH
@ Lep_WH
Definition: TVar.hh:78
SpinOneCouplings::reset
void reset()
Definition: TCouplings.cc:285
gHIGGS_VV_2
@ gHIGGS_VV_2
Definition: raw_couplings.txt:10
TVar::SelfDefine_spin1
@ SelfDefine_spin1
Definition: TVar.hh:180
SpinZeroCouplings::Hb4b4coupl
double Hb4b4coupl[SIZE_HQQ][2]
Definition: TCouplings.hh:34
__modjhugenmela_MOD_getpdfconstants
void __modjhugenmela_MOD_getpdfconstants(double *pdfzmass, int *pdfnloops, int *pdfnf)
SpinTwoCouplings::Gggcoupl
double Gggcoupl[SIZE_GGG][2]
Definition: TCouplings.hh:105
SpinZeroCouplings::Hggcoupl
double Hggcoupl[SIZE_HGG][2]
Definition: TCouplings.hh:28
TVar::MatrixElementName
TString MatrixElementName(TVar::MatrixElement temp)
Definition: TVar.cc:109
TVar::H0hplus
@ H0hplus
Definition: TVar.hh:127
gGRAVITON_VV_7
@ gGRAVITON_VV_7
Definition: raw_couplings.txt:141
TUtil::JHUGenMatEl
double JHUGenMatEl(const TVar::Process &process, const TVar::Production &production, const TVar::MatrixElement &matrixElement, TVar::event_scales_type *event_scales, MelaIO *RcdME, const double &EBEAM, TVar::VerbosityLevel verbosity)
Definition: TUtil.cc:4887
gHIGGS_KAPPA
@ gHIGGS_KAPPA
Definition: raw_couplings.txt:2
TVar::bbH
@ bbH
Definition: TVar.hh:68
TUtil::VHiggsMatEl
double VHiggsMatEl(const TVar::Process &process, const TVar::Production &production, const TVar::MatrixElement &matrixElement, TVar::event_scales_type *event_scales, MelaIO *RcdME, const double &EBEAM, bool includeHiggsDecay, TVar::VerbosityLevel verbosity)
Definition: TUtil.cc:6843
SpinZeroCouplings::H2qqcoupl
double H2qqcoupl[SIZE_HQQ][2]
Definition: TCouplings.hh:44
TVar::JJQCD_S
@ JJQCD_S
Definition: TVar.hh:82
SIZE_GQQ
@ SIZE_GQQ
Definition: raw_couplings.txt:121
TEvtProb::ResetIORecord
void ResetIORecord()
Definition: TEvtProb.cc:317
TEvtProb::ResetQuarkMasses
void ResetQuarkMasses()
Definition: TEvtProb.cc:323
gZPRIME_VV_1
@ gZPRIME_VV_1
Definition: raw_couplings.txt:111
TVar::Lep_ZH_TU
@ Lep_ZH_TU
Definition: TVar.hh:99
SpinZeroCouplings::H2ggcoupl
double H2ggcoupl[SIZE_HGG][2]
Definition: TCouplings.hh:43
gHIGGS_ZA_1_PRIME2
@ gHIGGS_ZA_1_PRIME2
Definition: raw_couplings.txt:46
TEvtProb::process
TVar::Process process
Definition: TEvtProb.hh:142
TUtil::ResetAmplitudeIncludes
void ResetAmplitudeIncludes()
Definition: TUtil.cc:3412