Difference between revisions of "Using SIMC"

From Hall A Wiki
Jump to: navigation, search
(SIMC repository and documentation)
 
Line 57: Line 57:
  
 
   h666 -> Draw("hsdelta - hsdeltai")
 
   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 13: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]