#include <AMNuGammaGenerator.h>

Inheritance diagram for genie::AMNuGammaGenerator:

Public Member Functions
	AMNuGammaGenerator ()

	AMNuGammaGenerator (string config)

	~AMNuGammaGenerator ()

void	ProcessEventRecord (GHepRecord *event_rec) const

Public Member Functions inherited from genie::EventRecordVisitorI
virtual	~EventRecordVisitorI ()

Public Member Functions inherited from genie::Algorithm
virtual	~Algorithm ()

virtual void	Configure (const Registry &config)

virtual void	Configure (string config)

virtual void	FindConfig (void)

virtual const Registry &	GetConfig (void) const

Registry *	GetOwnedConfig (void)

virtual const AlgId &	Id (void) const
	Get algorithm ID. More...

virtual AlgStatus_t	GetStatus (void) const
	Get algorithm status. More...

virtual bool	AllowReconfig (void) const

virtual AlgCmp_t	Compare (const Algorithm *alg) const
	Compare with input algorithm. More...

virtual void	SetId (const AlgId &id)
	Set algorithm ID. More...

virtual void	SetId (string name, string config)

const Algorithm *	SubAlg (const RgKey &registry_key) const

void	AdoptConfig (void)

void	AdoptSubstructure (void)

virtual void	Print (ostream &stream) const
	Print algorithm info. More...

Private Member Functions
void	AddPhoton (GHepRecord *event_rec) const

void	AddFinalStateNeutrino (GHepRecord *event_rec) const

void	AddTargetRemnant (GHepRecord *event_rec) const

void	AddRecoilNucleon (GHepRecord *event_rec) const

Additional Inherited Members
Static Public Member Functions inherited from genie::Algorithm
static string	BuildParamVectKey (const std::string &comm_name, unsigned int i)

static string	BuildParamVectSizeKey (const std::string &comm_name)

Protected Member Functions inherited from genie::EventRecordVisitorI
	EventRecordVisitorI ()

	EventRecordVisitorI (string name)

	EventRecordVisitorI (string name, string config)

Protected Member Functions inherited from genie::Algorithm
	Algorithm ()

	Algorithm (string name)

	Algorithm (string name, string config)

void	Initialize (void)

void	DeleteConfig (void)

void	DeleteSubstructure (void)

Registry *	ExtractLocalConfig (const Registry &in) const

Registry *	ExtractLowerConfig (const Registry &in, const string &alg_key) const
	Split an incoming configuration Registry into a block valid for the sub-algo identified by alg_key. More...

template<class T >
bool	GetParam (const RgKey &name, T &p, bool is_top_call=true) const

template<class T >
bool	GetParamDef (const RgKey &name, T &p, const T &def) const

template<class T >
int	GetParamVect (const std::string &comm_name, std::vector< T > &v, bool is_top_call=true) const
	Handle to load vectors of parameters. More...

int	GetParamVectKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const

int	AddTopRegistry (Registry *rp, bool owns=true)
	add registry with top priority, also update ownership More...

int	AddLowRegistry (Registry *rp, bool owns=true)
	add registry with lowest priority, also update ownership More...

int	MergeTopRegistry (const Registry &r)

int	AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
	Add registries with top priority, also udated Ownerships. More...

Protected Attributes inherited from genie::Algorithm
bool	fAllowReconfig

bool	fOwnsSubstruc
	true if it owns its substructure (sub-algs,...) More...

AlgId	fID
	algorithm name and configuration set More...

vector< Registry * >	fConfVect

vector< bool >	fOwnerships
	ownership for every registry in fConfVect More...

AlgStatus_t	fStatus
	algorithm execution status More...

AlgMap *	fOwnedSubAlgMp
	local pool for owned sub-algs (taken out of the factory pool) More...

Detailed Description

Author: Costas Andreopoulos <constantinos.andreopoulos cern.ch> University of Liverpool & STFC Rutherford Appleton Laboratory

Feb 15, 2008

Definition at line 25 of file AMNuGammaGenerator.h.

Constructor & Destructor Documentation

AMNuGammaGenerator::AMNuGammaGenerator ( )

Definition at line 28 of file AMNuGammaGenerator.cxx.

                                        :
 EventRecordVisitorI("genie::AMNuGammaGenerator")
 {
 
 }

AMNuGammaGenerator::AMNuGammaGenerator ( string config )

Definition at line 34 of file AMNuGammaGenerator.cxx.

                                                     :
 EventRecordVisitorI("genie::AMNuGammaGenerator", config)
 {
 
 }

AMNuGammaGenerator::~AMNuGammaGenerator ( )

Definition at line 40 of file AMNuGammaGenerator.cxx.

 {
 
 }

Member Function Documentation

void AMNuGammaGenerator::AddFinalStateNeutrino ( GHepRecord * event_rec ) const

private

Definition at line 112 of file AMNuGammaGenerator.cxx.

 {
 // Adding the final state neutrino
 // Just use 4-momentum conservation (init_neutrino = photon + final_neutrino)
 
   LOG("AMNuGammaGenerator", pINFO) << "Adding final state neutrino";
 
   GHepParticle * nu    = evrec->Probe();                       // incoming v
   GHepParticle * gamma = evrec->Particle(nu->FirstDaughter()); // gamma
   assert(nu);
   assert(gamma);
 
   const TLorentzVector & vtx = *(nu->X4()); // vtx
 
   const TLorentzVector & p4nu_lab    = *(nu->P4());
   const TLorentzVector & p4gamma_lab = *(gamma->P4());
   TLorentzVector p4_lab = p4nu_lab - p4gamma_lab;
 
   GHepParticle p(nu->Pdg(), kIStStableFinalState, 0,-1,-1,-1, p4_lab, vtx);
   evrec->AddParticle(p);
 }

void AMNuGammaGenerator::AddPhoton ( GHepRecord * event_rec ) const

private

Definition at line 53 of file AMNuGammaGenerator.cxx.

 {
 // Adding the final state photon
 //
   LOG("AMNuGammaGenerator", pINFO) << "Adding final state photon";
 
   RandomGen * rnd = RandomGen::Instance();
 
   // Get boost vector for transforms between LAB <-> NRF (nucleon rest frame)
   GHepParticle * nuc = evrec->HitNucleon();
   const TLorentzVector & p4nuc_lab = *(nuc->P4());
   TVector3 beta = p4nuc_lab.BoostVector();
 
   // Get the neutrino 4-momentum at the LAB
   GHepParticle * nu = evrec->Probe();
   const TLorentzVector & p4v_lab = *(nu->P4());
 
   // Get the neutrino 4-momentum at the NRF
   TLorentzVector p4v_nrf = p4v_lab;
   p4v_nrf.Boost(-1.*beta);
 
   // Generate the photon cos(theta) with respect to the neutrino direction
   // (at NRF) and a uniform azimuthal angle phi
   double costheta_gamma = -1.0 + 2.0 * rnd->RndKine().Rndm();
   double phi_gamma      =  2.0 * kPi * rnd->RndKine().Rndm();
 
   // Generate the fraction of the neutrino energy taken by the photon
   double efrac_gamma =  rnd->RndKine().Rndm();
 
   // Calculate the photon energy at the NRF
   double Ev_nrf = p4v_nrf.Energy();
   double Egamma_nrf = Ev_nrf * efrac_gamma;
 
   // Calculate the photon momentum components at a rotated NRF (NRF')
   // where z is along the neutrino direction
   double sintheta_gamma = TMath::Sqrt(1-TMath::Power(costheta_gamma,2));
   double pgamma_nrf_p = Egamma_nrf * costheta_gamma; // p(//)
   double pgamma_nrf_t = Egamma_nrf * sintheta_gamma; // p(-|)
   double px = pgamma_nrf_t * TMath::Sin(phi_gamma);
   double py = pgamma_nrf_t * TMath::Cos(phi_gamma);
   double pz = pgamma_nrf_p;
 
   // Take a unit vector along the neutrino direction @ the NRF
   TVector3 unit_nudir = p4v_nrf.Vect().Unit();
 
   // Rotate the photon momentum vector from NRF' -> NRF
   TVector3 p3gamma_nrf(px,py,pz);
   p3gamma_nrf.RotateUz(unit_nudir);
 
   // Get the photon 4-momentum back at the LAB
   TLorentzVector p4gamma_lab(p3gamma_nrf, Egamma_nrf);
   p4gamma_lab.Boost(beta);
 
   // Add the photon at the event record
   const TLorentzVector & vtx = *(nu->X4());
   GHepParticle p(kPdgGamma,kIStStableFinalState,0,-1,-1,-1,p4gamma_lab,vtx);
   evrec->AddParticle(p);
 }

void AMNuGammaGenerator::AddRecoilNucleon ( GHepRecord * event_rec ) const

private

p.SetRemovalEnergy(w);

Definition at line 134 of file AMNuGammaGenerator.cxx.

 {
 // Adding the recoil nucleon.
 // Recoil is negligible at the H3 model. However, for nuclear targets, the
 // hit nucleon was off-the-mass-shell (bound). Doing some magic here to bring
 // it back on-the-mass-shell so that it can be INTRANUKE'ed and appear in
 // the final state. The nucleon will keep its original (Fermi) 3-momentum but
 // take some energy back from the remnant nucleus. No such tweaking takes
 // place for free nucleon targets.
 
   LOG("AMNuGammaGenerator", pINFO) << "Adding recoil nucleon";
 
   // Get the interaction target
   GHepParticle * tgt_nucleus = evrec->TargetNucleus();
   bool is_nuclear_target = (tgt_nucleus!=0);
 
   // Get the hit nucleon
   GHepParticle * hitnuc = evrec->HitNucleon();
   assert(hitnuc);
 
   // Get the hit nucleon pdg code (= recoil nucleon pdg code)
   int pdgc = hitnuc->Pdg();
 
   // Get the hit nucleon 4-momentum (LAB)
   const TLorentzVector & p4n = *(hitnuc->P4());
   TLorentzVector p4(p4n);
 
   // Tweak the 4-momentum to bring the recoil nucleon on-the-mass-shell
   // (for nuclear targets only)
   if (is_nuclear_target) {
     double p = p4.Vect().Mag();
     double m = PDGLibrary::Instance()->Find(pdgc)->Mass();
     double E = TMath::Sqrt(m*m+p*p);
     p4.SetE(E);
   }
 
   // Get the vtx position
   GHepParticle * neutrino  = evrec->Probe();
   const TLorentzVector & vtx = *(neutrino->X4());
 
   // Add the recoil nucleon at the event record
   GHepStatus_t ist = (is_nuclear_target) ?
                       kIStHadronInTheNucleus : kIStStableFinalState;
   int mom = evrec->HitNucleonPosition();
 
   LOG("AMNuGammaGenerator", pINFO)
                   << "Adding recoil baryon [pdgc = " << pdgc << "]";
 
   GHepParticle p(pdgc, ist, mom,-1,-1,-1, p4, vtx);
 //  double w = hitnuc->RemovalEnergy();
 ///  p.SetRemovalEnergy(w);
   evrec->AddParticle(p);
 }

void AMNuGammaGenerator::AddTargetRemnant ( GHepRecord * event_rec ) const

private

Definition at line 188 of file AMNuGammaGenerator.cxx.

 {
 // Add the remnant nuclear target at the GHEP record
 
   LOG("AMNuGammaGenerator", pINFO) << "Adding final state nucleus";
 
   // Skip for non nuclear targets
   GHepParticle * nucleus = evrec->TargetNucleus();
   if (!nucleus) {
     LOG("AMNuGammaGenerator", pDEBUG)
      << "No nucleus in the initial state - no remnant nucleus to add in the f/s";
     return;
   }
 
   // Compute A,Z for final state nucleus & get its PDG code and its mass
   //GHepParticle * nucleon = evrec->HitNucleon();
 
   GHepParticle * hit_nucleon = evrec->HitNucleon(); // hit
   GHepParticle * rec_nucleon = evrec->Particle(hit_nucleon->FirstDaughter()); // recoil
 
   assert(rec_nucleon);
   int  npdgc = rec_nucleon->Pdg();
   bool is_p  = pdg::IsProton(npdgc);
   int A = nucleus->A();
   int Z = nucleus->Z();
   if (is_p) Z--;
   A--;
   int ipdgc = pdg::IonPdgCode(A, Z);
   TParticlePDG * remnant = PDGLibrary::Instance()->Find(ipdgc);
   if(!remnant) {
       LOG("AMNuGammaGenerator", pFATAL)
           << "No particle with [A = " << A << ", Z = " << Z
                             << ", pdgc = " << ipdgc << "] in PDGLibrary!";
       assert(remnant);
   }
 
   // Figure out the remnant 4-momentum
   double Mf  = remnant->Mass();
   double Mf2 = TMath::Power(Mf,2);
   double px  = -1.* rec_nucleon->Px();
   double py  = -1.* rec_nucleon->Py();
   double pz  = -1.* rec_nucleon->Pz();
   double E   = TMath::Sqrt(Mf2 + rec_nucleon->P4()->Vect().Mag2());
   E += (hit_nucleon->P4()->E() - rec_nucleon->P4()->E());
 
   //-- Add the nucleus to the event record
   LOG("AMNuGammaGenerator", pINFO)
        << "Adding nucleus [A = " << A << ", Z = " << Z
                                            << ", pdgc = " << ipdgc << "]";
   int mom = evrec->TargetNucleusPosition();
   evrec->AddParticle(
            ipdgc,kIStStableFinalState, mom,-1,-1,-1, px,py,pz,E, 0,0,0,0);
 }

void AMNuGammaGenerator::ProcessEventRecord ( GHepRecord * event_rec ) const

virtual

Implements genie::EventRecordVisitorI.

Definition at line 45 of file AMNuGammaGenerator.cxx.

 {
   this->AddPhoton(evrec);
   this->AddFinalStateNeutrino(evrec);
   this->AddRecoilNucleon(evrec);
 //this->AddTargetRemnant(evrec);
 }

The documentation for this class was generated from the following files:

Generator/src/Physics/AnomalyMediatedNuGamma/EventGen/AMNuGammaGenerator.h
Generator/src/Physics/AnomalyMediatedNuGamma/EventGen/AMNuGammaGenerator.cxx

Public Member Functions

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation