LArSoft  v06_85_00
Liquid Argon Software toolkit - http://larsoft.org/
evgen::NucleonDecay Class Reference
Inheritance diagram for evgen::NucleonDecay:
art::EDProducer art::ProducerBase art::Consumer art::EngineCreator art::ProductRegistryHelper

Public Types

using ModuleType = EDProducer
 
using WorkerType = WorkerT< EDProducer >
 
template<typename UserConfig , typename KeysToIgnore = void>
using Table = ProducerBase::Table< UserConfig, KeysToIgnore >
 

Public Member Functions

 NucleonDecay (fhicl::ParameterSet const &p)
 
 NucleonDecay (NucleonDecay const &)=delete
 
 NucleonDecay (NucleonDecay &&)=delete
 
NucleonDecayoperator= (NucleonDecay const &)=delete
 
NucleonDecayoperator= (NucleonDecay &&)=delete
 
void produce (art::Event &e) override
 
void beginJob () override
 
void beginRun (art::Run &run) override
 
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_tcreateEngine (seed_t seed)
 
base_engine_tcreateEngine (seed_t seed, std::string const &kind_of_engine_to_make)
 
base_engine_tcreateEngine (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< Consumernon_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

const genie::EventRecordVisitorI * mcgen
 
genie::NucleonDecayMode_t gOptDecayMode = genie::kNDNull
 
int dpdg = 0
 

Detailed Description

Definition at line 52 of file NucleonDecay_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::NucleonDecay::NucleonDecay ( fhicl::ParameterSet const &  p)
explicit

Definition at line 80 of file NucleonDecay_module.cc.

References art::EngineCreator::createEngine(), dpdg, fhicl::ParameterSet::get(), gOptDecayMode, mcgen, and seed.

81 {
82  genie::PDGLibrary::Instance(); //Ensure Messenger is started first in GENIE.
83 
84  string sname = "genie::EventGenerator";
85  string sconfig = "NucleonDecay";
86  genie::AlgFactory * algf = genie::AlgFactory::Instance();
87  mcgen =
88  dynamic_cast<const genie::EventRecordVisitorI *> (algf->GetAlgorithm(sname,sconfig));
89  if(!mcgen) {
90  throw cet::exception("NucleonDecay") << "Couldn't instantiate the nucleon decay generator";
91  }
92  int fDecayMode = p.get<int>("DecayMode");
93  gOptDecayMode = (genie::NucleonDecayMode_t) fDecayMode;
94 
95  if (p.get<int>("DecayedNucleon") > 0 ){
96  dpdg = p.get<int>("DecayedNucleon");
97  }
98  else{
99  dpdg = genie::utils::nucleon_decay::DecayedNucleonPdgCode(gOptDecayMode);
100  }
101 
102  produces< std::vector<simb::MCTruth> >();
103  produces< sumdata::RunData, art::InRun >();
104 
105  // create a default random engine; obtain the random seed from NuRandomService,
106  // unless overridden in configuration with key "Seed"
108  ->createEngine(*this, p, "Seed");
109 
110  unsigned int seed = art::ServiceHandle<rndm::NuRandomService>()->getSeed();
111  genie::utils::app_init::RandGen(seed);
112 }
genie::NucleonDecayMode_t gOptDecayMode
const genie::EventRecordVisitorI * mcgen
base_engine_t & createEngine(seed_t seed)
long seed
Definition: chem4.cc:68
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
evgen::NucleonDecay::NucleonDecay ( NucleonDecay const &  )
delete
evgen::NucleonDecay::NucleonDecay ( NucleonDecay &&  )
delete

Member Function Documentation

void evgen::NucleonDecay::beginJob ( )
overridevirtual

Reimplemented from art::EDProducer.

Definition at line 217 of file NucleonDecay_module.cc.

References DEFINE_ART_MODULE.

218 {
219  // Implementation of optional member function here.
220 }
void evgen::NucleonDecay::beginRun ( art::Run run)
overridevirtual

Reimplemented from art::EDProducer.

Definition at line 204 of file NucleonDecay_module.cc.

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

205 {
206 
207  // grab the geometry object to see what geometry we are using
209  std::unique_ptr<sumdata::RunData> runcol(new sumdata::RunData(geo->DetectorName()));
210 
211  run.put(std::move(runcol));
212 
213  return;
214 }
art::ProductID put(std::unique_ptr< PROD > &&)
Definition: Run.h:148
std::string DetectorName() const
Returns a string with the name of the detector, as configured.
Namespace collecting geometry-related classes utilities.
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().

148 {
149  if (!moduleContext_)
150  return ProductToken<T>::invalid();
151 
152  consumables_[BT].emplace_back(ConsumableType::Product,
153  TypeID{typeid(T)},
154  it.label(),
155  it.instance(),
156  it.process());
157  return ProductToken<T>{it};
158 }
static ProductToken< T > invalid()
Definition: ProductToken.h:47
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
template<typename T , art::BranchType BT>
void art::Consumer::consumesMany ( )
inherited

Definition at line 162 of file Consumer.h.

163 {
164  if (!moduleContext_)
165  return;
166 
167  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
168 }
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
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().

173 {
174  if (!moduleContext_)
175  return ViewToken<T>::invalid();
176 
177  consumables_[BT].emplace_back(ConsumableType::ViewElement,
178  TypeID{typeid(T)},
179  it.label(),
180  it.instance(),
181  it.process());
182  return ViewToken<T>{it};
183 }
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
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().

34 {
35  return rng()->createEngine(
36  placeholder_schedule_id(), seed, kind_of_engine_to_make);
37 }
long seed
Definition: chem4.cc:68
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
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().

43 {
44  return rng()->createEngine(
45  placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
46 }
long seed
Definition: chem4.cc:68
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
CurrentProcessingContext const * art::EDProducer::currentContext ( ) const
protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

121  {
122  return current_context_.get();
123  }
CPC_exempt_ptr current_context_
Definition: EDProducer.h:116
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_().

52 {
53  auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
54  return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
55 }
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_.

124  {
125  return ProducerBase::getProductID<PROD, B>(moduleDescription_,
126  instanceName);
127  }
ModuleDescription moduleDescription_
Definition: EDProducer.h:115
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().

58  {
59  auto const& pd =
60  get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
61  return pd.productID();
62  }
Int_t B
Definition: plot.C:25
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().

191 {
192  if (!moduleContext_)
193  return ProductToken<T>::invalid();
194 
195  consumables_[BT].emplace_back(ConsumableType::Product,
196  TypeID{typeid(T)},
197  it.label(),
198  it.instance(),
199  it.process());
200  return ProductToken<T>{it};
201 }
static ProductToken< T > invalid()
Definition: ProductToken.h:47
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
template<typename T , art::BranchType BT>
void art::Consumer::mayConsumeMany ( )
inherited

Definition at line 205 of file Consumer.h.

206 {
207  if (!moduleContext_)
208  return;
209 
210  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
211 }
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
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().

216 {
217  if (!moduleContext_)
218  return ViewToken<T>::invalid();
219 
220  consumables_[BT].emplace_back(ConsumableType::ViewElement,
221  TypeID{typeid(T)},
222  it.label(),
223  it.instance(),
224  it.process());
225  return ViewToken<T>{it};
226 }
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
bool art::ProducerBase::modifiesEvent ( ) const
inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

41  {
42  return true;
43  }
NucleonDecay& evgen::NucleonDecay::operator= ( NucleonDecay const &  )
delete
NucleonDecay& evgen::NucleonDecay::operator= ( NucleonDecay &&  )
delete
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().

90 {
91  if (!moduleContext_)
92  return;
93 
94  pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
95  for (auto& consumablesPerBranch : consumables_) {
96  cet::sort_all(consumablesPerBranch);
97  }
98 }
bool requireConsumes_
Definition: Consumer.h:137
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
void evgen::NucleonDecay::produce ( art::Event e)
overridevirtual

Implements art::EDProducer.

Definition at line 114 of file NucleonDecay_module.cc.

References simb::MCTruth::Add(), dpdg, art::RandomNumberGenerator::getEngine(), gOptDecayMode, simb::kUnknown, LOG_DEBUG, geo::BoxBoundedGeo::MaxX(), geo::BoxBoundedGeo::MaxY(), geo::BoxBoundedGeo::MaxZ(), mcgen, geo::BoxBoundedGeo::MinX(), geo::BoxBoundedGeo::MinY(), geo::BoxBoundedGeo::MinZ(), geo::GeometryCore::NTPC(), part, art::Event::put(), art::EngineCreator::rng(), simb::MCTruth::SetOrigin(), target, and geo::GeometryCore::TPC().

115 {
116  // Implementation of required member function here.
117  genie::EventRecord * event = new genie::EventRecord;
118  int target = 1000180400; //Only use argon target
119  int decay = (int)gOptDecayMode;
120  genie::Interaction * interaction = genie::Interaction::NDecay(target,decay,dpdg);
121  event->AttachSummary(interaction);
122 
123  // Simulate decay
124  mcgen->ProcessEventRecord(event);
125 
126 // genie::Interaction *inter = event->Summary();
127 // const genie::InitialState &initState = inter->InitState();
128 // std::cout<<"initState = "<<initState.AsString()<<std::endl;
129 // const genie::ProcessInfo &procInfo = inter->ProcInfo();
130 // std::cout<<"procInfo = "<<procInfo.AsString()<<std::endl;
131  LOG_DEBUG("NucleonDecay")
132  << "Generated event: " << *event;
133 
134  std::unique_ptr< std::vector<simb::MCTruth> > truthcol(new std::vector<simb::MCTruth>);
135  simb::MCTruth truth;
136 
139  CLHEP::HepRandomEngine &engine = rng->getEngine();
140  CLHEP::RandFlat flat(engine);
141 
142  // Find boundary of active volume
143  double minx = 1e9;
144  double maxx = -1e9;
145  double miny = 1e9;
146  double maxy = -1e9;
147  double minz = 1e9;
148  double maxz = -1e9;
149  for (size_t i = 0; i<geo->NTPC(); ++i){
150  const geo::TPCGeo &tpc = geo->TPC(i);
151  if (minx>tpc.MinX()) minx = tpc.MinX();
152  if (maxx<tpc.MaxX()) maxx = tpc.MaxX();
153  if (miny>tpc.MinY()) miny = tpc.MinY();
154  if (maxy<tpc.MaxY()) maxy = tpc.MaxY();
155  if (minz>tpc.MinZ()) minz = tpc.MinZ();
156  if (maxz<tpc.MaxZ()) maxz = tpc.MaxZ();
157  }
158 
159  // Assign vertice position
160  double X0 = flat.fire( minx, maxx );
161  double Y0 = flat.fire( miny, maxy );
162  double Z0 = flat.fire( minz, maxz );
163 
164  TIter partitr(event);
165  genie::GHepParticle *part = 0;
166  // GHepParticles return units of GeV/c for p. the V_i are all in fermis
167  // and are relative to the center of the struck nucleus.
168  // add the vertex X/Y/Z to the V_i for status codes 0 and 1
169  int trackid = 0;
170  std::string primary("primary");
171 
172  while( (part = dynamic_cast<genie::GHepParticle *>(partitr.Next())) ){
173 
174  simb::MCParticle tpart(trackid,
175  part->Pdg(),
176  primary,
177  part->FirstMother(),
178  part->Mass(),
179  part->Status());
180 
181  TLorentzVector pos(X0, Y0, Z0, 0);
182  TLorentzVector mom(part->Px(), part->Py(), part->Pz(), part->E());
183  tpart.AddTrajectoryPoint(pos,mom);
184  if(part->PolzIsSet()) {
185  TVector3 polz;
186  part->GetPolarization(polz);
187  tpart.SetPolarization(polz);
188  }
189  tpart.SetRescatter(part->RescatterCode());
190  truth.Add(tpart);
191 
192  ++trackid;
193  }// end loop to convert GHepParticles to MCParticles
194  truth.SetOrigin(simb::kUnknown);
195  truthcol->push_back(truth);
196  //FillHistograms(truth);
197  e.put(std::move(truthcol));
198 
199  delete event;
200  return;
201 
202 }
genie::NucleonDecayMode_t gOptDecayMode
void SetOrigin(simb::Origin_t origin)
Definition: MCTruth.h:78
cout<< "-> Edep in the target
Definition: analysis.C:54
const genie::EventRecordVisitorI * mcgen
double MinX() const
Returns the world x coordinate of the start of the box.
Definition: BoxBoundedGeo.h:90
Geometry information for a single TPC.
Definition: TPCGeo.h:37
double MaxX() const
Returns the world x coordinate of the end of the box.
Definition: BoxBoundedGeo.h:93
void Add(simb::MCParticle &part)
Definition: MCTruth.h:77
ProductID put(std::unique_ptr< PROD > &&product)
Definition: Event.h:102
base_engine_t & getEngine() const
TString part[npart]
Definition: Style.C:32
double MinZ() const
Returns the world z coordinate of the start of the box.
double MaxY() const
Returns the world y coordinate of the end of the box.
unsigned int NTPC(unsigned int cstat=0) const
Returns the total number of TPCs in the specified cryostat.
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
double MaxZ() const
Returns the world z coordinate of the end of the box.
#define LOG_DEBUG(id)
TPCGeo const & TPC(unsigned int const tpc=0, unsigned int const cstat=0) const
Returns the specified TPC.
Event generator information.
Definition: MCTruth.h:30
Namespace collecting geometry-related classes utilities.
double MinY() const
Returns the world y coordinate of the start of the box.
Event finding and building.
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().

126 {
127  if (!moduleContext_)
128  return;
129 
130  // If none of the branches have missing consumes statements, exit early.
131  if (std::all_of(cbegin(missingConsumes_),
132  cend(missingConsumes_),
133  [](auto const& perBranch) { return perBranch.empty(); }))
134  return;
135 
136  constexpr cet::HorizontalRule rule{60};
137  mf::LogPrint log{"MTdiagnostics"};
138  log << '\n'
139  << rule('=') << '\n'
140  << "The following consumes (or mayConsume) statements are missing from\n"
141  << module_context(moduleDescription_) << '\n'
142  << rule('-') << '\n';
143 
144  cet::for_all_with_index(
145  missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
146  for (auto const& pi : perBranch) {
147  log << " "
148  << assemble_consumes_statement(static_cast<BranchType>(i), pi)
149  << '\n';
150  }
151  });
152  log << rule('=');
153 }
cet::exempt_ptr< ModuleDescription const > moduleDescription_
Definition: Consumer.h:140
constexpr T pi()
Returns the constant pi (up to 35 decimal digits of precision)
bool moduleContext_
Definition: Consumer.h:136
ConsumableProductSets missingConsumes_
Definition: Consumer.h:139
void art::Consumer::validateConsumedProduct ( BranchType const  bt,
ProductInfo const &  pi 
)
protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

103 {
104  // Early exits if consumes tracking has been disabled or if the
105  // consumed product is an allowed consumable.
106  if (!moduleContext_)
107  return;
108 
109  if (cet::binary_search_all(consumables_[bt], pi))
110  return;
111 
112  if (requireConsumes_) {
114  "Consumer: an error occurred during validation of a "
115  "retrieved product\n\n")
116  << "The following consumes (or mayConsume) statement is missing from\n"
117  << module_context(moduleDescription_) << ":\n\n"
118  << " " << assemble_consumes_statement(bt, pi) << "\n\n";
119  }
120 
121  missingConsumes_[bt].insert(pi);
122 }
cet::exempt_ptr< ModuleDescription const > moduleDescription_
Definition: Consumer.h:140
bool requireConsumes_
Definition: Consumer.h:137
cet::coded_exception< errors::ErrorCodes, ExceptionDetail::translate > Exception
Definition: Exception.h:66
constexpr T pi()
Returns the constant pi (up to 35 decimal digits of precision)
ConsumableProducts consumables_
Definition: Consumer.h:138
bool moduleContext_
Definition: Consumer.h:136
ConsumableProductSets missingConsumes_
Definition: Consumer.h:139

Member Data Documentation

int evgen::NucleonDecay::dpdg = 0
private

Definition at line 76 of file NucleonDecay_module.cc.

Referenced by NucleonDecay(), and produce().

genie::NucleonDecayMode_t evgen::NucleonDecay::gOptDecayMode = genie::kNDNull
private

Definition at line 75 of file NucleonDecay_module.cc.

Referenced by NucleonDecay(), and produce().

const genie::EventRecordVisitorI* evgen::NucleonDecay::mcgen
private

Definition at line 74 of file NucleonDecay_module.cc.

Referenced by NucleonDecay(), and produce().


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