icecube.lilliput.segments module

icecube.lilliput.segments.I3IterativePandelFitter(tray, name, fitname='', pulses='OfflinePulses', n_iterations=4, minimizer=None, parametrization=None, seeds=[], domllh='SPE1st', tstype='TFirst', noiserate=1e-08, If=None)

Run an iterative Pandel Fit.

Use standard settings as for IC40-IC79 and bulk ice properties.

Parameters:
  • fitname (str) – Fit name

  • pulses (str) – Name of pulse map

  • n_iterations (int, optional) – Number of iterations

  • minimizer (I3MinimizerBase or str, optional) – Minimizer service

  • parametrization (I3ParametrizationBase or str, optional) – Track parametrization service

  • seeds (list) – Seed names

  • domllh ({"SPE1st", "SPEAll", "SPEqAll", "MPE"}, optional) – Likelihood function

  • tstype (str, optional) – Use "TNone" if seed has a reliable vertex time.

  • noiserate (float, optional) – Noise rate

Returns:

Names of services and fit for re-usage convenience

Return type:

tuple

icecube.lilliput.segments.I3ParaboloidPandelFitter(tray, name, fitname='', pulses='OfflinePulses', minimizer=None, parametrization=None, inputtrack=None, domllh='SPE1st', input_tstype='TFirst', grid_tstype='TFirst', noiserate=1e-08, If=None)

Run a Paraboloid fit with a Pandel likelihood.

Use standard settings as used for OnlineL2 and bulk ice properties This tray segment only works for gulliver Pandel track fits.

Parameters:
  • fitname (str) – Fit name

  • pulses (str) – Name of pulse map as for input track

  • minimizer (I3MinimizerBase or str, optional) – Minimizer service

  • parametrization (I3ParametrizationBase or str, optional) – Track parametrization service

  • inputtrack (str) – Name of input track

  • domllh ({"SPE1st", "SPEAll", "SPEqAll", "MPE"}, optional) – Likelihood as for input track

  • input_tstype (str, optional) – Vertex time strategy as for input track

  • grid_tstype (str, optional) – Use "TFirst" or "TChargeFraction" for grid.

  • noiserate (float) – Noise rate as for input track

Returns:

Names of services and fit for re-usage convenience

Return type:

tuple

Notes

The default vertex time strategy of "TFirst" for the input track was chosen because it is the standard for online level 2 and level 3. Supposing the input track has already an excellent timing, "TNone" should be the logical and optimal choice. According to some incomplete testing there should not be any dramatic differences between the two strategies.

icecube.lilliput.segments.I3SinglePandelFitter(tray, name, fitname='', pulses='OfflinePulses', seeds=[], minimizer=None, parametrization=None, domllh='SPE1st', tstype='TFirst', noiserate=1e-08, If=None)

Run a single Pandel fit.

Use standard settings as for IC40-IC79 and bulk ice properties.

Log-likelihood functions:

  • SPE1st: use time of first pulse in each DOM.

  • SPEAll: use times of all pulses in each DOM.

  • SPEqAll: use charge-weighted times of all pulses in each DOM.

  • MPE: use time of first pulse and total charge in each DOM.

Parameters:
  • fitname (str) – Fit name

  • pulses (str) – Name of pulse map

  • minimizer (I3MinimizerBase or str, optional) – Minimizer service

  • parametrization (I3ParametrizationBase or str, optional) – Track parametrization service

  • seeds (list) – Seed names

  • domllh ({"SPE1st", "SPEAll", "SPEqAll", "MPE"}, optional) – Likelihood function

  • tstype (str, optional) – Use "TNone" if seed has a reliable vertex time.

  • noiserate (float, optional) – Noise rate

Returns:

Names of services and fit for re-usage convenience

Return type:

tuple

Warning

Although SPEAll and SPEqAll were thought to be fundamentally correct, they work surprisingly worse than SPE1st.

Notes

The default noise rate of 10Hz does not have any physical basis. It was set at the time when we used only HLC pulses and based on flawed logic (the HLC dark noise rate was estimated to be of order 1Hz and then 10Hz seemed sort of conservative). Studies on low filtering level indicate that higher noise level reduce the rate of mis-reconstructed background, while at high cut level the effect of this setting on well reconstructed neutrino events seems to be minor. New studies and developments are underway (as of November 2011).

icecube.lilliput.segments.add_minuit_simplex_minimizer_service(tray, minimizer=None)

Add Minuit (Simplex) minimizer service to tray.

We need only one Minuit minimizer service in the tray, this can be used by all gulliver fitters that need it.

Parameters:

minimizer (I3MinimizerBase or str, optional) – If not None, use external minimizer service.

Returns:

Minimizer service; default: default_simplex

Return type:

I3MinimizerBase or str

icecube.lilliput.segments.add_pandel_likelihood_service(tray, pulses, domllh='SPE1st', noiserate=1e-08)

Add Pandel service to tray.

For every desired Pandel configuration we need at most one instance of the Pandel service in the tray. This service can be used by all gulliver fitters that need it.

With this function you can do some configuration, but not all. If you want to do more fancy configuration, then you’ll need to add the service yourself instead of using this convenience function.

Parameters:
  • pulses (str) – Name of pulse map

  • domllh ({"SPE1st", "SPEAll", "SPEqAll", "MPE"}, optional) – Likelihood function

  • noiserate (float) – Noise rate

Returns:

Name of Pandel service

Return type:

str

icecube.lilliput.segments.add_seed_service(tray, pulses, seeds=[], tstype='TFirst')

Add seed service to tray.

For every desired seeding configuration we need at most one instance of the basic seed service in the tray. This service can be used by all gulliver modules that need it.

Parameters:
  • pulses (str) – Name of pulse map

  • seeds (list) – Seed names

  • tstype (str, optional) – Use "TNone" if seed has a reliable vertex time.

Returns:

Name of seed service

Return type:

str

icecube.lilliput.segments.add_simple_track_parametrization_service(tray, parametrization=None)

Add track parametrization service to tray.

We need only one track parametrization service in the tray, this can be used by all gulliver fitters that need it.

Parameters:

parametrization (I3ParametrizationBase or str, optional) – If not None, use external track parametrization service.

Returns:

Track parametrization service; default: default_simpletrack

Return type:

I3ParametrizationBase or str