genie::DMEOutgoingDarkGenerator Class Reference

Generates the final state primary lepton in neutrino-electron events. Is a concrete implementation of the EventRecordVisitorI interface. More...

#include <DMEOutgoingDarkGenerator.h>

Inheritance diagram for genie::DMEOutgoingDarkGenerator:

Collaboration diagram for genie::DMEOutgoingDarkGenerator:

Public Member Functions
	DMEOutgoingDarkGenerator ()
	DMEOutgoingDarkGenerator (string config)
	~DMEOutgoingDarkGenerator ()
void	ProcessEventRecord (GHepRecord *event_rec) const
Public Member Functions inherited from genie::OutgoingDarkGenerator
virtual void	SetPolarization (GHepRecord *ev) const
virtual TVector3	NucRestFrame2Lab (GHepRecord *ev) const
virtual void	AddToEventRecord (GHepRecord *ev, int pdgc, const TLorentzVector &p4) const
void	Configure (const Registry &config)
void	Configure (string config)
Public Member Functions inherited from genie::EventRecordVisitorI
virtual	~EventRecordVisitorI ()
Public Member Functions inherited from genie::Algorithm
virtual	~Algorithm ()
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...

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)
static string	BuildParamMatKey (const std::string &comm_name, unsigned int i, unsigned int j)
static string	BuildParamMatRowSizeKey (const std::string &comm_name)
static string	BuildParamMatColSizeKey (const std::string &comm_name)
Protected Member Functions inherited from genie::OutgoingDarkGenerator
	OutgoingDarkGenerator ()
	OutgoingDarkGenerator (string name)
	OutgoingDarkGenerator (string name, string config)
virtual	~OutgoingDarkGenerator ()
void	LoadConfig (void)
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
template<class T >
int	GetParamMat (const std::string &comm_name, TMatrixT< T > &mat, bool is_top_call=true) const
	Handle to load matrix of parameters. More...
template<class T >
int	GetParamMatSym (const std::string &comm_name, TMatrixTSym< T > &mat, bool is_top_call=true) const
int	GetParamMatKeys (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::OutgoingDarkGenerator
bool	fApplyCoulombCorrection
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

Generates the final state primary lepton in neutrino-electron events. Is a concrete implementation of the EventRecordVisitorI interface.

Author

Costas Andreopoulos <c.andreopoulos cern.ch> University of Liverpool

Created:

July 13, 2005

License:

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

Definition at line 27 of file DMEOutgoingDarkGenerator.h.

Constructor & Destructor Documentation

DMEOutgoingDarkGenerator::DMEOutgoingDarkGenerator

(

)

Definition at line 37 of file DMEOutgoingDarkGenerator.cxx.

                                                   :
OutgoingDarkGenerator("genie::DMEOutgoingDarkGenerator")
{

}

DMEOutgoingDarkGenerator::DMEOutgoingDarkGenerator

(

string

config

)

Definition at line 43 of file DMEOutgoingDarkGenerator.cxx.

                                                               :
OutgoingDarkGenerator("genie::DMEOutgoingDarkGenerator", config)
{

}

DMEOutgoingDarkGenerator::~DMEOutgoingDarkGenerator

(

)

Definition at line 49 of file DMEOutgoingDarkGenerator.cxx.

{

}

Member Function Documentation

void DMEOutgoingDarkGenerator::ProcessEventRecord ( GHepRecord *  event_rec ) const

virtual

Reimplemented from genie::OutgoingDarkGenerator.

Definition at line 54 of file DMEOutgoingDarkGenerator.cxx.

References genie::OutgoingDarkGenerator::AddToEventRecord(), genie::Interaction::FSPrimLepton(), genie::Interaction::FSPrimLeptonPdg(), genie::Interaction::InitState(), genie::RandomGen::Instance(), genie::constants::kElectronMass, genie::Interaction::Kine(), genie::constants::kPi, genie::kRfLab, LOG, genie::GHepParticle::P4(), pINFO, pNOTICE, genie::GHepRecord::Probe(), genie::InitialState::ProbeE(), pWARN, genie::utils::kinematics::Q2(), genie::RandomGen::RndLep(), genie::OutgoingDarkGenerator::SetPolarization(), genie::GHepRecord::Summary(), and genie::Kinematics::y().

{
// This method generates the final state primary lepton for DME events

  Interaction * interaction = evrec->Summary();
  const InitialState & init_state = interaction->InitState();

  // Get selected kinematics
  double y = interaction->Kine().y(true);
  assert(y>0 && y<1);

  // Final state primary lepton PDG code
  int pdgc = interaction->FSPrimLeptonPdg();
  assert(pdgc!=0);

  // Compute the neutrino and muon energy
  double Ev  = init_state.ProbeE(kRfLab); 
  double El  = y*Ev;

  LOG("LeptonicVertex", pINFO)
       << "Ev = " << Ev << ", y = " << y << ", -> El = " << El;

  // Compute the momentum transfer and scattering angle
  double Ev2   = TMath::Power(Ev,2);
  double El2   = TMath::Power(El,2);
  double me    = kElectronMass;
  double ml    = interaction->FSPrimLepton()->Mass();
  double ml2   = TMath::Power(ml,2);
  double pl    = TMath::Sqrt(El2-ml2);   
  double pv    = TMath::Sqrt(Ev2-ml2);
  
  assert(El2>=ml2&&Ev2>=ml2);

  double Q2    = 2*(Ev-El)*me;
  double costh = (El*Ev - ml2 -0.5*Q2)/pl/pv;
  double sinth = TMath::Sqrt( TMath::Max(0., 1-TMath::Power(costh,2.)) );

  LOG("LeptonicVertex", pNOTICE)
       << "Q2 = " << Q2 << ", cos(theta) = " << costh;

  //warn about overflow in costheta and ignore it if it is small or abort
  if( TMath::Abs(costh)>1 ) {
     LOG("LeptonicVertex", pWARN)
       << "El = " << El << ", Ev = " << Ev << ", cos(theta) = " << costh;
     //if(TMath::Abs(costh)-1<0.3) costh = 1.0; //why?
  }
  assert(TMath::Abs(costh)<=1);

  // Compute the p components along and perpendicular the v direction 
  double plp = pl * costh; // p(//)
  double plt = pl * sinth; // p(-|)

  LOG("LeptonicVertex", pNOTICE)
        << "fsl: E = " << El << ", |p//| = " << plp << "[pT] = " << plt;

  // Randomize transverse components
  RandomGen * rnd = RandomGen::Instance();
  double phi  = 2*kPi * rnd->RndLep().Rndm();
  double pltx = plt * TMath::Cos(phi);
  double plty = plt * TMath::Sin(phi);

  // Take a unit vector along the neutrino direction
  TVector3 unit_nudir = evrec->Probe()->P4()->Vect().Unit();

  // Rotate lepton momentum vector from the reference frame (x'y'z') where 
  // {z':(neutrino direction), z'x':(theta plane)} to the LAB
  TVector3 p3l(pltx,plty,plp);
  p3l.RotateUz(unit_nudir);

  // Lepton 4-momentum in the LAB
  TLorentzVector p4l(p3l,El);

  // Create a GHepParticle and add it to the event record
  this->AddToEventRecord(evrec, pdgc, p4l);

  // Set final state lepton polarization
  this->SetPolarization(evrec);
}

---

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

• DMEOutgoingDarkGenerator.h
• DMEOutgoingDarkGenerator.cxx