Difference between revisions of "Using SIMC"

From Hall A Wiki
Jump to: navigation, search
 
(2 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 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.rzdat ''' will now be generated. 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.
+
* 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")



[back He3/H3 e,e'p experiment]

[back to Hall A wiki]