This entry will serve to summarize what tasks we've completed and still need to complete after last week's GENETIS hackathon.
Alex
My focus was on the XF scripting for making the PUEO antenna. Before last week, we had a script which made the geometry of a PUEO-like antenna. The xmacro script has been modified in the following ways:
- Now includes power sources connecting the ridges that are opposite each other.
- Generalized to use variables without hardcoding.
- Reads in data from the output files of the PUEO GA.
- Loops over all read in antennas to set up all of the simulations at once.
At this point, the script for setting up the simulations of the PUEO loop is just about ready for use. The only change to be made is to the frequency list. The script for printing out the gain data is also ready, though it only prints out the gain in theta and phi, not the cross polarization.
Here are some outstanding points that need to be hammered out, though they impede the quality and accuracy of a run rather than the actual functioning of the loop.
- The shape of the ridges in the generated antennas from the GA are considerably different from the PUEO antenna.
- This is likely something that can be fixed by adjusting the constraints--the subspace of the parameter space that leads to PUEO-like ridges may be very small as-is.
- You can see an example of an antenna from the GA attached (viewed from slightly below). It looks dramatically different. This seems to be due to constraints on the height and width being too small, but the ridges are also clearly misshapen.
- The curve on the ridges looks like it needs some work.
- For very PUEO-like antennas, the slopes look reasonable, but they do flare in or out slightly. This is more dramatic for ridges that deviate from PUEO's significantly.
- My best guess is that the parametric form of the curves is placing the curves in a plane that is different from the rest of the ridge. This should be fixable by somehow constraining the curve to a given plane.
- The simulation still only outputs polarization data for theta and phi, not the cross polarization.
- It looks like XF naturally gives us the theta and phi polarization. It's not clear how to get the cross polarization gain (which we think means the gain when the signal is at 45 degrees).
- One idea is to redo the simulation but with the antenna rotated 45 degrees in the azimuth. But I'm not certain that that iwll work.
- Amy indicated that the correct way to do this is by simulating the antenna with just one power source at a time.
- Leaving the VPol source would produce the VV and VH data and leaving the HPol would produce the HH and HV data (HV and VH could be flipped!).
- Currently, the antennas have a fixed opening angle. We'd like to add this as a parameter, though at the moment it might be ok to forego that.
- We need to make sure that the power sources aren't too close together (fixed!).
To do for this week:
1. Modify the script to make two simulations for each antenna, each with only one power source.
2. Adjust the constriants on the parameters to fix the curvature.
3. When we begin a run, we will want to just evolve the height and the inner side length (these determine both the inner side length and the outer side length).
Dylan and Jacob
Dylan and Jacob focused on getting PUEOsim/icemc working, which will be the simulation software needed for the PUEO loop. Here's a summary from Dylan:
- Plotting and bash scripts should be mostly good to go.
- Currently have 4/5 of the PUEO prerequisites for pueoBuilder, and the GENETIS-specific code is ready to implement and test when it's built. I’m trying to build the correct root version one more time before today’s meeting, and if it fails again I’m going to ask Will.
- The main thing left is to get the correct version of root installed, which while annoying, shouldn't take too long
To do for this week:
1. Once PUEOsim is working, run it with the usual settings to get the default performance. It will need to be submitted as a job array (ask Alex for help with this if you need; ask Will for advice on how many neutrinos to run).
|