This package generates the MC simulation for the Hall A HRSs using the Single Arm version. 1. phase_space_fall16 Generates events from the spectrometer and reconstruct them back after passing through the magnets and the detectors. As it is right now, it only works for the LHRS ###Input files: input/file.inp where file means the target to be analized (d2,h3,he3): has the information of the target. input/runplan.inp : contains the kinematic information runid, E0, p0, theta, specx_off, specy_off, specz_off where: id: run number or de kinematic number. E0: the Beam Energy p0: the Electron Momentum theta: spectrometer angle specx_off: Spectrometer x offset specy_off: Spectrometer y offset specz_off: Spectrometer z offset ##Compiling: # In the phase_space_fall16 folder make clean make #### Executing: # In the phase_space_fall16 folder ./mc_hrs_single file runid tarid file: name in the input folder (d2,h3,he3¦) runid: runid number in the runplan.inp file tarid: target identification number (1=d2, 2=h3, 3=he3) #### Output i. Summary of the simulation: outfiles/t2_"runid"_"tarid".hist ii. ntuple with all the information of the sumulation worksim/t2_"runid"_"tarid".rzdat 2. Weighting Take the ntuples generated by the single arm montecarlo along with the cross section table, and use them to generate the physics and the spectrometer variables that can be compared directly with the data. ###Input files: i) input.dat Name : (useless for now) ngentot : (the total number of events generated in the MC) maxev : Maximum numbers of events to be analilzed (In the most unlikely scenario it will be up to the number of events generated) dxp : Theta reconstruction cut (mRad) dyp : Phi reconstruction cut (mRad) delup : dP/P reconstruction cut (%), upper limit deldown : dP/P reconstruction cut (%), lower limit cs_flag : Not sure what is for (not being use at the moment) rc_flag : Radiative Corrections Flag (0:no, 1:yes) ii) targetdata.dat (1) tarid (2) Z (3) Target Length (cm) (4) Targar Mass (uma) (5) A (6) Target density (gm/cm2) (7) Target Name iii) table.out."runid"."tarid" It is the output of the XEMC code. (1) A (2) Z (3) E (Beam Energy) (4) E' Scatter electron momentum (5) Scatter Angle (6) Born Crooss Section (nb/sr/GeV) (7) Radiative Crorrection Factor (8) Xbj (9) Q2 This table will be used to assign the physics weight for each point simulated by the phase space simulation. iv) The output file of the phase space simulation ~/phase_space_fall16/t2_"runid"_"tarid".rzdat ##Compiling: # In the weighting folder make clean make #### Executing: # In the weighting folder ./recon_MC runid tarid #### Output The output file is output/mc_"runid"_"tarid".rzdat To convert to root use: h2root mc_"runid"_"tarid".rzdat