opdet::wvfAna Class Reference

Inheritance diagram for opdet::wvfAna:

Public Member Functions
	wvfAna (fhicl::ParameterSet const &p)
	wvfAna (wvfAna const &)=delete
	wvfAna (wvfAna &&)=delete
wvfAna &	operator= (wvfAna const &)=delete
wvfAna &	operator= (wvfAna &&)=delete
void	analyze (art::Event const &e) override
void	beginJob () override
void	endJob () override

Public Attributes
opdet::sbndPDMapAlg	pdMap

Private Attributes
size_t	fEvNumber
size_t	fChNumber
double	fSampling
double	fSampling_Daphne
double	fStartTime
double	fEndTime
std::string	fInputModuleName
std::vector< std::string >	fOpDetsToPlot
std::stringstream	histname
std::string	opdetType
std::string	opdetElectronics

Detailed Description

Definition at line 52 of file wvfAna_module.cc.

Constructor & Destructor Documentation

opdet::wvfAna::wvfAna ( fhicl::ParameterSet const &  p )

explicit

Definition at line 92 of file wvfAna_module.cc.

    :
    EDAnalyzer(p)  // ,
    // More initializers here.
  {
    fInputModuleName = p.get< std::string >("InputModule" );
    fOpDetsToPlot    = p.get<std::vector<std::string> >("OpDetsToPlot");

    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
    fSampling = clockData.OpticalClock().Frequency(); // MHz
    fSampling_Daphne = p.get<double>("DaphneFrequency" );
  }

opdet::wvfAna::wvfAna ( wvfAna const &  )

delete

opdet::wvfAna::wvfAna ( wvfAna &&  )

delete

Member Function Documentation

void opdet::wvfAna::analyze ( art::Event const &  e )

override

Definition at line 110 of file wvfAna_module.cc.

  {
    // Implementation of required member function here.
    std::cout << "My module on event #" << e.id().event() << std::endl;

    art::ServiceHandle<art::TFileService> tfs;
    fEvNumber = e.id().event();

    art::Handle< std::vector< raw::OpDetWaveform > > waveHandle;
    e.getByLabel(fInputModuleName, waveHandle);

    if(!waveHandle.isValid()) {
      std::cout << Form("Did not find any G4 photons from a producer: %s", "largeant") << std::endl;
    }

    // // example of usage for pdMap.getCollectionWithProperty()
    // //
    // // define a container
    // auto inBoxTwo = pdMap.getCollectionWithProperty("pds_box", 2);
    // // you can cout the whole json object
    // std::cout << "inBoxTwo:\t" << inBoxTwo << "\n";
    // // traverse its components in a loop
    // for (auto const &e: inBoxTwo) {
    //   std::cout << e["pd_type"] << " " << e["channel"] << ' ' << "\n";
    // }

    // // example of usage for pdMap.getCollectionFromCondition()
    // // define a lambda function with the conditions
    // auto subsetCondition = [](auto const& e)->bool
    //   // modify conditions as you want in the curly braces below
    //   {return e["pd_type"] == "pmt_uncoated" && e["tpc"] == 0;};
    // // get the container that satisfies the conditions
    // auto uncoatedsInTPC0 = pdMap.getCollectionFromCondition(subsetCondition);
    // std::cout << "uncoatedsInTPC0.size():\t" << uncoatedsInTPC0.size() << "\n";
    // for(auto const& e:uncoatedsInTPC0){
    //   std::cout << "e:\t" << e << "\n";
    // }

    std::cout << "Number of waveforms: " << waveHandle->size() << std::endl;

    std::cout << "fOpDetsToPlot:\t";
    for (auto const& opdet : fOpDetsToPlot){std::cout << opdet << " ";}
    std::cout << std::endl;

    int hist_id = 0;
    for(auto const& wvf : (*waveHandle)) {
      fChNumber = wvf.ChannelNumber();
      opdetType = pdMap.pdType(fChNumber);
      opdetElectronics = pdMap.electronicsType(fChNumber);
      if (std::find(fOpDetsToPlot.begin(), fOpDetsToPlot.end(), opdetType) == fOpDetsToPlot.end()) {continue;}
      histname.str(std::string());
      histname << "event_" << fEvNumber
               << "_opchannel_" << fChNumber
               << "_" << opdetType
               << "_" << hist_id;

      fStartTime = wvf.TimeStamp(); //in us
      if (opdetElectronics == "daphne"){
                                fEndTime = double(wvf.size()) / fSampling_Daphne + fStartTime;
                        } //in us
                        else{
        fEndTime = double(wvf.size()) / fSampling + fStartTime;
      } //in us

      //Create a new histogram
      TH1D *wvfHist = tfs->make< TH1D >(histname.str().c_str(), TString::Format(";t - %f (#mus);", fStartTime), wvf.size(), fStartTime, fEndTime);
      for(unsigned int i = 0; i < wvf.size(); i++) {
        wvfHist->SetBinContent(i + 1, (double)wvf[i]);
      }
      hist_id++;
    }
  }

void opdet::wvfAna::beginJob ( )

override

Definition at line 105 of file wvfAna_module.cc.

  {

  }

void opdet::wvfAna::endJob ( )

override

Definition at line 183 of file wvfAna_module.cc.

  {
  }

wvfAna& opdet::wvfAna::operator= ( wvfAna const &  )

delete

wvfAna& opdet::wvfAna::operator= ( wvfAna &&  )

delete

Member Data Documentation

size_t opdet::wvfAna::fChNumber

private

Definition at line 76 of file wvfAna_module.cc.

double opdet::wvfAna::fEndTime

private

Definition at line 80 of file wvfAna_module.cc.

size_t opdet::wvfAna::fEvNumber

private

Definition at line 75 of file wvfAna_module.cc.

std::string opdet::wvfAna::fInputModuleName

private

Definition at line 84 of file wvfAna_module.cc.

std::vector<std::string> opdet::wvfAna::fOpDetsToPlot

private

Definition at line 85 of file wvfAna_module.cc.

double opdet::wvfAna::fSampling

private

Definition at line 77 of file wvfAna_module.cc.

double opdet::wvfAna::fSampling_Daphne

private

Definition at line 78 of file wvfAna_module.cc.

double opdet::wvfAna::fStartTime

private

Definition at line 79 of file wvfAna_module.cc.

std::stringstream opdet::wvfAna::histname

private

Definition at line 86 of file wvfAna_module.cc.

std::string opdet::wvfAna::opdetElectronics

private

Definition at line 88 of file wvfAna_module.cc.

std::string opdet::wvfAna::opdetType

private

Definition at line 87 of file wvfAna_module.cc.

opdet::sbndPDMapAlg opdet::wvfAna::pdMap

Definition at line 71 of file wvfAna_module.cc.

---

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

• wvfAna_module.cc