fcs:eic:intro
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
fcs:eic:intro [2016/09/23 20:10] – [Reconstructed events] asc | fcs:eic:intro [2016/11/06 16:31] – [Looking at track resolution] asc | ||
---|---|---|---|
Line 77: | Line 77: | ||
====== Reconstructed events ====== | ====== Reconstructed events ====== | ||
Each detector in HepSim has " | Each detector in HepSim has " | ||
- | Click this link: [[http:// | + | Click this link: [[http:// |
- | The reconstruction tags are shown in the info page of this sample [[http:// | + | The reconstruction tags are shown in the info page of this sample [[http:// |
<code bash> | <code bash> | ||
Line 93: | Line 93: | ||
You can run such scripts using " | You can run such scripts using " | ||
- | <note tip>Look at the tutorial [[community: | + | <note tip>Look at the tutorial [[community: |
Here is an example of plotting the distribution of pseudorapidity for truth and reconstructed electrons: | Here is an example of plotting the distribution of pseudorapidity for truth and reconstructed electrons: | ||
Line 109: | Line 109: | ||
{{: | {{: | ||
- | 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 |
+ | |||
+ | ====== 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 {{: | ||
+ | for the MCparticle container inside LCIO file. We will use the "flat pT" sample called [[http:// | ||
+ | |||
+ | <code bash> | ||
+ | hs-get pgun_elec_eta4_pt30gev%rfull052 data 2 100 # get 100 files in 2 threads in the directory " | ||
+ | fpad mc_truth_flat.py | ||
+ | </ | ||
+ | |||
+ | The resulting image is: | ||
+ | |||
+ | {{: | ||
+ | |||
+ | Now, let's calculate track resolution. We will use this script {{: | ||
+ | makes fit of pT(reco)/ | ||
+ | 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 " | ||
+ | fpad track_resolution.py | ||
+ | </ | ||
+ | |||
+ | The final plot is shown below: | ||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | ====== How to continue ====== | ||
Now you can use this template to study: | Now you can use this template to study: | ||
Line 119: | Line 148: | ||
* 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 |
fcs/eic/intro.txt · Last modified: 2016/12/09 16:34 by asc