==========
Simulation
==========

-----------------
Detector modeling
-----------------

Related doc. for antenna model
`DB1591 <https://aradocs.wipac.wisc.edu/cgi-bin/DocDB/ShowDocument?docid=1591>`_,
`DB1693 <https://aradocs.wipac.wisc.edu/cgi-bin/DocDB/ShowDocument?docid=1693>`_,
`DB1838 <https://aradocs.wipac.wisc.edu/cgi-bin/DocDB/ShowDocument?docid=1838>`_

Related doc. for signal chain gain
`DB1179 <http://ara.physics.wisc.edu/docs/0011/001179/005/IceTemperatureRev3.pdf>`_,
`DB1188 <http://ara.physics.wisc.edu/docs/0011/001188/001/SystematicError.pdf>`_,
`DB2622 <https://aradocs.wipac.wisc.edu/cgi-bin/DocDB/ShowDocument?docid=2622>`_ ~ 2626,

In order to simulate realistic detector response,  detector model per station `A2 <https://github.com/MyoungchulK/MF_filters/tree/main/sim/ARA02/sim_table_full>`_ 
and `A3 <https://github.com/MyoungchulK/MF_filters/tree/main/sim/ARA03/sim_table_full>`_ is developed from data itself

Antenna model is based on pulser lift run performed at the South pole

Antenna model can be found at AraSim. 
`In_situ_VPol_Model.txt <https://github.com/ara-software/AraSim/blob/master/data/antennas/In_situ_VPol_Model.txt>`_, 
`In_situ_HPol_Model.txt <https://github.com/ara-software/AraSim/blob/master/data/antennas/In_situ_HPol_Model.txt>`_

Signal chain gain (electronic response) model is obtained from forced (software) trigger data and in-ice noise estimation

Detector model is produced by delow command::

    source ../setup.sh
    python3 script_executor.py -k rayl -s <station ID> -r <run number> -b 1              # produces noise and gain model

.. _fig-sim_hist0_0:
.. figure:: ../../figures/sim_hist0_0.png
    :align: center

    in-situ antenna model in dB scale

.. _fig-sim_hist0_1:
.. figure:: ../../figures/sim_hist0_1.png
    :align: center

    in-situ antenna model in linear scale

.. _fig-sim_hist3:
.. figure:: ../../figures/sim_hist3.png
    :align: center

    software trigger data to signal chain gain model

-------------------------
Signal / Noise generation
-------------------------

User can find setup file from here: `A2 <https://github.com/MyoungchulK/MF_filters/tree/main/sim/ARA02>`_ 
and `A3 <https://github.com/MyoungchulK/MF_filters/tree/main/sim/ARA03>`_

User can launch the simulation like executing below commands from :code:`AraSim/` path::

  source /home/mkim/analysis/MF_filters/setup.sh
  ./AraSim /home/mkim/analysis/MF_filters/sim/ARA02/sim_signal_setup_full/signal_E16_F1_A2_R1.txt 0 /data/ana/ARA/ARA02/sim_signal_full

Signal

- E^-1 spectrum between 10^7 ~ 10^12 GeV

- 8k events per 1) flavor, 2) decade, and 3) configuration 

- total # of event. A2: 840k, A3: 1080k

- Injection volume of each decade is following `previous analysis <https://github.com/ara-software/a23_4yr_diffuse_simconfigs>`_

- Each simulation set is tuned to follow each configuration in data sample section

Noise

- 1 mil. noise events per configuration 

- Recreating thermal noise event based on in-situ signal chain gain model and random distribution

.. _fig-sim_hist4:
.. figure:: ../../figures/sim_hist44.png
    :align: center

    sim summary plots

Signal generation supplement

Simulation is set to stop based on number of ‘triggered’ event
    
- It will continue to throw event until it gets desired number of ‘triggered’ event

- Thats why it has different thrown per decade and flavors, but similar number of triggered event

Simulation was run in E^-1 for each decade
    
- If It runs in E^-1 for all energy range at once, I will get only few low energy event

- Thats why triggered events have discontinuing shape, but I get decent amount of low energy event

Estd shape means estimated shape of thrown events.

- (Unfortunately) Sim didn’t save the energy of all thrown event. It only saved number of thrown events. So, I estimated E^-1 shape based on number of thrown events

.. _fig-sim_hist4_0:
.. figure:: ../../figures/sim_hist44_0.png
    :align: center

    A2 sim RZ plot in each energy decade and radius

.. _fig-sim_hist4_1:
.. figure:: ../../figures/sim_hist44_1.png
    :align: center

    A3 sim RZ plot in each energy decade and radius

--------------------
Effective Area Check
--------------------

.. _fig-sim_hist5:
.. figure:: ../../figures/sim_hist5.png
    :align: center

    A23 effective area in trigger level 

Effective area by triggered event

Based on new gain and threshold value, low energy has an improvement compare to previous analysis’s Aeff

A2 config 7 and A3 config 8 / 9 have lower area based on less # of trigger antenna at the time period