Particle detection

Particle detection infers a particle's passage from a physical change in matter. A cloud chamber makes the spatial pattern of ionization visible; an ionization chamber instead collects electric charge from ionization. Choosing the readout determines what can be inferred from a single event.

The Richter-Ting independent J/psi discovery experiments show how different readouts can establish the same particle resonance: reconstruct daughter invariant masses after a beam hits a stationary target, or tune the initial center-of-mass energy of two colliding beams. Both require an ensemble of events and calibrated particle detectors.

A transition-edge sensor detects deposited energy using the strong change in electrical resistance near a superconducting transition temperature.

Inferring a quark flavor from the charges and identities of observed decay products. A mistaken assignment dilutes a flavour-dependent decay rate difference.

The reconstructed location at which a particle decays, inferred from the trajectories of its decay products. Comparing two decay vertices can determine a decay-time difference when the motion of the parent system is known.

Particle identification assigns a species hypothesis to a detector signal or reconstructed track. Curvature in a magnetic field supplies magnetic rigidity; additional responses, such as those of a Cherenkov detector or a particle shower, help distinguish an Electron from a hadron. The identity hypothesis enters the energy–momentum relation and therefore the reconstructed invariant mass. An apparatus can observe a peak before it can identify every decay channel, as the SPEAR psi energy-scan experiment illustrates.

The elapsed time for a particle to traverse a specified path. For constant speed $v$ and known path length $L$, $t=L/v$; varying speed instead requires integrating $dt=ds/v(s)$.

Inferring a particle trajectory from recorded detector positions or a continuous trail. The inferred trajectory depends on spatial resolution, detector geometry and the assumed particle-motion model; it is not a direct picture of the particle itself.

The mean energy lost by a charged particle per unit path length through a material, with the particle species, material and energy specified. It is distinct from the particle range obtained by following the particle until it stops.

For a fast ion traversing a fixed material thickness with only a small fractional energy loss, the dipole excitation probability per atomic encounter scales as $Z^2/v^2$ at fixed impact parameter. Summing transferred excitation energies and integrating the encounter measure $2\pi bdb$ gives the leading electronic stopping power scaling, with slower logarithmic dependence from the allowed excitation energies and impact-parameter range.

At fixed ion kinetic energy $E=M{v}^2/2$, a heavier isotope is slower, so the leading loss scales as $Z^{2}M/E$. At fixed velocity, isotopes of the same element instead have approximately the same electronic loss; plotting against energy per nucleon removes the leading isotope mass dependence. The charge should be the effective projectile charge if its bound electrons screen the nucleus. At lower speeds, excitation thresholds, slow-passage suppression, charge exchange and nuclear energy transfer invalidate a simple extrapolation of the high-speed law.

A cascade producing multiple secondary particles from an initial energetic interaction. Its observed multiplicity depends on which particles reach the detector and on detection efficiency.

A particle detector converts an interaction involving a particle into a readable signal. A cloud chamber is a particularly direct example: tiny liquid drops turn an invisible ion trail into a photographable line.

An anticoincidence detector vetoes a candidate event when another detector fires during the specified coincidence interval, helping reject penetrating charged particles.

A multiwire proportional chamber is a proportional counter: parallel anode wires between cathodes in a gas. Charged particle ionization frees Electrons; an electric field drives them toward wires, triggering an electron avalanche read by electronic amplifiers.^[ref]

Georges Charpak's electronic readout avoided film scanning and full spark chamber discharge.^[ref]

Opposite polarities identify the avalanche wire. Positive ion motion dominates integrated signal; detection precedes complete collection. Calibrated Electron travel times enable drift chambers.^[ref]

The original apparatus used $40\mu \mathrm{m}$ stainless-steel wires between meshes $7.5\mathrm{mm}$ away, flowing argon with organic vapor, negative mesh voltage and individual electronic amplifiers. A $5.9\mathrm{keV}$ X-ray source calibrated pulse height while voltage was varied.^[ref]

A coincidence circuit defined a narrow proton or pion beam. Scanning it across wires tested localization; the response edge matched the beam width, limiting the inference. Distant counts were attributed to accidental coincidences. Separate Electron scans measured pulse delay versus position, comparing gas mixtures and pulse-height-dependent electronics.^[ref]

Thresholds rejected small pulses, reducing detection efficiency. Wire-independence tests did not establish unlimited separation of overlapping tracks.^[ref]

The centre-to-centre spacing $p$ of adjacent wires. A wire address alone specifies a cell, rather than a continuous position.

Consider an original idealized track reconstruction model. Let equally spaced wires have positions $x_k=kp$ and assign each event to the nearest wire, reporting $\hat{x}=x_k$. Assume exactly one correct wire address, perfect detection efficiency, exact wire positions and a uniform distribution of true positions within each cell. These assumptions explicitly exclude missed pulses, false addresses and several unresolved tracks.

Writing $u=x-x_k$, the measurement error is $e=\hat{x}-x=-u$, with $u$ uniform on $[-p/2,p/2]$. Symmetry gives zero expected value. Its variance and standard deviation follow by integration:For chosen wire pitches $p=2\mathrm{mm}$ and $1\mathrm{mm}$ this gives $0.577\mathrm{mm}$ and $0.289\mathrm{mm}$ respectively. These are model standard deviations, not measured chamber spatial resolutions. The largest possible absolute error in this model is $p/2$, which is a different quantity.

A useful subtlety is that the formula depends on the illumination. For an unresolved pencil beam fixed at $x=x_k+0.4p$, every reconstructed position is $x_k$. The repeated estimates have zero variance, but every error is $-0.4p$. Small repeat-to-repeat spread can coexist with a biased position. The uniform-cell result therefore cannot be transferred to an arbitrary beam profile or detector selection.

The wire address supplies one transverse coordinate in this model. Combining several separated measurements supports track reconstruction; a single address does not specify a complete trajectory. If two charged particles enter the same cell during one accepted interval, an address-only record can be identical to that from one particle. No reconstruction of that record can recover the lost multiplicity without additional information.

Can two readout values locate a hit between wires? Only with a specified response model. Place two readout wires at $0$ and $p$, put $u=x/p\in [0,1]$, and stipulate calibrated, noiseless electric chargesThese are assumptions for a toy sensor response, not the natural response law of adjacent anode wires. Under these assumptions the centroidrecovers position exactly. The apparent gain in precision comes from adding calibrated analogue information. Taking the mean of two fired wire addresses without their electric charges would instead give the fixed midpoint $p/2$.

For a chosen noise model, replace the two electric charges by $q_L+n_L$ and $q_R+n_R$, with independent random variables $n_L,n_R$ of zero expected value and common standard deviation $s$. Assume $s/Q\ll 1$ and no threshold censoring. First-order expansion of the ratio givesFor stipulated $p=2\mathrm{mm}$ and $s/Q=0.02$, this approximation gives $0.0283\mathrm{mm}$ at the midpoint and $0.0400\mathrm{mm}$ at either endpoint. The endpoint values describe uncensored signed channel readouts, including small negative values after baseline subtraction. They do not assume that both channels would pass a positive threshold. A ratio with a fluctuating denominator also need not have exactly zero bias; this is a local small-noise approximation.

Three original counterexamples show why calibration and selection matter:

Likewise, two equal hits at $u=1/4$ and $u=3/4$ have exactly the same total two-channel response as one hit at $u=1/2$ with twice the electric charge. A centroid can hide multiplicity even when each charge measurement is exact.

Diffusion enabled the lecture's Orsay interpolation, without establishing this model's linear response.^[ref]

In 1962, Leon Max Lederman, Melvin Schwartz, Jack Steinberger and collaborators distinguished muon neutrinos from electron neutrinos using a neutrino beam.^[ref]

The Alternating Gradient Synchrotron delivered $15\mathrm{GeV}$ protons onto beryllium. Particle decay of pions and kaons supplied neutrinos; a $13.5\mathrm{m}$ iron shield preceded a ten-ton aluminium spark chamber.^[ref]

Production and detection are separate. In modern neutrino flavour notation,Shielding suppresses the daughter muon; a newly produced muon inside the spark chamber's aluminium marks an interaction.^[ref]

Beam-timed coincidence circuits and veto counters rejected entering tracks. Machine-off runs measured cosmic ray contamination; shield-thickness changes tested neutrons. Placing lead near the target shortened the meson decay path and reduced events. Cosmotron Electron calibration established particle shower recognition.^[ref]

The paper compared $34$ single-muon candidates, including $5$ estimated cosmic ray events, with six shower candidates unlike calibrated Electron showers. Kaon decays could supply genuine electron neutrinos. The deficit favored distinct neutrino flavours under lepton universality.^[ref]

The inference concerns production and detection channels. It does not measure neutrino mass or establish permanently fixed neutrino flavour. Later neutrino oscillation experiments address changes during propagation.^[ref]

How can the relative abundance of tracks and showers distinguish hypotheses? Consider an original teaching model, with chosen inputs throughout. Condition on a sample of $N$ selected events, each assigned exactly one of two labels: track or shower. Multitrack events and unclassified events lie outside this model. A shower label need not identify an Electron; background can imitate it. Conversely, a true Electron interaction need not receive a shower label.

Let $b$ be the background fraction among selected events. Among the selected signal events, let $f$ be the fraction whose associated charged lepton is an Electron, with the remainder associated with a muon. These are fractions after event selection, rather than fractions of the incident neutrino particle flux. Define the conditional probabilitiesHere $a$ and $r$ describe classification within the selected sample. They are not measured historical detection efficiencies. The law of total probability gives the overall shower probabilityThis formula is the bridge between an interaction hypothesis and a visible label. It keeps beam contamination, background and mistaken classification conceptually distinct: genuine electron neutrino interactions belong in $f$, while non-signal events belong in $b$.

Assume independent random variables for the event labels, all with this same $q$. If $K$ counts shower labels, every specified pattern with $k$ showers has probability $q^k(1-q{)}^{N-k}$. Counting those patterns with a binomial coefficient yields the binomial distributionThe track count is $N-K$, so its expected value is $N(1-q)$. At fixed $N$, the two counts cannot fluctuate independently. The variance above measures sampling fluctuations under the stipulated model, not uncertainty in $a$, $r$, $b$ or $q_b$.

Choose $N=40$, $a=0.90$, $r=0.05$, $b=0.10$ and $q_b=0.40$. Compare a signal with only muon-associated interactions ($f=0$) against one with equal selected Electron- and muon-associated interactions ($f=1/2$). These give $q=0.085$ and $q=0.4675$, respectively. Their expected shower counts are $3.4$ and $18.7$; their expected track counts are $36.6$ and $21.3$. A background-only model with the specified $q_b$ gives $16$ showers and $24$ tracks on average. These are expectations, so they need not be integers.

For a hypothetical observation $K=4$, the likelihood ratio of the two signal models isThe binomial coefficient cancels because the same hypothetical sample is used in both likelihoods. This ratio is not the probability that either hypothesis is true. A separate original check gives $\mathrm{Pr}(K\le 5)\simeq 0.880$, $4.57\times 10^{-6}$ and $1.44\times 10^{-4}$ for the muon-associated, equal-mixture and specified background-only models, respectively.

Why are the controls essential? If the background-only model may choose any $q_b$, taking $q_b=K/N$ can fit the observed shower fraction. Counts alone then fail to distinguish it from a signal model. Independent background measurements and calibrated particle identification constrain that freedom. Changing the decay path additionally tests where the candidate-producing particles originate. The power of the experiment comes from combining apparatus interventions with event classification; the small binomial distribution model isolates only the latter part of that reasoning.

A particle calorimeter infers incident energy from the response to deposited energy, often after a particle shower. Segmentation supplies directional information relative to the collision position. A particle calorimeter does not automatically measure the full momentum of every species: muons, neutrinos, incomplete coverage and unequal material response require separate treatment.

A particle calorimeter designed to measure Electron and photon energies through particle showers. Track reconstruction helps distinguish charged Electrons from neutral photons.

A drift chamber uses the drift times of ionization Electrons in a gas to infer positions along charged particle tracks. In a known magnetic field, track reconstruction can estimate momentum and the sign of electric charge.

How can a detector establish that an unseen particle changes identity? The Super-Kamiokande atmospheric-neutrino experiment compares different travel paths; the Sudbury Neutrino Observatory compares different interaction channels. Both infer neutrino flavor from reaction products and calibrated responses. Their complementarity connects particle detection to neutrino oscillation, neutrino mass and neutrino astronomy.

Takaaki Kajita and Arthur B. McDonald shared the 2015 physics Nobel Prize. These were collaboration achievements: Yoichiro Suzuki's Super-Kamiokande work included the solar neutrino analysis, and George Ewan helped establish and lead the Sudbury Neutrino Observatory.^[ref][ref]

The older Kamiokande account explains water Cherenkov neutrino reconstruction; Raymond Davis Jr.'s Homestake experiment explains why a deficit in an electron neutrino channel alone leaves a source-versus-propagation ambiguity. Complementary channels can break that ambiguity.

Super-Kamiokande used $535$ days ($33.0$ kiloton-years) and $22.5$ fiducial kilotons of water. Inner photomultiplier times and charges reconstructed interactions; the outer anticoincidence detector vetoed entering particles and marked exits. Contained single-ring events used likelihood classification: Electron particle showers blur rings; muons give sharper edges. The $1.33\mathrm{GeV}$ visible-energy sample division is not a universal detection threshold. Simulation estimated classification efficiency above about $98\%$, distinct from total detection efficiency; stopping muons and decay Electrons checked classification.^[ref]

Downward neutrinos travel about $15\mathrm{km}$; upward ones about $13000\mathrm{km}$. The multi-GeV muon-like up/down asymmetry was $-0.296\pm 0.048\pm 0.010$; Electron-like events lacked a comparable effect. Calibration constrained directional bias. This was disappearance evidence; tau neutrino appearance was untagged and sterile conversion remained possible.^[ref]

Writing $U$ and $D$ for the selected upward and downward populations, the asymmetry is $A=(U-D)/(U+D)$: a negative value means fewer upward events. An absolute source-normalization error alone rescales both populations and leaves this ratio unchanged. A direction-dependent response error need not cancel.

For a deliberately simplified neutrino oscillation calculation, take two neutrino flavors related to two mass eigenstates by a real rotation matrix with angle $\theta$. A quantum superposition acquires different phases during propagation; projecting back onto the initial neutrino flavor gives the probabilityHere $L$ is path length and $E$ is neutrino energy; the numerical coefficient incorporates unit conversion.

The paper states the corresponding conversion formula.^[ref]

The path dependence has a geometrical origin. In an original spherical-Earth model, let $R$ be the radius and $h$ a fixed production height, neglecting detector depth. Let $\Theta$ be the zenith distance of the direction in the sky from which the neutrino arrives, so $c=\cos \Theta =+1$ means downward propagation. The backward ray from detector to production obeys $R^2+L^2+2RLc=(R+h{)}^2$. Its positive solution isThus $L(+1)=h$ and $L(-1)=2R+h$. The calculation explains why opposite directions test very different baselines, rather than attributing every observed direction to the same travel distance. Actual production heights and the angle between the outgoing charged lepton and incoming neutrino smear that relationship.

For an original chosen example, set $\theta =\pi /4$ and $\Delta m^2=2.5\times 10^{-3}{\mathrm{eV}}^2$. At $E=1\mathrm{GeV}$, $L=20\mathrm{km}$ gives $P_{\mu \mu }\simeq 0.9960$, while $L=500\mathrm{km}$ gives $P_{\mu \mu }\simeq 0.000279$. The first minimum is at $L/E\simeq 494.7km/GeV$. These are stipulated inputs, not a fit to the 1998 data. A broad distribution of true $L/E$, together with energy and directional uncertainty, smooths the sharp minima.

A direct amplitude check makes the mechanism explicit. DefineThen $S_{11}={\cos }^2\theta +{\sin }^2\theta e^{-2i\phi }$. Taking its squared absolute value yields the displayed survival formula because $1-\cos (2\phi )=2{\sin }^2\phi$. Unitarity gives $|S_{11}{|}^2+|S_{21}{|}^2=1$: disappearance from one channel becomes conversion to the other within this closed model.

This illustration does not describe all solar neutrino propagation. Coherent forward scattering in matter changes the effective propagation operator, and three active neutrino flavors require a larger mixing description. Even in the simplified model, neutrino oscillation measures a squared-mass difference; it does not determine the absolute neutrino mass.

The SNO Collaboration used $1000$ tonnes of heavy water in a $12\mathrm{m}$ acrylic sphere, surrounded by ordinary water shielding and $9456$ photomultipliers. Weak charged current (CC) breakup of a deuteron selects electron neutrinos; weak neutral current (NC) breakup measures all active neutrino flavors; Electron elastic scattering (ES) has reduced non-electron sensitivity.^[ref]

In 2001, $240.95$ live days, reconstructed radius below $5.50\mathrm{m}$ and effective Electron kinetic energy above $6.75\mathrm{MeV}$ gave a CC neutrino flux of $1.75\times 10^6{\mathrm{cm}}^{-2}{\mathrm{s}}^{-1}$ for the assumed boron-8 spectrum. Source calibrations measured $1.4\%$ instrumental-cut signal loss. SNO's own CC–ES difference was $1.6\sigma$; comparison with Super-Kamiokande's precise ES flux gave $3.3\sigma$ evidence for non-electron active neutrinos. This was not an NC measurement.^[ref]

In 2002 the pure-heavy water analysis used $306.4$ live days and a $5\mathrm{MeV}$ reconstructed kinetic-energy threshold. NC neutrons captured on deuterons emitted $6.25\mathrm{MeV}$ gamma rays, producing secondary Electron light. The NC reaction threshold was $2.2\mathrm{MeV}$, a different quantity. Californium-252 calibrated $29.9\%$ capture efficiency, but only $14.4\%$ detection efficiency after fiducial and energy cuts. Assays constrained backgrounds from radioactive decay and photodisintegration.^[ref]

Energy, radius and direction relative to the Sun statistically separated channels. Assuming the standard boron-8 spectral shape, the NC total active neutrino flux was $5.09_{-0.43}^{+0.44}{(stat.)}_{-0.43}^{+0.46}(syst.)\times 10^6{\mathrm{cm}}^{-2}{\mathrm{s}}^{-1}$, consistent with stellar structure predictions. The non-electron active component differed from zero by $5.3\sigma$. This established active flavour conversion, without individually tagging muon neutrinos and tau neutrinos. It did not exclude every sterile admixture.^[ref]

Why compare inferred neutrino fluxes rather than raw event counts? For channel $j$, exposure time $T$, target population $N_j$, cross section ${\sigma }_{j\alpha }$, and detection efficiency ${ϵ}_{j\alpha }$ including reconstruction and selection, a forward model isHere $B_j$ is expected background; $N_j$ on the left counts events, while $N_j^{\mathrm{target}}$ counts targets. Distinct cross sections, target numbers, thresholds and detection efficiencies prevent a direct CC/NC count ratio from being an electron neutrino fraction. Channel populations themselves can be correlated estimates from overlapping event patterns.

An original reduced linear algebra model isolates the inference. Assume specified spectral shapes and divide each background-corrected channel rate by its own calibrated electron neutrino response. Put ${\varphi }_x={\varphi }_{\mu }+{\varphi }_{\tau }$ and let $r$ be the spectrum- and response-weighted non-electron/electron ES sensitivity ratio. ThenThe equal $\mu$ and $\tau$ responses determine their sum, not their separate neutrino fluxes. Changing the assumed spectra generally changes the integrated coefficients; with arbitrary unknown spectra this small matrix is insufficient.

Choose illustrative inputs, not historical measurements: $r=0.15$, ${\varphi }_e=2$ and ${\varphi }_x=3$, with every neutrino flux in the same arbitrary unit. The forward prediction is $(y_{\mathrm{CC}},y_{\mathrm{NC}},y_{\mathrm{ES}})=(2,5,2.45)$. CC and NC independently invert to ${\varphi }_e=y_{\mathrm{CC}}=2$, ${\varphi }_x=y_{\mathrm{NC}}-y_{\mathrm{CC}}=3$; substituting into ES returns $2+0.15(3)=2.45$. Alternatively, CC and ES give ${\varphi }_x=(y_{\mathrm{ES}}-y_{\mathrm{CC}})/r=3$. Thus the third channel tests the response model and the inferred mixture.

CC alone admits any nonnegative ${\varphi }_x$. CC and ES lose the second independent constraint when $r=0$; for small positive $r$, an ES-minus-CC error is amplified by $1/r$. With chosen independent errors ${\sigma }_{\mathrm{ES}}={\sigma }_{\mathrm{CC}}=0.1$, their propagated standard deviation is $\sqrt{0.1^2+0.1^2}/0.15\simeq 0.943$. If CC and NC instead have those independent errors, the difference has standard deviation $\sqrt{0.1^2+0.1^2}\simeq 0.141$. Actual analyses must retain covariances, including those introduced by statistical channel separation. Complementarity supplies information; calibration establishes what information each channel carries.

IceCube turns Antarctic ice into a Cherenkov detector: $5160$ digital optical modules on $86$ strings span depths of $1450$–$2450\mathrm{m}$. Hot-water drilling allowed deployment before refreezing. Neutrino interactions produce charged particles; their Cherenkov radiation, rather than the incoming neutrino, supplies the optical readout.^[ref]

Neutrino astronomy asks where the incoming neutrino originated. A Cherenkov detector asks how its reaction products make measurable light. Connecting these questions requires track reconstruction, detection efficiency and background tests. An astronomical messenger and the light carrying information inside the apparatus are different stages of one inference.

A glass-enclosed photomultiplier with local digitization: pulse times and charges enter event reconstruction. Reciprocal electrical pulses synchronize module clocks to the surface time base.^[ref]

Weak charged current interactions of muon neutrinos usually give extended muon tracks. Compact particle showers occur in other channels, including weak neutral current interactions of every neutrino flavor. Tracks generally constrain direction better; escaping particles complicate energy inference.^[ref]

In the 2013 calibration study, six horizontal light-emitting diodes per module on string 63 flashed sequentially; other strings recorded charges and photon arrival times. A global fit inferred depth-dependent scattering of light and absorption coefficients. Dust-logger measurements checked layer structure; independent muon events tested simulated time patterns. Photomultiplier response and refrozen ice near modules remained limitations. Thus a late hit need not imply a late interaction: scattering of light changes the path.^[ref]

An original limiting calculation isolates this issue. Suppose a compact source emits at time $t_0$ into homogeneous, nondispersive ice with stipulated refractive index $n=1.32$. At distance $r$, the unscattered light arrives atwhere $c$ is the vacuum speed of light. For chosen $r=100\mathrm{m}$, the delay is $440.3\mathrm{ns}$. A scattered path $30\mathrm{m}$ longer adds $132.1\mathrm{ns}$ even with perfect clocks. These chosen values are not calibration measurements. Real ice requires wavelength-dependent propagation, module response and a spatially varying light-transport model; replacing that response by one arrival time discards information.

The 2013 selection rejected early boundary light. Tagged muons checked leakage; air shower companion muons also vetoed downward atmospheric neutrinos.^[ref]

A simple original conditional probability model explains why a veto changes the selected mixture. Let $A$ mean that an atmospheric neutrino interaction would pass all other cuts, $M$ that an accompanying muon reaches the veto, and $V$ that the veto fires. Set $a=P(M\mid A)$ and $ϵ=P(V\mid M,A)$, and neglect unrelated veto hits. The law of total probability givesFor chosen $a=0.8$ and $ϵ=0.95$, the passing fraction is $0.24$. Neither this value nor either input is an IceCube measurement. With equal incident candidate populations, this mechanism alone reduces the atmospheric contribution relative to a hypothetical population having $a=0$. It does not make every surviving event astronomical: $a=0$ or $ϵ=0$ gives no suppression, while $a=ϵ=1$ removes the entire chosen atmospheric population. If unrelated light vetoes a surviving event independently with probability $u$, its passing fraction becomes $(1-u)(1-aϵ)$; independence is an additional assumption.

In $662$ live days during 2010–2012, $28$ events exceeded the expected $10.6_{-3.6}^{+5.0}$ atmospheric events. Energy, direction and topology supported an extraterrestrial component. Combined searches gave $4.1\sigma$ with additional benchmark charm background; uncertain charm production affected significance. No source clustering was established.^[ref]

A diffuse excess and an identified source answer different questions. The first compares populations against a background hypothesis; the second needs an association in direction, perhaps also in time.

Escaping muons or neutrinos can make deposited energy smaller than incoming neutrino energy.^[ref]

How much does a small directional uncertainty help? Consider an original idealized model: a source direction fixed in advance, uniform background directions over the full sphere, and circular acceptance radius $\psi$. The accepted solid angle and single-event background probability areThe derivation integrates the area element on a unit sphere. For small $\psi$ in radians, $q\simeq {\psi }^2/4$: halving the radius quarters the accidental acceptance, to leading order. At $psi=0$ acceptance vanishes; at $\psi =\pi$ it covers the full sphere. The solid angle has units of steradians; $q$ is dimensionless.

For $0<\psi \le \pi$, take a deliberately crude point spread function: stipulate that every signal direction lies uniformly inside that cap. The signal and background probability density functions per steradian are $S=1/\Omega$ inside the cap and zero outside, and $B=1/(4\pi )$ everywhere. Both integrate to one. For an event inside the cap their likelihood ratio isThis ratio compares two specified directional models; it is not the probability that the event came from the source. Such a conditional probability would additionally require prior source odds and the other selection information. For fixed source fraction $f$ among selected events, a normalized mixture distribution has density $(1-f)B+fS$. For independent events the likelihood is the product of those densities; its ratio to background alone isOutside-cap events contribute factors $1-f$ rather than being silently discarded. At $f=0$, $\Lambda =1$; for a full-sky cap, $S=B$ and every $f$ gives $\Lambda =1$. A directional model that cannot distinguish source from background carries no directional evidence.

As an original numerical example, radii $1^{\circ }$ and $10^{\circ }$ give $q\simeq 7.615\times 10^{-5}$ and $7.596\times 10^{-3}$, respectively: nearly a hundredfold difference. For a fixed number $N$ of independent background events, the binomial distribution givesWith chosen $N=100$, those radii give approximately $0.00759$ and $0.5335$. A small single-event acceptance can become a common coincidence after many opportunities. This probability concerns a fixed cap, rather than a scan for the most favorable cap anywhere. For $m$ independent searches with common local tail probability $p$, the corresponding global tail probability is $1-(1-p{)}^m$; overlapping sky searches generally violate that independence. Actual source tests must include directional response, nonuniform exposure, energy information and the look-elsewhere effect in their null calibration.

The 2018 collaboration paper associated the 22 September 2017 neutrino alert with the flaring blazar TXS 0506+056 through directional and multiwavelength observations. Its approximately $3\sigma$ evidence included trial corrections and assumed relationships between neutrino and gamma ray emission. A single association did not establish a production mechanism or explain the entire diffuse neutrino flux.^[ref]

A particle detector in which an ionizing event triggers a gas discharge and a counted electrical pulse.

A detector trigger decides which candidate events cause a readout or recording. It can concentrate a limited recording capacity on rare interactions, but the accepted events then form a selected sample. A coincidence circuit uses a timing relation between several detectors; the Blackett-Occhialini triggered cloud-chamber experiment made that relation operate an entire camera-and-expansion sequence.

A coincidence circuit responds when detector signals occur within a specified time window. Two counters placed on opposite sides of a chamber can preferentially select a charged particle that crosses both. Simultaneity within the electronics' window alone does not establish that two signals have the same physical cause: unrelated detections produce accidental coincidences.

Suppose two independent background streams follow Poisson processes with rates $r_1$ and $r_2$. Define a coincidence by $|t_2-t_1|<w/2$: $w$ is the full accepted time-difference width. Around each first-counter signal, the expected number of second-counter signals is $r_{2}w$. The expected pair rate is therefore $r_1{r}_{2}w$. A circuit producing at most one output per first signal instead givesThe distinction matters when multiple detections occur in one window. If a different convention accepts $|t_2-t_1|<\tau$, its full width is $w=2\tau$ and the dilute expression is $2r_1{r}_2\tau$.

For a chosen numerical model, take background rates $r_1=100{\mathrm{s}}^{-1}$, $r_2=80{\mathrm{s}}^{-1}$ and $w=20\mu \mathrm{s}$. The dilute pair rate is $0.16{\mathrm{s}}^{-1}$. These are stipulated teaching inputs, not measurements of the historical chamber.

Making $w$ smaller also rejects some genuine events. Suppose true pairs arrive at rate $S_0$ before the timing cut, and their time difference has a centered normal distribution with standard deviation ${\sigma }_t$. Its accepted fraction, using the error function, isFor chosen $S_0=2.0{\mathrm{s}}^{-1}$ and ${\sigma }_t=5\mu \mathrm{s}$, the $20\mu \mathrm{s}$ window retains about $95.45\%$ of true pairs. In the dilute background model the fraction of accepted events that are true is $S/(S+R_{\mathrm{acc}})\simeq 92.27\%$. A narrow window trades detection efficiency against background; it does not make either quantity equal to one. Delaying one counter's signal far beyond genuine event correlations gives a proposed background check under stationarity, without assuming that the zero-delay excess itself is accidental.

A trigger changes which particle populations are common in the recorded sample. Let $N$ be the number of particles in an event, and stipulate an independent Bernoulli distribution for each particle's ability to produce an accepted telescope signal, with success probability $ϵ$. If any successful particle suffices, the event acceptance isThis is a deliberately simplified model. Its $ϵ$ includes the joint geometrical acceptance and response of the whole telescope; it is not an assumed universal single-counter efficiency. Actual shower directions, spatial correlations and shared electronics can violate the independence assumption.

The conditional probability formula givesFor chosen underlying fractions $\mathrm{Pr}(N=1)=0.90$, $\mathrm{Pr}(N=5)=0.10$ and $ϵ=0.15$, the two acceptances are $0.15$ and $0.5563$. Five-particle events consequently form about $29.18\%$ of recorded events. A population that is only $10\%$ of incoming events can occupy nearly $30\%$ of the photographs. This is selection bias, not evidence that the source population changed.

The interval during which a detector or readout cannot register another event while recovering from a response. A detector trigger may be fast while the camera, chamber or data acquisition is still busy.

In a simple nonextending model, every accepted event starts a fixed busy interval ${\tau }_d$; arrivals during it are lost and do not extend it. With Poisson process input rate $R$, the average cycle lasts ${\tau }_d+1/R$, giving observed rateThis is a specified mathematical model, not a universal correction. An instrument whose busy interval is extended by further arrivals behaves differently. Correcting the total rate also does not remove the multiplicity bias from a detector trigger.

A cloud chamber contains supersaturated vapor. Ionization along a charged particle's path supplies sites for nucleation; the resulting liquid drops scatter light and reveal the trail. The observed objects are drops, not magnified particles. Its beauty is that a spatially distributed microscopic interaction becomes a visible geometrical record.

Two ways of making supersaturated vapor must be distinguished. An expansion cloud chamber rapidly increases the volume of moist gas, producing a temporary sensitive interval through cooling by adiabatic expansion. A diffusion cloud chamber maintains a cold base and a warmer vapor supply, producing a continuously replenished sensitive layer through diffusion and a temperature gradient.

The central mechanism is a competition between two falling pressures: expansion reduces the actual water partial pressure, but cooling can reduce the equilibrium vapor pressure still faster. The ratio can therefore rise above one even though the gas becomes more dilute. The cloud-chamber expansion model below makes that comparison explicit.

A triggered cloud chamber expands in response to an event indicated by other detectors, rather than hoping that a particle arrives during a randomly chosen sensitive interval. The ionization trail is produced before expansion. The apparatus must respond while the ions remain sufficiently localized to make a useful photograph.

Patrick Maynard Stuart Blackett and Giuseppe Occhialini began their cosmic-ray collaboration in 1931. They used the coincidence methods of Walther Bothe and Bruno Rossi to make a particle initiate its own recording. Blackett later reported tracks on about $80\%$ of their counter-controlled photographs, compared with roughly $2$–$5\%$ for random expansions in small chambers. Those figures describe the reported apparatus and conditions, not a universal detection efficiency.^[ref]

In their 1933 account, Geiger counters above and below the chamber fed a coincidence circuit. Its pulse operated a thyratron, magnet and spring-release valve, letting a light piston expand the gas. A mercury-vapor lamp then illuminated the drops. A water-cooled solenoid supplied the magnetic field; two oblique cameras allowed track reconstruction. Moving the lower counter changed which trajectories could be recorded.^[ref]

Two timing problems must be separated. The coincidence circuit's window distinguishes nearly simultaneous counter signals. The mechanical expansion latency determines how much the existing ionization trail spreads before droplets form. For a chosen one-coordinate diffusion model with initial variance ${\sigma }_0^2$, diffusion coefficient $D$ and expansion delay $t_e$,Taking ${\sigma }_0=0$, $D=2.0\times 10^{-6}{\mathrm{m}}^2{\mathrm{s}}^{-1}$ and $t_e=20\mathrm{ms}$ gives ${\sigma }_{\mathrm{track}}\simeq 0.283\mathrm{mm}$. This is an original illustrative calculation, not the historical track width. It shows why a fast electronic pulse does not by itself guarantee a sharp photograph: the subsequent mechanical response also matters.

The photographs revealed groups containing positive and negative electron-like tracks, complementing the Anderson positron track experiment. Their interpretation depended on track curvature, direction, ionization and particle range, not simply on the presence of a coincidence.^[ref] Blackett explicitly noted that the counter control favored larger particle showers.^[ref] The multiplicity bias from a detector trigger gives a small model of that effect: a recording method can improve the study of rare interactions while making its photographs unrepresentative of the incident population.

Charles Thomson Rees Wilson's 1927 lecture recalls Electron cloud trails produced with X-rays in 1911, then alpha particle trails from a radium-tipped source. The improved 1912 chamber had a wide, shallow space with a suddenly falling floor to limit stirring. A roof–floor clearing electric field removed unwanted ions. Sharp photography required expansion, then the ionizing particle, then illumination.^[ref]

The piston mechanism connects the space beneath the moving floor to an evacuated vessel: opening the valve lets the floor drop. Water in a trough supplies the vapor. The Smithsonian's surviving demonstration object is explicitly a 1961 replica of the 1912 apparatus.^[ref]

After dust removal, X-rays increased droplet formation above a lower expansion threshold; larger expansions also produced dense background cloud.^[ref]

The control logic separates a localized ionization trail from fog. Delays permit ion displacement and drop growth; stirring deforms the trail. A photograph alone supplies neither an absolute event time nor a calibrated particle energy.

Consider an illustrative water–air mixture, initially saturated at $T_0=293.15\mathrm{K}$. Let $\lambda =V/V_0$ be the expansion ratio and approximate the carrier gas by a uniform ideal gas with constant heat capacity ratio $\gamma =1.40$. Before significant condensation, assume fixed numbers of carrier and vapor molecules, negligible incoming heat, and the reversible adiabatic expansion relation. The first law of thermodynamics, $n{C}_{V}dT=-pdV$, with the ideal gas law, givesThe vapor follows the same temperature, so its partial pressure scales as $T/V$. A rapid real piston motion need not be reversible; the formula is an idealized bulk estimate after neglecting flow and losses.

Use a physically motivated Integrated Clausius-Clapeyron relation for the equilibrium water vapor pressure, with stipulated rounded constants $p_0=2339\mathrm{Pa}$ at $T_0$, molar latent heat $L_m=45000\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}$ and molar gas constant $R=8.314\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$:These constants define the example, without fitting any historical measurements. Throughout the plotted interval $1\le \lambda \le 1.17$, the temperature stays above the freezing point of water. InitiallyThus the decrease in vapor pressure initially beats dilution. At $\lambda =1.15$, the model gives $T=277.21\mathrm{K}$, $p_v=1923.33\mathrm{Pa}$, $p_{\mathrm{sat}}=809.07\mathrm{Pa}$ and $S=2.3772$.

Why does $S>1$ not instantly make every molecule a visible drop? In the capillarity approximation of classical nucleation theory, a neutral spherical drop of radius $r$ has a positive surface energy cost and a negative bulk Gibbs free energy contribution. For molecular liquid volume $v_l$, surface tension $\sigma$ and Boltzmann constant $k_B$,Differentiating $\Delta G$ shows that the critical nucleus at $r_{\ast }$ is a maximum: smaller neutral drops tend to shrink, whereas larger ones can grow. With stipulated $\sigma =0.075\mathrm{J}{\mathrm{m}}^{-2}$, $v_l=2.99\times 10^{-29}{\mathrm{m}}^3$ and $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$, the $\lambda =1.15$ example gives $r_{\ast }=1.353\mathrm{nm}$ and $\Delta G_{\ast }/(k_{B}T)=150.3$. A substantial neutral-drop barrier can coexist with supersaturated vapor.

This is deliberately a neutral-drop comparison. An ion changes the cluster free energy through electrical interactions and molecular hydration; the formula does not calculate the ion-assisted nucleation threshold. At nanometre radii, using bulk surface tension is itself approximate. There is no predicted droplet count, growth time or detection efficiency here. Once many drops grow, latent heat release, vapor depletion and heat transfer from the chamber walls invalidate the dry adiabatic expansion trajectory. The example therefore explains the direction of the effect and the need for nucleation, without extrapolating its curves to Wilson's cloud limits.

A diffusion cloud chamber continuously supplies vapor from a warmer reservoir toward a colder base. A temperature gradient permits supersaturated vapor in a layer near the base, where new ionization trails can continually form drops. There is no repeated bulk piston expansion cycle. The two instruments share ion-assisted condensation, but their sensitive regions and timing differ.

CERN's modern teaching apparatus uses alcohol-soaked felt, a sealed transparent container, a dark metal base cooled with dry solid carbon dioxide, and a side light. The view is a thin sensitive region near the cold base, rather than a chamber-wide expansion photograph.^[ref][ref]

A bubble chamber photographs vapor bubbles along charged particle paths in a transparent superheated liquid. The liquid serves as target and recording medium.^[ref]

Compare the two phase changes: a cloud chamber forms liquid drops inside supersaturated vapor, whereas a bubble chamber forms vapor cavities inside a superheated liquid. Both use metastability to amplify a microscopic disturbance. The photons reaching the camera reveal the drops or bubbles. The photographed trail's width is not the particle's width: it includes bubble growth, illumination and spatial resolution.

Donald Arthur Glaser developed the bubble chamber in 1952 during experimental particle research at Michigan.^[ref]

His 1960 lecture, Elementary particles and bubble chambers, describes evacuated capillary bulbs containing pure diethyl ether. Cooling one bulb reduced pressure; the other remained hot. A cobalt-60 source initiated boiling. Small glass chambers subsequently used hot-oil baths, hand-cranked pistons and high-speed photography. Automatic cycling with counter-controlled flashes produced finer tracks. Local ionization deposits energy that can initiate bubble nucleation.^[ref]

Later metal-chamber photographs compared operation with and without a nearby radioactive source: wall boiling occurred without interior tracks; irradiation produced tracks before disruptive flow ended sensitivity.^[ref]

These observations separate two questions: can radiation start boiling, and can a detector preserve the geometry of a passing charged particle? Whole-vessel boiling alone answers only the first. Wall background, finite sensitivity and growth before photography constrain the second.

Reduce pressure before beam arrival; photograph bubbles from several views, then recompress.^[ref]

Unlike the Blackett-Occhialini triggered cloud-chamber experiment, Glaser used counter-controlled flashes after expansion. A detector trigger can select recording in sensitive liquid; subsequent expansion cannot recover a vanished latent track.^[ref]

The diagram deliberately separates the preparation through thermodynamics from the readout decision. Its green interval is a stipulated operating window, not a prediction from pressure alone. Excess superheating can increase unwanted nucleation; insufficient superheating can suppress tracks. A useful proposed sensitivity check holds the readout settings fixed while varying pressure or temperature, with matched background exposures. A proposed timing check moves the beam within the cycle; a proposed optical check varies flash delay to distinguish growing bubbles from changes in spatial resolution. These are experimental design examples, not additional historical measurements.

Why can hot liquid wait for a charged particle rather than immediately turn into vapor? A small cavity costs surface energy. Consider an original macroscopic capillarity model at fixed temperature with surface tension $\sigma >0$, imposed liquid pressure $p_l$, and driving difference $\Delta p=p_{\mathrm{sat}}(T)-p_l>0$, where $p_{\mathrm{sat}}$ is the planar-interface vapor pressure. Neglect dissolved gas, interface thickness, heat transfer and curvature dependence of $\sigma$. The reversible formation-work approximation isThe first term charges for surface area; the second rewards the new vapor volume. This is a macroscopic free energy model, not the energy that an actual detector must receive as localized heat. It shares the capillarity balance used in the cloud-chamber expansion model, but the new phase is now a bubble.

Differentiating givesAt the critical nucleus, $d^2\Delta G/d{r}^2=-8\pi \sigma <0$. The barrier is a maximum: cavities smaller than $r_{\ast }$ tend to shrink, while larger ones can grow in this model. Doubling $\Delta p$ halves $r_{\ast }$ and quarters $\Delta G_{\ast }$ if $\sigma$ stays fixed. This explains the strong sensitivity to operating conditions without claiming a universal detector threshold.

Choose $\sigma =0.010\mathrm{J}{\mathrm{m}}^{-2}$ and $\Delta p=2.0\times 10^5\mathrm{Pa}$ solely for a teaching example. ThenFor the same chosen $\sigma$ and twice the driving pressure, $r_{\ast }=50\mathrm{nm}$ and $\Delta G_{\ast }\simeq 0.654\mathrm{keV}$. These are calculated formation-work scales, not measured bubble thresholds for diethyl ether or any other substance.

A simplified way to discuss localization is to stipulate deposited energy $E_{\mathrm{loc}}$ within a segment of length $\ell$ and a trial conversion fraction $0<\eta \le 1$. A toy criterion $\eta E_{\mathrm{loc}}\gtrsim \Delta G_{\ast }$ gives an energy per unit length scale $\Delta G_{\ast }/(\eta \ell )$. It is only a bookkeeping estimate: a real bubble also needs latent heat, must exchange heat, and is affected by the spatial pattern and time scale of energy deposition. Consequently neither a single stopping power value nor this barrier alone predicts detection efficiency or the number of bubbles. Two tracks with the same mean stopping power can distribute their energy differently on the scale of a critical nucleus.

Stereoscopy enables three-dimensional track reconstruction.^[ref] Calibrated camera positions and optical geometry are needed before a photographed arc becomes a physical radius of curvature. Neutral particles do not directly leave the same continuous ionization trail as charged particles; an interaction or decay can instead reveal them through charged secondary tracks. A photon's inferred path between interactions is not a photographed luminous line.

In a uniform magnetic field, with negligible stopping power and elastic scattering over the fitted segment, the transverse projection is a circle. The Lorentz force turns transverse momentum at angular rate $v_{\perp }/r$, givingHere $r$ is measured in the plane perpendicular to the magnetic field. The relation remains valid with relativistic physical momentum; it does not require $p=mv$. Thus curvature determines magnetic rigidity, not mass or kinetic energy alone. For chosen $|q|=e$, $B=1.5\mathrm{T}$ and $r=0.40\mathrm{m}$, $p_{\perp }\simeq 0.180\mathrm{GeV}/c$. This is an invented example, not a fit to a historical photograph. Reconstructing a three-dimensional helix also determines the component parallel to the magnetic field; the transverse radius of curvature alone does not. Determining the sign of electric charge additionally requires the field direction and the direction of travel, as in the charge-sign ambiguity of a curved track.

A bubble chamber using liquid hydrogen as both target and recording medium. For a beam interaction on a proton, the target's initial momentum can usually be approximated as zero on the scale of the beam momentum. This makes four-momentum conservation especially useful: the initial state does not include the unknown motion of a proton bound inside a heavier atomic nucleus. The approximation still neglects molecular binding and thermal motion.

Luis Walter Alvarez's program extended the Glaser bubble-chamber experiment into large hydrogen bubble chambers and systematic analysis. His lecture credits John G. Wood, an accelerator technician, with early hydrogen tracks, and emphasizes technicians' responsibility for chamber development.^[ref]

Alvarez's 1956 report describes a ten-inch hydrogen bubble chamber in a steady magnetic field, exposed to a $700\mathrm{MeV}$ negative-pion beam for proton scattering. Liquid hydrogen was the target. Horizontal windows admitted illumination and stereo viewing; a shielded light source let bubble-scattered light reach the lenses. A temperature gradient initially restricted sensitivity near the top. Tilting and fast recompression improved uniformity.^[ref]

The same report compares tracks without a magnetic field against straightness expectations, testing spurious curvature. Multiple Coulomb scattering remained a limitation. This separates an instrumental check from the physical momentum inferred with the field on.^[ref]

The scientific opportunity is to turn a visible geometry into a test of an invisible intermediate state. The bubble-chamber pressure cycle prepares the sensitive target; bubble-chamber track readout supplies geometric measurements. Neither step alone identifies a particle resonance.

Stereoscopy supplies depth because different cameras see different projections of the same track. A functional reconstruction first matches the same interaction vertex and track segments across views. A calibrated lens and window model maps each measured image coordinate back to a ray inside the liquid. Refraction at the windows must be included. Intersecting the rays, or finding their closest agreement when measurements have errors, estimates three-dimensional coordinates. Several points then determine a trajectory, rather than only a projected arc.

An original ideal-camera model makes the geometry explicit. Take two parallel pinhole cameras with centers $(0,0,0)$ and $(b,0,0)$, focal length $f$, and optical axes along $z$. Use virtual image planes in front of the pinholes so that image coordinates have the same orientation as the scene. A point $(X,Y,Z)$ then givesSubtracting the horizontal image coordinates gives disparity $d=u_L-u_R=fb/Z$, henceFor stipulated $b=0.20\mathrm{m}$, $f=0.050\mathrm{m}$, $u_L=0.010\mathrm{m}$, $u_R=0.005\mathrm{m}$ and $v_L=0.002\mathrm{m}$, the reconstructed point is $(0.40,0.08,2.0)\mathrm{m}$. These are ideal optical inputs, not dimensions of Alvarez's apparatus. This model omits refraction and lens distortion; actual calibration replaces its simple ray geometry. Since $|dZ/dd|=fb/d^2$, the same disparity error gives worse depth precision for distant points.

The Lorentz force then supplies a momentum scale. With magnetic field $B$ uniform and negligible stopping power over a fitted segment, bubble-chamber track readout gives $p_{\perp }=|q|Br$. The local three-dimensional helix tangent supplies the missing direction: if $\alpha$ is its angle to the field, $p=p_{\perp }/\sin \alpha$. In practical units,Here $e$ is the elementary charge and $c$ the speed of light. A short, almost straight arc poorly constrains its radius of curvature; a visually sharp track can therefore have poor momentum precision.

Alvarez describes a progression from measured coordinates to momenta and angles, then conservation-constrained fits testing alternative track identities. Computer programs made the repeated tests practical; an insufficiently constrained event could retain several hypotheses.^[ref]

A mathematical formulation explains what such a fit can establish. Let $y$ collect measured momenta and angles, with positive-definite covariance matrix $V$. For a chosen reaction hypothesis, adjust them to $x$ by minimizing a generalized least squares objective,The constraints $g$ express four-momentum conservation and any independently known daughter masses. For each proposed identity, calculate energy from the energy–momentum relation, $E_i^2=m_i^2{c}^4+c^2|{\mathbf{p}}_i{|}^2$. A missing neutral particle adds unknowns; a measured decay vertex can add geometric constraints. Compatibility depends on the measurement errors and the independent constraints, not just the smallest raw ${\chi }^2$ among unrelated fits.