SpectralFunction2p2h.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
 or see $GENIE/LICENSE

 Author: Afroditi Papadopoulou <apapadop \at mit.edu>
         Massachusetts Institute of Technology - October 04, 2019

 @ October 4, 2019 - Afroditi Papadopoulou (AP)
   Created this new module that controls the addition of the recoil nucleon in the event record
   and extracts its kinematics
*/
//____________________________________________________________________________

#include <cstdlib>

#include <TLorentzVector.h>
#include <TVector3.h>
#include <TParticlePDG.h>
#include <TMath.h>

#include "Framework/Algorithm/AlgFactory.h"
#include "Framework/Algorithm/AlgConfigPool.h"
#include "Framework/Conventions/Constants.h"
#include "Framework/Conventions/Units.h"
#include "Physics/NuclearState/SpectralFunction2p2h.h"

#include "Physics/NuclearState/NuclearModel.h"
#include "Physics/NuclearState/NuclearModelI.h"
#include "Framework/EventGen/EVGThreadException.h"
#include "Framework/GHEP/GHepRecord.h"
#include "Framework/GHEP/GHepParticle.h"
#include "Framework/GHEP/GHepStatus.h"
#include "Framework/GHEP/GHepFlags.h"
#include "Framework/Interaction/Interaction.h"
#include "Framework/Messenger/Messenger.h"
#include "Physics/NuclearState/FermiMomentumTablePool.h"
#include "Physics/NuclearState/FermiMomentumTable.h"
#include "Framework/Numerical/RandomGen.h"
#include "Framework/ParticleData/PDGLibrary.h"
#include "Framework/ParticleData/PDGUtils.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/Utils/KineUtils.h"
#include "Physics/NuclearState/NuclearUtils.h"

using namespace genie;
using namespace genie::constants;

//___________________________________________________________________________
SpectralFunction2p2h::SpectralFunction2p2h() :
SecondNucleonEmissionI("genie::SpectralFunction2p2h")
{

}
//___________________________________________________________________________
SpectralFunction2p2h::SpectralFunction2p2h(string config) :
  SecondNucleonEmissionI("genie::SpectralFunction2p2h", config)
{

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

}
//___________________________________________________________________________
void SpectralFunction2p2h::ProcessEventRecord(GHepRecord * evrec) const
{

    Interaction *  interaction = evrec       -> Summary();
    InitialState * init_state  = interaction -> InitStatePtr();
    Target *       tgt         = init_state  -> TgtPtr();

    if ( tgt -> A() <= 2 ) return ;
    if ( tgt -> Z() < 2 ) return ;

    FermiMoverInteractionType_t interaction_type = fNuclModel->GetFermiMoverInteractionType();
    
    if ( interaction_type == kFermiMoveEffectiveSF2p2h_eject ) {

        GHepParticle * nucleon = evrec->HitNucleon();
        int second_nucleon_pdg = nucleon->Pdg() == kPdgProton ? kPdgNeutron : kPdgProton ;
        SecondNucleonEmissionI::EmitSecondNucleon( evrec, second_nucleon_pdg );

    }

}
//____________________________________________________________________________
void SpectralFunction2p2h::Configure(const Registry & config)
{
  Algorithm::Configure(config);
  this->LoadConfig();
}
//____________________________________________________________________________
void SpectralFunction2p2h::Configure(string config)
{
  Algorithm::Configure(config);
  this->LoadConfig();
}
//____________________________________________________________________________
void SpectralFunction2p2h::LoadConfig(void)
{
  SecondNucleonEmissionI::LoadConfig() ;
}
//____________________________________________________________________________