evgen::NucleonDecay Class Reference

Inheritance diagram for evgen::NucleonDecay:

Public Member Functions
	NucleonDecay (fhicl::ParameterSet const &p)
	NucleonDecay (NucleonDecay const &)=delete
	NucleonDecay (NucleonDecay &&)=delete
NucleonDecay &	operator= (NucleonDecay const &)=delete
NucleonDecay &	operator= (NucleonDecay &&)=delete
void	produce (art::Event &e) override
void	beginRun (art::Run &run) override

Private Attributes
const genie::EventRecordVisitorI *	mcgen
genie::NucleonDecayMode_t	gOptDecayMode = genie::kNDNull
int	dpdg = 0
CLHEP::RandFlat	flatDist

Detailed Description

Definition at line 52 of file NucleonDecay_module.cc.

Constructor & Destructor Documentation

evgen::NucleonDecay::NucleonDecay ( fhicl::ParameterSet const &  p )

explicit

Definition at line 80 of file NucleonDecay_module.cc.

  : art::EDProducer{p}
  // create a default random engine; obtain the random seed from NuRandomService,
  // unless overridden in configuration with key "Seed"
  , flatDist{art::ServiceHandle<rndm::NuRandomService>{}->createEngine(*this, p, "Seed")}

evgen::NucleonDecay::NucleonDecay ( NucleonDecay const &  )

delete

evgen::NucleonDecay::NucleonDecay ( NucleonDecay &&  )

delete

Member Function Documentation

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

override

Definition at line 201 of file NucleonDecay_module.cc.

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

NucleonDecay& evgen::NucleonDecay::operator= ( NucleonDecay const &  )

delete

NucleonDecay& evgen::NucleonDecay::operator= ( NucleonDecay &&  )

delete

void evgen::NucleonDecay::produce ( art::Event &  e )

override

Definition at line 116 of file NucleonDecay_module.cc.

{
  // Implementation of required member function here.
  genie::EventRecord * event = new genie::EventRecord;
  int target = 1000180400;  //Only use argon target
  int decay  = (int)gOptDecayMode;
  genie::Interaction * interaction = genie::Interaction::NDecay(target,decay,dpdg);
  event->AttachSummary(interaction);

  // Simulate decay
  mcgen->ProcessEventRecord(event);

//  genie::Interaction *inter = event->Summary();
//  const genie::InitialState &initState = inter->InitState();
//  std::cout<<"initState = "<<initState.AsString()<<std::endl;
//  const genie::ProcessInfo &procInfo = inter->ProcInfo();
//  std::cout<<"procInfo = "<<procInfo.AsString()<<std::endl;
  MF_LOG_DEBUG("NucleonDecay")
    << "Generated event: " << *event;

  std::unique_ptr< std::vector<simb::MCTruth> > truthcol(new std::vector<simb::MCTruth>);
  simb::MCTruth truth;

  art::ServiceHandle<geo::Geometry const> geo;

  // Find boundary of active volume
  double minx = 1e9;
  double maxx = -1e9;
  double miny = 1e9;
  double maxy = -1e9;
  double minz = 1e9;
  double maxz = -1e9;
  for (size_t i = 0; i<geo->NTPC(); ++i){
    const geo::TPCGeo &tpc = geo->TPC(i);
    if (minx>tpc.MinX()) minx = tpc.MinX();
    if (maxx<tpc.MaxX()) maxx = tpc.MaxX();
    if (miny>tpc.MinY()) miny = tpc.MinY();
    if (maxy<tpc.MaxY()) maxy = tpc.MaxY();
    if (minz>tpc.MinZ()) minz = tpc.MinZ();
    if (maxz<tpc.MaxZ()) maxz = tpc.MaxZ();
  }

  // Assign vertice position
  double X0 = flatDist.fire( minx, maxx );
  double Y0 = flatDist.fire( miny, maxy );
  double Z0 = flatDist.fire( minz, maxz );

  TIter partitr(event);
  genie::GHepParticle *part = 0;
  // GHepParticles return units of GeV/c for p.  the V_i are all in fermis
  // and are relative to the center of the struck nucleus.
  // add the vertex X/Y/Z to the V_i for status codes 0 and 1
  int trackid = 0;
  std::string primary("primary");

  while( (part = dynamic_cast<genie::GHepParticle *>(partitr.Next())) ){

    simb::MCParticle tpart(trackid,
                           part->Pdg(),
                           primary,
                           part->FirstMother(),
                           part->Mass(),
                           part->Status());

    TLorentzVector pos(X0, Y0, Z0, 0);
    TLorentzVector mom(part->Px(), part->Py(), part->Pz(), part->E());
    tpart.AddTrajectoryPoint(pos,mom);
    if(part->PolzIsSet()) {
      TVector3 polz;
      part->GetPolarization(polz);
      tpart.SetPolarization(polz);
    }
    tpart.SetRescatter(part->RescatterCode());
    truth.Add(tpart);

    ++trackid;
  }// end loop to convert GHepParticles to MCParticles
  truth.SetOrigin(simb::kUnknown);
  truthcol->push_back(truth);
  //FillHistograms(truth);
  e.put(std::move(truthcol));

  delete event;
}

Member Data Documentation

int evgen::NucleonDecay::dpdg = 0

private

Definition at line 75 of file NucleonDecay_module.cc.

CLHEP::RandFlat evgen::NucleonDecay::flatDist

private

Definition at line 76 of file NucleonDecay_module.cc.

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

private

Definition at line 74 of file NucleonDecay_module.cc.

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

private

Definition at line 73 of file NucleonDecay_module.cc.

---

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

• NucleonDecay_module.cc