Ticket #756: test_pol_sm.sin

printf "Starting polarization test: t tbar events"

openmp_num_threads = 2

model = SM

process eett  = e1, E1 => t, T

#-----------------------------------------------------------------
# Monoenergetic beam

beams = e1, E1

sqrts = 500 GeV

seed = 1

integrate (eett) { iterations = 2:2000 }

# Loop over e+ / e- helicity
scan int hel_em = (-1, 1) {
  scan int hel_ep = (-1, 1) {

    beams_pol_density = @(hel_em), @(hel_ep)

    integrate (eett) { iterations = 2:2000 }
  }
}




#-----------------------------------------------------------------
# ILC spectra from CIRCE1 + ISR

beams = e1, E1 => circe1  => isr

sqrts = 500 GeV
circe1_ver = 10
$circe1_acc = "ILC"

seed = 1

integrate (eett) { iterations = 2:2000 }

# Loop over e+ / e- helicity
scan int hel_em = (-1, 1) {
  scan int hel_ep = (-1, 1) {

    beams_pol_density = @(hel_em), @(hel_ep)

    integrate (eett) { iterations = 2:2000 }
  }
}

#-----------------------------------------------------------------
# Test spectra file

beams = e1, E1 => beam_events => isr
$beam_events_file = "uniform_spread_2.5%.dat"

sqrts = 500 GeV

seed = 1

integrate (eett) { iterations = 2:2000 }

# Loop over e+ / e- helicity
scan int hel_em = (-1, 1) {
  scan int hel_ep = (-1, 1) {

    beams_pol_density = @(hel_em), @(hel_ep)

    integrate (eett) { iterations = 2:2000 }
  }
}