Run 68 Processing and Analysis

This is meant to be a landing page to aid in coordinating processing and analysis tasks with R68 data.

Analysis planning document

Analysis task list table

Series summary spreadsheet: dataserieslist.xlsx

Notes:

• first-pass SiHV+NaI RQ data: TTree zip1 contains the SiHV channels, mapped to their proper Side 1 locations. TTree zip4 contains the two NaI waveforms in PA and PB. OF processing is not done on the NaI data so if zip1 is friended to zip4, OF-RQs will refer to SiHV without ambiguity. But RTFTwalk RQs need a specifier for detector, e.g. zip1.PAWKr50 versus zip4.PAWKRr50.
• Second Pass, n125V changes: Use empirical templates for HV detector, turn on OptimalFilterPhonon1X2 for PT, turn off PS1*, PS2* calculations, adjust P_PEAK_WINDOW_MIN/MAX for DETECTOR 1(4) to 75e-6/25e-6(44.8e-6/358.4) to reduce noise wall(ignore first 8 bins), turn on P_OF(1×2)_DELAY_INTERPOLATE, change P_RTFT_PULLBACK from 150 to 1 for both detectors, update P_RELATIVE_CALIBRATION to these values, set P_BIAS DETECTOR 4 values negative for inverted NaI pulses, change P_BASELINE_MIN from 1 to 9 to avoid 8-bin glitch
  • 2-Template OF is enabled for PT with the standard PTOF template as the “primary” template and a linear template as “residual” as explained for v9 templates here and here. The relevant RQs are the primary amplitude “PTOF1X2Pamps” which should be interpreted as the pulse amplitude and the residual amplitude “PTOF1X2Ramps”, which is proportional to the baseline slope.
  • A second parallel calculation of the single template OF which results in symmetric noise distributions as discussed here is also present. The relevant rqs are “PTOFSamps”, “PTOFSchisq”, etc.

Run 68 Templates

The following datasets have been defined:

*Note that from the beginning of the run until 09/26, the NaI PMT bias was mistakenly set to 1.2 kV, even though the run diary records it as being set to 1.1 kV. This table lists the correct biases.

You can use these files to load data into a ROOT session using something like this:

//Load the series given a list of series names
TString seriesName
std::ifstream infile("/home/mast/cdms/analysis/run068/datasets/p40_bkg_2.txt"); 
int Nser_Off=0, Nser_Off_max=NserMax
while( (infile >> seriesName) && (Nser_Off < Nser_Off_max{
  z_Off[Nser_Off]->Add(Form("%s%s/umn_*.root",RQDir.Data(),seriesName.Data()));
  e_Off[Nser_Off]->Add(Form("%s%s/umn_*.root",RQDir.Data(),seriesName.Data()));
  Nser_Off++;
}
infile.close();

These can vary for different data types and biases. Links to cuts and calibrations developed for specific datasets can be found below.

here (Defunct)

TCut crand = "EventCategory";
TCut cBShi="PFbs-PFbsmean>1500";
TCut cFstd="PFstd<2.5e-7";
TCut cofdel="(PTOFdelay<23e-6 && PTOFdelay>-74e-6)";
TCut cofwktiming = "(PTOFdelay*1e6<PTWKr50*1e6-390 && PTOFdelay*1e6>PTWKr50*1e6-470)";
TCut cglitch = "1e6*(PTWKf20-PTWKr100)<12";
TCut cspikey = "1e6*(PTWKf20-PTWKr100)<35";
TCut cchit = "PTOFchisq<(PTOFamps>6e-6)*((PTOFamps-6e-6)/3e-8)**2+PTOFamps/2e-9+5e3";

z->SetAlias("PTINTall_PTdbscorr","PTINTall-(PTdbs<0)*0.5*(PTdbs*4096./PAnorm)");
z->SetAlias("PTINTall_PTdbscorr_bscorr","PTINTall_PTdbscorr*(1.+0.18*(PFbs-PFbsmean)/2e3)");

TCut cofintHi_bkg="PTINTall_PTdbscorr_bscorr> (PTOFamps>1.09042e-05)*6.21097e+06*(PTOFamps-1.09042e-05)**2+170.476*PTOFamps+1e-4";
TCut cofintHi_Na="PTINTall_PTdbscorr_bscorr> (PTOFamps>7.96642e-06)*1.22109e+07*(PTOFamps-7.96642e-06)**2+169.406*PTOFamps+1e-4";
TCut cofintHi_PuBe="PTINTall_PTdbscorr_bscorr> (PTOFamps>9.59954e-06)*1.22032e+07*(PTOFamps-9.59954e-06)**2+171.253*PTOFamps+1e-4";
TCut cnopulsepus[3]={cofdel&&!cofintHi_bkg&&cofwktiming, cofdel&&!cofintHi_Na&&cofwktiming, cofdel&&!cofintHi_PuBe&&cofwktiming};

z->SetAlias("N100","nPTWKmax_N_N100");
z->SetAlias("n50_200_N100","nPTWKmax_50nA200nA_N100");
TCut cLERate_bkg = "N100==100&&n50_200_N100<101";
TCut cLERate_Na = "N100==100&&n50_200_N100<5";
TCut cLERate_PuBe = "N100==100&&n50_200_N100<4";
TCut cLERate[3]={cLERate_bkg, cLERate_Na, cLERate_PuBe};

TCut cgood[3]={!crand&&cFstd&&!cBShi&&cnopulsepus[0]&&!cspikey&&cchit&&cLERate[0], !crand&&cFstd&&!cBShi&&cnopulsepus[1]&&!cspikey&&cchit&&cLERate[1], !crand&&cFstd&&!cBShi&&cnopulsepus[2]&&!cspikey&&cchit&&cLERate[2]};

Full cut and calibration development can be found in the following links and the notes referenced therein.

• Reprocessed -125V Cuts
• Am Calibration