250 likes | 407 Views
Aesop An End-to-End Time-Domain ANITA Event Simulator. Stephen Hoover ANITA Collaboration Meeting December 2007 University of Hawaii. Aesop. Aesop is the time domain end-to-end ANITA event simulator developed at UCLA Available for general use – see ANITA note #370 for info
E N D
AesopAn End-to-End Time-Domain ANITA Event Simulator Stephen Hoover ANITA Collaboration Meeting December 2007 University of Hawaii
Aesop • Aesop is the time domain end-to-end ANITA event simulator developed at UCLA • Available for general use – see ANITA note #370 for info • Designed to be modular and easy to modify • Easy to update for ANITA-II • Easy to plug into other code Stephen Hoover, ANITA Collaboration Meeting 2007
Using Aesop in Standalone Mode • Separate package to interface with Aesop (EventSimulation program) • Can generate Askaryan signals, ground pulser signals, pure noise, all based on command-line inputs • Outputs to text file, displayed by ROOT script • Can read list of many event parameters from input file (e.g. from triggered events in icemc) • Output full events to ROOT file Stephen Hoover, ANITA Collaboration Meeting 2007
Using Aesop with Other Code • Create an “Aesop” object • Set event parameters (neutrino energy, payload orientation, etc.) • Out comes an “AnitaEventData” object • “AnitaEventData” extends Ryan’s “UsefulAnitaEvent” • So any analysis code that uses Ryan’s event reader can take output directly from Aesop • Can plug into Monte Carlos to get time-domain representations of events Stephen Hoover, ANITA Collaboration Meeting 2007
Aesop Features Askaryan pulse from theory Ground pulser signal Your signal Frequency-dependent changes to impulse response based on incident angle. From UH data taken by Predrag Miocinovic (ELOG #71). Antenna model System noise from measurements in ELOG #119. Also uses thermal noise from antenna and has option for non-thermal RFI. + Noise x 64 channels System impulse response from ELOG #121 and #119, SURF gain from ELOG #137. Convolve with impulse responses Output as ANITA event Useable by any code that uses Ryan’s EventReader. Stephen Hoover, ANITA Collaboration Meeting 2007
Simulation of Askaryan Signal • Class “AskaryanPulse” (extends “Signal”) • Uses parameterization from Jaime Alvarez-Muñiz in astro-ph/0512337 • Code from icemc • Phases +90° at pos. freq., -90° at neg. freq. Askaryan Signal 1° off-cone Stephen Hoover, ANITA Collaboration Meeting 2007
Simulation of Ground Pulser Signals • Class “GPSignal” (extends “Signal”) • Uses input files with lab-measured ground pulser signals • Lab measurements combined with Seavey antenna response to get E-field (see next slide) McMurdo Seavey Impulse Stephen Hoover, ANITA Collaboration Meeting 2007
Antenna Transmission • From a paper by Farr & Baum (used by Ped): (hN,TX is the antenna impulse response) • This works perfectly for Seavey tonebursts, but I have to multiply Vsrc by ~2.75 for Seavey impulses • Problem with measurement of Seavey impulse? • Not yet resolved Stephen Hoover, ANITA Collaboration Meeting 2007
ANITA Antenna Model • Class “AnitaAntenna” • Applies measured antenna impulse response • Modifies impulse response based on incident angle of signal Seavey Impulse Response Stephen Hoover, ANITA Collaboration Meeting 2007
Noise • Thermal / System noise in “AnitaAntenna” • Nonthermal sources added as signals from class “Signal” Stephen Hoover, ANITA Collaboration Meeting 2007
Signal (V) Signal Chain • Class “AnitaAntenna” • Measured impulse through all RF components from back of antenna to SURF input + attenuation in SURF • Slightly different impulse responses used for each channel Stephen Hoover, ANITA Collaboration Meeting 2007
Full-Payload Effects • Class “AnitaPayload” • Geometric time delay • Orientation of payload • Direction of incoming RF signal • Random jitter between and within phi sectors • Currently 35 ps within a phi sector • 65 ps between phi sectors Stephen Hoover, ANITA Collaboration Meeting 2007
Aesop Askaryan Signal 3x1019 eV neutrino, 40% of energy in hadronic shower, viewed from 1 degree off-cone. Signal traveled through 1 attenuation length of ice and 120 km of air. Stephen Hoover, ANITA Collaboration Meeting 2007
Trigger Simulation • Under construction! • Currently works by • Convert v-pol and h-pol SURF signals to LCP and RCP • Putting LCP & RCP signals through square bandpass filters • Applying simplistic exponential diode model • Checking for time coincidences & applying trigger logic • Downside of present implementation: • Slow (process ~5 events per minute) • Doesn’t quite work yet Stephen Hoover, ANITA Collaboration Meeting 2007
Verification • Three ways to compare Aesop to reality: • Thermal noise • Shape • Magnitude • Ground pulser (McMurdo Seavey) • Toneburst good • Problem with PCD signal? • SLAC • Work by David Goldstein • Quality of comparison currently limited by knowledge of geometry Stephen Hoover, ANITA Collaboration Meeting 2007
Red: Ch. 12H Black: Aesop, 160 K ant. noise Red: Ch. 15H Black: Aesop, 220 K ant. noise Comparing Thermal Noise Averaging FFTs of noise from events #6390500 to #6392400 in the 10% data (a quiet period) Stephen Hoover, ANITA Collaboration Meeting 2007
Final Words • Latest versions of Aesop available at http://home.physics.ucla.edu/~hoover/ • Hopefully on SVN server soonish • Please write your own classes and methods, and/or let me know what you’d find useful • No need to reinvent the wheel • (It’s a round sort of a thing on a stick, good for moving things around…) Stephen Hoover, ANITA Collaboration Meeting 2007
Additional Slides Stephen Hoover, ANITA Collaboration Meeting 2007
Extras • Included at no extra charge: • Deconvolution – using Aesop’s ANITA model to remove ANITA from the data • Analysis tools • AnitaEventData • Analytic signal, FFTs, Resampled waveforms • Software filter • Much, much more! Stephen Hoover, ANITA Collaboration Meeting 2007
Deconvolution Left: Askaryan pulse (as on p. 13 and p. 6) with no noise. Right: Waveform on left, deconvolved Stephen Hoover, ANITA Collaboration Meeting 2007
Deconvolution (Envelopes) Left: Envelope of Askaryan pulse from previous page. Right: Envelope of deconvolution Rise time (50% to 100%): 0.7 ns Rise time (50% to 100%): 1.05 ns Stephen Hoover, ANITA Collaboration Meeting 2007
Total Impulse Response Stephen Hoover, ANITA Collaboration Meeting 2007
Toneburst peak agrees to within 0.5 dBm/MHz! McM CW RFI Toneburst Comparison 300 MHz toneburst from McMurdo Seavey. Red: ANITA-I data Black: Aesop simulation McMurdo CW noise has been added to Aesop by hand, but there are no other manual adjustments Stephen Hoover, ANITA Collaboration Meeting 2007
Testing Aesop, part 1 • Aesop is the ANITA time-domain event simulator developed at UCLA • Signal of known magnitude and shape from ground pulser allows for check of gains and impulse responses in the simulation • Scale of simulated signal set matches scale of data without need for tuning! 600 MHz Toneburst ANITA-I Data Aesop Simulation Stephen Hoover, ANITA Collaboration Meeting 2007
Testing Aesop, part 2 Red: ANITA-I Event #266103, 23V: McM Seavey impulse, 10 dB attenuation Black: Aesop simulation of this event (but multiplied by fudge factor of 2.75!) Stephen Hoover, ANITA Collaboration Meeting 2007