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)
      noise.jpg

  • 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 :
    peakpos.jpg

  • cluster center, cluster size, cluster number each event, down: after remove bad channels:
    cluster.jpg

  • 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
      peakval.jpg

position resolution:

  • left SM1: signal = data - pedestal
  • right SM2: signal = data - pedestal - 3*RMS
  • if ask peak value>105, resolution improve 5 micro. reso_SM1.jpg

  • reason of the 2 pump round the main peak:
    • left original singal
    • middle SM1
    • right SM2
      reso3.jpg

  • 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
      reso2.jpg

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
    a.gif

sigal and noise re-analysis:

1,f(0-95), s(0-95) 共192道, IsChNumberRepeated表示该事例存在信号道道号重复的次数,原始数据为所有记录的数据去掉道号重复次数多于1次的事例后。(S0和f0分布比其他的信号道宽)

pdf:

2,噪音分析,

  • 没有信号的时候:
    • 第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
    • = 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.

Topic attachments
I Attachment History Action Size Date Who Comment
PowerPointppt PTGCana.ppt r1 manage 371.0 K 2010-02-26 - 11:06 ChengguangZhu  
GIFgif a.gif r1 manage 25.7 K 2010-03-05 - 09:20 ChengguangZhu  
JPEGjpg cluster.jpg r1 manage 117.8 K 2010-02-26 - 10:30 ChengguangZhu  
JPEGjpg noise.jpg r1 manage 90.7 K 2010-02-26 - 10:36 ChengguangZhu  
JPEGjpg peakpos.jpg r1 manage 134.0 K 2010-02-26 - 10:29 ChengguangZhu  
JPEGjpg peakval.jpg r1 manage 137.1 K 2010-02-26 - 10:35 ChengguangZhu  
JPEGjpg reso2.jpg r1 manage 77.0 K 2010-02-26 - 10:52 ChengguangZhu  
JPEGjpg reso3.jpg r1 manage 27.7 K 2010-02-26 - 10:52 ChengguangZhu  
JPEGjpg reso_SM1.jpg r1 manage 97.1 K 2010-02-26 - 10:51 ChengguangZhu  
Edit | Attach | Watch | Print version | History: r4 < r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r4 - 2011-01-26 - ChengguangZhu
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2024 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback