The first analysis
The first data with one PTGC
- up side channel 18 missing, downside channel 45 47 missing
- remove signals whose peak id locate at channel 0, 1, 2, 15, 16, 17, 18, 19, 20, 21, 42, 43, 44, 45, 46, 47 for up and down side
- after this cut 9367 events
- both cluster size (up and down) >1 8956 events
- red: upside, black downside for the following plot
noise
- pedestal jumping simultaniously for all channel
- first calculate average value of all channel
- noise of each channel is to get the RMS after subtract the average
- pedestal of each channel is to using the mean plus the average.
- RMS of each channel (unit is ADC point, one point 0.1fC)
- signal (if signal > pedestal + 3* RMS)
- SM1(signal method 1): signal = data - pedestal
- SM2(signal method 2): signal = data - pedestal - 3*RMS
Data
- peak, second peak, third peak positon (strip number), down: remove bad channels :
- cluster center, cluster size, cluster number each event, down: after remove bad channels:
- peak, second peak, third peak value (x=1000 mean 100fC) , down: remove bad channels.
- Downside signal is smaller, wire to cathode is 1.5mm, upside 1.4mm
position resolution:
- left SM1: signal = data - pedestal
- right SM2: signal = data - pedestal - 3*RMS
- if ask peak value>105, resolution improve 5 micro.
- reason of the 2 pump round the main peak:
- left original singal
- middle SM1
- right SM2
- SM1
- left: ask cluster size of upside = downside:
- middle: cluster size of upside = downside
- right: peak value of events in left plot, black of events in middle plot
Use gaus fit to fit the cluster signals
- Other than charge gravity methods, this methods doesn't introduce much bias measurement
- Signal is obtained by minus only pedestal
- Fit the following hist with gaus from -0.1 to 0.1
sigal and noise re-analysis:
1,f(0-95), s(0-95) 共192道,
IsChNumberRepeated表示该事例存在信号道道号重复的次数,原始数据为所有记录的数据去掉道号重复次数多于1次的事例后。(S0和f0分布比其他的信号道宽)
pdf:
2,噪音分析,
- 没有信号的时候:
- 第j个事例第i个信号道的测量值
- 48个信号道的平均值 1/48 * Σ(m_i) = 1/48*ΣP_i + 1/48*ΣN_i + 1/48*ΣΔ_i
- = P + n_j + Δ_i n<N(if no big noisy channel)
- M_i = m_i - 1/48 * Σ(m_i) = P_i + N_i + Δ_i - ( P + n_j + Δ_i )
- = (P_i - P) + (N_i - n_j) if (N_i > 3 n_j) N_i - n_j ≈N_i
- = (P_i - P) + N_i
-
- distribuion of M_i centers at (P_i - P), σ = N_i
- If there is no signal: S_i = M_i - (P_i - P) = N_i
- 有信号的时候(evtnumber: j):
- 第j个事例第i个信号道的测量值 m_i = P_i (pedestal) + N_i (noise) + Δ_i (disturbence)
- 48个信号道的平均值 1/48 * Σ(m_i) = 1/48*ΣP_i + 1/48*ΣN_i + 1/48*ΣΔ_i + 1/48*S_j
- = P + n_j + Δ_i + 1/48*S_j n<N(if no big noisy channel)
- M_i = m_i - 1/48 * Σ(m_i) = P_i + N_i + Δ_i - ( P + n_j + Δ_i +1/48*S_j)
- = (P_i - P) + (N_i - n_j) - 1/48*S_j if (N_i > 3 n_j) N_i - n_j ≈N_i
- = (P_i - P) + N_i - 1/48*S_j
-
- distribuion of M_i centers at (P_i - P - Avr(1/48*S_j) ), σ = σ(N_i - 1/48*S_j)
- if there is signal, S_i = M_i - (P_i - P - Avr(1/48*S_j) ) = S + Gaus(N_i - 1/48*S_j)
- 去除信号方法:
- 每个事例:分为4组:f(0-47),f(48-95),s(0-47),s(48-95)分别对应一个前放电路板,每组48个数据求平均1/48 * Σ(m_i),从测量数据中扣除:m_i - 1/48 * Σ(m_i)。
- 对每道求RMS:σ(N_i - 1/48*S_j) 平均值Average:(P_i - P - Avr(1/48*S_j) ) 。
- 判断信号:如果测量值m_i > ( (P_i - P - Avr(1/48*S_j) ) + σ(N_i - 1/48*S_j) ) 为信号,否则为噪音
- 将判定为信号的道设为0,然后剔除
- 剔除信号后,噪音按照上述步骤重新计算得到没有信号影响的P_i - P 和每个信号道的噪音 σ(N_i )
- 共192道:n1u(0-47), n1d(0-47), n2(0-47)n3(0-47) 其中n1d_0分布较宽,n2_0分布较宽,且有等间距的几个峰结构。 n1d_0和n2_0是否对应上面的f0,s0?
- 去掉信号之后的噪音二维图(噪音强度m_i v.s eventnumber), project to eventnumber给出基值P_i值随着时间的变化趋势
- 去掉扰动以后的噪音二维散点图(每一道实际噪音分布)
3 ,pdf, 图一: 噪音signa vs channelnumber, 图二: 噪音平均值,vs channnel number
4, pdf:X2measure - X2 calculate ( 每个信号道的测量值=信号-pedestal if 信号-pedestal > 3 signa )
Latex rendering error!! dvi file was not created.