Difference between revisions of "Setup for Batch Job With Varying Experimental CLAS12 Quantities"
(One intermediate revision by the same user not shown) | |||
Line 2: | Line 2: | ||
<pre> | <pre> | ||
− | + | mv LUND_Spread_LH2.LUND LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut.LUND</pre> | |
− | Spliting the LUND file into 38 parts of | + | Spliting the LUND file into 38 parts of 1000 events, or 3000 lines |
<pre> | <pre> | ||
− | split -d -l | + | split -d -l 3000 -a 3 LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut.LUND LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut_ |
</pre> | </pre> | ||
Line 13: | Line 13: | ||
<pre> | <pre> | ||
− | prename 's/( | + | prename 's/(LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut_\d{3})/$1.LUND/' LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut_* |
</pre> | </pre> | ||
− | Placing each of these files into its own directory, within a directory named | + | Placing each of these files into its own directory, within a directory named LH2_0Sol_p100Tor_11GeV_Phi0deg_ShieldOut |
<pre>find . -name "*.LUND" -exec sh -c 'mkdir "${1%.*}" ; mv "$1" "${1%.*}" ' _ {} \;</pre> | <pre>find . -name "*.LUND" -exec sh -c 'mkdir "${1%.*}" ; mv "$1" "${1%.*}" ' _ {} \;</pre> | ||
Line 43: | Line 43: | ||
The weighting found in the LUND file is applied to the scattering angle theta in the Center of Mass frame to produce the Moller Differential Cross Section. We should be able to use the information from the detector hits such as wire numbers and sectors to find a distribution of hits with respect to the scattering angle theta in the Lab Frame. The Theta distribution will need to undergo a Lorentz transformation to the Center of Mass Frame to apply the corresponding weight found in the LUND file. | The weighting found in the LUND file is applied to the scattering angle theta in the Center of Mass frame to produce the Moller Differential Cross Section. We should be able to use the information from the detector hits such as wire numbers and sectors to find a distribution of hits with respect to the scattering angle theta in the Lab Frame. The Theta distribution will need to undergo a Lorentz transformation to the Center of Mass Frame to apply the corresponding weight found in the LUND file. | ||
+ | |||
+ | |||
+ | =Plotting detector grids= | ||
+ | |||
+ | Using ~/src/CLAS/coatjava-2.4/bin/evio-dump -i LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_00/LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_00.evio, we can test for the first particle hits. We show 9 hits, all on wire 1. | ||
+ | <pre> | ||
+ | *********************** EVENT # 2 *********************** | ||
+ | |||
+ | +------------------------------------------------+------------+------------+ | ||
+ | | bank| nrows| ncols| | ||
+ | +------------------------------------------------+------------+------------+ | ||
+ | | DC::dgtz| 9| 11| | ||
+ | | DC::true| 9| 24| | ||
+ | | GenPart::true| 2| 7| | ||
+ | | HTCC::dgtz| 2| 6| | ||
+ | | PCAL::dgtz| 7| 7| | ||
+ | +------------------------------------------------+------------+------------+ | ||
+ | |||
+ | |||
+ | Press Enter for Next Event or Bank Name: DC::dgtz | ||
+ | *****>>>>> BANK DC::dgtz >>>> SIZE = 11 | ||
+ | sector (int) : 1 1 1 1 1 1 1 1 1 | ||
+ | superlayer (int) : 1 1 1 1 2 2 2 2 2 | ||
+ | layer (int) : 1 2 3 4 1 2 3 4 6 | ||
+ | wire (int) : 1 1 1 1 1 1 1 1 1 | ||
+ | LR (int) : -1 -1 -1 -1 1 1 -1 -1 -1 | ||
+ | |||
+ | </pre> | ||
+ | |||
+ | |||
+ | This is confirmed by the plot: | ||
+ | |||
+ | [[File:Region1DChits.png]] |
Latest revision as of 18:29, 27 December 2016
Using CLAS12 isotropic distribution LUND files for LH2 and NH3 targets and varying incident electron energy.
mv LUND_Spread_LH2.LUND LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut.LUND
Spliting the LUND file into 38 parts of 1000 events, or 3000 lines
split -d -l 3000 -a 3 LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut.LUND LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut_
Naming these files LUND files
prename 's/(LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut_\d{3})/$1.LUND/' LH2_0Sol_n100Tor_11GeV_Phi0deg_ShieldOut_*
Placing each of these files into its own directory, within a directory named LH2_0Sol_p100Tor_11GeV_Phi0deg_ShieldOut
find . -name "*.LUND" -exec sh -c 'mkdir "${1%.*}" ; mv "$1" "${1%.*}" ' _ {} \;
GCard Building
A program,File:Gcard builder.C is written to install a LUND specific gcard in each directory.
Commands Building
Similarly, File:Commands builder.C is written to install a specific commands file into each directory for testing.
Batch Conditions
The batch job will need to
cd /home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut cd /LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_00 ~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_00.LUND" -N=3000 eg12.gcard>/dev/null - - - cd /LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_37 ~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_37.LUND" -N=3000 eg12.gcard>/dev/null cd .. ~/src/CLAS/evioRootDV/bin/evio2root LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut.root 114000
The weighting found in the LUND file is applied to the scattering angle theta in the Center of Mass frame to produce the Moller Differential Cross Section. We should be able to use the information from the detector hits such as wire numbers and sectors to find a distribution of hits with respect to the scattering angle theta in the Lab Frame. The Theta distribution will need to undergo a Lorentz transformation to the Center of Mass Frame to apply the corresponding weight found in the LUND file.
Plotting detector grids
Using ~/src/CLAS/coatjava-2.4/bin/evio-dump -i LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_00/LH2_0Sol_n100Tor_11GeV_Phi10deg_ShieldOut_00.evio, we can test for the first particle hits. We show 9 hits, all on wire 1.
*********************** EVENT # 2 *********************** +------------------------------------------------+------------+------------+ | bank| nrows| ncols| +------------------------------------------------+------------+------------+ | DC::dgtz| 9| 11| | DC::true| 9| 24| | GenPart::true| 2| 7| | HTCC::dgtz| 2| 6| | PCAL::dgtz| 7| 7| +------------------------------------------------+------------+------------+ Press Enter for Next Event or Bank Name: DC::dgtz *****>>>>> BANK DC::dgtz >>>> SIZE = 11 sector (int) : 1 1 1 1 1 1 1 1 1 superlayer (int) : 1 1 1 1 2 2 2 2 2 layer (int) : 1 2 3 4 1 2 3 4 6 wire (int) : 1 1 1 1 1 1 1 1 1 LR (int) : -1 -1 -1 -1 1 1 -1 -1 -1
This is confirmed by the plot: