a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
sigma_tt_diff_mtt_LO_CMS_181106625 Class Reference

Differential cross section of pp->ttbar w.r.t. the invariant mass of the ttbar pair with the binning of 1811.06625. More...

#include <TopQuarkObservables.h>

+ Inheritance diagram for sigma_tt_diff_mtt_LO_CMS_181106625:

Detailed Description

Differential cross section of pp->ttbar w.r.t. the invariant mass of the ttbar pair with the binning of 1811.06625.

Definition at line 119 of file TopQuarkObservables.h.

Public Member Functions

double computeThValue ()
 
double ewgc (const std::string name) const
 
double ewgc (const std::string name, int i, int j) const
 
double ewgc (const std::string name, int i, int j, int k, int l) const
 
 sigma_tt_diff_mtt_LO_CMS_181106625 (const StandardModel &SM_i)
 FL constructor. More...
 
- Public Member Functions inherited from ThObservable
double getBinMax ()
 A get method to get the maximum value of the bin. More...
 
double getBinMin ()
 A get method to get the minimum value of the bin. More...
 
const StandardModelgetModel ()
 A get method to get the model. More...
 
const std::vector< std::string > getParametersForObservable ()
 A get method to get the parameters for the specific observable. More...
 
void setBinMax (double max)
 A set method to set the maximum value of the bin. More...
 
void setBinMin (double min)
 A set method to set the minimum value of the bin. More...
 
void setParametersForObservable (std::vector< std::string > parametersForObservable_i)
 A set method to get the parameters for the specific observable. More...
 
 ThObservable (const StandardModel &SM_i)
 Constructor. More...
 
 ThObservable (const ThObservable &orig)
 The copy constructor. More...
 
virtual ~ThObservable ()
 The default destructor. More...
 

Private Attributes

double b_max = 0.
 
double b_min = 0.
 
const TopQuarkObservables mytopobs
 

Additional Inherited Members

- Protected Attributes inherited from ThObservable
double max
 the bin maximum. More...
 
double min
 The bin minimum. More...
 
std::vector< std::string > parametersForObservable
 a vector of parameter namesfor the specific observable More...
 
const StandardModelSM
 A reference to an object of StandardMode class. More...
 

Constructor & Destructor Documentation

◆ sigma_tt_diff_mtt_LO_CMS_181106625()

sigma_tt_diff_mtt_LO_CMS_181106625::sigma_tt_diff_mtt_LO_CMS_181106625 ( const StandardModel SM_i)

FL constructor.

Definition at line 284 of file TopQuarkObservables.cpp.

285: ThObservable(SM_i), mytopobs(static_cast<const NPSMEFTd6General&> (SM))
286{
287
288 setParametersForObservable(make_vector<std::string>() << "SM_sigma_tt_bin_mtt_300_380_CMS_181106625"
289 << "SM_sigma_tt_bin_mtt_380_470_CMS_181106625" << "SM_sigma_tt_bin_mtt_470_620_CMS_181106625"
290 << "SM_sigma_tt_bin_mtt_620_820_CMS_181106625" << "SM_sigma_tt_bin_mtt_820_1100_CMS_181106625"
291 << "SM_sigma_tt_bin_mtt_1100_1500_CMS_181106625" << "SM_sigma_tt_bin_mtt_1500_2500_CMS_181106625");
292
293
294}
A model class for new physics in the form of the dimension-six effective Lagrangian.
void setParametersForObservable(std::vector< std::string > parametersForObservable_i)
A set method to get the parameters for the specific observable.
Definition: ThObservable.h:109
ThObservable(const StandardModel &SM_i)
Constructor.
Definition: ThObservable.h:32
const StandardModel & SM
A reference to an object of StandardMode class.
Definition: ThObservable.h:121
const TopQuarkObservables mytopobs

Member Function Documentation

◆ computeThValue()

double sigma_tt_diff_mtt_LO_CMS_181106625::computeThValue ( )
virtual
Returns
The value of Test_direct

Implements ThObservable.

Definition at line 296 of file TopQuarkObservables.cpp.

297{
298
299 b_min = getBinMin();
300 b_max = getBinMax();
301
302 bool flag_Quadratic= false; //Needs to be properly defined
303
304
305
306 if(b_min == 300 && b_max == 380){
307
308 double SM_sigma_tt_bin_300_380 = SM.getOptionalParameter("SM_sigma_tt_bin_mtt_300_380_CMS_181106625");
309 double sigma_tt_bin_300_380_madgraph = 54.003200;//pb
310
311 if(flag_Quadratic){
312
313 return 0.;
314
315 }
316 else{
317
318 double total = SM_sigma_tt_bin_300_380 + (+-1.080720*ewgc("CG")+0.015822*ewgc("Cqd1R",2,2,0,0)
319 +0.027498*ewgc("Cqd1R",2,2,1,1)+0.239796*ewgc("Cqd8R",2,2,0,0)+0.073740*ewgc("Cqd8R",2,2,1,1)
320 +0.089304*ewgc("Cqq1R",0,0,2,2)+1.256970*ewgc("Cqq1R",0,2,2,0)+0.033486*ewgc("Cqq1R",1,1,2,2)
321 +0.108318*ewgc("Cqq1R",1,2,2,1)+0.230184*ewgc("Cqq3R",0,0,2,2)+2.777490*ewgc("Cqq3R",0,2,2,0)
322 +0.067932*ewgc("Cqq3R",1,1,2,2)+0.457110*ewgc("Cqq3R",1,2,2,1)+0.071850*ewgc("Cqu1R",0,0,2,2)
323 +0.029322*ewgc("Cqu1R",1,1,2,2)+0.080706*ewgc("Cqu1R",2,2,0,0)+0.034674*ewgc("Cqu1R",2,2,1,1)
324 +0.535692*ewgc("Cqu8R",0,0,2,2)+0.091128*ewgc("Cqu8R",1,1,2,2)+0.347364*ewgc("Cqu8R",2,2,0,0)
325 +0.057720*ewgc("Cqu8R",2,2,1,1)+0.014940*ewgc("Cud1R",2,2,0,0)+0.030534*ewgc("Cud1R",2,2,1,1)
326 +0.241590*ewgc("Cud8R",2,2,0,0)+0.076716*ewgc("Cud8R",2,2,1,1)+-18.341500*ewgc("CuGR",2,2)
327 +0.081258*ewgc("CuuR",0,0,2,2)+1.238740*ewgc("CuuR",0,2,2,0)+0.036636*ewgc("CuuR",1,1,2,2)
328 +0.108228*ewgc("CuuR",1,2,2,1)
329 )*(SM_sigma_tt_bin_300_380/sigma_tt_bin_300_380_madgraph);
330
331 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
332
333 return total;
334
335
336 }
337
338 } else if(b_min == 380 && b_max == 470){
339
340 double SM_sigma_tt_bin_380_470 = SM.getOptionalParameter("SM_sigma_tt_bin_mtt_380_470_CMS_181106625");
341 double sigma_tt_bin_380_470_madgraph = 173.423000;//pb
342
343 if(flag_Quadratic){
344
345 return 0.;
346
347 }
348 else{
349
350 double total = SM_sigma_tt_bin_380_470 + (+-8.359080*ewgc("CG")+-0.042078*ewgc("Cqd1R",2,2,0,0)+0.004590*ewgc("Cqd1R",2,2,1,1)
351 +0.649842*ewgc("Cqd8R",2,2,0,0)+0.144468*ewgc("Cqd8R",2,2,1,1)+0.260772*ewgc("Cqq1R",0,0,2,2)+4.262270*ewgc("Cqq1R",0,2,2,0)
352 +-0.029244*ewgc("Cqq1R",1,1,2,2)+0.230286*ewgc("Cqq1R",1,2,2,1)+0.780282*ewgc("Cqq3R",0,0,2,2)+9.255620*ewgc("Cqq3R",0,2,2,0)
353 +0.084678*ewgc("Cqq3R",1,1,2,2)+1.257790*ewgc("Cqq3R",1,2,2,1)+0.071382*ewgc("Cqu1R",0,0,2,2)+0.007008*ewgc("Cqu1R",1,1,2,2)
354 +0.097896*ewgc("Cqu1R",2,2,0,0)+0.018918*ewgc("Cqu1R",2,2,1,1)+1.653840*ewgc("Cqu8R",0,0,2,2)+0.194502*ewgc("Cqu8R",1,1,2,2)
355 +1.000630*ewgc("Cqu8R",2,2,0,0)+0.057138*ewgc("Cqu8R",2,2,1,1)+-0.038910*ewgc("Cud1R",2,2,0,0)+0.004434*ewgc("Cud1R",2,2,1,1)
356 +0.648138*ewgc("Cud8R",2,2,0,0)+0.144918*ewgc("Cud8R",2,2,1,1)+-51.556900*ewgc("CuGR",2,2)+0.308838*ewgc("CuuR",0,0,2,2)
357 +4.176160*ewgc("CuuR",0,2,2,0)+0.030888*ewgc("CuuR",1,1,2,2)+0.226812*ewgc("CuuR",1,2,2,1)
358 )*(SM_sigma_tt_bin_380_470/sigma_tt_bin_380_470_madgraph);
359
360 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
361
362 return total;
363
364 }
365
366 } else if(b_min == 470 && b_max == 620){
367
368 double SM_sigma_tt_bin_470_620 = SM.getOptionalParameter("SM_sigma_tt_bin_mtt_470_620_CMS_181106625");
369 double sigma_tt_bin_470_620_madgraph = 152.893000;//pb
370
371 if(flag_Quadratic){
372
373 return 0.;
374
375 }
376 else{
377
378 double total = SM_sigma_tt_bin_470_620 + (+-12.015500*ewgc("CG")+-0.014436*ewgc("Cqd1R",2,2,0,0)+-0.010944*ewgc("Cqd1R",2,2,1,1)
379 +0.738569*ewgc("Cqd8R",2,2,0,0)+0.127001*ewgc("Cqd8R",2,2,1,1)+0.294306*ewgc("Cqq1R",0,0,2,2)+4.783330*ewgc("Cqq1R",0,2,2,0)
380 +-0.041070*ewgc("Cqq1R",1,1,2,2)+0.199307*ewgc("Cqq1R",1,2,2,1)+1.030540*ewgc("Cqq3R",0,0,2,2)+10.384200*ewgc("Cqq3R",0,2,2,0)
381 +0.084312*ewgc("Cqq3R",1,1,2,2)+1.174270*ewgc("Cqq3R",1,2,2,1)+0.134514*ewgc("Cqu1R",0,0,2,2)+-0.007254*ewgc("Cqu1R",1,1,2,2)
382 +0.095286*ewgc("Cqu1R",2,2,0,0)+0.000101*ewgc("Cqu1R",2,2,1,1)+1.890170*ewgc("Cqu8R",0,0,2,2)+0.177748*ewgc("Cqu8R",1,1,2,2)
383 +1.147850*ewgc("Cqu8R",2,2,0,0)+0.048700*ewgc("Cqu8R",2,2,1,1)+-0.042132*ewgc("Cud1R",2,2,0,0)+-0.013392*ewgc("Cud1R",2,2,1,1)
384 +0.739715*ewgc("Cud8R",2,2,0,0)+0.126949*ewgc("Cud8R",2,2,1,1)+-42.458000*ewgc("CuGR",2,2)+0.326724*ewgc("CuuR",0,0,2,2)
385 +4.680410*ewgc("CuuR",0,2,2,0)+0.007338*ewgc("CuuR",1,1,2,2)+0.194285*ewgc("CuuR",1,2,2,1)
386 )*(SM_sigma_tt_bin_470_620/sigma_tt_bin_470_620_madgraph);
387
388 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
389
390 return total;
391
392 }
393
394 } else if(b_min == 620 && b_max == 820){
395
396 double SM_sigma_tt_bin_620_820 = SM.getOptionalParameter("SM_sigma_tt_bin_mtt_620_820_CMS_181106625");
397 double sigma_tt_bin_620_820_madgraph = 66.954300;//pb
398
399 if(flag_Quadratic){
400
401 return 0.;
402
403 }
404 else{
405
406 double total = SM_sigma_tt_bin_620_820 + (+-6.708550*ewgc("CG")+-0.017060*ewgc("Cqd1R",2,2,0,0)+0.003028*ewgc("Cqd1R",2,2,1,1)
407 +0.511089*ewgc("Cqd8R",2,2,0,0)+0.076736*ewgc("Cqd8R",2,2,1,1)+0.224652*ewgc("Cqq1R",0,0,2,2)+3.446780*ewgc("Cqq1R",0,2,2,0)
408 +-0.023712*ewgc("Cqq1R",1,1,2,2)+0.115517*ewgc("Cqq1R",1,2,2,1)+0.804056*ewgc("Cqq3R",0,0,2,2)+7.363740*ewgc("Cqq3R",0,2,2,0)
409 +0.059872*ewgc("Cqq3R",1,1,2,2)+0.658982*ewgc("Cqq3R",1,2,2,1)+0.068678*ewgc("Cqu1R",0,0,2,2)+0.006047*ewgc("Cqu1R",1,1,2,2)
410 +0.049867*ewgc("Cqu1R",2,2,0,0)+0.005234*ewgc("Cqu1R",2,2,1,1)+1.331380*ewgc("Cqu8R",0,0,2,2)+0.102676*ewgc("Cqu8R",1,1,2,2)
411 +0.823487*ewgc("Cqu8R",2,2,0,0)+0.029953*ewgc("Cqu8R",2,2,1,1)+-0.042581*ewgc("Cud1R",2,2,0,0)+-0.001945*ewgc("Cud1R",2,2,1,1)
412 +0.511913*ewgc("Cud8R",2,2,0,0)+0.077032*ewgc("Cud8R",2,2,1,1)+-18.574200*ewgc("CuGR",2,2)+0.262504*ewgc("CuuR",0,0,2,2)
413 +3.359280*ewgc("CuuR",0,2,2,0)+0.010662*ewgc("CuuR",1,1,2,2)+0.112864*ewgc("CuuR",1,2,2,1)
414 )*(SM_sigma_tt_bin_620_820/sigma_tt_bin_620_820_madgraph);
415
416 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
417
418 return total;
419
420 }
421
422 } else if(b_min == 820 && b_max == 1100){
423
424 double SM_sigma_tt_bin_820_1100 = SM.getOptionalParameter("SM_sigma_tt_bin_mtt_820_1100_CMS_181106625");
425 double sigma_tt_bin_820_1100_madgraph = 24.163600;//pb
426
427 if(flag_Quadratic){
428
429 return 0.;
430
431 }
432 else{
433
434 double total = SM_sigma_tt_bin_820_1100 + (-2.802570*ewgc("CG")+-0.006901*ewgc("Cqd1R",2,2,0,0)+0.000237*ewgc("Cqd1R",2,2,1,1)
435 +0.333880*ewgc("Cqd8R",2,2,0,0)+0.037791*ewgc("Cqd8R",2,2,1,1)+0.161210*ewgc("Cqq1R",0,0,2,2)+2.286240*ewgc("Cqq1R",0,2,2,0)
436 +-0.012073*ewgc("Cqq1R",1,1,2,2)+0.059007*ewgc("Cqq1R",1,2,2,1)+0.566288*ewgc("Cqq3R",0,0,2,2)+4.802780*ewgc("Cqq3R",0,2,2,0)
437 +0.033808*ewgc("Cqq3R",1,1,2,2)+0.333648*ewgc("Cqq3R",1,2,2,1)+0.042755*ewgc("Cqu1R",0,0,2,2)+0.002326*ewgc("Cqu1R",1,1,2,2)
438 +0.035514*ewgc("Cqu1R",2,2,0,0)+0.001443*ewgc("Cqu1R",2,2,1,1)+0.875141*ewgc("Cqu8R",0,0,2,2)+0.051098*ewgc("Cqu8R",1,1,2,2)
439 +0.542868*ewgc("Cqu8R",2,2,0,0)+0.013234*ewgc("Cqu8R",2,2,1,1)+-0.023700*ewgc("Cud1R",2,2,0,0)+-0.001051*ewgc("Cud1R",2,2,1,1)
440 +0.332700*ewgc("Cud8R",2,2,0,0)+0.037501*ewgc("Cud8R",2,2,1,1)-6.912790*ewgc("CuGR",2,2)+0.178817*ewgc("CuuR",0,0,2,2)
441 +2.227770*ewgc("CuuR",0,2,2,0)+0.006347*ewgc("CuuR",1,1,2,2)+0.057377*ewgc("CuuR",1,2,2,1)
442 )*(SM_sigma_tt_bin_820_1100/sigma_tt_bin_820_1100_madgraph);
443
444 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
445
446 return total;
447
448 }
449
450 } else if(b_min == 1100 && b_max == 1500){
451
452 double SM_sigma_tt_bin_1100_1500 = SM.getOptionalParameter("SM_sigma_tt_bin_mtt_1100_1500_CMS_181106625");
453 double sigma_tt_bin_1100_1500_madgraph = 6.983800;//pb
454
455 if(flag_Quadratic){
456
457 return 0.;
458
459 }
460 else{
461
462 double total = SM_sigma_tt_bin_1100_1500 + (+-0.863200*ewgc("CG")+-0.004510*ewgc("Cqd1R",2,2,0,0)+0.000022*ewgc("Cqd1R",2,2,1,1)
463 +0.191513*ewgc("Cqd8R",2,2,0,0)+0.016414*ewgc("Cqd8R",2,2,1,1)+0.103366*ewgc("Cqq1R",0,0,2,2)+1.351660*ewgc("Cqq1R",0,2,2,0)
464 +-0.005690*ewgc("Cqq1R",1,1,2,2)+0.026285*ewgc("Cqq1R",1,2,2,1)+0.355618*ewgc("Cqq3R",0,0,2,2)+2.795780*ewgc("Cqq3R",0,2,2,0)
465 +0.014960*ewgc("Cqq3R",1,1,2,2)+0.147160*ewgc("Cqq3R",1,2,2,1)+0.024812*ewgc("Cqu1R",0,0,2,2)+0.001067*ewgc("Cqu1R",1,1,2,2)
466 +0.017626*ewgc("Cqu1R",2,2,0,0)+0.000632*ewgc("Cqu1R",2,2,1,1)+0.509667*ewgc("Cqu8R",0,0,2,2)+0.022537*ewgc("Cqu8R",1,1,2,2)
467 +0.318022*ewgc("Cqu8R",2,2,0,0)+0.006561*ewgc("Cqu8R",2,2,1,1)+-0.016528*ewgc("Cud1R",2,2,0,0)+-0.001045*ewgc("Cud1R",2,2,1,1)
468 +0.191209*ewgc("Cud8R",2,2,0,0)+0.016353*ewgc("Cud8R",2,2,1,1)+-2.105990*ewgc("CuGR",2,2)+0.112485*ewgc("CuuR",0,0,2,2)
469 +1.312690*ewgc("CuuR",0,2,2,0)+0.002397*ewgc("CuuR",1,1,2,2)+0.025537*ewgc("CuuR",1,2,2,1)
470 )*(SM_sigma_tt_bin_1100_1500/sigma_tt_bin_1100_1500_madgraph);
471
472 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
473
474 return total;
475
476 }
477
478 } else if(b_min == 1500 && b_max == 2500){
479
480 double SM_sigma_tt_bin_1500_2500 = SM.getOptionalParameter("SM_sigma_tt_bin_mtt_1500_2500_CMS_181106625");
481 double sigma_tt_bin_1500_2500_madgraph = 1.821290;//pb
482
483 if(flag_Quadratic){
484
485 return 0.;
486
487 }
488 else{
489
490 double total = SM_sigma_tt_bin_1500_2500 + (+-0.231926*ewgc("CG")+-0.003089*ewgc("Cqd1R",2,2,0,0)+-0.000047*ewgc("Cqd1R",2,2,1,1)
491 +0.130652*ewgc("Cqd8R",2,2,0,0)+0.007679*ewgc("Cqd8R",2,2,1,1)+0.078637*ewgc("Cqq1R",0,0,2,2)+0.983935*ewgc("Cqq1R",0,2,2,0)
492 +-0.002654*ewgc("Cqq1R",1,1,2,2)+0.014376*ewgc("Cqq1R",1,2,2,1)+0.252179*ewgc("Cqq3R",0,0,2,2)+1.968040*ewgc("Cqq3R",0,2,2,0)
493 +0.007838*ewgc("Cqq3R",1,1,2,2)+0.071898*ewgc("Cqq3R",1,2,2,1)+0.017883*ewgc("Cqu1R",0,0,2,2)+0.000657*ewgc("Cqu1R",1,1,2,2)
494 +0.011678*ewgc("Cqu1R",2,2,0,0)+0.000359*ewgc("Cqu1R",2,2,1,1)+0.363022*ewgc("Cqu8R",0,0,2,2)+0.011088*ewgc("Cqu8R",1,1,2,2)
495 +0.232613*ewgc("Cqu8R",2,2,0,0)+0.003450*ewgc("Cqu8R",2,2,1,1)+-0.011469*ewgc("Cud1R",2,2,0,0)+-0.000549*ewgc("Cud1R",2,2,1,1)
496 +0.130711*ewgc("Cud8R",2,2,0,0)+0.007725*ewgc("Cud8R",2,2,1,1)+-0.607749*ewgc("CuGR",2,2)+0.081595*ewgc("CuuR",0,0,2,2)
497 +0.955959*ewgc("CuuR",0,2,2,0)+0.001335*ewgc("CuuR",1,1,2,2)+0.014025*ewgc("CuuR",1,2,2,1)
498 )*(SM_sigma_tt_bin_1500_2500/sigma_tt_bin_1500_2500_madgraph);
499
500 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
501
502 return total;
503
504 }
505
506 }
507
508 else {
509 throw std::runtime_error("\nERROR: Please specify a correct binning range for sigma_tt_diff_mtt_LO_CMS_181106625.\n");
510 }
511
512}
const double getOptionalParameter(std::string name) const
A method to get parameters that are specific to only one set of observables.
Definition: QCD.h:450
double getBinMin()
A get method to get the minimum value of the bin.
Definition: ThObservable.h:82
double getBinMax()
A get method to get the maximum value of the bin.
Definition: ThObservable.h:91
double ewgc(const std::string name) const

◆ ewgc() [1/3]

double sigma_tt_diff_mtt_LO_CMS_181106625::ewgc ( const std::string  name) const
inline

Definition at line 136 of file TopQuarkObservables.h.

137 {
138 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
139 }
double getSMEFTCoeffEW(const std::string name) const
const NPSMEFTd6General & GetmyNPSMEFTd6General() const

◆ ewgc() [2/3]

double sigma_tt_diff_mtt_LO_CMS_181106625::ewgc ( const std::string  name,
int  i,
int  j 
) const
inline

Definition at line 141 of file TopQuarkObservables.h.

142 {
143 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
144 }

◆ ewgc() [3/3]

double sigma_tt_diff_mtt_LO_CMS_181106625::ewgc ( const std::string  name,
int  i,
int  j,
int  k,
int  l 
) const
inline

Definition at line 146 of file TopQuarkObservables.h.

147 {
148 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
149 }

Member Data Documentation

◆ b_max

double sigma_tt_diff_mtt_LO_CMS_181106625::b_max = 0.
private

Definition at line 154 of file TopQuarkObservables.h.

◆ b_min

double sigma_tt_diff_mtt_LO_CMS_181106625::b_min = 0.
private

Definition at line 153 of file TopQuarkObservables.h.

◆ mytopobs

const TopQuarkObservables sigma_tt_diff_mtt_LO_CMS_181106625::mytopobs
private

Definition at line 157 of file TopQuarkObservables.h.


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