GhostMerge.h

Go to the documentation of this file.

//
// Class to take in parallel run fit outputs and merge them  
//
//  - Writes out text and root files with overall surfaces, 
//    best fit spectra etc.
//

#ifndef GHOSTMERGE_H
#define GHOSTMERGE_H

#include "GhostInputs.h"
#include "GhostSample.h"
#include "NtupleUtils/NuCutImps.h"

#include "TObject.h"
#include "TFile.h"
#include "TH2D.h"
#include "TH1D.h"
#include "TString.h"
#include "TTree.h"
#include "TCanvas.h"
#include "TChain.h"
#include "TGraph.h"

#include "NtupleUtils/NuCutImps.h"

#include <vector>
#include <iostream>

using namespace std;

class GhostMerge : public TObject
{
  friend class GhostInputs;
  friend class GhostFakeData;
  friend class GhostFramework;
  friend class GhostSample;
 
 public:
  
  GhostMerge(GhostInputs*);
  ~GhostMerge();


  void AddFile(Char_t* filename){
    vFile.push_back(new TFile(filename,"read"));
  }
  void SetNumberInputFile(Int_t num) { fNumInputFiles = num; }
  void Pause();
  void AddSample(GhostSample* sample) { vSamples.push_back(sample); }
  void RunMerge();


 private:

  GhostInputs* gInputs;


  Int_t fNumInputFiles;
  vector<TFile*> vFile;
  vector<GhostSample*> vSamples;
  TFile* fOutputRootFile;


 //Grid variables
 Double_t fSin2MinOverall;
 Double_t fSin2MaxOverall;
 Double_t fDm2MinOverall;
 Double_t fDm2MaxOverall;

 Double_t fDm2BinWidth;
 Double_t fSin2BinWidth;
   
 Int_t fNumSin2BinsOverall;
 Int_t fNumDm2BinsOverall;
 Int_t fNumSin2Bins;
 Int_t fNumDm2Bins;


 //Best fit variables
 Double_t fMinLik;
 Int_t fMinIndexDm2;
 Int_t fMinIndexSin2;
 Double_t fMinLikSin2;
 Double_t fMinLikSin2Error;
 Double_t fMinLikDm2;
 Double_t fMinLikDm2Error;
 TGraph* gBestFit;

 //Fit Containers
 TH2D* hOverallSurface;
 TH2D* hDeltaOverallSurface;
 
 TH2D* hCCDetNormalisation;
 TH2D* hCCDetNormalisationError;
 TH2D* hCCRockNormalisation;
 TH2D* hCCRockNormalisationError;
 TH2D* hNuMuBarNormalisation;
 TH2D* hNuMuBarNormalisationError;
 TH2D* hNCNormalisation;
 TH2D* hNCNormalisationError;
 
 vector<TH2D*> v2DSystematicSurfaces;
 vector<TH2D*> v2DSystematicErrorSurfaces;
 vector<TH2D*> v1DSystematicSurfaces;
 vector<TH2D*> v1DSystematicErrorSurfaces;

 TH2D* hDecoherenceSurface;
 TH2D* hDecoherenceSurfaceError;

 ClassDef(GhostMerge,1)

};

#endif

---