Add getting beam energy from gpbeam parameter

Read gpbeam in THcRaster::ReadDatabase
Add jan05.kinematics
Modified DBASE/raster_test.database to fill g_ctp_kinematics_filename
with DBASE/jan05.kinematics
Modified raster_test.C to load g_ctp_kinematics_filename

examples/raster_test.C

@@ -18,6 +18,7 @@
   
   // Open and load parameter files
   //
+  gHcParms->Load(gHcParms->GetString("g_ctp_kinematics_filename"), RunNumber);
   gHcParms->Load(gHcParms->GetString("g_ctp_parm_filename"));
   // parameters not found in usual engine parameter files
   gHcParms->Load("PARAM/hcana.param");

src/THcRaster.cxx

@@ -96,6 +96,7 @@ Int_t THcRaster::ReadDatabase( const TDatime& date )
     {"fr_cal_mom",&fFrCalMom, kDouble},
     {"frx_adcpercm",&fFrXADCperCM, kDouble},
     {"fry_adcpercm",&fFrYADCperCM, kDouble},
+    {"pbeam",&fgpbeam, kDouble},
     {0}
   };
 
@@ -301,12 +302,13 @@ Int_t THcRaster::Process( ){
     gfry = (gfry_adc/gfry_adcpercm)*(gfr_cal_mom/ebeam)
   */
   // bpw- Iam not sure how to access the global variables. the below code doesnt work. so for now, hard code the beam energy.
-  if(gHcParms->Find("ebeam")){
-    tmp=*(Double_t *)gHcParms->Find("ebeam")->GetValuePointer();
-    cout <<"e Beam  = "<< tmp << endl;
+  eBeam=0.0001;
+  if(gHcParms->Find("gpbeam")){
+    eBeam=*(Double_t *)gHcParms->Find("gpbeam")->GetValuePointer();
   }
+  //  cout <<"e Beam  = "<< eBeam << " " << fgpbeam << endl;
   
-  eBeam = 4.02187;
+  //  eBeam = 4.02187;
   fXpos = (fXADC/fFrXADCperCM)*(fFrCalMom/eBeam);
   fYpos = (fYADC/fFrYADCperCM)*(fFrCalMom/eBeam);
 

src/THcRaster.h

@@ -40,6 +40,8 @@ class THcRaster : public THaBeamDet, public THcHitList {
   Int_t  ReadDatabase( const TDatime& date );
   Int_t DefineVariables( EMode mode );
 
+  Double_t       fgpbeam;     //
+
   Double_t       fXADC;     // X current
   Double_t       fYADC;     // Y current
   Double_t       fXpos;     // X position