Inheritance diagram for CRT::CRTSimValidation:

Public Member Functions
	CRTSimValidation (fhicl::ParameterSet const &p)

	CRTSimValidation (CRTSimValidation const &)=delete

	CRTSimValidation (CRTSimValidation &&)=delete

CRTSimValidation &	operator= (CRTSimValidation const &)=delete

CRTSimValidation &	operator= (CRTSimValidation &&)=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
art::InputTag	fCRTLabel

TH1D *	fDetsWithHits

TH1D *	fAllADCs

TH1D *	fDevWithinTrigger

TH1D *	fDeltaTAssoc

TH1D *	fTriggerDeltaT

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 46 of file CRTSimValidation_module.cc.

Constructor & Destructor Documentation

CRT::CRTSimValidation::CRTSimValidation ( fhicl::ParameterSet const & p )

explicit

Definition at line 85 of file CRTSimValidation_module.cc.

   :
   EDAnalyzer(p), fCRTLabel(p.get<art::InputTag>("CRTLabel"))
 {
   consumes<std::vector<CRT::Trigger>>(fCRTLabel);
   consumes<std::vector<art::Assns<sim::AuxDetSimChannel, CRT::Trigger>>>(fCRTLabel);
 }

CRT::CRTSimValidation::CRTSimValidation ( CRTSimValidation const & )

delete

CRT::CRTSimValidation::CRTSimValidation ( CRTSimValidation && )

delete

Member Function Documentation

void CRT::CRTSimValidation::analyze ( art::Event const & e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 93 of file CRTSimValidation_module.cc.

 {
   //Get the data products I plan to work with
   const auto& triggers = e.getValidHandle<std::vector<CRT::Trigger>>(fCRTLabel);
 
   art::FindManyP<sim::AuxDetSimChannel> trigToSim(triggers, e, fCRTLabel);
 
   //Get a handle to the Geometry service to look up AuxDetGeos from module numbers
   art::ServiceHandle<geo::Geometry> geom;
 
   //Mapping from channel to previous Trigger time
   std::unordered_map<size_t, double> prevTimes;
   for(const auto& trigger: *triggers)
   {
     fDetsWithHits->Fill(geom->AuxDet(trigger.Channel()).Name().c_str(), 1.0);
     const auto& hits = trigger.Hits();
     for(const auto& hit: hits) fAllADCs->Fill(hit.ADC());
 
     const auto found = prevTimes.find(trigger.Channel());
     if(found != prevTimes.end())
     {
       fTriggerDeltaT->Fill(trigger.Timestamp() - found->second);
       found->second = trigger.Timestamp();
     } 
     else prevTimes[trigger.Channel()] = trigger.Timestamp();
   }
 
   for(size_t trigIt = 0; trigIt < trigToSim.size(); ++trigIt)
   {
     const auto& trig = (*triggers)[trigIt];
     const auto& simChans = trigToSim.at(trigIt);
 
     size_t nIDEs = 0;
     const auto sumSq = std::accumulate(simChans.begin(), simChans.end(), 0., 
                                        [this, &trig, &nIDEs](double sum, const auto& channel)
                                        {
                                          const auto& ides = channel->AuxDetIDEs();
                                          nIDEs += ides.size();
                                          return sum + std::accumulate(ides.begin(), ides.end(), 0., 
                                                                       [this, &trig](double subSum, const auto& ide)
                                                                       {
                                                                         const auto diff = trig.Timestamp() - 
                                                                                           (ide.exitT + ide.entryT)/2.;
                                                                         fDeltaTAssoc->Fill(diff);
                                                                         if(fabs(diff) > 200) //TODO: Read pset from original module to use readout 
                                                                                              //      time here
                                                                         {
                                                                           mf::LogInfo("LargeDeltaT") << "Found large deltaT: " << diff << "\n"
                                                                                                      << "Timestamp is " << trig.Timestamp() << "\n"
                                                                                                      << "True time is " << (ide.exitT + 
                                                                                                                             ide.entryT)/2. << ".\n";
                                                                         }
 
                                                                         return subSum + diff*diff;
                                                                       });
                                        });
     fDevWithinTrigger->Fill(std::sqrt(sumSq/(nIDEs-1))); //Root of the mean of squares
   }
 }

void CRT::CRTSimValidation::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 153 of file CRTSimValidation_module.cc.

 {
   //Register ROOT objects I want to write with the TFileService
   art::ServiceHandle<art::TFileService> tfs;
   fDetsWithHits     = tfs->make<TH1D>("DetsWithHits", "AuxDetGeo Names that Had CRT::Triggers;AuxDetGeo Name;Events", 1, 0, -1); 
   //Trigger automatic binning since am plotting strings.
   fAllADCs          = tfs->make<TH1D>("AllADCs", "ADC Values from CRT::Hits on All Channels;ADC Value;Hits", 500, 0, 2000);
   fDevWithinTrigger = tfs->make<TH1D>("DevWithinTrigger", "Deviation of True Times from Associated Trigger Timestamp;"
                                                           "True Time #sigma;Energy Deposits", 100, 0, 300);
   fDeltaTAssoc      = tfs->make<TH1D>("DeltaTAssoc", "Time Difference Between Associated IDEs and Trigger Timestamp;#DeltaT_{IDE};Energy Deposits", 
                                       200, -100, 100);
   fTriggerDeltaT    = tfs->make<TH1D>("TriggerDeltaT", "#DeltaT Between Consecutive Triggers in a Module;#DeltaT;Trigger Pairs", 20, 0, 30);
 }

CRTSimValidation& CRT::CRTSimValidation::operator= ( CRTSimValidation const & )

delete

CRTSimValidation& CRT::CRTSimValidation::operator= ( CRTSimValidation && )

delete

Member Data Documentation

TH1D* CRT::CRTSimValidation::fAllADCs

private

Definition at line 73 of file CRTSimValidation_module.cc.

art::InputTag CRT::CRTSimValidation::fCRTLabel

private

Definition at line 67 of file CRTSimValidation_module.cc.

TH1D* CRT::CRTSimValidation::fDeltaTAssoc

private

Definition at line 79 of file CRTSimValidation_module.cc.

TH1D* CRT::CRTSimValidation::fDetsWithHits

private

Definition at line 72 of file CRTSimValidation_module.cc.

TH1D* CRT::CRTSimValidation::fDevWithinTrigger

private

Definition at line 77 of file CRTSimValidation_module.cc.

TH1D* CRT::CRTSimValidation::fTriggerDeltaT

private

Definition at line 80 of file CRTSimValidation_module.cc.

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

duneprototypes/duneprototypes/Protodune/singlephase/CRT/CRTSimValidation_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation