CorrelationFactory.cxx

#include <fstream>
#include <iostream>
#include "MagicDisplayConventions.h"
#include "CorrelationFactory.h"
#include "AnitaConventions.h"
#include "AnitaGeomTool.h"
#include "WaveformGraph.h"
#include "FFTtools.h"

#include "TString.h"
#include "TObjArray.h"
#include "TObjString.h"
#include "TVector3.h"
#include "TROOT.h"
#include "TPaveText.h"
#include "TPad.h"
#include "TText.h"
#include "TLatex.h"
#include "TGraph.h"
#include "WaveformGraph.h"
#include "TStyle.h"
#include "TCanvas.h"
#include "TAxis.h"
#include "TH1.h"
#include "THStack.h"
#include "TList.h"
#include "TObject.h"
#include "TTimeStamp.h"



CorrelationFactory*  CorrelationFactory::fgInstance = 0;


CorrelationFactory::CorrelationFactory()
{
  //Default constructor
  fDeltaT=1./2.6;  
  fgInstance=this;
  fGrWave[0]=NULL;
  fGrWave[1]=NULL;
  
}

CorrelationFactory::~CorrelationFactory()
{
  //Default destructor
}



//______________________________________________________________________________
CorrelationFactory*  CorrelationFactory::Instance()
{
  //static function
  return (fgInstance) ? (CorrelationFactory*) fgInstance : new CorrelationFactory();
}

void CorrelationFactory::addWaveformToCorrelation(WaveformGraph *inputWave)
{
  char graphTitle[180];
  if(fGrWave[0]==NULL)
    fGrWave[0]=inputWave;
  else {
    fGrWave[1]=inputWave;
    //Now we can correlate
    TGraph *grCor = FFTtools::getInterpolatedCorrelationGraph(fGrWave[0],fGrWave[1],fDeltaT);
    if(fGrWave[0]->fChan!=8) {
      sprintf(graphTitle,"Correlation Ant %d %c -%d %c",fGrWave[0]->fAnt,
          AnitaPol::polAsChar(fGrWave[0]->fPol),fGrWave[1]->fAnt,
          AnitaPol::polAsChar(fGrWave[1]->fPol));
    }
    else {
      sprintf(graphTitle,"Correlation SURF %d Chan %d - SURF %d Chan %d",
          fGrWave[0]->fSurf+1,fGrWave[0]->fChan+1,
          fGrWave[1]->fSurf+1,fGrWave[1]->fChan+1);
    }    
    new TCanvas();
    grCor->SetTitle(graphTitle);
    grCor->Draw("al");
    fGrWave[0]=0;
    fGrWave[1]=0;
  }

}