gEvGenDM.cxx File Reference

#include <cstdlib>
#include <cassert>
#include <sstream>
#include <string>
#include <vector>
#include <map>
#include <TFile.h>
#include <TTree.h>
#include <TSystem.h>
#include <TVector3.h>
#include <TH1.h>
#include <TF1.h>
#include "Framework/Algorithm/AlgConfigPool.h"
#include "Framework/Conventions/XmlParserStatus.h"
#include "Framework/Conventions/GBuild.h"
#include "Framework/Conventions/Controls.h"
#include "Framework/Conventions/Constants.h"
#include "Framework/Conventions/Units.h"
#include "Framework/EventGen/EventRecord.h"
#include "Framework/EventGen/GFluxI.h"
#include "Framework/EventGen/GEVGDriver.h"
#include "Framework/EventGen/GMCJDriver.h"
#include "Framework/EventGen/GMCJMonitor.h"
#include "Framework/Interaction/Interaction.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/Ntuple/NtpWriter.h"
#include "Framework/Ntuple/NtpMCFormat.h"
#include "Framework/Numerical/RandomGen.h"
#include "Framework/Numerical/Spline.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGLibrary.h"
#include "Framework/ParticleData/PDGUtils.h"
#include "Framework/Utils/AppInit.h"
#include "Framework/Utils/RunOpt.h"
#include "Framework/Utils/XSecSplineList.h"
#include "Framework/Utils/StringUtils.h"
#include "Framework/Utils/PrintUtils.h"
#include "Framework/Utils/SystemUtils.h"
#include "Framework/Utils/CmdLnArgParser.h"

Include dependency graph for gEvGenDM.cxx:

Go to the source code of this file.

Functions
void	GetCommandLineArgs (int argc, char **argv)
void	Initialize (void)
void	PrintSyntax (void)
bool	CheckUnitarityLimit (void)
void	GenerateEventsAtFixedInitState (void)
int	main (int argc, char **argv)

Variables
int	kDefOptNevents = 0
NtpMCFormat_t	kDefOptNtpFormat = kNFGHEP
Long_t	kDefOptRunNu = 0
int	gOptNevents
double	gOptDMEnergy
double	gOptDMEnergyRange
double	gOptDMMass
double	gOptZpCoupling
map< int, double >	gOptTgtMix
double	gOptMedRatio
Long_t	gOptRunNu
string	gOptFlux
bool	gOptWeighted
bool	gOptUsingFluxOrTgtMix = false
long int	gOptRanSeed
string	gOptInpXSecFile
string	gOptOutFileName
string	gOptStatFileName

Function Documentation

bool CheckUnitarityLimit

(

void

)

Definition at line 876 of file gEvGenDM.cxx.

References genie::units::cm2, genie::AlgConfigPool::CommonList(), genie::Registry::Get(), gOptDMEnergy, gOptDMMass, gOptMedRatio, genie::AlgConfigPool::Instance(), genie::constants::kNucleonMass, genie::constants::kPi, LOG, and pWARN.

{
  // Before generating the events, perform a simple sanity check
  // We estimate the leading divergent piece of the cross-section
  // We make sure it does not exceed the unitarity limit
  double gzp;
  Registry * r = AlgConfigPool::Instance()->CommonList("Param", "BoostedDarkMatter");
  r->Get("ZpCoupling", gzp);
  double gzp4 = TMath::Power(gzp,4);
  double Mzp  = gOptMedRatio * gOptDMMass;
  double Mzp2 = Mzp*Mzp;
  // The leading, forward-dominated piece is the same for both DM models
  double xsec_est = gzp4 / (4. * kPi * Mzp2);
  double ml   = gOptDMMass;
  double ml2  = ml*ml;
  double M    = kNucleonMass;
  double M2   = M*M;
  double Ed   = gOptDMEnergy;
  double Ed2  = Ed*Ed;
  double pcm2 = M2 * (Ed2 - ml2) / (ml2 + M2 + 2.*M*Ed);
  double xsec_lim = kPi / pcm2;
  bool unitary = xsec_lim > xsec_est;
  if (!unitary) {
    LOG("gevgen_dm", pWARN)
      << "Estimated a cross-section " << xsec_est/cm2 << " cm^2";
    LOG("gevgen_dm", pWARN)
      << "Unitarity limit set to " << xsec_lim/cm2 << " cm^2";
  }
  return unitary;
}

void GenerateEventsAtFixedInitState

(

void

)

void GetCommandLineArgs	(	int	argc,
		char **	argv
	)

void Initialize

(

void

)

int main	(	int	argc,
		char **	argv
	)

Definition at line 235 of file gEvGenDM.cxx.

References genie::PDGLibrary::AddDarkMatter(), genie::AlgConfigPool::CommonList(), GenerateEventsAtFixedInitState(), GetCommandLineArgs(), gOptDMMass, gOptMedRatio, gOptUsingFluxOrTgtMix, gOptZpCoupling, Initialize(), genie::PDGLibrary::Instance(), genie::AlgConfigPool::Instance(), genie::Registry::Lock(), LOG, pERROR, genie::Registry::Set(), and genie::Registry::UnLock().

{
  GetCommandLineArgs(argc,argv);
  PDGLibrary::Instance()->AddDarkMatter(gOptDMMass,gOptMedRatio);
  if (gOptZpCoupling > 0.) {
      Registry * r = AlgConfigPool::Instance()->CommonList("Param", "BoostedDarkMatter");
      r->UnLock();
      r->Set("ZpCoupling", gOptZpCoupling);
      r->Lock();
  }
  Initialize();


  // throw on NaNs and Infs...
#if defined(HAVE_FENV_H) && defined(HAVE_FEENABLEEXCEPT)
  feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
#endif
  //
  // Generate dark matter events
  //

  if(gOptUsingFluxOrTgtMix) {
#ifdef __CAN_GENERATE_EVENTS_USING_A_FLUX_OR_TGTMIX__
        GenerateEventsUsingFluxOrTgtMix();
#else
  LOG("gevgen_dm", pERROR)
    << "\n   To be able to generate dark matter events from a flux and/or a target mix"
    << "\n   you need to add the following config options at your GENIE installation:"
    << "\n   --enable-flux-drivers  --enable-geom-drivers \n" ;
#endif
  } else {
     GenerateEventsAtFixedInitState();
  }
  return 0;
}

void PrintSyntax

(

void

)

Variable Documentation

double gOptDMEnergy

Definition at line 219 of file gEvGenDM.cxx.

Referenced by CheckUnitarityLimit().

double gOptDMEnergyRange

Definition at line 220 of file gEvGenDM.cxx.

double gOptDMMass

Definition at line 221 of file gEvGenDM.cxx.

Referenced by CheckUnitarityLimit(), and main().

string gOptFlux

Definition at line 226 of file gEvGenDM.cxx.

string gOptInpXSecFile

Definition at line 230 of file gEvGenDM.cxx.

double gOptMedRatio

Definition at line 224 of file gEvGenDM.cxx.

Referenced by CheckUnitarityLimit(), and main().

int gOptNevents

Definition at line 218 of file gEvGenDM.cxx.

string gOptOutFileName

Definition at line 231 of file gEvGenDM.cxx.

long int gOptRanSeed

Definition at line 229 of file gEvGenDM.cxx.

Long_t gOptRunNu

Definition at line 225 of file gEvGenDM.cxx.

string gOptStatFileName

Definition at line 232 of file gEvGenDM.cxx.

map<int,double> gOptTgtMix

Definition at line 223 of file gEvGenDM.cxx.

bool gOptUsingFluxOrTgtMix = false

Definition at line 228 of file gEvGenDM.cxx.

bool gOptWeighted

Definition at line 227 of file gEvGenDM.cxx.

double gOptZpCoupling

Definition at line 222 of file gEvGenDM.cxx.

Referenced by main().

int kDefOptNevents = 0

Definition at line 213 of file gEvGenDM.cxx.

NtpMCFormat_t kDefOptNtpFormat = kNFGHEP

Definition at line 214 of file gEvGenDM.cxx.

Long_t kDefOptRunNu = 0

Definition at line 215 of file gEvGenDM.cxx.