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

Inheritance diagram for evgen::CosmicsGen:

Public Types
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

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

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

virtual	~CosmicsGen ()

void	produce (art::Event &evt)

void	beginJob ()

void	beginRun (art::Run &run)

void	reconfigure (fhicl::ParameterSet const &p)

template<typename PROD , BranchType B = InEvent>
ProductID	getProductID (std::string const &instanceName={}) const

template<typename PROD , BranchType B>
ProductID	getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const

bool	modifiesEvent () const

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

template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)

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

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

template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)

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

template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)

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

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

template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)

base_engine_t &	createEngine (seed_t seed)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)

seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const

void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)

void	prepareForJob (fhicl::ParameterSet const &pset)

void	showMissingConsumes () const

Private Attributes
evgb::CRYHelper *	fCRYHelp
	CRY generator object. More...

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

Detailed Description

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

Definition at line 46 of file CosmicsGen_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 43 of file EDProducer.h.

using art::EDProducer::WorkerType = WorkerT<EDProducer>

inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 105 of file CosmicsGen_module.cc.

References art::EngineCreator::createEngine(), fbuffbox, fhicl::ParameterSet::get(), and reconfigure().

     : fCRYHelp(0)
   {
     // create a default random engine; obtain the random seed from NuRandomService,
     // unless overridden in configuration with key "Seed"
     art::ServiceHandle<rndm::NuRandomService>()
       ->createEngine(*this, pset, "Seed");
     
     //the buffer box bounds specified here will extend on the cryostat boundaries
     fbuffbox = pset.get< std::vector<double> >("BufferBox",{0.0, 0.0, 0.0, 0.0, 0.0, 0.0});
 
     this->reconfigure(pset);
     
     produces< std::vector<simb::MCTruth> >();
     produces< sumdata::RunData, art::InRun >();    
   }

evgen::CosmicsGen::~CosmicsGen ( )

virtual

Definition at line 123 of file CosmicsGen_module.cc.

References fCRYHelp.

   {
     if(fCRYHelp) delete fCRYHelp;
   }

Member Function Documentation

void evgen::CosmicsGen::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 148 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, fPhotonsPerSample, and art::TFileDirectory::make().

   {
     art::ServiceHandle<art::TFileService> 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 )

virtual

Reimplemented from art::EDProducer.

Definition at line 185 of file CosmicsGen_module.cc.

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

   {
 
     // grab the geometry object to see what geometry we are using
     art::ServiceHandle<geo::Geometry> geo;
 
     std::unique_ptr<sumdata::RunData> runcol(new sumdata::RunData(geo->DetectorName()));
 
     run.put(std::move(runcol));
 
     return;
   }

template<typename T , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const & it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const & it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

 {
   return rng()->createEngine(placeholder_schedule_id(), seed);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make
	)

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make,
		label_t const &	engine_label
	)

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
 }

CurrentProcessingContext const * art::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

   {
     return current_context_.get();
   }

EngineCreator::seed_t EngineCreator::get_seed_value	(	fhicl::ParameterSet const &	pset,
		char const	key[] = `"seed"`,
		seed_t const	implicit_seed = `-1`
	)

inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

 {
   auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
   return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
 }

template<typename PROD , BranchType B>

ProductID art::EDProducer::getProductID ( std::string const & instanceName = {} ) const

inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

   {
     return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                                instanceName);
   }

template<typename PROD , BranchType B>

ProductID art::ProducerBase::getProductID	(	ModuleDescription const &	moduleDescription,
		std::string const &	instanceName
	)		const

inherited

Definition at line 56 of file ProducerBase.h.

References B, and art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

   {
     auto const& pd =
       get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
     return pd.productID();
   }

template<typename T , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const & it )

inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const & it )

inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

     {
       return true;
     }

cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )

staticinherited

Definition at line 76 of file Consumer.cc.

 {
   static Consumer invalid{InvalidTag{}};
   return &invalid;
 }

void art::Consumer::prepareForJob ( fhicl::ParameterSet const & pset )

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

 {
   if (!moduleContext_)
     return;
 
   pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
   for (auto& consumablesPerBranch : consumables_) {
     cet::sort_all(consumablesPerBranch);
   }
 }

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

virtual

Implements art::EDProducer.

Definition at line 199 of file CosmicsGen_module.cc.

References simb::MCTruth::Add(), bounds(), geo::GeometryCore::CryostatBoundaries(), 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(), simb::kCosmicRay, simb::MCParticle::Momentum(), geo::GeometryCore::Ncryostats(), 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> 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;
       // for c2: remove unused variables
       //int tpcPhotons   = 0;
       //int tpcElectrons = 0;
       //int tpcMuons     = 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     = 0;
         TH2F* hAngles   = 0;
         TH2F* hAnglesLo = 0;
         TH2F* hAnglesMi = 0;
         TH2F* hAnglesHi = 0;
         TH1F* hEnergy   = 0;
         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(unsigned int c = 0; c < geom->Ncryostats(); ++c){
         
           double bounds[6] = {0.};
           geom->CryostatBoundaries(bounds, c);
           
           //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;
 
             //The following code no longer works now that we require intersection with the cryostat boundary
             //For example, the particle could intersect this cryostat but miss its TPC, but intersect a TPC 
             //in another cryostat
             /*try{
               unsigned int tpc   = 0;
               unsigned int cstat = 0;
               geom->PositionToTPC(x2, tpc, cstat);
               if      (std::abs(particle.PdgCode())==13) ++tpcMuons;
               else if (std::abs(particle.PdgCode())==22) ++tpcPhotons;
               else if (std::abs(particle.PdgCode())==11) ++tpcElectrons;
             }
             catch(cet::exception &e){
               LOG_DEBUG("CosmicsGen") << "current particle does not go through any tpc";
             }*///
 
             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);
 
       /*fPhotonsInTPC  ->Fill(tpcPhotons);
       fElectronsInTPC->Fill(tpcElectrons);
       fMuonsInTPC    ->Fill(tpcMuons);*/
     }
 
     truthcol->push_back(truth);
     evt.put(std::move(truthcol));
   
     return;
   }// end produce

void evgen::CosmicsGen::reconfigure ( fhicl::ParameterSet const & p )

Definition at line 129 of file CosmicsGen_module.cc.

References fCRYHelp, art::RandomNumberGenerator::getEngine(), geo::GeometryCore::GetWorldVolumeName(), and art::EngineCreator::rng().

Referenced by CosmicsGen().

   {
     if(fCRYHelp){
       delete fCRYHelp; 
       fCRYHelp = 0;
     }
 
     // get the random number generator service and make some CLHEP generators
     art::ServiceHandle<art::RandomNumberGenerator> rng;
     CLHEP::HepRandomEngine& engine = rng->getEngine();
 
     art::ServiceHandle<geo::Geometry> geo;
 
     fCRYHelp = new evgb::CRYHelper(p, engine, geo->GetWorldVolumeName());
 
     return;
   }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

 {
   if (!moduleContext_)
     return;
 
   // If none of the branches have missing consumes statements, exit early.
   if (std::all_of(cbegin(missingConsumes_),
                   cend(missingConsumes_),
                   [](auto const& perBranch) { return perBranch.empty(); }))
     return;
 
   constexpr cet::HorizontalRule rule{60};
   mf::LogPrint log{"MTdiagnostics"};
   log << '\n'
       << rule('=') << '\n'
       << "The following consumes (or mayConsume) statements are missing from\n"
       << module_context(moduleDescription_) << '\n'
       << rule('-') << '\n';
 
   cet::for_all_with_index(
     missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
       for (auto const& pi : perBranch) {
         log << "  "
             << assemble_consumes_statement(static_cast<BranchType>(i), pi)
             << '\n';
       }
     });
   log << rule('=');
 }

void art::Consumer::validateConsumedProduct	(	BranchType const	bt,
		ProductInfo const &	pi
	)

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

 {
   // Early exits if consumes tracking has been disabled or if the
   // consumed product is an allowed consumable.
   if (!moduleContext_)
     return;
 
   if (cet::binary_search_all(consumables_[bt], pi))
     return;
 
   if (requireConsumes_) {
     throw Exception(errors::ProductRegistrationFailure,
                     "Consumer: an error occurred during validation of a "
                     "retrieved product\n\n")
       << "The following consumes (or mayConsume) statement is missing from\n"
       << module_context(moduleDescription_) << ":\n\n"
       << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
   }
 
   missingConsumes_[bt].insert(pi);
 }

Member Data Documentation

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

private

Definition at line 61 of file CosmicsGen_module.cc.

Referenced by CosmicsGen(), and produce().

evgb::CRYHelper* evgen::CosmicsGen::fCRYHelp

private

CRY generator object.

Definition at line 59 of file CosmicsGen_module.cc.

Referenced by produce(), reconfigure(), and ~CosmicsGen().

TH2F* evgen::CosmicsGen::fElectronAngles

private

Electron rate vs angle.

Definition at line 75 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fElectronAnglesHi

private

Electron rate vs angle, high momenta.

Definition at line 78 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fElectronAnglesLo

private

Electron rate vs angle, low momenta.

Definition at line 76 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fElectronAnglesMi

private

Electron rate vs angle, middle momenta.

Definition at line 77 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fElectronCosQ

private

Electron rate vs cos(Q)

Definition at line 79 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fElectronEnergy

private

Electron energy (GeV)

Definition at line 80 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 82 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 84 of file CosmicsGen_module.cc.

Referenced by beginJob().

TH1F* evgen::CosmicsGen::fElectronsPerSample

private

number of electrons in the sampled time window

Definition at line 81 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fMuonAngles

private

Muon rate vs angle.

Definition at line 87 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fMuonAnglesHi

private

Muon rate vs angle, high momenta.

Definition at line 90 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fMuonAnglesLo

private

Muon rate vs angle, low momenta.

Definition at line 88 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fMuonAnglesMi

private

Muon rate vs angle, middle momenta.

Definition at line 89 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fMuonCosQ

private

Muon rate vs cos(Q)

Definition at line 91 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fMuonEnergy

private

Muon energy (GeV)

Definition at line 92 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 94 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 96 of file CosmicsGen_module.cc.

Referenced by beginJob().

TH1F* evgen::CosmicsGen::fMuonsPerSample

private

number of muons in the sampled time window

Definition at line 93 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fPhotonAngles

private

Photon rate vs angle.

Definition at line 63 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fPhotonAnglesHi

private

Photon rate vs angle, high momenta.

Definition at line 66 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fPhotonAnglesLo

private

Photon rate vs angle, low momenta.

Definition at line 64 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH2F* evgen::CosmicsGen::fPhotonAnglesMi

private

Photon rate vs angle, middle momenta.

Definition at line 65 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fPhotonCosQ

private

Photon rate vs cos(Q)

Definition at line 67 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

TH1F* evgen::CosmicsGen::fPhotonEnergy

private

Photon energy (GeV)

Definition at line 68 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 70 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 72 of file CosmicsGen_module.cc.

Referenced by beginJob().

TH1F* evgen::CosmicsGen::fPhotonsPerSample

private

number of photons in the sampled time window

Definition at line 69 of file CosmicsGen_module.cc.

Referenced by beginJob(), and produce().

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

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

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation