Difference between revisions of "Cosmic Data Analysis"
(→Root Analysis Script) |
(→Root Analysis Script) |
||
Line 48: | Line 48: | ||
Access directory "analyzer" which contains directories "replay" and "rootfiles". Directory "replay" contains script for writing raw data files into root files which are saved in directory "rootfiles". In directory "rootfiles" there is the plot.C script. | Access directory "analyzer" which contains directories "replay" and "rootfiles". Directory "replay" contains script for writing raw data files into root files which are saved in directory "rootfiles". In directory "rootfiles" there is the plot.C script. | ||
− | |||
− | |||
To run plot.C type commands: | To run plot.C type commands: |
Revision as of 15:08, 3 July 2017
Contents
Motivation for Cosmic Ray Analysis
Cosmic rays provide a good proxy for high energy charged particles that will be found in Hall A during the time of experimentation. The robust nature of cosmic rays allow for vertical track testing at any point of the day/year, which makes them ideal for continuous testing/debugging methodology used in CDet.
DAQ Computer Access
The current computer setup for the DAQ is sbs2 which is located in the test lab, under the user adaq. it is recommended when setting up coda, to do so at the test lab setup, as remote ssh connections can be terminated periodically from the jlab host server.
To boot CODA (Data acquisition system)
1. Set the High voltage crates to On and Local
2. Login to sbs2 under user adaq
3 On CODA desktop
- open terminal and type msqld - on a separate terminal, type starcoda - wait for the main CODA window to appear before changing desktops
4. Change to desktop ROC
- open terminal, type telnet cdetts2 2006 - once the telnet is connected, it will tell you there is an escape character (^[) just press enter ONCE - an "->" prompt will appear, type reboot - The terminal will say "callbackControl: Do not understand the command: Session/control/setSession"
5. Change Desktop to HV (high Voltage)
- open terminal, ssh into highv@hvsrv3 - change directory into slowc - type './hvs TEST - refer to High Voltage Control to operate.
Hall A Analyzer
After a completed CODA run, the raw data must be converted into a usable root file.
On the computer sbs2 under user adaq, in directory analyzer/replay
1. Type analyzer
ASCII text will appear that resembles the root startup.
2. To analyze a file, type .x setup.C
3. When prompted, type the event number that is displayed on the CODA window
4. The program will then ask you to type the amount of events you wish to analyze (for all events type -1)
The root file will be saved in the directory analyzer/replay/rootfiles
Root Analysis Script
Access directory "analyzer" which contains directories "replay" and "rootfiles". Directory "replay" contains script for writing raw data files into root files which are saved in directory "rootfiles". In directory "rootfiles" there is the plot.C script.
To run plot.C type commands:
cd rootfiles
root -l
.x plot.C(####)
The #### is the run number to analyse. All events analyzed using the hall a analyser and the setup.C script will load for analysis.
The script plot.C contains many pointer methods with descriptions and directions for running as follows. All pointer methods must be run after executing the plot.C script for the specified run number.
Plot ADC Fit
To run plot adc fit type command:
plot_adc_fit(#)
The # is the number of the scintillating bar (1-14) to analyse.
The pointer method plot_adc_fit plots a 4 by 4 canvas of histograms for each of the fourteen scintillating paddles within the selected scintillating bars. ADC cuts are placed on neighbouring paddles at less than value 10 amplitude to ensure that the data used for the selected paddle has a near perfectly vertical track through the paddle.
The method also fits the data after cuts are applied using a Gaussian fit. This allows for further analysis of parameters such as means, standard deviations, amplitudes, and total event values.
Plot ADC
To run plot adc type command:
plot_adc(#)
The # is the number of the scintillating bar (1-14) to analyse.
The pointer method plot_adc plots a 4 by 4 canvas of histograms for each of the fourteen scintillating paddles within the selected scintillating bars. Each histogram plots two graphs on top of each other showing the ADC spectrum with just TDC cuts and then also with ADC cuts used for the selected paddle.
Plot TDC
To run plot tdc type command:
plot_tdc(#)
The # is the number of the scintillating bar (1-14) to analyse.
The pointer method plot_tdc plots a 4 by 4 canvas of histograms for each of the fourteen scintillating paddles within the selected scintillating bars. The histograms plot the TDC spectrum for each of the fourteen scintillating paddles in the selected bar.
Plot TDC and ADC
To run plot tdc adc type command:
plot_tdc_adc(#)
The # is the number of the scintillating bar (1-14) to analyse.
The pointer method plot_tdc_adc plots a 4 by 4 canvas of histograms for each of the fourteen scintillating paddles within the selected scintillating bars. The histograms plot the TDC spectrum and ADC spectrum for each of the fourteen scintillating paddles in the selected bar.
Plot Occupancy
To run plot occupancy type command:
plot_occupancy()
The pointer method plot_occupancy plots a canvas containing a three graphs. The main graph plots the occupancy for ever channel for every scintillating paddle in the half detector set-up. This allows the user to see if certain channels are providing sufficient data. The lower left graph plots the multiplicity. The lower right graph plots a heat map for all of channels in the half module set-up.
Plot Occupancy Single
To run plot occupancy single type command:
plot_occupancy_single(#)
The # is the number of the scintillating bar (1-14) to analyse.
The pointer method plot_occupancy_single plots a canvas containing graphs showing the occupancy for ever channel for every scintillating paddle in the half detector set-up. This allows the user to see if certain channels in and individual scintillating bar are providing sufficient data.