//////////////////////////////////////////////////////////////////////////
//
//  Right arm detector replay script
//
//////////////////////////////////////////////////////////////////////////
//	
//	Author : Vincent A. Sulkosky & Jin Huang
//
//////////////////////////////////////////////////////////////////////////


#include "def.h"

using namespace std;

void replay_det_R(Int_t runnumber=0,Int_t all=0,Int_t fstEvt=0,Bool_t QuietRun = kFALSE)
{

  THaApparatus* B = new THaIdealBeam("B","Idea Beam, for Test Only");
  gHaApps->Add( B );
  
  //cout<<"replay: Adding R-arm Helicity"<< endl;
  B->AddDetector( new THaADCHelicity("adchel.R","Beam helicity R-arm") );
  B->AddDetector( new THaG0Helicity("g0hel.R","Right arm G0 helicity") );	
  
  cout<<"replay: Adding UnRastered and Rastered Beam ..."<<endl;
  gHaApps->Add(new THaUnRasteredBeam("Rurb","Unrastered beam"));
  gHaApps->Add(new THaRasteredBeam("Rrb","Rastered Beam"));
  
  cout<<"replay: Adding Decoder Data ..."<<endl;
   gHaApps->Add(new THaDecData("D","Misc. Decoder Data"));
  
  cout<<"replay: Adding HRS-R ..."<<endl;
  THaApparatus* HRSR = new THaHRS("R","Right HRS");
  gHaApps->Add( HRSR );
  
  // add detectors that are not in the default config
  HRSR->AddDetector( new THaCherenkov("cer", "Gas Cherenkov counter" ));
  HRSR->AddDetector( new THaShower("ps", "Pre-shower" ));
  HRSR->AddDetector( new THaShower("sh", "Shower" ));


  // add neutron detector
  THaNA* NA=new THaNA("NA","NeutronDetector"); 
  gHaApps->Add(NA);
  
/*  cout<<"replay: adding Physics modules ..."<<endl;
  
  gHaPhysics->Add( new THaGoldenTrack("R.gold","Golden track for RHRS", "R") );
  
  THaPhysicsModule *Rpt_r = new THaReactionPoint(
						 "ReactPt_R","Reaction vertex for Right","R","B");
  gHaPhysics->Add( Rpt_r );
  
  // Correct for using an Extended target
  // This needs information about the Point of interaction (hence a THaVertexModule)
  THaPhysicsModule* TgC_r = new THaExtTarCor("ExTgtCor_R",
					     "Corrected for extended target, HRS-R",
					     "R","ReactPt_R");
  gHaPhysics->Add( TgC_r );
  
  cout<<"replay: Adding Target ..."<<endl;
  THaHe3Target* pT=new THaHe3Target("he3R","Hall A Polarized he3 target for RHRS DAQ");
  gHaApps->Add(pT);
  
  // add scalers
  THaScalerGroup* RightScalers = new THaScalerGroup("Right");
  gHaScalers->Add(RightScalers);
  
  ////////////////////////////////////////
  //    Norm Ana
  ////////////////////////////////////////
  //char buff[1000]="";
  //sprintf(buff,"./summaryfiles/RNormAna%d.log",runnumber);

  //cout<<"replay: Adding NormAna & saving result to "<<buff<<endl;
  //NormAna* norm = new NormAna("N","Normalization Analysis",buff,8);
  //gHaPhysics->Add(norm);

*/
 ////////////////////////////////////////
  //   Do replay
  ////////////////////////////////////////
  THaAnalyzer* analyzer = THaAnalyzer::GetInstance();
  if( !analyzer ) {
  analyzer = new THaAnalyzer;
  }
  analyzer->EnableBenchmarks();
  analyzer->EnableHelicity();
  analyzer->EnableScalers();


  ReplayCore(
	    runnumber,            //run #
	    all,                  //-1=replay all;0=ask for a number
	    100000,                   //defaut replay event num
	    "%s/e05102_det_R_%d.root", //output file format
	    "replay_det_R.odef",	//out define
	    "replay_det_R.cdef", 	//empty cut define
	    kTRUE,                 	//replay scalar?
	    fstEvt,					//First Event To Replay
		QuietRun				//whether ask user for inputs
	    );

  cout << endl;
  cout << "Look at histograms for Right arm" << endl
       << "\t online(\"detectorR\"," << runnumber << ")" << endl;
  cout << "And Neutron Detector.. " << endl
       << "\t online(\"detectorND\"," << runnumber << ")" << endl << endl;

  cout << "Exit ROOT session before performing another analysis:" << endl
       << "  .q" << endl << endl;

}