-
Notifications
You must be signed in to change notification settings - Fork 35
Configurations Enabling ICARUS NuMI Analysis #391
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
Conversation
Co-authored-by: Gianluca Petrillo <petrillo@slac.stanford.edu>
SBNEventWeight differences. Also make some adjustments to get it to build against geant4reweight v01_00_03e -q e20:s117:prof
* Placeholder for sampled parameter storage, for CAF compatibility. * Renamed fcl for consistency
…calorimetry variables
…e/gputnam-configs
… feature/gputnam-configs
jzennamo
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good, I might base your PR off my other PR
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This change is probably best left for this PR:
#394
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This change is probably best left for this PR:
#394
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This change is probably best left for this PR:
#394
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This change is probably best left for this PR:
#394
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This looks good, I might add a portion which corrects for the minbias gates directly to this accounting, something like this:
if(evtCount != 0 && totalMinBias != 0)
scale_factor = 1. - (evtCount/totalMinBias);
from https://github.com/SBNSoftware/sbncode/pull/394/files
| //Get the weight from the G4ReweightTraj | ||
| w = theReweighter->GetWeight( &theTraj ); | ||
| // Total weight is the product of track weights in the event | ||
| weight[j] *= std::max((float)0.0, w); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
should this be the max between 1 and w?
|
|
||
| // Do the same for elastic weight (should be 1 unless set to non-nominal ) | ||
| el_w = theReweighter->GetElasticWeight( &theTraj ); | ||
| weight[j] *= std::max((float)0.0,el_w); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
should this be the max between 1 and w?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Only comment is that maybe the random seeds should be a bit more....random ;)
| Name: "fPiPlusReacLow" | ||
| Cut: "reac" | ||
| Range: [10., 200.] | ||
| Nominal: 1 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
| Nominal: 1 | |
| Nominal: 1.0 |
| Name: "fProtonReac" | ||
| Cut: "reac" | ||
| Range: [10., 2005.] | ||
| Nominal: 1 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
| Nominal: 1 | |
| Nominal: 1.0 |
|
This PR has now been superseded by a combination of other PRs, most recently: #417, https://github.com/SBNSoftware/sbncode/pull/418/files, and #419. Closing this in favor of those. |
Gathers together a number of updates:
Depends on SBNSoftware/sbnanaobj#104.
Marking as a draft, for now.