12,457 research outputs found
A likelihood method to cross-calibrate air-shower detectors
We present a detailed statistical treatment of the energy calibration of
hybrid air-shower detectors, which combine a surface detector array and a
fluorescence detector, to obtain an unbiased estimate of the calibration curve.
The special features of calibration data from air showers prevent unbiased
results, if a standard least-squares fit is applied to the problem. We develop
a general maximum-likelihood approach, based on the detailed statistical model,
to solve the problem. Our approach was developed for the Pierre Auger
Observatory, but the applied principles are general and can be transferred to
other air-shower experiments, even to the cross-calibration of other
observables. Since our general likelihood function is expensive to compute, we
derive two approximations with significantly smaller computational cost. In the
recent years both have been used to calibrate data of the Pierre Auger
Observatory. We demonstrate that these approximations introduce negligible bias
when they are applied to simulated toy experiments, which mimic realistic
experimental conditions.Comment: 10 pages, 7 figure
Muon Track Reconstruction and Data Selection Techniques in AMANDA
The Antarctic Muon And Neutrino Detector Array (AMANDA) is a high-energy
neutrino telescope operating at the geographic South Pole. It is a lattice of
photo-multiplier tubes buried deep in the polar ice between 1500m and 2000m.
The primary goal of this detector is to discover astrophysical sources of high
energy neutrinos. A high-energy muon neutrino coming through the earth from the
Northern Hemisphere can be identified by the secondary muon moving upward
through the detector. The muon tracks are reconstructed with a maximum
likelihood method. It models the arrival times and amplitudes of Cherenkov
photons registered by the photo-multipliers. This paper describes the different
methods of reconstruction, which have been successfully implemented within
AMANDA. Strategies for optimizing the reconstruction performance and rejecting
background are presented. For a typical analysis procedure the direction of
tracks are reconstructed with about 2 degree accuracy.Comment: 40 pages, 16 Postscript figures, uses elsart.st
The wavefront of the radio signal emitted by cosmic ray air showers
Analyzing measurements of the LOPES antenna array together with corresponding
CoREAS simulations for more than 300 measured events with energy above
eV and zenith angles smaller than , we find that the radio
wavefront of cosmic-ray air showers is of approximately hyperbolic shape. The
simulations predict a slightly steeper wavefront towards East than towards
West, but this asymmetry is negligible against the measurement uncertainties of
LOPES. At axis distances m, the wavefront can be approximated by
a simple cone. According to the simulations, the cone angle is clearly
correlated with the shower maximum. Thus, we confirm earlier predictions that
arrival time measurements can be used to study the longitudinal shower
development, but now using a realistic wavefront. Moreover, we show that the
hyperbolic wavefront is compatible with our measurement, and we present several
experimental indications that the cone angle is indeed sensitive to the shower
development. Consequently, the wavefront can be used to statistically study the
primary composition of ultra-high energy cosmic rays. At LOPES, the
experimentally achieved precision for the shower maximum is limited by
measurement uncertainties to approximately g/cm. But the simulations
indicate that under better conditions this method might yield an accuracy for
the atmospheric depth of the shower maximum, , better than
g/cm. This would be competitive with the established air-fluorescence
and air-Cherenkov techniques, where the radio technique offers the advantage of
a significantly higher duty-cycle. Finally, the hyperbolic wavefront can be
used to reconstruct the shower geometry more accurately, which potentially
allows a better reconstruction of all other shower parameters, too.Comment: accepted by JCA
Modelling of the Surface Emission of the Low-Magnetic Field Magnetar SGR 0418+5729
We perform a detailed modelling of the post-outburst surface emission of the
low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this
source, B=6x10^12 G estimated from its spin-down rate, is in the observed range
of magnetic fields for normal pulsars. The source is further characterized by a
high pulse fraction and a single-peak profile. Using synthetic temperature
distribution profiles, and fully accounting for the general-relativistic
effects of light deflection and gravitational redshift, we generate synthetic
X-ray spectra and pulse profiles that we fit to the observations. We find that
asymmetric and symmetric surface temperature distributions can reproduce
equally well the observed pulse profiles and spectra of SGR 0418. Nonetheless,
the modelling allows us to place constraints on the system geometry (i.e. the
angles and that the rotation axis makes with the line of sight and
the dipolar axis, respectively), as well as on the spot size and temperature
contrast on the neutron star surface. After performing an analysis iterating
between the pulse profile and spectra, as done in similar previous works, we
further employed, for the first time in this context, a Markov-Chain
Monte-Carlo approach to extract constraints on the model parameters from the
pulse profiles and spectra, simultaneously. We find that, to reproduce the
observed spectrum and flux modulation: (a) the angles must be restricted to
or ; (b) the
temperature contrast between the poles and the equator must be at least a
factor of , and (c) the size of the hottest region ranges between
0.2-0.7 km (including uncertainties on the source distance). Last, we interpret
our findings within the context of internal and external heating models.Comment: 13 pages, 10 figures. Accepted for publication in MNRA
Prospects for the measurement of muon-neutrino disappearance at the FNAL-Booster
Neutrino physics is nowadays receiving more and more attention as a possible
source of information for the long-standing problem of new physics beyond the
Standard Model. The recent measurement of the mixing angle in the
standard mixing oscillation scenario encourages us to pursue the still missing
results on leptonic CP violation and absolute neutrino masses. However,
puzzling measurements exist that deserve an exhaustive evaluation. The NESSiE
Collaboration has been setup to undertake conclusive experiments to clarify the
muon-neutrino disappearance measurements at small , which will be able to
put severe constraints to models with more than the three-standard neutrinos,
or even to robustly measure the presence of a new kind of neutrino oscillation
for the first time. To this aim the use of the current FNAL-Booster neutrino
beam for a Short-Baseline experiment has been carefully evaluated. This
proposal refers to the use of magnetic spectrometers at two different sites,
Near and Far. Their positions have been extensively studied, together with the
possible performances of two OPERA-like spectrometers. The proposal is
constrained by availability of existing hardware and a time-schedule compatible
with the CERN project for a new more performant neutrino beam, which will
nicely extend the physics results achievable at the Booster. The possible FNAL
experiment will allow to clarify the current disappearance tension
with appearance and disappearance at the eV mass scale. Instead, a new
CERN neutrino beam would allow a further span in the parameter space together
with a refined control of systematics and, more relevant, the measurement of
the antineutrino sector, by upgrading the spectrometer with detectors currently
under R&D study.Comment: 76 pages, 52 figure
Compton Imaging of MeV Gamma-Rays with the Liquid Xenon Gamma-Ray Imaging Telescope (LXeGRIT)
The Liquid Xenon Gamma-Ray Imaging Telescope (LXeGRIT) is the first
realization of a liquid xenon time projection chamber for Compton imaging of
MeV gamma-ray sources in astrophysics. By measuring the energy deposit and the
three spatial coordinates of individual gamma-ray scattering points, the
location of the source in the sky is inferred with Compton kinematics
reconstruction. The angular resolution is determined by the detector's energy
and spatial resolutions, as well as by the separation in space between the
first and second scattering. The imaging response of LXeGRIT was established
with gamma-rays from radioactive sources, during calibration and integration at
the Columbia Astrophysics Laboratory, prior to the 2000 balloon flight mission.
In this paper we describe in detail the various steps involved in imaging
sources with LXeGRIT and present experimental results on angular resolution and
other parameters which characterize its performance as a Compton telescope.Comment: 22 pages, 20 figures, submitted to NIM
3D reconstruction of cerebral blood flow and vessel morphology from x-ray rotational angiography
Three-dimensional (3D) information on blood
flow and vessel morphology is important when
assessing cerebrovascular disease and when monitoring interventions. Rotational angiography
is nowadays routinely used to determine the geometry of the cerebral vasculature. To this end,
contrast agent is injected into one of the supplying arteries and the x-ray system rotates around
the head of the patient while it acquires a sequence of x-ray images. Besides information on the
3D geometry, this sequence also contains information on blood flow, as it is possible to observe
how the contrast agent is transported by the blood. The main goal of this thesis is to exploit
this information for the quantitative analysis of blood flow.
I propose a model-based method, called
flow map fitting, which determines the blood flow
waveform and the mean volumetric flow rate in the large cerebral arteries. The method uses a
model of contrast agent transport to determine the
flow parameters from the spatio-temporal
progression of the contrast agent concentration, represented by a flow map. Furthermore, it
overcomes artefacts due to the rotation (overlapping vessels and foreshortened vessels at some
projection angles) of the c-arm using a reliability map.
For the flow quantification, small changes to the clinical protocol of rotational angiography
are desirable. These, however, hamper the standard 3D reconstruction. Therefore, a new method
for the 3D reconstruction of the vessel morphology which is tailored to this application is also presented.
To the best of my knowledge, I have presented the first quantitative results for blood flow
quantification from rotational angiography. Additionally, the model-based approach overcomes
several problems which are known from flow quantification methods for planar angiography.
The method was mainly validated on images from different phantom experiments. In most
cases, the relative error was between 5% and 10% for the volumetric mean flow rate and between
10% and 15% for the blood flow waveform. Additionally, the applicability of the flow model was shown on clinical images from planar angiographic acquisitions. From this, I conclude that the method has the potential to give quantitative estimates of blood flow parameters during
cerebrovascular interventions
Advanced capabilities for planar X-ray systems
Mención Internacional en el título de doctorThe past decades have seen a rapid evolution towards the use of digital detectors
in radiology and a more flexible robotized movement of the system components,
X-ray tube and detector. This evolution opened the possibility for incorporating
advanced capabilities in these planar X-ray systems, and for providing new valuable
diagnostic information compared to the previous technology. Some of the current
challenges for radiography are to obtain more quantitative images and to reduce the
inherent superposition of tissues because of the 2D nature of the technique.
Dual energy radiography, based on the acquisition of two images at different
source voltages, enables a separate characterization of soft tissue and bone structures.
Its benefits over conventional radiography have been proven in different applications,
since it improves information content without adding significant extra
acquisition time or radiation dose.
In a different direction, a really disruptive advance would be to obtain 3D imaging
with systems designed just for planar images. The incorporation of tomographic
capabilities into these systems would have to deal with the acquisition of a limited
number of projections, with non-standard geometrical configurations.
This thesis presents original contributions in these two directions: dual energy
radiography and 3D imaging with X-ray systems designed for planar imaging. The
work is framed in a line of research of the Biomedical Imaging and Instrumentation
Group from the Bioengineering and Aerospace Department of University Carlos III
de Madrid working jointly with the University Hospital Gregorio Marañón, focused
on the advance of radiology systems. This research line is carried out in collaboration
with the group of Computer Architecture, Communications and Systems (ARCOS),
from the same university, the Imaging Research Laboratory (IRL) of the University
of Washington and the research center CREATIS, France. The research has a clear
focus on technology transfer to the industry through the company Sedecal, a Spanish
multinational among the 10 best world companies in the medical imaging field.
The first contribution of this thesis is a complete novel protocol to incorporate
dual energy capabilities that enable quantitative planar studies. The proposal is
based on the use of a preliminary calibration with a very simple and low-cost phantom
formed by two parts that represent soft tissue and bone equivalent materials.
This calibration is performed automatically with no strict placement requirements.
Compared to current Dual-energy X-ray Absorptiometry (DXA) systems, 1) it provides
real mass-thickness values directly, enabling quantitative planar studies instead
of relative comparisons, and 2) it is based on an automatic preliminary calibration without the need of interaction of an experienced technician.
The second contribution is a novel protocol for the incorporation of tomographic
capabilities into X-ray systems originally intended for planar imaging. For this purpose,
we faced three main challenges.
First, the geometrical trajectory of equipment follows non-standard circular orbits,
thus posing severe difficulties for reconstruction. To handle this, the proposed
protocol comprises a new geometrical calibration procedure that estimates all the
system parameters per-projection.
Second, the reconstruction of a limited number of projections from a reduced angular
span leads to severe artifacts when using conventional reconstruction methods.
To deal with these limited-view data, the protocol includes a novel advanced reconstruction
method that incorporates the surface information of the sample, which
can be extracted with a 3D light surface scanner. These data are introduced as an
imposed constraint following the Split Bregman formulation. The restriction of the
search space by exploiting the surface-based support becomes crucial for a complete
recovery of the external contour of the sample and surroundings when the angular
span is extremely reduced. The modular, efficient and flexible design followed for its
implementation allows for the reconstruction of limited-view data with non-standard
trajectories.
Third, the optimization of the acquisition protocols has not yet explored with
these systems. This thesis includes a study of the optimum acquisition protocols
that allowed us to identify the possibilities and limitations of these planar systems.
Using the surface-constrained method, it is possible to reduce the total number of
projections up to 33% and the angular span down to 60 degrees.
The contributions of this thesis open the way to provide depth and quantitative
information very valuable for the improvement of radiological diagnosis. This could
impact considerably the clinical practice, where conventional radiology is still the
imaging modality most used, accounting for 80-90% of the total medical imaging
exams. These advances open the possibility of new clinical applications in scenarios
where 1) the reduction of the radiation dose is key, such as lung cancer screening or
Pediatrics, according to the ALARA criteria (As Low As Reasonably Achievable),
2) a CT system is not usable due to movement limitations, such as during surgery
or in an ICU and 3) where costs issues complicate the availability of CT systems,
such as rural areas or underdeveloped countries.
The results of this thesis has a clear application in the industry, since it is part
of a proof of concept of the new generation of planar X-ray systems that will be
commercialized worldwide by the company SEDECAL (Madrid, Spain).Los últimos años están viendo un rápido avance de los sistemas de radiología hacia el
uso de detectores digitales y a una mayor flexibilidad de movimientos de los principales
componentes del sistema, el tubo de rayos X y el detector. Esta evolución abre
la posibilidad de incorporar capacidades avanzadas en sistemas de imagen plana por
rayos X proporcionando nueva información valiosa para el diagnóstico. Dos retos en
radiografía son obtener imágenes cuantitativas y reducir la superposición de tejidos
debida a la naturaleza proyectiva de la técnica.
La radiografía de energía dual, basada en la adquisición de dos imágenes a diferente
kilovoltaje, permite obtener imágenes de tejido blando y hueso por separado.
Los beneficios de esta técnica que aumenta la cantidad de información sin añadir
un tiempo de adquisición o de dosis de radiación extra significativos frente al uso de
radiografía convencional, han sido demostrados en diferentes aplicaciones.
En otra dirección, un avance realmente disruptivo sería la obtención de imagen
3D con sistemas diseñados únicamente para imagen plana. La incorporación de capacidades
tomográficas en estos sistemas tendría que lidiar con la adquisición de un
número limitado de proyecciones siguiendo trayectorias no estándar.
Esta tesis presenta contribuciones originales en esas dos direcciones: radiografía
de energía dual e imagen 3D con sistemas de rayos X diseñados para imagen plana.
El trabajo se encuadra en una línea de investigación del grupo de Imagen Biomédica
e Instrumentación del Departamento de Bioingeniería e Ingeniería Aerospacial de
la Universidad Carlos III de Madrid junto con el Hospital Universitario Gregorio
Marañon, centrada en el avance de sistemas de radiología. Esta línea de investigación
se desarollada en colaboración con el grupo Computer Architecture, Communications
and Systems (ARCOS), de la misma universidad, el grupo Imaging Research Laboratory
(IRL) de la Universidad de Washington y el centro de investigación CREATIS,
de Francia. Se trata de una línea de investigación con un claro enfoque de transferencia
tecnológica a la industria a través de la compañía SEDECAL, una multinacional
española de entre las 10 líderes del mundo en el campo de la radiología.
La primera contribución de esta tesis es un protocolo completo para incorporar
capacidades de energía dual que permitan estudios cuantitativos de imagen plana.
La propuesta se basa en una calibración previa con un maniquí simple y de bajo coste
formado por dos materiales equivalentes de tejido blando y hueso respectivamente.
Comparado con los sistemas actuales DXA (Dual-energy X-ray Absorptiometry),
1) proporciona valores reales de tejido atravesado, 2) se basa en una calibración
automática que no requiere la interacción de un técnico con gran experiencia. La segunda contribución es un protocolo nuevo para la incorporación de capacidades
tomográficas en sistemas de rayos X originariamente diseñados para imagen
plana. Para ello, nos enfrentamos a tres principales dificultades.
En primer lugar, las trayectorias que pueden seguir la fuente y el detector en
estos sistemas no constituyen órbitas circulares estándares, lo que plantea retos importantes
en la caracterización geométrica. Para solventarlo, el protocolo propuesto
incluye una calibración geométrica que estima todos los parámetros geométricos del
sistema para cada proyección.
En segundo lugar, la reconstrucción de un número limitado de proyecciones
adquiridas en un rango angular reducido da lugar a artefactos graves cuando se
reconstruye con algoritmos convencionales. Para lidiar con estos datos de ángulo
limitado, el protocolo incluye un nuevo método avanzado de reconstrucción que incorpora
la información de superficie de la muestra, que se puede se obtener con un
escáner 3D. Esta información se impone como una restricción siguiendo la formulación
de Split Bregman, para compensar la falta de datos. La restricción del espacio
de búsqueda a través de la explotación del soporte basado en superficie, es crucial
para una recuperación completa del contorno externo de la muestra cuando el rango
angular es extremadamente pequeño. El diseño modular, eficiente y flexible de la
implementación propuesta permite reconstruir datos de ángulo limitado obtenidos
con posiciones de fuente y detector no estándar.
En tercer lugar, hasta la fecha, no se ha explorado la optimización del protocolo
de adquisición con estos sistemas. Esta tesis incluye un estudio de los protocolos
óptimos de adquisición que permitió identificar las posibilidades y limitaciones de
estos sistemas de imagen plana. Gracias al método de reconstrucción basado en
superficie, es posible reducir el número total de proyecciones hasta el 33% y el rango
angular hasta 60 grados.
Las contribuciones de esta tesis abren la posibilidad de proporcionar información
de profundidad y cuantitativa muy valiosa para la mejora del diagnóstico radiológico.
Esto podría impactar considerablemente en la práctica clínica, donde la radiología
convencional es todavía la modalidad de imagen más utilizada, abarcando el 80-
90% del total de los exámenes de imagen médica. Estos avances abren la posibilidad
de nuevas aplicaciones clínicas en escenarios donde 1) la reducción de la dosis de
radiación es clave, como en screening de cáncer de pulmón, de acuerdo con el criterio
ALARA (As Low As Reasonably Achievable), 2) no se puede usar un sistema
TAC por limitaciones de movimiento como en cirugía o UCI, o 3) el coste limita la
disponibilidad de sistemas TAC, como en zonas rurales o en países subdesarrollados.
Los resultados de esta tesis presentan una clara aplicación industrial, ya que
son parte de un prototipo de la nueva generación de sistemas planos de rayos X que
serán distribuidos mundialmente por la compañía SEDECAL.This thesis has been developed as part of several research projects with public funding:
- DPI2016-79075-R. ”Nuevos escenarios de tomografía por rayos X”, IP: Mónica
Abella García, Ministerio de Economía y Competitividad, 01/01/2017-31/12/2019,
147.620 e.
- ”Nuevos escenarios de tomografía por rayos X (NEXT) DPI2016-79075-R.
Ministerio de Economía”, Industria y Competitividad. (Universidad Carlos
III de Madrid). 30/12/2016-29/12/2019. 147.620 e.
(…)
- FP7-IMI-2012 (GA-115337), ”PreDict-TB: Model-based preclinical development
of anti-tuberculosis drug combinations”. FP7-IMI - Seventh Framework
Programme (EC-EFPIA). Unión Europea. (Universidad Carlos III de Madrid).
01/05/2012-31/10/2017.
(…)
- TEC2013-47270-R, ”Avances en Imagen Radiológica (AIR)”, Ministerio de
Economía y Competitividad”, 01/01/2014-31/12/2016. IP: Mónica Abella Garcia
and Manuel Desco Menéndez. 160.204 e
(…)
- RTC-2014-3028-1, ”Nuevos Escenarios Clínicos con Radiología Avanzada (NECRA)”,
Ministerio de Economía y Competitividad, 01/06/2014-31/12/2016 IP: Mónica
Abella García. 2014-2016. 219.458,96 e
- IDI-20130301, ”Nuevo sistema integral de radiografía (INNPROVE: INNovative
image PROcessing in medicine and VEterinary)”, IP: Mónica Abella García
and Manuel Desco Menéndez. Ministerio de Economía y Competitividad.
Subcontratación CDTI, 14/01/2013-31/03/2015. Total: 1.860.629e (UC3M:
325.000e). (Art. 83)
- IPT-2012-0401-300000 INNPACTO 2012, ”Tecnologías para Procedimientos
Intraoperatorios Seguros y Precisos. XIORT. MINECO. (Universidad Carlos
III de Madrid). 01/01/2013-31/12/2015.Programa Oficial de Doctorado en Ingeniería MatemáticaPresidente: Doménec Ros Puig.- Secretario: Cyril Riddell.- Vocal: Yannick Boursie
- …