ex::HelloAuxDet Class Reference

Inheritance diagram for ex::HelloAuxDet:

Public Member Functions
	HelloAuxDet (fhicl::ParameterSet const &p)
	HelloAuxDet (HelloAuxDet const &)=delete
	HelloAuxDet (HelloAuxDet &&)=delete
HelloAuxDet &	operator= (HelloAuxDet const &)=delete
HelloAuxDet &	operator= (HelloAuxDet &&)=delete
void	analyze (art::Event const &e) override
void	beginJob () override
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
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::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
fhicl::ParameterSetID	selectorConfig () const
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 Attributes
std::vector< art::InputTag >	fDigitLabels
std::vector< art::InputTag >	fSimLabels
art::InputTag	fPartLabel
TH1D *	fXDistToTrueTrackFront
TH1D *	fYDistToTrueTrackFront

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
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

Definition at line 99 of file HelloAuxDet_module.cc.

Constructor & Destructor Documentation

ex::HelloAuxDet::HelloAuxDet ( fhicl::ParameterSet const &  p )

explicit

Definition at line 130 of file HelloAuxDet_module.cc.

  :
  EDAnalyzer(p)  // ,
 // More initializers here.
{
  fDigitLabels = p.get<std::vector<art::InputTag>>("DigitLabels", {});
  fSimLabels = p.get<std::vector<art::InputTag>>("SimLabels", {});
  fPartLabel = p.get<art::InputTag>("PartLabel", "largeant");
}

ex::HelloAuxDet::HelloAuxDet ( HelloAuxDet const &  )

delete

ex::HelloAuxDet::HelloAuxDet ( HelloAuxDet &&  )

delete

Member Function Documentation

void ex::HelloAuxDet::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 140 of file HelloAuxDet_module.cc.

{
  //Print() AuxDetDigits
  //::PrintAllIfEmpty<raw::AuxDetDigit>(fDigitLabels, e);

  //Print() AuxDetSimChannels
  //::PrintAllIfEmpty<sim::AuxDetSimChannel>(fSimLabels, e);

  //Produce a map from TrackId to MCParticle
  const auto parts = e.getValidHandle<std::vector<simb::MCParticle>>(fPartLabel);
  std::map<int, simb::MCParticle> idToPart; //Map of Geant TrackId to particle
  for(const auto& part: *parts) idToPart[part.TrackId()] = part;

  idToPart[-1] = idToPart[0]; //TODO: Remove this hack that makes sure the primary particle is recognized
  mf::LogInfo("Number of Particles") << "There are " << parts->size() << " simb::MCParticles in this event.\n";

  //Make histograms showing that channel mapping works for AuxDetSimChannels
  art::ServiceHandle<geo::Geometry> geom;
  for(const auto& label: fSimLabels)
  {
    const auto channels = e.getValidHandle<std::vector<sim::AuxDetSimChannel>>(label); 
    for(const auto& channel: *channels)
    {
      const auto& auxDet = geom->AuxDet(channel.AuxDetID());
      const auto& sens = auxDet.SensitiveVolume(channel.AuxDetSensitiveID());
      const auto center = sens.GetCenter();
      TLorentzVector centerV(center.X(), center.Y(), center.Z(), 0.); //TODO: Use GenVector when I have time to figure it out
      for(const auto& ide: channel.AuxDetIDEs())
      {
        const auto found = idToPart.find(ide.trackID);
        if(found != idToPart.end())
        {
          const auto& part = found->second;
          const auto closestPt = std::min_element(part.Trajectory().begin(), part.Trajectory().end(), 
                                                  [&centerV](const auto& first, const auto& second)
                                                  {
                                                    //SimulationBase coordinates are in mm, but Geometry service coordinates are in cm!  
                                                    //Convert to cm.
                                                    return (first.first*0.1-centerV).Vect().Mag() < (second.first*0.1-centerV).Vect().Mag();
                                                  });

          if(closestPt != part.Trajectory().end())
          {
            //Figure out the x and y positions where the CRT was hit
            //TODO: This is a hack to partially get around my lack of knowledge about the CRT channel mapping.  
            //      I'd like to figure this out by identifying the AuxDet.  
            if(sens.toWorldCoords(geo::AuxDetSensitiveGeo::LocalPoint_t(0, 0, sens.HalfLength())).Y() - center.Y() 
               > sens.toWorldCoords(geo::AuxDetSensitiveGeo::LocalPoint_t(0, 0, sens.HalfLength())).X() - center.X()) 
            //If a CRT hit that gives x position
            {
              fXDistToTrueTrackFront->Fill(centerV.X() - closestPt->first.X());
            }
            else //This is a CRT hit that gives y position
            {
              fYDistToTrueTrackFront->Fill(centerV.Y() - closestPt->first.Y());
            }
          }
          else mf::LogWarning("No Closest Point") << "No point found that is considered \"closest\" to AuxDetSensitive!\n";
        }
        else mf::LogWarning("No Such TrackID") << "No map entry for TrackID " << ide.trackID << "\n";
      }
    }
  }
}

void ex::HelloAuxDet::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 205 of file HelloAuxDet_module.cc.

{
  // Print all of the AuxDetGeos that the Geometry service knows about.  
  const auto geom = lar::providerFrom<geo::Geometry>();

  //Get AuxDetGeos and print their names.  
  for(size_t det = 0; det < geom->NAuxDets(); ++det) 
  {
    const auto& geo = geom->AuxDet(det);
    mf::LogInfo("AuxDetGeometry") << "AuxDetGeo number " << det << ", " << geo.Name() << ", is centered at " << geo.GetCenter() 
                                  << " and has sensitive volumes:\n";
    for(size_t sens = 0; sens < geo.NSensitiveVolume(); ++sens)
    {
      const auto& sensGeo = geo.SensitiveVolume(sens);
      const auto center = sensGeo.GetCenter();
      const bool xPos = (sensGeo.toWorldCoords(geo::AuxDetSensitiveGeo::LocalPoint_t(0, 0, sensGeo.HalfLength())).Y() - center.Y() 
                         > sensGeo.toWorldCoords(geo::AuxDetSensitiveGeo::LocalPoint_t(0, 0, sensGeo.HalfLength())).X() - center.X());
      mf::LogInfo("AuxDetGeometry") << "Sensitive volume " << sens << ": " << sensGeo.TotalVolume()->GetName() << " with center " 
                                    << sensGeo.GetCenter() << " which gives information about " << std::string(xPos?"X":"Y") << " coordinates.\n";
    }
  }

  //Tell the framework that I intend to consume AuxDetDigits and AuxDetSimChannels
  for(const auto& label: fDigitLabels) consumes<std::vector<raw::AuxDetDigit>>(label);
  for(const auto& label: fSimLabels) consumes<std::vector<sim::AuxDetSimChannel>>(label);
  consumes<std::vector<simb::MCParticle>>(fPartLabel);

  //Create histograms owned by "The Framework" that I will fill
  art::ServiceHandle<art::TFileService> tfs;
  fXDistToTrueTrackFront = tfs->make<TH1D>("XDistToTrueTrackFront", "X Distance from True Trajectory to CRT Hits"
                                                                    ";X[mm];Events", 150, -30, 30);
  fYDistToTrueTrackFront = tfs->make<TH1D>("YDistToTrueTrackFront", "Y Distance from True Trajectory to CRT Hits"
                                                                    ";Y[mm];Events", 150, -30, 30);
}

HelloAuxDet& ex::HelloAuxDet::operator= ( HelloAuxDet const &  )

delete

HelloAuxDet& ex::HelloAuxDet::operator= ( HelloAuxDet &&  )

delete

Member Data Documentation

std::vector<art::InputTag> ex::HelloAuxDet::fDigitLabels

private

Definition at line 120 of file HelloAuxDet_module.cc.

art::InputTag ex::HelloAuxDet::fPartLabel

private

Definition at line 122 of file HelloAuxDet_module.cc.

std::vector<art::InputTag> ex::HelloAuxDet::fSimLabels

private

Definition at line 121 of file HelloAuxDet_module.cc.

TH1D* ex::HelloAuxDet::fXDistToTrueTrackFront

private

Definition at line 125 of file HelloAuxDet_module.cc.

TH1D* ex::HelloAuxDet::fYDistToTrueTrackFront

private

Definition at line 126 of file HelloAuxDet_module.cc.

---

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

• duneprototypes/duneprototypes/Protodune/singlephase/CRT/HelloAuxDet_module.cc