HFI detection chain

Detection chain

(Francesco Piacentini)

This entire section has been moved to this page.

Bolometers

The heart of the HFI - the detectors - are bolometers, solid-state devices in which the incoming radiation dissipates its energy as heat that increases the temperature of a thermometer. The instrument Flight Model total number of bolometers is 52, split into 6 channels at central frequencies of 100, 143, 217, 353, 545, and 857GHz. Two extra bolometers not optically coupled to the telescope are added to the focal plane to monitor thermal noise (Dark Bolometers).

Thirty two of these bolometers are polarization sensitive allowing a map of the CMB polarisation to be built.

Bolometers consist of (i) an absorber that transforms the in-coming radiation into heat; (ii) a thermometer that is thermally linked to the absorber and measures the temperature changes; and (iii) a weak thermal link to a thermal sink, to which the bolometer is attached.

The detectors are semi-conducting NTD thermistor bolometers. There are two kinds of detector modules:

• SWB, spider web bolometers, are mounted on an absorbing spider-web of metallised silicon nitride
• PSB, polarization sensitive bolometers, are mounted or on a parallel absorbing grid of metallised silicon nitride.

In the spider-web bolometers, or SWBs (Bock et al. 1995; Mauskopf et al. 1997), the absorbers consist of metallic grids deposited on a Si3 N4 substrate in the shape of a spider web. The mesh design and the impedance of the metallic layer are adjusted to match vacuum impedance and maximize the absorption of millimeter waves, while minimizing the cross section to particles. The absorber is designed to offer equal impedance to any linearly polarised radiation. Nevertheless, the thermometer and its electrical leads define a privileged orientation (Fig. SWB) that makes the SWBs slightly sensitive to polarisation, as detailed in (Rosset et al. 2010). The thermometers are made of neutron transmuted doped (NTD) germanium (Haller et al. 1996), chosen to have an impedance of about 10 MΩ at their operating temperature.

Picture of a spider web bolometer. This a 143 GHz module, The temperature sensor is at the center of the absorbing grid

The absorber of PSBs is a rectangular grid (Fig. PSB) with metallization in one direction (Jones et al. 2003). Electrical fields parallel to this direction develop currents and then deposit some power in the grid, while perpendicular electrical fields propagate through the grid without significant interaction. A second PSB perpendicular to the first one absorbs the other polarisation. Such a PSB pair measures two polarisations of radiation collected by the same horns and filtered by the same devices, which minimises the systematic effects: differences between polarised beams collected by a given horn are typically less than −30 dB of the peak. The differences in the spectral responses of a PSB pair also proved to be a few percent in the worst case. Each pair of PSBs sharing the same horn is able to measure the intensity Stokes parameter and the Q parameter associated with its local frame. An associated pair of PSBs rotated by 45◦ scans exactly the same line (if the geometrical alignment is perfect), providing the U Stokes parameter.

Picture of a polarization sensitive bolometer. This a 217 GHz module, The temperature sensor is at the an edge of the absorbing grid

The detectors operate at a temperature close to 100 mK, while the filters are distributed on the 100mK, 1.6 K, and 4 K stages in such a way that the heat load on the coldest stages is minimized to limit the heat load on the detectors and to decrease the heat lift requirement and thus enhance the mission lifetime. The self-emission of the 4K stage is minimised to limit the photon noise contribution on the detectors from the instrument. The HFI Focal Plane Unit accommodates sub-millimiter absorbing material in order to decrease the scattering inside it.

to be change into a real table

Focal plane layout

The layout of the detectors in the focal plane is defined to cope with the scanning strategy. The HFI horns are positioned at the centre of the focal plane, where the optical quality is good enough for the high frequencies. The curvature of rows compensates for the distortion of images by the telescope. A pair of identical SWB will scan the same circle on the sky to provide additional redundancy. Similar horns feeding PSBs are also aligned so that two pairs of PSBs rotated by 45◦ with respect to each other scan the same line. This will provide the Q and U Stokes parameters with minimal correction for the pointing (Rosset et al. 2010). Residual systematics will come from the differences between the beam shapes of the two horns. In all cases except for the 100 GHz horns, a measurement is also done by a pair of similar channels shifted by 1.25 arcminutes in the cross-scan direction, to ensure adequate sampling. In the figures focal plane layout is reported.

Focal plane layout as seen from outside the celestial sphere. Each spot represents an horn in the focal plane. Coordinates are in pointing direction with respect to the telescope boresight. PSBs detectors are indicated with a cross oriented as the two polarimeters axes. Scan direction is from left to right.

Focal plane layout as seen from outside the celestial sphere. In this case PSBs orientations are reported. Red rod is for a elements of PSBs, blue rod is for b elements.

Readout

The AC bias readout electronics of the HFI instrument (Gaertner et al. 1997) includes a number of original features proposed by several laboratories (CRTBT in Grenoble, CESR in Toulouse and IAS in Orsay), which were validated on the Diabolo experiment and on the balloon-borne Archeops experiment. It was developed for space by the CESR in Toulouse.

The particular features of the HFI AC bias readout are mainly

• i)that the cold load resistors were replaced by capacitors because they have no Johnson noise;
• ii)that the detectors are biased by applying a triangular voltage to the load capacitors

which produces a square current [Ibias ] in the capacitors, and a square voltage [T bias ] that compensates for the stray capacitance of the wiring (producing a nearly square bias current into the bolometer;

• iii)that a square offset com-pensation signal is subtracted to the bolometric signal to minimise the amplitude of the signal that has to be amplified and digitized;
• iv)that the electronic scheme is symmetrical and uses a differential amplification scheme to optimize the immunity to electromagnetic interferences;
• v)and finally that every parameter of the REU can be set by commands, which

is made possible by using digital-to-analog and analog-to-digital converters extensively.

The readout electronics consist of 72 channels designed to perform low noise measurements of the impedance of 52 bolometers, two blind bolometers, 16 accurate low temperature thermometers, all in the 10 MΩ range, one resistor of 10 MΩ and one capacitor of 100 pF. The semiconductor bolometers and thermometers of Planck-HFI operate at cryogenic temperature around 100 mK on the focal plane, with impedance of about 10 MΩ when biased at the optimal current. The readout electronics on the contrary have to operate at “room” temperature (300 K). The distance between the two extremities of the readout chain is about 10 m and could represent a point of extreme susceptibility to electromagnetic interference. The readout electronic chain was split into three boxes. These are the JFET box, located on the 50 K stage of the satellite at 2.2 m from the focal plane, the pre-amplifier unit (PAU), located 1.8m further at 300 K, and the REU, located on the opposite side of the satellite, 5 m away. Each of the three boxes (JFET, PAU and REU) consists of 12 belts of six channels. The first nine belts are dedicated to the bolometers, and the three last ones to the accurate thermometers, the resistor and the capacitor (see figure Organization of the HFI readout).

Principles of readout electronics. The three modules of the chain are shown: JFet Box, Pre-Amplifier Unit (PAU) and Readout Electronic Unit (REU)

Organization of the HFI readout. Each row represents a belt. Each belt has 6 channels. (to be change into a real table ?)

Principles of the readout electronics.

See figure [Principles of readout electronics]. The bolometer is biased by a square wave AC current obtained by the differentiation of a triangular voltage through a load capacitance, in a completely differential architecture. The presence of the stray capacitance due to losses of charge in the wiring requires a correction of the shape of the square bias current by a transient voltage. Thus the bias voltage generation is controlled by the two parameters I-bias and T-bias that express the amplitude of the triangular and transient voltage. The compensation voltage added to the bolometric signal to optimize the dynamic of the chain is controlled by the V-bal parameter.

Parameters of the Readout Unit can be set to optimize the detectors performance.

The modulation frequency of the AC bias system, fmod of the square bias current can be tuned from 70 Hz to 112 Hz by the telecommand parameters: Nsample, which defines the number of samples per half period of modulated signal, fdiv which determines the sampling frequency of the ADC.

The optimal frequencies are around 90 Hz.

Each channel has its own settings for the following parameters: I-bias, amplitude of the triangular bias voltage;

T-bias, amplitude of the transient bias voltage;

V-bal, amplitude of the square compensation voltage;

G-amp, value of the programmable gain of the REU [1/3, 1, 3, 7.6];

N-blank, number of blanked samples at the beginning of halfperiod not taken into account during integration of the signal;

S-phase, phase shift when computing the integrated signal.

All these parameters influence the effective response of the detection chains, and were optimized during the calibration campaigns and confirmed during the calibration and performance verification (CPV) phase following the launch of Planck. The scientific signal is provided by the integral of the signal on each half-period, between limits fixed by the S-phase and N-blank parameters.

The interaction of modulated readout electronics with semiconductor bolometers is rather different from that of a classical DC bias readout (Jones 1953). The differences were seen during the calibration of the HFI, although the readout electronics was designed to mimic the operation of a DC biased bolometric system as far as possible. With the AC readout the maximum of responsivity is lower and is obtained for higher bias current in the bolometer with respect to the DC model. This is caused by the stray capacitance in the wiring which has negligible effects for a DC bias and a major effect for an AC bias. In our case, a stray capacitance of 150 pF induces increases of NEP ranging from 4% (857 GHz bolometers) to 10% (100 GHz bolometers) and also affects the HFI time response. Details of the effect of the HFI AC bias system into the time response of the detectors are discussed in the Time Response Section.

JFETs

Given the high impedance of the bolometers and the length of the connecting cables, it it is essential that the impedance of the signal is lowered as close as possible to the detectors. In our system this is accomplished by means of JFET source followers, located in boxes connected to the 50 K stage. The JFET box has been designed, developed and tested in the Observational Cosmology group in the Physics Department of the University of Rome "La Sapienza" (Brienza D. et al 2006). There are two JFETs per channel, since the readout is fully differential, and they provide a current amplification of the signal while keeping the voltage amplification very close to unity.

Inside the box, the JFETs are mounted on a thermally insulated plate with an active temperature control to keep them at the optimal temperature of 110 K. With a dissipated power lower than 240 mW, mainly produced by the JFETs and the source resistors, we obtained a noise power spectral density of less than 3 nV Hz1/2 for the frequency range of interest. This increases the total noise ofall bolometer channels by less than 5%.

Data compression

The output of the readout electronics unit (REU) consists of one value for each of the 72 science channels for each half-period of modulation. This number, [math]S_{REU}[/math] , is the sum of the 40 16-bit ADC signal values obtained within the given half-period. The data processor unit (DPU) performs a lossy quantization of [math]S_{REU}[/math].

We define a compression slice of 254 [math]S_{REU}[/math] values, corresponding to about 1.4 s of observation for each detector and to a strip of sky about 8 degrees long. The mean [math]\langle S_{REU} \rangle[/math] of the data within each compression slice is computed, and data are demodulated using this mean:

[math]S_{demod,i} = (S_{REU,i} − \langle S_{REU} \rangle ) ∗ (−1)^{i}[/math]

where [math]1 \lt i \lt 254[/math] is the running index within the compression slice.

The mean [math]\langle S_{demod} \rangle[/math] of the demodulated data [math]S_{demod,i}[/math] is computed and subtracted, and the resulting data slice is quantized according to a step size Q that is fixed per detector:

[math] S_{DPU,i} = \mbox{round} \left[( S_{demod,i} − \langle S_{demod} \rangle) /Q \right ] [/math]

This is the lossy part of the algorithm: the required compression factor, obtained through the tuning of the quantization step Q, adds a measure of noise to the data. Assuming Gaussian white noise with standard deviation [math]\sigma[/math], a quantization setting of [math]\sigma[/math]/Q = 2 adds 1% to the noise (Pajot et al. 2010; Pratt 1978). The value of [math]\sigma[/math] was determined at the end of the CPV phase after subtraction of the signal from the timeline.

The two means [math]\langle S_{REU} \rangle[/math] and [math]\langle S_{demod} \rangle[/math] are computed as 32-bit words and sent through the telemetry, together with the [math]S_{DPU,i}[/math] values. Variable-length encoding of the [math]S_{DPU,i}[/math] values is performed on board, and the inverse decoding is applied on ground.

Optimal use of the bandpass available for the downlink (75 kb/s average for HFI science) was obtained initially by using a value of Q = [math]\sigma[/math]/2.5 for all bolometer signals. After 12 December 2009, only for the 857 GHz detectors, the value was reset to Q = [math]\sigma[/math]/2.0 to avoid data loss due to exceeding the limit of downlink rate.

Time_response

The HFI bolometers and readout electronics have a finite response time to changes in incident optical power. The bolometers are thermal detectors of radiation whose response time is determined by the thermal circuit defined by the heat capacity of the detector and thermal conductivity.

Due to Planck's nearly constant scan rate, the time response is degenerate with the optical beam. However, because of the long time scale effects present in the time response, the time response is deconvolved from the data in the processing of the HFI data (see TOI processing).

The time response of the HFI bolometers and readout electronics is modeled as a Fourier domain transfer function (called the LFER4 model) consisting of the product of an bolometer thermal response [math]F(\omega)[/math] and an electronics response [math]H'(\omega)[/math].

[math]\label{LFER4def}TF^{LFER4}(\omega) = F(\omega) H'(\omega)[/math]

LFER4 model

If we write the input signal (power) on a bolometer as [math]\label{bol_in} s_0(t)=e^{i\omega t} [/math] the bolometer physical impedance can be written as: [math]\label{bol_out} s(t)=e^{i\omega t}F(\omega) [/math] where [math]\omega[/math] is the angular frequency of the signal and [math]F(\omega)[/math] is the complex intrinsic bolometer transfer function. For HFI the bolometer transfer function is modelled as the sum of 4 single pole low pass filters: [math]\label{bol_tf} F(\omega) = \sum_{i=0,4} \frac{a_i}{1 + i\omega\tau_i} [/math] The modulation of the signal is done with a square wave, written here as a composition of sine waves of decreasing amplitude: [math]\label{sigmod} s'(t)=e^{i\omega t}F(\omega)\sum_{k=0}^{\infty} \frac{e^{i\omega_r(2k+1)t}-e^{-i\omega_r(2k+1)t}}{2i(2k+1)} [/math] where we have used the Euler relation [math]\sin x=(e^{ix}-e^{-ix})/2i[/math] and [math]\omega_r[/math] is the angular frequency of the square wave. The modulation frequency is [math]f_{mod} = \omega_r/2\pi[/math] and was set to [math]f_{mod} = 90.18759 [/math]Hz in flight. This signal is then filtered by the complex electronic transfer function [math]H(\omega)[/math]. Setting: [math]\omega_k^+=\omega+(2k+1)\omega_r[/math] [math]\omega_k^-=\omega-(2k+1)\omega_r[/math] we have: [math]\label{sigele} \Sigma(t)=\sum_{k=0}^\infty\frac{F(\omega)}{2i(2k+1)}\left[H(\omega_k^+)e^{i\omega_k^+t}-H(\omega_k^-)e^{i\omega_k^-t}\right] [/math] This signal is then sampled at high frequency ([math]2 f_{mod} NS[/math]). [math]NS[/math] is one of the parameters of the HFI electronics and corresponds to the number of high frequency samples in each modulation semi-period. In order to obtain an output signal sampled every [math]\pi/\omega_r[/math] seconds, we must integrate on a semiperiod, as done in the HFI readout. To also include a time shift [math]\Delta t[/math], the integral is calculated between [math]n\pi/\omega_r+\Delta t[/math] and [math](n+1)\pi/\omega_r+\Delta t[/math] (with [math]T=2 \pi/\omega_r[/math] period of the modulation). The time shift [math]\Delta t[/math] is encoded in the HFI electronics by the parameter [math]S_{phase}[/math], with the relation [math]\Delta t = S_{phase}/NS/f_{mod} [/math].

After integration, the [math]n[/math]-sample of a bolometer can be written as [math]\label{eqn:output} Y(t_n) = (-1)^n F(\omega) H'(\omega) e^{i t_n \omega} [/math] where [math]\label{tfele} H'(\omega) = \frac 12 \sum_{k=0}^\infty e^{-i(\frac{\pi\omega}{2\omega_r}+\omega\Delta t)} \Bigg[ \frac{H(\omega_k^+)e^{i\omega_k^+ \Delta t}}{(2k+1)\omega_k^+} \left(1-e^{\frac{i\omega_k^+\pi}{\omega_r}}\right) \\ - \frac{H(\omega_k^-)e^{i\omega_k^- \Delta t}}{(2k+1)\omega_k^-} \left(1-e^{\frac{i\omega_k^-\pi}{\omega_r}}\right) \Bigg] [/math]

The output signal in equation eqn:output can be demodulated (thus removing the [math](-1)^n[/math]) and compared to the input signal in equation bol_in. The overall transfer function is composed of the bolometer transfer function and the effective electronics transfer function, [math]H'(\omega)[/math]: [math]TF(\omega) = F(\omega) H'(\omega) [/math]

The shape of [math]H(\omega)[/math] is obtained combining low and high-pass filters with Sallen Key topologies (with their respective time constants) and accounting also for the stray capacitance low pass filter given by the bolometer impedance combined with the stray capacitance of the cables. The sequence of filters that define the electronic band-pass function [math]H(\omega) = h_0*h_1*h_2*h_3*h_4*h_{5}[/math] are listed in table table:readout_electronics_filters.

Parameters of LFER4 model

The LFER4 model has are a total of 10 parameters([math]A_1[/math],[math]A_2[/math],[math]A_3[/math],[math]A_4[/math],[math]\tau_1[/math],[math]\tau_2[/math],[math]\tau_3[/math],[math]\tau_4[/math],[math]S_{phase}[/math],[math]\tau_{stray}[/math]) 9 of which are independent, for each bolometer. The free parameters of the LFER4 model are determined using in-flight data in the following ways:

• [math]S_{phase}[/math] is fixed at the value of the REU setting.
• [math]\tau_{stray}[/math] is measured during the QEC test during CPV.
• [math]A_1[/math], [math]\tau_1[/math], [math]A_2[/math], [math]\tau_2[/math] are fit forcing the compactness of the scanning beam.
• [math]A_3[/math], [math]\tau_3[/math], [math]A_{4}[/math] [math]\tau_4[/math] are fit by forcing agreement of survey 2 and survey 1 maps.
• The overall normalization of the LFER4 model is forced to be 1.0 at the signal frequency of the dipole.

The details of determining the model parameters are given in (reference P03c paper) and the best-fit parameters listed here in table table:LFER4pars.

HFI electronics filter sequence

HFI electronics filter sequence. We define $s = i \omega$

Filter	Parameters	Function
0. Stray capacitance low pass filter	[math]\tau_{stray}= R_{bolo} C_{stray}[/math]	[math]h_0 = \frac{1}{1.0+\tau_{stray}*s}[/math]
1. Low pass filter	[math]R_1=1[/math]k[math]\Omega[/math] [math]C_1=100[/math]nF	[math]h_1 = \frac{2.0+R_1*C_1*s}{2.0*(1.0+R_1*C_1*s)}[/math]
2. Sallen Key high pass filter	[math]R_2=51[/math]k[math]\Omega[/math] [math]C_2=1\mu[/math]	[math]h_2= \frac{(R_2*C_2*s)^2}{(1.0+R_2*C_2*s)^2}[/math]3
3. Sign reverse with gain		[math]h_3=-5.1[/math]
4. Single pole low pass filter with gain	[math]R_4=10[/math]k[math]\Omega[/math] [math]C_4=10[/math]nF	[math]h_4= \frac{1.5}{1.0+R_4*C_4*s}[/math]
5. Single pole high pass filter coupled to a Sallen Key low pass filter	[math]R_9=18.7[/math]k[math]\Omega[/math] [math]R_{12}=37.4[/math]k[math]\Omega[/math] [math]C=10.0[/math]nF [math]R_{78}=510[/math]k[math]\Omega[/math] [math]C_{18}=1.0\mu[/math]F [math]K_3 = R_9^2*R_{78}*R_{12}^2*C^2*C_{18}[/math] [math]K_2 = R_9*R_{12}^2*R_{78}*C^2+R_{9}^2*R_{12}^2*C^2+R_9*R_{12}^2*R_{78}*C_{18}*C[/math] [math]K_1 =R_9*R_{12}^2*C+R_{12}*R_{78}*R_9*C_{18}[/math]	[math]h_{5} = \frac{2.0*R_{12}*R_9*R_{78}*C_{18}*s}{s^3*K_3 + s^2*K_2+ s*K_1 + R_{12}*R_9 } [/math]

Parameters for LFER4 model.

Bolometer	[math]A_1[/math]	[math]\tau_1[/math] (s)	[math]A_2[/math]	[math]\tau_2[/math] (s)	[math]A_3[/math]	[math]\tau_3[/math] (s)	[math]A_4[/math]	[math]\tau_4[/math] (s)	[math]\tau_{stray}[/math] (s)	[math]S_{phase}[/math] (s)
100-1a	0.392	0.01	0.534	0.0209	0.0656	0.0513	0.00833	0.572	0.00159	0.00139
100-1b	0.484	0.0103	0.463	0.0192	0.0451	0.0714	0.00808	0.594	0.00149	0.00139
100-2a	0.474	0.00684	0.421	0.0136	0.0942	0.0376	0.0106	0.346	0.00132	0.00125
100-2b	0.126	0.00584	0.717	0.0151	0.142	0.0351	0.0145	0.293	0.00138	0.00125
100-3a	0.744	0.00539	0.223	0.0147	0.0262	0.0586	0.00636	0.907	0.00142	0.00125
100-3b	0.608	0.00548	0.352	0.0155	0.0321	0.0636	0.00821	0.504	0.00166	0.00125
100-4a	0.411	0.0082	0.514	0.0178	0.0581	0.0579	0.0168	0.37	0.00125	0.00125
100-4b	0.687	0.0113	0.282	0.0243	0.0218	0.062	0.00875	0.431	0.00138	0.00139
143-1a	0.817	0.00447	0.144	0.0121	0.0293	0.0387	0.0101	0.472	0.00142	0.00125
143-1b	0.49	0.00472	0.333	0.0156	0.134	0.0481	0.0435	0.27	0.00149	0.00125
143-2a	0.909	0.0047	0.0763	0.017	0.00634	0.1	0.00871	0.363	0.00148	0.00125
143-2b	0.912	0.00524	0.0509	0.0167	0.0244	0.0265	0.0123	0.295	0.00146	0.00125
143-3a	0.681	0.00419	0.273	0.00956	0.0345	0.0348	0.0115	0.317	0.00145	0.00125
143-3b	0.82	0.00448	0.131	0.0132	0.0354	0.0351	0.0133	0.283	0.00161	0.000832
143-4a	0.914	0.00569	0.072	0.0189	0.00602	0.0482	0.00756	0.225	0.00159	0.00125
143-4b	0.428	0.00606	0.508	0.00606	0.0554	0.0227	0.00882	0.084	0.00182	0.00125
143-5	0.491	0.00664	0.397	0.00664	0.0962	0.0264	0.0156	0.336	0.00202	0.00139
143-6	0.518	0.00551	0.409	0.00551	0.0614	0.0266	0.0116	0.314	0.00153	0.00111
143-7	0.414	0.00543	0.562	0.00543	0.0185	0.0449	0.00545	0.314	0.00186	0.00139
217-5a	0.905	0.00669	0.0797	0.0216	0.00585	0.0658	0.00986	0.342	0.00157	0.00111
217-5b	0.925	0.00576	0.061	0.018	0.00513	0.0656	0.0094	0.287	0.00187	0.00125
217-6a	0.844	0.00645	0.0675	0.0197	0.0737	0.0316	0.0147	0.297	0.00154	0.00125
217-6b	0.284	0.00623	0.666	0.00623	0.0384	0.024	0.0117	0.15	0.00146	0.00111
217-7a	0.343	0.00548	0.574	0.00548	0.0717	0.023	0.0107	0.32	0.00152	0.00139
217-7b	0.846	0.00507	0.127	0.0144	0.0131	0.0479	0.0133	0.311	0.00151	0.00139
217-8a	0.496	0.00722	0.439	0.00722	0.0521	0.0325	0.0128	0.382	0.00179	0.00111
217-8b	0.512	0.00703	0.41	0.00703	0.0639	0.0272	0.0139	0.232	0.00173	0.00125
217-1	0.0136	0.00346	0.956	0.00346	0.0271	0.0233	0.00359	1.98	0.00159	0.00111
217-2	0.978	0.00352	0.014	0.0261	0.00614	0.042	0.00194	0.686	0.0016	0.00125
217-3	0.932	0.00355	0.0336	0.00355	0.0292	0.0324	0.00491	0.279	0.00174	0.00125
217-4	0.658	0.00135	0.32	0.00555	0.0174	0.0268	0.00424	0.473	0.00171	0.00111
353-3a	0.554	0.00704	0.36	0.00704	0.0699	0.0305	0.0163	0.344	0.0017	0.00125
353-3b	0.219	0.00268	0.671	0.00695	0.0977	0.0238	0.0119	0.289	0.00157	0.00111
353-4a	0.768	0.00473	0.198	0.00993	0.0283	0.0505	0.00628	0.536	0.00181	0.00125
353-4b	0.684	0.00454	0.224	0.0108	0.0774	0.08	0.0149	0.267	0.00166	0.00111
353-5a	0.767	0.00596	0.159	0.0124	0.0628	0.0303	0.0109	0.357	0.00156	0.00111
353-5b	0.832	0.00619	0.126	0.0111	0.0324	0.035	0.0096	0.397	0.00166	0.00111
353-6a	0.0487	0.00176	0.855	0.006	0.0856	0.0216	0.0105	0.222	0.00199	0.00125
353-6b	0.829	0.00561	0.127	0.00561	0.0373	0.0252	0.00696	0.36	0.00228	0.00111
353-1	0.41	0.000743	0.502	0.00422	0.0811	0.0177	0.0063	0.329	0.00132	0.00097
353-2	0.747	0.00309	0.225	0.00726	0.0252	0.0447	0.00267	0.513	0.00154	0.00097
353-7	0.448	0.0009	0.537	0.0041	0.0122	0.0273	0.00346	0.433	0.00178	0.00125
353-8	0.718	0.00223	0.261	0.00608	0.0165	0.038	0.00408	0.268	0.00177	0.00111
545-1	0.991	0.00293	0.00743	0.026	0.00139	2.6	0	0	0.00216	0.00111
545-2	0.985	0.00277	0.0128	0.024	0.00246	2.8	0	0	0.00187	0.00097
545-4	0.972	0.003	0.0277	0.025	0.000777	2.5	0	0	0.00222	0.00111
857-1	0.974	0.00338	0.0229	0.025	0.00349	2.2	0	0	0.00176	0.00111
857-2	0.84	0.00148	0.158	0.00656	0.00249	3.2	0	0	0.0022	0.00125
857-3	0.36	4.22e-05	0.627	0.0024	0.0111	0.017	0.002	1.9	0.00152	0.00126
857-4	0.278	0.0004	0.719	0.00392	0.00162	0.09	0.00152	0.8	0.00149	0.000558