nearline::NearlineAna Class Reference

Inheritance diagram for nearline::NearlineAna:

Public Member Functions
	NearlineAna (fhicl::ParameterSet const &p)
	NearlineAna (NearlineAna const &)=delete
	NearlineAna (NearlineAna &&)=delete
NearlineAna &	operator= (NearlineAna const &)=delete
NearlineAna &	operator= (NearlineAna &&)=delete
void	reconfigure (fhicl::ParameterSet const &p)
void	printConfig ()
void	analyze (art::Event const &e) override
void	beginJob () override
void	endJob () override
size_t	getRawDigits (art::Event const &e, art::Handle< std::vector< raw::RawDigit >> &digitHandle)
size_t	getHits (art::Event const &e, art::Handle< std::vector< recob::Hit >> &hitsHandle)
void	makeHitsPerEventPlots ()
void	fillHitsPerEventPlots (art::Event const &e)
void	makePedestalPerEventPlots ()
void	fillPedestalPerEventPlots (art::Event const &e)
void	writePedestalPerEventSummaryFile (std::string fileName)
void	makePedestalPerTickPlots ()
void	fillPedestalPerTickPlots (art::Event const &e)
void	writePedestalPerTickSummaryFile (std::string fileName)
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 ()
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 &)
std::string const &	processName () const
bool	wantAllEvents () const
bool	wantEvent (Event const &e)
fhicl::ParameterSetID	selectorConfig () const
Handle< TriggerResults >	getTriggerResults (Event const &e) 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
bool	fVerboseOutput
bool	fUseOnlineChannels
std::string	fChannelMapFile
std::map< int, int >	fChannelMap
bool	fMakeHitsPerEventPlots
std::vector< unsigned int >	fHitsPerEventChannels
std::vector< TH1I * >	fVecHitsPerEventPlots
art::InputTag	fHitsTag
bool	fMakePedestalPerEventPlots
bool	fWritePedestalPerEventFile
std::string	fPedestalPerEventFileName
std::vector< unsigned int >	fPedestalPerEventChannels
std::vector< TH1I * >	fVecPedestalPerEventPlots
art::InputTag	fRawDigitsTag
bool	fMakePedestalPerTickPlots
bool	fWritePedestalPerTickFile
std::string	fPedestalPerTickFileName
std::vector< unsigned int >	fPedestalPerTickChannels
std::vector< TH1I * >	fVecPedestalPerTickPlots
TTree *	fHeader
unsigned int	fRun
unsigned int	fSubrun
int	fFirstEvent
int	fLastEvent
int	fNevents
unsigned int	fStartYear
unsigned int	fEndYear
unsigned int	fStartMonth
unsigned int	fEndMonth
unsigned int	fStartDay
unsigned int	fEndDay
double	fStartHour
double	fEndHour
unsigned long long int	fStartTime
unsigned long long int	fEndTime
TH1I *	fHistNearlineVersion

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &paths, fhicl::ParameterSet const &config)
detail::ProcessAndEventSelectors &	processAndEventSelectors ()
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 NearlineAna_module.cc.

Constructor & Destructor Documentation

nearline::NearlineAna::NearlineAna ( fhicl::ParameterSet const &  p )

explicit

Definition at line 127 of file NearlineAna_module.cc.

  :
  EDAnalyzer(p),
  fHeader(0),
  fRun(0),
  fSubrun(0),
  fFirstEvent(1e9),
  fLastEvent(-1),
  fNevents(0),
  fStartYear(0),
  fEndYear(0),
  fStartMonth(0),
  fEndMonth(0),
  fStartDay(0),
  fEndDay(0),
  fStartHour(0.0),
  fEndHour(0.0),
  fStartTime(-1), // this is an unsigned int so it will default to a huge number
  fEndTime(0)
{

  {
    mf::LogInfo logInfo("NearlineAna::NearlineAna");
    logInfo << "NearlineAna" << "\n";
  }
  reconfigure(p);

  printConfig();

}

nearline::NearlineAna::NearlineAna ( NearlineAna const &  )

delete

nearline::NearlineAna::NearlineAna ( NearlineAna &&  )

delete

Member Function Documentation

void nearline::NearlineAna::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 345 of file NearlineAna_module.cc.

{

  //
  // Extract event info for the header...
  //
  unsigned int           run    = e.run();
  unsigned int           subrun = e.subRun();
  unsigned int           event  = e.id().event();
  unsigned long long int time   = e.time().value();

  fNevents++;
  fRun    = run;
  fSubrun = subrun;

  // Don't assume first/last events are coorelated with start/end times...
  if(time < fStartTime && e.time() != art::Timestamp::invalidTimestamp())        fStartTime = time;
  if((int)event < fFirstEvent) fFirstEvent = event;
  if(time > fEndTime)          fEndTime = time;
  if((int)event > fLastEvent)  fLastEvent = event;

  //
  // Fill the desired plots.
  //
  if(fMakeHitsPerEventPlots) fillHitsPerEventPlots(e);
  if(fMakePedestalPerEventPlots) fillPedestalPerEventPlots(e);
  if(fMakePedestalPerTickPlots) fillPedestalPerTickPlots(e);

}

void nearline::NearlineAna::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 247 of file NearlineAna_module.cc.

                                  {

  mf::LogInfo logInfo("NearlineAna::beginJob");
  logInfo << "beginJob" << "\n";

  // book the header object
  art::ServiceHandle<art::TFileService> tfs;
  fHeader = tfs->make<TTree>("Header","Subrun Information");
  
  fHeader->Branch("Run",&fRun);
  fHeader->Branch("Subrun",&fSubrun);
  fHeader->Branch("FirstEvent",&fFirstEvent);
  fHeader->Branch("LastEvent",&fLastEvent);
  fHeader->Branch("Nevents",&fNevents);
  fHeader->Branch("StartYear",&fStartYear);
  fHeader->Branch("StartMonth",&fStartMonth);
  fHeader->Branch("StartDay",&fStartDay);
  fHeader->Branch("StartHour",&fStartHour);
  fHeader->Branch("EndYear",&fEndYear);
  fHeader->Branch("EndMonth",&fEndMonth);
  fHeader->Branch("EndDay",&fEndDay);
  fHeader->Branch("EndHour",&fEndHour);

  // Set Nearline Version Number
  fHistNearlineVersion = tfs->make<TH1I>("hist_nearline_version", "hist_nearline_version", 2, 0, 2);
  fHistNearlineVersion->GetXaxis()->SetBinLabel(1,"NearlineMinorVersion");
  fHistNearlineVersion->GetXaxis()->SetBinLabel(2,"NearlineMajorVersion");
  fHistNearlineVersion->SetBinContent(1, NearlineMinorVersion);
  fHistNearlineVersion->SetBinContent(2, NearlineMajorVersion);
  
  if(fMakeHitsPerEventPlots) makeHitsPerEventPlots();
  if(fMakePedestalPerEventPlots) makePedestalPerEventPlots();
  if(fMakePedestalPerTickPlots) makePedestalPerTickPlots();

}

void nearline::NearlineAna::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 285 of file NearlineAna_module.cc.

                                {

  //
  // Compute header info.
  //

  //
  // DISECTING the time from evt.time().value() into "human readable" format to display the date & time
  //
  unsigned int hour, minute, second;
  int          nano;

  // Get the time stamp.  art::Timestamp::value() returns a TimeValue_t which is a typedef to unsigned long long.
  // The conventional use is for the upper 32 bits to have the seconds since 1970 epoch and the lower 32 bits to be
  // the number of nanoseconds with the current second.
  //
  // NOTE: It seems that the above is NOT the convention for the 35t events. They only seem to use the lower 32 bits
  //       for the year/month/day/hour/second. For now, I have reversed the values of lup and llo to get the time right.
  //
  //       THESE VARIABLES WILL NEED TO BE SWITCHED BACK IF USING THE LOWER 32 BITS FOR NANOSECONDS IS IMPLEMENTED!!!


  const unsigned long int mask32 = 0xFFFFFFFFUL;

  // taking start time apart

  // unsigned long int lup = ( fStartTime >> 32 ) & mask32;
  // unsigned long int llo = fStartTime & mask32;
  unsigned long int llo = ( fStartTime >> 32 ) & mask32; // reversed value (see above comment)
  unsigned long int lup = fStartTime & mask32;           // reversed value (see above comment)
  TTimeStamp ts1(lup, (int)llo);
  ts1.GetDate(kTRUE,0,&fStartYear,&fStartMonth,&fStartDay);
  ts1.GetTime(kTRUE,0,&hour,&minute,&second);
  nano = ts1.GetNanoSec();
  double sec = ((double)second + (double)nano/1.0e9);
  fStartHour = (double)hour + (double)minute/60.0 + sec/3600.0;

  // taking end time apart
  // lup = ( fEndTime >> 32 ) & mask32;
  // llo = fEndTime & mask32;
  llo = ( fEndTime >> 32 ) & mask32; // reversed value (see above comment)
  lup = fEndTime & mask32;           // reversed value (see above comment)
  TTimeStamp ts2(lup, (int)llo);
  ts2.GetDate(kTRUE,0,&fEndYear,&fEndMonth,&fEndDay);
  ts2.GetTime(kTRUE,0,&hour,&minute,&second);
  nano = ts2.GetNanoSec();
  sec = ((double)second + (double)nano/1.0e9);
  fEndHour = (double)hour + (double)minute/60.0 + sec/3600.0;

  fHeader->Fill();

  //Write out pedestal plots summary to files
  if(fWritePedestalPerEventFile) writePedestalPerEventSummaryFile(fPedestalPerEventFileName);
  if(fWritePedestalPerTickFile) writePedestalPerTickSummaryFile(fPedestalPerTickFileName);


}

void nearline::NearlineAna::fillHitsPerEventPlots

(

art::Event const &

e

)

Definition at line 523 of file NearlineAna_module.cc.

                                                                 {


  mf::LogInfo logInfo("NearlineAna::fillHistPerEventPlots");
  art::Handle<std::vector<recob::Hit>> hitHandle;
  size_t numHits = getHits(e, hitHandle);

  //Count the number of hits on channels we are watching
  std::vector<unsigned int> numHitsPerChannel(fHitsPerEventChannels.size(), 0);

  if(numHits==0) return;

  for(size_t hitIter=0;hitIter<numHits;hitIter++){
    art::Ptr<recob::Hit> hitVec(hitHandle,hitIter);
    raw::ChannelID_t channel = hitVec->Channel();
    // geo::View_t      view = hitVec->View();
    // geo::SigType_t   signalType = hitVec->SignalType();
    // geo::WireID      wireID = hitVec->WireID();
    // logInfo << "Hit Channel: " << channel 
    //         << " View: " << view 
    //         << " SignalType: " << signalType 
    //         << " WireID: " << wireID << "\n";

    for(size_t index=0;index<fHitsPerEventChannels.size();index++){
      auto this_channel = (fUseOnlineChannels ? fChannelMap.at(fHitsPerEventChannels.at(index)) : fHitsPerEventChannels.at(index));
      if(this_channel!=channel) continue;
      numHitsPerChannel.at(index) += 1;
    }//channel index
  }//hitIter

  //Now fill histogram
  for(size_t index=0;index<fHitsPerEventChannels.size();index++){
    TH1I* histTemp = fVecHitsPerEventPlots.at(index);
    histTemp->Fill(numHitsPerChannel.at(index));
  }//channel index

  if(fVerboseOutput){
    logInfo << "Channel (Number of Hits): \n";
    for(size_t index=0;index<fHitsPerEventChannels.size();index++){
      logInfo << fHitsPerEventChannels.at(index) << " (" << numHitsPerChannel.at(index) << ") ";
    }
    logInfo << "\n";
  }


}

void nearline::NearlineAna::fillPedestalPerEventPlots

(

art::Event const &

e

)

Definition at line 572 of file NearlineAna_module.cc.

                                                                     {

  mf::LogInfo logInfo("NearlineAna::fillPedestalPerEventPlots");

  // for(size_t index=0;index<fVecPedestalPerEventPlots.size();index++){
  //   auto channel = fPedestalPerEventChannels.at(index);
  //   auto offline_channel = (fUseOnlineChannels ? fChannelMap.at(channel) : channel );
  //   auto hist = fVecPedestalPerEventPlots.at(index);
  //   logInfo << "channel " << (fUseOnlineChannels ? "(online/offline): " + std::to_string(channel) + "/" + std::to_string(offline_channel) : "(offline): " + std::to_string(offline_channel) );
  //   logInfo << " hist_title: " << hist->GetTitle() << "\n";
  // }


  art::Handle<std::vector<raw::RawDigit> > digitHandle;
  size_t numDigitChans = getRawDigits(e, digitHandle);

  //Loop through the vector of rawDigits and pick out channels that match our list of channels

  for(size_t rdIter=0;rdIter<numDigitChans;rdIter++){
    art::Ptr<raw::RawDigit> digitVec(digitHandle, rdIter);
    auto channel =  digitVec->Channel();    
    for(size_t index=0;index<fPedestalPerEventChannels.size();index++){
      auto this_channel = (fUseOnlineChannels ? fChannelMap.at(fPedestalPerEventChannels.at(index)) : fPedestalPerEventChannels.at(index));

      //Only proceed if this is a channel of interest
      if(this_channel!=channel) continue; 

      //DEBUG 
      if(fUseOnlineChannels && fVerboseOutput) logInfo << "this_channel (online/offline): " << this_channel << " (" << fPedestalPerEventChannels.at(index) << "/" << fChannelMap.at(fPedestalPerEventChannels.at(index)) << ")\n";

      auto numSamples = digitVec->Samples();
      auto compression = digitVec->Compression();

      //Only proceed if there is a non-zero number of samples
      if(numSamples==0) continue;

      //We should uncompress the samples in case there is some form of compression applied
      raw::RawDigit::ADCvector_t ADCsUncompressed(numSamples);
      raw::Uncompress(digitVec->ADCs(), ADCsUncompressed, compression);

      //Calculate the channel ADC average
      double averageADC=0;
      for(unsigned int sample=0;sample<numSamples;sample++){
        averageADC+=ADCsUncompressed[sample];
      }//sample

      averageADC/=numSamples; 

      TH1I* histTemp = fVecPedestalPerEventPlots.at(index);
      histTemp->Fill(averageADC);

    }//index
  }//rdIter

}

void nearline::NearlineAna::fillPedestalPerTickPlots

(

art::Event const &

e

)

Definition at line 631 of file NearlineAna_module.cc.

                                                                    {

  mf::LogInfo logInfo("NearlineAna::fillPedestalPerTickPlots");

  // for(size_t index=0;index<fVecPedestalPerTickPlots.size();index++){
  //   auto channel = fPedestalPerTickChannels.at(index);
  //   auto offline_channel = (fUseOnlineChannels ? fChannelMap.at(channel) : channel );
  //   auto hist = fVecPedestalPerTickPlots.at(index);
  //   logInfo << "channel " << (fUseOnlineChannels ? "(online/offline): " + std::to_string(channel) + "/" + std::to_string(offline_channel) : "(offline): " + std::to_string(offline_channel) );
  //   logInfo << " hist_title: " << hist->GetTitle() << "\n";
  // }


  art::Handle<std::vector<raw::RawDigit> > digitHandle;
  size_t numDigitChans = getRawDigits(e, digitHandle);

  //Loop through the vector of rawDigits and pick out channels that match our list of channels

  for(size_t rdIter=0;rdIter<numDigitChans;rdIter++){
    art::Ptr<raw::RawDigit> digitVec(digitHandle, rdIter);
    auto channel =  digitVec->Channel();    
    for(size_t index=0;index<fPedestalPerTickChannels.size();index++){
      auto this_channel = (fUseOnlineChannels ? fChannelMap.at(fPedestalPerTickChannels.at(index)) : fPedestalPerTickChannels.at(index));

      //Only proceed if this is a channel of interest
      if(this_channel!=channel) continue; 

      //DEBUG 
      if(fUseOnlineChannels && fVerboseOutput) logInfo << "this_channel (online/offline): " << this_channel << " (" << fPedestalPerTickChannels.at(index) << "/" << fChannelMap.at(fPedestalPerTickChannels.at(index)) << ")\n";

      auto numSamples = digitVec->Samples();
      auto compression = digitVec->Compression();

      //Only proceed if there is a non-zero number of samples
      if(numSamples==0) continue;

      //We should uncompress the samples in case there is some form of compression applied
      raw::RawDigit::ADCvector_t ADCsUncompressed(numSamples);
      raw::Uncompress(digitVec->ADCs(), ADCsUncompressed, compression);

      TH1I* histTemp = fVecPedestalPerTickPlots.at(index);

      for(unsigned int sample=0;sample<numSamples;sample++){
        histTemp->Fill(ADCsUncompressed[sample]);
      }//sample
    }//index
  }//rdIter

}

size_t nearline::NearlineAna::getHits	(	art::Event const &	e,
		art::Handle< std::vector< recob::Hit >> &	hitsHandle
	)

Definition at line 405 of file NearlineAna_module.cc.

                                                                                                  {

  bool retVal = e.getByLabel(fHitsTag, hitsHandle);
  if(retVal!=true){
    mf::LogWarning("NearlineAna::getHits") << "Getting Hits FAIL: " << fHitsTag << std::endl;
      return 0;
  }
  try { hitsHandle->size(); }
  catch(std::exception const&) {
    mf::LogError("NearlineAna::getHits") << "WARNING: Issue with hitsHandle for Hits" << std::endl;
    return 0;
  }

  if(!hitsHandle.isValid()){
    mf::LogError("NearlineAna::getHits") << "Run: " << e.run()
                            << ", SubRun: " << e.subRun()
                            << ", Event: " << e.event()
                            << " is NOT VALID";
    throw cet::exception("Hit NOT VALID");
    return 0;
  }

  return hitsHandle->size();

}

size_t nearline::NearlineAna::getRawDigits	(	art::Event const &	e,
		art::Handle< std::vector< raw::RawDigit >> &	digitHandle
	)

Definition at line 377 of file NearlineAna_module.cc.

                                                                                                           {

  bool retVal = e.getByLabel(fRawDigitsTag, digitHandle);
  if(retVal!=true){
    mf::LogWarning("NearlineAna::getRawDigits") << "Getting RawDigits FAIL: " << fRawDigitsTag << std::endl;
    return 0;
  }
  
  try { digitHandle->size(); }
  catch(std::exception const&) {
    mf::LogError("NearlineAna::getRawDigits") << "WARNING: Issue with digitHandle for RawDigits" << std::endl;
    return 0;
  }

  if(!digitHandle.isValid()){
    mf::LogError("NearlineAna::getRawDigits") << "Run: " << e.run()
                            << ", SubRun: " << e.subRun()
                            << ", Event: " << e.event()
                            << " is NOT VALID";
    throw cet::exception("RawDigit NOT VALID");
    return 0;
  }

  return digitHandle->size();
}

void nearline::NearlineAna::makeHitsPerEventPlots

(

)

Definition at line 433 of file NearlineAna_module.cc.

                                               {

  mf::LogInfo logInfo("NearlineAna::makeHitsPerEventPlots");
  if(fVerboseOutput) logInfo << "fHitsPerEventChannels:" << "\n";

  art::ServiceHandle<art::TFileService> tfs;

  for(auto channel: fHitsPerEventChannels){
    int numBins = 100;
    int xmin = 0;
    int xmax = 3200;
    std::string hist_name = "hhits_per_event_chan_" + std::to_string(channel);
    std::string hist_title = "Hits Per Event - Channel " 
    + (fUseOnlineChannels ? 
       "(online/offline) " + std::to_string(channel) + "/" + std::to_string(fChannelMap.at(channel)) : 
       "(offline) " + std::to_string(channel));

    TH1I* histTemp = tfs->make<TH1I>(hist_name.c_str(), hist_title.c_str(), numBins, xmin, xmax);
    histTemp->GetXaxis()->SetTitle("Hits per Event");
    histTemp->GetYaxis()->SetTitle("Events");
    fVecHitsPerEventPlots.push_back(histTemp);
    if(fVerboseOutput) logInfo << "channel: " << channel << " hist_name: " << hist_name << " hist_title: " << hist_title << "\n";
  }
  if(fVerboseOutput) logInfo << "\n";

}

void nearline::NearlineAna::makePedestalPerEventPlots

(

)

Definition at line 462 of file NearlineAna_module.cc.

                                                   {


  mf::LogInfo logInfo("NearlineAna::makePedestalPerEventPlots");
  if(fVerboseOutput) logInfo << "fPedestalPerEventChannels:" << "\n";

  art::ServiceHandle<art::TFileService> tfs;

  for(auto channel: fPedestalPerEventChannels){
    int numBins = 128;
    int xmin = 0;
    int xmax = 4096;
    std::string hist_name = "hped_per_event_chan_" + std::to_string(channel);
    std::string hist_title = "Average ADC Per Event - Channel " 
    + (fUseOnlineChannels ? 
       "(online/offline) " + std::to_string(channel) + "/" + std::to_string(fChannelMap.at(channel)) : 
       "(offline) " + std::to_string(channel));

    TH1I* histTemp = tfs->make<TH1I>(hist_name.c_str(), hist_title.c_str(), numBins, xmin, xmax);
    histTemp->GetXaxis()->SetTitle("ADC");
    histTemp->GetYaxis()->SetTitle("Events");
    fVecPedestalPerEventPlots.push_back(histTemp);
    if(fVerboseOutput) logInfo << "channel: " << channel << " hist_name: " << hist_name << " hist_title: " << hist_title << "\n";
  }
  if(fVerboseOutput) logInfo << "\n";

}

void nearline::NearlineAna::makePedestalPerTickPlots

(

)

Definition at line 492 of file NearlineAna_module.cc.

                                                  {


  mf::LogInfo logInfo("NearlineAna::makePedestalPerTickPlots");
  if(fVerboseOutput) logInfo << "fPedestalPerTickChannels:" << "\n";

  art::ServiceHandle<art::TFileService> tfs;

  for(auto channel: fPedestalPerTickChannels){
    int numBins = 128;
    int xmin = 0;
    int xmax = 4096;
    std::string hist_name = "hped_per_tick_chan_" + std::to_string(channel);

    std::string hist_title = "ADC Per Tick - Channel " 
    + (fUseOnlineChannels ? 
       "(online/offline) " + std::to_string(channel) + "/" + std::to_string(fChannelMap.at(channel)) : 
       "(offline) " + std::to_string(channel));

    TH1I* histTemp = tfs->make<TH1I>(hist_name.c_str(), hist_title.c_str(), numBins, xmin, xmax);
    histTemp->GetXaxis()->SetTitle("ADC");
    histTemp->GetYaxis()->SetTitle("Events");
    fVecPedestalPerTickPlots.push_back(histTemp);
    if(fVerboseOutput) logInfo << "channel: " << channel << " hist_name: " << hist_name << " hist_title: " << hist_title << "\n";
  }
  if(fVerboseOutput) logInfo << "\n";

}

NearlineAna& nearline::NearlineAna::operator= ( NearlineAna const &  )

delete

NearlineAna& nearline::NearlineAna::operator= ( NearlineAna &&  )

delete

void nearline::NearlineAna::printConfig

(

)

Definition at line 192 of file NearlineAna_module.cc.

                                     {

  mf::LogInfo logInfo("NearlineAna::printConfig");

  logInfo << "fVerboseOutput: " << (fVerboseOutput ? "true" : "false") << "\n";
  logInfo << "fRawDigitsTag: " << fRawDigitsTag << "\n";
  logInfo << "fHitsTag: " << fHitsTag << "\n";
  logInfo << "fUseOnlineChannels: " << (fUseOnlineChannels ? "true" : "false") << "\n";

  logInfo << "fMakeHitsPerEventPlots: " << (fMakeHitsPerEventPlots ? "true" : "false") << "\n";
  if(fHitsPerEventChannels.size()){
    logInfo << "fHitsPerEventChannels";
    if(fUseOnlineChannels) logInfo << " (online/offline): ";
    else logInfo << "(offline): ";
    for(size_t i=0;i<fHitsPerEventChannels.size();i++){
      logInfo << " " << fHitsPerEventChannels.at(i);
      if(fUseOnlineChannels) logInfo << "/" << fChannelMap.at(fHitsPerEventChannels.at(i));
    }//fHitsPerEventChannels
    logInfo << "\n";
  }

  logInfo << "fMakePedestalPerEventPlots: " << (fMakePedestalPerEventPlots ? "true" : "false") << "\n";
  if(fPedestalPerEventChannels.size()){
    logInfo << "fPedestalPerEventChannels";
    if(fUseOnlineChannels) logInfo << " (online/offline): ";
    else logInfo << "(offline): ";
    for(size_t i=0;i<fPedestalPerEventChannels.size();i++){
      logInfo << " " << fPedestalPerEventChannels.at(i);
      if(fUseOnlineChannels) logInfo << "/" << fChannelMap.at(fPedestalPerEventChannels.at(i));
    }//fPedestalPerEventChannels
    logInfo << "\n";
  }
  logInfo << "fWritePedestalPerEventFile: " << (fWritePedestalPerEventFile ? "true" : "false") << "\n";
  if(fWritePedestalPerEventFile) logInfo << "fPedestalPerEventFileName: " << fPedestalPerEventFileName << "\n";


  logInfo << "fMakePedestalPerTickPlots: " << (fMakePedestalPerTickPlots ? "true" : "false") << "\n";
  if(fPedestalPerTickChannels.size()){
    logInfo << "fPedestalPerTickChannels";
    if(fUseOnlineChannels) logInfo << " (online/offline): ";
    else logInfo << "(offline): ";
    for(size_t i=0;i<fPedestalPerTickChannels.size();i++){
      logInfo << " " << fPedestalPerTickChannels.at(i);
      if(fUseOnlineChannels) logInfo << "/" << fChannelMap.at(fPedestalPerTickChannels.at(i));
    }//fPedestalPerTickChannels
    logInfo << "\n";
  }
  logInfo << "fWritePedestalPerTickFile: " << (fWritePedestalPerTickFile ? "true" : "false") << "\n";
  if(fWritePedestalPerTickFile) logInfo << "fPedestalPerTickFileName: " << fPedestalPerTickFileName << "\n";


}

void nearline::NearlineAna::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 160 of file NearlineAna_module.cc.

                                                                {

  mf::LogInfo logInfo("NearlineAna::reconfigure");
  logInfo << "reconfigure" << "\n";

  fVerboseOutput = p.get<bool>("VerboseOutput", false);

  fUseOnlineChannels = p.get<bool>("UseOnlineChannels", true);
  fChannelMapFile = p.get<std::string>("TPCChannelMapFile");

  if(fUseOnlineChannels) BuildTPCChannelMap(fChannelMapFile, fChannelMap);

  fMakeHitsPerEventPlots = p.get<bool>("MakeHitsPerEventPlots", true);
  fHitsPerEventChannels = p.get<std::vector<unsigned int>>("HitsPerEventChannels", {1,2,3,4});

  fHitsTag = p.get<art::InputTag>("HitsTag", "a:b:c");

  fMakePedestalPerEventPlots = p.get<bool>("MakePedestalPerEventPlots", true);
  fPedestalPerEventChannels = p.get<std::vector<unsigned int>>("PedestalPerEventChannels", {1,2,3,4});
  fWritePedestalPerEventFile = p.get<bool>("WritePedestalPerEventFile", false);
  fPedestalPerEventFileName = p.get<std::string>("PedestalPerEventFileName", "PedestalPerEvent.txt");

  fMakePedestalPerTickPlots = p.get<bool>("MakePedestalPerTickPlots", true);
  fPedestalPerTickChannels = p.get<std::vector<unsigned int>>("PedestalPerTickChannels", {1,2,3,4});
  fWritePedestalPerTickFile = p.get<bool>("WritePedestalPerTickFile", false);
  fPedestalPerTickFileName = p.get<std::string>("PedestalPerTickFileName", "PedestalPerTick.txt");

  fRawDigitsTag = p.get<art::InputTag>("RawDigitsTag", "a:b:c");
}

void nearline::NearlineAna::writePedestalPerEventSummaryFile

(

std::string

fileName

)

Definition at line 683 of file NearlineAna_module.cc.

                                                                            {
  std::ostringstream my_ostream;
  my_ostream << "online_channel " << "offline_channel " << "pedestal_mean " << "pedestal_rms " << "\n";

  for(size_t index=0;index<fPedestalPerEventChannels.size();index++){
    auto online_channel = -1;
    auto offline_channel = fPedestalPerEventChannels.at(index);
    if(fUseOnlineChannels){
      online_channel = fPedestalPerEventChannels.at(index);
      offline_channel = fChannelMap.at(fPedestalPerEventChannels.at(index));
    } 
    TH1I* histTemp = fVecPedestalPerEventPlots.at(index);
    double mean = histTemp->GetMean();
    double rms = histTemp->GetRMS();
    my_ostream << online_channel << " " << offline_channel << " " << mean << " " << rms << "\n";
  }//index

  std::ofstream outFile(fileName);
  if(outFile.is_open()) outFile << my_ostream.str();
  else mf::LogWarning("writePedestalPerEventSummaryFile") << "FAILED to open file: " << fileName;
  outFile.close();

  //DEBUG
  mf::LogInfo("writePedestalPerEventSummaryFile") << my_ostream.str();
}

void nearline::NearlineAna::writePedestalPerTickSummaryFile

(

std::string

fileName

)

Definition at line 711 of file NearlineAna_module.cc.

                                                                           {
  std::ostringstream my_ostream;
  my_ostream << "online_channel " << "offline_channel " << "pedestal_mean " << "pedestal_rms " << "\n";

  for(size_t index=0;index<fPedestalPerTickChannels.size();index++){
    auto online_channel = -1;
    auto offline_channel = fPedestalPerTickChannels.at(index);
    if(fUseOnlineChannels){
      online_channel = fPedestalPerTickChannels.at(index);
      offline_channel = fChannelMap.at(fPedestalPerTickChannels.at(index));
    } 
    TH1I* histTemp = fVecPedestalPerTickPlots.at(index);
    double mean = histTemp->GetMean();
    double rms = histTemp->GetRMS();
    my_ostream << online_channel << " " << offline_channel << " " << mean << " " << rms << "\n";
  }//index

  std::ofstream outFile(fileName);
  if(outFile.is_open()) outFile << my_ostream.str();
  else mf::LogWarning("writePedestalPerEventSummaryFile") << "FAILED to open file: " << fileName;
  outFile.close();

  //DEBUG
  mf::LogInfo("writePedestalPerTickSummaryFile") << my_ostream.str();

}

Member Data Documentation

std::map<int,int> nearline::NearlineAna::fChannelMap

private

Definition at line 81 of file NearlineAna_module.cc.

std::string nearline::NearlineAna::fChannelMapFile

private

Definition at line 80 of file NearlineAna_module.cc.

unsigned int nearline::NearlineAna::fEndDay

private

Definition at line 114 of file NearlineAna_module.cc.

double nearline::NearlineAna::fEndHour

private

Definition at line 116 of file NearlineAna_module.cc.

unsigned int nearline::NearlineAna::fEndMonth

private

Definition at line 112 of file NearlineAna_module.cc.

unsigned long long int nearline::NearlineAna::fEndTime

private

Definition at line 118 of file NearlineAna_module.cc.

unsigned int nearline::NearlineAna::fEndYear

private

Definition at line 110 of file NearlineAna_module.cc.

int nearline::NearlineAna::fFirstEvent

private

Definition at line 106 of file NearlineAna_module.cc.

TTree* nearline::NearlineAna::fHeader

private

Definition at line 103 of file NearlineAna_module.cc.

TH1I* nearline::NearlineAna::fHistNearlineVersion

private

Definition at line 121 of file NearlineAna_module.cc.

std::vector<unsigned int> nearline::NearlineAna::fHitsPerEventChannels

private

Definition at line 84 of file NearlineAna_module.cc.

art::InputTag nearline::NearlineAna::fHitsTag

private

Definition at line 86 of file NearlineAna_module.cc.

int nearline::NearlineAna::fLastEvent

private

Definition at line 107 of file NearlineAna_module.cc.

bool nearline::NearlineAna::fMakeHitsPerEventPlots

private

Definition at line 83 of file NearlineAna_module.cc.

bool nearline::NearlineAna::fMakePedestalPerEventPlots

private

Definition at line 88 of file NearlineAna_module.cc.

bool nearline::NearlineAna::fMakePedestalPerTickPlots

private

Definition at line 95 of file NearlineAna_module.cc.

int nearline::NearlineAna::fNevents

private

Definition at line 108 of file NearlineAna_module.cc.

std::vector<unsigned int> nearline::NearlineAna::fPedestalPerEventChannels

private

Definition at line 91 of file NearlineAna_module.cc.

std::string nearline::NearlineAna::fPedestalPerEventFileName

private

Definition at line 90 of file NearlineAna_module.cc.

std::vector<unsigned int> nearline::NearlineAna::fPedestalPerTickChannels

private

Definition at line 98 of file NearlineAna_module.cc.

std::string nearline::NearlineAna::fPedestalPerTickFileName

private

Definition at line 97 of file NearlineAna_module.cc.

art::InputTag nearline::NearlineAna::fRawDigitsTag

private

Definition at line 93 of file NearlineAna_module.cc.

unsigned int nearline::NearlineAna::fRun

private

Definition at line 104 of file NearlineAna_module.cc.

unsigned int nearline::NearlineAna::fStartDay

private

Definition at line 113 of file NearlineAna_module.cc.

double nearline::NearlineAna::fStartHour

private

Definition at line 115 of file NearlineAna_module.cc.

unsigned int nearline::NearlineAna::fStartMonth

private

Definition at line 111 of file NearlineAna_module.cc.

unsigned long long int nearline::NearlineAna::fStartTime

private

Definition at line 117 of file NearlineAna_module.cc.

unsigned int nearline::NearlineAna::fStartYear

private

Definition at line 109 of file NearlineAna_module.cc.

unsigned int nearline::NearlineAna::fSubrun

private

Definition at line 105 of file NearlineAna_module.cc.

bool nearline::NearlineAna::fUseOnlineChannels

private

Definition at line 79 of file NearlineAna_module.cc.

std::vector<TH1I*> nearline::NearlineAna::fVecHitsPerEventPlots

private

Definition at line 85 of file NearlineAna_module.cc.

std::vector<TH1I*> nearline::NearlineAna::fVecPedestalPerEventPlots

private

Definition at line 92 of file NearlineAna_module.cc.

std::vector<TH1I*> nearline::NearlineAna::fVecPedestalPerTickPlots

private

Definition at line 99 of file NearlineAna_module.cc.

bool nearline::NearlineAna::fVerboseOutput

private

Definition at line 77 of file NearlineAna_module.cc.

bool nearline::NearlineAna::fWritePedestalPerEventFile

private

Definition at line 89 of file NearlineAna_module.cc.

bool nearline::NearlineAna::fWritePedestalPerTickFile

private

Definition at line 96 of file NearlineAna_module.cc.

---

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

• dunetpc/dune/NearlineMonitor/NearlineAna_module.cc