a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
Public Member Functions | Private Attributes | List of all members
sigma_tta_diff_LO_CMS_semileptonic_210701508 Class Reference

#include <TopQuarkObservables.h>

+ Inheritance diagram for sigma_tta_diff_LO_CMS_semileptonic_210701508:

Detailed Description

Definition at line 676 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_tta_diff_LO_CMS_semileptonic_210701508 (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 StandardModel & getModel ()
 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 StandardModel & SM
 A reference to an object of StandardMode class. More...
 

Constructor & Destructor Documentation

◆ sigma_tta_diff_LO_CMS_semileptonic_210701508()

sigma_tta_diff_LO_CMS_semileptonic_210701508::sigma_tta_diff_LO_CMS_semileptonic_210701508 ( const StandardModel &  SM_i)

FL constructor.

Definition at line 2359 of file TopQuarkObservables.cpp.

2360: ThObservable(SM_i), mytopobs(static_cast<const NPSMEFTd6General&> (SM))
2361{
2362
2363 setParametersForObservable(make_vector<std::string>() << "SM_sigma_tta_bin_20_35_CMS_semileptonic"
2364 << "SM_sigma_tta_bin_35_50_CMS_semileptonic" << "SM_sigma_tta_bin_50_80_CMS_semileptonic"
2365 << "SM_sigma_tta_bin_80_120_CMS_semileptonic" << "SM_sigma_tta_bin_120_160_CMS_semileptonic"
2366 << "SM_sigma_tta_bin_160_200_CMS_semileptonic" << "SM_sigma_tta_bin_200_260_CMS_semileptonic"
2367 << "SM_sigma_tta_bin_260_320_CMS_semileptonic" << "SM_sigma_tta_bin_320_380_CMS_semileptonic");
2368
2369
2370}
NPSMEFTd6General
A model class for new physics in the form of the dimension-six effective Lagrangian.
Definition: NPSMEFTd6General.h:593
ThObservable::setParametersForObservable
void setParametersForObservable(std::vector< std::string > parametersForObservable_i)
A set method to get the parameters for the specific observable.
Definition: ThObservable.h:109
ThObservable::ThObservable
ThObservable(const StandardModel &SM_i)
Constructor.
Definition: ThObservable.h:32
ThObservable::SM
const StandardModel & SM
A reference to an object of StandardMode class.
Definition: ThObservable.h:121
sigma_tta_diff_LO_CMS_semileptonic_210701508::mytopobs
const TopQuarkObservables mytopobs
Definition: TopQuarkObservables.h:713

Member Function Documentation

◆ computeThValue()

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

Implements ThObservable.

Definition at line 2373 of file TopQuarkObservables.cpp.

2374{
2375
2376 b_min = getBinMin();
2377 b_max = getBinMax();
2378
2379 bool flag_Quadratic= false; //Needs to be properly defined
2380
2381
2382
2383
2384
2385 if(b_min == 20 && b_max == 35){
2386
2387 double SM_sigma_tta_bin_20_35 = SM.getOptionalParameter("SM_sigma_tta_bin_20_35_CMS_semileptonic");
2388 double sigma_tta_bin_20_35_madgraph = 0.583675; //pb
2389
2390
2391 if(flag_Quadratic){
2392
2393 return 0.;
2394
2395 }
2396 else{
2397
2398
2399 double total = SM_sigma_tta_bin_20_35 +(-0.034106*ewgc("CG")+0.004794*ewgc("Cqd8R",2,2,0,0)
2400 +0.000605*ewgc("Cqd8R",2,2,1,1)+0.082493*ewgc("Cqq1R",0,2,2,0)
2401 +0.002529*ewgc("Cqq1R",1,2,2,1)+0.146282*ewgc("Cqq3R",0,2,2,0)+0.010782*ewgc("Cqq3R",1,2,2,1)
2402 +0.036042*ewgc("Cqu8R",0,0,2,2)+0.001578*ewgc("Cqu8R",1,1,2,2)+0.031190*ewgc("Cqu8R",2,2,0,0)
2403 +0.000906*ewgc("Cqu8R",2,2,1,1)+-0.007204*ewgc("CuBR",2,2)+0.008516*ewgc("Cud8R",2,2,0,0)
2404 +0.001057*ewgc("Cud8R",2,2,1,1)+-0.154829*ewgc("CuGR",2,2)
2405 +0.080748*ewgc("CuuR",0,2,2,0)+0.002498*ewgc("CuuR",1,2,2,1)+-0.005006*ewgc("CuWR",2,2)
2406 +-0.061811*ewgc("CHD")+-0.035571*ewgc("CHl3R",0,0)+-0.035571*ewgc("CHl3R",1,1)+-0.132541*ewgc("CHWB")
2407 +0.035506*ewgc("CllR",0,1,1,0)+-0.000081*ewgc("Cqd1R",2,2,0,0)+-0.000027*ewgc("Cqd1R",2,2,1,1)
2408 +0.006399*ewgc("Cqq1R",0,0,2,2)+-0.000210*ewgc("Cqq1R",1,1,2,2)+0.016394*ewgc("Cqq3R",0,0,2,2)
2409 +0.000986*ewgc("Cqq3R",1,1,2,2)+0.002978*ewgc("Cqu1R",0,0,2,2)+0.000069*ewgc("Cqu1R",1,1,2,2)
2410 +0.002181*ewgc("Cqu1R",2,2,0,0)+0.000072*ewgc("Cqu1R",2,2,1,1)+-0.000689*ewgc("Cud1R",2,2,0,0)
2411 +-0.000057*ewgc("Cud1R",2,2,1,1)+0.005823*ewgc("CuuR",0,0,2,2)+0.000189*ewgc("CuuR",1,1,2,2))
2412 *(SM_sigma_tta_bin_20_35/sigma_tta_bin_20_35_madgraph);
2413
2414 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2415
2416 return total;
2417
2418 }
2419 } else if(b_min == 35 && b_max == 50){
2420
2421
2422 double SM_sigma_tta_bin_35_50 = SM.getOptionalParameter("SM_sigma_tta_bin_35_50_CMS_semileptonic");
2423 double sigma_tta_bin_35_50_madgraph = 0.305601; //pb
2424
2425
2426 if(flag_Quadratic){
2427
2428 return 0.;
2429
2430 }
2431 else{
2432
2433 double total = SM_sigma_tta_bin_35_50 +(-0.018555*ewgc("CG")+0.002547*ewgc("Cqd8R",2,2,0,0)
2434 +0.000295*ewgc("Cqd8R",2,2,1,1)+0.040351*ewgc("Cqq1R",0,2,2,0)
2435 +0.001256*ewgc("Cqq1R",1,2,2,1)+0.075386*ewgc("Cqq3R",0,2,2,0)+0.005401*ewgc("Cqq3R",1,2,2,1)
2436 +0.018788*ewgc("Cqu8R",0,0,2,2)+0.000807*ewgc("Cqu8R",1,1,2,2)+0.016262*ewgc("Cqu8R",2,2,0,0)
2437 +0.000470*ewgc("Cqu8R",2,2,1,1)+-0.006535*ewgc("CuBR",2,2)+0.004622*ewgc("Cud8R",2,2,0,0)
2438 +0.000557*ewgc("Cud8R",2,2,1,1)-0.080445*ewgc("CuGR",2,2)
2439 +0.039496*ewgc("CuuR",0,2,2,0)+0.001217*ewgc("CuuR",1,2,2,1)+-0.004069*ewgc("CuWR",2,2)
2440 +-0.032332*ewgc("CHD")+-0.018575*ewgc("CHl3R",0,0)+-0.018575*ewgc("CHl3R",1,1)+-0.069327*ewgc("CHWB")
2441 +0.018597*ewgc("CllR",0,1,1,0)+-0.000057*ewgc("Cqd1R",2,2,0,0)+-0.000002*ewgc("Cqd1R",2,2,1,1)
2442 +0.003004*ewgc("Cqq1R",0,0,2,2)+-0.000132*ewgc("Cqq1R",1,1,2,2)+0.008326*ewgc("Cqq3R",0,0,2,2)
2443 +0.000523*ewgc("Cqq3R",1,1,2,2)+0.001386*ewgc("Cqu1R",0,0,2,2)+0.000057*ewgc("Cqu1R",1,1,2,2)
2444 +0.001019*ewgc("Cqu1R",2,2,0,0)+0.000055*ewgc("Cqu1R",2,2,1,1)+-0.000289*ewgc("Cud1R",2,2,0,0)
2445 +-0.000021*ewgc("Cud1R",2,2,1,1)+0.002900*ewgc("CuuR",0,0,2,2)+0.000133*ewgc("CuuR",1,1,2,2))
2446 *(SM_sigma_tta_bin_35_50/sigma_tta_bin_35_50_madgraph);
2447
2448 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2449
2450 return total;
2451
2452 }
2453 } else if(b_min == 50 && b_max == 80){
2454
2455 double SM_sigma_tta_bin_50_80 = SM.getOptionalParameter("SM_sigma_tta_bin_50_80_CMS_semileptonic");
2456 double sigma_tta_bin_50_80_madgraph = 0.319505; //pb
2457
2458
2459 if(flag_Quadratic){
2460
2461 return 0.;
2462
2463 }
2464 else{
2465
2466
2467 double total = SM_sigma_tta_bin_50_80 +(-0.020248*ewgc("CG")+0.002577*ewgc("Cqd8R",2,2,0,0)
2468 +0.000287*ewgc("Cqd8R",2,2,1,1)+0.039600*ewgc("Cqq1R",0,2,2,0)
2469 +0.001092*ewgc("Cqq1R",1,2,2,1)+0.077011*ewgc("Cqq3R",0,2,2,0)
2470 +0.005369*ewgc("Cqq3R",1,2,2,1)+0.019650*ewgc("Cqu8R",0,0,2,2)
2471 +0.000769*ewgc("Cqu8R",1,1,2,2)+0.017049*ewgc("Cqu8R",2,2,0,0)
2472 +0.000467*ewgc("Cqu8R",2,2,1,1)+-0.010718*ewgc("CuBR",2,2)
2473 +0.004924*ewgc("Cud8R",2,2,0,0)+0.000555*ewgc("Cud8R",2,2,1,1)
2474 -0.082799*ewgc("CuGR",2,2)+0.038642*ewgc("CuuR",0,2,2,0)
2475 +0.001064*ewgc("CuuR",1,2,2,1)+-0.006289*ewgc("CuWR",2,2)
2476 +-0.033790*ewgc("CHD")+-0.019417*ewgc("CHl3R",0,0)+-0.019417*ewgc("CHl3R",1,1)
2477 +-0.072465*ewgc("CHWB")+0.019405*ewgc("CllR",0,1,1,0)+-0.000134*ewgc("Cqd1R",2,2,0,0)
2478 +-0.000033*ewgc("Cqd1R",2,2,1,1)+0.002483*ewgc("Cqq1R",0,0,2,2)+-0.000227*ewgc("Cqq1R",1,1,2,2)
2479 +0.008481*ewgc("Cqq3R",0,0,2,2)+0.000442*ewgc("Cqq3R",1,1,2,2)+0.001386*ewgc("Cqu1R",0,0,2,2)
2480 +0.000011*ewgc("Cqu1R",1,1,2,2)+0.000965*ewgc("Cqu1R",2,2,0,0)+0.000015*ewgc("Cqu1R",2,2,1,1)
2481 +-0.000494*ewgc("Cud1R",2,2,0,0)+-0.000044*ewgc("Cud1R",2,2,1,1)+0.002770*ewgc("CuuR",0,0,2,2)
2482 +0.000056*ewgc("CuuR",1,1,2,2))
2483 *(SM_sigma_tta_bin_50_80/sigma_tta_bin_50_80_madgraph);
2484
2485 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2486
2487 return total;
2488
2489
2490 }
2491 } else if(b_min == 80 && b_max == 120){
2492
2493
2494 double SM_sigma_tta_bin_80_120 = SM.getOptionalParameter("SM_sigma_tta_bin_80_120_CMS_semileptonic");
2495 double sigma_tta_bin_80_120_madgraph = 0.198070; //pb
2496
2497
2498 if(flag_Quadratic){
2499
2500 return 0.;
2501
2502 }
2503 else{
2504
2505 double total = SM_sigma_tta_bin_80_120 +(-0.013527*ewgc("CG")+0.001628*ewgc("Cqd8R",2,2,0,0)
2506 +0.000164*ewgc("Cqd8R",2,2,1,1)+0.023160*ewgc("Cqq1R",0,2,2,0)
2507 +0.000587*ewgc("Cqq1R",1,2,2,1)+0.047833*ewgc("Cqq3R",0,2,2,0)
2508 +0.003240*ewgc("Cqq3R",1,2,2,1)+0.012405*ewgc("Cqu8R",0,0,2,2)
2509 +0.000452*ewgc("Cqu8R",1,1,2,2)+0.010800*ewgc("Cqu8R",2,2,0,0)
2510 +0.000273*ewgc("Cqu8R",2,2,1,1)+-0.009514*ewgc("CuBR",2,2)
2511 +0.003221*ewgc("Cud8R",2,2,0,0)+0.000351*ewgc("Cud8R",2,2,1,1)
2512 -0.050365*ewgc("CuGR",2,2)+0.022607*ewgc("CuuR",0,2,2,0)
2513 +0.000576*ewgc("CuuR",1,2,2,1)+-0.005448*ewgc("CuWR",2,2)
2514 +-0.020954*ewgc("CHD")+-0.012043*ewgc("CHl3R",0,0)+-0.012043*ewgc("CHl3R",1,1)
2515 +-0.044910*ewgc("CHWB")+0.012023*ewgc("CllR",0,1,1,0)+-0.000046*ewgc("Cqd1R",2,2,0,0)
2516 +-0.000006*ewgc("Cqd1R",2,2,1,1)+0.001338*ewgc("Cqq1R",0,0,2,2)+-0.000141*ewgc("Cqq1R",1,1,2,2)
2517 +0.005279*ewgc("Cqq3R",0,0,2,2)+0.000310*ewgc("Cqq3R",1,1,2,2)+0.000893*ewgc("Cqu1R",0,0,2,2)
2518 +0.000039*ewgc("Cqu1R",1,1,2,2)+0.000660*ewgc("Cqu1R",2,2,0,0)+0.000008*ewgc("Cqu1R",2,2,1,1)
2519 +-0.000280*ewgc("Cud1R",2,2,0,0)+-0.000026*ewgc("Cud1R",2,2,1,1)+0.001709*ewgc("CuuR",0,0,2,2)
2520 +0.000051*ewgc("CuuR",1,1,2,2))
2521 *(SM_sigma_tta_bin_80_120/sigma_tta_bin_80_120_madgraph);
2522
2523 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2524
2525 return total;
2526
2527 }
2528 } else if(b_min == 120 && b_max == 160){
2529
2530
2531 double SM_sigma_tta_bin_120_160 = SM.getOptionalParameter("SM_sigma_tta_bin_120_160_CMS_semileptonic");
2532 double sigma_tta_bin_120_160_madgraph = 0.097725; //pb
2533
2534
2535 if(flag_Quadratic){
2536
2537 return 0.;
2538
2539 }
2540 else{
2541
2542 double total = SM_sigma_tta_bin_120_160 +(-0.007834*ewgc("CG")+0.000808*ewgc("Cqd8R",2,2,0,0)
2543 +0.000071*ewgc("Cqd8R",2,2,1,1)+0.011182*ewgc("Cqq1R",0,2,2,0)
2544 +0.000245*ewgc("Cqq1R",1,2,2,1)+0.024260*ewgc("Cqq3R",0,2,2,0)
2545 +0.001549*ewgc("Cqq3R",1,2,2,1)+0.006409*ewgc("Cqu8R",0,0,2,2)
2546 +0.000204*ewgc("Cqu8R",1,1,2,2)+0.005568*ewgc("Cqu8R",2,2,0,0)
2547 +0.000122*ewgc("Cqu8R",2,2,1,1)+-0.005632*ewgc("CuBR",2,2)
2548 +0.001708*ewgc("Cud8R",2,2,0,0)+0.000162*ewgc("Cud8R",2,2,1,1)
2549 -0.024540*ewgc("CuGR",2,2)+0.010926*ewgc("CuuR",0,2,2,0)
2550 +0.000235*ewgc("CuuR",1,2,2,1)+-0.003185*ewgc("CuWR",2,2)
2551 +-0.010353*ewgc("CHD")+-0.005952*ewgc("CHl3R",0,0)+-0.005952*ewgc("CHl3R",1,1)+-0.022160*ewgc("CHWB")
2552 +0.005918*ewgc("CllR",0,1,1,0)+-0.000019*ewgc("Cqd1R",2,2,0,0)+-0.000016*ewgc("Cqd1R",2,2,1,1)
2553 +0.000566*ewgc("Cqq1R",0,0,2,2)+-0.000079*ewgc("Cqq1R",1,1,2,2)+0.002746*ewgc("Cqq3R",0,0,2,2)
2554 +0.000133*ewgc("Cqq3R",1,1,2,2)+0.000456*ewgc("Cqu1R",0,0,2,2)+-0.000005*ewgc("Cqu1R",1,1,2,2)
2555 +0.000346*ewgc("Cqu1R",2,2,0,0)+-0.000009*ewgc("Cqu1R",2,2,1,1)+-0.000154*ewgc("Cud1R",2,2,0,0)
2556 +-0.000030*ewgc("Cud1R",2,2,1,1)+0.000819*ewgc("CuuR",0,0,2,2)+0.000004*ewgc("CuuR",1,1,2,2))
2557 *(SM_sigma_tta_bin_120_160/sigma_tta_bin_120_160_madgraph);
2558
2559 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2560
2561 return total;
2562
2563 }
2564 } else if(b_min == 160 && b_max == 200){
2565
2566
2567 double SM_sigma_tta_bin_160_200 = SM.getOptionalParameter("SM_sigma_tta_bin_160_200_CMS_semileptonic");
2568 double sigma_tta_bin_160_200_madgraph = 0.053162; //pb
2569
2570
2571 if(flag_Quadratic){
2572
2573 return 0.;
2574
2575 }
2576 else{
2577
2578 double total = SM_sigma_tta_bin_160_200 +(-0.005116*ewgc("CG")+0.000482*ewgc("Cqd8R",2,2,0,0)
2579 +0.000047*ewgc("Cqd8R",2,2,1,1)+0.006230*ewgc("Cqq1R",0,2,2,0)
2580 +0.000137*ewgc("Cqq1R",1,2,2,1)+0.013887*ewgc("Cqq3R",0,2,2,0)
2581 +0.000853*ewgc("Cqq3R",1,2,2,1)+0.003756*ewgc("Cqu8R",0,0,2,2)
2582 +0.000123*ewgc("Cqu8R",1,1,2,2)+0.003272*ewgc("Cqu8R",2,2,0,0)
2583 +0.000076*ewgc("Cqu8R",2,2,1,1)+-0.003224*ewgc("CuBR",2,2)
2584 +0.001020*ewgc("Cud8R",2,2,0,0)+0.000101*ewgc("Cud8R",2,2,1,1)
2585 -0.013391*ewgc("CuGR",2,2)+0.006063*ewgc("CuuR",0,2,2,0)
2586 +0.000133*ewgc("CuuR",1,2,2,1)+-0.001792*ewgc("CuWR",2,2)
2587 +-0.005617*ewgc("CHD")+-0.003223*ewgc("CHl3R",0,0)+-0.003223*ewgc("CHl3R",1,1)
2588 +-0.012030*ewgc("CHWB")+0.003235*ewgc("CllR",0,1,1,0)+-0.000003*ewgc("Cqd1R",2,2,0,0)
2589 +0.000004*ewgc("Cqd1R",2,2,1,1)+0.000283*ewgc("Cqq1R",0,0,2,2)+-0.000035*ewgc("Cqq1R",1,1,2,2)
2590 +0.001579*ewgc("Cqq3R",0,0,2,2)+0.000085*ewgc("Cqq3R",1,1,2,2)+0.000238*ewgc("Cqu1R",0,0,2,2)
2591 +0.000009*ewgc("Cqu1R",1,1,2,2)+0.000208*ewgc("Cqu1R",2,2,0,0)+0.000010*ewgc("Cqu1R",2,2,1,1)
2592 +-0.000083*ewgc("Cud1R",2,2,0,0)+-0.000005*ewgc("Cud1R",2,2,1,1)+0.000471*ewgc("CuuR",0,0,2,2)
2593 +0.000014*ewgc("CuuR",1,1,2,2))
2594 *(SM_sigma_tta_bin_160_200/sigma_tta_bin_160_200_madgraph);
2595
2596 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2597
2598 return total;
2599
2600
2601 }
2602 } else if(b_min == 200 && b_max == 260){
2603
2604
2605 double SM_sigma_tta_bin_200_260 = SM.getOptionalParameter("SM_sigma_tta_bin_200_260_CMS_semileptonic");
2606 double sigma_tta_bin_200_260_madgraph = 0.040676; //pb
2607
2608
2609 if(flag_Quadratic){
2610
2611 return 0.;
2612
2613 }
2614 else{
2615
2616 double total = SM_sigma_tta_bin_200_260 +(-0.004958*ewgc("CG")+0.000418*ewgc("Cqd8R",2,2,0,0)
2617 +0.000036*ewgc("Cqd8R",2,2,1,1)+0.005096*ewgc("Cqq1R",0,2,2,0)
2618 +0.000101*ewgc("Cqq1R",1,2,2,1)+0.011849*ewgc("Cqq3R",0,2,2,0)
2619 +0.000680*ewgc("Cqq3R",1,2,2,1)+0.003215*ewgc("Cqu8R",0,0,2,2)
2620 +0.000093*ewgc("Cqu8R",1,1,2,2)+0.002804*ewgc("Cqu8R",2,2,0,0)
2621 +0.000057*ewgc("Cqu8R",2,2,1,1)+-0.002405*ewgc("CuBR",2,2)
2622 +0.000885*ewgc("Cud8R",2,2,0,0)+0.000078*ewgc("Cud8R",2,2,1,1)
2623 -0.010460*ewgc("CuGR",2,2)+0.004970*ewgc("CuuR",0,2,2,0)
2624 +0.000098*ewgc("CuuR",1,2,2,1)+-0.001356*ewgc("CuWR",2,2)
2625 +-0.004304*ewgc("CHD")+-0.002471*ewgc("CHl3R",0,0)+-0.002471*ewgc("CHl3R",1,1)
2626 +-0.009223*ewgc("CHWB")+0.002470*ewgc("CllR",0,1,1,0)+-0.000011*ewgc("Cqd1R",2,2,0,0)
2627 +0.000001*ewgc("Cqd1R",2,2,1,1)+0.000251*ewgc("Cqq1R",0,0,2,2)+-0.000039*ewgc("Cqq1R",1,1,2,2)
2628 +0.001359*ewgc("Cqq3R",0,0,2,2)+0.000063*ewgc("Cqq3R",1,1,2,2)+0.000204*ewgc("Cqu1R",0,0,2,2)
2629 +0.000004*ewgc("Cqu1R",1,1,2,2)+0.000151*ewgc("Cqu1R",2,2,0,0)+0.000004*ewgc("Cqu1R",2,2,1,1)
2630 +-0.000074*ewgc("Cud1R",2,2,0,0)+-0.000005*ewgc("Cud1R",2,2,1,1)+0.000416*ewgc("CuuR",0,0,2,2)
2631 +0.000007*ewgc("CuuR",1,1,2,2))
2632 *(SM_sigma_tta_bin_200_260/sigma_tta_bin_200_260_madgraph);
2633
2634
2635 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2636
2637 return total;
2638
2639 }
2640 } else if(b_min == 260 && b_max == 320){
2641
2642
2643 double SM_sigma_tta_bin_260_320 = SM.getOptionalParameter("SM_sigma_tta_bin_260_320_CMS_semileptonic");
2644 double sigma_tta_bin_260_320_madgraph = 0.018994; //pb
2645
2646
2647 if(flag_Quadratic){
2648
2649 return 0.;
2650
2651 }
2652 else{
2653
2654 double total = SM_sigma_tta_bin_260_320 +(-0.003098*ewgc("CG")+0.000227*ewgc("Cqd8R",2,2,0,0)
2655 +0.000013*ewgc("Cqd8R",2,2,1,1)+0.002712*ewgc("Cqq1R",0,2,2,0)
2656 +0.000044*ewgc("Cqq1R",1,2,2,1)+0.006445*ewgc("Cqq3R",0,2,2,0)
2657 +0.000337*ewgc("Cqq3R",1,2,2,1)+0.001768*ewgc("Cqu8R",0,0,2,2)
2658 +0.000042*ewgc("Cqu8R",1,1,2,2)+0.001534*ewgc("Cqu8R",2,2,0,0)
2659 +0.000025*ewgc("Cqu8R",2,2,1,1)+-0.001038*ewgc("CuBR",2,2)
2660 +0.000492*ewgc("Cud8R",2,2,0,0)+0.000035*ewgc("Cud8R",2,2,1,1)
2661 -0.005064*ewgc("CuGR",2,2)+0.002652*ewgc("CuuR",0,2,2,0)
2662 +0.000042*ewgc("CuuR",1,2,2,1)+-0.000592*ewgc("CuWR",2,2)
2663 +-0.002010*ewgc("CHD")+-0.001156*ewgc("CHl3R",0,0)+-0.001156*ewgc("CHl3R",1,1)
2664 +-0.004307*ewgc("CHWB")+0.001150*ewgc("CllR",0,1,1,0)+-0.000011*ewgc("Cqd1R",2,2,0,0)
2665 +-0.000005*ewgc("Cqd1R",2,2,1,1)+0.000094*ewgc("Cqq1R",0,0,2,2)+-0.000021*ewgc("Cqq1R",1,1,2,2)
2666 +0.000730*ewgc("Cqq3R",0,0,2,2)+0.000028*ewgc("Cqq3R",1,1,2,2)+0.000112*ewgc("Cqu1R",0,0,2,2)
2667 +-0.000003*ewgc("Cqu1R",1,1,2,2)+0.000081*ewgc("Cqu1R",2,2,0,0)+-0.000003*ewgc("Cqu1R",2,2,1,1)
2668 +-0.000049*ewgc("Cud1R",2,2,0,0)+-0.000008*ewgc("Cud1R",2,2,1,1)+0.000214*ewgc("CuuR",0,0,2,2)
2669 +-0.000002*ewgc("CuuR",1,1,2,2))
2670 *(SM_sigma_tta_bin_260_320/sigma_tta_bin_260_320_madgraph);
2671
2672 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2673
2674 return total;
2675
2676
2677 }
2678 } else if(b_min == 320 && b_max == 380){
2679
2680
2681 double SM_sigma_tta_bin_320_380 = SM.getOptionalParameter("SM_sigma_tta_bin_320_380_CMS_semileptonic");
2682 double sigma_tta_bin_320_380_madgraph = 0.009391; //pb
2683
2684
2685 if(flag_Quadratic){
2686
2687 return 0.;
2688
2689 }
2690 else{
2691
2692 double total = SM_sigma_tta_bin_320_380 +(-0.001921*ewgc("CG")+0.000136*ewgc("Cqd8R",2,2,0,0)
2693 +0.000007*ewgc("Cqd8R",2,2,1,1)+0.001575*ewgc("Cqq1R",0,2,2,0)
2694 +0.000023*ewgc("Cqq1R",1,2,2,1)+0.003819*ewgc("Cqq3R",0,2,2,0)
2695 +0.000186*ewgc("Cqq3R",1,2,2,1)+0.001048*ewgc("Cqu8R",0,0,2,2)
2696 +0.000022*ewgc("Cqu8R",1,1,2,2)+0.000910*ewgc("Cqu8R",2,2,0,0)
2697 +0.000013*ewgc("Cqu8R",2,2,1,1)+-0.000451*ewgc("CuBR",2,2)
2698 +0.000293*ewgc("Cud8R",2,2,0,0)+0.000019*ewgc("Cud8R",2,2,1,1)
2699 -0.002605*ewgc("CuGR",2,2)+0.001534*ewgc("CuuR",0,2,2,0)
2700 +0.000023*ewgc("CuuR",1,2,2,1)+-0.000262*ewgc("CuWR",2,2)
2701 +-0.000996*ewgc("CHD")+-0.000574*ewgc("CHl3R",0,0)+-0.000574*ewgc("CHl3R",1,1)
2702 +-0.002131*ewgc("CHWB")+0.000567*ewgc("CllR",0,1,1,0)+-0.000004*ewgc("Cqd1R",2,2,0,0)
2703 +-0.000003*ewgc("Cqd1R",2,2,1,1)+0.000057*ewgc("Cqq1R",0,0,2,2)+-0.000012*ewgc("Cqq1R",1,1,2,2)
2704 +0.000442*ewgc("Cqq3R",0,0,2,2)+0.000016*ewgc("Cqq3R",1,1,2,2)+0.000069*ewgc("Cqu1R",0,0,2,2)
2705 +-0.000001*ewgc("Cqu1R",1,1,2,2)+0.000046*ewgc("Cqu1R",2,2,0,0)+-0.000001*ewgc("Cqu1R",2,2,1,1)
2706 +-0.000029*ewgc("Cud1R",2,2,0,0)+-0.000005*ewgc("Cud1R",2,2,1,1)+0.000125*ewgc("CuuR",0,0,2,2)
2707 +-0.0000005*ewgc("CuuR",1,1,2,2))
2708 *(SM_sigma_tta_bin_320_380/sigma_tta_bin_320_380_madgraph);
2709
2710 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2711
2712 return total;
2713
2714 }
2715 } else{
2716 throw std::runtime_error("\nERROR: Please specify a correct binning range for sigma_tta_diff_LO_CMS_semileptonic.\n");
2717 }
2718
2719
2720
2721}
QCD::getOptionalParameter
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
ThObservable::getBinMin
double getBinMin()
A get method to get the minimum value of the bin.
Definition: ThObservable.h:82
ThObservable::getBinMax
double getBinMax()
A get method to get the maximum value of the bin.
Definition: ThObservable.h:91
sigma_tta_diff_LO_CMS_semileptonic_210701508::b_min
double b_min
Definition: TopQuarkObservables.h:709
sigma_tta_diff_LO_CMS_semileptonic_210701508::ewgc
double ewgc(const std::string name) const
Definition: TopQuarkObservables.h:692
sigma_tta_diff_LO_CMS_semileptonic_210701508::b_max
double b_max
Definition: TopQuarkObservables.h:710

◆ ewgc() [1/3]

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

Definition at line 692 of file TopQuarkObservables.h.

693 {
694 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
695 }
NPSMEFTd6General::getSMEFTCoeffEW
double getSMEFTCoeffEW(const std::string name) const
Definition: NPSMEFTd6General.h:6790
TopQuarkObservables::GetmyNPSMEFTd6General
const NPSMEFTd6General & GetmyNPSMEFTd6General() const
Definition: TopQuarkObservables.h:27

◆ ewgc() [2/3]

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

Definition at line 697 of file TopQuarkObservables.h.

698 {
699 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
700 }

◆ ewgc() [3/3]

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

Definition at line 702 of file TopQuarkObservables.h.

703 {
704 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
705 }

Member Data Documentation

◆ b_max

double sigma_tta_diff_LO_CMS_semileptonic_210701508::b_max = 0.
private

Definition at line 710 of file TopQuarkObservables.h.

◆ b_min

double sigma_tta_diff_LO_CMS_semileptonic_210701508::b_min = 0.
private

Definition at line 709 of file TopQuarkObservables.h.

◆ mytopobs

const TopQuarkObservables sigma_tta_diff_LO_CMS_semileptonic_210701508::mytopobs
private

Definition at line 713 of file TopQuarkObservables.h.

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