GLAST ADASS 2005 Oct 4 1/21 GLAST Large Area GLAST Large Area Telescope: Telescope: LAT Science Analysis Software’s Use of Gaudi & ROOT Heather Kelly NASA Goddard Space Flight Center Science Systems and Applications Inc. For the GLAST LAT Collaboration SAS [email protected]Gamma-ray Large Gamma-ray Large Area Space Area Space Telescope Telescope http://www-glast.slac.stanford.edu/software
21
Embed
GLAST Large Area Telescope: LAT Science Analysis Software’s Use of Gaudi & ROOT Heather Kelly
Gamma-ray Large Area Space Telescope. GLAST Large Area Telescope: LAT Science Analysis Software’s Use of Gaudi & ROOT Heather Kelly NASA Goddard Space Flight Center Science Systems and Applications Inc. For the GLAST LAT Collaboration SAS [email protected]. - PowerPoint PPT Presentation
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
GLAST ADASS 2005 Oct 4
1/21
GLAST Large Area Telescope:GLAST Large Area Telescope:
LAT Science Analysis Software’s Use of Gaudi & ROOT
Heather KellyNASA Goddard Space Flight CenterScience Systems and Applications Inc.For the GLAST LAT Collaboration [email protected]
Gamma-ray Large Gamma-ray Large Area Space Area Space TelescopeTelescope
http://www-glast.slac.stanford.edu/software
GLAST ADASS 2005 Oct 4
2/21
GLAST MissionGLAST Mission
GLAST measures the direction, energy and arrival time of celestial gamma rays
-LAT measures gamma-rays in the energy range ~20 MeV - >300 GeV
- There is no telescope now covering this range!!
- GBM provides correlative observations of transient events in the energy range ~20 keV – 20 MeV
Launch: August 2007 Florida
Orbit: 565 km, 28.5o inclination
Lifetime: 5 years (minimum; 10 yrs goal)
GLAST ADASS 2005 Oct 4
3/21
GLAST ParticipationGLAST Participation
NASA - DoE Partnership on LAT
LAT is being built by an international teamStanford University (SLAC & HEPL, Physics)Goddard Space Flight CenterNaval Research LaboratoryUniversity of California, Santa CruzUniversity of WashingtonOhio State UniversityCEA/Saclay & IN2P3 (France)INFN & ASI (Italy)Hiroshima University, ISAS, RIKEN (Japan)Royal Inst. of Technology & Stockholm Univ. (Sweden)
GBM is being built by US and GermanyMPE, Garching (Germany)Marshall Space Flight Center
Single Photon Angular Resolution3.5o @ 100 MeV0.15o @ 10 GeV
Good Energy ResolutionE/E ~ 10%; 100 MeV – 10 GeV ~ < 20%; 10 GeV – 300 GeV
GLAST ADASS 2005 Oct 4
6/21
GLAST Concept
• Low profile for wide f.o.v.
• Segmented anti-detector to minimize self-veto at high E.
• Finely segmented calorimeter for enhanced background rejection and shower leakage correction.
• High-efficiency, precise track detectors located close to the conversions foils to minimize multiple-scattering errors.
• Modular, redundant design.
• No consumables.
• Low power consumption (650 W)
Calorimeter (energy measurement)
Particle tracking detectors
Conversion foils (W)
Charged particle anticoincidence shield
e+ e-
Pair production is the dominant photon interaction in our energy range
detection – pair conversion telescope
GLAST ADASS 2005 Oct 4
7/21
A Special NoteA Special Note
• This talk discusses the LAT data format and software.• FITS and FTOOLs will be delivered to the user community.
• Science End-Users view GLAST data as photon lists, tend to be taken in aggregate.
• LAT team looks at data on an event by event basis.– Upon triggering GLAST freezes its detectors and writes out
its buffers.– Once on the ground, the game is to separate gammas from
background. (1000:1)• Determine best estimate for direction and energy.
– Part of the LAT collaboration’s job is to provide the highest quality photon lists to the science user community.
GLAST ADASS 2005 Oct 4
8/21
Some BackgroundSome Background
• ~25 developers distributed across 9 time zones– VRVS (www.vrvs.org) and Instant Messaging (ICQ/Gaim)– Face-to-Face meetings twice a year– Close-knit team
• Support Windows, Red Hat Linux, (perhaps Mac OS X)– Why Windows??? Visual Studio Development Environment
• Currently ~35 MB of source code for:– Detailed Monte Carlo Simulation– Reconstruction and Calibration algorithms– I/O infrastructure
• GLAST LAT software is open and available to anyone who wants it.– Brookhaven group - KOPIO
• C++ was adopted for coding efforts – early on
• By 1999, we had detailed Monte Carlo simulation and prototype reconstruction code.– This was before the proposal!!!
GLAST ADASS 2005 Oct 4
9/21
Why Use a Framework?Why Use a Framework?
• By 1999, we had plenty of code – but it was unmaintainable.• Our MC generator was custom built, planned move to Geant4
– Looked to re-use existing code as much as possible• HEP offered a number of packages: CLHEP, Geant4
• Small team, of distributed programmers, with varying amounts of OO C++ programming experience.
• Provide one code system for simulation, test data analysis, and flight operations.
• Needed organization and clear lines of division between the components of our code.
• Flexible• Extensible
GLAST ADASS 2005 Oct 4
10/21
Gaudi FrameworkGaudi Framework
• An Object Oriented C++ Framework– “The implementation of an architecture which defines a structure
flexible enough to support all types of physics data processing needs..from simulation to analysis to visualization.”
• http://proj-gaudi.web.cern.ch/proj-gaudi/welcome.html• Promotes the creation of small maintainable components which can
be loaded at runtime – Components = algorithms, services, I/O, etc
• Standard Gaudi interface for each type of component• Embraces the division of data from algorithms• Gaudi provides a number of basic services such as:
– Event Data Service (TDS) – a fancy OO common block!– Messaging and Logging Services– Infrastructure for I/O where various output formats may be
supported• Provides standard event loop• Job parameters are handled via an input ASCII file• Optional Python interface
GLAST ADASS 2005 Oct 4
11/21
Example of Using Gaudi ToolsExample of Using Gaudi Tools
• http://root.cern.ch (circa 1995)• A floor wax and a dessert topping
– Machine independent file format– An analysis toolkit written in C++– Created to address the needs of large scale data analysis
• Object I/O• Self-describing files• Utilizes gzip algorithm for compression (user controlled)• Typically store event data in trees• C/C++ interpreter for command line analysis• PyROOT – python interface
GLAST ADASS 2005 Oct 4
16/21
GLAST ADASS 2005 Oct 4
17/21
Pretty PicturesPretty Pictures
GLAST ADASS 2005 Oct 4
18/21
Why ROOT?Why ROOT?
• As of 1999, our output was tab-delimited ASCII ntuples• GLAST code is Object-Oriented C++• A OO friendly file format seemed better suited for our needs
– Store our objects “as is” in ROOT• Either read them back into our sim/analysis framework• Or use ROOT/JAS or any ROOT capable analysis toolkit
to perform analysis.• Open-source, widely used within the HEP community• Well supported and good user documentation• Pledged backward compatibility• Development is actively extending and improving the system
– Planning support for multi-core – Adding context-sensitive help
GLAST ADASS 2005 Oct 4
19/21
Persistent ROOT OutputPersistent ROOT Output
• Four types of ROOT files: MC, Digitization, Reconstruction and ntuple.
• Two Levels of Output at end of Reconstruction:– Root Trees
• Basically, all the output of all steps of reconstruction• Enough information to read back in and continue
reconstruction from that point• Detailed offline analysis for reconstruction algorithm
improvements• Main component of System Tests
– Output Ntuple• A detailed branch which contains enough information
for checking of reconstruction performance• The analysis data which is ultimately FITisfied.
GLAST ADASS 2005 Oct 4
20/21
Advantages and DisadvantagesAdvantages and Disadvantages