Differences

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

--- fcs:eic:intro [2016/10/04 17:30] – [Reconstructed events] asc
+++ fcs:eic:intro [2016/12/09 16:34] (current) – [Hit timing information] asc
@@ Line 110: / Line 110: @@
 Note that the plot is done for all electrons starting from 1 GeV (a region where the efficiency is quite small).
+====== Hit timing information ======
+You can also look at timing of silicon tracker using this script: {{:fcs:eic:mc_hits_time.py|eic:mc_hits_time.py}}
+<code bash>
+hs-get pgun_elec_eta4_pt30gev%rfull052 data # 2 100  get 100 files in 2 threads in the directory "data"
+fpad eic:mc_hits_time.py
+</code>
+The plot is shown here:
+{{:fcs:eic:electrons_hit_timing.png?400|}}
+====== 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 162: @@
   * 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 128: / Line 171: @@
 <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>