gar::rosim::RawDigitAna Class Reference

Inheritance diagram for gar::rosim::RawDigitAna:

Public Member Functions
	RawDigitAna (fhicl::ParameterSet const &pset)
	Standard constructor and destructor for an FMWK module. More...
virtual	~RawDigitAna ()
void	analyze (const ::art::Event &evt)
void	beginJob ()
void	reconfigure (fhicl::ParameterSet const &pset)
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::string	fReadoutModuleLabel
	module label for the Geant More...
TTree *	fChannelTree
	Tree to keep track of sanity check info. More...
EventInfo	fEvt
	Struct containing event identification. More...
EnergyDep	fEDep
	Struct containing energy deposition info. More...
ChannelInfo	fChannelInfo
	Struct containing channel info. More...

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 76 of file RawDigitAna_module.cc.

Constructor & Destructor Documentation

gar::rosim::RawDigitAna::RawDigitAna ( fhicl::ParameterSet const &  pset )

explicit

Standard constructor and destructor for an FMWK module.

Definition at line 102 of file RawDigitAna_module.cc.

    : EDAnalyzer(pset)
    {
      this->reconfigure(pset);
    }

gar::rosim::RawDigitAna::~RawDigitAna ( )

virtual

Definition at line 110 of file RawDigitAna_module.cc.

    {
    }

Member Function Documentation

void gar::rosim::RawDigitAna::analyze

(

const ::art::Event &

evt

)

Definition at line 144 of file RawDigitAna_module.cc.

    {
      // set the event id info
      fEvt.run    = evt.run();
      fEvt.subrun = evt.subRun();
      fEvt.event  = evt.id().event();
      
      // get the list of particles from this event
      auto bt = gar::providerFrom<cheat::BackTracker>();
      
      // get the raw digits for this event
      auto digCol = evt.getValidHandle<std::vector<gar::raw::RawDigit> >(fReadoutModuleLabel);
      
      if( digCol.failedToGet() ){
        MF_LOG_VERBATIM("RawDigitAna")
        << "failed to get raw digits, bail";
        return;
      }
 
      // get the find many for the energy depositions
      ::art::FindMany<gar::sdp::EnergyDeposit> fmEnergyDep(digCol, evt, fReadoutModuleLabel);
      
      if( !fmEnergyDep.isValid() ){
        MF_LOG_WARNING("RawDigitAna")
        << "Unable to find valid association between RawDigits and "
        << "energy deposits, no analysis in this module is possible";
        
        return;
      }

      auto geo = gar::providerFrom<geo::GeometryGAr>();

      float        xyz[3]     = {0.};
      float        xyzChan[3] = {0.};
      unsigned int chan       = 0;
      
      for(size_t d = 0; d < digCol->size(); ++d){
        std::vector<const gar::sdp::EnergyDeposit*> edepsCol;
        fmEnergyDep.get(d, edepsCol);
        if(edepsCol.size() < 1) continue;
        
        MF_LOG_DEBUG("RawDigitAna")
        << "There are "
        << edepsCol.size()
        << " energy depositions for channel "
        << (*digCol)[d].Channel();
        
        geo->ChannelToPosition((*digCol)[d].Channel(), xyz);
        
        // fill the channel information
        fChannelInfo.channel = (int)(*digCol)[d].Channel();
        fChannelInfo.x       = xyz[0];
        fChannelInfo.y       = xyz[1];
        fChannelInfo.z       = xyz[2];
        
        // loop over the energy deposition collections to fill the tree
        for(auto edep : edepsCol){
        
          // get the MCParticle for this track ID
          auto part = bt->TrackIDToParticle(edep->TrackID());
          
          fEDep.trackID = edep->TrackID();
          fEDep.pdg     = part->PdgCode();
          fEDep.x       = edep->X();
          fEDep.y       = edep->Y();
          fEDep.z       = edep->Z();
          fEDep.dX      = edep->dX();
          fEDep.t       = edep->Time();
          fEDep.e       = edep->Energy();
          
          MF_LOG_DEBUG("RawDigitAna")
          << "pos: ("
          << fEDep.x
          << ", "
          << fEDep.y
          << ", "
          << fEDep.z
          << ") e: "
          << fEDep.e
          << " t: "
          << fEDep.t
          << " dX: "
          << fEDep.dX;
          
          // closure test for the channel map
          xyz[0] = fEDep.x;
          xyz[1] = fEDep.y;
          xyz[2] = fEDep.z;
          
          chan = geo->NearestChannel(xyz);
          if (chan != geo->GapChannelNumber())
            {
               geo->ChannelToPosition(chan, xyzChan);
          
//          if(std::abs(xyz[1] - fChannelInfo.y) > 0.33)
//            MF_LOG_VERBATIM("RawDigitAna")
//            << "Channel "
//            << fChannelInfo.channel
//            << " / "
//            << chan
//            << " is off from the energy deposit in y: ("
//            << fChannelInfo.y
//            << ", "
//            << fChannelInfo.z
//            << ") vs ("
//            << xyz[1]
//            << ", "
//            << xyz[2]
//            << ") "
//            << std::abs(xyz[1] - fChannelInfo.y);
//
//          if(std::abs(xyz[2] - fChannelInfo.z) > 10.33)
//            MF_LOG_VERBATIM("RawDigitAna")
//            << "Channel "
//            << fChannelInfo.channel
//            << " / "
//            << chan
//            << " is off from the energy deposit in z: ("
//            << fChannelInfo.y
//            << ", "
//            << fChannelInfo.z
//            << ") vs ("
//            << xyz[1]
//            << ", "
//            << xyz[2]
//            << ") "
//            << std::abs(xyz[2] - fChannelInfo.z);

            }

          // make the tree flat in terms of the energy depositions
          fChannelTree->Fill();
          
        } // end loop over collection of EnergyDeposits
      } // end loop over RawDigit collection
      
      return;
    } // analyze

void gar::rosim::RawDigitAna::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 115 of file RawDigitAna_module.cc.

    {
      ::art::ServiceHandle<::art::TFileService> tfs;
      
      fChannelTree = tfs->make<TTree>("ChannelTree", "ChannelTree");
      
      std::string description("run/I:subrun/I:event/I");
      
      fChannelTree->Branch("info", &fEvt, description.c_str());
      
      description = "trackID/I:pdg/I:x/F:y/F:z/F:e/F:dX/F:t/F";
      
      fChannelTree->Branch("edep", &fEDep, description.c_str());
      
      description = "channel/I:x/F:y/F:z/F";
      
      fChannelTree->Branch("chan", &fChannelInfo, description.c_str());
      
    }

void gar::rosim::RawDigitAna::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 136 of file RawDigitAna_module.cc.

    {
      fReadoutModuleLabel = p.get< std::string >("ReadoutModuleLabel");
      
      return;
    }

Member Data Documentation

ChannelInfo gar::rosim::RawDigitAna::fChannelInfo

private

Struct containing channel info.

Definition at line 93 of file RawDigitAna_module.cc.

TTree* gar::rosim::RawDigitAna::fChannelTree

private

Tree to keep track of sanity check info.

Definition at line 90 of file RawDigitAna_module.cc.

EnergyDep gar::rosim::RawDigitAna::fEDep

private

Struct containing energy deposition info.

Definition at line 92 of file RawDigitAna_module.cc.

EventInfo gar::rosim::RawDigitAna::fEvt

private

Struct containing event identification.

Definition at line 91 of file RawDigitAna_module.cc.

std::string gar::rosim::RawDigitAna::fReadoutModuleLabel

private

module label for the Geant

Definition at line 89 of file RawDigitAna_module.cc.

---

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

• garsoft/ReadoutSimulation/RawDigitAna_module.cc