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/06 16:26]
asc [Looking at track resolution]
fcs:eic:intro [2016/12/09 16:33]
asc [Hit timing information]
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). 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 ====== ====== 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. 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":+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> <code bash>
Line 126: Line 138:
  
 Now, let's calculate track resolution. We will use this script {{:fcs:eic:track_resolution.py|track_resolution.py}} which 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): +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> <code bash>
 hs-get pgun_elec_eta4_pt30gev%rfull052 data  2 20 # get 20 files in 2 threads in the directory "data" hs-get pgun_elec_eta4_pt30gev%rfull052 data  2 20 # get 20 files in 2 threads in the directory "data"
Line 133: Line 145:
 </code> </code>
  
 +The final plot is shown below:
  
 {{:fcs:eic:pt_resolution.png?400|}} {{:fcs:eic:pt_resolution.png?400|}}
Line 156: 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