Difference between revisions of "Using SIMC"
From Hall A Wiki
(4 intermediate revisions by the same user not shown) | |||
Line 1: | Line 1: | ||
== SIMC repository and documentation == | == SIMC repository and documentation == | ||
− | * [[SIMC short documentation]] | + | * [[https://hallaweb.jlab.org/wiki/images/0/07/Simc.pdf SIMC short documentation]] |
* [https://github.com/JeffersonLab/simc_gfortran SIMC on github] | * [https://github.com/JeffersonLab/simc_gfortran SIMC on github] | ||
− | |||
== Brief Overview == | == Brief Overview == | ||
Line 10: | Line 9: | ||
* SIMC main inputs are | * SIMC main inputs are | ||
− | Beam energy , Spectrometers angles , Determination of which spectrometer detects electrons and which detects hadrons , Target parameters (A,Z,density etc.) , Generator acceptance and resolution | + | Beam energy , Spectrometers angles , Determination of which spectrometer detects electrons and which detects hadrons , |
+ | Target parameters (A,Z,density etc.) , Generator acceptance and resolution | ||
* The generator will generate events with the reaction 4 physical variables | * The generator will generate events with the reaction 4 physical variables | ||
in-plane angle '''y'''' , out-of-plane angle '''x'''' , '''y'''(target) , momentum difference from central spectrometer momentum - '''delta''' | in-plane angle '''y'''' , out-of-plane angle '''x'''' , '''y'''(target) , momentum difference from central spectrometer momentum - '''delta''' | ||
− | |||
* COSY simulation tool will generate the '''hrs_cosy''' file which include the transportation parameters for the particle along the spectrometer | * COSY simulation tool will generate the '''hrs_cosy''' file which include the transportation parameters for the particle along the spectrometer | ||
Line 44: | Line 43: | ||
worksim/hrsr_eep_h.rzdat (paw ntuple) | worksim/hrsr_eep_h.rzdat (paw ntuple) | ||
− | * convert the paw file into root | + | * convert the paw file into root |
+ | cd worksim/ | ||
h2root hrsr_eep_h.rzdat | h2root hrsr_eep_h.rzdat | ||
− | * a root file '''hrsr_eep_h. | + | * a root file '''hrsr_eep_h.root ''' will now be generated in '''/worksim/'''. |
+ | |||
+ | * The tree is called '''h666''', and the most important variables are the reconstructed 4 variables (hs-electron arm, ss-hadron arm) and those generates (add an i), e.g. | ||
hsdelta , hsdeltai ,hsytar , hsytari ,hsyptar , hsyptari , hsxptar , hsxptari | hsdelta , hsdeltai ,hsytar , hsytari ,hsyptar , hsyptari , hsxptar , hsxptari | ||
+ | |||
+ | * To draw, e.g. the difference between the generated and reconstructed momentum, use | ||
+ | |||
+ | h666 -> Draw("hsdelta - hsdeltai") | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | [[https://hallaweb.jlab.org/wiki/index.php/He3_/_H3_e%2Ce%27p_experiment_E12-13-012 back He3/H3 e,e'p experiment]] | ||
+ | |||
+ | [[http://hallaweb.jlab.org/wiki back to Hall A wiki]] |
Latest revision as of 14:34, 25 February 2015
SIMC repository and documentation
Brief Overview
- SIMC main inputs are
Beam energy , Spectrometers angles , Determination of which spectrometer detects electrons and which detects hadrons , Target parameters (A,Z,density etc.) , Generator acceptance and resolution
- The generator will generate events with the reaction 4 physical variables
in-plane angle y' , out-of-plane angle x' , y(target) , momentum difference from central spectrometer momentum - delta
- COSY simulation tool will generate the hrs_cosy file which include the transportation parameters for the particle along the spectrometer
- The output of the simulation would be the 4 physical variables reconstructed at the focal plane of each spectrometer
How to run SIMC
- compile SIMC directory: go to the SIMC directory, change what ever necessary in Makefile and type
make
- create an input file under the directory infiles following the repository examples, e.g.
infiles/hrsr_eep_h.inp
- execute SIMC with an input file, e.g.
simc (enter input file)> hrsr_eep_h
- after the run was completed 4 output files are generated
outfiles/hrsr_eep_h.geni (information) outfiles/hrsr_eep_h.hist (run settings) outfiles/hrsr_eep_h.gen (output text file) worksim/hrsr_eep_h.rzdat (paw ntuple)
- convert the paw file into root
cd worksim/ h2root hrsr_eep_h.rzdat
- a root file hrsr_eep_h.root will now be generated in /worksim/.
- The tree is called h666, and the most important variables are the reconstructed 4 variables (hs-electron arm, ss-hadron arm) and those generates (add an i), e.g.
hsdelta , hsdeltai ,hsytar , hsytari ,hsyptar , hsyptari , hsxptar , hsxptari
- To draw, e.g. the difference between the generated and reconstructed momentum, use
h666 -> Draw("hsdelta - hsdeltai")