NormGenerator.cxx

Go to the documentation of this file.

//____________________________________________________________________________
/*
 Copyright (c) 2003-2024, The GENIE Collaboration
 For the full text of the license visit http://copyright.genie-mc.org

 Igor Kakorin <kakorin@jinr.ru>
 Joint Institute for Nuclear Research
*/
//____________________________________________________________________________

#include <cstdlib>

#include <TLorentzVector.h>


#include "Framework/Algorithm/AlgConfigPool.h"
#include "Framework/Conventions/GBuild.h"
#include "Framework/EventGen/EventGeneratorI.h"
#include "Framework/EventGen//RunningThreadInfo.h"
#include "Framework/GHEP/GHepRecord.h"
#include "Framework/GHEP/GHepParticle.h"
#include "Physics/Common/NormGenerator.h"

using namespace genie;

//___________________________________________________________________________
NormGenerator::NormGenerator():
EventRecordVisitorI("genie::NormGenerator")
{

}
//___________________________________________________________________________
NormGenerator::NormGenerator(string config):
EventRecordVisitorI("genie::NormGenerator")
{

}
//___________________________________________________________________________
NormGenerator::~NormGenerator()
{

}
//___________________________________________________________________________
void NormGenerator::ProcessEventRecord(GHepRecord * evrec) const
{
  Interaction * interaction = evrec->Summary();
  const InitialState & init_state = interaction -> InitState();

  // Access cross section algorithm for running thread
  RunningThreadInfo * rtinfo = RunningThreadInfo::Instance();
  const EventGeneratorI * evg = rtinfo->RunningThread();
  fXSecModel = evg->CrossSectionAlg();

  //bool isHeavyNucleus = tgt->A()>=3;
  GHepParticle * probe  = evrec->Probe();
  int pdg_probe = probe->Pdg();
  int iprobe = evrec->ProbePosition();
  TLorentzVector * p4v = probe->GetP4();
  TLorentzVector * vtx = probe->X4();
  
  GHepParticle * nucltgt = evrec->TargetNucleus();
  int pdg_tgt = nucltgt->Pdg();
  int inucltgt = evrec->TargetNucleusPosition();
  TLorentzVector * p4_tgt = nucltgt->GetP4();
  TLorentzVector * vtx_tgt = nucltgt->X4();

  
  evrec->AddParticle(pdg_probe, kIStStableFinalState, iprobe,-1,-1,-1, *p4v, *vtx);
  evrec->AddParticle(pdg_tgt, kIStStableFinalState, inucltgt,-1,-1,-1, *p4_tgt, *vtx_tgt);
  
   //-- Determine the status code
  //const Target & tgt = interaction->InitState().Tgt();


  // update the interaction summary
  evrec->Summary()->KinePtr()->SetFSLeptonP4(*p4v);
  
  double xsec = fXSecModel->XSec(interaction, kPSfE);
  
  evrec->SetDiffXSec(xsec, kPSfE);
  
  return;
}
//___________________________________________________________________________
void NormGenerator::Configure(const Registry & config)
{
  Algorithm::Configure(config);
  this->LoadConfig();
}
//____________________________________________________________________________
void NormGenerator::Configure(string config)
{
  Algorithm::Configure(config);
  this->LoadConfig();
}
//____________________________________________________________________________
void NormGenerator::LoadConfig(void)
{
}
//____________________________________________________________________________