Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Verify that systematic universes utilize same random numbers between ICARUS and SBND #186

Open
jzennamo opened this issue Nov 8, 2021 · 2 comments
Assignees

Comments

@jzennamo
Copy link
Contributor

jzennamo commented Nov 8, 2021

When drawing a systematic universe in either the flux or cross-section uncertainties we would want them be correlated throws between the near and far detector. We can check this by verifying that the random numbers used for the systematic weight calculations are the same for the near and far detector for the same systematic universe.

e.g. the pion production flux uncertainty will create N universes, universe k, should use the same random numbers for ICARUS and SBND

@cjbacchus
Copy link
Contributor

One easy way to check this would be to dump the universe configurations from the resulting CAF files and check they are identical. The information is stored within the globalTree tree, but I don't think there is a simple tool to dump it.

@Ooohu
Copy link
Contributor

Ooohu commented Nov 24, 2021

  • [ x ] Using cafmaker* in sbncode
  • [ x ] Using cafmaker* in icaruscode
    Both jobs give the same random number sequences as the following:
SBNEventWeight Flux: configure horncurrent
SBNEventWeight Flux: `parameter_sigma` was not set; now it is set as 1
SBNEventWeight Flux: Calculator type: Unisim
Check 0th Random # -0.837599
Check 1th Random # 0.874643
Check 2th Random # -0.776626
Check 3th Random # -1.0661
Check 4th Random # -1.13575
Check 5th Random # 1.48533
Check 6th Random # -1.38929
Check 7th Random # -1.38041
Check 8th Random # 1.50705
Check 9th Random # -1.35532
Check 10th Random # -1.214
Check 11th Random # 1.58389
Check 12th Random # 0.672717
Check 13th Random # -0.667711
Check 14th Random # 1.36864
Check 15th Random # 0.0516781
Check 16th Random # 0.486945
Check 17th Random # 1.98229
Check 18th Random # 0.118562
Check 19th Random # 0.219084
...
...
GENIE weight calculator genie_sbnd_multisim_fake will generate 20 multisim universes
%MSG
CHECK  mode·
Check 0th Random # 1.24676
Check 1th Random # -0.893461
Check 2th Random # -2.36704
Check 3th Random # 0.758384
Check 4th Random # -0.657734
Check 5th Random # 1.67073
Check 6th Random # -0.229847
Check 7th Random # 0.05604
Check 8th Random # 0.938183
Check 9th Random # 0.896552
Check 10th Random # -0.659323
Check 11th Random # -0.418972
Check 12th Random # 0.785706
Check 13th Random # 0.258532
Check 14th Random # -0.224273
Check 15th Random # -1.1992
Check 16th Random # 1.87169
Check 17th Random # 2.27482
Check 18th Random # 0.061628
Check 19th Random # -2.09423
...
Check 0th Random # 0.782663
Check 1th Random # 2.07384
Check 2th Random # -1.47801
Check 3th Random # 0.127152
Check 4th Random # -0.607407
Check 5th Random # -0.274356
Check 6th Random # 0.271645
Check 7th Random # 1.25048
Check 8th Random # 0.136712
Check 9th Random # -0.49169
Check 10th Random # 1.34227
Check 11th Random # -0.179574
Check 12th Random # -0.166795
Check 13th Random # 0.0767447
Check 14th Random # -0.664041
Check 15th Random # 0.265925
Check 16th Random # -1.57082
Check 17th Random # -1.99057
Check 18th Random # -0.0507365
Check 19th Random # -1.27382

Here I use random seed 7 for first 20 universes multisim for one flux uncertainty; random seed 100 for first 20 universes for 2 genie uncertainties. They all use the same FHiCL listed under sbncode/sbncode/SBNEventWeight/jobs

@Ooohu Ooohu closed this as completed Nov 24, 2021
@Ooohu Ooohu reopened this Nov 24, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants