View source for Spot++

<br>
[[Main_Page|<B>PREX Main</B>]]<< [[Information for Shift Takers]] << [[Shift Leader]] << [[HOW TOs]]

== Description ==

Performing a spot++ allows us to "image" the target or raster size by plotting the hits as a function of raster x and y current
* <font color="blue"> CASE 1</font> If you run with T1 or T2 trigger and target in place, spot gives a view of the target density
* <font color="blue"> CASE 2</font> If you run with pulser trigger T8 and high current, spot gives a (misleading but consistent) view of the rastered spot size using the BPMs

The steps involved in doing this measurement are:

# Take the spot++ run
# Analyze the run
# Make the plots and post to the logbook

== How To Take a High Current spot++ run ==
# This will give only raster size. Monitor BPM positions for target locks.
# the high current spot run can be taken at any current above 1 uA or so.  (Below this, the bpms won't work). This does not interfere with parity production.
# Set the trigger prescale factor for PS8 =1 for RHRS and PS8 = 100 for LHRS, and all other triggers =0.
# Start the counting DAQs while beam is on.
# End the run after 50k events are collected. Run spot++
# If you log the results (and you should) make a follow up entry to label this run as a "high-current spot" run.
# The 2D spot plot should look uniform.

== How To Take a Low Current spot++ run ==

Detailed instructions from [https://logbooks.jlab.org/entry/3711297 halog 3711297]<br>
To center the beam on the Carbon hole center [https://prex.jlab.org/wiki/index.php/Low_current_running_beam_setup#Carbon_hole_target_.28spot.2B.2B.29]
# Stop Parity DAQ, turn off beam modulation and Aq feedback.
#* Ask MCC for 0.5uA current, raster requirement and target lock for a particular target (see white board).
# Set up the detectors for taking the spot data
#* Turn S3 and/or S0 scintillator HVs on for both arms. These are connected to channels L0.2 (S0 tab) in the HRS HV control panels.
#* Turn on other scintillator paddles and/or VDCs if requested
#* Observe the S3 and/or S0 trigger rates in the scaler GUI. Aim for less than 100 kHz.
#* If S3 and/or S0 rates are too high, ask for a lower current.
# While the shift leader was performing the above, the target operator should prepare the [[CODA#Starting Counting Mode DAQ | HRS DAQ]]
#* Set prescale factor PS1=1,PS2=1 and PS6=1 and all other prescale factor should be zero.
#* Take a test run to make sure things are working while the beam is being set up
#* Take at least 50k events for a quick spot, but for a target uniformity spot ~1M is better
#* End CODA runs
# Back out of the spot run
#* Turn S3 and S0 HVs off
# Unless otherwise instructed, go back to production running
#* Start a parity run

== How to analyse the run ==

# Log into the a-onl account
#* <font face="Courier">ssh –X a-onl@aonl2</font># Then type either:
#* <font face="Courier"> spot_L ''run#'' </font> for the LHRS (Recommended) <p> <font face="Courier"> spot_R ''run#'' </font> for the RHRS </p>
#*If no run # is given, analyzes last run taken
# Plots should pop up after a few seconds
#* The most important plot is on the first page, titled "beam on target y vs x", see, for example [https://logbooks.jlab.org/system/temporary/FADC_prex_21079.pdf RHRS run 21079]
# When prompted, type <font face="Courier"> y </font> to upload the plots to the halog or <font face="Courier"> n </font> to skip this step

[https://hallaweb.jlab.org/equipment/daq/spot.html More Info]

== More information ==

==== Making plots for long Runs ====

Most of the time we just want the normal length of spot++ run, which is 50k events.  But occasionally we want to take a higher statistics run, and the spot default is to only analyse 50k events.  However, a command like:

<font face="Courier"> spot_L ''run#'' ''firstEvent'' ''lastEvent''</font>

will analyze the events in the specified range. 

This will output the .pdf to the logbook in the same way as above.

==== raster_check.C ====

An additional way to make plots is to use the raster_check script.

* go to the <font face="Courier"> /adaqfs/home/a-onl/rastersize/</font> directory
* open root
** <code> .L raster_check.C</code>
** <code> raster_check(<1 or 0>,run#)</code> (1 for LHRS, 0 for RHRS)
*** This version of the function creates a single plot of the map of the raster at the target. This plot as well as the rest of the macro uses the Raster2.rawcur
*** If this is done on the carbon hole target, you should see the hole on the plot. Make a rough visual estimate of what the bounds of the hole are on both the left and y axis
** <code> raster_check(<1 or 0>,run #,xmin,xmax,ymin,ymax) </code> (1 for LHRS, 0 for RHRS)
*** <code>xmin</code>, <code>xmax</code>, <code>ymin</code>, <code>ymax</code> are the bounds on the carbon hole target from the 2D plot above
*** This function creates two plots, profile plots of the 2D raster map. Each plot has three different fiducial cuts applied in order to estimate the hole position better.
*** By zooming in on the x-axis of each plot and using the Event Statusbar, estimate the point on the x-axis of each plot where all three histograms overlap each other, on both sides of the hole.
*** To estimate raster size from the profile plots calculate 2*(max bin - min bin)/(max histogram overlap point - min histogram overlap point). Do this independently for both x and y.
** <code> raster_check(<1 or 0>,run#,0,100e3,0,100e3)</code> (1 for LHRS, 0 for RHRS)
*** This function creates the profile plots like above but applies no fiducial cut to the histograms. This is best when looking at spot runs on the production target.
* Post these plots to the halog if desired

==== More information about spot++ ====

* Hint: Type "spot help" to learn more options for spot++
* spot++ output files go to:
 /chafs1/work1/spot/
* For experts - macros, source codes, DB files are at:
 /adaqfs/home/a-onl/rastersize/ 

* You can get "strange printouts" of errors if you only type spot++ and the last run was for R-HRS because spot++ is aliased to spot_L

=== What do the plots mean? ===
Need to put examples of output when we have them

'''BPMA y vs x:''' Shows the beam location at BPM A

'''BPMA x:''' Shows the number of events at each x value in BPM A

'''BPMA y:''' Shows the number of events at each y value in BPM A

'''BPMB y vs x, BPMB x, BPMB y:''' Same as above but for BPM B

'''Beam at target y vs x, …x, …y:''' 	Same as above for at the target

'''BPMA Antenna 1, 2, 3, 4:''' Signal at the four wires of BPM A, one for X+, X-, Y+, Y-, respectively

'''BPMB Antenna 1, 2, 3, 4:''' Same as above but for BPM B

=== Cavity BCMs === 
* The cavities are used in <0.5 uA running territory (for commissioning a lead target with spot ++ check)
* The cavity calibration is not quite reliable.  Positions are best found by running at 1uA to set the position using the stripline bpms 1H04A and 1H04E, and then record the corresponding cavity positions.
* Bob's EPICs variable printing to terminal script: 
 ssh apar@adaq
 ~/scripts/showCav.C