Chiquita data: iterative fit - Jan. 29

The goal here is to compare the various data sets taken from the chiquita data. It is probably most interesting to first compare the original set of 500 showers with the set of showers found by analysiscut2 during the same time period.

The aperture histogram for the first set is the most recent one, (using analysiscut2) and takes into account a cutoff of theta <= 45 degrees.

Aperture from simulations:

Iterative fits:

• Original 500 showers.
• Same time interval, 5295 showers found by analysiscut2.
• Total chiquita data set, 16067 showers found by analysiscut2.

Results:

Data histogram:

• 500 showers
• 5000 showers
• 15000 showers

Data x E^3 histogram:

• 500 showers
• 5000 showers
• 15000 showers

Data x E^3, iterated:

• 500 showers
• 5000 showers
• 15000 showers

The aperture histogram for the second set is taken from the set of simulations used for the spectrum in the ICRC paper, but in order to be consistent with the above, I have added a cutoff of theta <= 45 degrees. I wouldn't trust the low end of rhist because the simulated energies only extend down to 10^16.4.

Aperture from simulations:

Iterative fits:

• Original 500 showers.
• Same time interval, 5295 showers found by analysiscut2.
• Total chiquita data set, 16067 showers found by analysiscut2.

Results:

Data histogram:

• 500 showers
• 5000 showers
• 15000 showers

Data x E^3 histogram:

• 500 showers
• 5000 showers
• 15000 showers

Data x E^3, iterated:

• 500 showers
• 5000 showers
• 15000 showers

• 500 showers, temperature-corrected, no cut on theta, older version of the iteration code, hand-filtered for accidentals.
• All chiquita data, not temperature-corrected, no cut on theta, older code, not filtered for accidentals. <- has errors
• 500 showers, temperature-corrected, with a cut on theta <= 45 degrees, older code, hand-filtered.
• All chiquita data (post 06/10/04), not temperature-corrected, not filtered, older code, includes cut on theta.

One possible problem is that chiquita showers are being filtered into the chicos list... I looked through part of the analysiscut2 list and found quite a few that would qualify as chiquita showers (most are not very high energy, though).