opdet::OpHitAna Class Reference

Inheritance diagram for opdet::OpHitAna:

Public Member Functions
	OpHitAna (const fhicl::ParameterSet &)
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 Attributes
std::string	fInputModule
float	fSampleFreq
float	fTimeBegin
float	fTimeEnd
bool	fMakeHistPerChannel
bool	fMakeHistAllChannels
bool	fMakeHitTree
TTree *	fHitTree
Int_t	fEventID
Int_t	fOpChannel
Float_t	fPeakTime
Float_t	fNPe

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 33 of file OpHitAna_module.cc.

Constructor & Destructor Documentation

opdet::OpHitAna::OpHitAna

(

const fhicl::ParameterSet &

)

Definition at line 87 of file OpHitAna_module.cc.

    : EDAnalyzer(pset)
  {

    // Indicate that the Input Module comes from .fcl
    fInputModule = pset.get<std::string>("InputModule");

    art::ServiceHandle<OpDigiProperties> odp;
    fTimeBegin  = odp->TimeBegin();
    fTimeEnd    = odp->TimeEnd();
    fSampleFreq = odp->SampleFreq();


    fMakeHistPerChannel     = pset.get<bool>("MakeHistPerChannel");
    fMakeHistAllChannels    = pset.get<bool>("MakeHistAllChannels");
    fMakeHitTree     = pset.get<bool>("MakeHitTree");

    // If required, make TTree for output

    if(fMakeHitTree)
      {
        art::ServiceHandle<art::TFileService const> tfs;
        fHitTree = tfs->make<TTree>("HitTree","HitTree");
        fHitTree->Branch("EventID",   &fEventID,   "EventID/I");
        fHitTree->Branch("OpChannel", &fOpChannel, "OpChannel/I");
        fHitTree->Branch("PeakTime",  &fPeakTime,  "PeakTime/F");
        fHitTree->Branch("NPe",       &fNPe,       "NPe/F");
      }


  }

Member Function Documentation

void opdet::OpHitAna::analyze

(

const art::Event &

evt

)

Definition at line 120 of file OpHitAna_module.cc.

  {

    // Create a handle for our vector of pulses
    art::Handle< std::vector< recob::OpHit > > HitHandle;

    // Create string for histogram name
    char HistName[50];


    // Read in HitHandle
    evt.getByLabel(fInputModule, HitHandle);

    // Access ART's TFileService, which will handle creating and writing
    // histograms for us.
    art::ServiceHandle<art::TFileService const> tfs;

    art::ServiceHandle<geo::Geometry const> geom;
    int NOpChannels = geom->NOpChannels();

    std::vector<TH1D*> HitHist;

    sprintf(HistName, "Event %d AllOpDets", evt.id().event());

    TH1D * AllHits = nullptr;
    if(fMakeHistAllChannels)
      {
        AllHits = tfs->make<TH1D>(HistName, ";t (ns);",
                                  int((fTimeEnd - fTimeBegin) * fSampleFreq),
                                  fTimeBegin * 1000.,
                                  fTimeEnd * 1000.);
      }

    for(int i=0; i!=NOpChannels; ++i)
      {

        sprintf(HistName, "Event %d OpDet %i", evt.id().event(), i);
        if(fMakeHistPerChannel) HitHist.push_back ( tfs->make<TH1D>(HistName, ";t (ns);",
                                                              int((fTimeEnd - fTimeBegin) * fSampleFreq),
                                                              fTimeBegin * 1000.,
                                                              fTimeEnd * 1000.));

      }


    fEventID = evt.id().event();


    // For every OpHit in the vector
    for(unsigned int i = 0; i < HitHandle->size(); ++i)
      {
        // Get OpHit
        art::Ptr< recob::OpHit > TheHitPtr(HitHandle, i);
        recob::OpHit TheHit = *TheHitPtr;

        fOpChannel = TheHit.OpChannel();
        fNPe       = TheHit.PE();
        fPeakTime  = TheHit.PeakTime();

        if(fMakeHitTree)
          fHitTree->Fill();

        if(fMakeHistPerChannel)
          {
            if(fOpChannel>int(HitHist.size()))
              {
                mf::LogError("OpHitAna")<<"Error : Trying to fill channel out of range, " << fOpChannel;
              }
            HitHist[fOpChannel]->Fill(fPeakTime,fNPe);
          }

        if(fMakeHistAllChannels) AllHits->Fill(fPeakTime, fNPe);

      }

  }

Member Data Documentation

Int_t opdet::OpHitAna::fEventID

private

Definition at line 61 of file OpHitAna_module.cc.

TTree* opdet::OpHitAna::fHitTree

private

Definition at line 60 of file OpHitAna_module.cc.

std::string opdet::OpHitAna::fInputModule

private

Definition at line 48 of file OpHitAna_module.cc.

bool opdet::OpHitAna::fMakeHistAllChannels

private

Definition at line 56 of file OpHitAna_module.cc.

bool opdet::OpHitAna::fMakeHistPerChannel

private

Definition at line 55 of file OpHitAna_module.cc.

bool opdet::OpHitAna::fMakeHitTree

private

Definition at line 57 of file OpHitAna_module.cc.

Float_t opdet::OpHitAna::fNPe

private

Definition at line 64 of file OpHitAna_module.cc.

Int_t opdet::OpHitAna::fOpChannel

private

Definition at line 62 of file OpHitAna_module.cc.

Float_t opdet::OpHitAna::fPeakTime

private

Definition at line 63 of file OpHitAna_module.cc.

float opdet::OpHitAna::fSampleFreq

private

Definition at line 49 of file OpHitAna_module.cc.

float opdet::OpHitAna::fTimeBegin

private

Definition at line 50 of file OpHitAna_module.cc.

float opdet::OpHitAna::fTimeEnd

private

Definition at line 51 of file OpHitAna_module.cc.

---

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

• larana/larana/OpticalDetector/OpHitAna_module.cc