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


#include "def.h"

using namespace std;

void replay_det_L(Int_t runnumber=0,Int_t all=-1,Int_t fstEvt=0,Bool_t QuietRun = kTRUE)
{

    THaApparatus* B = new THaIdealBeam("B","Idea Beam, for Test Only");
    gHaApps->Add( B );
  
    cout<<"replay: Adding UnRastered and Rastered Beam ..."<<endl;
    gHaApps->Add(new THaUnRasteredBeam("urb","Unrastered beam"));
    gHaApps->Add(new THaRasteredBeam("rb","Rastered Beam"));
  
    cout<<"replay: Adding Decoder Data ..."<<endl;
    gHaApps->Add(new THaDecData("DL","Misc. Decoder Data"));
  
    cout<<"replay: Adding HRS-L ..."<<endl;
    THaApparatus* HRSL = new THaHRS("L","Left HRS");
    gHaApps->Add( HRSL );
  
  // add detectors that are not in the default config
    HRSL->AddDetector( new THaCherenkov("cer", "Gas Cherenkov counter" ));
    HRSL->AddDetector( new THaCherenkov("a1", "A1 Cherenkov counter" ));
    HRSL->AddDetector( new THaShower("prl1", "Pre-shower pion rej." ));
    HRSL->AddDetector( new THaShower("prl2", "Show pion rej." ));

    cout<<"replay: adding Physics modules ..."<<endl;
    gHaPhysics->Add( new THaGoldenTrack("L.gold","Golden track for LHRS", "L") );
  
	THaPhysicsModule *Rpt_l = new THaReactionPoint(
			"ReactPt_L","Reaction vertex for Left","L","rb");
	gHaPhysics->Add( Rpt_l );
  
	THaPhysicsModule *Rpt_l = new THaReactionPoint(
			"ReactPt_Lurb","Reaction vertex for Left","L","urb");
	gHaPhysics->Add( Rpt_l );
  
  // Correct for using an Extended target
  // This needs information about the Point of interaction (hence a THaVertexModule)
	THaPhysicsModule* TgC_l = new THaExtTarCor("ExTgtCor_L",
											   "Corrected for extended target, HRS-L",
			  "L","ReactPt_L");
	gHaPhysics->Add( TgC_l );
	THaPhysicsModule* TgC_l = new THaExtTarCor("ExTgtCor_Lurb",
											   "Corrected for extended target, HRS-L",
			  "L","ReactPt_Lurb");
	gHaPhysics->Add( TgC_l );

	cout<<"replay: adding Physics modules ..."<<endl;
	THaPrimaryKine *PriKine=new THaPrimaryKine("PriKine","kinematics of scattering of electron to BB","ExTgtCor_L","rb",.938272013);
	gHaPhysics->Add(PriKine);
  
    
	cout<<"replay: adding Physics modules ..."<<endl;
	THaPrimaryKine *PriKine=new THaPrimaryKine("PriKine_C","kinematics of scattering of electron to BB","ExTgtCor_L","rb",12*.931494028-.000511*6);
	gHaPhysics->Add(PriKine);
    
	cout<<"replay: adding Physics modules ..."<<endl;
	THaPrimaryKine *PriKine=new THaPrimaryKine("PriKine_He3","kinematics of scattering of electron to BB","ExTgtCor_L","rb",3.0160293*.931494028-.000511*2);
	gHaPhysics->Add(PriKine);
  
  
    ReplayCore(
               runnumber,            //run #
               all,                  //-1=replay all;0=ask for a number
               -1,                   //defaut replay event num
               "%s/e06010_det_L_%d.root", //output file format
               "replay_det_L.odef",	//out define
               "replay_det_L.cdef", 	//empty cut define
               kFALSE,                 	//replay scalar?
               fstEvt,					//First Event To Replay
               QuietRun				//whether ask user for inputs
              );
}