#define process_cxx

#include "process.h"
#include <iostream>

//using namespace std;

void process::Begin(TTree * /*tree*/)
{
   TString option = GetOption();
   
   //OutputFile = new TFile("testout.root","recreate");
   
   //DeclareHistos();
}

void process::DeclareHistos()
{

   CalibratedSpec = new TH1D("calib", "Calibrated Spectrum", 11500, 0, 11500);
   CalibratedSpec->GetXaxis()->SetTitle("Energy");
   CalibratedSpec->GetYaxis()->SetTitle("Counts");

}

void process::SlaveBegin(TTree * /*tree*/)
{
   TString option = GetOption();
   
   OutputFile = new TFile("testout.root","recreate");

   DeclareHistos();
}

Bool_t process::Process(Long64_t entry)
{
   fChain->GetTree()->GetEntry(entry);

   int i;
   double calibvalue;

   for(i = 0; i< GeB_hits; i++)
   {
      if(GeB_value[i] > 0)
      {
         if(GeB_channel[i]==50)
         {
            calibvalue=(double(GeB_value[i])*double(GeB_value[i]))*G_GeBlueCore[5][0]+double(GeB_value[i])*G_GeBlueCore[5][1]+G_GeBlueCore[5][2];
            if(calibvalue>5)
            {
               CalibratedSpec->Fill(calibvalue);
            }
         }
      }
   }
   
   return kTRUE;
}


void process::SlaveTerminate()
{
	OutputFile->Write();
   //OutputFile->Close();
}

void process::Terminate()
{
	//OutputFile->Write();
	//OutputFile->Close();
}