This page documents the location of MC sample for the ILC angular analysis.

JHU

Spin 0 signals

If the absolute link of the lhe file is not there, take the location by the root file.

sqrts g1 g2 g4 fa2(phase) fa3(phase) location root file lhe file nEvents created by
250 1 0 0 0 0 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130618/ unweighted_unpol_g1_1M_false.root unweighted_unpol_g1_1M.lhe 1M Yaofu
250 0 1 0 1 0 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130618/ unweighted_unpol_g2_1M_false.root unweighted_unpol_g2_1M.lhe 1M Yaofu
250 0 0 1 0 1 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130618/ unweighted_unpol_g4_1M_false.root unweighted_unpol_g4_1M.lhe 1M Yaofu
250 1 0.54071 0 0.908844 0 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130618/ unweighted_unpol_g1_p_g2_1M_false.root unweighted_unpol_g1_p_g2_1M.lhe 1M Yaofu
250 1 i(0.54071) 0 0.908844 (pi/2) 0 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130618/ unweighted_unpol_g1_p_ig2_1M_false.root unweighted_unpol_g1_p_ig2_1M.lhe 1M Yaofu
250 1 0 0.83265 0 0.848419 export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130626/ unweighted_unpol_fa3_250GeV_5M_false.root /export/d1/scratch/yaofu/Events_20130626/unweighted_unpol_fa3_250GeV_5M.lhe 5M Yaofu
250 1 0 i0.83265 0 0.848419(pi/2) /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130618/ unweighted_unpol_g1_p_ig4_1M_false.root unweighted_unpol_g1_p_ig4_1M.lhe 1M Yaofu
250 1 0.54071 0.83265 0.601801 0.33784 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130618/ unweighted_unpol_g1_p_g2_p_g4_1M_false.root unweighted_unpol_g1_p_g2_p_g4_1M.lhe 1M Yaofu
250 1 i(0.54071) 0.83265 0.601801(pi/2) 0.33784 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130618/ unweighted_unpol_g1_p_ig2_p_g4_1M_false.root unweighted_unpol_g1_p_ig2_p_g4_1M.lhe 1M Yaofu
250 1 0.54071 i0.83265 0.601801 0.33784(pi/2) /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130618/ unweighted_unpol_g1_p_g2_ig4_1M_false.root unweighted_unpol_g1_p_g2_ig4_1M.lhe 1M Yaofu
250 1 0 0.117316 0 0.1 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130626/ unweighted_unpol_f_3_250GeV_5M_false.root /export/d1/scratch/yaofu/Events_20130626/unweighted_unpol_f_3_250GeV_5M.lhe 5M Yaofu
250 1 0 0.0807427 0 0.05 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130730/ unweighted_unpol_f_3_p05_5M_false.root unweighted_unpol_f_3_p05_5M.lhe 5M yaofu
250 1 0 0.0618949 0 0.03 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130730/ unweighted_unpol_f_3_p03_10M_false.root unweighted_unpol_f_3_p03_10M.lhe 10M yaofu
250 1 0 0.35194 0 0.5 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130829/ ee_ZH_f3_p5_phi_0_250GeV_1M_false.root ee_ZH_f3_p5_phi_0_250GeV_1M.lhe 1M yaofu
250 1 0 i0.35194 0 0.5 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130829/ ee_ZH_f3_p5_phi_pi2_250GeV_1M_false.root ee_ZH_f3_p5_phi_pi2_250GeV_1M.lhe 1M yaofu
|*350* | 1 | 0 | 0 | 1 | 0 | /export/d1/scratch/snowmass2013/ilcsample
350 0 0 1 0 1 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130701/ unweighted_unpol_g4_350GeV_2M_false.root unweighted_unpol_g4_350GeV_2M.lhe 2M yaofu
350 1 0 4.6863E-02 0 0.1 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130701/ unweighted_unpol_model8_2M_false.root unweighted_unpol_model8_2M.lhe 2M yaofu
500 1 0 2.62636E-02 0 0.1 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130626/ unweighted_unpol_f_3_500GeV_5M_false.root /export/d1/scratch/yaofu/Events_20130626/unweighted_unpol_f_3_500GeV_5M.lhe 5M Yaofu
500 1 0 0 0 0 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130701/ unweighted_unpol_g1_500GeV_2M_false.root unweighted_unpol_g1_500GeV_2M.lhe 2M yaofu
500 0 0 1 0 1 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130701/ unweighted_unpol_g4_500GeV_2M_false.root unweighted_unpol_g4_500GeV_2M.lhe 2M yaofu
1000 0 0 1 0 1 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130701/ unweighted_unpol_g4_1TeV_1M_false.root unweighted_unpol_g4_1TeV_1M.lhe 1M yaofu
1000 1 0 0 0 0 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130701/ unweighted_unpol_g1_1TeV_1M_false.root unweighted_unpol_g1_1TeV_1M.lhe 1M yaofu
1000 1 0 0.011296 0 0.1 /export/d1/scratch/snowmass2013/ilcsamples/signal/spin0/Events_20130701/ unweighted_unpol_model9_1M_false.root unweighted_unpol_model9_1M.lhe 1M yaofu

f2/f3 calcualtions

  • Calculate the fractions in terms of the couplings for ILC/LHC

void ilc_fa3()
{
// ILC numbers at 1000 GeV at mH= 125 GeV (narrow Z width approximation)
   Double_t sigma1_v = 8.95721; 
   Double_t sigma2_v = 8208.91;  
   Double_t sigma4_v = 7800.2;
// ILC numbers at 500 GeV at mH= 125 GeV (narrow Z width approximation)
   Double_t sigma1_h = 2.57246; 
   Double_t sigma2_h = 516.556;  
   Double_t sigma4_h = 414.378;
// ILC numbers at 350 GeV at mH= 125 GeV (narrow Z width approximation)
   Double_t sigma1_m = 1.48872; 
   Double_t sigma2_m = 125.387;  
   Double_t sigma4_m = 75.3199;
// ILC numbers at 250 GeV at mH= 125 GeV (narrow Z width approximation)
   Double_t sigma1_e = 0.981396;  // was 0.946962 at 126 GeV
   Double_t sigma2_e = 33.4674;  // was 32.1981 at 126 GeV
   Double_t sigma4_e = 7.9229; // was 7.45502 at 126 GeV
// LHC numbers at mH= 125 GeV (JHUGen)
   Double_t sigma1_p = 1.860351; // was 2.03971 at 126 GeV
   Double_t sigma2_p = 0.672859; // was 0.77517 at 126 GeV
   Double_t sigma4_p = 0.284353; // was 0.32689 at 126 GeV
//
   Double_t g1 = 1.;
   Double_t g2 = 0.;
   Double_t g3 = 0.;
   Double_t g4 = 0.1;
//
   Double_t fg2_v = sigma2_v*g2*g2 / (  sigma1_v*g1*g1 + sigma2_v*g2*g2 + sigma4_v*g4*g4 );
   Double_t fg4_v = sigma4_v*g4*g4 / (  sigma1_v*g1*g1 + sigma2_v*g2*g2 + sigma4_v*g4*g4 );
   Double_t fg2_h = sigma2_h*g2*g2 / (  sigma1_h*g1*g1 + sigma2_h*g2*g2 + sigma4_h*g4*g4 );
   Double_t fg4_h = sigma4_h*g4*g4 / (  sigma1_h*g1*g1 + sigma2_h*g2*g2 + sigma4_h*g4*g4 );
   Double_t fg2_m = sigma2_m*g2*g2 / (  sigma1_m*g1*g1 + sigma2_m*g2*g2 + sigma4_m*g4*g4 );
   Double_t fg4_m = sigma4_m*g4*g4 / (  sigma1_m*g1*g1 + sigma2_m*g2*g2 + sigma4_m*g4*g4 );
   Double_t fg2_e = sigma2_e*g2*g2 / (  sigma1_e*g1*g1 + sigma2_e*g2*g2 + sigma4_e*g4*g4 );
   Double_t fg4_e = sigma4_e*g4*g4 / (  sigma1_e*g1*g1 + sigma2_e*g2*g2 + sigma4_e*g4*g4 );
   Double_t fg2_p = sigma2_p*g2*g2 / (  sigma1_p*g1*g1 + sigma2_p*g2*g2 + sigma4_p*g4*g4 );
   Double_t fg4_p = sigma4_p*g4*g4 / (  sigma1_p*g1*g1 + sigma2_p*g2*g2 + sigma4_p*g4*g4 );
//
   cout << " " << endl;   
   cout << " g1= "<< g1 << endl;   
   cout << " g2= "<< g2 << endl;   
   cout << " g3= "<< g3 << endl;   
   cout << " g4= "<< g4 << endl;   
   cout << " " << endl;   
   cout << " LHC (mH=125 GeV):   fa2 = " << fg2_p << "   fa3 = "<< fg4_p << endl;   
   cout << " ILC (250 GeV):      f_2 = " << fg2_e << "   f_3 = "<< fg4_e << endl;   
   cout << " ILC (350 GeV):      f_2 = " << fg2_m << "   f_3 = "<< fg4_m << endl;   
   cout << " ILC (500 GeV):      f_2 = " << fg2_h << "   f_3 = "<< fg4_h << endl;   
   cout << " ILC (1000 GeV):     f_2 = " << fg2_v << "   f_3 = "<< fg4_v << endl;   
   cout << " " << endl;  
}

  • Calculate the fractions in terms of fa3 converting from f_3 on ILC

void converter_fa3()
{
//
   Double_t f_3 = 0.1;
   Double_t energy = 350.;
   cout << " fa3 = "<< convert_to_fa3(f_3, energy) << endl;  
   cout << " " << endl;  
//   
} 
//
Double_t convert_to_fa3(Double_t f_3, Double_t energy)
{
// LHC numbers at mH= 125 GeV (JHUGen)
   Double_t sigma1_p = 1.860351; // was 2.03971 at 126 GeV
   Double_t sigma2_p = 0.672859; // was 0.77517 at 126 GeV
   Double_t sigma4_p = 0.284353; // was 0.32689 at 126 GeV
// ILC numbers at 250 GeV at mH= 125 GeV (narrow Z width approximation)
   Double_t sigma1_e = 0.981396;  // was 0.946962 at 126 GeV
   Double_t sigma2_e = 33.4674;  // was 32.1981 at 126 GeV
   Double_t sigma4_e = 7.9229; // was 7.45502 at 126 GeV
// ILC numbers at 350 GeV at mH= 125 GeV (narrow Z width approximation)
   if(fabs(energy-350.)<1.){
     sigma1_e = 1.48872; 
     sigma2_e = 125.387;  
     sigma4_e = 75.3199;
   }
// ILC numbers at 500 GeV at mH= 125 GeV (narrow Z width approximation)
   if(fabs(energy-500.)<1.){
     sigma1_e = 2.57246; 
     sigma2_e = 516.556;  
     sigma4_e = 414.378;
   }
//------------
   Double_t fg2_e = 0.;
   Double_t fg4_e = f_3;
//------------
   Double_t g1 = 1.;
   Double_t g2 = sqrt(fg2_e/(1-fg2_e-fg4_e))*sqrt(sigma1_e/sigma2_e);
   Double_t g4 = sqrt(fg4_e/(1-fg2_e-fg4_e))*sqrt(sigma1_e/sigma4_e);
//------------
   Double_t fg2_p = sigma2_p*g2*g2 / (  sigma1_p*g1*g1 + sigma2_p*g2*g2 + sigma4_p*g4*g4 );
   Double_t fg4_p = sigma4_p*g4*g4 / (  sigma1_p*g1*g1 + sigma2_p*g2*g2 + sigma4_p*g4*g4 );
//------------
   return fg4_p;
}

  • Andrei's model generation notes

Using LHC cross-sections:

   Double_t sigma1_p = 1.860351; // was 2.03971 at 126 GeV
   Double_t sigma2_p = 0.672859; // was 0.77517 at 126 GeV
  
Edit | Attach | Watch | Print version | History: r30 < r29 < r28 < r27 < r26 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r30 - 2013-08-31 - YanyanGao
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback