Ticket #161: Run.LHCB_T_Asym.tt_Zprime.dat

#export RIVET_ANALYSIS_PATH=$PWD

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

  ANALYSIS = Rivet
  ANALYSIS_OUTPUT = Asym

  BATCH_MODE 0; OUTPUT 2;

  ## full description of all parameters with
  ## Sherpa SHOW_MODEL_SYNTAX=1
  MODEL      = SM+U1_B
  ## mass and coupling of leptophobic Z'
  M_Z' = 160;
  g'_1 = 0.5;
  MIX  = 0;
  MASSIVE[5] 1  
  YUYKAWA[5] 0
  ACTIVE[23] 0
}(run);
(analysis){
  BEGIN_RIVET {
    -a LHCB_TOP_ASYMMETRY
  } END_RIVET
}(analysis)
(processes){
  Process 93 93 -> 6[a] -6[b]
  Decay  6[a] -> -13  14  5
  Decay -6[b] ->  93  93 -5
  Order_EW  2;
  Order_QCD 0;
  End process;

  Process 93 93 -> 6[a] -6[b]
  Decay  6[a] ->  93  93  5
  Decay -6[b] ->  13 -14 -5
  Order_EW  2;
  Order_QCD 0;
  End process;

  Process 2 2 -> 6[a] 6[b]
  Decay  6[a] -> -13  14  5
  Decay  6[b] ->  93  93  5
  Order_EW  2;
  Order_QCD 0;
  End process;

  Process 93 93 -> -6[a] -6[b]
  Decay -6[a] ->  93  93 -5
  Decay -6[b] ->  13 -14 -5
  Order_EW  2;
  Order_QCD 0;
  End process;
}(processes);
(selector){
  Mass  13 -14 10 E_CMS;
  Mass -13  14 10 E_CMS;
  PT     5     10 E_CMS;
  PT    -5     10 E_CMS;
}(selector)
(me){
  ME_SIGNAL_GENERATOR Amegic;
  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);