Inheritance diagram for opdet::CalibrationAnalysis:

Public Member Functions
	CalibrationAnalysis (fhicl::ParameterSet const &)

virtual	~CalibrationAnalysis ()

void	analyze (art::Event const &) 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 Member Functions
void	beginJob () override

void	endJob () override

Private Attributes
std::string	fInputModule

std::string	fInstanceName

std::string	fOpHitModule

double	fSampleFreq

double	fTimeBegin

std::map< int, TH1D * >	averageWaveforms

std::map< int, int >	waveformCount

std::map< int, int >	OpHitCount

std::map< int, double >	FirstHitTimePerChannel

TH1D *	fPedestalMeanPerChannel

TH1D *	fPedestalSigmaPerChannel

TH1D *	fIntegratedSignalMeanPerChannel

TH1D *	fFractionSamplesNearMaximumPerChannel

TH1D *	fNumberOfWaveformsProcessedPerChannel

TH1D *	fFirstOpHitTimeMean

TH1D *	fSecondOpHitTimeMean

TH1D *	fFirstSecondDiffOpHitTimeMean

TH1D *	fNumberOfOpHitsPerChannelPerEvent

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 45 of file CalibrationAnalysis_module.cc.

Constructor & Destructor Documentation

opdet::CalibrationAnalysis::CalibrationAnalysis ( fhicl::ParameterSet const & pset )

Definition at line 108 of file CalibrationAnalysis_module.cc.

         : EDAnalyzer(pset)
     {
 
         // Read the fcl-file
         fInputModule  = pset.get< std::string >("InputModule");
         fInstanceName = pset.get< std::string >("InstanceName");
         fOpHitModule = pset.get<std::string>("OpHitModule");
         // Obtain parameters from TimeService
         auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
         fSampleFreq = clockData.OpticalClock().Frequency();
 
         // Assume starting at 0
         fTimeBegin  = 0;
 
 
     }

opdet::CalibrationAnalysis::~CalibrationAnalysis ( )

virtual

Definition at line 128 of file CalibrationAnalysis_module.cc.

129 {

130 }

Member Function Documentation

void opdet::CalibrationAnalysis::analyze ( art::Event const & evt )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 207 of file CalibrationAnalysis_module.cc.

     {
 
 
         // Get OpDetWaveforms from the event
         art::InputTag itag1(fInputModule, fInstanceName);
         auto waveformHandle = evt.getHandle< std::vector< raw::OpDetWaveform > >(itag1);
 
         // Access ART's TFileService, which will handle creating and writing
         // histograms for us
         art::ServiceHandle< art::TFileService > tfs;
 
         for (size_t i = 0; i < waveformHandle->size(); i++)
         {
 
             // This was probably required to overcome the "const" problem 
             // with OpDetPulse::Waveform()
             art::Ptr< raw::OpDetWaveform > waveformPtr(waveformHandle, i);
             raw::OpDetWaveform pulse = *waveformPtr;
             int channel = pulse.ChannelNumber();
 
             // Create the TH1 if it doesn't exist
             auto waveform = averageWaveforms.find( channel );
             if ( waveform == averageWaveforms.end() ) {
                 TString histName = TString::Format("avgwaveform_channel_%03i", channel);
                 averageWaveforms[channel] =  tfs->make< TH1D >(histName, ";t (us);", pulse.size(), 0, double(pulse.size()) / fSampleFreq);
             }
 
             double PedestalMean = 0;
             double PedestalVariance = 0;
             double IntegratedSignalMean = 0;
             // double IntegratedSignalSigma = 0;
             
 
             // Add this waveform to this histogram
             for (size_t tick = 0; tick < pulse.size(); tick++) {
                 averageWaveforms[channel]->Fill(double(tick)/fSampleFreq, pulse[tick]);
 
                 if(tick < 30)
                   PedestalMean += pulse[tick];
 
                 if(pulse[tick] > 4000)
                   fFractionSamplesNearMaximumPerChannel->Fill(channel,1.0/pulse.size());
             }
 
             PedestalMean /= 30;
 
 
             
             for (size_t tick = 0; tick < pulse.size(); tick++) {
               
               PedestalVariance += ( pulse[tick] - PedestalMean)*( pulse[tick] - PedestalMean);
               
               IntegratedSignalMean += pulse[tick] - PedestalMean;
             }
             
 
             PedestalVariance /= (pulse.size() - 1);
             double PedestalSigma = sqrt(PedestalVariance);
             
             fPedestalMeanPerChannel->Fill(channel,PedestalMean);
             fPedestalSigmaPerChannel->Fill(channel,PedestalSigma);
 
             fIntegratedSignalMeanPerChannel->Fill(channel,IntegratedSignalMean);
 
             // Count number of waveforms on each channel
             waveformCount[channel]++;
 
             
         }
 
         auto OpHitHandle = evt.getHandle< std::vector< recob::OpHit > >(fOpHitModule);
         
         // Map to store how many OpHits are on one optical channel per event
         std::map< int, int   > OpHitCountPerEvent;
 
 
 
 
         if(OpHitHandle->size() > 0)
           std::cout << "OpHitHandle->size() = " << OpHitHandle->size() << std::endl;
 
         for(size_t i=0; i!=OpHitHandle->size(); ++i)
           {
             
             int channel = OpHitHandle->at(i).OpChannel();
             //std::cout << "hit number = " << i << ", online channel number = " << OpHitHandle->at(i).OpChannel() << ", offline channel number = " << channel << std::endl;
             ++OpHitCount[channel];
             ++OpHitCountPerEvent[channel];
 
             double firstophitpeaktime;
             double secondophitpeaktime;
 
             if(OpHitCountPerEvent[channel] == 1){
               firstophitpeaktime = OpHitHandle->at(i).PeakTime();
               FirstHitTimePerChannel[channel] = firstophitpeaktime;
               fFirstOpHitTimeMean->Fill(channel, firstophitpeaktime);
             
             }
             else if(OpHitCountPerEvent[channel] == 2){
               secondophitpeaktime = OpHitHandle->at(i).PeakTime();
               fSecondOpHitTimeMean->Fill(channel, secondophitpeaktime);
               fFirstSecondDiffOpHitTimeMean->Fill(channel, secondophitpeaktime-FirstHitTimePerChannel[channel]);
             }
 
           }
         
         for (auto iter = OpHitCountPerEvent.begin(); iter != OpHitCountPerEvent.end(); iter++){
           fNumberOfOpHitsPerChannelPerEvent->Fill(iter->second);
         }
 
     }

void opdet::CalibrationAnalysis::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 134 of file CalibrationAnalysis_module.cc.

     {
      
       // Access ART's TFileService, which will handle creating and writing
       // histograms for us
       art::ServiceHandle< art::TFileService > tfs;
 
 
       fPedestalMeanPerChannel =  tfs->make< TH1D >("PedestalMeanPerChannel", "Pedestal Means;Channel Number;Pedestal Mean (ADC)", 711, 0.5, 711.5);
       fPedestalSigmaPerChannel =  tfs->make< TH1D >("PedestalSigmaPerChannel", "Pedestal Sigma;Channel Number;Pedestal Sigma (ADC)", 711, 0.5, 711.5);
       fIntegratedSignalMeanPerChannel =  tfs->make< TH1D >("IntegratedSignalMeanPerChannel", "Integrated Signal Means;Channel Number; Integrated SignalMean (ADC)", 711, 0.5, 711.5);
       fIntegratedSignalMeanPerChannel->Sumw2();
       fFractionSamplesNearMaximumPerChannel =  tfs->make< TH1D >("FractionSamplesNearMaximumPerChannel", "Fraction of Samples Near Maximum;Channel Number;Fraction Near Maximum", 711, 0.5, 711.5);
 
       fNumberOfWaveformsProcessedPerChannel =  tfs->make< TH1D >("NumberOfWaveformsProcessedPerChannel", "Number of Waveforms Processed per Channel;Channel Number;Fraction Near Maximum", 711, 0.5, 711.5);
 
       fFirstOpHitTimeMean =  tfs->make< TH1D >("FirstOpHitTimeMean", "Mean of first OpHit time per Channel;Channel Number;First OpHit time mean (ticks)", 711, 0.5, 711.5);
       //fFirstOpHitTimeSigma = tfs->make< TH1D >("FirstOpHitTimeSigma","Sigma of first OpHit time per Channel;Channel Number;First OpHit time sigma (ticks)", 711, 0.5, 711.5);
       fSecondOpHitTimeMean =  tfs->make< TH1D >("SecondOpHitTimeMean", "Mean of second OpHit time per Channel;Channel Number;Second OpHit time mean (ticks)", 711, 0.5, 711.5);
       //fSecondOpHitTimeSigma = tfs->make< TH1D >("SecondOpHitTimeSigma","Sigma of second OpHit time per Channel;Channel Number;Second OpHit time sigma (ticks)", 711, 0.5, 711.5);
       fFirstSecondDiffOpHitTimeMean =  tfs->make< TH1D >("FirstSecondDiffOpHitTimeMean", "Mean of first-second OpHit time difference per Channel;Channel Number;(second - first) OpHit time mean (ticks)", 711, 0.5, 711.5);
       //fFirstSecondDiffOpHitTimeSigma = tfs->make< TH1D >("FirstSecondDiffOpHitTimeSigma","Sigma of first-second OpHit time difference per Channel;Channel Number;(second - first) OpHit time sigma (ticks)", 711, 0.5, 711.5);
      fNumberOfOpHitsPerChannelPerEvent =  tfs->make< TH1D >("NumberOfOpHitsPerChannelPerEvent", "Number of OpHits in one channel per event;Number of OpHits;", 16, -0.5, 15.5);
       
     }

void opdet::CalibrationAnalysis::endJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 161 of file CalibrationAnalysis_module.cc.

     {
 
         for (auto iter = averageWaveforms.begin(); iter != averageWaveforms.end(); iter++)
         {
             mf::LogInfo("Scaling down channel ") << iter->first << " by 1./" << waveformCount[iter->first] << std::endl;
             iter->second->Scale(1./waveformCount[iter->first]);
 
             //Also scale down pedestal mean, sigma, and integrated signal mean and sigma histograms
 
             double adjustedPedestalMean = fPedestalMeanPerChannel->GetBinContent(iter->first)/waveformCount[iter->first];
             fPedestalMeanPerChannel->SetBinContent(iter->first,adjustedPedestalMean);
             double adjustedPedestalSigma = fPedestalSigmaPerChannel->GetBinContent(iter->first)/waveformCount[iter->first];
             fPedestalSigmaPerChannel->SetBinContent(iter->first,adjustedPedestalSigma);
 
             double adjustedFractionSamplesNearMaximumPerChannel = fFractionSamplesNearMaximumPerChannel->GetBinContent(iter->first)/waveformCount[iter->first];
             fFractionSamplesNearMaximumPerChannel->SetBinContent(iter->first,adjustedFractionSamplesNearMaximumPerChannel);
            
             fNumberOfWaveformsProcessedPerChannel->SetBinContent(iter->first,waveformCount[iter->first]);
 
             double adjustedIntegratedSignalMeanPerChannel = fIntegratedSignalMeanPerChannel->GetBinContent(iter->first)/waveformCount[iter->first];
             //double adjustedIntegratedSignalErrorPerChannel = fIntegratedSignalMeanPerChannel->GetBinError(iter->first)/waveformCount[iter->first];
 
             fIntegratedSignalMeanPerChannel->SetBinContent(iter->first,adjustedIntegratedSignalMeanPerChannel);
             //fIntegratedSignalMeanPerChannel->SetBinError(iter->first,adjustedIntegratedSignalErrorPerChannel);
 
         }
 
         for (auto iter = OpHitCount.begin(); iter != OpHitCount.end(); iter++){
           double adjustedFirstOpHitTimeMean = fFirstOpHitTimeMean->GetBinContent(iter->first)/(iter->second);
           fFirstOpHitTimeMean->SetBinContent(iter->first,adjustedFirstOpHitTimeMean);
 
           double adjustedSecondOpHitTimeMean = fSecondOpHitTimeMean->GetBinContent(iter->first)/(iter->second);
           fSecondOpHitTimeMean->SetBinContent(iter->first,adjustedSecondOpHitTimeMean);
 
           double adjustedFirstSecondDiffOpHitTimeMean = fFirstSecondDiffOpHitTimeMean->GetBinContent(iter->first)/(iter->second);
           fFirstSecondDiffOpHitTimeMean->SetBinContent(iter->first,adjustedFirstSecondDiffOpHitTimeMean);
 
         }
 
 
 
     }