gtestInteraction.cxx File Reference

#include "Framework/Conventions/Constants.h"
#include "Framework/Interaction/Interaction.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/ParticleData/PDGCodes.h"

Include dependency graph for gtestInteraction.cxx:

Go to the source code of this file.

Functions
int	main (int, char **)

Function Documentation

int main	(	int	,
		char **
	)

Definition at line 28 of file gtestInteraction.cxx.

References genie::Kinematics::ClearRunningValues(), genie::Interaction::InitState(), genie::Interaction::InitStatePtr(), genie::Interaction::KinePtr(), genie::kIntWeakCC, genie::constants::kNucleonMass, genie::kPdgNuMu, genie::kPdgProton, genie::kPdgUQuark, genie::kScDeepInelastic, LOG, pINFO, genie::Target::SetHitNucP4(), genie::Target::SetHitNucPdg(), genie::Target::SetHitQrkPdg(), genie::InitialState::SetProbeP4(), genie::Kinematics::SetQ2(), genie::Kinematics::SetW(), genie::Kinematics::Setx(), genie::Kinematics::Sety(), genie::InitialState::Tgt(), genie::InitialState::TgtPtr(), and genie::Kinematics::UseSelectedKinematics().

{
  //-- build an initial state
  Target         Fe(26,56);
  int            numu(kPdgNuMu); 
  TLorentzVector pnu(0,0,8,8);
  InitialState   init_state(Fe, numu);

  init_state.SetProbeP4(pnu);

  //-- build process info
  ProcessInfo proc(kScDeepInelastic, kIntWeakCC);

  //-- create an interaction & print it
  Interaction interaction(init_state, proc);

  LOG("test", pINFO) << "Printing an interaction object";
  LOG("test", pINFO) << interaction;

  //-- set struck nucleon & quark info in the initial state's Target
  //   note: using methods ending in Ptr -> they return a 'writable'
  //         object that can be modified

  TLorentzVector pnucl(0,0,0,kNucleonMass);

  interaction.InitStatePtr()->TgtPtr()->SetHitNucP4(pnucl);
  interaction.InitStatePtr()->TgtPtr()->SetHitNucPdg(kPdgProton);
  interaction.InitStatePtr()->TgtPtr()->SetHitQrkPdg(kPdgUQuark);

  //-- get a 'read-only' InitialState and print it (check that struck nucleon
  //   and quark were set) 
  //   note: using the methods not ending in Ptr to get a 'read-only' object

  const InitialState & cinit = interaction.InitState();

  LOG("test", pINFO) << "Printing initial state after setting struck nucl/qrk";
  LOG("test", pINFO) << "\n" << cinit;

  //-- take just the Target from the initial state and print it

  LOG("test", pINFO) << "Printing target after setting struck nucl/qrk";
  const Target & ctgt = interaction.InitState().Tgt();
  LOG("test", pINFO) << "\n" <<ctgt;

  //-- change the struck nucleon
  //-- instead of using the long syntax above, get a writable Target object first

  Target * wtgt = interaction.InitStatePtr()->TgtPtr();
  wtgt->SetHitNucPdg(kPdgProton);

  LOG("test", pINFO) << "Printing target after changing struck nucl";
  LOG("test", pINFO) << "\n" << *wtgt;

  //-- take the Kinematics object and set some

  Kinematics * wkine = interaction.KinePtr();
 
  wkine->Setx(0.1781);
  wkine->Sety(0.6892);
  wkine->SetQ2(3.2218);

  LOG("test", pINFO) << "Printing kinematics after setting x,y,Q2";
  LOG("test", pINFO) << "\n" <<*wkine;

  //-- modify some & add a new
  wkine->Setx(0.2);
  wkine->Sety(0.2);
  wkine->SetW(1.9219);

  LOG("test", pINFO) << "Printing kinematics after modifying x,y and adding W";
  LOG("test", pINFO) << "\n" << *wkine;

  //-- now set "selected" kinematics
  wkine->Setx(0.25, true);
  wkine->Sety(0.21, true);
  wkine->SetW(2.89, true);

  LOG("test", pINFO) << "Printing kinematics after setting 'selected'";
  LOG("test", pINFO) << "\n" << *wkine;

  //-- now delete 'running' kinematics
  wkine->ClearRunningValues();

  LOG("test", pINFO) << "Printing kinematics after deleting 'running'";
  LOG("test", pINFO) << "\n" << *wkine;

  //-- copy the 'selected' kinematics to the 'running' ones
  wkine->UseSelectedKinematics();

  LOG("test", pINFO) << "Printing kinematics after copying 'selected'";
  LOG("test", pINFO) << "\n" << *wkine;

  //-- see that the interaction was updated

  LOG("test", pINFO) << "Printing the interaction object after all changes";
  LOG("test", pINFO) << interaction;
}