CHICOS Analysiscut2

Running analysiscut2 on old data:

I am modifying the program simShowers3 to throw simulated showers over the whole CHICOS array. This program writes a file of showers which pass analysiscut2. This file can then be given to autoFit4 which reconstructs the showers and further filters the results. The output of autoFit can be read in root with the macro read_chicos_sim.C.

I ran a test that threw showers between 10^19 eV and 10^21 eV. It throws enough showers that 10 pass analysiscut2 in each bin. It threw 1974 showers alltogether, and took about 1/2 hour. Some plots:

This is the energy spectrum of showers thrown on the array.

This is the simulated energy of showers which passed all cuts. 10 showers per bin passed analysiscut2, but fewer than that passed all cuts.

This is the reconstructed energy of showers which passed all cuts.

This is a map of all showers thrown on the array with a map of the array (12/31/03) overlaid. Note that before 10/31/03 the array is not dynamic. All showers before that date use that site information file.

Map of input core locations of showers which pass the cuts:

Same, overlaid on map of showers thrown. May need to adjust the area over which showers are thrown. All but 2 showers with energy less than 10^20 were in the SG valley.

Map of output (reconstructed) core locations:

Reconstructed vs. simulated parameters:

Notes:
- The range of theta generated was 0.4 < costheta < 1, or 0 < theta < 66 degrees. This is larger than the allowed range of theta in the reconstructor.
- AutoFit4 should be changed to ignore unhealthy sites even if hit. Same should apply when reconstructing real showers, although it will mean we lose a few.
- A .health file will be written by autoFit4 for each shower generated at a specific time. Add line to delete them?