Differences

This shows you the differences between two versions of the page.

--- fcs:eic:intro [2016/11/04 19:50]
asc [Reconstructed events]
+++ fcs:eic:intro [2016/11/06 19:18]
asc [Displaying reconstructed events]
@@ Line 111: / Line 111: @@
 Note that the plot is done for all electrons starting from 1 GeV (a region where the efficiency is quite small).
-====== Looking at other samples ======
+====== Looking at track resolution ======
 You can also look at single particles with flat pT, Eta and Phi. Here is an example this using a simple validation script {{:fcs:eic:mc_truth_flat.py|}}
-for the MCparticle container inside LCIO file:
+for the MCparticle container inside LCIO file. We will use the "flat pT" sample called [[http://atlaswww.hep.anl.gov/hepsim/info.php?item=242|pgun_elec_eta4_pt30gev]] reconstructed with the tag "rfull052" (SiEIC with 1T solenoid):
-<code>
+<code bash>
 hs-get pgun_elec_eta4_pt30gev%rfull052 data  2 100 # get 100 files in 2 threads in the directory "data"
 fpad mc_truth_flat.py
@@ Line 123: / Line 123: @@
 The resulting image is:
-{{:fcs:eic:mc_truth_pt_ele.png?400|}}
+{{:fcs:eic:mc_truth_pt_ele.png?500|}}
+Now, let's calculate track resolution. We will use this script {{:fcs:eic:track_resolution.py|track_resolution.py}} which
+makes fit of pT(reco)/pT(truth) and plots the resolution as a function of pT(truth).
+This command shown below makes resolution plots, fits with a Gaussian function, and plot resolution as a function of pT:
+<code bash>
+hs-get pgun_elec_eta4_pt30gev%rfull052 data  2 20 # get 20 files in 2 threads in the directory "data"
+fpad track_resolution.py
+</code>
+The final plot is shown below:
+{{:fcs:eic:pt_resolution.png?400|}}
 ====== How to continue ======
@@ Line 135: / Line 148: @@
   * fake rate (cannot be done with this sample!)
   * calculate Q2 and x_bj, check efficiency and resolutions for these variable
-===== Displaying reconstructed events =====
+====== Displaying reconstructed events ======
@@ Line 144: / Line 157: @@
 <code bash>
-wget http://mc.hep.anl.gov/asc/hepsim/events/misc/pgun/pgun_eta45_p30_elec//rfull051/pgun_elec30gev_001_pandora.slcio
+wget http://mc.hep.anl.gov/asc/hepsim/events/misc/pgun/pgun_eta45_p30_elec/rfull051/pgun_elec30gev_001_pandora.slcio
 </code>

ATLAS Support Center

User Tools

Site Tools

Differences

Page Tools