Odd-Low-Pass Filter for A Wave-Corrected Assessment of the Local Midplane

Software for the paper titled "A wave-corrected assessment of the local midplane", published in the Astrophysical Journal (ApJ). The paper proposes an Odd-Low-Pass filter (OLPF) that removes small, odd-parity structures (like the vertical waves) from the observation of the Milky Way, allowing for more accurate analyses of the Galaxy's midplane and vertical structures.

The Python code here encompasses the OLPF, Verlet Integration Simulation, and an analysis without the OLPF, as discussed in the paper. 

• R_cuts_and_CMD_trim.py

This script performs cuts on R ranges and color-magnitude diagrams (CMDs), generating _trimmed.csv files for each bin of R = 0.1 kpc. It outputs both the trimmed.csv files and an initial_cuts.csv. The latter includes data from all of the trimmed.csv files before dividing into separate bins. Pass the _trimmed.csv files generated here to the ''fileList'' in the OLPF.py. This should be the first step of the entire analysis.

• OLPF.py (errorbar.py)

File that filters the input data with the OLPF and produces the average vertical height of each radial bin with error bars in a figure. It also includes a method to investigate the effects of OLPF by subtracting the OLPF'd data from the raw data.

• Verlet Integration Simulation.py

Simulation using Verlet Integration that plots the motion of a selected number of stars in the vertical direction after every timestep of choice, considering gravity and background potentials.

• analysis without OLPF.py

Midplane analysis without OLPF that produces the average vertical height of each radial bin with error bars in a figure. The midplane analysis is the same as in the OLPF.py except the raw data is not OLPF'd.