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

Inheritance diagram for evgen::CosmicsGen:

Public Types
using	ModuleType = EDProducer

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

Public Member Functions
	CosmicsGen (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

Private Attributes
std::vector< double >	fbuffbox

TH2F *	fPhotonAngles
	Photon rate vs angle. More...

TH2F *	fPhotonAnglesLo
	Photon rate vs angle, low momenta. More...

TH2F *	fPhotonAnglesMi
	Photon rate vs angle, middle momenta. More...

TH2F *	fPhotonAnglesHi
	Photon rate vs angle, high momenta. More...

TH1F *	fPhotonCosQ
	Photon rate vs cos(Q) More...

TH1F *	fPhotonEnergy
	Photon energy (GeV) More...

TH1F *	fPhotonsPerSample
	number of photons in the sampled time window More...

TH1F *	fPhotonsInCStat

TH1F *	fPhotonsInTPC

TH2F *	fElectronAngles
	Electron rate vs angle. More...

TH2F *	fElectronAnglesLo
	Electron rate vs angle, low momenta. More...

TH2F *	fElectronAnglesMi
	Electron rate vs angle, middle momenta. More...

TH2F *	fElectronAnglesHi
	Electron rate vs angle, high momenta. More...

TH1F *	fElectronCosQ
	Electron rate vs cos(Q) More...

TH1F *	fElectronEnergy
	Electron energy (GeV) More...

TH1F *	fElectronsPerSample
	number of electrons in the sampled time window More...

TH1F *	fElectronsInCStat

TH1F *	fElectronsInTPC

TH2F *	fMuonAngles
	Muon rate vs angle. More...

TH2F *	fMuonAnglesLo
	Muon rate vs angle, low momenta. More...

TH2F *	fMuonAnglesMi
	Muon rate vs angle, middle momenta. More...

TH2F *	fMuonAnglesHi
	Muon rate vs angle, high momenta. More...

TH1F *	fMuonCosQ
	Muon rate vs cos(Q) More...

TH1F *	fMuonEnergy
	Muon energy (GeV) More...

TH1F *	fMuonsPerSample
	number of muons in the sampled time window More...

TH1F *	fMuonsInCStat

TH1F *	fMuonsInTPC

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

evgb::CRYHelper	fCRYHelp
	CRY generator object. More...

Detailed Description

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

Definition at line 37 of file CosmicsGen_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::CosmicsGen::CosmicsGen ( fhicl::ParameterSet const & pset )

explicit

Definition at line 91 of file CosmicsGen_module.cc.

References art::detail::EngineCreator::createEngine(), fbuffbox, fCRYHelp, and fEngine.

     : art::EDProducer{pset}
     , fbuffbox{pset.get<std::vector<double>>("BufferBox", {0.0, 0.0, 0.0, 0.0, 0.0, 0.0})}
     // 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"))
     , fCRYHelp{pset, fEngine, art::ServiceHandle<geo::Geometry const>{} -> GetWorldVolumeName()}
   {
     produces<std::vector<simb::MCTruth>>();
     produces<sumdata::RunData, art::InRun>();
   }

Member Function Documentation

void evgen::CosmicsGen::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDProducer.

Definition at line 106 of file CosmicsGen_module.cc.

References fElectronAngles, fElectronAnglesHi, fElectronAnglesLo, fElectronAnglesMi, fElectronCosQ, fElectronEnergy, fElectronsInCStat, fElectronsInTPC, fElectronsPerSample, fMuonAngles, fMuonAnglesHi, fMuonAnglesLo, fMuonAnglesMi, fMuonCosQ, fMuonEnergy, fMuonsInCStat, fMuonsInTPC, fMuonsPerSample, fPhotonAngles, fPhotonAnglesHi, fPhotonAnglesLo, fPhotonAnglesMi, fPhotonCosQ, fPhotonEnergy, fPhotonsInCStat, fPhotonsInTPC, and fPhotonsPerSample.

   {
     art::ServiceHandle<art::TFileService const> tfs;
 
     fPhotonAngles =
       tfs->make<TH2F>("fPhotonAngles", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fPhotonAnglesLo =
       tfs->make<TH2F>("fPhotonAnglesLo", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fPhotonAnglesMi =
       tfs->make<TH2F>("fPhotonAnglesMi", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fPhotonAnglesHi =
       tfs->make<TH2F>("fPhotonAnglesHi", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fPhotonCosQ = tfs->make<TH1F>("fPhotonCosQ", ";cos#theta;tracks", 50, -1.0, 1.0);
     fPhotonEnergy = tfs->make<TH1F>("fPhotonEnergy", ";E (GeV)", 5000, 0.0, 1000.0);
     fPhotonsPerSample =
       tfs->make<TH1F>("fPhotonsPerSample", ";Number Photons;Samples", 100, 0, 1000);
     fPhotonsInCStat =
       tfs->make<TH1F>("fPhotonsInCryostat", ";Number Photons;Samples", 100, 0, 1000);
     fPhotonsInTPC = tfs->make<TH1F>("fPhotonsInTPC", ";Number Photons;Samples", 100, 0, 1000);
 
     fElectronAngles =
       tfs->make<TH2F>("fElectronAngles", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fElectronAnglesLo =
       tfs->make<TH2F>("fElectronAnglesLo", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fElectronAnglesMi =
       tfs->make<TH2F>("fElectronAnglesMi", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fElectronAnglesHi =
       tfs->make<TH2F>("fElectronAnglesHi", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fElectronCosQ = tfs->make<TH1F>("fElectronCosQ", ";cos#theta;tracks", 50, -1.0, 1.0);
     fElectronEnergy = tfs->make<TH1F>("fElectronEnergy", ";E (GeV)", 5000, 0.0, 1000.0);
     fElectronsPerSample =
       tfs->make<TH1F>("fElectronsPerSample", ";Number Electrons;Samples", 100, 0, 1000);
     fElectronsInCStat =
       tfs->make<TH1F>("fElectronsInCryotat", ";Number Electrons;Samples", 100, 0, 1000);
     fElectronsInTPC = tfs->make<TH1F>("fElectronsInTPC", ";Number Electrons;Samples", 100, 0, 1000);
 
     fMuonAngles =
       tfs->make<TH2F>("fMuonAngles", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fMuonAnglesLo =
       tfs->make<TH2F>("fMuonAnglesLo", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fMuonAnglesMi =
       tfs->make<TH2F>("fMuonAnglesMi", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fMuonAnglesHi =
       tfs->make<TH2F>("fMuonAnglesHi", ";#phi;cos#theta", 36, -180.0, 180.0, 50, -1.0, 1.0);
     fMuonCosQ = tfs->make<TH1F>("fMuonCosQ", ";cos#theta;tracks", 50, -1.0, 1.0);
     fMuonEnergy = tfs->make<TH1F>("fMuonEnergy", ";E (GeV)", 5000, 0.0, 1000.0);
     fMuonsPerSample = tfs->make<TH1F>("fMuonsPerSample", ";Number Muons;Samples", 100, 0, 1000);
     fMuonsInCStat = tfs->make<TH1F>("fMuonsInCryostat", ";Number Muons;Samples", 100, 0, 1000);
     fMuonsInTPC = tfs->make<TH1F>("fMuonsInTPC", ";Number Muons;Samples", 100, 0, 1000);
   }

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

overrideprivatevirtual

Reimplemented from art::EDProducer.

Definition at line 158 of file CosmicsGen_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 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";
   }

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

overrideprivatevirtual

Implements art::EDProducer.

Definition at line 165 of file CosmicsGen_module.cc.

References util::abs(), simb::MCTruth::Add(), DEFINE_ART_MODULE, geo::GeometryCore::DetLength(), fbuffbox, fCRYHelp, fElectronAngles, fElectronAnglesHi, fElectronAnglesLo, fElectronAnglesMi, fElectronCosQ, fElectronEnergy, fElectronsInCStat, fElectronsPerSample, fMuonAngles, fMuonAnglesHi, fMuonAnglesLo, fMuonAnglesMi, fMuonCosQ, fMuonEnergy, fMuonsInCStat, fMuonsPerSample, fPhotonAngles, fPhotonAnglesHi, fPhotonAnglesLo, fPhotonAnglesMi, fPhotonCosQ, fPhotonEnergy, fPhotonsInCStat, fPhotonsPerSample, simb::MCTruth::GetParticle(), geo::GeometryCore::Iterate(), simb::kCosmicRay, simb::MCParticle::Momentum(), simb::MCTruth::NParticles(), simb::MCParticle::PdgCode(), simb::MCParticle::Position(), art::Event::put(), evgb::CRYHelper::Sample(), simb::MCTruth::SetOrigin(), and geo::GeometryCore::SurfaceY().

   {
     std::unique_ptr<std::vector<simb::MCTruth>> truthcol(new std::vector<simb::MCTruth>);
 
     // fill some histograms about this event
     art::ServiceHandle<geo::Geometry const> geom;
 
     int nCrossCryostat = 0;
 
     simb::MCTruth truth;
 
     while (nCrossCryostat < 1) {
 
       simb::MCTruth pretruth;
       truth.SetOrigin(simb::kCosmicRay);
       fCRYHelp.Sample(pretruth, geom->SurfaceY(), geom->DetLength(), 0);
 
       int numPhotons = 0;
       int numElectrons = 0;
       int numMuons = 0;
       int allPhotons = 0;
       int allElectrons = 0;
       int allMuons = 0;
 
       // loop over particles in the truth object
       for (int i = 0; i < pretruth.NParticles(); ++i) {
         simb::MCParticle particle = pretruth.GetParticle(i);
         const TLorentzVector& v4 = particle.Position();
         const TLorentzVector& p4 = particle.Momentum();
         double x0[3] = {v4.X(), v4.Y(), v4.Z()};
         double dx[3] = {p4.Px(), p4.Py(), p4.Pz()};
 
         if (std::abs(particle.PdgCode()) == 13)
           ++allMuons;
         else if (std::abs(particle.PdgCode()) == 22)
           ++allPhotons;
         else if (std::abs(particle.PdgCode()) == 11)
           ++allElectrons;
 
         TH1F* hCosQ = nullptr;
         TH2F* hAngles = nullptr;
         TH2F* hAnglesLo = nullptr;
         TH2F* hAnglesMi = nullptr;
         TH2F* hAnglesHi = nullptr;
         TH1F* hEnergy = nullptr;
         if (std::abs(particle.PdgCode()) == 13) {
           hCosQ = fMuonCosQ;
           hAngles = fMuonAngles;
           hAnglesLo = fMuonAnglesLo;
           hAnglesMi = fMuonAnglesMi;
           hAnglesHi = fMuonAnglesHi;
           hEnergy = fMuonEnergy;
         }
         else if (std::abs(particle.PdgCode()) == 22) {
           hCosQ = fPhotonCosQ;
           hAngles = fPhotonAngles;
           hAnglesLo = fPhotonAnglesLo;
           hAnglesMi = fPhotonAnglesMi;
           hAnglesHi = fPhotonAnglesHi;
           hEnergy = fPhotonEnergy;
         }
         else if (std::abs(particle.PdgCode()) == 11) {
           hCosQ = fElectronCosQ;
           hAngles = fElectronAngles;
           hAnglesLo = fElectronAnglesLo;
           hAnglesMi = fElectronAnglesMi;
           hAnglesHi = fElectronAnglesHi;
           hEnergy = fElectronEnergy;
         }
 
         // now check if the particle goes through any cryostat in the detector
         // if so, add it to the truth object.
         for (auto const& cryostat : geom->Iterate<geo::CryostatGeo>()) {
 
           double bounds[6] = {0.};
           cryostat.Boundaries(bounds);
 
           //add a buffer box around the cryostat bounds to increase the acceptance
           //(geometrically) at the CRY level to make up for particles we will loose
           //due to multiple scattering effects that pitch in during GEANT4 tracking
           //By default, the buffer box has zero size
           for (unsigned int cb = 0; cb < 6; cb++)
             bounds[cb] = bounds[cb] + fbuffbox[cb];
 
           //calculate the intersection point with each cryostat surface
           bool intersects_cryo = false;
           for (int bnd = 0; bnd != 6; ++bnd) {
             if (bnd < 2) {
               double p2[3] = {bounds[bnd],
                               x0[1] + (dx[1] / dx[0]) * (bounds[bnd] - x0[0]),
                               x0[2] + (dx[2] / dx[0]) * (bounds[bnd] - x0[0])};
               if (p2[1] >= bounds[2] && p2[1] <= bounds[3] && p2[2] >= bounds[4] &&
                   p2[2] <= bounds[5]) {
                 intersects_cryo = true;
                 break;
               }
             }
             else if (bnd >= 2 && bnd < 4) {
               double p2[3] = {x0[0] + (dx[0] / dx[1]) * (bounds[bnd] - x0[1]),
                               bounds[bnd],
                               x0[2] + (dx[2] / dx[1]) * (bounds[bnd] - x0[1])};
               if (p2[0] >= bounds[0] && p2[0] <= bounds[1] && p2[2] >= bounds[4] &&
                   p2[2] <= bounds[5]) {
                 intersects_cryo = true;
                 break;
               }
             }
             else if (bnd >= 4) {
               double p2[3] = {x0[0] + (dx[0] / dx[2]) * (bounds[bnd] - x0[2]),
                               x0[1] + (dx[1] / dx[2]) * (bounds[bnd] - x0[2]),
                               bounds[bnd]};
               if (p2[0] >= bounds[0] && p2[0] <= bounds[1] && p2[1] >= bounds[2] &&
                   p2[1] <= bounds[3]) {
                 intersects_cryo = true;
                 break;
               }
             }
           }
 
           if (intersects_cryo) {
             truth.Add(particle);
 
             if (std::abs(particle.PdgCode()) == 13)
               ++numMuons;
             else if (std::abs(particle.PdgCode()) == 22)
               ++numPhotons;
             else if (std::abs(particle.PdgCode()) == 11)
               ++numElectrons;
 
             if (hCosQ != 0) {
               double cosq = -p4.Py() / p4.P();
               double phi = std::atan2(p4.Pz(), p4.Px());
               phi *= 180 / M_PI;
               hCosQ->Fill(cosq);
               hAngles->Fill(phi, cosq);
               if (p4.E() < 1.0)
                 hAnglesLo->Fill(phi, cosq);
               else if (p4.E() < 10.0)
                 hAnglesMi->Fill(phi, cosq);
               else
                 hAnglesHi->Fill(phi, cosq);
               hEnergy->Fill(p4.E());
             }      //end if there is a cos(theta) histogram
             break; //leave loop over cryostats to avoid adding particle multiple times
           }        // end if particle goes into a cryostat
         }          // end loop over cryostats in the detector
 
       } // loop on particles
 
       nCrossCryostat = truth.NParticles();
 
       fPhotonsPerSample->Fill(allPhotons);
       fElectronsPerSample->Fill(allElectrons);
       fMuonsPerSample->Fill(allMuons);
 
       fPhotonsInCStat->Fill(numPhotons);
       fElectronsInCStat->Fill(numElectrons);
       fMuonsInCStat->Fill(numMuons);
     }
 
     truthcol->push_back(truth);
     evt.put(std::move(truthcol));
   } // end produce

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);
   }

Member Data Documentation

std::vector<double> evgen::CosmicsGen::fbuffbox

private

Definition at line 46 of file CosmicsGen_module.cc.

Referenced by CosmicsGen(), and produce().

evgb::CRYHelper evgen::CosmicsGen::fCRYHelp

private

CRY generator object.

Definition at line 84 of file CosmicsGen_module.cc.

Referenced by CosmicsGen(), and produce().

TH2F* evgen::CosmicsGen::fElectronAngles

private

Electron rate vs angle.

Definition at line 60 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fElectronAnglesHi

private

Electron rate vs angle, high momenta.

Definition at line 63 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fElectronAnglesLo

private

Electron rate vs angle, low momenta.

Definition at line 61 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fElectronAnglesMi

private

Electron rate vs angle, middle momenta.

Definition at line 62 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fElectronCosQ

private

Electron rate vs cos(Q)

Definition at line 64 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fElectronEnergy

private

Electron energy (GeV)

Definition at line 65 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fElectronsInCStat

private

number of electrons in the cryostat during the sampled time window

Definition at line 67 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fElectronsInTPC

private

number of electrons in the tpc during the sampled time window

Definition at line 69 of file CosmicsGen_module.cc.

Referenced by beginJob().

TH1F* evgen::CosmicsGen::fElectronsPerSample

private

number of electrons in the sampled time window

Definition at line 66 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

CLHEP::HepRandomEngine& evgen::CosmicsGen::fEngine

private

art-managed random-number engine

Definition at line 83 of file CosmicsGen_module.cc.

Referenced by CosmicsGen().

TH2F* evgen::CosmicsGen::fMuonAngles

private

Muon rate vs angle.

Definition at line 72 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fMuonAnglesHi

private

Muon rate vs angle, high momenta.

Definition at line 75 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fMuonAnglesLo

private

Muon rate vs angle, low momenta.

Definition at line 73 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fMuonAnglesMi

private

Muon rate vs angle, middle momenta.

Definition at line 74 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fMuonCosQ

private

Muon rate vs cos(Q)

Definition at line 76 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fMuonEnergy

private

Muon energy (GeV)

Definition at line 77 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fMuonsInCStat

private

number of muons in the cryostat during the sampled time window

Definition at line 79 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fMuonsInTPC

private

number of muons in the tpc during the sampled time window

Definition at line 81 of file CosmicsGen_module.cc.

Referenced by beginJob().

TH1F* evgen::CosmicsGen::fMuonsPerSample

private

number of muons in the sampled time window

Definition at line 78 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fPhotonAngles

private

Photon rate vs angle.

Definition at line 48 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fPhotonAnglesHi

private

Photon rate vs angle, high momenta.

Definition at line 51 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fPhotonAnglesLo

private

Photon rate vs angle, low momenta.

Definition at line 49 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fPhotonAnglesMi

private

Photon rate vs angle, middle momenta.

Definition at line 50 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fPhotonCosQ

private

Photon rate vs cos(Q)

Definition at line 52 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fPhotonEnergy

private

Photon energy (GeV)

Definition at line 53 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fPhotonsInCStat

private

number of photons in the cryostat during the sampled time window

Definition at line 55 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fPhotonsInTPC

private

number of photons in the tpc during the sampled time window

Definition at line 57 of file CosmicsGen_module.cc.

Referenced by beginJob().

TH1F* evgen::CosmicsGen::fPhotonsPerSample

private

number of photons in the sampled time window

Definition at line 54 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

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

larsim/v09_43_00/source/larsim/EventGenerator/CRY/CosmicsGen_module.cc

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation