- THnSparse *fOriginalPrior; // This is the original prior distribution : will not be modified
- THnSparse *fInverseResponse; // Inverse response matrix
- THnSparse *fMeasuredEstimate; // Estimation of the measured (M) spectrum given the a priori (T) distribution
- THnSparse *fConditional; // Matrix holding the conditional probabilities P(M|T)
- THnSparse *fProjResponseInT; // Projection of the response matrix on TRUE axis
- THnSparse *fUnfolded; // Unfolded spectrum
- Int_t *fCoordinates2N; // Coordinates in 2N (measured,true) space
- Int_t *fCoordinatesN_M; // Coordinates in measured space
- Int_t *fCoordinatesN_T; // Coordinates in true space
-
+ //
+ Double_t fMaxConvergence; // Convergence criterion in case of correlated error calculation
+ Int_t fNRandomIterations; // Number of random distributed measured spectra
+ THnSparse *fResponse; // Copy of the original response matrix (modified)
+ THnSparse *fPrior; // Copy of the original prior spectrum (modified)
+ THnSparse *fEfficiency; // Copy of original efficiency (modified)
+ THnSparse *fMeasured; // Copy of the original measureed spectrum (modified)
+ THnSparse *fInverseResponse; // Inverse response matrix
+ THnSparse *fMeasuredEstimate; // Estimation of the measured (M) spectrum given the a priori (T) distribution
+ THnSparse *fConditional; // Matrix holding the conditional probabilities P(M|T)
+ THnSparse *fUnfolded; // Unfolded spectrum (modified before and during error calculation)
+ THnSparse *fUnfoldedFinal; // Final unfolded spectrum
+ Int_t *fCoordinates2N; // Coordinates in 2N (measured,true) space
+ Int_t *fCoordinatesN_M; // Coordinates in measured space
+ Int_t *fCoordinatesN_T; // Coordinates in true space
+
+
+ /* correlated error calculation */
+ THnSparse *fRandomResponse; // Randomized distribution for each bin of the response matrix to calculate correlated errors
+ THnSparse *fRandomEfficiency; // Randomized distribution for each bin of the efficiency spectrum to calculate correlated errors
+ THnSparse *fRandomMeasured; // Randomized distribution for each bin of the measured spectrum to calculate correlated errors
+ TRandom3 *fRandom3; // Object to get random number following Poisson distribution
+ THnSparse *fDeltaUnfoldedP; // Profile of the delta-unfolded distribution
+ THnSparse *fDeltaUnfoldedN; // Entries of the delta-unfolded distribution (count for each bin)
+ Short_t fNCalcCorrErrors; // Book-keeping to prevend infinite loop
+ UInt_t fRandomSeed; // Random seed
+