opdet::OpDigiAna Class Reference

Inheritance diagram for opdet::OpDigiAna:

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

Private Attributes
std::string	fInputModule
std::string	fInstanceName
float	fSampleFreq
float	fTimeBegin
float	fTimeEnd
float	fZeroSupThresh
double	fSPEAmp
bool	fMakeBipolarHist
bool	fMakeUnipolarHist

Detailed Description

Definition at line 29 of file OpDigiAna_module.cc.

Constructor & Destructor Documentation

opdet::OpDigiAna::OpDigiAna

(

const fhicl::ParameterSet &

pset

)

Definition at line 67 of file OpDigiAna_module.cc.

    : EDAnalyzer(pset)
  {

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


    // Indicate that the Input Module comes from .fcl
    fInputModule = pset.get<std::string>("InputModule");
    fInstanceName = pset.get<std::string>("InstanceName");
    fZeroSupThresh = pset.get<double>("ZeroSupThresh") * fSPEAmp;

    fMakeBipolarHist    = pset.get<bool>("MakeBipolarHist");
    fMakeUnipolarHist = pset.get<bool>("MakeUnipolarHist");


  }

Member Function Documentation

void opdet::OpDigiAna::analyze

(

const art::Event &

evt

)

Definition at line 90 of file OpDigiAna_module.cc.

  {

    // Create a handle for our vector of pulses
    art::Handle< std::vector< raw::OpDetPulse > > WaveformHandle;

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


    // Read in WaveformHandle
    evt.getByLabel(fInputModule, fInstanceName, WaveformHandle);

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


    std::vector<std::string> histnames;

    // For every OpDet waveform in the vector given by WaveformHandle
    for(unsigned int i = 0; i < WaveformHandle->size(); ++i)
      {
        // Get OpDetPulse
        art::Ptr< raw::OpDetPulse > ThePulsePtr(WaveformHandle, i);
        raw::OpDetPulse ThePulse = *ThePulsePtr;

        // Make an instance of histogram and its pointer, changing all units to ns
        // Notice that histogram axis is in ns but binned by 1/64MHz;
        //   appropriate conversions are made from beginning and end time
        //   in us, and frequency in MHz.

        // Make sure the histogram name is unique since there can be multiple pulses
        // per event and photo detector
        unsigned int pulsenum = 0;
        while (true) {
            sprintf(HistName, "Event_%d_OpDet_%i_Pulse_%i", evt.id().event(), ThePulse.OpChannel(), pulsenum);
            auto p = std::find(histnames.begin(), histnames.end(), HistName);
            if ( p != histnames.end() ) {
                // Found a duplicate
                pulsenum++;
            }
            else {
                // Found a unique name
                histnames.push_back(HistName);
                break;
            }
        }



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

        sprintf(HistName, "Event_%d_uni_OpDet_%i_Pulse_%i", evt.id().event(), ThePulse.OpChannel(), pulsenum);
        TH1D* UnipolarHist = nullptr;
        if(fMakeUnipolarHist)
          {
            UnipolarHist= tfs->make<TH1D>(HistName, ";t (ns);",
                                          int((fTimeEnd - fTimeBegin) * fSampleFreq),
                                          fTimeBegin * 1000.,
                                          fTimeEnd * 1000.);
          }

        for(unsigned int binNum = 0; binNum < ThePulse.Waveform().size(); ++binNum)
          {
            // Fill histogram with waveform

            if(fMakeBipolarHist) PulseHist->SetBinContent( binNum,
                                                         (double) ThePulse.Waveform()[binNum] );
            if((binNum>0)&&(fMakeUnipolarHist))
              {
                double BinContent =     (double) ThePulse.Waveform()[binNum] +
                  (double) UnipolarHist->GetBinContent(binNum-1);
                if(BinContent>fZeroSupThresh)
                  UnipolarHist->SetBinContent(binNum,BinContent);
                else
                  UnipolarHist->SetBinContent(binNum,0);

              }

          }

      }
  }

Member Data Documentation

std::string opdet::OpDigiAna::fInputModule

private

Definition at line 46 of file OpDigiAna_module.cc.

std::string opdet::OpDigiAna::fInstanceName

private

Definition at line 47 of file OpDigiAna_module.cc.

bool opdet::OpDigiAna::fMakeBipolarHist

private

Definition at line 56 of file OpDigiAna_module.cc.

bool opdet::OpDigiAna::fMakeUnipolarHist

private

Definition at line 57 of file OpDigiAna_module.cc.

float opdet::OpDigiAna::fSampleFreq

private

Definition at line 48 of file OpDigiAna_module.cc.

double opdet::OpDigiAna::fSPEAmp

private

Definition at line 54 of file OpDigiAna_module.cc.

float opdet::OpDigiAna::fTimeBegin

private

Definition at line 49 of file OpDigiAna_module.cc.

float opdet::OpDigiAna::fTimeEnd

private

Definition at line 50 of file OpDigiAna_module.cc.

float opdet::OpDigiAna::fZeroSupThresh

private

Definition at line 52 of file OpDigiAna_module.cc.

---

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

• OpDigiAna_module.cc