Table of Contents
Run 2
Summary
- Upgrades from last run: added a calibrated RuOx thermometer on the detector housing and a thermal link between the 600mK floor and still.
- Detector housing base T: ~25mK. This is from extrapolation. The calibration lower limit is 43.8mK.
- Installed cdmsbat and analyzed background data (with help from Lauren, Bruno, and Jeter). Results shown below.
Setup
| Detector | G9 (former T1Z5) |
|---|---|
| Tower | T1 from Soudan |
| Face | #5 |
| Side coax | (the shortest one) |
| SQUET | ZC15 |
| Stripline | #20 |
| Readout | DCRC v1 |
A thermal link between the 600mK tower floor and the fridge still is added. We also mounted a calibrated RuOx thermometer from Lakeshore on the detector housing.
Fridge Performance
The base temperature of the thermometer on the detector housing is around 25mK this time. In the last run, the base temperature of the thermometer on the mixing chamber was about 38mK. The base temperature is also relatively stable.
Tc Measurement
Send 1mVpp 100Hz triangle signal, apply different power to the mixing chamber, and watch at which resistance (temperature) the triangle wave disappear or reappear. The resistances were read when the temperatures were changing fast, there were some lag effect.
| MC power (uW) | 50 | 0 | 20 | 0 | 10 | 0 | 10 | 6 | 5 |
|---|---|---|---|---|---|---|---|---|---|
| Temperature change | ↑ | ↓ | ↑ | ↓ | ↑ | ↓ | ↑ | ↑ | ↓ |
| phA (kΩ) | missed | 26.96 | 31.6 | 26.9 | 29.6 | 27.67 | 28.8 | 29.4 | 28 |
| phA (mK) | missed | 81.3 | 73.8 | 81.4 | 76.7 | 80.0 | 78.0 | 77.0 | 79.4 |
| phD (kΩ) | missed | 26.75 | 30.8 | 26.4 | 28.8 | 27.6 | 28.1 | 28.5 | 27.7 |
| phD (mK) | missed | 81.7 | 74.9 | 82.4 | 78.0 | 80.1 | 79.2 | 78.5 | 79.9 |
- G9 after Implantation (Run 14 results)
- SQUET # 8, Sidecoax ZB5-1, Stripline #25, ZIP Board #5, RTF Board #1.7
| Sensor | Transition | R Superconducting | Rsc (75 Run 270) | R Normal at 150 mK | R Normal from IbIs | Rn (75 Run 270) |
|---|---|---|---|---|---|---|
| A | 79 mK-85 mK | 12.1 mOhms | 9.1 mOhms | 1.48 Ohms | 1.50 Ohms | 1.55 Ohms |
| B | 83 mK-87 mK | 10.9 mOhms | 11.2 mOhms | 2.17 Ohms | 2.80 Ohms | 1.85 Ohms |
| C | 81 mK-85 mK | 16.5 mOhms | 13.0 mOhms | 2.17 Ohms | 2.12 Ohms | 2.0 Ohms |
| D | 82 mK-88 mK | 21.1 mOhms | 9.5 mOhms | 1.53 Ohms | 1.51 Ohms | 1.66 Ohms |
We obtained consistent Tc values with that from CWRU at June 5, 2001.
Ic vs T Measurement
| Ic vs T | ||||||
|---|---|---|---|---|---|---|
| HP3325A (mVpp) | Ic(uA) | Tower T Up (mK) | Tower T Down (mK) | |||
| A | D | A | D | A | D | |
| 2 | 0.63 | 0.63 | 75.7 | 76.6 | 78.8 | 81 |
| 4 | 2.72 | 2.74 | 73.6 | 75.7 | 78.9 | 79.6 |
| 6 | 4.04 | 4.12 | 75.4 | 77.1 | 78.7 | 79.3 |
| 8 | 5.4 | 5.45 | 73.7 | 77.2 | 78.8 | 79.8 |
| 10 | 6.35 | 6.45 | 74.8 | 77.1 | 78.8 | 79.5 |
| 20 | 12.5 | 12.7 | 73.8 | 75.2 | 78.1 | 78.5 |
| 30 | 20.6 | 20.8 | 73.7 | 74 | 77.8 | 78 |
| 40 | 27.6 | 28 | 77.3 | 77.6 | ||
| 50 | 35 | 35.4 | 72.5 | 72.7 | 77.6 | 78.1 |
| 60 | 41.6 | 41.8 | 74 | 74 | 77.3 | 78.1 |
| 80 | 55.5 | 56.5 | 72.6 | 72.7 | 76.9 | 77 |
| 120 | 62 | 63.5 | 73.6 | 73.7 | 76 | 76.3 |
| 140 | 75 | 75 | 71.6 | 72.6 | ||
| 160 | 81 | 81 | 70 | 71.4 | ||
| 180 | 96 | 95 | 73 | 73 | ||
| 200 | 108 | 108 | 67 | 69 | ||
| 240 | 138 | 141 | 54 | 61 | ||
| 280 | 161 | 168 | 32.9 | 36.9 | ||
IbIs
| Senser | RSC (mΩ) | RN (Ω) |
|---|---|---|
| A | 4.0 | 1.32 |
| B | 2.2 | 1.66 |
| C | 3.1 | 1.87 |
| D | 1.9 | 1.48 |
Noise
The noise spectrum from the DAQ software and the analysis package cdmsbat are plotted. Note that the psd from the cdmsbat are normalized to roughly match that from the DAQ software.
DCRC and DAQ vi
DCRC
- DCRC v1 is used. Trigger thresh hold can not be adjusted with the resolution as needed.
- Charge: ADS1602, 2.5MHz 16-bit digitizer
- Phonon: ADS1601, 1.25MHz 16-bit digitizer
- Charge gain after Vf: 10.
DAC vi
- The glitches on the pulses caused by the lookback time are removed.
- Trace length 1.6ms, lookback time 400us
- Charge trace length 3998, pre-pulse 998, post pulse 3000.
- Phonon trace length 1998, pre-pulse 498, post pulse 1500.
- Flashing during data taking feature is added to the vi.
Features good to add
- Stop the run after a given amount time/number of events.
- Save hardware settings in a file to be used in data processing
cdmsbat
cdmsbat is installed and configured to run successfully with the help of Bruno, Lauren, and Jeter. Some information can be found here
Data Summary
| Data Run # | Run nickname | Data Start | Data End | Comment |
|---|---|---|---|---|
| 201 | f3md3m | 07100312_1454_F0001 | 07100312_2110_F0008 | 100 Files. DCRC: Preset 2; LED1: I=400uA, ON=100us, OFF=2500us; Data taking: 3min, Flash: 3min |
| 202 | f3md15m | 07100312_2258_F0001 | 07100313_0108_F0035 | 137 Files. DCRC: Preset 2; LED1: I=400uA, ON=100us, OFF=2500us; Data taking: 15min, Flash: 3min |
| 203 | f3md30m | 07100313_1619_F0001 | 07100313_1817_F0070 | 211 files. DCRC: Preset 2; LED1: I=400uA, ON=100us, OFF=2500us; Data taking: 30min, Flash: 3min |
| 204 | f3md60m | 07100314_1819_F0001 | 07100314_1819_F0142 | 142 Files. DCRC: Preset 2. LED1: I=400uA, ON=100us, OFF=2500us; Data taking: 60min, Flash: 3min |
| 205 | f12sd5m | 07100315_1212_F0001 | 07100315_1329_F0013 | 78 files. DCRC : Preset 2. LED1: I=1000uA, ON=100us, OFF=2500us; Data taking: 5min, Flash: 0.2min |
| 206 | f12md60m | 07100317_1014_F0001 | 07100317_1124_F0056 | 200 Files. DCRC: Preset 2. LED1: I=1000uA, ON=100us, OFF=2500us; Data taking: 60min, Flash: 0.2min |
Run 201 (f3md3m) Data Analysis
The last two digits '01' number the sub run in the main run. The name 'f3md3m' means flashing 3 min and them taking data for 3 min. Run201 analysis is in this page. Others can be found in the links below.
f3md60m
Used the same calibration parameters as f3md3m. f3md60m plots
Charge calibration
Energy calibration
We didn't have chance to take data with a radioactive source in this run. Instead the energy scale is determined by figuring out all the gains in the data acquisition and processing chain. The amplitude of the charge pulse Vf on the feedback capacitor Cf is Vf = Q/Cf. The gain from Vf to the digitizer is 10 (from Jeter). The resolution of the digitizer ADS1602 is VREF(3V)/(215-1) = 9.16×10-5V. The data recorded in the raw files is in units of digitizer bins. The corresponding RQs to the charge pulse amplitudes are QIOFvolts and QOOFvolts. The conversion factor from the pulse height to the RQ is calculated by a fit to the data in the following table.
| Table 1 QIOFvolts and the raw pulse heights for selected pulses | |||
|---|---|---|---|
| Series # | Event # | Pulse height | QIOFvolts |
| 071003121454 | 30003 | 2507 | 2.068234317712845e+03 |
| 071003121454 | 30005 | 1726.66 | 1.337342215483828e+03 |
| 071003121454 | 30008 | 1910.29 | 1.693592437027741e+03 |
| 071003121454 | 30009 | 5595.45 | 4.909761309521775e+03 |
| 071003121454 | 30011 | 11665.6 | 1.026814996876608e+04 |
| 071003121454 | 30013 | 1384.3 | 1.346075181854962e+03 |
| 071003121454 | 30015 | 8205.68 | 6.719766858661344e+03 |
| 071003121454 | 30018 | 2174.3 | 1.826612958114279e+03 |
| 071003121454 | 30020 | 2224.76 | 1.951766024551855e+03 |
| 071003121454 | 30025 | 10337.5 | 8.926799698376875e+03 |
| 071003121454 | 30026 | 1315.1 | 1.132223032238271e+03 |
| 071003121454 | 30028 | 1958.51 | 1.707466613788928e+03 |
| 071003121454 | 30030 | 4379.22 | 3.745395850027360e+03 |
| 071003121454 | 30031 | 1750.67 | 1.616621127118419e+03 |
| 071003121454 | 30032 | 6153.69 | 5.533173664480933e+03 |
So we have the relationships:
Vf = Q/CfVf*g = Ad*R; g=10, Ad is the pulse height, and R=9.16×10-5VAd*b = QIOFvoltsEr = N*ε = ε*Q/e; ε = 3eV and e=1.6×10-19C
From above we get
Er = Cf*R/(g*b*e)*ε*QIOFvolts = 198.56*QIOFvolts (eV) = qa**QIOFvolts (eV)
Crosstalk correction
Using a method similar to the standard one to to the crosstalk correction. The transformation is
- Qi = a*QIOFvolts + b*QOOFvolts
- Qo = c*QIOFvolts + d*QOOFvolts
By charge conservation we have
- Qi + Qo = a*QIOFvolts + b*QOOFvolts + c*QIOFvolts + d*QOOFvolts = QIOFvolts + QOOFvolts ⇒
- a + c = 1
- b + d = 1
To fully determine the transformation matrix, we need to know the slopes of the QIOFvolts and QOOFvolts bands in the (QIOFvolts, QIOFvolts) plot. The slopes can be obtained by fits to the bands directly. But the QOOFvolts band doesn't have so many events to define it clearly. We use a different technique to determine it. In the zeroth order, define
- Qpart = (QIOFvolts - QOOFvlots)/(QIOFvolts + QOOFvlots).
On the (QOOFvolts, QIOFvolts) plot, the QIOFvolts and QOOFvolts bands are represented by two straight lines QOOFvolts = ki * QIOFvolts and QOOFvolts = ko * QIOFvolts respectively. On the Qpart histogram, the two bands are two peaks at
- Qipart = (1 - ki)/(1 + ki)
- Qopart = (1 - ko)/(1 + ko)
The slope ko can be determined from the value of Qopart. ki, however, is very small, and Qipart is very close to 1. This method can not be used. We have to fit the QIOFvolts band to get ki.
Once we get ki and ko, redefine the charge partition as
- Qpart = (Qi - Qo)/(Qi + Qo)
Since Qi and Qo are orthogonal, we know Qpart = 1 for Qi band and Qpart = -1 for Qo band
- (a + b*ki - c - d*ki)/(a + b*ki + c + d*ki) = 1 ⇒
- c + d*ki = 0
- (a + b*ko - c - d*ko)/(a + b*ko + c + d*ko) = -1 ⇒
- a + b*ko = 0
Now the transformation can be completely determined.
Cuts and plots
After energy correction
- Qi = qa * QIOFvolts
- Qo = qa * QOOFvolts
The following cuts are used in the plots
| Plot | Cut |
|---|---|
| Qpart = (Qi - Qo)/(Qi + Qo) | No cut |
| (Qi, Qo) | c1 = Pt > 5 & inrange(Xdel, -50, 50) & inrange(Ydel, -50, 50) |
| Qi band fit | c2 = c1 & Qo > -10 & Qo < 10 |
| (Qsum, QSOFchisq) | 0.93 < Qpart & Qpart < 1.07 & 5 < Qsum & Qsum < 200 |
| Energy and crosstalk correction parameters | ||||
|---|---|---|---|---|
| qa | a | b | c | d |
| 198.56 eV | 0.999 | 0.072 | 0.001 | 0.928 |
Plots after energy and crosstalk corrections
Position dependence
There is no spectral line I can use to do the position correction. I just plot the events before the correction.
Phonon calibration
Relative amplitude correction
Make the following transformation:
- Pa = r(1) * PAOFev
- Pb = r(2) * PBOFev
- Pc = r(3) * PCOFev
- Pd = r(4) * PDOFev
Phonon energies are adjusted to let the four channels have the same distribution. The normalization is chosen to be
- r(1) + r(2) + r(2) + r(4) = 4
For Run 'f3md3m', the parameters are
- r = [1.2706 1.0191 1.0751 0.6353]
The phonon energy distribution after adjustment is
Overall calibration
The over all calibration minimizes the χ2 between the phonon energy Pt and the charge energy Qsum in the energy range 5 < Qsum < 100keV. This gives an overall factor of 2.6853e-08.
(Qsum, Pt), (Pt, Qsum/Pt)
Neutralization plot (Time, Qsum/Pt)
The yield Y = Qsum/Pt is plotted versus the time after flashing to see the neutralization change as time goes on.





















