Design Systematic uncertanties code

[jtenavidal]

What we need:

• Plotting code: Must be able to add systematics on top of statistical error
• Effects of fiducial on subtracted particles. Acceptance should compensate for signal particle fiducial effects. Changing 3 degrees each side (L and R). Must re-calculate with +shift and -shift. Need to understand where to modify the sub-routine in the code. Add a configurable parameter (angle shift) and pass it as an option to the fiducial routine in Substraction only. Must be able to turn off and on the cut on the app, and change sign. In configuration we only have the absolute value. Add files in automatic script calculation.

Can add shift here: https://github.com/e4nu/e4nuanalysiscode/blob/bddb3402c61d30c71b25982a4a8b9e1b4bfed20d/src/analysis/BackgroundI.h#L39

GetFiducialCut(Optionalshift)

• Sector to sector : we already have the plots we only need to take the difference between the sectors. Take maximum difference as uncertainty?

• change beta cut in GetCharge app.

• Model uncertainty on acceptance correction. Modelling FSI won't change much. We can change overall RES vs DIS scaling. Reweight?

https://github.com/e4nu/e4nuanalysiscode/blob/bddb3402c61d30c71b25982a4a8b9e1b4bfed20d/src/analysis/BackgroundI.h#L152 Could have a set event error:

temp_event.SetEventError( (pow(temp_event.GetEventError(),2)+pow(rot_error,2)  ) ; 

• Uncertainty due to number of rotations. Can I return also the error?

[jtenavidal]

Developing in develop/Systematics

[jtenavidal]

Fiducial cut:

Bool_t Fiducial::FiducialCut( const int pdg, const double beam_en, TVector3 momentum, const bool is_data, const bool apply_fiducial ) {