• To Do
  • Since pythia 6 wont work on the umd cluster or this macbook i need to put it on a computer in the office
  • Do decied what large event runs for pythia 6 and pythia 8 to do
  • learn to use plot digitizer for article plots
  • add the data points gotten from plot difitizer to completed code for plotting article and flux.

• to move files scp file dagherni@pa-pub.umd.edu:folder

• I want to run the corsika file through hawc sim in order to see how far i can throw the core and still trigger on the array

Outriggers

• Need to make plots of radius vs #triggered showeres/#totalshowers (Where triggered showeres are showeres that hit ~30pmts on the HAWC300 array) for different energys. In order to find the max radius our experiment can detect for different energy.

• Need to see the most optimal method of arranging the outriggers, Going to use the sunflower sprial to start but may mask tanks off to create a denser array near HAWC300 and lesser dense array farther away from it.

• use files that dirk sent : NoBadChannels.xml HAWC300_v1.0.xml and plot-detector.cc to create a file that is titled what ever you choose when you compiled plot detector.

• Now i am trying to use a file titled create_spiral_outriggers.py but i am getting the following thing appears

Breaking Homer

Outriggers

• Need to make plots of radius vs #triggered showeres/#totalshowers (Where triggered showeres are showeres that hit ~30pmts on the HAWC300 array) for different energys. In order to find the max radius our experiment can detect for different energy.
• Need to see the most optimal method of arranging the outriggers, Going to use the sunflower sprial to start but may mask tanks off to create a denser array near HAWC300 and lesser dense array farther away from it.
• use files that dirk sent : NoBadChannels.xml HAWC300_v1.0.xml and plot-detector.cc to create a file that is titled what ever you choose when you compiled plot detector.
• Now i am trying to use a file titled create_spiral_outriggers.py but i am getting the following thing appears

  • nick@nick-Satellite-P755:~/workspace/dagherni/outriggers$ python create_spiral_outriggers.py
    This script requires: numpy, matplotlib, argparse, xml and StringIO

  • i beleive i have numpy and matplotlib but i need to check for the others

• to fill histogram ande declare it

  Hist ZpT( "Z0 pT spectrum", 100, 0., 100.); 

• ZpT.fill( 22.7, 1.)

• • creat a new decay channel to go to 100% of the time for the z' to go to either bbbar,ttbat,ww,zz
  • Initialize collidpor and run 10000 events and colect energy of photons after every event

• fos wiki link to a student who worked on z' in pythia 8 https://hep.pa.msu.edu/wiki/bin/view/ZPrime/E6Generation

To do with pat

• go over thought proccess of replicating the dark matter collision
• go over the pythia 6 code in order to connect the thought proccess and code
• what parts of the code are specific to dark matter

• the key points of the file are:
  • a positron and electron annihalate into a z, z' and a gamma. We only want the proccess that goes to z'
  • turn off inital-state radiation,photonic part of pdf for e+ e- interactions.
  • Allow muons,pions,kaons, and k_long to decay. (try and add nutron decay to pythia 8)
  • in muon decay:mu -> nu_ebar + e- + nu_mu + gamma is 1.4% and mu -> nu_ebar + e- + nu_mu is 98.6%. Need to set these decays to there per centage and creat a decay channel to store the rare decay.
  • turn off all pre-existing Z' decay channels
  • give z' parameters of mass=com, no mass width, no broadening, and no lifetime
  • 
    

• pythia 8 html manual http://home.thep.lu.se/~torbjorn/pythia81html/Welcome.html

Pythia 6

code structure

• The code calls for 100000 events, and a mass of 1 TeV
• "Basically, this means that there are three ways by which you c

  ommunicate with the

  programs. First, by setting common-block variables, you sp

  ecify the details of how the pro-

  grams should perform specific tasks, e.g. which subprocesse

  s should be generated, which

  particle masses should be assumed, which coupling constant

  s used, which fragmentation

  scenarios, and so on with hundreds of options and parameters

  . Second, by calling subrou-

  tines you tell the programs to generate events according to t

  he rules established above.

  Normally there are few subroutine arguments, and those are u

  sually related to details of

  the physical situation, such as what c.m. energy to assume fo

  r events. Third, you can

  either look at the common block

  PYJETS

  to extract information on the generated event, or

  you can call on various functions and subroutines to analyse

  the event further for you" -Pythia 6.4 manual
• Appendix B is an index to all subprogrograms and common-block variables

• Use "xcdf-utility info" and "xcdf-utility dump" in a triggered file and a reco file.

• Week of june 8th
  • Do aerie hw 1
  • how does pythia 8 work
  • how do you make plot the dark matter in pythia 8

• ssh -Y yourname@pa-pub.umd.edu and then into ssh yourname@sequoia

Tutorial

• RootClass2014.pdf

Weekly goals

• Week of June 1st
  • have pats plots of Pythia 6 working
  • how pythia8 and root work together
  • read "Reliability of Monte Carlo event generators for gamma-ray dark matter searches"
  • get mary cluster and start aerie homework

Getting Acquainted

• "Sensitivity of the high altitude water cherenkov detector to sources of multi-TeV gamma rays"
  • Sensitivity_HAWC.pdf
• "Notes on Calculating Dark Matter Fluxes"
  • HAWC_Harding_note.pdf
• "WIMP Annihilation Limits from Segue 1 with Early HAWC-30 Data"
  • HAWC_Pretz_note.pdf
• "Prepared for submission to JHEP

  Reliability of Monte Carlo event generators for

  gamma{ray dark matter searches"

  • MCdiscrepancyDM.pdf
• "VERITAS Deep Observations of the Dwarf Spheroidal Galaxy

  Segue 1"

  • Veritas_Segue1.pdf

Aerie

Homework

• http://private.hawc-observatory.org/wiki/index.php/Analysis_Workshop_2015_06_15

Root

Browser

Installing root

Plotting functions

Plotting Histograms

pointers

Analyze

Plotting N-tuples/trees

• • • <del>scrolling the graph from left to right is a little glitchy</del>

Analyze

• setup:define the histogram
• loop section: fill the histogram with values
• wrapupsection: draw the histogram

Readings

• Sensitivity of the high altitude water cherenkov detector to sources of multi-TeV gamma rays
  • https://hep.pa.msu.edu/wiki/pub/HAWC/JenniferRanta/Sensitivity_HAWC.pdf

• • • if you right click the bottom graph it dissapears

• • • the graphs on the right are getting cut off the screen unless you make the page large
    • right clicking the graph changes the the y axis
    • zooming in changes the y axis but when you zoom back out the axis doesnt change back
  • VAC monitoring tab
    • the hovering information bubble is not showing the point the cursor is pointing at, it is off by inches
  • UPS Monitoring tab
    • nothing that i havn't mentioned earlier regarding the control of the plot
  • Reconstruction Monitoring tab
    • Can not access, still waiting for hawc wiki account
  • Residuals monitoring tab
    • Can not access, still waiting for hawc wiki account
  • Webcams tab
    • Can not access, still waiting for hawc wiki account

• • • when hovering over a tank it displays the last time updated but not the range in time the graph displays. maybe add the range and make the graph bigger if your hovering over, it is hard to read unless i click it to go to the expanded page

• • GTC Status tab
    • scrolling the graph from left to right is a little glitchy
    • it seems like the y range is too high, all the points look like a line at the bottom, so maybe a zoom option or an ability to adjust the range
    • being able to select the range in time that is displayed may be helpfull as well
  • EMS Monitoring tab
    • when hovering over the graphs the bubble of information doesnt always display what the crusor is on, sometimes it is off by ~inch
    • I tried to scroll left and right like i did for the GTC tab but it actually zoomed me in
    • If you click plot with out a start date, an end state or either then no plots appear
    • if you make the start and end date the same then a plot still appears but there is nothing on the plot
    • I think it would be useful to seperate the zoom and scroll so you could zoom in multiple times and still be able to scroll left and right
  • Lighting Monitoring tab
    • The start date and end date have the same problems as before.
    • im trying to drag the image to my desktop to save it but i keep zooming/scrolling on the graph
  • Computer Status tab
    • MY comments for this page are the same for the previous tabs regarding the start/end date and the scroll/zoom feature. Also still trouble trying to drag it instead of zooming/scrolling
    • For the bottom graph i clicked on a colored square in the key and instead of highlighting that color in the graph it removed that color.
  • Network Monitoring
    • same comments on zoom/scroll and draging to save
    • while i was clicking around on the bottom 2 graphs, i was able to drag each of them too the side and delete them from the page. I was only able to do it twice and i cant figure out how i did it haha.
  • HV monitoring
    • Draging the image off the screen does not save or do anything.
  • LV monitoring
    • the graphs on the right are getting cut off the screen so you cant see them all

Browser

• TBrowser b to open a browser
• TFile myfile("filename")

• TH1D h1("rootname",title,#ofbins,binmin,binmax)

pointers

• TH1D *hist1 = new TH1D ("h1","a histogram",100,-3,3)
• this creates a pointer

Plotting N-tuples/trees

• tree1->Scan() runs like 25 events
• tree1->Print() it displays the names of the variables and the size of the tree
• tree1->Draw("variableyouwant") makes a histogram for a variable in the tree

Breaking Homer

• HW scalers tab
  • If you click anywhere on the color changing key on the right or one of the 4 example circles labeled 0,20,40,60 then you need to refresh the page in order to click on anything else.
  • If you click on the tank and not the dots in the tank then you need to refresh in order to click on anything else
  • when hovering over a tank it displays the last time updated but not the range in time the graph displays. maybe add the range and make the graph bigger if your hovering over
  • Once you click on a dot and the graph displays it is kind of small so a zoom option would be usefull

Plotting functions

• TF2 f1("func1",thefunction,xmin,xmax,ymin,ymax)
• then next line type f2.Draw() use "surf1" for a surface plot

Plotting Histograms

• TH1D h1("hist1",title,#ofbins,binmin,binmax)
  • This creats an empty histogram, we must fill it
  • h1 is the c++ name and hist1 is the root name
  • use h1.FillRandom("gauss",10000) to fill it with 10000 points of data fraom the functions "guass"
• TF1 myfunc("myfunc",
• hist1 is the root name of this histogram
• h1.GetBinWidth(0) to get the bin width
• TF1 myfunc("myfunc","gaus",0,3)
  • myfunc.SetParameters(10.,1.0,0.5) to set P0,P1,P2 the parameters of the gaussian function "Constant" "Mean and "sigma"
  • *

Root

• Installing root

• • went to the root website and went to downloads tab and clicked patches v5-34-00. Downloaded the source link * after downloading the source link and extracting the files run "./configure" and install any needed required packages, use the website for assistance
  • After you run configure after installing necessary packages. run "make" in the terminal and then run "make install"

Pythia 8

• In order to run a file type( from the worksheet file, can be found on the pythia website)
  
  • make mainNN
  • ./mainNN.exe > mainNN.out

Running and compiling fortran code

• For help with fortran code http://wwwcdf.pd.infn.it/localdoc/f77_sun.pdf
• type gfortran -static filetocompile.f pythia.o
  • the file to run is saved as a.out (./a.out to open)

November 18

• After having trouble dool booting linux i have decieded to use a virtual macheine software to launch linux
• I have downloaded and am able to succesfull access linux now
• 
  

Dual booting Linux with windows 8

• In order to use pythia I am going to dual boot ubuntu on to my computer. I have followed the instructions from this website
  • http://www.everydaylinuxuser.com/2014/05/install-ubuntu-1404-alongside-windows.html
• In order to do this you need 2 usbs, (I used 16 gb to make sure i had enough space) and an external hard drive (you really just need anything that will back up your personal files)
• The main steps that I took were
  • create a recovery drive with one of the usbs
  • back up all personal files onto an external hard drive
  • Install ubuntu from their website
  • Install the driver that will put ubuntu on too the 2nd usb
  • turn off fast boot option
  • turn off secure boot, for help see https://neosmart.net/wiki/disabling-secure-boot/

• Read up on pythia