Differences

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

--- fcs:eic:intro [2016/10/04 17:29]
asc
+++ fcs:eic:intro [2016/11/06 16:31]
asc [Looking at track resolution]
@@ Line 109: / Line 109: @@
 {{:fcs:eic:mc_electrons.png?400|}}
-Note that the plot is done for all electrons starting from 1 GeV (where efficiency is low).
+Note that the plot is done for all electrons starting from 1 GeV (a region where the efficiency is quite small).
+====== 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. 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 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
+</code>
+The resulting image is:
+{{: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 ======
 Now you can use this template to study:
@@ Line 119: / 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 ======

ATLAS Support Center

User Tools

Site Tools

Differences

Page Tools