User Tools

Site Tools


fcs:eic:intro

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Last revision Both sides next revision
fcs:eic:intro [2016/11/04 19:50]
asc [Reconstructed events]
fcs:eic:intro [2016/12/09 16:33]
asc [Hit timing information]
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). 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 ======+ 
 +===== 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|}} 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" 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 fpad mc_truth_flat.py
Line 123: Line 135:
 The resulting image is: 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 ====== ====== How to continue ======
Line 135: Line 160:
   * fake rate (cannot be done with this sample!)   * fake rate (cannot be done with this sample!)
   * calculate Q2 and x_bj, check efficiency and resolutions for these variable   * calculate Q2 and x_bj, check efficiency and resolutions for these variable
-===== Displaying reconstructed events =====+====== Displaying reconstructed events ======
  
  
Line 144: Line 169:
  
 <code bash> <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> </code>
  
fcs/eic/intro.txt · Last modified: 2016/12/09 16:34 by asc