Ticket #166: Run.LHCB_T_Asym.tt.HL.dat

#export RIVET_ANALYSIS_PATH=$PWD

(run){
  EVENTS 1000000;
  RESULT_DIRECTORY = Results_TT_SM
  GENERATE_RESULT_DIRECTORY = 1

  #ANALYSIS = Rivet
  #ANALYSIS_OUTPUT = Asym

  BATCH_MODE 0; OUTPUT 2;
}(run);
(analysis){
  BEGIN_RIVET {
    -a TTBAR_ASYM
  } END_RIVET
}(analysis)
(processes){
  Process 93 93 ->   6[a] -6[b] 93{1}
  Decay    6[a] ->  24[c]     5
  Decay   24[c] ->    -13    14
  Decay   -6[b] -> -24[d]    -5
  Decay  -24[d] ->     93    93
  CKKW sqr(20/E_CMS)
  Order_EW 4;
  End process;

  #Process 93 93 -> 6[a]  -6[b] 93{1}
  #Decay    6[a] ->  93  93   5
  #Decay   -6[b] ->  13 -14  -5
  #CKKW sqr(20/E_CMS)
  #Order_EW 0;
  #End process;

  #Process 93 93 -> 6[a]  93 93{1}
  #Decay    6[a] -> -13  14   5
  #CKKW sqr(20/E_CMS)
  #Order_EW 2;
  #End process;

  #Process 93 93 -> 93 -6[a] 93{1}
  #Decay   -6[a] ->  13 -14   5
  #CKKW sqr(20/E_CMS)
  #Order_EW 2;
  #End process;
}(processes);
(selector){     
  Mass  93   93   10 E_CMS;     
  Mass  13  -14   10 E_CMS;     
  Mass -13   14   10 E_CMS;     
  PT     6         2 E_CMS;
  PT    -6         2 E_CMS;
  PT     5         2 E_CMS;
  PT    -5         2 E_CMS;
}(selector)
(me){
  ME_SIGNAL_GENERATOR Comix
  #SCALES VAR{Abs2(p[0]+p[1])};
}(me);
(beam){
  BEAM_1  2212; BEAM_ENERGY_1 3500;
  BEAM_2  2212; BEAM_ENERGY_2 3500;
}(beam);
(integration){
  INTEGRATOR 7; ERROR 0.1;
  FINISH_OPTIMIZATION On;
}(integration);