Simulation Tutorial

Define (import) amplitude (function) to simulate

The function will be use by the rejection method to "carve" a\ new distribution into the input simulated data read in next lines.

Read input (flat) simulated data (in condense format)

Read data full (from gamp files)

The format of the input data will depend on the Amplitudes: In this example the standard HEL angles, polarization angle (alpha) and other information neccesary are given for PWA (see below)

Produce simulation mask

A boolean file of (False and True)

Check on the waves and resonances that will be produced

This is a matrix with the weigths of each wave on each resonance

Check on how many events will be kept through the masking

Apply mask to input data

new_data will contain the simulated data in the same format that data/datag

Mask full data format

Write new_data into a Pandas Dateframe

Plot simulated data intensity versus mass

Plot phi_GJ vs costheta_GJ) of simulated data (with 4 different contracts(gamma))

Plot cos(theta_GJ) vs mass for simulated data (with 4 different contrasts)

Plot phiGJ vs mass for simulated data (with 4 different contrasts)

Histogram of alpha/Phi

alpha/Phi is the polarization angle

Plot mass versus alpha/Phi