View source for Search for CP violation in the gamma-ray sky

== Schema ==

[[File:illustration.jpg|500px]]




Illustration of the cut-sky with gamma rays distributed on it. Patches of radius R
degrees are centered on the highest energy gamma rays. In those patches we
test if the lower energy photons are distributed along left- or right-handed spirals. The
statistics Q measures the handedness averaged over all patches in the sky. For
more details, refer to "Search for CP violation in the gamma-ray sky", Hiroyuki Tashiro, 
Wenlei Chen, Francesc Ferrer, and Tanmay Vachaspati, [http://arxiv.org/abs/1310.4826 '''arxiv:1310.4826'''].

== Fermi Data ==

[[Fermi Data]]

== Galactic Latitude and Source Cuts ==

[[Galactic latitude and Source Cuts]]

== Evaluation of CP Odd Statistic (Q) ==


[[Q Evaluation routines]]

== Monte Carlo Simulations ==

[[MC Routines]]

== Results ==


[[File: Q_Rad_80_clean_jan.jpg|200px]]
[[File:Combinedjan14z100bgt80.jpg|600px]]

CP odd statistics, Q x 10<sup>6</sup>, versus patch radius R in degrees for the 6 different energy (E<sub>1</sub>,E<sub>2</sub>) combinations. E<sub>3</sub> is always taken to be 50 GeV. On the left are results of Tashiro's routines; on the right are consistent results obtained using Ferrer's routines. 1σ spreads of the Monte Carlo simulations are shown as magenta error bars. The data points are the results from Fermi data together with 1σ standard errors. Points that deviate by more than 2σ, where σ = max{Monte Carlo spread, standard error}, are colored in red.

== Progression ==

[[File:Cp-slide-sept-jan.pdf|800px]]

== Exposure ==

To account for the non-uniform exposure of Fermi-LAT we construct the time exposure map in the following way:

First create a livetime cube with ''gltcube'', then use ''gtexpcube2'' to obtain full sky exposure maps. These are used to generate Monte Carlo samples. Figure 2 in 1412.3171 shows exposure maps over weeks 9-328 using ''v9r33p0' Fermi tools.

The following bash script generates the exposure maps: [[File:Timecube.sh]]. The script calls a python routine that makes additional cuts on the exposure like e.g. galactic latitude > 50: [[File:Cutblat.py]] and [[File:Bcut.py]].

Some of the files generated by the shell script can be found in the file listings (accessed from 'Special pages'), but not the ''exp...fits'' ones, since they are too large.

For example, in the 10-20 bin, the script first generates [[File:Cube10_20GeVp202.fits]], which is used by gtexpcube to generate the (missing) exp file. These tools need the time information file [[File:Gtmktime10_20GeVp202.fits]] which is generated by ''gtmktime'' from the event file [[File:Gtseclect10_20GeVp202.fits]]. You also need the spacecraft file that contains information on where was the satellite at a particular time and that is provided by Fermi. With this exposure file one can select the relevant Fermi data as done in the ''cutblat.py'' script that creates e.g. [[File:Jan14z100e10_20GeVp202.fits]]

With this information you can generate Monte Carlos by weighting each direction according to the exposure as shown in this script [[File:Exposuremonte.py]]

== Contact Information ==

Wenlei Chen: wenleichen AT wustl.edu

Francesc Ferrer: ferrer AT physics.wustl.edu

Hiroyuki Tashiro: hiroyuki.tashiro AT nagoya-u.jp

Tanmay Vachaspati: tvachasp AT asu.edu