A module to check the results from the Monte Carlo generator. More...

Inheritance diagram for evgen::NDKGen:

Public Types
using	ModuleType = EDProducer

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Member Functions
	NDKGen (fhicl::ParameterSet const &pset)

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

void	fillProductDescriptions ()

void	registerProducts (ProductDescriptions &productsToRegister)

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Private Member Functions
void	produce (art::Event &evt) override

void	beginJob () override

void	beginRun (art::Run &run) override

void	endJob () override

std::string	ParticleStatus (int StatusCode)

std::string	ReactionChannel (int ccnc, int mode)

void	FillHistograms (simb::MCTruth const &mc)

Private Attributes
std::string	fNdkFile

std::ifstream	fEventFile

std::string	fNDKModuleLabel

CLHEP::HepRandomEngine &	fEngine
	art-managed random-number engine More...

TH1F *	fGenerated [6]
	Spectra as generated. More...

TH1F *	fVertexX
	vertex location of generated events in x More...

TH1F *	fVertexY
	vertex location of generated events in y More...

TH1F *	fVertexZ
	vertex location of generated events in z More...

TH2F *	fVertexXY
	vertex location in xy More...

TH2F *	fVertexXZ
	vertex location in xz More...

TH2F *	fVertexYZ
	vertex location in yz More...

TH1F *	fDCosX
	direction cosine in x More...

TH1F *	fDCosY
	direction cosine in y More...

TH1F *	fDCosZ
	direction cosine in z More...

TH1F *	fMuMomentum
	momentum of outgoing muons More...

TH1F *	fMuDCosX
	direction cosine of outgoing mu in x More...

TH1F *	fMuDCosY
	direction cosine of outgoing mu in y More...

TH1F *	fMuDCosZ
	direction cosine of outgoing mu in z More...

TH1F *	fEMomentum
	momentum of outgoing electrons More...

TH1F *	fEDCosX
	direction cosine of outgoing e in x More...

TH1F *	fEDCosY
	direction cosine of outgoing e in y More...

TH1F *	fEDCosZ
	direction cosine of outgoing e in z More...

TH1F *	fCCMode
	CC interaction mode. More...

TH1F *	fNCMode
	CC interaction mode. More...

TH1F *	fECons
	histogram to determine if energy is conserved in the event More...

Detailed Description

A module to check the results from the Monte Carlo generator.

Definition at line 48 of file NDKGen_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 17 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 26 of file Producer.h.

Constructor & Destructor Documentation

evgen::NDKGen::NDKGen ( fhicl::ParameterSet const & pset )

explicit

Definition at line 103 of file NDKGen_module.cc.

References art::detail::EngineCreator::createEngine(), fEngine, fEventFile, and fNdkFile.

     : EDProducer{pset}
     , fNdkFile{pset.get<std::string>("NdkFile")}
     , fEventFile{fNdkFile}
     // create a default random engine; obtain the random seed from NuRandomService,
     // unless overridden in configuration with key "Seed"
     , fEngine(art::ServiceHandle<rndm::NuRandomService>()->registerAndSeedEngine(createEngine(0),
                                                                                  pset,
                                                                                  "Seed"))
   {
     produces<std::vector<simb::MCTruth>>();
     produces<sumdata::RunData, art::InRun>();
 
     if (!fEventFile.good()) {
       throw cet::exception("NDKGen") << "Could not open file: " << fNdkFile << '\n';
     }
   }

Member Function Documentation

void evgen::NDKGen::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDProducer.

Definition at line 123 of file NDKGen_module.cc.

References geo::GeometryCore::DetHalfHeight(), geo::GeometryCore::DetHalfWidth(), geo::GeometryCore::DetLength(), fCCMode, fDCosX, fDCosY, fDCosZ, fECons, fEDCosX, fEDCosY, fEDCosZ, fEMomentum, fGenerated, fMuDCosX, fMuDCosY, fMuDCosZ, fMuMomentum, fNCMode, fVertexX, fVertexXY, fVertexXZ, fVertexY, fVertexYZ, fVertexZ, x, y, and z.

   {
     art::ServiceHandle<art::TFileService const> tfs;
 
     fGenerated[0] = tfs->make<TH1F>("fGenerated_necc", "", 100, 0.0, 20.0);
     fGenerated[1] = tfs->make<TH1F>("fGenerated_nebcc", "", 100, 0.0, 20.0);
     fGenerated[2] = tfs->make<TH1F>("fGenerated_nmcc", "", 100, 0.0, 20.0);
     fGenerated[3] = tfs->make<TH1F>("fGenerated_nmbcc", "", 100, 0.0, 20.0);
     fGenerated[4] = tfs->make<TH1F>("fGenerated_nnc", "", 100, 0.0, 20.0);
     fGenerated[5] = tfs->make<TH1F>("fGenerated_nbnc", "", 100, 0.0, 20.0);
 
     fDCosX = tfs->make<TH1F>("fDCosX", ";dx/ds", 200, -1., 1.);
     fDCosY = tfs->make<TH1F>("fDCosY", ";dy/ds", 200, -1., 1.);
     fDCosZ = tfs->make<TH1F>("fDCosZ", ";dz/ds", 200, -1., 1.);
 
     fMuMomentum = tfs->make<TH1F>("fMuMomentum", ";p_{#mu} (GeV/c)", 500, 0., 50.);
     fMuDCosX = tfs->make<TH1F>("fMuDCosX", ";dx/ds;", 200, -1., 1.);
     fMuDCosY = tfs->make<TH1F>("fMuDCosY", ";dy/ds;", 200, -1., 1.);
     fMuDCosZ = tfs->make<TH1F>("fMuDCosZ", ";dz/ds;", 200, -1., 1.);
 
     fEMomentum = tfs->make<TH1F>("fEMomentum", ";p_{e} (GeV/c)", 500, 0., 50.);
     fEDCosX = tfs->make<TH1F>("fEDCosX", ";dx/ds;", 200, -1., 1.);
     fEDCosY = tfs->make<TH1F>("fEDCosY", ";dy/ds;", 200, -1., 1.);
     fEDCosZ = tfs->make<TH1F>("fEDCosZ", ";dz/ds;", 200, -1., 1.);
 
     fCCMode = tfs->make<TH1F>("fCCMode", ";CC Interaction Mode;", 5, 0., 5.);
     fCCMode->GetXaxis()->SetBinLabel(1, "QE");
     fCCMode->GetXaxis()->SetBinLabel(2, "Res");
     fCCMode->GetXaxis()->SetBinLabel(3, "DIS");
     fCCMode->GetXaxis()->SetBinLabel(4, "Coh");
     fCCMode->GetXaxis()->SetBinLabel(5, "kInverseMuDecay");
     fCCMode->GetXaxis()->CenterLabels();
 
     fNCMode = tfs->make<TH1F>("fNCMode", ";NC Interaction Mode;", 5, 0., 5.);
     fNCMode->GetXaxis()->SetBinLabel(1, "QE");
     fNCMode->GetXaxis()->SetBinLabel(2, "Res");
     fNCMode->GetXaxis()->SetBinLabel(3, "DIS");
     fNCMode->GetXaxis()->SetBinLabel(4, "Coh");
     fNCMode->GetXaxis()->SetBinLabel(5, "kNuElectronElastic");
     fNCMode->GetXaxis()->CenterLabels();
 
     fECons = tfs->make<TH1F>("fECons", ";#Delta E(#nu,lepton);", 500, -5., 5.);
 
     art::ServiceHandle<geo::Geometry const> geo;
     double x = 2.1 * geo->DetHalfWidth();
     double y = 2.1 * geo->DetHalfHeight();
     double z = 2. * geo->DetLength();
     int xdiv = TMath::Nint(2 * x / 5.);
     int ydiv = TMath::Nint(2 * y / 5.);
     int zdiv = TMath::Nint(2 * z / 5.);
 
     fVertexX = tfs->make<TH1F>("fVertexX", ";x (cm)", xdiv, -x, x);
     fVertexY = tfs->make<TH1F>("fVertexY", ";y (cm)", ydiv, -y, y);
     fVertexZ = tfs->make<TH1F>("fVertexZ", ";z (cm)", zdiv, -0.2 * z, z);
 
     fVertexXY = tfs->make<TH2F>("fVertexXY", ";x (cm);y (cm)", xdiv, -x, x, ydiv, -y, y);
     fVertexXZ = tfs->make<TH2F>("fVertexXZ", ";z (cm);x (cm)", zdiv, -0.2 * z, z, xdiv, -x, x);
     fVertexYZ = tfs->make<TH2F>("fVertexYZ", ";z (cm);y (cm)", zdiv, -0.2 * z, z, ydiv, -y, y);
   }

void evgen::NDKGen::beginRun ( art::Run & run )

overrideprivatevirtual

Reimplemented from art::EDProducer.

Definition at line 184 of file NDKGen_module.cc.

References geo::GeometryCore::DetectorName(), art::fullRun(), and art::Run::put().

   {
     art::ServiceHandle<geo::Geometry const> geo;
     run.put(std::make_unique<sumdata::RunData>(geo->DetectorName()), art::fullRun());
   }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::consumes ( InputTag const & tag )

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

   {
     collector_.consumesMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::consumesView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const & tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

   {
     return collector_.consumesView<T, BT>(tag);
   }

void art::detail::Producer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Producer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

   {
     auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

   {
     auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

void art::detail::Producer::doEndJob ( )

inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Producer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto r = rp.makeRun(mc, rp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto sr = srp.makeSubRun(mc, srp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

bool art::detail::Producer::doEvent	(	EventPrincipal &	ep,
		ModuleContext const &	mc,
		std::atomic< std::size_t > &	counts_run,
		std::atomic< std::size_t > &	counts_passed,
		std::atomic< std::size_t > &	counts_failed
	)

inherited

Definition at line 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

   {
     auto e = ep.makeEvent(mc);
     ++counts_run;
     ProcessingFrame const frame{mc.scheduleID()};
     produceWithFrame(e, frame);
     e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
     ++counts_passed;
     return true;
   }

void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void evgen::NDKGen::endJob ( )

overrideprivatevirtual

Reimplemented from art::EDProducer.

Definition at line 191 of file NDKGen_module.cc.

References fEventFile.

   {
     fEventFile.close();
   }

void evgen::NDKGen::FillHistograms ( simb::MCTruth const & mc )

private

look for the outgoing lepton in the particle stack just interested in the first one

Definition at line 383 of file NDKGen_module.cc.

   {
     // Decide which histograms to put the spectrum in
     int id = -1;
     if (mc.GetNeutrino().CCNC() == simb::kCC) {
       fCCMode->Fill(mc.GetNeutrino().Mode());
       if (mc.GetNeutrino().Nu().PdgCode() == 12)
         id = 0;
       else if (mc.GetNeutrino().Nu().PdgCode() == -12)
         id = 1;
       else if (mc.GetNeutrino().Nu().PdgCode() == 14)
         id = 2;
       else if (mc.GetNeutrino().Nu().PdgCode() == -14)
         id = 3;
       else
         return;
     }
     else {
       fNCMode->Fill(mc.GetNeutrino().Mode());
       if (mc.GetNeutrino().Nu().PdgCode() > 0)
         id = 4;
       else
         id = 5;
     }
     if (id == -1) abort();
 
     // Fill the specta histograms
     fGenerated[id]->Fill(mc.GetNeutrino().Nu().E());
 
     //< fill the vertex histograms from the neutrino - that is always
     //< particle 0 in the list
     simb::MCNeutrino mcnu = mc.GetNeutrino();
     const simb::MCParticle nu = mcnu.Nu();
 
     fVertexX->Fill(nu.Vx());
     fVertexY->Fill(nu.Vy());
     fVertexZ->Fill(nu.Vz());
 
     fVertexXY->Fill(nu.Vx(), nu.Vy());
     fVertexXZ->Fill(nu.Vz(), nu.Vx());
     fVertexYZ->Fill(nu.Vz(), nu.Vy());
 
     double mom = nu.P();
     if (std::abs(mom) > 0.) {
       fDCosX->Fill(nu.Px() / mom);
       fDCosY->Fill(nu.Py() / mom);
       fDCosZ->Fill(nu.Pz() / mom);
     }
 
     std::cout << "REACTION:  " << ReactionChannel(mc.GetNeutrino().CCNC(), mc.GetNeutrino().Mode())
               << std::endl;
     std::cout << "-----------> Particles in the Stack = " << mc.NParticles() << std::endl;
     std::cout << std::setiosflags(std::ios::left) << std::setw(20) << "PARTICLE"
               << std::setiosflags(std::ios::left) << std::setw(32) << "STATUS" << std::setw(18)
               << "E (GeV)" << std::setw(18) << "m (GeV/c2)" << std::setw(18) << "Ek (GeV)"
               << std::endl
               << std::endl;
 
     const TDatabasePDG* databasePDG = TDatabasePDG::Instance();
 
     // Loop over the particle stack for this event
     for (int i = 0; i < mc.NParticles(); ++i) {
       simb::MCParticle part(mc.GetParticle(i));
       std::string name;
       if (part.PdgCode() == 18040)
         name = "Ar40 18040";
       else if (part.PdgCode() != -99999) {
         name = databasePDG->GetParticle(part.PdgCode())->GetName();
       }
 
       int code = part.StatusCode();
       std::string status = ParticleStatus(code);
       double mass = part.Mass();
       double energy = part.E();
       double Ek = (energy - mass); // Kinetic Energy (GeV)
 
       std::cout << std::setiosflags(std::ios::left) << std::setw(20) << name
                 << std::setiosflags(std::ios::left) << std::setw(32) << status << std::setw(18)
                 << energy << std::setw(18) << mass << std::setw(18) << Ek << std::endl;
     }
 
     if (mc.GetNeutrino().CCNC() == simb::kCC) {
 
       for (int i = 0; i < mc.NParticles(); ++i) {
         simb::MCParticle part(mc.GetParticle(i));
         if (abs(part.PdgCode()) == 11) {
           fEMomentum->Fill(part.P());
           fEDCosX->Fill(part.Px() / part.P());
           fEDCosY->Fill(part.Py() / part.P());
           fEDCosZ->Fill(part.Pz() / part.P());
           fECons->Fill(nu.E() - part.E());
           break;
         }
         else if (abs(part.PdgCode()) == 13) {
           fMuMomentum->Fill(part.P());
           fMuDCosX->Fill(part.Px() / part.P());
           fMuDCosY->Fill(part.Py() / part.P());
           fMuDCosZ->Fill(part.Pz() / part.P());
           fECons->Fill(nu.E() - part.E());
           break;
         }
       } // end loop over particles
     }   //end if CC interaction
   }

void art::Modifier::fillProductDescriptions ( )

inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

   {
     ProductRegistryHelper::fillDescriptions(moduleDescription());
   }

std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

   {
     return collector_.getConsumables();
   }

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const & wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

   {
     return doMakeWorker(wp);
   }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::mayConsume ( InputTag const & tag )

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

   {
     collector_.mayConsumeMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const & tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

std::string evgen::NDKGen::ParticleStatus ( int StatusCode )

private

Definition at line 342 of file NDKGen_module.cc.

Referenced by FillHistograms().

   {
     switch (StatusCode) {
     case 0: return "kIStInitialState";
     case 1: return "kIStFinalState";
     case 11: return "kIStNucleonTarget";
     default: return "Status Unknown";
     }
   }

void evgen::NDKGen::produce ( art::Event & evt )

overrideprivatevirtual

Implements art::EDProducer.

Definition at line 197 of file NDKGen_module.cc.

References simb::MCTruth::Add(), simb::MCParticle::AddTrajectoryPoint(), E, art::EventID::event(), fEngine, fEventFile, art::Event::id(), in, geo::GeometryCore::Iterate(), simb::kUnknown, simb::MCTruth::NParticles(), art::Event::put(), and simb::MCTruth::SetOrigin().

   {
 
     std::cout << std::endl;
     std::cout << "------------------------------------------------------------------------------"
               << std::endl;
     std::cout << "event  : " << evt.id().event() << std::endl;
 
     std::string name, k, dollar;
 
     // event dump format on file output by the two commands ....
     // gevgen_ndcy  -g 1000180400  -m 8 -n 400 -o ndk
     // gevdump -f ndk.1000.ghep.root > ndk.out
     std::string Name;
     int Idx, Ist, PDG, Mother1, Mother2, Daughter1, Daughter2;
     double Px, Py, Pz, E, m;
     std::string p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12, p13, p14;
 
     int trackid = -1; // set track id to -i as these are all primary particles and have id <= 0
     std::string primary("primary");
     int FirstMother = -1;
     double Mass = -9999;
     int Status = -9999;
 
     double P; // momentum of MCParticle IN GeV/c
 
     TLorentzVector Neutrino;
     TLorentzVector Lepton;
     TLorentzVector Target;
     TLorentzVector q;
     TLorentzVector Hadron4mom;
 
     TLorentzVector Tpos;
 
     std::unique_ptr<std::vector<simb::MCTruth>> truthcol(new std::vector<simb::MCTruth>);
     simb::MCTruth truth;
 
     art::ServiceHandle<geo::Geometry const> geo;
     CLHEP::RandFlat flat(fEngine);
 
     double const fvCut{5.0}; // force vtx to be this far from any wall.
 
     // Find boundary of active volume
     double minx = 1e9;
     double maxx = -1e9;
     double miny = 1e9;
     double maxy = -1e9;
     double minz = 1e9;
     double maxz = -1e9;
     for (auto const& tpc : geo->Iterate<geo::TPCGeo>(geo::CryostatID{0})) {
       auto const world = tpc.GetCenter();
       if (minx > world.X() - tpc.HalfWidth()) minx = world.X() - tpc.HalfWidth();
       if (maxx < world.X() + tpc.HalfWidth()) maxx = world.X() + tpc.HalfWidth();
       if (miny > world.Y() - tpc.HalfHeight()) miny = world.Y() - tpc.HalfHeight();
       if (maxy < world.Y() + tpc.HalfHeight()) maxy = world.Y() + tpc.HalfHeight();
       if (minz > world.Z() - tpc.Length() / 2.) minz = world.Z() - tpc.Length() / 2.;
       if (maxz < world.Z() + tpc.Length() / 2.) maxz = world.Z() + tpc.Length() / 2.;
     }
 
     // Assign vertice position
     double X0 = 0.0 + flat.fire(minx + fvCut, maxx - fvCut);
     double Y0 = 0.0 + flat.fire(miny + fvCut, maxy - fvCut);
     double Z0 = 0.0 + flat.fire(minz + fvCut, maxz - fvCut);
 
     std::cout << "NDKGen_module: X, Y, Z of vtx: " << X0 << ", " << Y0 << ", " << Z0 << std::endl;
 
     int GenieEvt = -999;
 
     if (!fEventFile.good()) std::cout << "NdkFile: Problem reading Ndk file" << std::endl;
 
     while (getline(fEventFile, k)) {
 
       if (k.find("** Event:") != std::string::npos) {
         std::istringstream in;
         in.clear();
         in.str(k);
         std::string dummy;
         in >> dummy >> dummy >> dummy >> dummy >> dummy >> dummy >> dummy >> dummy >> dummy >>
           dummy >> GenieEvt;
         std::cout << "Genie Evt " << GenieEvt << " art evt " << evt.id().event() << "\n";
       }
 
       if (GenieEvt + 1 != static_cast<int>(evt.id().event()))
         continue;
       else {
 
         if (!k.compare(0, 25, "GENIE Interaction Summary")) // testing for new event.
           break;
         if (k.compare(0, 1, "|") || k.compare(1, 2, "  "))
           continue; // uninteresting line if it doesn't start with "|" and if second and third characters aren't spaces.
         if (k.find("Fin-Init") != std::string::npos) continue; // Meh.
         if (k.find("Ar") != std::string::npos) continue;       // Meh.
         if (k.find("Cl") != std::string::npos) continue; // ignore chlorine nucleus in nnbar events
         if (k.find("HadrBlob") != std::string::npos) continue; // Meh.
         if (k.find("NucBindE") != std::string::npos) continue; // Meh. atmo
         if (k.find("FLAGS") != std::string::npos) break;       // Event end. gevgen_ndcy
         if (k.find("Vertex") != std::string::npos) break;      // Event end. atmo
 
         if (!k.compare(26, 1, "1")) // New event or stable particles.
         {
 
           std::istringstream in;
           in.clear();
           in.str(k);
 
           in >> p1 >> Idx >> p2 >> Name >> p3 >> Ist >> p4 >> PDG >> p5 >> Mother1 >> p6 >>
             Mother2 >> p7 >> Daughter1 >> p8 >> Daughter2 >> p9 >> Px >> p10 >> Py >> p11 >> Pz >>
             p12 >> E >> p13 >> m >> p14;
           if (Ist != 1) continue;
 
           std::cout << "PDG = " << PDG << std::endl;
 
           const TDatabasePDG* databasePDG = TDatabasePDG::Instance();
           const TParticlePDG* definition = databasePDG->GetParticle(PDG);
           Mass = definition->Mass();      // GeV
           if (E - Mass < 0.001) continue; // KE is too low.
 
           Status = 1;
 
           simb::MCParticle mcpart(trackid, PDG, primary, FirstMother, Mass, Status);
 
           P = std::sqrt(pow(E, 2.) - pow(Mass, 2.)); // GeV/c
           std::cout << "Momentum = " << P << std::endl;
 
           TLorentzVector pos(X0, Y0, Z0, 0);
 
           Tpos = pos; // for target
 
           TLorentzVector mom(Px, Py, Pz, E);
 
           mcpart.AddTrajectoryPoint(pos, mom);
           truth.Add(mcpart);
 
         } // loop over particles in an event
         truth.SetOrigin(simb::kUnknown);
       }
     } // end while loop
 
     std::cout << "NDKGen.cxx: Putting " << truth.NParticles() << " tracks on stack." << std::endl;
     truthcol->push_back(truth);
     evt.put(std::move(truthcol));
   }

std::string evgen::NDKGen::ReactionChannel	(	int	ccnc,
		int	mode
	)

private

Definition at line 353 of file NDKGen_module.cc.

Referenced by FillHistograms().

   {
     std::string ReactionChannelName = " ";
 
     if (ccnc == 0)
       ReactionChannelName = "kCC";
     else if (ccnc == 1)
       ReactionChannelName = "kNC";
     else
       std::cout << "Current mode unknown!! " << std::endl;
 
     if (mode == 0)
       ReactionChannelName += "_kQE";
     else if (mode == 1)
       ReactionChannelName += "_kRes";
     else if (mode == 2)
       ReactionChannelName += "_kDIS";
     else if (mode == 3)
       ReactionChannelName += "_kCoh";
     else if (mode == 4)
       ReactionChannelName += "_kNuElectronElastic";
     else if (mode == 5)
       ReactionChannelName += "_kInverseMuDecay";
     else
       std::cout << "interaction mode unknown!! " << std::endl;
 
     return ReactionChannelName;
   }

void art::Modifier::registerProducts ( ProductDescriptions & productsToRegister )

inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

   {
     ProductRegistryHelper::registerProducts(productsToRegister,
                                             moduleDescription());
   }

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void art::ModuleBase::sortConsumables ( std::string const & current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }