Inheritance diagram for opdet::LEDCalibrationAna:

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

void	endJob ()

void	analyze (const art::Event &)

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
uint32_t	ShaperToChannel (uint32_t Shaper)

Private Attributes
std::string	fInputModule

uint32_t	fTriggerChannel

float	fCoincThreshold

float	fMaxTimeMean

float	fMaxTimeThresh

float	fAreaDivs

float	fAreaMin

float	fAreaMax

float	fTriggerDelay

pmtana::PulseRecoManager	fPulseRecoMgr

pmtana::AlgoThreshold	fThreshAlg

pmtana::PedAlgoEdges	fPedAlg

TTree *	fPulseTree

TTree *	fPulseTreeNonCoinc

Float_t	fPeak

Float_t	fPedRMS

Float_t	fPedMean

Float_t	fArea

Float_t	fTBegin

Float_t	fTEnd

Float_t	fTMax

Float_t	fOffset

Int_t	fEventID

Int_t	fRunID

Int_t	fChannel

Int_t	fShaper

bool	fMakeNonCoincTree

std::map< uint32_t, std::vector< double > >	fAreas

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 42 of file LEDCalibrationAna_module.cc.

Constructor & Destructor Documentation

opdet::LEDCalibrationAna::LEDCalibrationAna ( const fhicl::ParameterSet & )

Definition at line 95 of file LEDCalibrationAna_module.cc.

     : EDAnalyzer(pset), fPulseRecoMgr(), fThreshAlg(), fPedAlg()
   {
     // Indicate that the Input Module comes from .fcl
     fInputModule = pset.get<std::string>("InputModule");
     fTriggerChannel = pset.get<uint32_t>("TriggerChannel");
     fTriggerDelay = pset.get<uint32_t>("TriggerDelay");
     fCoincThreshold = pset.get<float>("CoincThreshold");
     fMaxTimeThresh = pset.get<float>("MaxTimeThresh");
     fMaxTimeMean = pset.get<float>("MaxTimeMean");
     fAreaMin = pset.get<float>("AreaMin");
     fAreaMax = pset.get<float>("AreaMax");
     fAreaDivs = pset.get<float>("AreaDivs");
     fMakeNonCoincTree = pset.get<bool>("MakeNonCoincTree");
 
     fPulseRecoMgr.AddRecoAlgo(&fThreshAlg);
     fPulseRecoMgr.SetDefaultPedAlgo(&fPedAlg);
 
     art::ServiceHandle<art::TFileService const> tfs;
 
     fPulseTree = tfs->make<TTree>("fPulseTree", "fPulseTree");
     fPulseTree->Branch("Area", &fArea, "Area/F");
     fPulseTree->Branch("Peak", &fPeak, "Peak/F");
     fPulseTree->Branch("TBegin", &fTBegin, "TBegin/F");
     fPulseTree->Branch("TEnd", &fTEnd, "TEnd/F");
     fPulseTree->Branch("TMax", &fTMax, "TMax/F");
     fPulseTree->Branch("Channel", &fChannel, "Channel/I");
     fPulseTree->Branch("Shaper", &fShaper, "Shaper/I");
     fPulseTree->Branch("PedMean", &fPedMean, "PedMean/F");
     fPulseTree->Branch("PedRMS", &fPedRMS, "PedRMS/F");
     fPulseTree->Branch("Offset", &fOffset, "Offset/F");
     fPulseTree->Branch("EventID", &fEventID, "EventID/I");
     fPulseTree->Branch("RunID", &fRunID, "RunID/I");
 
     fPulseTreeNonCoinc = tfs->make<TTree>("fPulseTreeNonCoinc", "fPulseTreeNonCoinc");
     fPulseTreeNonCoinc->Branch("Area", &fArea, "Area/F");
     fPulseTreeNonCoinc->Branch("Peak", &fPeak, "Peak/F");
     fPulseTreeNonCoinc->Branch("TBegin", &fTBegin, "TBegin/F");
     fPulseTreeNonCoinc->Branch("TEnd", &fTEnd, "TEnd/F");
     fPulseTreeNonCoinc->Branch("TMax", &fTMax, "TMax/F");
     fPulseTreeNonCoinc->Branch("Channel", &fChannel, "Channel/I");
     fPulseTreeNonCoinc->Branch("Shaper", &fShaper, "Shaper/I");
     fPulseTreeNonCoinc->Branch("PedMean", &fPedMean, "PedMean/F");
     fPulseTreeNonCoinc->Branch("PedRMS", &fPedRMS, "PedRMS/F");
     fPulseTreeNonCoinc->Branch("EventID", &fEventID, "EventID/I");
     fPulseTreeNonCoinc->Branch("RunID", &fRunID, "RunID/I");
   }

Member Function Documentation

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

Definition at line 215 of file LEDCalibrationAna_module.cc.

   {
     auto const clock_data =
       art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
 
     fRunID = evt.run();
     fEventID = evt.event();
 
     // Create a handle for our vector of pulses
     art::Handle<std::vector<raw::OpDetWaveform>> OpDetWaveformHandle;
 
     // Read in WaveformHandle
     evt.getByLabel(fInputModule, OpDetWaveformHandle);
 
     std::map<uint32_t, std::vector<int>> OrgOpDigitByChannel;
 
     for (size_t i = 0; i != OpDetWaveformHandle->size(); ++i) {
       OrgOpDigitByChannel[ShaperToChannel(OpDetWaveformHandle->at(i).ChannelNumber())].push_back(i);
     }
 
     std::vector<uint32_t> FrameNumbersForTrig;
     std::vector<uint32_t> TimeSlicesForTrig;
 
     for (size_t i = 0; i != OrgOpDigitByChannel[fTriggerChannel].size(); ++i) {
       double TimeStamp =
         OpDetWaveformHandle->at(OrgOpDigitByChannel[fTriggerChannel][i]).TimeStamp();
       uint32_t Frame = clock_data.OpticalClock().Frame(TimeStamp);
       uint32_t TimeSlice = clock_data.OpticalClock().Sample(TimeStamp);
       FrameNumbersForTrig.push_back(Frame);
       TimeSlicesForTrig.push_back(TimeSlice);
     }
 
     for (size_t i = 0; i != OpDetWaveformHandle->size(); ++i) {
       double TimeStamp = OpDetWaveformHandle->at(i).TimeStamp();
       uint32_t Frame = clock_data.OpticalClock().Frame(TimeStamp);
       uint32_t TimeSlice = clock_data.OpticalClock().Sample(TimeStamp);
       fShaper = OpDetWaveformHandle->at(i).ChannelNumber();
       fChannel = ShaperToChannel(fShaper);
 
       if (uint32_t(fChannel) != fTriggerChannel) {
         for (size_t j = 0; j != FrameNumbersForTrig.size(); ++j) {
           if ((Frame == FrameNumbersForTrig.at(j)) &&
               (fabs(TimeSlice - TimeSlicesForTrig.at(j) - fTriggerDelay) < fCoincThreshold)) {
 
             const raw::OpDetWaveform& wf = OpDetWaveformHandle->at(i);
 
             //fPulseRecoMgr.RecoPulse(wf);
             fPulseRecoMgr.Reconstruct(wf);
 
             size_t NPulses = fThreshAlg.GetNPulse();
 
             fOffset = TimeSlice - TimeSlicesForTrig.at(j);
             //fPedMean = fThreshAlg.PedMean();
             //fPedRMS  = fThreshAlg.PedRms();
 
             for (size_t k = 0; k != NPulses; ++k) {
               if (fabs(fMaxTimeMean - fThreshAlg.GetPulse(k).t_max) < fMaxTimeThresh) {
 
                 fPeak = fThreshAlg.GetPulse(k).peak;
                 fArea = fThreshAlg.GetPulse(k).area;
                 fTBegin = fThreshAlg.GetPulse(k).t_start;
                 fTEnd = fThreshAlg.GetPulse(k).t_end;
                 fTMax = fThreshAlg.GetPulse(k).t_max;
                 fPedMean = fThreshAlg.GetPulse(k).ped_mean;
                 fPedRMS = fThreshAlg.GetPulse(k).ped_sigma;
 
                 fPulseTree->Fill();
 
                 fAreas[fChannel].push_back(fArea);
               }
               else if (fMakeNonCoincTree) {
                 fPeak = fThreshAlg.GetPulse(k).peak;
                 fArea = fThreshAlg.GetPulse(k).area;
                 fTBegin = fThreshAlg.GetPulse(k).t_start;
                 fTEnd = fThreshAlg.GetPulse(k).t_end;
                 fTMax = fThreshAlg.GetPulse(k).t_max;
 
                 fPulseTreeNonCoinc->Fill();
               }
             }
           }
         }
       }
     }
   }

void opdet::LEDCalibrationAna::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 145 of file LEDCalibrationAna_module.cc.

   {
     art::ServiceHandle<art::TFileService const> tfs;
 
     for (auto it = fAreas.begin(); it != fAreas.end(); ++it) {
       uint32_t Channel = it->first;
 
       std::stringstream histname;
       histname.flush();
       histname << "ch" << Channel << "area";
 
       TH1D* HistArea = tfs->make<TH1D>(
         histname.str().c_str(), histname.str().c_str(), fAreaDivs, fAreaMin, fAreaMax);
 
       for (size_t j = 0; j != it->second.size(); ++j) {
         HistArea->Fill(it->second.at(j));
       }
 
       std::stringstream fitname;
       fitname.flush();
       fitname << "ch" << Channel << "fit";
 
       double Max = HistArea->GetMaximum();
       double Mid = HistArea->GetBinContent(fAreaDivs / 2.);
 
       TF1* GausFit = new TF1(fitname.str().c_str(), "gaus(0)+gaus(3)+gaus(6)", fAreaMin, fAreaMax);
 
       GausFit->SetParameters(Mid,
                              (fAreaMin + fAreaMax) / 2.,
                              (fAreaMax - fAreaMin) / 2.,
                              Max,
                              0,
                              (fAreaMin + fAreaMax) / 8.,
                              Max / 5.,
                              0,
                              (fAreaMin + fAreaMax) / 4.);
 
       GausFit->SetParLimits(0, 0, 1.1 * Max);
       GausFit->SetParLimits(1, 0, fAreaMax);
       GausFit->SetParLimits(2, 0, fAreaMax);
 
       GausFit->SetParLimits(3, 0, 1.1 * Max);
       GausFit->FixParameter(4, 0);
       GausFit->SetParLimits(5, 0, (fAreaMin + fAreaMax) / 2.);
 
       GausFit->SetParLimits(6, 0, 1.1 * Max);
       GausFit->FixParameter(7, 0);
       GausFit->SetParLimits(8, 0, (fAreaMin + fAreaMax) / 2.);
 
       HistArea->Fit(GausFit);
 
       double Mean = GausFit->GetParameter(1);
       double Width = GausFit->GetParameter(2);
 
       double MeanErr = GausFit->GetParError(1);
       double WidthErr = GausFit->GetParError(2);
 
       double NPE = pow(Mean, 2) / pow(Width, 2);
       double SPEScale = Mean / NPE;
 
       double NPEError = NPE * pow(2. * (pow(MeanErr / Mean, 2) + pow(WidthErr / Width, 2)), 0.5);
       double SPEError = SPEScale * pow(2. * pow(WidthErr / Width, 2) + pow(MeanErr / Mean, 2), 0.5);
 
       std::cout << "Channel " << Channel << ":\tSPE Scale \t" << SPEScale << "\t +/- \t" << SPEError
                 << ",\t NPE \t" << NPE << "\t +/- \t" << NPEError << std::endl;
     }
   }

uint32_t opdet::LEDCalibrationAna::ShaperToChannel ( uint32_t Shaper )

private

Definition at line 303 of file LEDCalibrationAna_module.cc.

   {
     static std::map<uint32_t, uint32_t> ShaperToChannelMap;
     if (ShaperToChannelMap.size() == 0) {
 
       // temporary
       for (size_t i = 0; i != 40; ++i) {
         ShaperToChannelMap[i] = i;
       }
     }
 
     return ShaperToChannelMap[Shaper];
   }