CHICOS Analysiscut2

CTShower Code

Running analysiscut2 on old data:

New Showers

I wanted to look at how well the automatic reconstructor was doing and whether it would be useful to put in more cuts on quality of fit. Three things to look at are how many sites are fit with non-zero intensity, and how far off in intensity and time the worst-fit hits in the shower are.

I had one problem in getting this info because I wanted to loop over all hits in the shower object and check the fit, but when the only hits at a site are above the charge threshold of 1000 MIP, this breaks. I removed showers with this problem from the file just for the purpose of getting some statistics, but it would be good to figure out a way around this.

Number of sites in the shower fit with >= 0.1 MIP:

Distance between time average and fit for worst fit site in the shower:

Distance between intensity average and fit for worst fit site in the shower:

Correlation between intensity/time offsets and number of fit sites:

Having a hit far off in intensity fit seems to correlate with large showers. On the other hand, a hit far off in time is more common in showers with < 4 sites fit with non-zero intensity. I would propose making a cut on showers with < 4 hit sites AND one hit more than 5 us away from the fit time. This would remove 8 showers from the best candidate list:

This doesn't actually seem like a very helpful cut. It might be good to cut out all showers with < 3 sites fit AND all showers with a site off by > 5 us.