Data analysis

Basic steps how to evaluate data from test beam

1. Reading of raw data format from test beam (only for experts)
2. Basic preprocessing of data to prepare them for any physicist users and non-programmers
3. Generation of DST data files with including some special information inside
4. Prepare basic set of macros to help non experts to work with DST
5. To do some special analysis on telescopes to support simulation results of quality of test beam geometry
6. DUT analysis:
   1. Remove pedestal nad common noise of DUT
   2. Find hits/clusters on DUT
   3. Alignment of DUT and masking of bad channels
   4. Residual plots of DUT
   5. Efficiency of DUT, timing
   6. Noise occupancy of DUT
   7. Cluster size of DUT
   8. Bias of DUT
   9. Angle scans of DUT
   10. Interstrip analysis of all before of DUT
   11. Special analysis DAQ parameters of DUT
7. If it is possible to do analysis compare with simulation to find space resolution in DUT position
8. Laser tests

Raw data

Data preprocessing

• precise telescope intersection point using Eta correction alghoritm (R. Turchetta: Spatial resolution of silicon microstrip detectors, Nuclear Instruments and Methods in Physics Research A335 (1993) 44-58 North-Holland, S. Straulino at al.: Spatial resolution of double-sided silicon microstrip detectors for the PAMELA apparatus, Nuclear Instruments and Methods in Physics Research A 556 (2006) 100–114)
  • Plot of interstrip position of position of particle without Eta correction using method of centre of gravity (COG) from two highest or three highest strips 

    two highest strips give base on their signal S(R) and S(L) position:
    

    

    than Eta function for COG is:
    

    

  • Eta correction
  • Plot of interstrip position of position of particle after Eta correction

    New position with Eta correction will be:
    

    

  • Calculation of error of interstrip position after Eta correction - plot of average error value

    Error for Eta position is calculate by:
    

    

• telescope alignment was realized to find shift and rotation for each telescope plane, note that TMinuit class in ROOT was not working precisely if sharp cut for tracks was realized! Basic minimizing steps are:
• particle track find for each event where in each axis was exactly one hit in each telescope detected (independent for each axis)
• Chi2 for track was calculate
• intersection with detector under test (DUT) was calculate
• confidence and predicted region for estimation level 95% (2*sigma) was calculate in DUT position

DST production

Basic macros for evaluation

• UpgreatDST_depfet.cpp call .x UpgreatDST_depfet.cpp(RunNo) functions: macro show content of some event, all histograms, example how to write loop over some/all events, how to use tracking, how to find intersection point in telescope plane and make residuals include alignment information, how to open DST for append next histograms (danger for non-experts)
• LinearFit.cpp content some functions for find track, using of this function is also include: call .x LinearFit.cpp, or it is also in macro UpgreatDST_depfet.cpp

root [0] new TCanvas()      //to create panel for looking plots
(class TCanvas*)0x5a3a168
root [1] TFile *f = new TFile("./dst/DST2116.root")  //open DST file - in ROOT memory
root [2] DST->Draw("Channel2")          //to show histogram of all telescopes - mixture of all (wrong)
root [3] DST->Draw("Channel2[0][1][1]")  //to show histogram of first each hit in tel 1 axis 1
root [4] DST->Draw("Channel2[0][0][0]:Channel2[0][2][0]")  //scater plot - correlations btw tel0 and tel2
root [5] DST->Draw("Channel2[0][0][0]:Channel2[0][2][0]","flag==1")  //the same but only if exist track
root [6] DST->Draw("Channel2[0][0][0]:Channel2[0][0][1]","flag==1")  //correlation btw axis0 and axis1 (hitmap)
root [7] DST->Draw("Channel2[0][1][1]","Channel2[0][1][1]<400&Channel2[0][1][1]>200&Flag=1") // hitmap in smaller region
root [8] DST->Show()  //content of event 0
root [9] DST->Show(10)  //content of event 10

Space precision in DUT

1. telescope errors was set to value calculate base on theory of Eta correction
2. telescope errors was set to constant value 12 um

Results: Telescope errors was set to value calculate base on theory of Eta correction:

                                    Chi2Cut: 10 30 50 70 100 %
Tel 0 Resolution at Ifinity Energy, Axis 0 : 2.0 3.6 5.2 5.9 6.2
Tel 0 Resolution at Ifinity Energy, Axis 1 : 4.8 6.7 7.3 7.2 7.0

Tel 1 Resolution at Ifinity Energy, Axis 0 : 4.2 5.9 6.8 6.7 6.9
Tel 1 Resolution at Ifinity Energy, Axis 1 : 1.1 2.8 4.6 5.4 5.8

Tel 2 Resolution at Ifinity Energy, Axis 0 : 1.9 3.3 5.0 5.6 6.0
Tel 2 Resolution at Ifinity Energy, Axis 1 : 1.9 3.6 5.1 5.7 6.0

Tel 3 Resolution at Ifinity Energy, Axis 0 : 0.8 1.7 3.3 4.6 5.4
Tel 3 Resolution at Ifinity Energy, Axis 1 : 2.2 3.4 5.1 5.6 5.8

                                               Chi2Cut: 10 30 50 70 100 %
DUT: Confidence Region of Track Ifinity Energy, Tel 0 : 4.2 6.8 9.9 12.5 14.5
DUT: Confidence Region of Track Ifinity Energy, Tel 1 : 4.0 6.8 9.8 11.4 13.6
Chi2Cut: 10 30 50 70 100 %
DUT: Prediction of Intersection Ifinity Energy, Tel 0 : 9.3 15.1 22.2 28.3 32.1
DUT: Prediction of Intersection Ifinity Energy, Tel 1 : 8.9 15.2 22.2 27.1 30.6

• Summary_cut100_Tel.txt
• Summary_cut75_Tel.txt
• Summary_cut50_Tel.txt
• Summary_cut25_Tel.txt
• Summary_cut10_Tel.txt

• TelResiduals1_TelErr.gif - for each telescope each Chi2 cut SigmaResiduals vs Energy, similar for DUT mean of confidence and prediction region vs Energy
• TelResiduals2_TelErr.gif - for each telescope each Chi2 cut SigmaResiduals^2 vs 1/(Energy^2), similar for DUT mean of confidence and prediction region square vs 1/(Energy^2)
• R15_DUTResid2_Cut100_2153.gif - for DUT: hitmap, confidence and prediction region, chi2 for track, all for no Chi2 cut, energy 5GeV
• R15_AlignResid2_Cut100_2153.gif - for telescopes: residuals, all for no Chi2 cut, energy 5GeV
• R15_DUTResid2_Cut50_2153.gif - for DUT: hitmap, confidence and prediction region, chi2 for track, all for Chi2 cut 50% of tracks accept, energy 5GeV
• R15_AlignResid2_Cut50_2153.gif - for telescopes: residuals, all for Chi2 cut 50% of tracks accept, energy 5GeV

Telescope errors was set to constant value 12 um

                                    Chi2Cut: 10 30 50 70 100 %
Tel 0 Resolution at Ifinity Energy, Axis 0 : 1.4 1.8 3.3 4.3 5.9
Tel 0 Resolution at Ifinity Energy, Axis 1 : 1.2 2.0 2.8 4.2 5.5

Tel 1 Resolution at Ifinity Energy, Axis 0 : 1.6 2.6 3.9 5.8 7.1
Tel 1 Resolution at Ifinity Energy, Axis 1 : 1.6 2.2 3.8 5.1 6.7

Tel 2 Resolution at Ifinity Energy, Axis 0 : 1.5 2.3 3.6 5.5 6.7
Tel 2 Resolution at Ifinity Energy, Axis 1 : 1.2 2.4 3.9 5.3 6.4

Tel 3 Resolution at Ifinity Energy, Axis 0 : 1.4 1.8 3.2 4.1 5.6
Tel 3 Resolution at Ifinity Energy, Axis 1 : 1.1 1.9 3.0 4.2 5.4

                                               Chi2Cut: 10 30 50 70 100 %
DUT: Confidence Region of Track Ifinity Energy, Tel 0 : 3.0 5.0 7.1 9.4 11.4
DUT: Confidence Region of Track Ifinity Energy, Tel 1 : 2.6 4.4 7.4 9.1 10.2
Chi2Cut: 10 30 50 70 100 %
DUT: Prediction of Intersection Ifinity Energy, Tel 0 : 5.7 10.3 13.0 22.1 25.7
DUT: Prediction of Intersection Ifinity Energy, Tel 1 : 6.5 9.9 14.8 18.7 23.4

• Summary_cut100_12umTel.txt
• Summary_cut75_12umTel.txt
• Summary_cut50_12umTel.txt
• Summary_cut25_12umTel.txt
• Summary_cut10_12umTel.txt

• TelResiduals1_12umTelErr.gif - for each telescope each Chi2 cut SigmaResiduals vs Energy, similar for DUT mean of confidence and prediction region vs Energy
• TelResiduals2_12umTelErr.gif - for each telescope each Chi2 cut SigmaResiduals^2 vs 1/(Energy^2), similar for DUT mean of confidence and prediction region square vs 1/(Energy^2)
• R15_DUTResidConstErr2_Cut100_2153.gif - for DUT: hitmap, confidence and prediction region, chi2 for track, all for no Chi2 cut, energy 5GeV
• R15_AlignResidConstErr2_Cut100_2153.gif - for telescopes: residuals, all for no Chi2 cut, energy 5GeV
• R15_DUTResidConstErr2_Cut50_2153.gif - for DUT: hitmap, confidence and prediction region, chi2 for track, all for Chi2 cut 50% of tracks accept, energy 5GeV
• R15_AlignResidConstErr2_Cut50_2153.gif - for telescopes: residuals, all for Chi2 cut 50% of tracks accept, energy 5GeV

DUT analysis

// *********************************************************************************
// ***************** Here you add your code for DUT evaluation *********************
// *********************************************************************************

Remove pedestal nad common noise of DUT

Example how to use it is on: DUT1_DST.cpp

Example of real data acquised DUTMap.gif

Removing of pedestal DUTMap2.gif

Removing of pedestal - detail, common noise is visible DUTMap2a.gif

Removing of pedestal - detail, common noise is visible DUTMap2b.gif

Removing of pedestal and common noise - detail DUTMap3.gif

Removing of pedestal and common noise - detail DUTMap3a.gif

Histogram of noise distribution DUTMap3a.gif

Find hits/clusters on DUT

Example how to use it is on: DUT1_DST.cpp

Alignment of DUT and masking of bad channels

Residual plots of DUT

Efficiency of DUT and timing analysis

Noise occupancy of DUT

Cluster size of DUT

Bias of DUT

Angle scans of DUT

Interstrip analysis of all before of DUT

Special analysis

Analysis and simulations

Laser tests