Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F8723562
eeww_polarized.sin
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
10 KB
Subscribers
None
eeww_polarized.sin
View Options
# $Id: eeww_polarized.sin 7444 2016-02-17 15:37:20Z jr_reuter $
# Use WHIZARD's initial and final state polarization facilities to
# examine the angular distributions for all 36 helicity combinations
# in e+, e- -> W+, W-
########################################################################
#
# Copyright (C) 1999-2016 by
# Wolfgang Kilian <kilian@physik.uni-siegen.de>
# Thorsten Ohl <ohl@physik.uni-wuerzburg.de>
# Juergen Reuter <juergen.reuter@desy.de>
# with contributions from
# Christian Speckner <cnspeckn@googlemail.com>
#
# WHIZARD is free software; you can redistribute it and/or modify it
# under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2, or (at your option)
# any later version.
#
# WHIZARD is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
#
########################################################################
# Define and compile the process
model = SM
process eeww = "e+", "e-" => "W+", "W-"
compile
# Define the beams and event count (for simulation)
sqrts = 200 GeV
beams = "e+", "e-"
n_events = 100000
# Tell WHIZARD to retain the W polarization in the generated events
polarized "W+", "W-"
# Define title etc. as global variables, that will be used by PLOT
$description =
"A WHIZARD 2.2 Example. $W$ pair production in polarized ILC beams at 200 GeV. Analysis for polarized $W$s."
# Define histograms for all different final state helicies.
histogram cta_mm_mm (-1, 1, 0.1) { $title = "Polarization: -1 -1 -1 -1" }
histogram cta_pm_mm (-1, 1, 0.1) { $title = "Polarization: +1 -1 -1 -1" }
histogram cta_mp_mm (-1, 1, 0.1) { $title = "Polarization: -1 +1 -1 -1" }
histogram cta_pp_mm (-1, 1, 0.1) { $title = "Polarization: +1 +1 -1 -1" }
histogram cta_mm_m0 (-1, 1, 0.1) { $title = "Polarization: -1 -1 -1 0" }
histogram cta_pm_m0 (-1, 1, 0.1) { $title = "Polarization: +1 -1 -1 0" }
histogram cta_mp_m0 (-1, 1, 0.1) { $title = "Polarization: -1 +1 -1 0" }
histogram cta_pp_m0 (-1, 1, 0.1) { $title = "Polarization: +1 +1 -1 0" }
histogram cta_mm_mp (-1, 1, 0.1) { $title = "Polarization: -1 -1 -1 +1" }
histogram cta_pm_mp (-1, 1, 0.1) { $title = "Polarization: +1 -1 -1 +1" }
histogram cta_mp_mp (-1, 1, 0.1) { $title = "Polarization: -1 +1 -1 +1" }
histogram cta_pp_mp (-1, 1, 0.1) { $title = "Polarization: +1 +1 -1 +1" }
histogram cta_mm_0p (-1, 1, 0.1) { $title = "Polarization: -1 -1 0 +1" }
histogram cta_pm_0p (-1, 1, 0.1) { $title = "Polarization: +1 -1 0 +1" }
histogram cta_mp_0p (-1, 1, 0.1) { $title = "Polarization: -1 +1 0 +1" }
histogram cta_pp_0p (-1, 1, 0.1) { $title = "Polarization: +1 +1 0 +1" }
histogram cta_mm_00 (-1, 1, 0.1) { $title = "Polarization: -1 -1 0 0" }
histogram cta_pm_00 (-1, 1, 0.1) { $title = "Polarization: +1 -1 0 0" }
histogram cta_mp_00 (-1, 1, 0.1) { $title = "Polarization: -1 +1 0 0" }
histogram cta_pp_00 (-1, 1, 0.1) { $title = "Polarization: +1 +1 0 0" }
histogram cta_mm_0m (-1, 1, 0.1) { $title = "Polarization: -1 -1 0 -1" }
histogram cta_pm_0m (-1, 1, 0.1) { $title = "Polarization: +1 -1 0 -1" }
histogram cta_mp_0m (-1, 1, 0.1) { $title = "Polarization: -1 +1 0 -1" }
histogram cta_pp_0m (-1, 1, 0.1) { $title = "Polarization: +1 +1 0 -1" }
histogram cta_mm_pm (-1, 1, 0.1) { $title = "Polarization: -1 -1 +1 -1" }
histogram cta_pm_pm (-1, 1, 0.1) { $title = "Polarization: +1 -1 +1 -1" }
histogram cta_mp_pm (-1, 1, 0.1) { $title = "Polarization: -1 +1 +1 -1" }
histogram cta_pp_pm (-1, 1, 0.1) { $title = "Polarization: +1 +1 +1 -1" }
histogram cta_mm_p0 (-1, 1, 0.1) { $title = "Polarization: -1 -1 +1 0" }
histogram cta_pm_p0 (-1, 1, 0.1) { $title = "Polarization: +1 -1 +1 0" }
histogram cta_mp_p0 (-1, 1, 0.1) { $title = "Polarization: -1 +1 +1 0" }
histogram cta_pp_p0 (-1, 1, 0.1) { $title = "Polarization: +1 +1 +1 0" }
histogram cta_mm_pp (-1, 1, 0.1) { $title = "Polarization: -1 -1 +1 +1" }
histogram cta_pm_pp (-1, 1, 0.1) { $title = "Polarization: +1 -1 +1 +1" }
histogram cta_mp_pp (-1, 1, 0.1) { $title = "Polarization: -1 +1 +1 +1" }
histogram cta_pp_pp (-1, 1, 0.1) { $title = "Polarization: +1 +1 +1 +1" }
# Loop over e+ / e- helicity
scan int hel_ep = (-1, 1) {
scan int hel_em = (-1, 1) {
# Apply the helicity setup
beams_pol_density = @(hel_ep), @(hel_em)
# (Re)calculate the integral for each initial state polarization
integrate (eeww)
# The analysis setup fills the different histograms
analysis =
if (all Hel == -1 ["W+"] and all Hel == -1 ["W-"] ) then
if (all Hel == -1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_mm_mm (eval cos (Theta) ["W+"]) elsif
(all Hel == -1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_mp_mm (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_pm_mm (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_pp_mm (eval cos (Theta) ["W+"])
endif
endif; if (all Hel == -1 ["W+"] and all Hel == 0 ["W-"] ) then
if (all Hel == -1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_mm_m0 (eval cos (Theta) ["W+"]) elsif
(all Hel == -1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_mp_m0 (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_pm_m0 (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_pp_m0 (eval cos (Theta) ["W+"])
endif
endif; if (all Hel == -1 ["W+"] and all Hel == 1 ["W-"] ) then
if (all Hel == -1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_mm_mp (eval cos (Theta) ["W+"]) elsif
(all Hel == -1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_mp_mp (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_pm_mp (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_pp_mp (eval cos (Theta) ["W+"])
endif
endif; if (all Hel == 0 ["W+"] and all Hel == -1 ["W-"] ) then
if (all Hel == -1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_mm_0m (eval cos (Theta) ["W+"]) elsif
(all Hel == -1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_mp_0m (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_pm_0m (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_pp_0m (eval cos (Theta) ["W+"])
endif
endif; if (all Hel == 0 ["W+"] and all Hel == 0 ["W-"] ) then
if (all Hel == -1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_mm_00 (eval cos (Theta) ["W+"]) elsif
(all Hel == -1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_mp_00 (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_pm_00 (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_pp_00 (eval cos (Theta) ["W+"])
endif
endif; if (all Hel == 0 ["W+"] and all Hel == 1 ["W-"] ) then
if (all Hel == -1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_mm_0p (eval cos (Theta) ["W+"]) elsif
(all Hel == -1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_mp_0p (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_pm_0p (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_pp_0p (eval cos (Theta) ["W+"])
endif
endif; if (all Hel == 1 ["W+"] and all Hel == -1 ["W-"] ) then
if (all Hel == -1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_mm_pm (eval cos (Theta) ["W+"]) elsif
(all Hel == -1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_mp_pm (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_pm_pm (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_pp_pm (eval cos (Theta) ["W+"])
endif
endif; if (all Hel == 1 ["W+"] and all Hel == 0 ["W-"] ) then
if (all Hel == -1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_mm_p0 (eval cos (Theta) ["W+"]) elsif
(all Hel == -1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_mp_p0 (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_pm_p0 (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_pp_p0 (eval cos (Theta) ["W+"])
endif
endif; if (all Hel == 1 ["W+"] and all Hel == 1 ["W-"] ) then
if (all Hel == -1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_mm_pp (eval cos (Theta) ["W+"]) elsif
(all Hel == -1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_mp_pp (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == -1 [incoming "e-"]) then
record cta_pm_pp (eval cos (Theta) ["W+"]) elsif
(all Hel == +1 [incoming "e+"] and
all Hel == +1 [incoming "e-"]) then
record cta_pp_pp (eval cos (Theta) ["W+"])
endif
endif
# Generate events (keep the helicity information). As this takes
# some time, we request status information. The event sample has to
# be rebuilt as it is different for different initial state polarization.
simulate (eeww) {
?rebuild_events = true
?polarized_events = true
checkpoint = 10000
}
}
}
# Write out the generated histograms in eeww_polarized.ps / eeww_polarized.pdf
compile_analysis { $out_file = "eeww_polarized.dat" }
File Metadata
Details
Attached
Mime Type
text/plain
Expires
Mon, Jan 20, 8:41 PM (22 h, 51 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4242322
Default Alt Text
eeww_polarized.sin (10 KB)
Attached To
rWHIZARDSVN whizardsvn
Event Timeline
Log In to Comment