-
Notifications
You must be signed in to change notification settings - Fork 1
Data collection
Collecting the data is a multi-step process.
Lassen example:
cd /usr/workspace/coldqcd/c51/x_files/project_2/production/a12m310_a
python ../nucleon_elastic_FF/scripts/get_res_phi.py 300 5560 5
python ../nucleon_elastic_FF/scripts/avg_mres.py --cfgs 300 5560 5 --srcs srcs0-7
Individual data files from each config get placed here by the get script
/p/gpfs1/walkloud/c51/x_files/project_2/production/a12m310_a/data/a12m310_a_cfg.h5
The avg script then places the src averaged and concatenated file here
/p/gpfs1/walkloud/c51/x_files/project_2/production/a12m310_a/data/avg/a12m310_a_avg.h5
Lassen example:
cd /usr/workspace/coldqcd/c51/x_files/project_2/production/a12m310_a
python ../nucleon_elastic_FF/scripts/get_spec.py 300 5560 5
python ../nucleon_elastic_FF/scripts/avg_spec.py --cfgs 300 5560 5
Individual data files from each config get placed here by the get script
/p/gpfs1/walkloud/c51/x_files/project_2/production/a12m310_a/data/a12m310_a_cfg.h5
The avg script then places the src averaged and concatenated file here
/p/gpfs1/walkloud/c51/x_files/project_2/production/a12m310_a/data/avg/a12m310_a_avg.h5
Lassen example
cd /usr/workspace/coldqcd/c51/x_files/project_2/production/a12m310_a
python ../nucleon_elastic_FF/scripts/get_formfac.py 300 5560 5
python ../nucleon_elastic_FF/scripts/avg_formfac.py --cfgs 300 5560 5
These routines only collect the charges or zero-momentum current insertion results. Individual data files from each config get placed here by the get script
/p/gpfs1/walkloud/c51/x_files/project_2/production/a12m310_a/data/a12m310_a_cfg.h5
The avg script then places the src averaged and concatenated file here
/p/gpfs1/walkloud/c51/x_files/project_2/production/a12m310_a/data/avg/a12m310_a_avg.h5
spec_4D happens in two steps, tslice and then avg. Here is a Lassen example:
cd /usr/workspace/coldqcd/c51/x_files/project_2/production/a12m310_a
python ../nucleon_elastic_FF/scripts/tslice_4D.py spec --cfgs 300 5295 5 -t 0.5
This will take a few hours to a day. After the slicing is done, then, we average over sources:
python ../nucleon_elastic_FF/scripts/avg_4D.py spec --cfgs 300 5295 5
and finally, concatenate them all together to the final 4D file.
python ../nucleon_elastic_FF/scripts/concatenate_spec_4D_tslice.py 300 5560 5
This happens in three steps as well, tslice and then avg and then concatenate. The tslice step usually takes so long, I split it into various sets via screen.
cd /usr/workspace/coldqcd/c51/x_files/project_2/production/a12m310_a
screen -S a12m310_a_ff_1000
python ../nucleon_elastic_FF/scripts/tslice_4D.py formfac --cfgs 1000 1995 5
repeat for 2000, 3000, 4000, 5000 and 300. These will take about a day each. After they are finished, in the same screen, I then do
python ../nucleon_elastic_FF/scripts/avg_4D.py formfac --cfgs 1000 1995 5
This will take a little less time. Finally, concatenate all the 4D files together.
python ../nucleon_elastic_FF/scripts/concatenate_formfac_4D_tslice.py 300 5295 5