Updated To-Do List:

1. Add ARA constants in plotting code: ie add Veff, theta, R for ARA. This will sit as a constant line for us to compare each evolving parameter to. For example, in the plot Veff vs Generation, we will see a constant line for ARA's Veff to see how it is evolving in time. 

TASK 1: MAKE Veff PLOT. Note: One of our variables in this run is Energy. Unfortunately, we care about how Veff changes with Energy for ARA_actualBicone. This means that before plotting "Veff versus Generation" we need to generate Veff for ARA_actualBicone at the energy we are assigning the run. 

STEP 1: Make plotting code to plot "Veff vs Generation". 

• You can follow the way we make LRPlot.py.
  • However, in this case, we will need to take the files /fs/project/PAS0654/BiconeEvolutionOSC/BiconeEvolution/current_antenna_evo_build/XF_Loop/Evolutionary_Loop/Run_Ouputs/$RunName/AraOut${gen}_${i}.txt -- ie the AraSim outputs for each individual in ALL generations. ie Gen 0 - ${gen}  -- INCLUDING the one we created in Gen 1 for ARA_actualBicone -- and plot Veff from those. 
  • We also want to plot error bars with Veff. The error bars for each are labeled "and Veff (water eq)" at the bottom of these AraSim output files we are already reading in for Veff. 
    • Note, to know how to read in Veff, you can look at how we are reading that *exact* line from the AraSim output files in on the fitnessFunction_ARA.cpp code in the Antenna_Performance_Metric directory, since we are doing it there. 
    • Note that since we are plotting "Veff vs Generation", we will either have to read in all the files AraOut${gen}_${i}.txt (up to how may ${gen} we are in) every time -- including the ARA_actualBicone one-- or we can append them all to one .txt file. It's your choice. 
  • Make sure you label the units for each axis every time, and that you have a key! If not, Amy will hate the plot!

STEP2: Add V_eff for ARA_actualBicone as a constat in the newly made "Veff vs Gen" plotting code. Remember, as I mentioned above, we must first put the ARA_actualBicone through AraSim to see Veff at the same energy as the rest of the individuals. Here's what to do:

•  Edit the Bash script to submit a job that runs the gain patterns of ARA_actualBicone through AraSim in the initial gen (gen 0)! 
  • There's a current file that exits for ARA_actualBicone that has the gain patterns of the current antenna, in the correct format for AraSim to run. In the bash script (XF_Loop.sh) we need to create ONE more AraSim job submission for this antenna (as a comparison to what we see with our genetically evolved ones) just in gen 0 (the first generation).
    • This current AraSim input file for ARA_actualBicone exists in BiconeEvolutionOSC/AraSim/Ara_Bicone6in_output.txt
      • This needs to be the input for our .job script. 
• Add lines to your newly made "Veff vs Gen" plotting code that reads in the Veff output for this ARA_actualBicone AraSim output file (whatever you've chosen to name it when you wrote it into the job script. 

 STEP 3: Move the ARA_actualBicone output file into the $Run_Name directory so it is not overwritten ever. 

• We will want to make sure this is not overwritten after the run ends, and we start a new one. 
• We will also be using it for plotting "Veff vs Generation". 

STEP 4: Correct the "wait" command for the AraSim job submissions. 

• The "wait" command, on the first iteration, will have to wait for $NPOP+1 now-- not $NPOP.

TASK 2: MAKE "Theta vs Gen" Plot.  

STEP 1: Make plotting code to plot "Theta vs Generation". 

• You can follow the way we make LRPlot.py; however, MAKE SURE YOU LABEL AXIS WITH UNITS, AND HAVE A KEY.
• This will be exactly the same as this other code, except you will need to:
  1. Read in theta from your generation DNA files
  2. Add the ARA_actualBicone constant to the code as a reference point

TASK 3: Fix "L vs Gen" and "R vs Gen" plot.  

STEP 1: Correct both plots to have labeled the units for each axis every time, and a key. This is plotted in the "LRPlot.py" code. 

STEP 2: Fix our 'R' values. I do believe we are plotting the incorrect "R". Here is what I think is happening:

• We need to double-check and see what ARA defines "R" as on their current antennas. 
  • Check with Amy?
  • What does Brian Clark's thesis say?
  • Can you find a paper that shows a picture of how they define it?
• If we are wrong, this means we need to change the "R vs Gen" part in the "LRPlot.py" code so that we:
  1. Read in R, L, and Theta for each generation.
     • See how this is already done for reading in R, and repeat it for L and theta. 
  2. Do trig to calculate the same "R" as ARA_actualBicone defines. See picture below!
  3. Store that newly calculated "R" as an array, and plot THAT.
  4. Don't forget to add the ARA_actualBicone Radius as a constant line. 

TASK 4: See PDF for other items!