SBNDOpT0Finder Class Reference

Inheritance diagram for SBNDOpT0Finder:

Public Member Functions
	SBNDOpT0Finder (fhicl::ParameterSet const &p)
	SBNDOpT0Finder (SBNDOpT0Finder const &)=delete
	SBNDOpT0Finder (SBNDOpT0Finder &&)=delete
SBNDOpT0Finder &	operator= (SBNDOpT0Finder const &)=delete
SBNDOpT0Finder &	operator= (SBNDOpT0Finder &&)=delete
void	produce (art::Event &e) override

Private Member Functions
void	DoMatch (art::Event &e, int tpc, std::unique_ptr< std::vector< anab::T0 >> &t0_v, std::unique_ptr< art::Assns< recob::Slice, anab::T0 >> &slice_t0_assn_v, std::unique_ptr< art::Assns< recob::OpFlash, anab::T0 >> &flash_t0_assn_v)
	Performs the matching in a specified tpc. More...
bool	ConstructLightClusters (art::Event &e, unsigned int tpc)
	Constructs all the LightClusters (TPC Objects) in a specified TPC. More...
float	GetNPhotons (const float charge, const art::Ptr< recob::PFParticle > &pfp)
	Returns the number of photons given charge and PFParticle. More...
std::vector< int >	PDNamesToList (std::vector< std::string >)
	Convert from a list of PDS names to a list of op channels. More...
std::vector< int >	GetUncoatedPTMList (std::vector< int > ch_to_use)
	Returns a list of uncoated PMTs that are a subset of those in ch_to_use. More...

Private Attributes
std::unique_ptr < SemiAnalyticalModel >	_semi_model
fhicl::ParameterSet	_vuv_params
fhicl::ParameterSet	_vis_params
::flashmatch::FlashMatchManager	_mgr
	The flash matching manager. More...
std::vector < flashmatch::FlashMatch_t >	_result_v
	Matching result will be stored here. More...
std::vector< std::string >	_opflash_producer_v
	The OpFlash producers (to be set) More...
std::vector< unsigned int >	_tpc_v
	TPC number per OpFlash producer (to be set) More...
std::string	_slice_producer
	The Slice producer (to be set) More...
double	_flash_trange_start
	The time start from where to include flashes (to be set) More...
double	_flash_trange_end
	The time stop from where to stop including flashes (to be set) More...
float	_charge_to_n_photons_track
	The conversion factor betweeen hit integral and photons (to be set) More...
float	_charge_to_n_photons_shower
	The conversion factor betweeen hit integral and photons (to be set) More...
std::vector< std::string >	_photo_detectors
	The photodetector to use (to be set) More...
std::vector< int >	_opch_to_use
	List of opch to use (will be infered from _photo_detectors) More...
std::vector< int >	_uncoated_pmts
	List of uncoated opch to use (will be infered from _opch_to_use) More...
opdet::sbndPDMapAlg	_pds_map
	map for photon detector types More...
std::vector < flashmatch::QCluster_t >	_light_cluster_v
	Vector that contains all the TPC objects. More...
std::map< int, art::Ptr < recob::Slice > >	_clusterid_to_slice
std::map< int, art::Ptr < recob::OpFlash > >	_flashid_to_opflash
	Will contain map tpc object id -> Slice. More...
TTree *	_tree1
	Will contain map flash id -> OpFlash. More...
int	_run
int	_subrun
int	_event
int	_tpc
int	_matchid
int	_flashid
int	_tpcid
double	_t0
double	_score
double	_tpc_xmin
double	_qll_xmin
double	_hypo_pe
double	_flash_pe
std::vector< double >	_flash_spec
std::vector< double >	_hypo_spec
TTree *	_tree2
std::vector< float >	_dep_x
std::vector< float >	_dep_y
std::vector< float >	_dep_z
std::vector< float >	_dep_charge
std::vector< float >	_dep_n_photons
std::vector< int >	_dep_slice

Detailed Description

Definition at line 56 of file SBNDOpT0Finder_module.cc.

Constructor & Destructor Documentation

SBNDOpT0Finder::SBNDOpT0Finder ( fhicl::ParameterSet const &  p )

explicit

Definition at line 137 of file SBNDOpT0Finder_module.cc.

  : EDProducer{p}
{
  produces<std::vector<anab::T0>>();
  produces<art::Assns<recob::Slice, anab::T0>>();
  produces<art::Assns<recob::OpFlash, anab::T0>>();

  ::art::ServiceHandle<geo::Geometry> geo;

  _vuv_params = p.get<fhicl::ParameterSet>("VUVHits");
  _vis_params = p.get<fhicl::ParameterSet>("VIVHits");
  _semi_model = std::make_unique<SemiAnalyticalModel>(_vuv_params, _vis_params, true, false);

  _opflash_producer_v = p.get<std::vector<std::string>>("OpFlashProducers");
  _tpc_v = p.get<std::vector<unsigned int>>("TPCs");
  _slice_producer = p.get<std::string>("SliceProducer");

  _flash_trange_start = p.get<double>("FlashVetoTimeStart", 0);
  _flash_trange_end = p.get<double>("FlashVetoTimeEnd", 2);

  _photo_detectors = p.get<std::vector<std::string>>("PhotoDetectors");
  _opch_to_use = this->PDNamesToList(_photo_detectors);
  _uncoated_pmts = this->GetUncoatedPTMList(_opch_to_use);

  _charge_to_n_photons_track = p.get<float>("ChargeToNPhotonsTrack");
  _charge_to_n_photons_shower = p.get<float>("ChargeToNPhotonsShower");

  if (_tpc_v.size() != _opflash_producer_v.size()) {
    throw cet::exception("SBNDOpT0Finder")
      << "TPC vector and OpFlash producer vector don't have the same size, check your fcl params.";
  }

  _mgr.Configure(p.get<flashmatch::Config_t>("FlashMatchConfig"));

  _mgr.SetChannelMask(_opch_to_use);

  _mgr.SetUncoatedPMTs(_uncoated_pmts);

  _mgr.SetSemiAnalyticalModel(std::move(_semi_model));

  _flash_spec.resize(geo->NOpDets(), 0.);
  _hypo_spec.resize(geo->NOpDets(), 0.);

  art::ServiceHandle<art::TFileService> fs;

  _tree1 = fs->make<TTree>("slice_deposition_tree","");
  _tree1->Branch("run",             &_run,                             "run/I");
  _tree1->Branch("subrun",          &_subrun,                          "subrun/I");
  _tree1->Branch("event",           &_event,                           "event/I");
  _tree1->Branch("dep_slice", "std::vector<int>", &_dep_slice);
  _tree1->Branch("dep_x", "std::vector<float>", &_dep_x);
  _tree1->Branch("dep_y", "std::vector<float>", &_dep_y);
  _tree1->Branch("dep_z", "std::vector<float>", &_dep_z);
  _tree1->Branch("dep_charge", "std::vector<float>", &_dep_charge);
  _tree1->Branch("dep_n_photons", "std::vector<float>", &_dep_n_photons);

  _tree2 = fs->make<TTree>("flash_match_tree","");
  _tree2->Branch("run",             &_run,                             "run/I");
  _tree2->Branch("subrun",          &_subrun,                          "subrun/I");
  _tree2->Branch("event",           &_event,                           "event/I");
  _tree2->Branch("tpc",             &_tpc,                             "tpc/I");
  _tree2->Branch("matchid",         &_matchid,                         "matchid/I");
  _tree2->Branch("tpcid",           &_tpcid,                           "tpcid/I");
  _tree2->Branch("flashid",         &_flashid,                         "flashid/I");
  _tree2->Branch("tpc_xmin",        &_tpc_xmin,                        "tpc_xmin/D");
  _tree2->Branch("qll_xmin",        &_qll_xmin,                        "qll_xmin/D");
  _tree2->Branch("t0",              &_t0,                              "t0/D");
  _tree2->Branch("score",           &_score,                           "score/D");
  _tree2->Branch("hypo_pe",         &_hypo_pe,                         "hypo_pe/D");
  _tree2->Branch("flash_pe",        &_flash_pe,                        "flash_pe/D");
  _tree2->Branch("hypo_spec",       "std::vector<double>",             &_hypo_spec);
  _tree2->Branch("flash_spec",      "std::vector<double>",             &_flash_spec);
}

SBNDOpT0Finder::SBNDOpT0Finder ( SBNDOpT0Finder const &  )

delete

SBNDOpT0Finder::SBNDOpT0Finder ( SBNDOpT0Finder &&  )

delete

Member Function Documentation

bool SBNDOpT0Finder::ConstructLightClusters ( art::Event &  e, unsigned int  tpc  )

private

Constructs all the LightClusters (TPC Objects) in a specified TPC.

Definition at line 402 of file SBNDOpT0Finder_module.cc.

                                                                         {
  // One slice is one QCluster_t.
  // Start from a slice, get all the PFParticles, from there get all the spacepoints, from
  // there get all the hits on the collection plane.
  // Use the charge on the collection plane to estimate the light, and the 3D spacepoint
  // position for the 3D location.

  _light_cluster_v.clear();

  ::art::Handle<std::vector<recob::Slice>> slice_h;
  e.getByLabel(_slice_producer, slice_h);
  if(!slice_h.isValid() || slice_h->empty()) {
    mf::LogWarning("SBNDOpT0Finder") << "Don't have good Slices." << std::endl;
    return false;
  }

  ::art::Handle<std::vector<recob::PFParticle>> pfp_h;
  e.getByLabel(_slice_producer, pfp_h);
  if(!pfp_h.isValid() || pfp_h->empty()) {
    mf::LogWarning("SBNDOpT0Finder") << "Don't have good PFParticle." << std::endl;
    return false;
  }

  ::art::Handle<std::vector<recob::SpacePoint>> spacepoint_h;
  e.getByLabel(_slice_producer, spacepoint_h);
  if(!spacepoint_h.isValid() || spacepoint_h->empty()) {
    mf::LogWarning("SBNDOpT0Finder") << "Don't have good SpacePoint." << std::endl;
    return false;
  }

  // Construct the vector of Slices
  std::vector<art::Ptr<recob::Slice>> slice_v;
  art::fill_ptr_vector(slice_v, slice_h);

  // Get the associations between slice->pfp->spacepoint->hit
  art::FindManyP<recob::PFParticle> slice_to_pfps (slice_h, e, _slice_producer);
  art::FindManyP<recob::SpacePoint> pfp_to_spacepoints (pfp_h, e, _slice_producer);
  art::FindManyP<recob::Hit> spacepoint_to_hits (spacepoint_h, e, _slice_producer);

  // Loop over the Slices
  for (size_t n_slice = 0; n_slice < slice_h->size(); n_slice++) {

    flashmatch::QCluster_t light_cluster;

    _dep_slice.clear();
    _dep_x.clear();
    _dep_y.clear();
    _dep_z.clear();
    _dep_charge.clear();
    _dep_n_photons.clear();

    // Get the associated PFParticles
    std::vector<art::Ptr<recob::PFParticle>> pfp_v = slice_to_pfps.at(n_slice);

    for (size_t n_pfp = 0; n_pfp < pfp_v.size(); n_pfp++) {

      auto pfp = pfp_v[n_pfp];

      // Get the associated SpacePoints
      std::vector<art::Ptr<recob::SpacePoint>> spacepoint_v = pfp_to_spacepoints.at(pfp.key());

      for (size_t n_spacepoint = 0; n_spacepoint < spacepoint_v.size(); n_spacepoint++) {

        auto spacepoint = spacepoint_v[n_spacepoint];

        // Get the associated hits
        std::vector<art::Ptr<recob::Hit>> hit_v = spacepoint_to_hits.at(spacepoint.key());

        for (size_t n_hit = 0; n_hit < hit_v.size(); n_hit++) {

          auto hit = hit_v[n_hit];

          // Only select hits from the collection plane
          if (hit->View() != geo::kZ) {
            continue;
          }

          // Only use hits (and so spacepoints) that are in the specified TPC
          if (hit->WireID().TPC != tpc) {
            continue;
          }

          const auto &position(spacepoint->XYZ());
          const auto charge(hit->Integral());

          // Emplace this point with charge to the light cluster
          light_cluster.emplace_back(position[0],
                                     position[1],
                                     position[2],
                                     GetNPhotons(charge, pfp));

          // Also save the quantites for the output tree
          _dep_slice.push_back(_light_cluster_v.size());
          _dep_x.push_back(position[0]);
          _dep_y.push_back(position[1]);
          _dep_z.push_back(position[2]);
          _dep_charge.push_back(charge);
          _dep_n_photons.push_back(GetNPhotons(charge, pfp));
        }
      } // End loop over Spacepoints
    } // End loop over PFParticle

    _tree1->Fill();

    // Don't include clusters with zero points
    if (!light_cluster.size()) {
      continue;
    }

    // Save the light cluster, and remember the correspondance from index to slice
    _clusterid_to_slice[_light_cluster_v.size()] = slice_v.at(n_slice);
    _light_cluster_v.emplace_back(light_cluster);

  } // End loop over Slices

  return true;
}

void SBNDOpT0Finder::DoMatch ( art::Event &  e, int  tpc, std::unique_ptr< std::vector< anab::T0 >> &  t0_v, std::unique_ptr< art::Assns< recob::Slice, anab::T0 >> &  slice_t0_assn_v, std::unique_ptr< art::Assns< recob::OpFlash, anab::T0 >> &  flash_t0_assn_v  )

private

Performs the matching in a specified tpc.

Definition at line 250 of file SBNDOpT0Finder_module.cc.

                                                                                               {

  _flashid_to_opflash.clear();
  _clusterid_to_slice.clear();

  auto const & flash_h = e.getValidHandle<std::vector<recob::OpFlash>>(_opflash_producer_v[tpc]);
  if(!flash_h.isValid() || flash_h->empty()) {
    mf::LogInfo("SBNDOpT0Finder") << "Don't have good flashes from producer "
                                  << _opflash_producer_v[tpc] << std::endl;
    return;
  }

  // Construct the vector of OpFlashes
  std::vector<art::Ptr<recob::OpFlash>> flash_v;
  art::fill_ptr_vector(flash_v, flash_h);

  ::art::ServiceHandle<geo::Geometry> geo;

  int n_flashes = 0;
  std::vector<::flashmatch::Flash_t> all_flashes;

  for (size_t n = 0; n < flash_v.size(); n++) {

    auto const& flash = *flash_v[n];

    mf::LogDebug("SBNDOpT0Finder") << "Flash time from " << _opflash_producer_v[tpc] << ": " << flash.Time() << std::endl;

    if(flash.Time() < _flash_trange_start || _flash_trange_end < flash.Time()) {
      continue;
    }

    _flashid_to_opflash[n_flashes] = flash_v[n];

    n_flashes++;

    // Construct a Flash_t
    ::flashmatch::Flash_t f;
    f.x = f.x_err = 0;
    f.pe_v.resize(geo->NOpDets());
    f.pe_err_v.resize(geo->NOpDets());
    for (unsigned int op_ch = 0; op_ch < f.pe_v.size(); op_ch++) {
      unsigned int opdet = geo->OpDetFromOpChannel(op_ch);
      if (std::find(_opch_to_use.begin(), _opch_to_use.end(), op_ch) == _opch_to_use.end()) {
        f.pe_v[opdet] = 0;
        f.pe_err_v[opdet] = 0;
      } else {
        f.pe_v[opdet] = flash.PE(op_ch);
        f.pe_err_v[opdet] = sqrt(flash.PE(op_ch));
      }
    }
    f.y = flash.YCenter();
    f.z = flash.ZCenter();
    f.y_err = flash.YWidth();
    f.z_err = flash.ZWidth();
    f.time = flash.Time();
    f.idx = n_flashes-1;
    all_flashes.resize(n_flashes);
    all_flashes[n_flashes-1] = f;

  } // flash loop

  // Don't waste time if there are no flashes
  if (n_flashes == 0) {
    mf::LogInfo("SBNDOpT0Finder") << "Zero good flashes in this event." << std::endl;
    return;
  }

  // Get all the ligh clusters
  // auto light_cluster_v = GetLighClusters(e);
  if (!ConstructLightClusters(e, tpc)) {
    mf::LogInfo("SBNDOpT0Finder") << "Cannot construct Light Clusters." << std::endl;
    return;
  }

  // Don't waste time if there are no clusters
  if (!_light_cluster_v.size()) {
    mf::LogInfo("SBNDOpT0Finder") << "No slices to work with in TPC " << tpc << "." << std::endl;
    return;
  }

  // Emplace flashes to Flash Matching Manager
  for (auto f : all_flashes) {
    _mgr.Emplace(std::move(f));
  }

  // Emplace clusters to Flash Matching Manager
  for (auto lc : _light_cluster_v) {
    _mgr.Emplace(std::move(lc));
  }

  // Run the matching
  _result_v = _mgr.Match();

  // Loop over the matching results
  for(_matchid = 0; _matchid < (int)(_result_v.size()); ++_matchid) {

    auto const& match = _result_v[_matchid];

    _tpcid    = match.tpc_id;
    _flashid  = match.flash_id;
    _score    = match.score;
    _qll_xmin = match.tpc_point.x;

    mf::LogInfo("SBNDOpT0Finder") << "Matched TPC object " << _tpcid
                                  << " with flash number " << _flashid
                                  << " in TPC " << tpc
                                  << " -> score: " << _score
                                  << ", qll xmin: " << _qll_xmin << std::endl;

    // Get the minimum x position of the TPC Object
    _tpc_xmin = 1.e4;
    for(auto const& pt : _mgr.QClusterArray()[_tpcid]) {
      if(pt.x < _tpc_xmin) _tpc_xmin = pt.x;
    }

    // Get the matched flash time, the t0
    auto const& flash = _mgr.FlashArray()[_flashid];
    _t0 = flash.time;

    // Save the reconstructed flash and hypothesis flash PE spectrum
    if(_hypo_spec.size() != match.hypothesis.size()) {
      throw cet::exception("SBNDOpT0Finder") << "Hypothesis size mismatch!";
    }
    for(size_t pmt=0; pmt<_hypo_spec.size(); ++pmt) _hypo_spec[pmt]  = match.hypothesis[pmt];
    for(size_t pmt=0; pmt<_hypo_spec.size(); ++pmt) _flash_spec[pmt] = flash.pe_v[pmt];

    // Also save the total number of photoelectrons
    _flash_pe = 0.;
    _hypo_pe  = 0.;
    for(auto const& v : _hypo_spec) _hypo_pe += v;
    for(auto const& v : _flash_spec) _flash_pe += v;

    _tree2->Fill();

    // Construct the anab::T0 dataproduc to put in the Event
    auto t0 = anab::T0(_t0,        // "Time": The recontructed flash time, or t0
                       _flash_pe,  // "TriggerType": placing the reconstructed total PE instead
                       _tpcid,     // "TriggerBits": placing the tpc id instead
                       _flashid,   // "ID": placing the flash id instead
                       _score);    // "TriggerConfidence": Matching score

    t0_v->push_back(t0);
    util::CreateAssn(*this, e, *t0_v, _clusterid_to_slice[_tpcid], *slice_t0_assn_v);
    util::CreateAssn(*this, e, *t0_v, _flashid_to_opflash[_flashid], *flash_t0_assn_v);
  }

}

float SBNDOpT0Finder::GetNPhotons ( const float  charge, const art::Ptr< recob::PFParticle > &  pfp  )

private

Returns the number of photons given charge and PFParticle.

Definition at line 520 of file SBNDOpT0Finder_module.cc.

                                                                      {
  return charge * (::lar_pandora::LArPandoraHelper::IsTrack(pfp) ? _charge_to_n_photons_track
                                                                 : _charge_to_n_photons_shower);
}

std::vector< int > SBNDOpT0Finder::GetUncoatedPTMList ( std::vector< int >  ch_to_use )

private

Returns a list of uncoated PMTs that are a subset of those in ch_to_use.

Definition at line 539 of file SBNDOpT0Finder_module.cc.

                                                                          {
  std::vector<int> out_v;

  for (auto ch : ch_to_use) {
    if (_pds_map.isPDType(ch, "pmt_uncoated")) {
      out_v.push_back(ch);
    }
  }

  return out_v;
}

SBNDOpT0Finder& SBNDOpT0Finder::operator= ( SBNDOpT0Finder const &  )

delete

SBNDOpT0Finder& SBNDOpT0Finder::operator= ( SBNDOpT0Finder &&  )

delete

std::vector< int > SBNDOpT0Finder::PDNamesToList ( std::vector< std::string >  pd_names )

private

Convert from a list of PDS names to a list of op channels.

Definition at line 526 of file SBNDOpT0Finder_module.cc.

                                                                          {

  std::vector<int> out_ch_v;

  for (auto name : pd_names) {
    auto ch_v = _pds_map.getChannelsOfType(name);
    out_ch_v.insert(out_ch_v.end(), ch_v.begin(), ch_v.end());
  }

  return out_ch_v;

}

void SBNDOpT0Finder::produce ( art::Event &  e )

override

Definition at line 211 of file SBNDOpT0Finder_module.cc.

{
  std::unique_ptr<std::vector<anab::T0>> t0_v (new std::vector<anab::T0>);
  std::unique_ptr< art::Assns<recob::Slice, anab::T0>> slice_t0_assn_v (new art::Assns<recob::Slice, anab::T0>);
  std::unique_ptr< art::Assns<recob::OpFlash, anab::T0>> flash_t0_assn_v (new art::Assns<recob::OpFlash, anab::T0>);

  // _mgr.PrintConfig();

  _run    = e.id().run();
  _subrun = e.id().subRun();
  _event  = e.id().event();

  // Loop over the specified TPCs
  for (auto tpc : _tpc_v) {

    mf::LogInfo("SBNDOpT0Finder") << "Performing matching in TPC " << tpc << std::endl;

    // Reset the manager and the result vector
    _mgr.Reset();
    _result_v.clear();
    _tpc = tpc;

    // Tell the manager what TPC and cryostat we are going to be doing
    // the matching in. For SBND, the cryostat is always zero.
    _mgr.SetTPCCryo(tpc, 0);

    // Perform the matching in the specified TPC
    DoMatch(e, tpc, t0_v, slice_t0_assn_v, flash_t0_assn_v);

  }

  // Finally, place the anab::T0 vector and the associations in the Event
  e.put(std::move(t0_v));
  e.put(std::move(slice_t0_assn_v));
  e.put(std::move(flash_t0_assn_v));

  return;
}

Member Data Documentation

float SBNDOpT0Finder::_charge_to_n_photons_shower

private

The conversion factor betweeen hit integral and photons (to be set)

Definition at line 107 of file SBNDOpT0Finder_module.cc.

float SBNDOpT0Finder::_charge_to_n_photons_track

private

The conversion factor betweeen hit integral and photons (to be set)

Definition at line 106 of file SBNDOpT0Finder_module.cc.

std::map<int, art::Ptr<recob::Slice> > SBNDOpT0Finder::_clusterid_to_slice

private

Definition at line 118 of file SBNDOpT0Finder_module.cc.

std::vector<float> SBNDOpT0Finder::_dep_charge

private

Definition at line 132 of file SBNDOpT0Finder_module.cc.

std::vector<float> SBNDOpT0Finder::_dep_n_photons

private

Definition at line 132 of file SBNDOpT0Finder_module.cc.

std::vector<int> SBNDOpT0Finder::_dep_slice

private

Definition at line 133 of file SBNDOpT0Finder_module.cc.

std::vector<float> SBNDOpT0Finder::_dep_x

private

Definition at line 132 of file SBNDOpT0Finder_module.cc.

std::vector<float> SBNDOpT0Finder::_dep_y

private

Definition at line 132 of file SBNDOpT0Finder_module.cc.

std::vector<float> SBNDOpT0Finder::_dep_z

private

Definition at line 132 of file SBNDOpT0Finder_module.cc.

int SBNDOpT0Finder::_event

private

Definition at line 122 of file SBNDOpT0Finder_module.cc.

double SBNDOpT0Finder::_flash_pe

private

Definition at line 127 of file SBNDOpT0Finder_module.cc.

std::vector<double> SBNDOpT0Finder::_flash_spec

private

Definition at line 128 of file SBNDOpT0Finder_module.cc.

double SBNDOpT0Finder::_flash_trange_end

private

The time stop from where to stop including flashes (to be set)

Definition at line 104 of file SBNDOpT0Finder_module.cc.

double SBNDOpT0Finder::_flash_trange_start

private

The time start from where to include flashes (to be set)

Definition at line 103 of file SBNDOpT0Finder_module.cc.

int SBNDOpT0Finder::_flashid

private

Definition at line 124 of file SBNDOpT0Finder_module.cc.

std::map<int, art::Ptr<recob::OpFlash> > SBNDOpT0Finder::_flashid_to_opflash

private

Will contain map tpc object id -> Slice.

Definition at line 119 of file SBNDOpT0Finder_module.cc.

double SBNDOpT0Finder::_hypo_pe

private

Definition at line 127 of file SBNDOpT0Finder_module.cc.

std::vector<double> SBNDOpT0Finder::_hypo_spec

private

Definition at line 129 of file SBNDOpT0Finder_module.cc.

std::vector<flashmatch::QCluster_t> SBNDOpT0Finder::_light_cluster_v

private

Vector that contains all the TPC objects.

Definition at line 116 of file SBNDOpT0Finder_module.cc.

int SBNDOpT0Finder::_matchid

private

Definition at line 124 of file SBNDOpT0Finder_module.cc.

::flashmatch::FlashMatchManager SBNDOpT0Finder::_mgr

private

The flash matching manager.

Definition at line 96 of file SBNDOpT0Finder_module.cc.

std::vector<int> SBNDOpT0Finder::_opch_to_use

private

List of opch to use (will be infered from _photo_detectors)

Definition at line 110 of file SBNDOpT0Finder_module.cc.

std::vector<std::string> SBNDOpT0Finder::_opflash_producer_v

private

The OpFlash producers (to be set)

Definition at line 99 of file SBNDOpT0Finder_module.cc.

opdet::sbndPDMapAlg SBNDOpT0Finder::_pds_map

private

map for photon detector types

Definition at line 113 of file SBNDOpT0Finder_module.cc.

std::vector<std::string> SBNDOpT0Finder::_photo_detectors

private

The photodetector to use (to be set)

Definition at line 109 of file SBNDOpT0Finder_module.cc.

double SBNDOpT0Finder::_qll_xmin

private

Definition at line 126 of file SBNDOpT0Finder_module.cc.

std::vector<flashmatch::FlashMatch_t> SBNDOpT0Finder::_result_v

private

Matching result will be stored here.

Definition at line 97 of file SBNDOpT0Finder_module.cc.

int SBNDOpT0Finder::_run

private

Definition at line 122 of file SBNDOpT0Finder_module.cc.

double SBNDOpT0Finder::_score

private

Definition at line 125 of file SBNDOpT0Finder_module.cc.

std::unique_ptr<SemiAnalyticalModel> SBNDOpT0Finder::_semi_model

private

Definition at line 92 of file SBNDOpT0Finder_module.cc.

std::string SBNDOpT0Finder::_slice_producer

private

The Slice producer (to be set)

Definition at line 101 of file SBNDOpT0Finder_module.cc.

int SBNDOpT0Finder::_subrun

private

Definition at line 122 of file SBNDOpT0Finder_module.cc.

double SBNDOpT0Finder::_t0

private

Definition at line 125 of file SBNDOpT0Finder_module.cc.

int SBNDOpT0Finder::_tpc

private

Definition at line 123 of file SBNDOpT0Finder_module.cc.

std::vector<unsigned int> SBNDOpT0Finder::_tpc_v

private

TPC number per OpFlash producer (to be set)

Definition at line 100 of file SBNDOpT0Finder_module.cc.

double SBNDOpT0Finder::_tpc_xmin

private

Definition at line 126 of file SBNDOpT0Finder_module.cc.

int SBNDOpT0Finder::_tpcid

private

Definition at line 124 of file SBNDOpT0Finder_module.cc.

TTree* SBNDOpT0Finder::_tree1

private

Will contain map flash id -> OpFlash.

Definition at line 121 of file SBNDOpT0Finder_module.cc.

TTree* SBNDOpT0Finder::_tree2

private

Definition at line 131 of file SBNDOpT0Finder_module.cc.

std::vector<int> SBNDOpT0Finder::_uncoated_pmts

private

List of uncoated opch to use (will be infered from _opch_to_use)

Definition at line 111 of file SBNDOpT0Finder_module.cc.

fhicl::ParameterSet SBNDOpT0Finder::_vis_params

private

Definition at line 94 of file SBNDOpT0Finder_module.cc.

fhicl::ParameterSet SBNDOpT0Finder::_vuv_params

private

Definition at line 93 of file SBNDOpT0Finder_module.cc.

---

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

• SBNDOpT0Finder_module.cc