configuration macros - SpiRIT-Collaboration/BeamAnalysis GitHub Wiki

Here we discuss macros used to determine parameters for configuration.

• aoqByRun.C This macro looks over a range of reconstructed runs, and is used to determine fine corrections to the TOF offset supplied in ridf_events.csv. The goal AoQ (mass to charge ratio) must be supplied, as well as a range of AoQ that the isotope of interest should be in. Currently, this is only useful for fine-tuning the TOF offset, and a preliminary guess has to found manually to start the process. For 132Sn, I've used 280 ns as the initial TOF offset.

• makeTSUM.C This macro finds reasonable values to put in the "BigRIPSPPAC" configuration file, for the four columns: txsum_min, txsum_max, tysum_min, tysum_max. For TxSum and TySum, the minimum and maximum are chosen to be within three standard deviations of the main peak. We should expect that the peak with the highest TSum value should correspond to events with a single particle. Smaller values of TSum can indicate a second particle coming through the detector. If there are multiple peaks in the TSum spectra, there could be problems fitting the correct peak. For this reason, the macro saves png figures of the TSum spectra to allow manual inspection. In the event that the main peak is not contained within the histogram, the range should be broadened. In the event of multiple peaks, the range should be tightened around the main peak. The macro outputs a csv file with the results of the fit, including a suggestion for the max and min for tx and ty sum. This TSUM gate is important to set correctly, as a modified version of "TArtCalibPPAC.cc" uses these values to fill in missing information.

• Determining IC configuration file The Ion Chamber requires 4 values to be set, and this changes for each beam. We fit the energy lost in the ion chamber using a linear fit , with ch2mev_0 the offset and ch2mev_1 the slope: ICMeVSqSum=ch2mev_0+ch2mev_1*ADCSqSum

I used ICfitting1.C to find the ADCSqSum for different isotopes, and then used LISE++ to determine the actual energy lost. From that information I manually fit the information together using LINEST in Excel.

After the energy loss has been calibrated, we have to find Z as a linear fit to a function of ICMeVSqSum and beta(the function can be found in standaloneICCalib.C).

Z=zcoef_1+zcoef_0*f(beta, ICMeVSqSum)

standaloneICCalib.C can be used to determine these parameters, and it attempts to do so automatically. The user may need to do some changes to their macro, including adjusting the range of "histICMeV" to have correct peaks.

Within the 112 Sn setting, I have observed that early runs do not match later runs, but I have not created multiple IC files for this setting. The IC configuration should work for the majority of 112 Sn runs.