sim::GenericCRTUtility Class Reference

#include "GenericCRT.h"

Public Member Functions
	GenericCRTUtility (const std::string energyUnitsScale)
sim::AuxDetIDE	toAuxDetIDE (sim::AuxDetHit const &InputHit) const
std::vector< unsigned int >	GetAuxDetChannels (const std::vector< sim::AuxDetHit > &InputHitVector) const
std::vector< sim::AuxDetSimChannel >	GetAuxDetSimChannels (const std::vector< sim::AuxDetHit > &InputHitVector) const
sim::AuxDetSimChannel	GetAuxDetSimChannelByNumber (const std::vector< sim::AuxDetHit > &InputHitVector, unsigned int inputchannel) const

Private Attributes
art::ServiceHandle< geo::Geometry const >	fGeo
double	fEnergyUnitsScale

Detailed Description

Definition at line 25 of file GenericCRT.h.

Constructor & Destructor Documentation

sim::GenericCRTUtility::GenericCRTUtility

(

const std::string

energyUnitsScale

)

Title: GenericCRT Utility Class Author: Andrzej Szelc (andrz.nosp@m.ejs@.nosp@m.fnal..nosp@m.gov)

Description: Class with Algorithms to convert sim::AuxDetHits to sim::AuxDetSimChannels

Definition at line 19 of file GenericCRT.cxx.

References fEnergyUnitsScale.

{
  HepTool::Evaluator eval;
  eval.setStdMath();
  eval.setSystemOfUnits();

  const std::string scaleExpression = "MeV / " + energyUnitsScale;
  fEnergyUnitsScale = eval.evaluate(scaleExpression.c_str());

  if (eval.status() != 0) fEnergyUnitsScale = 1.;
}

Member Function Documentation

std::vector< unsigned int > sim::GenericCRTUtility::GetAuxDetChannels

(

const std::vector< sim::AuxDetHit > &

InputHitVector

)

const

Definition at line 52 of file GenericCRT.cxx.

References util::size().

Referenced by GetAuxDetSimChannels().

{

  std::vector<unsigned int> AuxDetChanNumber;
  AuxDetChanNumber.reserve(size(InputHitVector));

  for (auto const& hit : InputHitVector) {

    std::vector<unsigned int>::iterator Chanitr =
      std::find(AuxDetChanNumber.begin(), AuxDetChanNumber.end(), hit.GetID());

    if (Chanitr == AuxDetChanNumber.end()) { //If trackID is already in the map, update it
      //if channel ID is not in the set yet, add it
      AuxDetChanNumber.push_back(hit.GetID());
    } //
  }

  return AuxDetChanNumber;
}

sim::AuxDetSimChannel sim::GenericCRTUtility::GetAuxDetSimChannelByNumber	(	const std::vector< sim::AuxDetHit > &	InputHitVector,
		unsigned int	inputchannel
	)		const

Definition at line 73 of file GenericCRT.cxx.

References fGeo, geo::GeometryCore::FindAuxDetSensitiveAtPosition(), and toAuxDetIDE().

Referenced by GetAuxDetSimChannels().

{
  std::vector<sim::AuxDetIDE> IDEvector;
  //loop over sim::AuxDetHits and assign them to AuxDetSimChannels.

  size_t ad_id_no = 9999;
  size_t ad_sen_id_no = 9999;

  for (auto const& auxDetHit : InputHitVector) {

    double xcoordinate = (auxDetHit.GetEntryX() + auxDetHit.GetExitX()) / 2.0;
    double ycoordinate = (auxDetHit.GetEntryY() + auxDetHit.GetExitY()) / 2.0;
    double zcoordinate = (auxDetHit.GetEntryZ() + auxDetHit.GetExitZ()) / 2.0;
    geo::Point_t const worldPos{xcoordinate, ycoordinate, zcoordinate};

    if (auxDetHit.GetID() == inputchannel) // this is the channel we want.
    {
      // Find the IDs given the hit position
      fGeo->FindAuxDetSensitiveAtPosition(worldPos, ad_id_no, ad_sen_id_no, 0.0001);

      mf::LogDebug("GenericCRTUtility")
        << "Found an AuxDetHit with ID " << auxDetHit.GetID() << " for AuxDet ID " << ad_id_no
        << " Sens ID " << ad_sen_id_no << std::endl;

      auto tempIDE = toAuxDetIDE(auxDetHit);

      std::vector<sim::AuxDetIDE>::iterator IDEitr =
        std::find(IDEvector.begin(), IDEvector.end(), tempIDE);

      if (IDEitr != IDEvector.end()) { //If trackID is already in the map, update it
        //Andrzej's note - following logic from AuxDetReadout in Legacy, but why are the other paremeters getting overwritten like that?
        IDEitr->energyDeposited += tempIDE.energyDeposited;
        IDEitr->exitX = tempIDE.exitX;
        IDEitr->exitY = tempIDE.exitY;
        IDEitr->exitZ = tempIDE.exitZ;
        IDEitr->exitT = tempIDE.exitT;
        IDEitr->exitMomentumX = tempIDE.exitMomentumX;
        IDEitr->exitMomentumY = tempIDE.exitMomentumY;
        IDEitr->exitMomentumZ = tempIDE.exitMomentumZ;
      }
      else { //if trackID is not in the set yet, add it
        IDEvector.push_back(std::move(tempIDE));
      } //else

      // break;
    } // end if the AuxDetHit channel checks out.

  } // end main loop on AuxDetHit

  mf::LogDebug("GenericCRTUtility")
    << "Returning AuxDetSimChannel for ID " << ad_id_no << " " << ad_sen_id_no << ", with "
    << IDEvector.size() << " IDEs." << std::endl;

  //push back the AuxDetSimChannel Vector.
  //TODO check the last parameter values.
  return sim::AuxDetSimChannel(ad_id_no, std::move(IDEvector), ad_sen_id_no);
}

std::vector< sim::AuxDetSimChannel > sim::GenericCRTUtility::GetAuxDetSimChannels

(

const std::vector< sim::AuxDetHit > &

InputHitVector

)

const

Definition at line 133 of file GenericCRT.cxx.

References GetAuxDetChannels(), GetAuxDetSimChannelByNumber(), and util::size().

Referenced by sim::GenericCRT::produce().

{
  auto const auxDetChannels = GetAuxDetChannels(InputHitVector);
  std::vector<sim::AuxDetSimChannel> auxDetVector;
  auxDetVector.reserve(size(auxDetChannels));

  for (auto const channelNum : auxDetChannels) {
    auxDetVector.push_back(GetAuxDetSimChannelByNumber(InputHitVector, channelNum));
  }

  return auxDetVector;
}

sim::AuxDetIDE sim::GenericCRTUtility::toAuxDetIDE

(

sim::AuxDetHit const &

InputHit

)

const

Definition at line 31 of file GenericCRT.cxx.

References sim::AuxDetIDE::energyDeposited, sim::AuxDetIDE::entryT, sim::AuxDetIDE::entryX, sim::AuxDetIDE::entryY, sim::AuxDetIDE::entryZ, sim::AuxDetIDE::exitMomentumX, sim::AuxDetIDE::exitMomentumY, sim::AuxDetIDE::exitMomentumZ, sim::AuxDetIDE::exitT, sim::AuxDetIDE::exitX, sim::AuxDetIDE::exitY, sim::AuxDetIDE::exitZ, fEnergyUnitsScale, sim::AuxDetHit::GetEnergyDeposited(), sim::AuxDetHit::GetEntryT(), sim::AuxDetHit::GetEntryX(), sim::AuxDetHit::GetEntryY(), sim::AuxDetHit::GetEntryZ(), sim::AuxDetHit::GetExitMomentumX(), sim::AuxDetHit::GetExitMomentumY(), sim::AuxDetHit::GetExitMomentumZ(), sim::AuxDetHit::GetExitT(), sim::AuxDetHit::GetExitX(), sim::AuxDetHit::GetExitY(), sim::AuxDetHit::GetExitZ(), sim::AuxDetHit::GetTrackID(), and sim::AuxDetIDE::trackID.

Referenced by GetAuxDetSimChannelByNumber().

{
  sim::AuxDetIDE outputIDE;

  outputIDE.trackID = InputHit.GetTrackID();
  outputIDE.energyDeposited = InputHit.GetEnergyDeposited() * fEnergyUnitsScale;
  outputIDE.entryX = InputHit.GetEntryX();
  outputIDE.entryY = InputHit.GetEntryY();
  outputIDE.entryZ = InputHit.GetEntryZ();
  outputIDE.entryT = InputHit.GetEntryT();
  outputIDE.exitX = InputHit.GetExitX();
  outputIDE.exitY = InputHit.GetExitY();
  outputIDE.exitZ = InputHit.GetExitZ();
  outputIDE.exitT = InputHit.GetExitT();
  outputIDE.exitMomentumX = InputHit.GetExitMomentumX();
  outputIDE.exitMomentumY = InputHit.GetExitMomentumY();
  outputIDE.exitMomentumZ = InputHit.GetExitMomentumZ();

  return outputIDE;
}

Member Data Documentation

double sim::GenericCRTUtility::fEnergyUnitsScale

private

Definition at line 45 of file GenericCRT.h.

Referenced by GenericCRTUtility(), and toAuxDetIDE().

art::ServiceHandle<geo::Geometry const> sim::GenericCRTUtility::fGeo

private

Definition at line 43 of file GenericCRT.h.

Referenced by GetAuxDetSimChannelByNumber().

---

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

• larsim/v09_43_00/source/larsim/GenericCRT/GenericCRT.h
• larsim/v09_43_00/source/larsim/GenericCRT/GenericCRT.cxx