icecube.finiteReco.segments_photospline module

icecube.finiteReco.segments_photospline.advancedSplineLengthReco(tray, name, inputPulses, inputReco, cylinderRadius=200, geometry='IC86', AmplitudeTable='', MaxRadius=inf, TiltTableDir='', TimingSigma=0.0, TimingTable='', If=<function <lambda>>)

Advanced length reconstruction. The reconstructed vertex and stopping point are found by a likelihood minimization using gulliver. Only hits within a cylinder with radius cylinderRadius around inputReco are considered. The simple length reconstruction is run first, and its output is used for seeding the advanced reconstruction. Note that the pulse map inputPulses should in most cases have a soft hit cleaning, e.g. classicRT.

This reconstruction puts four objects into the frame (note that the names always contain the name name of the tray segment, to allow multiple instances of the segment being run in the same script):

  • inputReco_Contained_name: An I3Particle, which has the direction of inputReco, but the newly reconstructed length. The position is set to the reconstructed start point.

  • inputReco_Contained_name_FiniteCuts: Some cut values related to the length reconstruction.

  • inputReco_Contained_name_StartStopParams: The likelihood values for a starting, stopping and infinite track. These values are supposed to be used as likelihood ratios, i.e. as LLHInf-LLHStart and LLHInf-LLHStop. All other uses on your own risk!

  • inputReco_Contained_name_FitParams: The gulliver FitParams for the length reconstruction.

icecube.finiteReco.segments_photospline.simpleSplineLengthReco(tray, name, inputPulses, inputReco, cylinderRadius=200, doSimpleOnly=True, geometry='IC86', AmplitudeTable='', MaxRadius=inf, TiltTableDir='', TimingSigma=0.0, TimingTable='', If=<function <lambda>>)

Simple length reconstruction. The reconstructed vertex and stopping point are found by simply projecting back from the first and last hit of inputPulses on the track inputReco. Only hits within a cylinder with radius cylinderRadius around that track are considered. Note that the input pulse map inputPulses should in most cases have a soft hit cleaning, e.g. classicRT.

This reconstruction puts three objects into the frame (note that the names always contain the name name of the tray segment, to allow multiple instances of the segment being run in the same script):

  • inputReco_Contained_name: An I3Particle, which has the direction of inputReco, but the newly reconstructed length. The position is set to the reconstructed start point.

  • inputReco_Contained_name_FiniteCuts: Some cut values related to the length reconstruction.

  • inputReco_Contained_name_StartStopParams: The likelihood values for a starting, stopping and infinite track. These values are supposed to be used as likelihood ratios, i.e. as LLHInf-LLHStart and LLHInf-LLHStop. All other uses on your own risk!