magnetic monopoles are particles carrying a quantized magnetic charge and are predicted in various theories . in classical electrodynamics , their existence would symmetrize maxwell s equations with respect to the sources of the electromagnetic field .
quantum mechanically , the existence of magnetic monopoles implies that both electric charge and the hypothetical magnetic charge , are quantized , given that the associated electromagnetic fields still satisfy maxwell s equations @xcite . the resulting magnetic elementary charge , called the dirac charge @xmath5 , is    @xmath6    where @xmath7 is the electric elementary charge and @xmath8 is the fine structure constant .    in grand unified theories ( guts )
@xcite magnetic monopoles appear as stable , finite energy solutions of the field equations @xcite . the predicted
masses range from @xmath9 to @xmath10 @xcite and their magnetic charges are integer multiples of the dirac charge @xmath11 .
the lower part of the mass range up to @xmath12 refers to intermediate mass monopoles ( imms ) which arise from intermediate stages of symmetry breaking below the gut scale .
in contrast the superheavy monopoles with masses at the gut scale may have been created during the phase transition associated with the spontaneous breakdown of the unified gauge symmetry in the early universe at @xmath13 after the big bang @xcite . the monopole mass and
charge depend on the underlying gauge group , the symmetry breaking hierarchy , and the type and temperature of the phase transition in a particular gut
.    since magnetic monopoles are stable , they should still be present in cosmic rays .
the number density today depends on the existence of an inflationary epoch and on the time of creation , which could be before , during or after this epoch @xcite .
since then , monopoles have been accelerated by large - scale cosmic magnetic fields .
the kinetic - energy gain by passing through a magnetic field @xmath14 is given by    @xmath15    where @xmath16 is the magnetic charge .
the maximum kinetic energy of a magnetic monopole due to acceleration in cosmic magnetic fields is rather uncertain but can reach @xmath17 @xcite .
therefore , monopoles with masses at , or above , this energy scale should be non - relativistic . based on the propagation of magnetic monopoles in the galactic magnetic field an upper bound on the monopole flux
can be calculated , assuming the galactic magnetic field does not decrease faster than it can be regenerated .
this assumption constrains the monopole flux to be less than @xmath18 , which is called the parker bound @xcite .
taking into account the fields during galaxy formation , the limit was extended by adams et al . to be less than @xmath19 for monopoles with masses below @xmath10 @xcite .
many experiments have searched for relic magnetic monopoles , but there is no experimental proof for their existence .
the current best limits for magnetic monopoles constrain their flux to a level of @xmath20 depending on the monopole speed and interaction mechanism @xcite .
consequently , searches for magnetic monopoles require very large detectors .
the icecube neutrino observatory currently is the world s largest neutrino detector .
the primary goal is the detection of cherenkov light from electrically - charged secondary particles produced when high - energy astrophysical neutrinos interact in the surrounding matter @xcite .
however , icecube can also be used to search for magnetic monopoles .
depending on their speed monopoles have different signatures in icecube .
relativistic monopoles with a speed above the cherenkov threshold , e.g. @xmath21 in ice , can be detected by the cherenkov light they directly produce @xcite .
non - relativistic monopoles that catalyze the decay of nucleons in the detector medium can , in contrast , be detected by the cherenkov light from electrically charged secondary particles produced in subsequent nucleon decays along the monopole trajectory ( sec .
[ subsec : rubakov - callan_effect ] ) . therefore ,
different analysis strategies are needed in order to cover both detection channels .
this paper presents the results of a search for non - relativistic magnetic monopoles which would catalyze the proton decays via the rubakov - callan effect in icecube .
the icecube neutrino observatory consists of the in - ice detector , icecube , and the surface air shower detector , icetop .
it is located at the geographic south pole . for the in - ice
detector @xmath0 of the antartic ice , which is used as detection medium , has been instrumented .
the detector consists of 86 strings equipped with 60 digital optical modules ( doms ) each .
the dom , the sensor of the icecube detector , consists of a glass pressure housing enclosing a @xmath22 diameter hamamatsu photomultiplier tube ( pmt ) with the electronics needed for signal digitization , and a set of leds for calibration purposes @xcite .
signals that pass a threshold of about @xmath23 photo - electrons are digitized and recorded .
this process is called a dom launch or for simplicity a _ hit _ in the following .
two hits are labeled as hard local coincidences ( also called hlc pair ) , if their time difference is less than @xmath24 and the corresponding doms are nearest or next - to - nearest neighbors on the same string .
the recorded data is sent to the surface and a trigger algorithm evaluates the time and position of the hits and decides whether they form an event .
for example , for relativistic particles a simple multiplicity trigger requiring eight hlc hits within a time window of @xmath25 , called smt8 , is used .
the doms are deployed at depths between @xmath26 and @xmath27 @xcite . at depths below @xmath28 ,
eight inner strings are placed with smaller separations from each other and thus form a region of denser instrumentation . together with seven central standard strings they form deepcore , a low - energy sub - detector @xcite .
the construction of icecube was completed december 16 , 2010 but data taken during intermediate construction stages were already used for physics analyses during earlier years .
one of the two presented analyses uses data taken from may 2009 to may 2010 , when icecube was operating in its 59-string configuration ( ic-59 ) .
the other analysis uses the fully installed detector .
non - relativistic magnetic monopoles would themselves be too slow to emit cherenkov light when propagating through the icecube detector .
however , relativistic charged secondary particles , produced in monopole interactions with the surrounding matter , can produce cherenkov light and thus can be detected by the icecube detector .
the energy loss of a magnetic monopole due to ionization can be described by a modified bethe - bloch formula @xcite , which is valid for speeds @xmath29 . for lower speeds in the range from @xmath30 to @xmath31 ahlen and
kinoshita introduced a model to calculate the energy loss of magnetic monopoles @xcite .
later , ritson extended this model for speeds below @xmath30 @xcite . for magnetic monopoles with e.g. @xmath30
the energy loss is of the order of @xmath32 @xcite .
only electrons above the cherenkov threshold of @xmath33 kinetic energy emit detectable cherenkov light . however , the maximum transferred energy of a monopole with e.g. @xmath30 on an atomic electron is typically @xmath34 and no cherenkov light is produced .
rubakov @xcite and callan @xcite showed that some monopoles could catalyze nucleon decays along their trajectories ( rubakov - callan effect ) .
this effect depends on the gauge group of the respective gut theory @xcite and on assumptions , e.g. on the fermion masses or the relative velocity between the quarks and the monopole , used in the calculation @xcite .
furthermore , this process is not possible for intermediate mass monopoles with masses below @xmath12 @xcite .
therefore , the sensitivity of this analysis is contrained to heavier monopoles ( gut scale ) .
figure [ fig : protondecay ] illustrates the catalyzed decay of a proton by a gut monopole into a positron and a neutral pion :    @xmath35    for this decay channel almost the full rest mass energy of the proton is transferred to electromagnetic particles . because of the high light yield this channel is used as a benchmark in the analyses .
the catalysis cross section for nucleon decays @xmath36 depends not only on the cross section @xmath37 @xcite , but also on the monopole speed @xmath38 :    @xmath39    the correction @xmath40 takes into account an additional angular momentum of the monopole - nucleus - system and becomes relevant for speeds below the speed threshold @xmath41 .
depending on the sign of @xmath42 the catalysis cross section is enhanced or suppressed . both parameters @xmath42 and @xmath41 depend on the nucleus @xcite .
current estimates for the catalysis cross sections are of the order of @xmath43 to @xmath44 @xcite .
the rubakov - callan effect results in small electromagnetic or hadronic cascades from catalyzed nucleon decays along the monopole trajetory through the detector .
this is illustrated in fig .
[ fig : monopole_signature ] .
experimentally , the relevant parameter is the mean free path @xmath45 between two decays .
that is    @xmath46    where @xmath47 is the particle density of the medium through which the monopole propagates .    the energy of each cascade , and therefore the number of emitted cherenkov photons , depends on the decay channel ( e.g. eq . [ eqn : promising_decay_channel ] ) .
a general quantity is the track length , @xmath48 , that a relativistic particle carrying a single electric charge would have to travel in order to emit the same number of cherenkov photons as the average number expected from a proton decay , @xmath49 @xcite . using this track length per proton decay , @xmath48
, the monopole s mean free path @xmath45 can be converted into the light yield per monopole track length @xmath50 .
@xmath51    a monopole with @xmath52 will therefore produce the same number of cherenkov photons per track length as a single - electric - charge , relativistic particle without stochastic energy losses along its track @xcite .
this implies that @xmath50 can be used to express the resulting monopole flux limits without assuming a specific decay channel ( sec .
[ sec : flux_limits ] ) .
this ansatz is valid as long as the monopole s light emission can be approximated as being continuous .
this condition is satisfied for a mean free path much smaller than the detector spacing . from an experimental point of view the speed @xmath53 and
the mean free path @xmath45 are the characterizing parameters for the detection of such monopoles .
are symbolized by orange rays . ]
searches for slow monopoles based on the rubakov - callan effect have been pioneered with the underground detectors imb and kamiokande - ii @xcite and the underwater detectors in lake baikal @xcite .
a similar search has also been performed with amanda , the predecessor of icecube @xcite .    during the commissioning of the full detector ( ic-86 ) in may 2011 , a dedicated trigger for slow particle signatures ( slow - particle - trigger , sec .
[ sec : slop ] ) in deepcore was implemented .
the denser instrumentation of deepcore allows icecube to detect monopoles of low light emission , i.e. with rather large values of mean distances @xmath45 between induced catalysis points . in 2009 , the deployment of the first deepcore strings was still ahead . due to the larger spacing and the lack of an appropriate trigger , ic-59 was blind for large @xmath45 . for smaller @xmath45
the mentioned drawbacks were balanced by the larger geometrical area compared to deepcore .
the signal expectation was determined from monte carlo simulations of magnetic monopoles in icecube , while the background expectation was determined from experimental data itself , with only supplementary simulations .
icecube simulation includes particle injection and propagation , taking into account appropriate particle interactions , as well as the full detector response to the generated cherenkov photons .
the arrival directions of magnetic monopoles are assumed to be isotropic .
the starting points of simulated monopole tracks are generated randomly on a disc of fixed size .
the distance of the plane is fixed with respect to the deepcore detector but its orientation is random .
it is assumed that the magnetic monopoles are not substantially decelerated along their track and their velocity is constant @xcite .
the distances between the catalyzed nucleon decays are simulated as a poisson process with a mean free path @xmath45 along the monopole track .
each nucleon decay is simulated as an electromagnetic cascade with an energy of @xmath54 , corresponding to the benchmark detection channel ( eq . [ eqn : promising_decay_channel ] ) . the simulation and propagation of the cherenkov light from these cascades
is done with the software package _ photonics _ @xcite using the ice model described in @xcite for the ic-59 analysis and an improved version described in @xcite for the ic-86/deepcore analysis .
background noise in the doms has to be superimposed on the signal .
this noise consists of uncorrelated random noise , mostly from radioactive decays in the doms and correlated noise because of after pulses and signals from atmospheric muons . for the ic-59 analysis ,
the random noise is simulated as a poisson process and the atmospheric muons are simulated using the software package corsika @xcite based on a 5-component model for cosmic rays with the hadronic interaction model sibyll @xcite and the hrandel flux model @xcite . for the simulation of noise in the ic-86/deepcore analysis
the detector response of simulated monopole signals is superimposed with random and correlated noise hits from experimental data .
these noise hits were recorded with a fixed rate trigger ( frt ) that was implemented to measure and analyze background noise in the detector .
more details on the frt data are given in sec .
[ subsec : background_estimation ] .
figure [ fig : simulated_monopole ] shows a simulated monopole event with @xmath30 and @xmath55 . because of the low speed , the event duration for a monopole is typically a factor of 1000 longer than for muon events and a large number of noise hits are recorded .
however , the monopole also produces a large amount of cherenkov light in the detector .
therefore , its signature can be separated from the randomly distributed noise hits already by eye .     and @xmath55 with superimposed background noise .
the black line represents the monopole track .
the doms are shown as tiny black dots .
the color code illustrates the time scale from red for early times to blue for later times .
the radii of the colored spheres scale with the number of recorded photoelectrons . ]
the experimental data set was recorded between may 2011 and may 2012 with a dedicated slow particle trigger applied to deepcore . in this period
the livetime of the detector was 351 days , with a total number of approximately 50 million triggered events .
multiple icecube triggers are implemented in the software of the data acquisition system @xcite .
most of them are sensitive to signatures of relativistic particles , e.g. muons , so they have little sensitivity to non - relativistic magnetic monopoles . only for the case of very bright magnetic monopoles
the large amount of light can frequently prompt triggers for relativistic particles .
this case is described in sec .
[ sec : ic59set ] .    the slow - particle - trigger ( slop trigger )
was first implemented in may 2011 @xcite .
for the first year , the trigger operated only on the subdetector deepcore . since may 2012
, the trigger has been operating on the full icecube detector .
are removed .
c ) all combinations of three hlc pairs , called triplet , with a time difference @xmath56 $ ] between two pairs are built .
d ) the cuts on the quality criteria @xmath57 and @xmath58 remove two more triplets . if the remaining triplets overlap in time and fulfill @xmath59 , a trigger is generated and the full detector data within the time span from the first to the last hlc pair of the triplets is recorded .
@xcite ]    the slop - trigger searches for time isolated local coincidences in nearby doms caused by subsequent nucleon decays along the monopole trajetory .
these coincidences have to be consistent with a straight particle track of constant speed .
the slop - trigger is illustrated in fig .
[ fig : triplet_selection ] .
specific values for the different trigger parameters are listed in table [ tab : trigger_conditions ] .
it is based on local coincidences of hits ( hlcs , sec .
[ icecube ] ) . for the trigger ,
the position and time , defined by the first hit of the hlc pair , of all hlc pairs are stored in a list ( fig .
[ fig : triplet_selection]a ) .
since muons pass the detector within @xmath60 they produce several hlc pairs within a short time . by removing all hlc pairs with time differences @xmath61 from the list , muon hits are efficiently rejected ( fig .
[ fig : triplet_selection]b ) .
the remaining hlc pairs are searched for every combination of three hlc pairs , the triplets ( fig .
[ fig : triplet_selection]c ) .
the time difference between any two hlc pairs within a triplet has to be in the range @xmath62 $ ] . furthermore only triplets that match a track - like signature are kept .
therefore two quality criteria are required : the contributing hlc pairs have to be ordered along a line and the time differences have to be consistent with a constant speed ( fig . [
fig : triplet ] ) .
the first can be verified by the parameter @xmath63 . if @xmath64 all hlc pairs are located on a line .
the second can be checked by the parameter    @xmath65    where @xmath66 with @xmath67 corresponds to the speed between the j - th and the i - th hlc pair within a triplet . for a monopole with a constant speed
all hlc pairs should be connected by a constant speed and therefore @xmath68 should be valid .
all triplets not satisfying these quality criteria are removed from the set of triplets ( fig .
[ fig : triplet_selection]d ) .     and the time @xmath69 of the first hit of an hlc pair .
the trigger observables are the distances @xmath70 and time differences @xmath71 . ]    finally , if the number of triplets in the set overlapping in time , @xmath72 , is greater than a minimum number of triplets @xmath73 , the trigger is launched . when these conditions are met , the full detector data from the first to the last hlc pair in the list of triplets are stored , also including those dom signals not contributing to the trigger .
the maximum event duration of the trigger is restricted to @xmath74 .
ll parameter & value + @xmath75 & @xmath76 + @xmath77 & @xmath78 + @xmath79 & @xmath80 + @xmath57 & @xmath81 + @xmath58 & @xmath82 + @xmath73 & @xmath83 + @xmath84 & @xmath85 +      to investigate the characteristics of the slop events , we use an experimental data set of @xmath86 days of live - time .
this is sufficiently short to exclude a significant signal contamination given by current flux limits ( sec .
[ sec : intro ] ) and hence the data can be considered as background .    .
the maximum is at about @xmath87 . for comparison the event duration distributions of the generated background events ( black )
is superimposed .
the number of entries is normalized to one . ]
figure [ fig : trigger_length ] shows the distribution of event durations of slop triggered events .
typical durations are of the order of milliseconds , whereas other icecube triggers have typical durations of a few microseconds .
figure [ fig:2days_ntriplet ] compares the @xmath72 distribution of the experimental data sample with simulated monopoles of @xmath30 and @xmath55 .
while the background distribution decreases rapidly for larger @xmath72 , the signal distribution is almost flat .
therefore , the quantity @xmath72 discriminates well between signal and background events .
the exponential decrease of the background distribution indicates a possible poissonian random process for combinations of hlc pairs which result in a triplet .
distribution of the experimental test data sample ( blue ) in comparison to a distribution of simulated monopoles with @xmath30 and @xmath55 ( green ) .
in addition , an exponential function is fitted to the tail of the experimental distribution for @xmath88 ( red ) . ]    to understand the underlying random processes for background events we developed a method to generate a high statistics sample of background events by reshuffling experimental events recorded with the frt .
the frt fires at fixed time intervals ( e.g. , every thirty seconds ) , and dom data from the entire detector are recorded over a time interval of 10ms .
the resulting events contain all types of random and correlated backgrounds , and highly unlikely any signal .
the frt events of @xmath89 length were split into snippets of @xmath90 , which were then randomly re - ordered to form new @xmath89 events .
the newly assembled events are then passed to the slop trigger algorithm ( fig .
[ fig : reshuffling_method ] ) .
this way , a total of 400 seconds of frt data were re - shuffled to generate a background sample of about 25 days of live - time equivalent .
the generated sample closely resembles the experimental slop - triggered events .
figure [ fig : trigger_length ] compares the event duration of the generated background events to the slop - triggered events in two days of experimental data .
the method reproduces the measured event duration distribution resonably well over several orders of magnitude . for shorter event duration
the distribution of the generated data sample tends to be below the distribution of the experimental test data sample .
this is expected because this method can not correctly model noise hits that are correlated over time scales larger than the length of the @xmath90 snippets .
below @xmath90 the triplets are characterized by the same dom combinations due to the low statistics of the frt events .
the overall good agreement indicates that correlated noise is a subdominant effect and is only relevant for short time - scales .
we will presume later that different triplets due to correlated noise are themselves based on largely independent sets of hlc pulses .
therefore , for large values of @xmath72 the contribution from correlated noise triplets is added as a random process similar to the triplets from uncorrelated noise .    .
they are split into @xmath90 snippets .
the snippets are shuffled randomly to build new @xmath89 frames . then the slop trigger algorithm is applied . ]
figure [ fig : ntriplet_cleaned_frt ] compares the @xmath72 distributions of experimental data and generated background .
overall both distributions are similar and show an exponential decay .
the differences can be understood by two effects .
the first is the aforementioned effect that noise correlations over time scales longer than @xmath90 are not taken into account , which is expected to increase the number of triplets . by removing triplets which arise from hlc - pairs fulfilling the typical time - scale of correlated noise ( @xmath91 or @xmath92 )
the agreement improves .
however overall correlated noise has only a small effect on these distributions .
more importantly , the frt data and the slop test - data do not correspond to the same data taking period .
the dom noise rate shows slow slight drifts over long periods of time .
the chance probability of producing large n - triplet values depends on this random noise .
this effect is accounted for by the background fit described in the following section . in conclusion ,
the observed background is understood by the noise characteristics of the doms .
distributions of the experimental test data set ( green ) and the generated background events ( black ) .
triplets caused by hlc pairs fulfulling @xmath91 or @xmath92 are not taken into account ( cleaned ) . ]      as result of the findings in the previous section , the @xmath47-triplet distribution for background is estimated by fitting the experimental data with a simple probabilistic model .
the generic ansatz assumes that the probability to find a triplet ( 3 hlc - pairs ) can be described with a combinatorial model . for a number of @xmath93 hlc pairs the maximum number of possible triplets
is given by @xmath94 .
if the probability @xmath95 for any three out of @xmath93 hlc pairs to build a triplet is approximately constant , then the probability for @xmath47 triplets out of @xmath96 possible triplets follows the binomial distribution :    @xmath97    as the hlc - pairs themselves arrise from random noise , the probability to observe @xmath93 hlc pairs is given by a poissonian :    @xmath98    where @xmath99 is the mean expectation for the number of hlc pairs @xmath93 in a given time window . the total probability to observe the number @xmath47 triplets
is given by a sum over all binomial probabilities @xmath100 weighted with the probability to observe @xmath93 hlc pairs .
this results in    @xmath101    the sum starts with the minimum number of hlc - pairs @xmath102 which are combinatorially required for @xmath103 triplets .
this is given by the solution of the equation @xmath104 .
here , @xmath105 is the greatest integer less or equal to @xmath106 .
the parameter @xmath107 is the total normalisation of @xmath108 .    with this ansatz it is possible to describe the distribution of @xmath72 with only three parameters @xmath107 , @xmath99 and @xmath95 .
figure [ fig : ntriplet_fit ] shows the fit of this model to two normalized , experimental @xmath72 distributions which are based on slop data corresponding to different noise rates . since the distributions are normalized only @xmath99 and @xmath95 have to be fit .
the background model well describes the @xmath72 distributions over several orders of magnitude .
moreover the increase in the noise rate is reflected in an increase of the value of @xmath99 which depends on the noise rate . in summary
it can be confirmed that the background events from the slop trigger are dominated by random noise .
distributions of experimental slop data .
the blue distribution corresponds to a @xmath109 higher noise rate than the red one .
both distributions are normalized to one .
the solid lines show the results of the background model fit and the fit parameters @xmath99 and @xmath95 are shown in the boxes . ]
the analysis searches for monopoles from all directions . also the random background is largely isotropic and the event selection does not depend specifically on the direction of the monopole .
however , an important observable is the speed of the track , which can be estimated with the _ linefit _ @xcite .
this algorithm is based on a simple ansatz in which the geometry of the cherenkov cone and the optical properties of the medium are ignored and the particle is assumed to travel with a velocity @xmath110 along a straight line through the detector . a pseudo-@xmath111 is constructed with the positions @xmath112 and times @xmath113 of the hlc pairs of all @xmath93 selected triplets :    @xmath114    hlc pairs which participate in multiple triplets are taken into account multiple times .
this @xmath115 can be minimized analytically with respect to the speed @xmath116 and vertex @xmath117 .
note that @xmath118 is arbitrarily normalized and can not be interpreted statistically in terms of goodness of fit .
the following analysis uses only the estimated speed @xmath119 .    in fig .
[ fig : velocity_reco ] the distributions of the reconstructed speeds are shown .
the reconstructed speeds are a reasonable estimate of the true speed , in particular for faint monopoles ( see @xmath120 ) . for brighter monopoles ( see @xmath121 ) ,
the reconstructed speeds slightly underestimate the true speed .
this reconstruction algorithm is simple , robust and fast , while still yielding a sufficient accuracy .
it also allows us to approximate the monopole direction by the direction of @xmath122 . the mean difference between the true and reconstructed direction varies between @xmath123 and @xmath124 depending on the monopole speed and
the mean free path @xmath45 .     and at the bottom for monopoles with @xmath121
are shown .
for comparison , reconstructed experimental data corresponding to a live - time of 8h are plotted .
the three dotted black lines show the true speeds and the speed of light .
all distributions are normalized to one.,title="fig : " ] +   and at the bottom for monopoles with @xmath121 are shown . for comparison ,
reconstructed experimental data corresponding to a live - time of 8h are plotted .
the three dotted black lines show the true speeds and the speed of light .
all distributions are normalized to one.,title="fig : " ]      for this first icecube analysis of slop data a robust approach based on @xmath72 as single final selection criterion and the determination of the expected background from experimental data was chosen .
figure [ fig : pdf_ntriplet_final ] shows the probability density distributions of @xmath72 for events with a reconstructed speed of at least @xmath125 ( top ) and with a reconstructed speed less than @xmath125 ( bottom ) . while the signal expectation extends to very high @xmath72 , the distributions of the experimental data decrease rapidly .
the final cuts on @xmath72 were optimized for maximum sensitivity based on the _ model rejection factor _ @xcite .
the optimization resulted in the following criteria / cuts : @xmath126 for a reconstructed speed @xmath127 and @xmath128 for @xmath129 .     for events with larger reconstructed speed ( top ) and for events with smaller reconstructed speed ( bottom ) . in black the distributions of one year
experimental data are shown .
the signal distributions are shown with decreasing @xmath45 in blue , red and green .
the final cuts on @xmath72 are shown by the dashed black line.,title="fig : " ] +   for events with larger reconstructed speed ( top ) and for events with smaller reconstructed speed ( bottom ) . in black the distributions of one year experimental data are shown .
the signal distributions are shown with decreasing @xmath45 in blue , red and green .
the final cuts on @xmath72 are shown by the dashed black line.,title="fig : " ]    these selection cuts were defined before unblinding the full experimental data . here , an iterative two step procedure was chosen .
first , 10% experimental data was unblinded with the selection determined by the aforementioned experimental 2 days data sample . after no signal or unexpected background
was observed the same procedure was applied to the full experimental data .
figure [ fig : velocity_ntriplet_final ] shows the resulting @xmath72-speed distribution for the full year of experimental data . after the final selection only one experimental event with @xmath130 and @xmath131 remains , but not well separated from background .
closer inspection revealed no evidence for an obvious track - like signature , in particular most triplets would not have survived tighter causality requirements .
as this observation is consistent with the expected number of about @xmath132 background events ( see below ) , we do not consider this result as positive detection .
-speed distribution of one year of experimental data .
the final cuts on @xmath72 are shown by the dashed red lines . the boundary between the two speed regions
is shown by the dashed black line . ]      with no observed monopole signal we have derived an upper limit on the flux of non - relativistic magnetic monopoles ( sec .
[ sec : flux_limits ] ) . for this
, the background model is fit to the @xmath72 distributions for both speed ranges ( fig .
[ fig : dist_ntriplet_final ] ) .
it is found that the background model ( sec . [
sec : backgroundmodel ] ) well describes the @xmath72 distribution over several orders of magnitude .
distributions for events with larger reconstructed speed ( top ) and for events with a smaller reconstructed speed ( bottom ) .
the black data points show the distributions of the full experimental data .
the expected signal is shown with decreasing @xmath45 in blue , red and green .
the fitted functions @xmath108 are shown in purple and the final selections on @xmath72 as dashed black lines.,title="fig : " ] +    distributions for events with larger reconstructed speed ( top ) and for events with a smaller reconstructed speed ( bottom ) .
the black data points show the distributions of the full experimental data .
the expected signal is shown with decreasing @xmath45 in blue , red and green .
the fitted functions @xmath108 are shown in purple and the final selections on @xmath72 as dashed black lines.,title="fig : " ]    based on the fit results , the expected number of background events with @xmath133 can be calculated by the integral    @xmath134    the total expected number of background events is defined by the sum of the expectation of both speed regions @xmath135 . by varying the fit - parameters within their fitted uncertainty a pseudo - experiment with different expected numbers of background events can be performed .
figure [ fig : pdf_background ] shows the resulting probability density distribution of the total expected number of background events .
the median expected number of background events is @xmath136 . here
, the statistical uncertainty is approximated by the difference between the median and the quantiles @xmath137 and @xmath138 .     and
@xmath138 are shown by the dash - dotted red line . ]
the search for magnetic monopoles presented in this section uses data taken during the season 2009 - 2010 when icecube was running in its 59-string configuration .
this analysis used the data taken with the standard icecube triggers .
the standard trigger that is used for highly energetic relativistic particles is a simple multiplicity trigger ( smt ) , which requires at least eight hlc hits within a sliding time window of 5 @xmath139 ( smt-8 ) .
other triggers are optimized for relativistic particles with lower energies .
data are recorded over at least the time interval over which the trigger condition of any of the triggers is fulfilled . for hlc hits ,
the full pmt waveforms are digitized and recorded @xcite .
not all triggered events were transmitted to the northern hemisphere by satellite .
events of various categories ( e.g. track - like , cascade - like , very bright events , etc . ) have been selected by various online filters at the south pole @xcite .
although the filters are optimized for relativistic particles , they may accept bright monopole events if a sufficient number of doms are hit .
this analysis uses the cascade and high - energy filters , which have the best acceptance for non - relativistic monopoles .
the total livetime of this data set is 311.25 days , with an average rate of selected events of 85.5hz .
the efficiency of this filter selection with respect to the multiplicity trigger is above 75% for monopoles of @xmath140 and @xmath141 .
slow monopoles with a catalysis cross section @xmath36 much larger than @xmath142 appear as very bright tracks .
simulations of the detector response to such tracks show that the multiplicity condition is fulfilled over most of the monopole crossing time , or that successive triggers occur close enough in time for the recording intervals to overlap .
so , a large fraction of a monopole s catalysis signature would be captured in a single event , if @xmath36 is sufficiently high . for @xmath143 @xmath144 monopoles
still yield multiple triggers , but the triggers occur less frequently , so that the signature is often split up into several sub - events
. the smaller the cross section , the more the monopole event splits up and the larger are the gaps between the sub - events .
eventually , the signal becomes indistinguishable from the background .
therefore , this analysis focuses on catalysis cross sections above @xmath145 @xmath144 . for monopoles with such high @xmath36
, the ic-59 analysis achieves a better sensitivity than the analysis using the slop trigger .
this is simply because the ic-59 array had a much larger detection volume than the deepcore array available to the previously described analysis .
future monopole searches will use data taken after 2012 , when the slop trigger was operating on the full ic-86 array .
these analyses will take advantage of both the large detection volume of the full ic-86 array and the high efficiency of the slop trigger .
the high - energy and cascade filters provide a data sample with about @xmath146 events .
the vast majority of these events are down - going atmospheric muons .
this background is reduced using a set of straight cuts in a first step .
these cuts are based on the time and location of the detected cherenkov photons .
contrary to the ic-86/deepcore analysis , whose cut parameters where defined using the time and location of dom launches or hlc pairs , this analysis uses a feature extraction algorithm on the pmt waveforms , which reconstructs the constituent pmt pulses caused by individual photo electrons . in a second step a multivariate analysis is adopted to reduce the background further .
the variables used for background reduction are :    1 .   the event duration @xmath147 defined as the time difference between the last and first pulse registred by a dom in an event .
2 .   the reconstructed speed v from the linefit .
3 .   the number of clusters ( @xmath148 ) , which is defined by the reconstructed pulses on all doms sorted into groups of pulses which occur close in space and time .
each such group is called a cluster and the total number of these clusters in an event is used as a cut variable .
bright signal tracks tend to have a higher number of clusters than atmospheric muon background events .
4 .   the total number of photoelectrons collected in the whole detector divided by the event duration , @xmath149 .
median of the distance between clusters along the reconstructed track
the center of gravity ( cog ) of the event , defined as the average spatial coordinates of all hits .
+ straight cuts are applied to variables 1 - 5 , chosen to substantially reduce background while keeping the signal efficiency reasonably high ( table [ table : cuts_effic_ic59 ] ) .
the cut on variable 6 removes events that only traverse a corner of the detector .
+ after applying those cuts , a multivariate analysis is performed on the remaining data to define a final selection criterion .
in addition to variables 3 , 4 and 5 this multivariate analysis considers the following variables : 7 .
mean distance of the hit doms to the center of gravity ( cog ) of the event divided by the event duration .
number of clusters divided by the event duration . 9 .
number of simple multiplicity triggers divided by the number of strings with hit doms      data are divided into two sets according to the monopole track brightness ( i.e. the catalysis cross section ) .
the cross section values for which we optimized the analysis and derive flux limits are @xmath150 and @xmath151 , which correspond to @xmath152 and @xmath121 respectively .
figure [ fig : time_durations_velocity_1 ] compares event duration @xmath147 and reconstructed speed v of experimental data to those of bright monopoles with simulated @xmath45 = 1 mm and speeds @xmath53 of @xmath31 and @xmath153 .
the signal efficiencies are 57.0% and 75.4% at the filter level for @xmath154 and @xmath155 , respectively .
a cut @xmath156s reduces the data by a factor @xmath157 while keeping 43.5% of the signal for @xmath158 and 45.3% for @xmath30 at the filter level .
note that the average duration of triggered events is shorter for slower speeds than for faster , because slower monopole events are more likely to be split into multiple sub - events .    , for experimental data and simulated bright monopoles with @xmath45 = 1 mm ( i.e. @xmath159 ) , before applying a cut .
the green histograms represents monopoles with @xmath160 , the blue histogram with @xmath154 .
the gray histogram represents the data .
the red dashed line marks the value of the chosen cut which is set at @xmath161 .
bottom : the same for the reconstructed speed v with a cut at v < @xmath162m / ns .
histograms are normalized to 1.,title="fig : " ] + , for experimental data and simulated bright monopoles with @xmath45 = 1 mm ( i.e. @xmath159 ) , before applying a cut .
the green histograms represents monopoles with @xmath160 , the blue histogram with @xmath154 .
the gray histogram represents the data .
the red dashed line marks the value of the chosen cut which is set at @xmath161 .
bottom : the same for the reconstructed speed v with a cut at v < @xmath162m / ns .
histograms are normalized to 1.,title="fig : " ]    relativistic single - muon tracks have a reconstructed speed v around @xmath163 .
events having passed the preceding cut on the time duration @xmath147 are enriched with coincident muons from uncorrelated air showers , resulting in a lower v ( see fig .
[ fig : time_durations_velocity_1 ] , bottom ) . for monopoles , the speed v is close to the simulated values .
cutting at v
< @xmath162m / ns ( corresponding to @xmath164 ) reduces the background by another order of magnitude .
further cuts on variables 3 - 6 reduce the background by another factor two . in total after this first set of cuts the data rate is reduced by a factor @xmath165 while the signal efficiencies only drop to 33.6% and 34.3% for @xmath166 and @xmath155 , respectively .
data reduction factors , rates and signal efficiencies before and after each applied cut are presented in table [ table : cuts_effic_ic59 ] .
for a ten times lower @xmath36 , the monopole tracks are dimmer and the signal efficiency drops dramatically . before applying any cut , the efficiencies at the filter level are 41% and 43% for @xmath154 and @xmath160 , respectively .
figure [ fig : time_durations_velocity_2 ] compares the same variables presented in fig .
[ fig : time_durations_velocity_1 ] .
signal and background are much less separated than for @xmath152 . moreover
, none of the events with @xmath30 has a duration that exceeds 800@xmath99s , which is much less than the 3ms necessary to cross the full array ; i.e. most events are split into one or more sub - events which have a shorter event duration in comparison to @xmath152 .
the properties of these sub - events are determined by hits from muons and from noise falling in the time window of the monopole passage .
thus the cuts applied on @xmath147 and v had to be slightly relaxed compared to @xmath152 : @xmath147 > 28 @xmath139 and v < @xmath167m / ns . after excluding events with a reconstructed center of gravity ( cog ) at outer strings ,
the data rate is reduced by a factor @xmath168 .
table [ table : cuts_effic_ic59 ] shows the data reduction factors , rates and the final signal efficiencies before and after each applied cut .
the signal efficiencies drop to 13.9% and 3.1% for @xmath166 and @xmath155 respectively .
l|lllllll & & + & before the cut & cut 1 & cut 2 & cut 3 & cut 4 & cut 5 & cut 6 + @xmath166 & 57% & 43.5% & 42.3% & 41.9% & 41.9% & 41.8% & 33.6% + @xmath155 & 75.4% & 45.3% & 41.1% & 41% & 41% & 39.8% & 34.3% + experiment : reduction factor &  & @xmath157 & @xmath169 & @xmath170 & @xmath171 & @xmath172 & @xmath173 + experiment : rate[s@xmath174 & 85.5 & @xmath175 & @xmath176 & @xmath177 & @xmath178 & @xmath179 & @xmath180 +   + & & + & before the cut & cut 1 & cut 2 & cut 6 + @xmath166 & 41% & 17.2% & 17% & 13.9% + @xmath155 & 43% & 3.5% & 3.23% & 3.1% + experiment : reduction factor &  & @xmath181 & @xmath182 & @xmath183 + experiment : rates[s@xmath174 & 85.5 & @xmath184 & @xmath185 & @xmath186 +     but with @xmath121 , i.e. 10 times lower .
the red dashed line marks the value of the chosen cut which is set at @xmath187 for the event duration and v < @xmath167m / ns for the speed .
histograms are normalized to 1.,title="fig : " ] +   but with @xmath121 , i.e. 10 times lower .
the red dashed line marks the value of the chosen cut which is set at @xmath187 for the event duration and v < @xmath167m / ns for the speed .
histograms are normalized to 1.,title="fig : " ]      to optimize the sensitivity for bright monopoles a multivariate analysis is used .
it classifies each event by a boosted decision trees ( bdt ) score in the range @xmath188 $ ] @xcite .
a bdt score of @xmath189 characterizes a background - like event whereas a bdt score of @xmath190 characterizes a signal - like event .    for the analysis ,
a sample of 10% of all experimental data ( burn sample ) was divided into two equally - sized sets .
bdts have been trained on each combination ( @xmath53,@xmath45 ) of the signal monte carlo and on 50% of the corresponding burn sample , using combinations of the variables described above .
the sensitivity was estimated from the other 50% of the burn sample by fitting an exponential function to the tail of the bdt score distribution . over a large range of the bdt scores ,
the fit describes the data rather well .
still , its extension into the signal region has no strict physics justification .
the final cut on the bdt scores for each combination of ( @xmath53,@xmath45 ) is obtained by using the model rejection factor ( mrf ) method @xcite . for the chosen high catalysis cross sections the limits for three ( @xmath53 , @xmath45 ) combinations are significantly better or comparable to those of the ic-86/deepcore analysis .
the fourth combination ( @xmath155 , @xmath191 cm ) is not competitive because the optimal cut results in 42 expected background events for the full data sample .
figure [ fig : bdt_limits ] shows the bdt scores for data and signal ( @xmath192 mm ) for one year data taking ( 311.25 days live time ) after the unblinding .
the optimized cut on the bdt scores leaves only one event which merely passes the cut .
no events pass the cuts for @xmath193 mm and @xmath193 cm .
the one surviving event was inspected visually .
it contains two nearly vertical high - energy muons which subsequently cross the whole detector and trigger two neighboring strings .
it has a time duration @xmath194 and a reconstructed speed of v = @xmath195m / ns .
mm , and speed @xmath160 .
the dot dashed line shows the optimized cut on the bdt score obtained from the model rejection factor method .
one event survived the bdt cut and is compatible with background . ]    the expected number of background events after unblinding is calculated using the exponential fit to the bdt score distributions of the one year experimental data .
the numbers of expected background events and of observed events , as well as the cut values on the bdt score for each parameter combination ( @xmath53 , @xmath45 ) are shown in table [ tab : background_events ] .
the higher number of expected background events compared to the actually observed number of events suggests that the exponential fit over - estimates the background rather than under - estimates it .
c|lll & + & @xmath196 & @xmath197 & @xmath198 + @xmath166 & 0.46 & @xmath199 & 0 + @xmath155 & 0.48 & @xmath200 & 1 +   + & + & @xmath196 & @xmath197 & @xmath198 + @xmath166 & 0.5 & @xmath201 & 0 + @xmath155 & not sensitive & not sensitive & not sensitive +
the calculation of upper flux limits takes into account the statistical and systematic uncertainties in the background and signal predictions .
because the number of expected background events is estimated from experimental data , only the statistical uncertainties of the fit parameters of the background model are relevant .    for signal the imperfect detector description
is taken into account .
for example in case of the ic-86/deepcore search , the random noise leads to an increase of @xmath72 for signal events .
furthermore the optical light detection efficiency is important .
this efficiency takes into account the cumulative effect of the light yield of nucleon decays , where a single electromagnetic cascade is simulated instead of several daughter particles , the light propagation through the ice and its detection by the doms .
these effects result in an uncertainty of the detection efficiency for magnetic monopoles which is used to derive the upper limits .
the impact of these uncertainties on the flux limits is estimated by simulating monopoles with simulation parameters changed within their estimated uncertainties .
the uncertainties of the superimposed background noise , the light yield of nucleon decays and the light propagation through ice are estimated by their differences in the detection efficiencies of signal simulations taking into account different approaches ( sec . [ sec : simulation_monopole ] ) . for the ic-86/deepcore analysis the superimposed noise
can be described by random and correlated noise hits from experimental data or noise simulated as a poisson process and atmospheric muons simulated using the software package corsika . since for the ic-59 analysis
no unbiased experimental data exists the background noise can be simulated by a noise generator that also takes into account correlated noise hits . for reasons of simplification
the proton decay is simulated as a single electromagnetic cascade with an isotropic direction which is valid as long as the mean free path is much smaller than the icecube spacing . due to kinematics in the proton decay ( eq .
[ eqn : promising_decay_channel ] ) two back - to - back electromagnetic cascades with an isotropic direction have to be simulated .
the uncertainties due to this simplification are estimated by the differences between both approaches .
for the light propagation through ice the two ice models described in @xcite and @xcite are used .
the uncertainty of the optical efficiency of doms can be estimated as @xmath202 .
signal simulations based on optical dom efficiencies varied by @xmath202 are compared with simulations based on the default settings .
the differences in the detection efficiencies are used as an estimate for the uncertainty .
we quantify each systematic effect @xmath203 by the ratio of the resulting detection efficiency @xmath204 relative to the detection efficiency with baseline assumptions @xmath205 :    @xmath206    the resulting changes are displayed in table [ tab : sys_eff_area ] .
cr|ll|ll|ll|ll & & & & & + & @xmath207 $ ] & & & & & & & & + & 3.0 & + 29%@xmath208 - 2% & + 23% & + 36% & -16% & + 36%@xmath208 - 23% & + 46%@xmath208 - 29% & + 20%@xmath208 - 11% & + 72% + & 1.0 & + 52% & + 6%@xmath208 - 11% & + 27% & + 1%@xmath208 - 11% & + 14%@xmath208 - 7% & + 16%@xmath208 - 20% & -17% & + 12% + & 0.3 & + 31% & + 5%@xmath208 - 8% & + 15% & + 1%@xmath208 - 8% & @xmath20910% & @xmath20911% & -16% & + 8%@xmath208 - 2% + & 0.1 & + 19%@xmath208 - 1% & + 3%@xmath208 - 4% & @xmath2097% & + 1%@xmath208 - 6% & + 5%@xmath208 - 11% & + 8%@xmath208 - 6% & -15% & + 1%@xmath208 - 6% + & 0.03 & + 17%@xmath208 - 2% & + 1%@xmath208 - 5% & + 9%@xmath208 - 4% & + 1%@xmath208 - 4% & + 10%@xmath208 - 5% & + 6%@xmath208 - 3% & -9% & -7% + & 0.01 & + 15%@xmath208 - 4% & -4% & + 11%@xmath208 - 2% & @xmath2092% & + 12%@xmath208 - 1% & + 5%@xmath208 - 0.3% & @xmath2095% & -10% + & 0.001 & + 15%@xmath208 - 4% & -4% & + 11%@xmath208 - 2% & @xmath2092% & + 12%@xmath208 - 1% & + 5%@xmath208 - 0.3% & @xmath209 5% & -10% + & 0.01 & + 2% &  &
-5% &  & + 9%@xmath208 - 6% &  & + 7% &  + & 0.001 & @xmath210 & @xmath211 & -3% & -1% & -5% & + 9%@xmath208 - 2% & + 2% & + 4% +    note , that these calculations are limited by computing resources and correspondingly by the statistics of simulated events .
the resulting statistical uncertainties of effective areas are typically a few percent as shown in table [ tab : stat_eff_area ] and included in the total error .
rr|ll & & + & @xmath212 $ ] & @xmath158 & @xmath30 + & 3.0 & @xmath213 & @xmath214 + & 1.0 & 8% & 5% + & 0.3 & 8% & 2% + & 0.1 & 8% & 2% + & 0.03 & 9% & 2% + & 0.01 & 10% & 2% + & 0.001 & 12% & 2% + & 0.01 & 0.2% &  + & 0.001 & 0.1% & 0.1% +    for the calculation of the final flux limits we perform high statistics computer experiments . in each
we randomize the effect of each systematics effect @xmath215 according to its specific uncertainty . for each parameter combination @xmath53 and @xmath36
this results in the effective probability density distribution for the relative change of the detection efficiency @xmath216 taking into account all uncertainties .
an example is shown in fig .
[ fig : signal_systematics ]
. multiplying these distributions by the detection efficiency with baseline assumptions @xmath217 one gets the probability density distributions for the detection efficiencies .
monopoles for the ic-86/deepcore analysis . ]
the flux limits on non - relativistic magnetic monopoles are calculated assuming an isotropic flux and the proton decay channel @xmath218 ( eq . [ eqn : promising_decay_channel ] ) with the catalysis cross section @xmath36 which depends on the speed @xmath53 ( eq . [ eqn : sigma_cat_theo ] ) .
using the quantity @xmath50 ( eq . [ eqn : monopole_photon_yield ] ) the flux limits can also be expressed without assuming a specific decay channel .
the flux limits are calculated based on a generalization of the approach by rolke et al .
@xcite , which takes into account the uncertainties of the signal detection efficiency and the expected number of background events .
therefore , a three - dimensional likelihood fit is performed with the following parameters : expected number of signal events @xmath99 , signal detection efficiency @xmath219 and expected number of background events @xmath220 .
the likelihood function is defined by    @xmath221    where the number of observed events @xmath222 follows a poisson distribution with the expectation value @xmath223 . the functions @xmath224 and @xmath225 represent the probability density distributions of the signal detection efficiency and the expected number of background events ( figs .
[ fig : pdf_background ] and [ fig : signal_systematics ] ) .
the flux limits for each monopole speed @xmath53 and catalysis cross section @xmath36 are calculated by    @xmath226    where @xmath227 is the upper limit at the 90% confidence level for @xmath228 event . the upper limit is obtained from the profile likelihood function defined in @xcite . in order to exclude unphysical values of @xmath227 the expected number of signal events @xmath99 is constrained to be greater than or equal to zero . as a consequence of this method
the upper limits at the 90% confidence level have a slight over - coverage of about 5% .
the other parameters are the size of the signal generation disc @xmath229 , the solid angle @xmath230 and the livetime @xmath231 .
figure [ fig : limits ] shows the resulting direct detection limits on the flux of non - relativistic magnetic monopoles in comparison to the current best experimental flux limits by the macro experiment @xcite .     and catalysis cross section @xmath36 of the ic-59 analysis and ic-86/deepcore analysis .
the dashed lines are limits published by the macro experiment @xcite . here , macro 1 is an analysis developed for monopoles catalyzing the proton decay .
macro 2 is the standard - macro - analysis , which is sensitive to monopoles ionizing the surrounding matter .
additionally , the icecube limits are shown as a function of @xmath50 which is proportional to the averaged cherenkov photon yield per nuceleon decay ( not valid for macro limits ) . ]    above @xmath232 corresponding to @xmath233 the previous flux limits are improved by more than one order of magnitude . moreover , for such large cross sections the monopole flux can be constrained up to a level @xmath234 which is three orders of magnitude below the parker bound .    assuming monopoles are the dominant part of dark matter , i.e. the relic mass density of monopoles is similar to the dark matter mass density , our most stringent flux limits constrain the monopole mass to be at least of the order of the planck mass @xmath235 @xcite .
this implies that monopoles with masses significantly smaller than the planck mass do not contribute dominantly to the dark matter mass density .
indirect searches for monopole induced proton decays set very strong bounds on monopoles with non - relativistic speeds , e.g. the limits from super - kamiokande @xcite assuming gravitational trapping of monopoles in the sun .
also a variety of bounds based on observations of neutron stars , white dwarfs and gas giants have been obtained .
these bounds range from @xmath236 and depend on the catalysis cross sections as well as on details of the assumed astrophysical scenarios @xcite .
although the direct icecube searches are not as stringent as indirect searches the former are not affected by astrophysical uncertainties .
thus the direct icecube limits can be considered as a robust upper bound on the monopole flux , if the rubakov - callan effect is realized in nature .
data taken from may 2011 until may 2012 with a dedicated slow - particle trigger and for the brightest monopoles data taken from may 2009 until may 2010 with standard - icecube - triggers were analysed .
the analysis , which is based on data of the slow - particle trigger , was developed by using simulated monopole events and experimental data to estimate background properties .
for this first analysis of such a signal in icecube a robust approach based on a single final selection criterion and the comparison between the number of expected background events and observed experimental events is chosen . using experimental data ,
the number of expected background events can be estimated to @xmath237 .
the ic-59 analysis based on standard - icecube - triggers is sensitive only for bright monopoles with @xmath238 .
the analysis used boosted decision trees ( bdt ) to discriminate between monopole signal and background .
the expected number of background events is derived from a fit of the bdt scores tails with an exponential function for each ( @xmath53,@xmath45 ) . the number of observed events after unblinding is @xmath239 for an expected background of @xmath240 .
this event contains multiple coincident muons , which renders it compatible with a background event .
the obtained flux limits for @xmath166 and @xmath241 m , @xmath242 m from the ic-59 analysis are better than the ones from the ic-86/deepcore analysis because of the bigger effective area . for @xmath155
the limits are comparable since the standard icecube triggers are less sensitive to the monopole signal in comparison to the dedicated slow - particle trigger .    in both analyses
no monopole signal has been observed .
thus , the limits on the flux of non - relativistic magnetic monopoles  catalysing the proton decay  are improved by about more than one order of magnitude in comparison to macro @xcite for most of the investigated parameter space and reach down to about three orders of magnitude below the parker limit .
since may 2012 the dedicated slow - particle trigger has been updated to the full icecube detector . from this upgrade , we expect an improvement in sensitivity by roughly an order of magnitude @xcite .
this gain is supplemented by improvements of the data selection which have been developed after completion of this analysis .
examples are the implementation of a kalman - filter based hlc hit selection , which improves the angular and speed reconstruction , and the implementation of an event selection based on a boosted decision tree @xcite .
we acknowledge the support from the following agencies : u.s . national science foundation - office of polar programs , u.s . national science foundation - physics division , university of wisconsin alumni research foundation , the grid laboratory of wisconsin ( glow ) grid infrastructure at the university of wisconsin - madison , the open science grid ( osg ) grid infrastructure ; u.s .
department of energy , and national energy research scientific computing center , the louisiana optical network initiative ( loni ) grid computing resources ; natural sciences and engineering research council of canada , westgrid and compute / calcul canada ; swedish research council , swedish polar research secretariat , swedish national infrastructure for computing ( snic ) , and knut and alice wallenberg foundation , sweden ; german ministry for education and research ( bmbf ) , deutsche forschungsgemeinschaft ( dfg ) , l.b . was funded by the dfg sonderforschungsbereich 676 , helmholtz alliance for astroparticle physics ( hap ) , research department of plasmas with complex interactions ( bochum ) , germany ; fund for scientific research ( fnrs - fwo ) , fwo odysseus programme , flanders institute to encourage scientific and technological research in industry ( iwt ) , belgian federal science policy office ( belspo ) ; university of oxford , united kingdom ; marsden fund , new zealand ; australian research council ; japan society for promotion of science ( jsps ) ; the swiss national science foundation ( snsf ) , switzerland ; national research foundation of korea ( nrf ) ; danish national research foundation , denmark ( dnrf )                                                                                p.  byron , h.  yang , and j.  zhub , boosted decision trees , a powerful event classifier , in _ statistical problems in particle physics , astrophysics and cosmology : proceedings of phystat05 _ vol .
40 , p. 139
, imperial college pr , 2006 .