• View
• Edit
• History
• Print

Experimental Code (after)

General Strategy

Make sure you have a version of MyClass.h that is appropriate for the tuple you are analysing

1. If necessary .. .see section 12.21.1 Creating a Class with MakeClass
2. Use a folder for each data set, even a sub-folder for each code you are trying.

Make an elliptical cut on DWC0 where you look how the events that could be radiation in the BGO are mapped on DWC0.

1. Below is an example of quite tight cuts, leading to the core of the BGO. (reversed cut shown for clarity)

Look at DWC2 .... (depending on the B-field and how the DWC2 has been moved) .... you will see the position of the (deflected) positrons.

1. Note that positrons that have radiated photons have less energy and will be deflected more.

You can project from the hits in DWC0 and DWC1 to DWC2, and plot their difference.

1. You can see in the code how extrapolation is done by taking the line from the hits in DWC0 and DWC1 to the DWC2.
2. Below is an example of the difference of the projected hit on DWC2 from the actual hit.

1. The radiation results in a lower momentum positron leading to the positrons of interest missing the DWC2. These positrons are on the right of the DWC2_projected-actual. A cut can be made for these events

Now look at the target position for the photon energy region region of interest, with also the DWC0 and DWC2 cuts.

1. You can see in the code how position of the position at the target position is extracted by taking the line from the DWC0 and DWC1 at the target plane.
2. Try to position the elliptical cut where you think the target is (except for Si, which was bigger then the veto scintillator hole).

Now display the angle space variables and try to see an axis or plane .... make a graphical cut on these.

Now with all cuts, display again the full energy spectrum and see if you see any photon peaks.

An example for the Si target ... run 9552

1. Here we were at E_e+=6 GeV, expecting a 100 channeling radiation peak at 45 MeV, so the cut in E_e+=(30-60) MeV
2. The sample code is Analyse_9552.cpp with the MyClass_9552.h
3. ***** More recent sample code is Analyse_und_tilt.cpp with the MyClass_und_tilt.h
4. ***** Even More recent sample code with cut variables defined is Analyse_und_tilt2.cpp with the previous MyClass
5. You can run the code with cuts on/off or invert them to visualise the various cuts at various stages.
6. The stats are just very low, .. but there could be a remnant thats a photon peak.

Undulator_05Tilt Analysis

1. Analysis code here
2. Report is on the Google Drive under 'Run Analysis'