aboutsummaryrefslogtreecommitdiff
path: root/src
AgeCommit message (Expand)Author
2019-09-09add math header file to zdab_utils.c for fmin() and fmax()tlatorre
2019-09-09update get_event() to use the PMT types from the snoman.ratdb filetlatorre
2019-09-09use EHS instead of QHS in is_muon()tlatorre
2019-09-09lower QLX threshold from 4000 -> 1000 in is_flashertlatorre
2019-09-09reset the PMT_FLAG_CHARGE bit in get_event()tlatorre
2019-09-09update fit and zdab-cat to skip ZDAB recordstlatorre
2019-09-09add the --gtid command line option to zdab-cattlatorre
2019-09-09uncalibrate MC times and charges in get_event()tlatorre
2019-09-09switch to using ept instead of pt1 in is_flasher() and is_slot_early()tlatorre
2019-09-09only loop over normal PMTs in is_slot_early()tlatorre
2019-09-09fix a typo in is_slot_early()tlatorre
2019-08-28fix some error handling in zebra.ctlatorre
2019-08-27fix the fts cuttlatorre
2019-08-27add the in time channel data cleaning cut from SNOtlatorre
2019-08-27update neck cuttlatorre
2019-08-26fix a bug in get_{shower,delta_ray}_weights() causing crashtlatorre
2019-08-26sort particle combo arraytlatorre
2019-08-05update sno_charge to prevent nan when charge is negativetlatorre
2019-08-05add ability to specify a particle combo on the command linetlatorre
2019-07-29use standard PCA time for MCtlatorre
2019-07-29update run-fit to write output files with suffix hdf5tlatorre
2019-07-29write out the hdf5 file after every fittlatorre
2019-07-29flag PMT charges below qlotlatorre
2019-07-29don't need N_ACOS anymoretlatorre
2019-07-29add a faster version of fast_acos()tlatorre
2019-07-29avoid a division in get_expected_charge()tlatorre
2019-07-29fast_sqrt -> sqrttlatorre
2019-07-29switch to using the multiphoton PCA timetlatorre
2019-07-29fix a bug in theta0 min calculationtlatorre
2019-07-29ev.gtid -> bev.gtr_idtlatorre
2019-07-16fix bug introduced in ebe2799tlatorre
2019-07-16don't double count PMT pairs in fts cuttlatorre
2019-07-16update neck tube cut to include time difference changes due to cable changestlatorre
2019-07-16delete some unused #defines and add a comment to the PSUP radiustlatorre
2019-07-16use QLX if QHS is railedtlatorre
2019-07-16multiply rayleigh scattering by 2 to account for the 50% PMT coveragetlatorre
2019-07-12don't load DQXX file for run 10000 by defaulttlatorre
2019-07-12don't flag PMT hits based on the uncalibrated QHS for MCtlatorre
2019-07-11switch from YAML output to HDF5 to speed things uptlatorre
2019-07-08add vertex time field to zdab-cattlatorre
2019-07-05add MCVX time to the YAML filetlatorre
2019-06-25add ability to read gzipped zdab filestlatorre
2019-06-23update how neck events are flaggedtlatorre
2019-06-23update flasher cuttlatorre
2019-06-20update zdab-cat to output 50 MHz clock time and trigger typetlatorre
2019-06-20fix cat-grid-jobs againtlatorre
2019-06-20update zdab-cat to emit multiple YAML documentstlatorre
2019-06-20fix empty list at top of YAML output in zdab-cattlatorre
2019-06-19add data cleaning word and ftp, ftk, and rsp info to zdab-cat outputtlatorre
2019-06-19add FTP, RSP, and FTK results to the output filetlatorre
the range - now only fit to the end of the PSUP ^ add a time limit for the fit ^ multi-particle fit p optimize CHARGE_FRACTION ? find out why likelihood is sometimes returning nan ^ add pdf for cerenkov light from delta rays ^ tweak find peaks algorithm ^ update zebra to allow use of links p minimum energy for each particle type - I do already have a minimum energy which is the Cerenkov threshold, however I think what I meant here was to have something significantly higher than that since a particle just at the Cerenkov threshold won't produce many PMT hits. I think I was mostly concerned about this issue because many fits were always fitting better by including a low energy electron, however I have solved that problem by simply dropping multi particle fits which have an electron with an energy less than 20 MeV. ^ check all zebra links ^ fix zero logical record size bug in zebra.c p speed up fit - have optimized most of the likelihood function. I think further optimization will require a different technique for evaluating the integrals ^ nhit cut p add flasher cut p default max particles - submit-grid-jobs fits up to 3 particles by default ^ add environment variable to control where to look for pmt.txt and titles files p update charge code to handle highest values for qhs and qlx - don't have this coded into the charge probability code, but I now use QHL if QHS is railed and flag any PMT with a railed QHL ^ use qhs and qlx based on variable in PMT bank - Note: I don't use the best charge variable in the PMT bank, but I do select QHL when QHS is railed ^ fix bug in path_init() ^ load DQXX file based on run number ^ update rayleigh scattering lengths ^ fix delta ray calculation p figure out why electron energy bias is +10% - updated code to use a more accurate approximation for the number of shower photons which seems to bring the bias down to 5% p update find peaks algorithm to do single particle quick fits - I did update the find peaks algorithm to work *much* better, but it still doesn't do single particle quick fits. I think the next improvement would be to add the ability to actually determine the number of rings. ^ figure out how to combine SNO fitter data with my fitter for final analysis. For example, how to apply neutron follower cut? - double check that zdab-reprocess is correct for the D2O and salt phases since it appears to be from the NCD phase ? add code to compute expected deviation from nll_best to normalize psi - tried several different versions of this and nothing seemed to perform as well as psi/nhit. ^ add term to likelihood for probability that a channel is misccalibrated or didn't make it into the event - when calculating the first order statistic for a gaussian is alpha = pi/8 a better choice than 0.375? - speed up charge initialization - determine *real* mean and standard deviation of single PE charge distribution. TELLIE SNO+ data? ^ extend electron range table - extended the electron range table up to 1 TeV ? fix time PDF. Currently take the first order statistic of dark noise + indirect light + direct light all together, but this isn't correct. - thought more about this and it's not actually obvious if what I'm doing is wrong. Need to think more about this.