SBNDFlashAna Class Reference

Inheritance diagram for SBNDFlashAna:

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

Private Attributes
std::vector< std::string >	_flash_label_v
std::vector< std::string >	_ophit_label_v
std::vector< std::string >	_mcflash_label_v
std::vector< int >	_tpc_v
opdet::sbndPDMapAlg	_pds_map
TTree *	_flash_tree
int	_run
int	_subrun
int	_event
int	_tpc
double	_flash_time
double	_flash_total_pe
std::vector< double >	_flash_pe_v
double	_flash_y
double	_flash_yerr
double	_flash_z
double	_flash_zerr
std::vector< double >	_flash_ophit_time
std::vector< double >	_flash_ophit_amp
std::vector< double >	_flash_ophit_area
std::vector< double >	_flash_ophit_width
std::vector< double >	_flash_ophit_pe
std::vector< int >	_flash_ophit_ch
std::vector< std::string >	_flash_ophit_chname
TTree *	_mcflash_tree

Detailed Description

Definition at line 31 of file SBNDFlashAna_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 77 of file SBNDFlashAna_module.cc.

  : EDAnalyzer{p}
{
  _flash_label_v   = p.get<std::vector<std::string>>("FlashProducer",   _flash_label_v);
  _mcflash_label_v = p.get<std::vector<std::string>>("MCFlashProducer", _mcflash_label_v);
  _ophit_label_v   = p.get<std::vector<std::string>>("OpHitProducer", _ophit_label_v);
  _tpc_v           = p.get<std::vector<int>>("TPC", _tpc_v);

  art::ServiceHandle<art::TFileService> fs;
  _flash_tree = fs->make<TTree>("FlashTree","");
  _flash_tree->Branch("run", &_run, "run/I");
  _flash_tree->Branch("subrun", &_subrun, "subrun/I");
  _flash_tree->Branch("event", &_event, "event/I");
  _flash_tree->Branch("tpc", &_tpc, "tpc/I");
  _flash_tree->Branch("flash_time", &_flash_time, "flash_time/D");
  _flash_tree->Branch("flash_total_pe", &_flash_total_pe, "flash_total_pe/D");
  _flash_tree->Branch("flash_y", &_flash_y, "flash_y/D");
  _flash_tree->Branch("flash_yerr", &_flash_yerr, "flash_yerr/D");
  _flash_tree->Branch("flash_z", &_flash_z, "flash_z/D");
  _flash_tree->Branch("flash_zerr", &_flash_zerr, "flash_zerr/D");
  _flash_tree->Branch("flash_pe_v", "std::vector<double>", &_flash_pe_v);
  _flash_tree->Branch("flash_ophit_time", "std::vector<double>", &_flash_ophit_time);
  _flash_tree->Branch("flash_ophit_amp", "std::vector<double>", &_flash_ophit_amp);
  _flash_tree->Branch("flash_ophit_area", "std::vector<double>", &_flash_ophit_area);
  _flash_tree->Branch("flash_ophit_width", "std::vector<double>", &_flash_ophit_width);
  _flash_tree->Branch("flash_ophit_pe", "std::vector<double>", &_flash_ophit_pe);
  _flash_tree->Branch("flash_ophit_ch", "std::vector<int>", &_flash_ophit_ch);
  _flash_tree->Branch("flash_ophit_chname", "std::vector<std::string>", &_flash_ophit_chname);

  _mcflash_tree = fs->make<TTree>("MCFlashTree","");
  _mcflash_tree->Branch("run", &_run, "run/I");
  _mcflash_tree->Branch("subrun", &_subrun, "subrun/I");
  _mcflash_tree->Branch("event", &_event, "event/I");
  _mcflash_tree->Branch("tpc", &_tpc, "tpc/I");
  _mcflash_tree->Branch("flash_time", &_flash_time, "flash_time/D");
  _mcflash_tree->Branch("flash_total_pe", &_flash_total_pe, "flash_total_pe/D");
  _mcflash_tree->Branch("flash_y", &_flash_y, "flash_y/D");
  _mcflash_tree->Branch("flash_yerr", &_flash_yerr, "flash_yerr/D");
  _mcflash_tree->Branch("flash_z", &_flash_z, "flash_z/D");
  _mcflash_tree->Branch("flash_zerr", &_flash_zerr, "flash_zerr/D");
  _mcflash_tree->Branch("flash_pe_v", "std::vector<double>", &_flash_pe_v);
}

SBNDFlashAna::SBNDFlashAna ( SBNDFlashAna const &  )

delete

SBNDFlashAna::SBNDFlashAna ( SBNDFlashAna &&  )

delete

Member Function Documentation

void SBNDFlashAna::analyze ( art::Event const &  e )

override

Definition at line 120 of file SBNDFlashAna_module.cc.

{
  // if (_flash_label_v.size() != 2) {
  //   std::cout << "FlashProducer size has to be 2!" << std::endl;
  //   throw std::exception();
  // }
  // if (_mcflash_label_v.size() != 2) {
  //   std::cout << "MCFlashProducer size has to be 2!" << std::endl;
  //   throw std::exception();
  // }


  art::Handle<std::vector<recob::OpHit>> ophit_h;
  e.getByLabel("ophitpmt", ophit_h);
  if(!ophit_h.isValid()){
    std::cout << "Invalid producer for recob::OpHit: " << "ophitpmt" << ". Ignoring." << std::endl;
  }
  std::vector<art::Ptr<recob::OpHit>> ophit_v;
  art::fill_ptr_vector(ophit_v, ophit_h);
  // for (auto oh : ophit_v) {
  //   std::cout << "ophit time " << oh->PeakTime() << std::endl;
  // }



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

  for (size_t l = 0; l < _flash_label_v.size(); l++) {
    art::Handle<std::vector<recob::OpFlash>> flash_h;
    e.getByLabel(_flash_label_v[l], flash_h);
    if(!flash_h.isValid()){
      std::cout << "Invalid producer for recob::OpFlash: " << _flash_label_v[l] << ". Ignoring." << std::endl;
    }

    art::FindManyP<recob::OpHit> flashToOpHitAssns(flash_h, e, _ophit_label_v[l]);

    for (size_t i = 0; i < flash_h->size(); i++) {
      auto const& f = (*flash_h)[i];
      std::cout << "Flash " << i << ", time " << f.AbsTime() << std::endl;
      _tpc = _tpc_v[l];
      _flash_time = f.AbsTime();
      _flash_total_pe = f.TotalPE();
      _flash_pe_v = f.PEs();
      _flash_y = f.YCenter();
      _flash_yerr = f.YWidth();
      _flash_z = f.ZCenter();
      _flash_zerr = f.ZWidth();

      _flash_ophit_time.clear();
      _flash_ophit_amp.clear();
      _flash_ophit_area.clear();
      _flash_ophit_width.clear();
      _flash_ophit_pe.clear();
      _flash_ophit_ch.clear();

      std::vector<art::Ptr<recob::OpHit>> ophit_v = flashToOpHitAssns.at(i);
      for (auto ophit : ophit_v) {
        _flash_ophit_time.push_back(ophit->PeakTimeAbs());
        _flash_ophit_amp.push_back(ophit->Amplitude());
        _flash_ophit_area.push_back(ophit->Area());
        _flash_ophit_width.push_back(ophit->Width());
        _flash_ophit_pe.push_back(ophit->PE());
        _flash_ophit_ch.push_back(ophit->OpChannel());
        _flash_ophit_chname.push_back(_pds_map.pdType(ophit->OpChannel()));
        
      }

      _flash_tree->Fill();
    }
  }

  for (size_t l = 0; l < _mcflash_label_v.size(); l++) {
    art::Handle<std::vector<recob::OpFlash>> flash_h;
    e.getByLabel(_mcflash_label_v[l], flash_h);
    if(!flash_h.isValid()){
      std::cout << "Invalid producer for MC recob::OpFlash: " << _mcflash_label_v[l] << ". Ignoring." << std::endl;
    }
    for (size_t i = 0; i < flash_h->size(); i++) {
      auto const& f = (*flash_h)[i];
      _tpc = _tpc_v[l];
      _flash_time = f.Time();
      _flash_total_pe = f.TotalPE();
      _flash_pe_v = f.PEs();
      _flash_y = f.YCenter();
      _flash_yerr = f.YWidth();
      _flash_z = f.ZCenter();
      _flash_zerr = f.ZWidth();
      _mcflash_tree->Fill();
    }
  }

}

SBNDFlashAna& SBNDFlashAna::operator= ( SBNDFlashAna const &  )

delete

SBNDFlashAna& SBNDFlashAna::operator= ( SBNDFlashAna &&  )

delete

Member Data Documentation

int SBNDFlashAna::_event

private

Definition at line 56 of file SBNDFlashAna_module.cc.

std::vector<std::string> SBNDFlashAna::_flash_label_v

private

Definition at line 48 of file SBNDFlashAna_module.cc.

std::vector<double> SBNDFlashAna::_flash_ophit_amp

private

Definition at line 66 of file SBNDFlashAna_module.cc.

std::vector<double> SBNDFlashAna::_flash_ophit_area

private

Definition at line 67 of file SBNDFlashAna_module.cc.

std::vector<int> SBNDFlashAna::_flash_ophit_ch

private

Definition at line 70 of file SBNDFlashAna_module.cc.

std::vector<std::string> SBNDFlashAna::_flash_ophit_chname

private

Definition at line 71 of file SBNDFlashAna_module.cc.

std::vector<double> SBNDFlashAna::_flash_ophit_pe

private

Definition at line 69 of file SBNDFlashAna_module.cc.

std::vector<double> SBNDFlashAna::_flash_ophit_time

private

Definition at line 65 of file SBNDFlashAna_module.cc.

std::vector<double> SBNDFlashAna::_flash_ophit_width

private

Definition at line 68 of file SBNDFlashAna_module.cc.

std::vector<double> SBNDFlashAna::_flash_pe_v

private

Definition at line 60 of file SBNDFlashAna_module.cc.

double SBNDFlashAna::_flash_time

private

Definition at line 58 of file SBNDFlashAna_module.cc.

double SBNDFlashAna::_flash_total_pe

private

Definition at line 59 of file SBNDFlashAna_module.cc.

TTree* SBNDFlashAna::_flash_tree

private

Definition at line 55 of file SBNDFlashAna_module.cc.

double SBNDFlashAna::_flash_y

private

Definition at line 61 of file SBNDFlashAna_module.cc.

double SBNDFlashAna::_flash_yerr

private

Definition at line 62 of file SBNDFlashAna_module.cc.

double SBNDFlashAna::_flash_z

private

Definition at line 63 of file SBNDFlashAna_module.cc.

double SBNDFlashAna::_flash_zerr

private

Definition at line 64 of file SBNDFlashAna_module.cc.

std::vector<std::string> SBNDFlashAna::_mcflash_label_v

private

Definition at line 50 of file SBNDFlashAna_module.cc.

TTree* SBNDFlashAna::_mcflash_tree

private

Definition at line 73 of file SBNDFlashAna_module.cc.

std::vector<std::string> SBNDFlashAna::_ophit_label_v

private

Definition at line 49 of file SBNDFlashAna_module.cc.

opdet::sbndPDMapAlg SBNDFlashAna::_pds_map

private

Definition at line 53 of file SBNDFlashAna_module.cc.

int SBNDFlashAna::_run

private

Definition at line 56 of file SBNDFlashAna_module.cc.

int SBNDFlashAna::_subrun

private

Definition at line 56 of file SBNDFlashAna_module.cc.

int SBNDFlashAna::_tpc

private

Definition at line 57 of file SBNDFlashAna_module.cc.

std::vector<int> SBNDFlashAna::_tpc_v

private

Definition at line 51 of file SBNDFlashAna_module.cc.

---

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

• SBNDFlashAna_module.cc