Inheritance diagram for DAQToOffline::GoodRun:

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

virtual	~GoodRun ()

void	analyze (art::Event const &evt) override

void	beginJob () override

void	beginRun (const art::Run &r) override

void	endRun (const art::Run &r) override

void	endJob () override

void	printParameterSet ()

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
TTree *	fTree

int	RunNumber

int	TotEvents

int	NumberOfRCEs

double	FracRCEs

double	AvADCPedDiff

double	FracWaveforms

double	FracOpHits

int	nPTBTrigsOn110

int	nPTBTrigsOn111

int	nPTBTrigsOn112

int	nPTBTrigsOn113

int	nPTBTrigsOn114

int	nPTBTrigsOn115

int	EvJustSSPs

TTree *	FlatTree

int	Event

unsigned int	DigSize

unsigned int	NSamples

std::vector< int >	ChanNum

std::vector< float >	ADCVec

bool	fMakeFlatTree

std::string	fCounterModuleLabel

std::string	fWaveformModuleLabel

std::string	fRawDigitModuleLabel

std::string	fOpHitModuleLabel

short unsigned int	fMaxSamples

std::vector< int >	MyUsefulChans

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 47 of file GoodRun_module.cc.

Constructor & Destructor Documentation

DAQToOffline::GoodRun::GoodRun ( fhicl::ParameterSet const & p )

explicit

Definition at line 96 of file GoodRun_module.cc.

   : EDAnalyzer(pset)
   , fMakeFlatTree       (pset.get<bool>("MakeFlatTree"))
     //---------------------------RCE--------------------------------------
   , fCounterModuleLabel (pset.get<std::string>("CounterModuleLabel"))
   , fWaveformModuleLabel(pset.get<std::string>("WaveformModuleLabel"))
   , fRawDigitModuleLabel(pset.get<std::string>("RawDigitModuleLabel"))
   , fOpHitModuleLabel   (pset.get<std::string>("OpHitModuleLabel"))
   , fMaxSamples         (pset.get<short unsigned int>("MaxSamples", 2048))
     //---------------------------SSP--------------------------------------
 //  ,fSSPFragType        ( pset.get<std::string>("SSPFragType"))
 //  ,fSSPRawDataLabel    ( pset.get<std::string>("SSPRawDataLabel"))
     //---------------------------PTB--------------------------------------
 //  ,fPTBFragType        ( pset.get<std::string>("PTBFragType"))
 //  ,fPTBRawDataLabel    ( pset.get<std::string>("PTBRawDataLabel"))
     //--------------------------------------------------------------------
 {
 }

DAQToOffline::GoodRun::~GoodRun ( )

virtual

Definition at line 115 of file GoodRun_module.cc.

115 {

116 }

Member Function Documentation

void DAQToOffline::GoodRun::analyze ( art::Event const & evt )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 166 of file GoodRun_module.cc.

 {
   if (TotEvents==0) {
     RunNumber = evt.run();
     std::cout << "Got runNumber it is " << RunNumber << std::endl;
   }
   ++TotEvents;
 
   // get raw::RawDigits
   art::Handle< std::vector< raw::RawDigit> > externalRawDigitListHandle;
   std::vector< art::Ptr< raw::RawDigit> > digits;
   if (evt.getByLabel(fRawDigitModuleLabel, externalRawDigitListHandle) )
     art::fill_ptr_vector(digits,externalRawDigitListHandle);
   
   if (fMakeFlatTree) {
     Event    = evt.event();
     DigSize  = digits.size();
     NSamples = digits[0]->Samples();
     if (NSamples > fMaxSamples) NSamples = fMaxSamples;
     // skip this because types may be different: NSamples = std::min( digits[0]->Samples(), fMaxSamples );
 
     ADCVec.clear();
     for (unsigned int dig=0; dig<DigSize; ++dig ) {
       ChanNum.push_back( digits[dig]->Channel() );
       unsigned int Samp=0;
       while ( Samp < NSamples ) {
         //std::cerr << "Looking at digit " << dig << ", sample " << Samp << ", my value is " << digits[dig]->ADC(Samp) << std::endl;
         ADCVec.push_back ( digits[dig]->ADC(Samp) );
         ++Samp;
       }
       //std::cout << "Looking at Event " << Event << ", DigSize " << DigSize << ", ChanNum " << ChanNum[dig] << ", NSamples " << NSamples << ", ADCVec.size() " << ADCVec.size() << std::endl;
     }
     FlatTree->Fill();
   } else {
     // get raw::ExternalTriggers
     art::Handle< std::vector< raw::ExternalTrigger> > externalTriggerListHandle;
     std::vector< art::Ptr< raw::ExternalTrigger> > trigs;
     if (evt.getByLabel(fCounterModuleLabel, externalTriggerListHandle) )
       art::fill_ptr_vector(trigs,externalTriggerListHandle);
 
     // get raw::OpDetWaveforms
     art::Handle< std::vector< raw::OpDetWaveform> > externalWaveformListHandle;
     std::vector< art::Ptr< raw::OpDetWaveform> > waveforms;
     if (evt.getByLabel(fWaveformModuleLabel, externalWaveformListHandle) )
       art::fill_ptr_vector(waveforms,externalWaveformListHandle);
 
     // get recob::OpHits
     art::Handle< std::vector< recob::OpHit> > externalOpHitListHandle;
     std::vector< art::Ptr< recob::OpHit> > ophits;
     if (evt.getByLabel(fOpHitModuleLabel, externalOpHitListHandle) )
       art::fill_ptr_vector(ophits,externalOpHitListHandle);
 
     bool RCEsPresent = false;
     //bool PTBPresent  = false;
     bool WavePresent = false;
     //bool OHitPresent = false;
 
     // RCEs
     if (digits.size()) {
       RCEsPresent = true;
       ++FracRCEs;
       if (NumberOfRCEs == 0) {
         NumberOfRCEs = digits.size() / 128;
       
         //GET THE LIST OF BAD CHANNELS.
         lariov::ChannelStatusProvider const& channelStatus = art::ServiceHandle<lariov::ChannelStatusService>()->GetProvider();
         lariov::ChannelStatusProvider::ChannelSet_t const BadChannels = channelStatus.BadChannels();
       
         for (size_t dig=0; dig<digits.size(); ++dig) {
           double AvADC = 0;
           unsigned int Channel = digits[dig]->Channel();
           double Pedestal = digits[dig]->GetPedestal();
           for (size_t samp=0; samp<digits[dig]->Samples(); ++samp)
             AvADC += digits[dig]->ADC(samp);
           AvADC = AvADC / digits[dig]->Samples();
           // Do I want to use this channel? Both AvADC and Pedestal, plus not a 'bad' channel...
           bool UseChan = true;
           if ( AvADC == 0 && Pedestal == 0)
             UseChan = false;
           for (auto it = BadChannels.begin(); it != BadChannels.end(); it++) {
             if(Channel==*it) {
               UseChan = false;
               break;
             }
           } // Loop through bad chans.
           if ( UseChan ) {
             double ADCPedDiff = AvADC - Pedestal;
             //ADCDiff->Fill( Channel, ADCPedDiff );
             //ADCPed ->Fill( AvADC  , Pedestal   );
             AvADCPedDiff += fabs(ADCPedDiff);
           } // If have data for this tick, ie not turned off.
         } // Loop over digits
         AvADCPedDiff = AvADCPedDiff / digits.size();
       } // NumberOfRCEs == 0
     }
   
     // PTB
     if (trigs.size()) {
       //PTBPresent = true;
       for (size_t tr=0; tr<trigs.size(); ++tr) {
         if (trigs[tr]->GetTrigID() < 100) continue;
         if (trigs[tr]->GetTrigID() == 110) ++nPTBTrigsOn110;
         if (trigs[tr]->GetTrigID() == 111) ++nPTBTrigsOn111;
         if (trigs[tr]->GetTrigID() == 112) ++nPTBTrigsOn112;
         if (trigs[tr]->GetTrigID() == 113) ++nPTBTrigsOn113;
         if (trigs[tr]->GetTrigID() == 114) ++nPTBTrigsOn114;
         if (trigs[tr]->GetTrigID() == 115) ++nPTBTrigsOn115;
       }
     }
 
     // Waveforms
     if (waveforms.size()) {
       WavePresent = true;
       ++FracWaveforms;
       if (!RCEsPresent) ++EvJustSSPs;
     }
 
     // OpHits
     if (ophits.size()) {
       //OHitPresent = true;
       ++FracOpHits;
       if (!RCEsPresent && !WavePresent) ++EvJustSSPs;
     }
   }
 }

void DAQToOffline::GoodRun::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 118 of file GoodRun_module.cc.

                                    {
   
   art::ServiceHandle<art::TFileService> tfs;
     
   if (!fMakeFlatTree) {
     fTree = tfs->make<TTree>("RunList","RunList Information");
     fTree->Branch("RunNumber"     ,&RunNumber     );
     fTree->Branch("TotEvents"     ,&TotEvents     );
     fTree->Branch("NumberOfRCEs"  ,&NumberOfRCEs  );
     fTree->Branch("FracRCEs"      ,&FracRCEs      );
     fTree->Branch("AvADCPedDiff"  ,&AvADCPedDiff  );
     fTree->Branch("FracWaveforms" ,&FracWaveforms );
     fTree->Branch("FracOpHits"    ,&FracOpHits    );
     fTree->Branch("nPTBTrigsOn110",&nPTBTrigsOn110);
     fTree->Branch("nPTBTrigsOn111",&nPTBTrigsOn111);
     fTree->Branch("nPTBTrigsOn112",&nPTBTrigsOn112);
     fTree->Branch("nPTBTrigsOn113",&nPTBTrigsOn113);
     fTree->Branch("nPTBTrigsOn114",&nPTBTrigsOn114);
     fTree->Branch("nPTBTrigsOn115",&nPTBTrigsOn115);
     fTree->Branch("EvJustSSPs"    ,&EvJustSSPs    );
     //ADCDiff = tfs->make<TH2D>("ADCDiff","Difference in average ADC value and pedestal; Channel number; Difference (ADCs)", 2049,0,2048, 400,-50,50);
     //ADCPed  = tfs->make<TH2D>("ADCPed" ,"Comparison of Average ADC values, and pedestal; Average ADC; Pedestal", 1200, 400, 1000, 1200, 400, 1000);
 
   } else {
 
     // I need to choose which channels I want to look at.....
     // Runs 18407 - 18411 used RCEs 0, 4,     12,     15
     // Runs 18412 - 18422 used RCEs 0, 4, 11, 12, 13, 15
     // Runs 18423 - 18430 used RCEs 0, 4, 11
 
     // To convert from the RCE number to the channels I want to store, look at Mikes detailed channel map here:
     //   https://cdcvs.fnal.gov/redmine/projects/35ton/wiki/Channel_map
     FlatTree = tfs->make<TTree>("FlatDigitTree","FlatDigitTree");
     FlatTree->Branch("Event"   ,&Event   ,"Event/I"   );
     FlatTree->Branch("NSamples",&NSamples,"NSamples/I");
     FlatTree->Branch("DigSize" ,&DigSize ,"DigSize/I" );
     FlatTree->Branch("ChanNum" ,&ChanNum);
     FlatTree->Branch("ADCVec"  ,&ADCVec );
   }
 }

void DAQToOffline::GoodRun::beginRun ( const art::Run & r )

override

Definition at line 159 of file GoodRun_module.cc.

                                                    {
 
   RunNumber = TotEvents = NumberOfRCEs = FracRCEs = AvADCPedDiff = FracWaveforms = FracOpHits = nPTBTrigsOn110 = nPTBTrigsOn111 = nPTBTrigsOn112 = nPTBTrigsOn113 = nPTBTrigsOn114 = nPTBTrigsOn115 = EvJustSSPs = 0;
   std::cout << "At the start of the run....going to reset all the variables. " << std::endl;
 
 }

void DAQToOffline::GoodRun::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 306 of file GoodRun_module.cc.

306 {

307 }

void DAQToOffline::GoodRun::endRun ( const art::Run & r )

override

Definition at line 292 of file GoodRun_module.cc.

 {
   if (!fMakeFlatTree) {
     FracRCEs      = FracRCEs / TotEvents;
     FracWaveforms = FracWaveforms / TotEvents;
     FracOpHits    = FracOpHits / TotEvents;
     std::cout << "At the end of the run, I am putting the following into the TTree:\n"
               << RunNumber << " " << TotEvents << " " << NumberOfRCEs << " " << FracRCEs << " " << AvADCPedDiff << " " << FracWaveforms << " " << FracOpHits << " "
               << nPTBTrigsOn110 << " " << nPTBTrigsOn111 << " " << nPTBTrigsOn112 << " " << nPTBTrigsOn113 << " " << nPTBTrigsOn114 << " " << nPTBTrigsOn115 << " " << EvJustSSPs
               << std::endl;
     fTree->Fill();
   }
 }

void DAQToOffline::GoodRun::printParameterSet ( )

Member Data Documentation

std::vector<float> DAQToOffline::GoodRun::ADCVec

private

Definition at line 87 of file GoodRun_module.cc.

double DAQToOffline::GoodRun::AvADCPedDiff

private

Definition at line 67 of file GoodRun_module.cc.

std::vector<int > DAQToOffline::GoodRun::ChanNum

private

Definition at line 86 of file GoodRun_module.cc.

unsigned int DAQToOffline::GoodRun::DigSize

private

Definition at line 84 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::Event

private

Definition at line 83 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::EvJustSSPs

private

Definition at line 76 of file GoodRun_module.cc.

std::string DAQToOffline::GoodRun::fCounterModuleLabel

private

Definition at line 90 of file GoodRun_module.cc.

TTree* DAQToOffline::GoodRun::FlatTree

private

Definition at line 82 of file GoodRun_module.cc.

bool DAQToOffline::GoodRun::fMakeFlatTree

private

Definition at line 89 of file GoodRun_module.cc.

short unsigned int DAQToOffline::GoodRun::fMaxSamples

private

Definition at line 91 of file GoodRun_module.cc.

std::string DAQToOffline::GoodRun::fOpHitModuleLabel

private

Definition at line 90 of file GoodRun_module.cc.

double DAQToOffline::GoodRun::FracOpHits

private

Definition at line 69 of file GoodRun_module.cc.

double DAQToOffline::GoodRun::FracRCEs

private

Definition at line 66 of file GoodRun_module.cc.

double DAQToOffline::GoodRun::FracWaveforms

private

Definition at line 68 of file GoodRun_module.cc.

std::string DAQToOffline::GoodRun::fRawDigitModuleLabel

private

Definition at line 90 of file GoodRun_module.cc.

TTree* DAQToOffline::GoodRun::fTree

private

Definition at line 62 of file GoodRun_module.cc.

std::string DAQToOffline::GoodRun::fWaveformModuleLabel

private

Definition at line 90 of file GoodRun_module.cc.

std::vector<int> DAQToOffline::GoodRun::MyUsefulChans

private

Definition at line 93 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::nPTBTrigsOn110

private

Definition at line 70 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::nPTBTrigsOn111

private

Definition at line 71 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::nPTBTrigsOn112

private

Definition at line 72 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::nPTBTrigsOn113

private

Definition at line 73 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::nPTBTrigsOn114

private

Definition at line 74 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::nPTBTrigsOn115

private

Definition at line 75 of file GoodRun_module.cc.

unsigned int DAQToOffline::GoodRun::NSamples

private

Definition at line 85 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::NumberOfRCEs

private

Definition at line 65 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::RunNumber

private

Definition at line 63 of file GoodRun_module.cc.

int DAQToOffline::GoodRun::TotEvents

private

Definition at line 64 of file GoodRun_module.cc.

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

dunetpc/dune/daqinput35t/GoodRun_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation