gar::evgen::SingleGen Class Reference

module to produce single or multiple specified particles in the detector More...

Inheritance diagram for gar::evgen::SingleGen:

Public Member Functions
	SingleGen (fhicl::ParameterSet const &pset)
virtual	~SingleGen ()
void	produce (::art::Event &evt)
void	beginRun (::art::Run &run)
void	reconfigure (fhicl::ParameterSet const &p)
Public Member Functions inherited from art::EDProducer
	EDProducer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDProducer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Producer
virtual	~Producer () noexcept
	Producer (fhicl::ParameterSet const &)
	Producer (Producer const &)=delete
	Producer (Producer &&)=delete
Producer &	operator= (Producer const &)=delete
Producer &	operator= (Producer &&)=delete
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)
Public Member Functions inherited from art::Modifier
	~Modifier () noexcept
	Modifier ()
	Modifier (Modifier const &)=delete
	Modifier (Modifier &&)=delete
Modifier &	operator= (Modifier const &)=delete
Modifier &	operator= (Modifier &&)=delete
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
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)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Member Functions
void	SampleOne (unsigned int i, simb::MCTruth &mct)
void	Sample (simb::MCTruth &mct)
void	printVecs (std::vector< std::string > const &list)
bool	PadVector (std::vector< double > &vec)

Private Attributes
int	fMode
bool	fPadOutVectors
std::vector< int >	fPDG
	PDG code of particles to generate. More...
std::vector< double >	fP0
	Central momentum (GeV/c) to generate. More...
std::vector< double >	fSigmaP
	Variation in momenta (GeV/c) More...
int	fPDist
	How to distribute momenta (gaus or uniform) More...
std::vector< double >	fX0
	Central x position (cm) in world coordinates. More...
std::vector< double >	fY0
	Central y position (cm) in world coordinates. More...
std::vector< double >	fZ0
	Central z position (cm) in world coordinates. More...
std::vector< double >	fT0
	Central t position (s) in world coordinates. More...
std::vector< double >	fSigmaX
	Variation in x position (cm) More...
std::vector< double >	fSigmaY
	Variation in y position (cm) More...
std::vector< double >	fSigmaZ
	Variation in z position (cm) More...
std::vector< double >	fSigmaT
	Variation in t position (s) More...
int	fPosDist
	How to distribute xyz (gaus, or uniform) More...
int	fTDist
	How to distribute t (gaus, or uniform) More...
std::vector< double >	fTheta0XZ
	Angle in XZ plane (degrees) More...
std::vector< double >	fTheta0YZ
	Angle in YZ plane (degrees) More...
std::vector< double >	fSigmaThetaXZ
	Variation in angle in XZ plane. More...
std::vector< double >	fSigmaThetaYZ
	Variation in angle in YZ plane. More...
int	fAngleDist
	How to distribute angles (gaus, uniform) More...
CLHEP::HepRandomEngine &	fEngine

Static Private Attributes
static const int	kUNIF = 0
static const int	kGAUS = 1

Additional Inherited Members
Public Types inherited from art::EDProducer
using	ModuleType = EDProducer
using	WorkerType = WorkerT< EDProducer >
Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >
Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >
Static Public Member Functions inherited from art::EDProducer
static void	commitEvent (EventPrincipal &ep, Event &e)
Protected Member Functions inherited from art::ModuleBase
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 ()

Detailed Description

module to produce single or multiple specified particles in the detector

Definition at line 62 of file SingleGen_module.cc.

Constructor & Destructor Documentation

gar::evgen::SingleGen::SingleGen ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 116 of file SingleGen_module.cc.

                                                      :
      art::EDProducer{pset},
      fEngine(art::ServiceHandle<rndm::NuRandomService>()->createEngine(*this,pset,"Seed"))
    {

      this->reconfigure(pset);

      produces< std::vector<simb::MCTruth> >();
      produces< sumdata::RunData, ::art::InRun >();

    }

gar::evgen::SingleGen::~SingleGen ( )

virtual

Definition at line 129 of file SingleGen_module.cc.

    {
    }

Member Function Documentation

void gar::evgen::SingleGen::beginRun

(

::art::Run &

run

)

Definition at line 235 of file SingleGen_module.cc.

    {

      // grab the geometry object to see what geometry we are using
      auto geo = gar::providerFrom<geo::GeometryGAr>();
      std::string name(geo->DetectorName());
      std::unique_ptr<gar::sumdata::RunData> runcol(new gar::sumdata::RunData(name));

      run.put(std::move(runcol));

      return;
    }

bool gar::evgen::SingleGen::PadVector ( std::vector< double > &  vec )

private

Definition at line 208 of file SingleGen_module.cc.

    {
      // check if the vec has the same size as fPDG
      if( vec.size() != fPDG.size() ){
        // if not padding out the vectors always cause an
        // exception to be thrown if the vector in question
        // is not the same size as the fPDG vector
        // the exception is thrown in the reconfigure method
        // that calls this one
        if     (!fPadOutVectors) return false;
        else {
          // if padding of vectors is desired but the vector in
          // question has more than one entry it isn't clear
          // what the padded values should be so cause
          // an exception
          if(vec.size() != 1) return false;

          // pad it out
          vec.resize(fPDG.size(), vec[0]);

        }// end if padding out vectors
      }// end if the vector size is not the same as fPDG

      return true;
    }

void gar::evgen::SingleGen::printVecs ( std::vector< std::string > const &  list )

private

Definition at line 394 of file SingleGen_module.cc.

    {

      MF_LOG_INFO("SingleGen")
      << " You are using vector values for SingleGen configuration.\n   "
      << " Some of the configuration vectors may have been padded out ,"
      << " because they (weren't) as long as the pdg vector"
      << " in your configuration. \n"
      << " The new input particle configuration is:\n" ;

      std::string values;
      for(size_t i = 0; i < list.size(); ++i){

        values.append(list[i]);
        values.append(": [ ");

        for(size_t e = 0; e < fPDG.size(); ++e){
          std::stringstream buf;
          buf.width(10);
          if(i == 0 ) buf << fPDG[e]          << ", ";
          buf.precision(5);
          if(i == 1 ) buf << fP0[e]           << ", ";
          if(i == 2 ) buf << fSigmaP[e]             << ", ";
          if(i == 3 ) buf << fX0[e]                 << ", ";
          if(i == 4 ) buf << fY0[e]                 << ", ";
          if(i == 5 ) buf << fZ0[e]         << ", ";
          if(i == 6 ) buf << fSigmaX[e]             << ", ";
          if(i == 7 ) buf << fSigmaY[e]             << ", ";
          if(i == 8 ) buf << fSigmaZ[e]             << ", ";
          if(i == 9 ) buf << fTheta0XZ[e]     << ", ";
          if(i == 10) buf << fTheta0YZ[e]     << ", ";
          if(i == 11) buf << fSigmaThetaXZ[e] << ", ";
          if(i == 12) buf << fSigmaThetaYZ[e] << ", ";
          if(i == 13) buf << fT0[e]           << ", ";
          if(i == 14) buf << fSigmaT[e]       << ", ";
          values.append(buf.str());
        }

        values.erase(values.find_last_of(","));
        values.append(" ] \n");

      }// end loop over vector names in list

      MF_LOG_INFO("SingleGen") << values;

      return;
    }

void gar::evgen::SingleGen::produce

(

::art::Event &

evt

)

unique_ptr allows ownership to be transferred to the art::Event after the put statement

Definition at line 249 of file SingleGen_module.cc.

    {

      ///unique_ptr allows ownership to be transferred to the ::art::Event after the put statement
      std::unique_ptr< std::vector<simb::MCTruth> > truthcol(new std::vector<simb::MCTruth>);

      simb::MCTruth truth;
      truth.SetOrigin(simb::kSingleParticle);
      Sample(truth);

      MF_LOG_DEBUG("SingleGen") << truth;

      truthcol->push_back(truth);

      evt.put(std::move(truthcol));

      return;
    }

void gar::evgen::SingleGen::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 134 of file SingleGen_module.cc.

    {

      fPadOutVectors = p.get< bool                >("PadOutVectors");
      fMode          = p.get< int                 >("ParticleSelectionMode");
      fPDG           = p.get< std::vector<int>    >("PDG");
      fP0            = p.get< std::vector<double> >("P0");
      fSigmaP        = p.get< std::vector<double> >("SigmaP");
      fPDist         = p.get< int                 >("PDist");
      fTDist         = p.get< int                 >("TDist");
      fX0            = p.get< std::vector<double> >("X0");
      fY0            = p.get< std::vector<double> >("Y0");
      fZ0            = p.get< std::vector<double> >("Z0");
      fT0            = p.get< std::vector<double> >("T0");
      fSigmaX        = p.get< std::vector<double> >("SigmaX");
      fSigmaY        = p.get< std::vector<double> >("SigmaY");
      fSigmaZ        = p.get< std::vector<double> >("SigmaZ");
      fSigmaT        = p.get< std::vector<double> >("SigmaT");
      fPosDist       = p.get< int                 >("PosDist");
      fTheta0XZ      = p.get< std::vector<double> >("Theta0XZ");
      fTheta0YZ      = p.get< std::vector<double> >("Theta0YZ");
      fSigmaThetaXZ  = p.get< std::vector<double> >("SigmaThetaXZ");
      fSigmaThetaYZ  = p.get< std::vector<double> >("SigmaThetaYZ");
      fAngleDist     = p.get< int                 >("AngleDist");

      std::vector<std::string> vlist(15);
      vlist[0]  = "PDG";
      vlist[1]  = "P0";
      vlist[2]  = "SigmaP";
      vlist[3]  = "X0";
      vlist[4]  = "Y0";
      vlist[5]  = "Z0";
      vlist[6]  = "SigmaX";
      vlist[7]  = "SigmaY";
      vlist[8]  = "SigmaZ";
      vlist[9]  = "Theta0XZ";
      vlist[10] = "Theta0YZ";
      vlist[11] = "SigmaThetaXZ";
      vlist[12] = "SigmaThetaYZ";
      vlist[13] = "T0";
      vlist[14] = "SigmaT";


      // begin tests for multiple particle error possibilities
      std::string list;
      if( !this->PadVector(fP0          ) ) list.append(vlist[1] .append(", \n"));
      if( !this->PadVector(fSigmaP      ) ) list.append(vlist[2] .append(", \n"));
      if( !this->PadVector(fX0          ) ) list.append(vlist[3] .append(", \n"));
      if( !this->PadVector(fY0          ) ) list.append(vlist[4] .append(", \n"));
      if( !this->PadVector(fZ0          ) ) list.append(vlist[5] .append(", \n"));
      if( !this->PadVector(fSigmaX      ) ) list.append(vlist[6] .append(", \n"));
      if( !this->PadVector(fSigmaY      ) ) list.append(vlist[7] .append(", \n"));
      if( !this->PadVector(fSigmaZ      ) ) list.append(vlist[8] .append(", \n"));
      if( !this->PadVector(fTheta0XZ    ) ) list.append(vlist[9] .append(", \n"));
      if( !this->PadVector(fTheta0YZ    ) ) list.append(vlist[10].append(", \n"));
      if( !this->PadVector(fSigmaThetaXZ) ) list.append(vlist[11].append(", \n"));
      if( !this->PadVector(fSigmaThetaYZ) ) list.append(vlist[12].append("  \n"));
      if( !this->PadVector(fT0          ) ) list.append(vlist[13].append(", \n"));
      if( !this->PadVector(fSigmaT      ) ) list.append(vlist[14].append(", \n"));

      if(list.size() > 0)
        throw cet::exception("SingleGen") << "The "<< list
        << "\n vector(s) defined in the fhicl files has/have "
        << "a different size than the PDG vector "
        << "\n and it has (they have) more than one value, "
        << "\n disallowing sensible padding "
        << " and/or you have set fPadOutVectors to false. \n";

      if(fPDG.size() > 1 && fPadOutVectors) this->printVecs(vlist);

      return;
    }

void gar::evgen::SingleGen::Sample ( simb::MCTruth &  mct )

private

Definition at line 364 of file SingleGen_module.cc.

    {

      switch (fMode) {
        case 0: // List generation mode: every event will have one of each
                // particle species in the fPDG array
          for (unsigned int i=0; i<fPDG.size(); ++i) {
            SampleOne(i,mct);
          }//end loop over particles
          break;
        case 1: // Random selection mode: every event will exactly one particle
                // selected randomly from the fPDG array
        {
          CLHEP::RandFlat flat(fEngine);

          unsigned int i=flat.fireInt(fPDG.size());
          SampleOne(i,mct);
        }
          break;
        default:
          MF_LOG_WARNING("UnrecognizeOption")
          << "SingleGen does not recognize ParticleSelectionMode "
          << fMode;
          break;
      } // switch on fMode

      return;
    }

void gar::evgen::SingleGen::SampleOne ( unsigned int  i, simb::MCTruth &  mct  )

private

Definition at line 271 of file SingleGen_module.cc.

                                                {

      // get the random number generator service and make some CLHEP generators
      CLHEP::RandFlat   flat(fEngine);
      CLHEP::RandGaussQ gauss(fEngine);

        // Choose momentum
      double p = 0.0;
      double m = 0.0;
      if (fPDist == kGAUS) {
        p = gauss.fire(fP0[i], fSigmaP[i]);
      }
      else {
        p = fP0[i] + fSigmaP[i]*(2.0*flat.fire()-1.0);
      }

      static TDatabasePDG  pdgt;
      TParticlePDG* pdgp = pdgt.GetParticle(fPDG[i]);
      if (pdgp) m = pdgp->Mass();

        // Choose position
      TVector3 x;
      if (fPosDist == kGAUS) {
        x[0] = gauss.fire(fX0[i], fSigmaX[i]);;
        x[1] = gauss.fire(fY0[i], fSigmaY[i]);
        x[2] = gauss.fire(fZ0[i], fSigmaZ[i]);
      }
      else {
        x[0] = fX0[i] + fSigmaX[i]*(2.0*flat.fire()-1.0);
        x[1] = fY0[i] + fSigmaY[i]*(2.0*flat.fire()-1.0);
        x[2] = fZ0[i] + fSigmaZ[i]*(2.0*flat.fire()-1.0);
      }

      double t = 0.;
      if(fTDist==kGAUS){
        t = gauss.fire(fT0[i], fSigmaT[i]);
      }
      else{
        t = fT0[i] + fSigmaT[i]*(2.0*flat.fire()-1.0);
      }

      TLorentzVector pos(x[0], x[1], x[2], t);

      // Choose angles
      double thxz = 0;
      double thyz = 0;

      double thyzrads = 0;
      double thyzradsplussigma = 0;
      double thyzradsminussigma = 0;

      if (fAngleDist == kGAUS) {
        thxz = gauss.fire(fTheta0XZ[i], fSigmaThetaXZ[i]);
        thyz = gauss.fire(fTheta0YZ[i], fSigmaThetaYZ[i]);
      }
      else {

        // Choose angles flat in phase space, which is flat in theta_xz
        // and flat in sin(theta_yz).

        thxz = fTheta0XZ[i] + fSigmaThetaXZ[i]*(2.0*flat.fire()-1.0);

        //Taking asin of sin gives value between -Pi/2 and Pi/2 regardless of user input
        thyzrads           = std::asin(std::sin((M_PI/180.)*(fTheta0YZ[i])));
        thyzradsplussigma  = TMath::Min((thyzrads + ((M_PI/180.)*fabs(fSigmaThetaYZ[i]))), M_PI/2.);
        thyzradsminussigma = TMath::Max((thyzrads - ((M_PI/180.)*fabs(fSigmaThetaYZ[i]))), -M_PI/2.);

        double sinthyzmin = std::sin(thyzradsminussigma);
        double sinthyzmax = std::sin(thyzradsplussigma);
        double sinthyz    = sinthyzmin + flat.fire() * (sinthyzmax - sinthyzmin);
        thyz = (180. / M_PI) * std::asin(sinthyz);
      }

      double thxzrad = thxz*M_PI/180.0;
      double thyzrad = thyz*M_PI/180.0;

      TLorentzVector pvec(p*std::cos(thyzrad)*std::sin(thxzrad),
                          p*std::sin(thyzrad),
                          p*std::cos(thxzrad)*std::cos(thyzrad),
                          std::sqrt(p*p+m*m));

      // set track id to -i as these are all primary particles and have id <= 0
      int trackid = -1 * (i + 1);
      std::string primary("primary");

      simb::MCParticle part(trackid, fPDG[i], primary);
      part.AddTrajectoryPoint(pos, pvec);

      mct.Add(part);
    }

Member Data Documentation

int gar::evgen::SingleGen::fAngleDist

private

How to distribute angles (gaus, uniform)

Definition at line 109 of file SingleGen_module.cc.

CLHEP::HepRandomEngine& gar::evgen::SingleGen::fEngine

private

Definition at line 111 of file SingleGen_module.cc.

int gar::evgen::SingleGen::fMode

private

Particle Selection Mode 0–generate a list of all particles, 1–generate a single particle selected randomly from the list

Definition at line 84 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fP0

private

Central momentum (GeV/c) to generate.

Definition at line 92 of file SingleGen_module.cc.

bool gar::evgen::SingleGen::fPadOutVectors

private

Select to pad out configuration vectors if they are not of of the same length as PDG false: don't pad out - all values need to specified true: pad out - default values assumed and printed out

Definition at line 87 of file SingleGen_module.cc.

std::vector<int> gar::evgen::SingleGen::fPDG

private

PDG code of particles to generate.

Definition at line 91 of file SingleGen_module.cc.

int gar::evgen::SingleGen::fPDist

private

How to distribute momenta (gaus or uniform)

Definition at line 94 of file SingleGen_module.cc.

int gar::evgen::SingleGen::fPosDist

private

How to distribute xyz (gaus, or uniform)

Definition at line 103 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fSigmaP

private

Variation in momenta (GeV/c)

Definition at line 93 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fSigmaT

private

Variation in t position (s)

Definition at line 102 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fSigmaThetaXZ

private

Variation in angle in XZ plane.

Definition at line 107 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fSigmaThetaYZ

private

Variation in angle in YZ plane.

Definition at line 108 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fSigmaX

private

Variation in x position (cm)

Definition at line 99 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fSigmaY

private

Variation in y position (cm)

Definition at line 100 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fSigmaZ

private

Variation in z position (cm)

Definition at line 101 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fT0

private

Central t position (s) in world coordinates.

Definition at line 98 of file SingleGen_module.cc.

int gar::evgen::SingleGen::fTDist

private

How to distribute t (gaus, or uniform)

Definition at line 104 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fTheta0XZ

private

Angle in XZ plane (degrees)

Definition at line 105 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fTheta0YZ

private

Angle in YZ plane (degrees)

Definition at line 106 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fX0

private

Central x position (cm) in world coordinates.

Definition at line 95 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fY0

private

Central y position (cm) in world coordinates.

Definition at line 96 of file SingleGen_module.cc.

std::vector<double> gar::evgen::SingleGen::fZ0

private

Central z position (cm) in world coordinates.

Definition at line 97 of file SingleGen_module.cc.

const int gar::evgen::SingleGen::kGAUS = 1

staticprivate

Definition at line 82 of file SingleGen_module.cc.

const int gar::evgen::SingleGen::kUNIF = 0

staticprivate

Definition at line 81 of file SingleGen_module.cc.

---

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

• garsoft/EventGenerator/SingleGen_module.cc