QSigExChecks

class QSigExChecks

    public:

                      QSigExChecks(TFile* aFile = NULL)
                      QSigExChecks(const QSigExChecks&)
              virtual ~QSigExChecks()
                 void AddEvents(TFile* aFile = NULL, Char_t* aDir = "NULL", Char_t* aTreeName = "NULL")
                 void AddNCDEvents(TFile* aFile = NULL)
                 void AddPMTEvents(TFile* aFile = NULL)
               Bool_t CheckEigenValues(TMatrixDSym* aMatrix = NULL)
               Bool_t CheckSum(TMatrixDSym* aMatrix = NULL, Int_t aNoEvents = 0)
       static TClass* Class()
         TMatrixDSym* GetCovariance(TFile* aFile = NULL, Char_t* aDir = "NULL", Char_t* aName = "NULL")
                 void GetF2EMaps(TFile* aFile = NULL, Char_t* aDir = "NULL")
                 void GetFitResult(TFile* aFile = NULL, Char_t* aDir = "NULL")
                 void GetNCDF2EMaps(TFile* aFile = NULL)
                Int_t GetNoEvents()
                Int_t GetParActive(Int_t i)
             Double_t GetParameter(Int_t i)
                 void GetParameters(TFile* aFile = NULL, Char_t* aDir = "NULL")
             Double_t GetParMinusError(Int_t i)
             Double_t GetParPlusError(Int_t i)
                 void GetPMTF2EMaps(TFile* aFile = NULL)
             Double_t GetSum(TMatrixDSym* aMatrix = NULL)
      virtual TClass* IsA() const
        QSigExChecks& operator=(const QSigExChecks&)
             Double_t PlotLikelihoodFunction(Int_t ipar, Int_t jpar, Char_t* outbase = "NULL")
             Double_t PlotLikelihoodFunction(Int_t ipar, Char_t* outbase = "NULL")
                 void PrintCorrelationMatrix()
                 void PrintEventNumbers()
                 void PrintFlux2Events()
                 void PrintParameters()
         virtual void ShowMembers(TMemberInspector& insp, char* parent)
         virtual void Streamer(TBuffer& b)
                 void StreamerNVirtual(TBuffer& b)

Data Members

    private:

         static const Int_t fMaxFitPars           this is for the quantile calculation
      static const Double_t fQuantileProbability  
      static const Double_t fQuantile             
                     TFile* fFile                 
                    TString fPMTEventDir          
                    TString fPMTF2EDir            
                    TString fNCDEventDir          
                    TString fNCDF2EDir            
                    TString fFitNumbersDir        
                    TString fFitMinuitDir         
                    TString fFitParametersDir     
                    TString fCovarianceName       
                    TString fStatusName           
                   Long64_t fNoEvents             total number of events in data sample
               TMatrixDSym* fCovariance           
                      Int_t fnFluxes              
                      Int_t fnActiveFluxes        
                   TString* fFluxName             
                      Int_t fnF2EMaps             
                QSigExDis** fF2EMap               !
                  TString** fF2EFluxName          !
                     Int_t* fF2EFluxIndex         !
                  TString** fF2EEvtName           !
                      Int_t fnFitPars             
                  TString** fFitParName           !
                     Int_t* fFitParActive         !
                     Int_t* fFitParIndexOfActive  ! the index of the parameter as function of the active index 
                  Double_t* fFitPar               !
                  Double_t* fFitParMErr           !
                  Double_t* fFitParPErr           !
                  Double_t* fFitParErr            !

 this class implements a number of sanity checks to verify    
 a successful fit. The most useful check is the comparison of 
 the sum of the covariance matrix with the total number of    
 events in the input data.                                    
 the class interfaces with the Fit result via the output file //  
 from the main qsigex classes.                                

 Clean up memory

  cout << "<QSigExChecks::GetParameters>" << endl;

 read the flux 2 event mapping functions

 check whether the eigenvalues are all positive as they should be for a covariance matrix
 return true if all are positive (>=0), false if at least one is negative.

 plot and print the likelihood function for one parameter.
 presently this is done for one parameter

 plot and print the likelihood function for one parameter.
 presently this is done for one parameter