opdet::FullOpHitFinder Class Reference

Inheritance diagram for opdet::FullOpHitFinder:

Public Member Functions
	FullOpHitFinder (const fhicl::ParameterSet &)
virtual	~FullOpHitFinder ()
void	produce (art::Event &)
void	beginJob ()
void	endJob ()

Private Member Functions
std::map< int, int >	GetChannelMap ()
std::vector< double >	GetSPEScales ()
std::vector< double >	GetSPEShifts ()

Private Attributes
std::string	fInputModule
std::vector< std::string >	fInputLabels
std::set< unsigned int >	fChannelMasks
pmtana::PulseRecoManager	fPulseRecoMgr
pmtana::PMTPulseRecoBase *	fThreshAlg
pmtana::PMTPedestalBase *	fPedAlg
std::string	_output_filename
std::vector< double >	_tstart_v
std::vector< double >	_tmax_v
std::vector< double >	_tend_v
std::vector< double >	_tcross_v
std::vector< double >	_amp_v
std::vector< double >	_area_v
std::vector< double >	_ped_mean_v
std::vector< double >	_ped_sigma_v
std::vector< double >	_wf
double	_tstart
int	_run
int	_event
int	_ch
TTree *	_wftree
TTree *	_hittree
TTree *	_geotree
TFile *	_ofile

Detailed Description

Definition at line 48 of file FullOpHitFinder_module.cc.

Constructor & Destructor Documentation

opdet::FullOpHitFinder::FullOpHitFinder ( const fhicl::ParameterSet &  pset )

explicit

Definition at line 101 of file FullOpHitFinder_module.cc.

                                                                :
    EDProducer{pset},
    fPulseRecoMgr()
  {

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

    for (auto const& ch : pset.get< std::vector< unsigned int > >
                            ("ChannelMasks", std::vector< unsigned int >()))
      fChannelMasks.insert(ch);

    // Initialize the hit finder algorithm
    auto const hit_alg_pset = pset.get< fhicl::ParameterSet >("HitAlgoPset");
    std::string threshAlgName = hit_alg_pset.get< std::string >("Name");
    if      (threshAlgName == "Threshold")
      fThreshAlg = new pmtana::AlgoThreshold(hit_alg_pset);
    else if (threshAlgName == "SiPM")
      fThreshAlg = new pmtana::AlgoSiPM(hit_alg_pset);
    else if (threshAlgName == "SlidingWindow")
      fThreshAlg = new pmtana::AlgoSlidingWindow(hit_alg_pset);
    else if (threshAlgName == "FixedWindow")
      fThreshAlg = new pmtana::AlgoFixedWindow(hit_alg_pset);
    else if (threshAlgName == "CFD" )
      fThreshAlg = new pmtana::AlgoCFD(hit_alg_pset);
    else throw art::Exception(art::errors::UnimplementedFeature)
                    << "Cannot find implementation for "
                    << threshAlgName << " algorithm.\n";

    auto const ped_alg_pset = pset.get< fhicl::ParameterSet >("PedAlgoPset");
    std::string pedAlgName = ped_alg_pset.get< std::string >("Name");
    if      (pedAlgName == "Edges")
      fPedAlg = new pmtana::PedAlgoEdges(ped_alg_pset);
    else if (pedAlgName == "RollingMean")
      fPedAlg = new pmtana::PedAlgoRollingMean(ped_alg_pset);
    else if (pedAlgName == "UB"   )
      fPedAlg = new pmtana::PedAlgoUB(ped_alg_pset);
    else throw art::Exception(art::errors::UnimplementedFeature)
                    << "Cannot find implementation for "
                    << pedAlgName << " algorithm.\n";

    fPulseRecoMgr.AddRecoAlgo(fThreshAlg);
    fPulseRecoMgr.SetDefaultPedAlgo(fPedAlg);

    _output_filename = pset.get<std::string>("OutputFile","out.root");

  }

opdet::FullOpHitFinder::~FullOpHitFinder ( )

virtual

Definition at line 221 of file FullOpHitFinder_module.cc.

  {

    delete fThreshAlg;
    delete fPedAlg;

  }

Member Function Documentation

void opdet::FullOpHitFinder::beginJob

(

)

Definition at line 150 of file FullOpHitFinder_module.cc.

  {
    // analyzie tuple prep
    _ofile = TFile::Open(_output_filename.c_str(),"RECREATE");

    _wftree = new TTree("wftree","wftree");
    _wftree->Branch("run",&_run,"run/I");
    _wftree->Branch("event",&_event,"event/I");
    _wftree->Branch("ch",&_ch,"ch/I");
    _wftree->Branch("wf",&_wf);
    _wftree->Branch("ped_mean",&_ped_mean_v);
    _wftree->Branch("ped_sigma",&_ped_sigma_v);
    _wftree->Branch("tstart",&_tstart,"tstart/D");

    _hittree = new TTree("hittree","hittree");
    _hittree->Branch("run",&_run,"run/I");
    _hittree->Branch("event",&_event,"event/I");
    _hittree->Branch("ch",&_ch,"ch/I");
    _hittree->Branch("tstart",&_tstart_v);
    _hittree->Branch("tmax",&_tmax_v);
    _hittree->Branch("tend",&_tend_v);
    _hittree->Branch("tcross",&_tcross_v);
    _hittree->Branch("amp",&_amp_v);
    _hittree->Branch("area",&_area_v);
    _hittree->Branch("ped_mean",&_ped_mean_v);
    _hittree->Branch("ped_sigma",&_ped_sigma_v);

    _geotree = new TTree("geotree","geotree");
    std::vector<double> pmtX, pmtY, pmtZ;
    std::vector<double> minX, minY, minZ;
    std::vector<double> maxX, maxY, maxZ;
    auto const geop = lar::providerFrom<geo::Geometry>();
    double PMTxyz[3];
    for(size_t opch=0; opch<geop->NOpChannels(); ++opch) {
      geop->OpDetGeoFromOpChannel(opch).GetCenter(PMTxyz);
      pmtX.push_back(PMTxyz[0]);
      pmtY.push_back(PMTxyz[1]);
      pmtZ.push_back(PMTxyz[2]);
    }
    for(auto iter=geop->begin_TPC(); iter!=geop->end_TPC(); ++iter) {
      auto const& tpc = (*iter);
      minX.push_back(tpc.BoundingBox().MinX());
      minY.push_back(tpc.BoundingBox().MinY());
      minZ.push_back(tpc.BoundingBox().MinZ());
      maxX.push_back(tpc.BoundingBox().MaxX());
      maxY.push_back(tpc.BoundingBox().MaxY());
      maxZ.push_back(tpc.BoundingBox().MaxZ());
    }
    _geotree->Branch("pmtX",&pmtX);
    _geotree->Branch("pmtY",&pmtY);
    _geotree->Branch("pmtZ",&pmtZ);
    _geotree->Branch("minX",&minX);
    _geotree->Branch("minY",&minY);
    _geotree->Branch("minZ",&minZ);
    _geotree->Branch("maxX",&maxX);
    _geotree->Branch("maxY",&maxY);
    _geotree->Branch("maxZ",&maxZ);
    _geotree->Fill();
  }

void opdet::FullOpHitFinder::endJob

(

)

Definition at line 210 of file FullOpHitFinder_module.cc.

  {
    _ofile->cd();
    _wftree->Write();
    _hittree->Write();
    _geotree->Write();
    _ofile->Close();
  }

std::map< int, int > opdet::FullOpHitFinder::GetChannelMap ( )

private

std::vector< double > opdet::FullOpHitFinder::GetSPEScales ( )

private

std::vector< double > opdet::FullOpHitFinder::GetSPEShifts ( )

private

void opdet::FullOpHitFinder::produce

(

art::Event &

evt

)

Definition at line 230 of file FullOpHitFinder_module.cc.

  {

    _event = evt.id().event();
    _run   = evt.id().run();

    // These is the storage pointer we will put in the event
    std::unique_ptr< std::vector< recob::OpHit > >
                      HitPtr(new std::vector< recob::OpHit >);

    //
    // Get the pulses from the event
    //

    if(fInputLabels.empty()) fInputLabels.push_back("");
    for (auto label : fInputLabels) {

      art::Handle< std::vector< raw::OpDetWaveform > > wfHandle;
      evt.getByLabel(fInputModule, label, wfHandle);
      if (!wfHandle.isValid()) continue; // Skip non-existent collections

      for(auto const& waveform : (*wfHandle)) {
        
        _ch = static_cast< int >(waveform.ChannelNumber());
        _tstart = waveform.TimeStamp();
        
        _wf.clear();
        _wf.resize(waveform.size());
        for(size_t idx=0; idx<_wf.size(); ++idx)
          _wf[idx] = waveform[idx];
        
        fPulseRecoMgr.Reconstruct(waveform);
        
        // Record waveforms
        _ped_mean_v = fPedAlg->Mean();
        _ped_sigma_v = fPedAlg->Sigma();
        _wftree->Fill();
        
        // Record pulses
        auto const& pulses = fThreshAlg->GetPulses();
        size_t npulse = pulses.size();
        _tstart_v.resize(npulse); _tmax_v.resize(npulse); _tend_v.resize(npulse); _tcross_v.resize(npulse);
        _amp_v.resize(npulse); _area_v.resize(npulse);
        _ped_mean_v.resize(npulse); _ped_sigma_v.resize(npulse);
        
        for(size_t idx=0; idx<npulse; ++idx) {
          auto const& pulse = pulses[idx];
          _tstart_v[idx] = pulse.t_start;
          _tmax_v[idx]   = pulse.t_max;
          _tend_v[idx]   = pulse.t_end;
          _tcross_v[idx] = pulse.t_cfdcross;
          _amp_v[idx]    = pulse.peak;
          _area_v[idx]   = pulse.area;
          _ped_mean_v[idx]  = pulse.ped_mean;
          _ped_sigma_v[idx] = pulse.ped_sigma;
        }
        _hittree->Fill();
      }
    }
  }

Member Data Documentation

std::vector<double> opdet::FullOpHitFinder::_amp_v

private

Definition at line 78 of file FullOpHitFinder_module.cc.

std::vector<double> opdet::FullOpHitFinder::_area_v

private

Definition at line 78 of file FullOpHitFinder_module.cc.

int opdet::FullOpHitFinder::_ch

private

Definition at line 82 of file FullOpHitFinder_module.cc.

int opdet::FullOpHitFinder::_event

private

Definition at line 82 of file FullOpHitFinder_module.cc.

TTree* opdet::FullOpHitFinder::_geotree

private

Definition at line 86 of file FullOpHitFinder_module.cc.

TTree* opdet::FullOpHitFinder::_hittree

private

Definition at line 85 of file FullOpHitFinder_module.cc.

TFile* opdet::FullOpHitFinder::_ofile

private

Definition at line 87 of file FullOpHitFinder_module.cc.

std::string opdet::FullOpHitFinder::_output_filename

private

Definition at line 76 of file FullOpHitFinder_module.cc.

std::vector<double> opdet::FullOpHitFinder::_ped_mean_v

private

Definition at line 79 of file FullOpHitFinder_module.cc.

std::vector<double> opdet::FullOpHitFinder::_ped_sigma_v

private

Definition at line 79 of file FullOpHitFinder_module.cc.

int opdet::FullOpHitFinder::_run

private

Definition at line 82 of file FullOpHitFinder_module.cc.

std::vector<double> opdet::FullOpHitFinder::_tcross_v

private

Definition at line 77 of file FullOpHitFinder_module.cc.

std::vector<double> opdet::FullOpHitFinder::_tend_v

private

Definition at line 77 of file FullOpHitFinder_module.cc.

std::vector<double> opdet::FullOpHitFinder::_tmax_v

private

Definition at line 77 of file FullOpHitFinder_module.cc.

double opdet::FullOpHitFinder::_tstart

private

Definition at line 81 of file FullOpHitFinder_module.cc.

std::vector<double> opdet::FullOpHitFinder::_tstart_v

private

Definition at line 77 of file FullOpHitFinder_module.cc.

std::vector<double> opdet::FullOpHitFinder::_wf

private

Definition at line 80 of file FullOpHitFinder_module.cc.

TTree* opdet::FullOpHitFinder::_wftree

private

Definition at line 84 of file FullOpHitFinder_module.cc.

std::set< unsigned int > opdet::FullOpHitFinder::fChannelMasks

private

Definition at line 69 of file FullOpHitFinder_module.cc.

std::vector< std::string > opdet::FullOpHitFinder::fInputLabels

private

Definition at line 68 of file FullOpHitFinder_module.cc.

std::string opdet::FullOpHitFinder::fInputModule

private

Definition at line 67 of file FullOpHitFinder_module.cc.

pmtana::PMTPedestalBase* opdet::FullOpHitFinder::fPedAlg

private

Definition at line 73 of file FullOpHitFinder_module.cc.

pmtana::PulseRecoManager opdet::FullOpHitFinder::fPulseRecoMgr

private

Definition at line 71 of file FullOpHitFinder_module.cc.

pmtana::PMTPulseRecoBase* opdet::FullOpHitFinder::fThreshAlg

private

Definition at line 72 of file FullOpHitFinder_module.cc.

---

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

• FullOpHitFinder_module.cc