ELOG Updates and Results

Amy Connolly , Other, Simulation, Plots for QC sites,

Tue, 27 Feb 2018 23:13:15 -0500

Hi All,

We need to revive the icemcQC and AraSimQC pages, and I thought that in the meantime we could keep a list and/or associated code for plots we'd like to see go in there.

Here are a couple plots I just made code attached is forstevenanddave.cc

event weight (attenuation factor) vs. angle that neutrino makes wrt up when standing under the anita balloon

event weight (attenuation factor) vs. angle that rf makes wrt up from an observer on the anita balloon

BC: Adding old Stephen Hoover plots from ANITA elog: https://www.phys.hawaii.edu/elog/anita_notes/32

Brian Clark, Analysis, Analysis, Testbed Channel Mapping and Antenna Information, ARA

Tue, 10 Oct 2017 11:04:05 -0400

This is the Testbed polarization channel mapping. This is the polarization result if you use the getGraphfromRFChan Function:

/* Channel mappings for the testbed
Channel 0: H Pol
Channel 1: H Pol
Channel 2: V Pol
Channel 3: V Pol
Channel 4: V Pol
Channel 5: H Pol
Channel 6: V Pol
Channel 7: H Pol
Channel 8: V Pol
Channel 9: H Pol
Channel 10: V Pol
Channel 11: H Pol
Channel 12: H Pol
Channel 13: H Pol
Channel 14: Surface
Channel 15: Surface
*/

Also, the Testbed has a somewhat bizarre menagerie of antennas. Here's how to understand it:

Check out the table of antennas for the testbed (table 1 in both papers): https://arxiv.org/pdf/1105.2854.pdf , https://arxiv.org/pdf/1404.5285.pdf

Basically the testbed was weird. There are four bowtie slotted cylinders deployed at ~30 m (these are the "deep hpol") and four bicones deployed at ~30 m (these are the "deep vpol"). So 8 total there: four V, four H.

Quote:

I am trying to write a script that will plot a 2d histogram of effective volume versus two of ShelfMC's parameters.

The script prompts the user for which two parameters (out of five that we vary in our parameter space scan) to plot along the x- and y-axes, as well as what values to hold the other 3 parameters constant at. It then collects the necessary root files from simulation results, generates a plotting script, and runs the plotting script to produce a plot in pdf form.

After many struggles I have the script written to the point where it functions, but the plots don't look right. Some plots look like they could be actual data (like Veff_A_I), and others just look flat-out wrong (like Veff_R_S).

I have yet to pin down the cause of this, but hopefully will be able to sometime in the near future.

Spoorthi Nagasamudram, Modeling, Simulation, A different way to implement ray tracing in AraSim (possibly even other simulations),

Sun, 17 Sep 2017 20:10:21 -0400

Hi,

I've attached a copy of some of the work I did over the summer on ray tracing and how I did it. Please let me know if you have any questions.

PS The plots are sidewards for some reasons. Sorry about that.

Amy Connolly (for Suren), Modeling, General, Using spherical harmonics to search for ideal antenna beam patters, Other

Fri, 15 Sep 2017 23:25:07 -0400

Suren put a bunch of info about the work he did this summer on github:

https://github.com/osu-particle-astrophysics/Spherical-Harmonics

He added this in an email:

I added more spherical harmonics into the Chi^2 code so now we can test up to L=12. Adding the additional harmonics brought the chi^2 for the first frequency gain fit down about 30%. However, I also fit the first phase, and it seems to be much harder to fit, especially at the poles. The Chi^2 in in the 200s. Therefore, I am unsure how you wish to proceed with the phases. Perhaps just avoiding the 0 degree and 5 degree theta cones is the way to go. This would involve modifying the chi^2 code and the data to not have those 0 and 5 degree cones.

Amy Connolly, Analysis, Analysis, Info on generating pseudoexperiments, calculating likelihoods from them and finding p-values, Other

Wed, 13 Sep 2017 09:28:15 -0400

Will point to a bunch of papers and stuff here.

Julie Rolla, Modeling, Other, Jordan's Code Antenna Optimization/Evolution, Other

Fri, 08 Sep 2017 16:51:28 -0400

Below is the explaination regarding this code from Jordan. Attached is the code in its last form.

Hi Lucas,

I'll try to fit a bunch of stuff into this email, so hopefully some of it is useful/interesting.

As far as the evolutionary algorithm goes, I'll send you the code and let you look at it yourself. the paper that I am developing it from comes from https://www.researchgate.net/publication/282857432_INTEGRATION_OF_REMCOM_XFDTD_AND_OCTAVE_FOR_THE_AUTOMATED_OPTIMAL_DESIGN_OF_MICROWAVE_ANTENNAS. My code has a couple differences in that I didn't look at the past k iterations like the paper does, but instead I have 5 parallel designs being run and look at how many of those improve the output. This is subtly different and might not be good, because it considers each design at the same time, so the mean doesn't have time to adjust before it is reevaluated if that makes sense. Anyways, just something to think about.

The basic structure for the code is that it is run with command line arguments, so that you compile it with whatever name ( I usually call it evolved_dipole, but doesn't matter). So it is run as

$./evolved_dipole --start

for the first time in a run and then every subsequent time in a run

$./evolve_dipole --cont

The --start will create the initial population, and record the parameters in a .csv file called handshake.csv. The --cont will read in a file called handshook.csv that theoretically will have the output of the antenna simulations for each antenna.

The first obvious thing I can think of that is missing in this script is that it doesn't write

Abdullah Alhag, Analysis, Analysis, GP algorithms,

Fri, 28 Jul 2017 17:43:20 -0400

In this post, I will be pointing out the advantage and the disadvantage of the GP algorithms I came across, particular Eureqa and HeuristicLab.

Eureqa is by far the fastest genetic algorithm software I came across. It is over simplified and easy to use. It has some built-in fitness function and also with some playing with the function that is being solved for and some other feature, it is possible for one to write his/her own fitness function. Moreover, the software is available for free for academic use and for most platform. Other features come with the software is the ability to normalize the data in different ways and even handle outliers and missing data. The program support a large collection of functions including trig and more complex one.

One the other hand, HeuristicLab is much slower than Eureqa but still far faster than Karoo-GP. The latest version of the software was released a year ago, and the support for the software is fairly slow. It is only supported for Windows; however, there is plans to adopted to Linux systems. The software support way more feature than Eureqa or karoo and even different regression and classification algorithms. You could also get the function ready to use in many software such as MATLAB, Excel, Mathematica, and much more. Another cool feature is that it shows you a three of the function and the weight of each node (operation or operand), greener means the node has more weight, see attached. It should be noticed that the software has the tendency to grow large three which could be fixed by changing the default max three length and the max three depth. The software has a problem with the last update of windows 10, you will get the blue screen if you opened too many windows, so be careful.

In booth software you could change how much of the data set goes to training and how much goes to testing, be sure to shuffle the data in HeuristicLab as it will otherwise distribute the data as training and testing non-randomly. Booth software by

Brian Clark and Ian Best, Lab Measurement, Hardware, Mac Addresses,

Wed, 05 Jul 2017 19:22:52 -0400

This is a "bank" of mac addresses that we obtained for the lab. They were taken by Ian Best (one of Jim' students) in summer 2017. We purchased 25 AT24MAC402 EEPROMS (https://www.arrow.com/en/products/at24mac402-stum-t/microchip-technology) and used a SOT-23 breakout board and a bus pirate to retrieve their internal mac addresses.

If you take one, please note where you used it so no one tries to take the same one twice:

Serial Number	EUI (Mac) Address	Used?
0x0A70080064100461105CA000A0000000	FC:C2:3D:0D:A6:71	Spare ADAQ (ADAQF002)
0x0A70080064100460FC6DA000A0000000	FC:C2:3D:0D:A7:37	Spare ADAQ (ADAQF003)
0x0A700800641004611C8DA000A0000000	FC:C2:3D:0D:A8:7C	Spare ADAQ (ADAQF004)
0x0A70080064100461E47FA000A0000000	FC:C2:3D:0D:BC:BF	-- reserved for testing
0x0A700800641004611C9FA000A0000000	FC:C2:3D:0D:BE:02
0x0A70080064100461D8D9A000A0000000	FC:C2:3D:0D:C0:2D
0x0A700800641004611C6EA000A0000000	FC:C2:3D:0D:C6:62
0x0A70080064100460F0A0A000A0000000	FC:C2:3D:0D:C8:AF
0x0A70080064100461F870A000A0000000	FC:C2:3D:0D:DC:22
0x0A700800641004612CC7A000A0000000	FC:C2:3D:0D:E9:F4	ARA3 ADAQ (ADAQG003)
0x0A70080064100461EC43A000A0000000	FC:C2:3D:0D:EF:78
0x0A7008006410046228B7A000A0000000	FC:C2:3D:0D:FE:3D	&nbs Brian Clark, Analysis, Analysis, Estimate of ARA Station-Year/ Livetime, ARA Mon, 24 Apr 2017 22:27:29 -0400 In response to a request by Amy, I make an estimae of the number of deep "station-years" of data obtained by ARA so far. This means roughly (# deep stations) * (# months livetime). This is very approximate, and only counts days where ARA has data in the storage vault on cobalt. It doesn't verify that cal pulsers are running, or that we actually have data for very hour of every day, etc. Not accounting for 2013 A1 data, I get the following estimate. All I did was "ls \| wc -l" on all of the data directories to count the number of days. ARA1: 285 days 2012 + 124 days 2014 + 29 days 2015 + 117 days 2016 + 0 days 2017 = 555 days total ARA2: 211 days 2013 + 310 days 2014 + 345 days 2015 + 314 days 2016 + 109 days 2017 = 1289 days total ARA 3: 214 days 2013 + 303 days 2014 + 251 days 2015 + 292 days 2016 + 0 days 2017 = 1060 days total So for the three stations that is 2904 days total, or ~8 station years of data. The spreadsheet with the calculation is attached, including which directories I searched over to make the count. Hannah Hasan, Other, Simulation, Plotting ShelfMC Parameter Space, Other Sun, 23 Apr 2017 14:54:50 -0400 I am trying to write a script that will plot a 2d histogram of effective volume versus two of ShelfMC's parameters. The script prompts the user for which two parameters (out of five that we vary in our parameter space scan) to plot along the x- and y-axes, as well as what values to hold the other 3 parameters constant at. It then collects the necessary root files from simulation results, generates a plotting script, and runs the plotting script to produce a plot in pdf form. After many struggles I have the script written to the point where it functions, but the plots don't look right. Some plots look like they could be actual data (like Veff_A_I), and others just look flat-out wrong (like Veff_R_S). I have yet to pin down the cause of this, but hopefully will be able to sometime in the near future. Brian Clark and Patrick Allison, Other, Hardware, ARAFE Master Documentation, ARA Fri, 14 Apr 2017 12:51:53 -0400 Here is documentation on the ARAFE Master firmware and software design for the next generation of ARA stations. I include both the pdf and tex source code. The firmware is located at: https://github.com/ara-daq-hw/ArafeMasterSoftware The software is located at: https://github.com/ara-daq-hw/ArafeMasterSoftware Revision History 2017.04.26: Typo fixes, fault curve addition, and python hex preparation instructions. Suren Gourapura and Brian Clark, Other, Hardware, ARAFE master Python communication, ARA Fri, 14 Apr 2017 12:51:27 -0400 With Brian's help, I am writing the python serial commander code used to control and troubleshoot the ARAFE master board. I have worked on it for about 2 weeks so far, and progress is slow but measureable! Recently, we are powering on a channel and are able to measure the clock on an ARAFE board we attach to it. You can check out our progress here: https://github.com/ara-daq-hw/arafe-master/blob/master/python_serial_commander.py

ELOG Updates and Results

Amy Connolly , Other, Simulation, Plots for QC sites,

Brian Clark, Analysis, Analysis, Testbed Channel Mapping and Antenna Information, ARA

Spoorthi Nagasamudram, Modeling, Simulation, A different way to implement ray tracing in AraSim (possibly even other simulations),

Amy Connolly (for Suren), Modeling, General, Using spherical harmonics to search for ideal antenna beam patters, Other

Amy Connolly, Analysis, Analysis, Info on generating pseudoexperiments, calculating likelihoods from them and finding p-values, Other

Julie Rolla, Modeling, Other, Jordan's Code Antenna Optimization/Evolution, Other

Abdullah Alhag, Analysis, Analysis, GP algorithms,

Brian Clark and Ian Best, Lab Measurement, Hardware, Mac Addresses,

Brian Clark, Analysis, Analysis, Estimate of ARA Station-Year/ Livetime, ARA

Hannah Hasan, Other, Simulation, Plotting ShelfMC Parameter Space, Other

Brian Clark and Patrick Allison, Other, Hardware, ARAFE Master Documentation, ARA

Suren Gourapura and Brian Clark, Other, Hardware, ARAFE master Python communication, ARA