Differences

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

--- fcs:eic:intro [2016/09/19 22:41] – [Reconstructed events] asc
+++ fcs:eic:intro [2016/12/09 16:34] (current) – [Hit timing information] asc
@@ Line 77: / Line 77: @@
 ====== Reconstructed events ======
 Each detector in HepSim has "reconstruction" tag, i.e. a data folder where files after simulation and reconstruction are located. How to find the name of such folder?
-Click this link: [[http://atlaswww.hep.anl.gov/hepsim/lookup.php?info=sieic1|check sieic1 tag]]. It shows that this detector has the tag **rfull051**.
+Click this link: [[http://atlaswww.hep.anl.gov/hepsim/lookup.php?info=sieic1|check sieic1 tag]]. It shows that this detector has the [[http://atlaswww.hep.anl.gov/hepsim/list.php?find=rfull051|rfull051 tag]].
-The reconstruction tags are shown in the info page of this sample [[http://atlaswww.hep.anl.gov/hepsim/info.php?item=214|pgun_elec_eta45_30gev]]. List available files after the full simulation and reconstruction using **rfull051** tag:
+The reconstruction tags are shown in the info page of this sample [[http://atlaswww.hep.anl.gov/hepsim/info.php?item=214|pgun_elec_eta45_30gev]]. List available files after the full simulation and reconstruction using [[http://atlaswww.hep.anl.gov/hepsim/list.php?find=rfull051|rfull051 tag]]:
 <code bash>
@@ Line 93: / Line 93: @@
 You can run such scripts using "fpad" (Linux/Mac) or jaspp GUI (Windows).
-Look at the tutorial [[community:hepsim:usage_full|here]]. It also has a C++ package to read SLCIO files.
+<note tip>Look at the tutorial [[community:hepsim:usage_full|here]] which shows how  to use C++/ROOT  for  reading  SLCIO files. We show this Python example since it does not require compilation, and you can run it on Windows. For more sophisticated studies that involve FastJet and other tools, use C++/ROOT.</note>
 Here is an example of plotting the distribution of pseudorapidity for truth and reconstructed electrons:
@@ Line 109: / Line 109: @@
 {{:fcs:eic:mc_electrons.png?400|}}
-You can use this template to study:
+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:
   * electron efficiency as a function of pT
@@ Line 117: / 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 126: / 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>