Running analysiscut2 on old data:
New ShowersI 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.
Test 1: 10^18 - 10^21 eV, 10 showers per bin (0.2 in Log(E)). After running for 1 hour and throwing 10000 showers in the first bin, nothing had passed. This is surprising because showes of this energy should at least be picked up when they land near Chiquita.
Spectrum of showers thrown (just to check that the generation is correctly using a 1/E^3 spectrum in each bin):
Core location of showers thrown:
Core location of showers thrown near Caltech:
There were nearly 50 showers thrown near the Chiquita array and none of them triggered. Earlier simulations using only Chiquita gave a much higher rate, more like 50% at that energy. Maybe the difference comes from using the CHICOS LDF vs the Chiquita LDF? Also, the Chiquita sites were not operational until halfway through 2003, so that will cut the rate in half.
Plotting the Chiquita LDF and the CHICOS LDF together shows they are very similar at 10^18 eV:
Edited to add: The high rate of showers being accepted before were accepted as chiquita showers, not chicos showers. That must explain the difference.
Test 2: 10^18.8 - 10^21 eV, 100 showers per bin. About 50000 showers were thrown on the array. This took about 5.5 hours, plus about 4 hours to run autoFit on the resulting shower file. 1200 showers passed analysiscut 2 (this is the 100 per bin) and 910 passed all cuts.
The y:x distribution of core locations reveals that showers are being thrown more densely in the SG valley. This is my mistake - they are being thrown evenly in x, but y was being chosen according to the range of each valley, resulting in a denser coverage of the SG valley. Will be fixed in the next iteration.
Fixed version throws showers evenly:
Other plots:
Core location of showers which passed:
San Gabriel array could possibly be given larger area.
At lower energies almost all showers are in the SG Valley:
Reconstructor spreads for showers in each input bin:
