opdet::OpHitAna Class Reference

Inheritance diagram for opdet::OpHitAna:

Public Member Functions
	OpHitAna (const fhicl::ParameterSet &)
void	analyze (const art::Event &)

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

Detailed Description

Definition at line 33 of file OpHitAna_module.cc.

Constructor & Destructor Documentation

opdet::OpHitAna::OpHitAna

(

const fhicl::ParameterSet &

pset

)

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:

• OpHitAna_module.cc