• View
• Edit
• History
• Print

Experimental Code

Tuesday 15 September

Analysis of BGO data for calibration from Brandon's MBP 2:

1. Compatible with the beam setup and ntuple from the night run.
2. Comparison between noise, clearance and incrementation here
3. Visual representation of low-energy noise here

Wednesday 16 September

DWC Analysis from Brandon's MBP 2:

1. Plotting of the DWC mean averages here

MyClass.h file for afternoon ntuples from Brandon's MBP 2:

1. here

Thursday 17 September

MyClass.h file for Silicon run from Brandon's MBP 2:

1. here

Raw Theta Plot Analysis Code from Brandon's MBP 2:

1. here ... Note ... this code has an example of a disk shaped area cut for a DWC. It uses the formula √( ((x-a)/α)² + ((y-b)/β)² ) < 1.

You can download a new MyClass.h .... all the nTuples with a "y*" in the Google doc spreadsheet have a new nTuple which includes angle information. Here (MyClass.h)

 BGO fit code


 MyClass2015_17 .... latest MyClass.h .... but note the name ... rename  or change the #include as necessary

Friday 18 September

Here is a recent code with various angle stuff

 Code

Rough Analysis of Prime Theta Values with reference to BGO readings from Brandon's MBP 2:

1. File 1
2. File 2
3. File 3

How to chain several tuples together

1. Place yourself in the folder where the nTuples you want to chain together are.
2. Open a root session
   1. TChain chain("<TTree object>"); for us its tree
   2. chain.Add("<ntuple-1.root"); for us its 1442505676complete.root
   3. chain.Add("<ntuple-2.root"); for us its 1442499126complete.root
   4. chain.Add("<ntuple-3.root"); for us its 1442518129complete.root etc ... as often as necessary
   5. chain.Merge("<filename>.root"); this is your choice of output nTuple filename
   6. chain.Write();

18 Sept BGO fit code {{Attach:bgofit15.cpp|BGOfit }}

Saturday 19 September

1. Graphical cutting code from Brandon's MBP 2: here (Purely graphically cut at the moment)
2. We got the graphical cut code working ... but it was flaky ... it was hard to make the polygon ... it kept closing prematurely or not initiating. Anyway, the latest codes here with the MyClass.h. The TCutG parts are commented out now ... as they weren't working ... but there are still two loops over filling with the second loop having cuts on the incident angle space and also the target shape projected onto the DWC1. For Tim ... also the Makefile we are using.
3. Here is the last picture we got from the analysis at midnight ... maybe there is something in the channeling regime ... so we decide to collect more stats here and don't change anything ....

Wednesday 23 September

1. Trial Code
   1. Use with chained tuple to maximise the data
   2. Makes a cut on photon energy, and on DWC0 (pointing to the core of the BGO) and on the extrapolated target position (very tight), then displays the incident angle space to make a cut there with the graphical editor, and then puts a new energy spectrum up with these cuts.
   3. test.cpp