pdsp::CheckCRT Class Reference

Inheritance diagram for pdsp::CheckCRT:

Public Member Functions
	CheckCRT (fhicl::ParameterSet const &p)
	CheckCRT (CheckCRT const &)=delete
	CheckCRT (CheckCRT &&)=delete
CheckCRT &	operator= (CheckCRT const &)=delete
CheckCRT &	operator= (CheckCRT &&)=delete
void	analyze (art::Event const &e) override
void	beginJob () override
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
fhicl::ParameterSetID	selectorConfig () const
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Member Functions
std::vector< unsigned short >	GetModules (int pixel)

Private Attributes
TTree *	fTree
int	run
int	subrun
int	event
std::vector< uint32_t >	crtmask
std::vector< ULong64_t >	rdtime
std::vector< unsigned int >	crtmodule
std::vector< unsigned long long >	crttime
std::vector< double >	trkstartx_pandora
std::vector< double >	trkstarty_pandora
std::vector< double >	trkstartz_pandora
std::vector< double >	trkendx_pandora
std::vector< double >	trkendy_pandora
std::vector< double >	trkendz_pandora
std::vector< double >	t0_pandora
std::vector< double >	trkstartx_pmtrack
std::vector< double >	trkstarty_pmtrack
std::vector< double >	trkstartz_pmtrack
std::vector< double >	trkendx_pmtrack
std::vector< double >	trkendy_pmtrack
std::vector< double >	trkendz_pmtrack
std::vector< double >	t0_pmtrack

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 34 of file CheckCRT_module.cc.

Constructor & Destructor Documentation

pdsp::CheckCRT::CheckCRT ( fhicl::ParameterSet const &  p )

explicit

Definition at line 84 of file CheckCRT_module.cc.

  : EDAnalyzer{p}  // ,
  // More initializers here.
{
}

pdsp::CheckCRT::CheckCRT ( CheckCRT const &  )

delete

pdsp::CheckCRT::CheckCRT ( CheckCRT &&  )

delete

Member Function Documentation

void pdsp::CheckCRT::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 251 of file CheckCRT_module.cc.

{

  run = e.run();
  subrun = e.subRun();
  event = e.id().event();

  crtmask.clear();
  rdtime.clear();
  crtmodule.clear();
  crttime.clear();

  trkstartx_pandora.clear();
  trkstarty_pandora.clear();
  trkstartz_pandora.clear();
  trkendx_pandora.clear();
  trkendy_pandora.clear();
  trkendz_pandora.clear();
  t0_pandora.clear();

  trkstartx_pmtrack.clear();
  trkstarty_pmtrack.clear();
  trkstartz_pmtrack.clear();
  trkendx_pmtrack.clear();
  trkendy_pmtrack.clear();
  trkendz_pmtrack.clear();
  t0_pmtrack.clear();

  art::ValidHandle<std::vector<raw::RDTimeStamp>> timeStamps = e.getValidHandle<std::vector<raw::RDTimeStamp>>("timingrawdecoder:daq");

  int trigger = -1;
  // Check that we have good information
  if(timeStamps.isValid() && timeStamps->size() == 1){
    // Access the trigger information. Beam trigger flag = 0xc
    const raw::RDTimeStamp& timeStamp = timeStamps->at(0);
    trigger = timeStamp.GetFlags();
  }

  if (trigger!=13) return;

  auto const& pdspctbs = *e.getValidHandle<std::vector<raw::ctb::pdspctb>>("ctbrawdecoder:daq");
  if (!pdspctbs.empty()){
    const size_t npdspctbs = pdspctbs.size();
    for (size_t i=0; i<npdspctbs; ++i){
//      std::cout << "Number of triggers: " << pdspctbs[i].GetNTriggers() << std::endl;
//      std::cout << "Number of chstatuses: " << pdspctbs[i].GetNChStatuses() << std::endl;
//      std::cout << "Number of feedbacks: " << pdspctbs[i].GetNFeedbacks() << std::endl;
//      std::cout << "Number of miscs: " << pdspctbs[i].GetNMiscs() << std::endl;
//      
//      const std::vector<raw::ctb::Trigger> HLTriggers = pdspctbs[i].GetHLTriggers();
//      std::cout << "Number of HL Triggers: " << HLTriggers.size() << std::endl;
//      for (size_t j=0; j<HLTriggers.size(); ++j){
//        std::cout << "HL Trigger: " << j << " " << HLTriggers[j].word_type << " " << HLTriggers[j].trigger_word << " " << HLTriggers[j].timestamp << std::endl;
//      }
      
      const std::vector<raw::ctb::ChStatus> chs = pdspctbs[i].GetChStatusAfterHLTs();
      std::cout << "Number of CH Status after HLTs: " << chs.size() << std::endl;
      for (size_t j=0; j<chs.size(); ++j){
        //std::cout << " chs after HLT: " << j << " " << chs[j].word_type << " " << chs[j].pds << " " << chs[j].crt << " " << chs[j].beam_hi << " " << chs[j].beam_lo << std::endl;
        crtmask.push_back(chs[j].crt);
      }
//      
//      
//      const std::vector<raw::ctb::Trigger> LLTriggers = pdspctbs[i].GetLLTriggers();
//      std::cout << "Number of LL Triggers: " << LLTriggers.size() << std::endl;
//      for (size_t j=0; j<LLTriggers.size(); ++j){
//        std::cout << "LL Trigger: " << j << " " << LLTriggers[j].word_type << " " << LLTriggers[j].trigger_word << " " << LLTriggers[j].timestamp << std::endl;
//      }               
//      
//      const std::vector<raw::ctb::WordIndex> indexes = pdspctbs[i].GetIndexes();
//      std::cout << "Number of Indexes: " << indexes.size() << std::endl;
//      for (size_t j=0; j<indexes.size(); ++j){
//        std::cout << "Index: " << j << " " << indexes[j].word_type << " " << indexes[j].index << std::endl;
//      }               
      
    }
  }
  //Save CRT trigger information
  const auto & triggers = *e.getValidHandle < std::vector < CRT::Trigger >> ("crtrawdecoder");
  for (const auto& trigger: triggers){
    //std::cout<<"module: "<<trigger.Channel()<<" time: "<<trigger.Timestamp()<<std::endl;
    crtmodule.push_back(trigger.Channel());
    crttime.push_back(trigger.Timestamp());
  }
  for (const auto& time : *timeStamps){
    rdtime.push_back(time.GetTimeStamp());
  }
  /*
  for (auto const& crt_mask : crtmask){
    int pixel[2] = {-1, -1};
    for (size_t i = 0; i<32; ++i){
      if (crt_mask & (1<<i)){
        if (i<16){
          pixel[0] = i;
        }
        else {
          pixel[1] = i;
        }
      }
    }
    std::cout<<"Pixels "<<pixel[0]<<" "<<pixel[1]<<std::endl;
    for (size_t i = 0; i<2; ++i){
      if (pixel[i]!=-1){
        std::vector<unsigned short> modules = GetModules(pixel[i]);
        for (size_t j = 0; j<2; ++j){
          for (const auto& trigger: triggers){
            if (trigger.Channel() == modules[j]){
              std::cout<<"Pixel: "<<pixel[i]<<" module: "<<trigger.Channel()<<" time: "<<trigger.Timestamp()<<" "<<timeStamps->at(0).GetTimeStamp()<<std::endl;
            }
          }
        }
      }
    }
  }
  */

  std::vector< art::Ptr<recob::Track> > pandoratrks;
  auto pandoratrkHandle = e.getHandle< std::vector<recob::Track> >("pandoraTrack");
  if (pandoratrkHandle) art::fill_ptr_vector(pandoratrks, pandoratrkHandle);
  else return;
  
  auto pfpListHandle = e.getHandle< std::vector<recob::PFParticle> >("pandora");
  art::FindManyP<recob::PFParticle> fmpfp(pandoratrkHandle, e, "pandoraTrack");
  art::FindManyP<anab::T0> fmt0pandora(pfpListHandle, e, "pandora");

  for (size_t i = 0; i<pandoratrks.size(); ++i){
    auto & trk = pandoratrks[i];
    trkstartx_pandora.push_back(trk->Vertex().X());
    trkstarty_pandora.push_back(trk->Vertex().Y());
    trkstartz_pandora.push_back(trk->Vertex().Z());
    trkendx_pandora.push_back(trk->End().X());
    trkendy_pandora.push_back(trk->End().Y());
    trkendz_pandora.push_back(trk->End().Z());
    double t0 = 0;
    auto &pfps = fmpfp.at(trk.key());
    if (!pfps.empty()){
      auto &t0s = fmt0pandora.at(pfps[0].key());
      if (!t0s.empty()){
        t0 = t0s[0]->Time();
      }
    }
    t0_pandora.push_back(t0);
  }

  /*
  std::vector< art::Ptr<recob::Track> > pmtracktrks;
  auto pmtracktrkHandle = e.getHandle< std::vector<recob::Track> >("pmtrack");
  if (pmtracktrkHandle) art::fill_ptr_vector(pmtracktrks, pmtracktrkHandle);
  
  art::FindManyP<anab::T0> fmt0pmtrack(pmtracktrkHandle, e, "pmtrack");

  for (size_t i = 0; i<pmtracktrks.size(); ++i){
    auto & trk = pmtracktrks[i];
    trkstartx_pmtrack.push_back(trk->Vertex().X());
    trkstarty_pmtrack.push_back(trk->Vertex().Y());
    trkstartz_pmtrack.push_back(trk->Vertex().Z());
    trkendx_pmtrack.push_back(trk->End().X());
    trkendy_pmtrack.push_back(trk->End().Y());
    trkendz_pmtrack.push_back(trk->End().Z());
    double t0 = 0;
    auto &t0s = fmt0pmtrack.at(trk.key());
    if (!t0s.empty()){
      t0 = t0s[0]->Time();
    }
    t0_pmtrack.push_back(t0);
  }
  */
  fTree->Fill();
}

void pdsp::CheckCRT::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 222 of file CheckCRT_module.cc.

{
  art::ServiceHandle<art::TFileService> tfs;
  fTree = tfs->make<TTree>("crt","crt tree");
  fTree->Branch("run",&run,"run/I");
  fTree->Branch("subrun",&subrun,"subrun/I");
  fTree->Branch("event",&event,"event/I");
  fTree->Branch("crtmask", &crtmask);
  fTree->Branch("rdtime", &rdtime);
  fTree->Branch("crtmodule", &crtmodule);
  fTree->Branch("crttime", &crttime);
  fTree->Branch("trkstartx_pandora",&trkstartx_pandora);
  fTree->Branch("trkstarty_pandora",&trkstarty_pandora);
  fTree->Branch("trkstartz_pandora",&trkstartz_pandora);
  fTree->Branch("trkendx_pandora",&trkendx_pandora);
  fTree->Branch("trkendy_pandora",&trkendy_pandora);
  fTree->Branch("trkendz_pandora",&trkendz_pandora);
  fTree->Branch("t0_pandora",&t0_pandora);

  fTree->Branch("trkstartx_pmtrack",&trkstartx_pmtrack);
  fTree->Branch("trkstarty_pmtrack",&trkstarty_pmtrack);
  fTree->Branch("trkstartz_pmtrack",&trkstartz_pmtrack);
  fTree->Branch("trkendx_pmtrack",&trkendx_pmtrack);
  fTree->Branch("trkendy_pmtrack",&trkendy_pmtrack);
  fTree->Branch("trkendz_pmtrack",&trkendz_pmtrack);
  fTree->Branch("t0_pmtrack",&t0_pmtrack);

}

std::vector< unsigned short > pdsp::CheckCRT::GetModules ( int  pixel )

private

Definition at line 89 of file CheckCRT_module.cc.

                                                           {
  std::vector<unsigned short> modules(2,-1);
  if (pixel == 0){
    modules[0] = 1;
    modules[1] = 2;
  }
  else if (pixel == 1){
    modules[0] = 0;
    modules[1] = 2;
  }
  else if (pixel == 2){
    modules[0] = 15;
    modules[1] = 13;
  }
  else if (pixel == 3){
    modules[0] = 14;
    modules[1] = 13;
  }
  else if (pixel == 4){
    modules[0] = 14;
    modules[1] = 12;
  }
  else if (pixel == 5){
    modules[0] = 9;
    modules[1] = 11;
  }
  else if (pixel == 6){
    modules[0] = 9;
    modules[1] = 10;
  }
  else if (pixel == 7){
    modules[0] = 8;
    modules[1] = 10;
  }
  else if (pixel == 8){
    modules[0] = 7;
    modules[1] = 5;
  }
  else if (pixel == 9){
    modules[0] = 6;
    modules[1] = 5;
  }
  else if (pixel == 10){
    modules[0] = 6;
    modules[1] = 4;
  }
  else if (pixel == 11){
    modules[0] = 1;
    modules[1] = 3;
  }
  else if (pixel == 12){
    modules[0] = 0;
    modules[1] = 3;
  }
  else if (pixel == 13){
    modules[0] = 15;
    modules[1] = 12;
  }
  else if (pixel == 14){
    modules[0] = 8;
    modules[1] = 11;
  }
  else if (pixel == 15){
    modules[0] = 7;
    modules[1] = 4;
  }
  else if (pixel == 16){
    modules[0] = 17;
    modules[1] = 18;
  }
  else if (pixel == 17){
    modules[0] = 16;
    modules[1] = 18;
  }
  else if (pixel == 18){
    modules[0] = 31;
    modules[1] = 29;
  }
  else if (pixel == 19){
    modules[0] = 30;
    modules[1] = 29;
  }
  else if (pixel == 20){
    modules[0] = 30;
    modules[1] = 28;
  }
  else if (pixel == 21){
    modules[0] = 25;
    modules[1] = 27;
  }
  else if (pixel == 22){
    modules[0] = 25;
    modules[1] = 26;
  }
  else if (pixel == 23){
    modules[0] = 24;
    modules[1] = 26;
  }
  else if (pixel == 24){
    modules[0] = 23;
    modules[1] = 21;
  }
  else if (pixel == 25){
    modules[0] = 22;
    modules[1] = 21;
  }
  else if (pixel == 26){
    modules[0] = 22;
    modules[1] = 20;
  }
  else if (pixel == 27){
    modules[0] = 17;
    modules[1] = 19;
  }
  else if (pixel == 28){
    modules[0] = 16;
    modules[1] = 19;
  }
  else if (pixel == 29){
    modules[0] = 31;
    modules[1] = 28;
  }
  else if (pixel == 30){
    modules[0] = 24;
    modules[1] = 27;
  }
  else if (pixel == 31){
    modules[0] = 23;
    modules[1] = 20;
  }
  return modules;
}

CheckCRT& pdsp::CheckCRT::operator= ( CheckCRT const &  )

delete

CheckCRT& pdsp::CheckCRT::operator= ( CheckCRT &&  )

delete

Member Data Documentation

std::vector<uint32_t> pdsp::CheckCRT::crtmask

private

Definition at line 61 of file CheckCRT_module.cc.

std::vector<unsigned int> pdsp::CheckCRT::crtmodule

private

Definition at line 63 of file CheckCRT_module.cc.

std::vector<unsigned long long> pdsp::CheckCRT::crttime

private

Definition at line 64 of file CheckCRT_module.cc.

int pdsp::CheckCRT::event

private

Definition at line 60 of file CheckCRT_module.cc.

TTree* pdsp::CheckCRT::fTree

private

Definition at line 56 of file CheckCRT_module.cc.

std::vector<ULong64_t> pdsp::CheckCRT::rdtime

private

Definition at line 62 of file CheckCRT_module.cc.

int pdsp::CheckCRT::run

private

Definition at line 58 of file CheckCRT_module.cc.

int pdsp::CheckCRT::subrun

private

Definition at line 59 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::t0_pandora

private

Definition at line 71 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::t0_pmtrack

private

Definition at line 79 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkendx_pandora

private

Definition at line 68 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkendx_pmtrack

private

Definition at line 76 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkendy_pandora

private

Definition at line 69 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkendy_pmtrack

private

Definition at line 77 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkendz_pandora

private

Definition at line 70 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkendz_pmtrack

private

Definition at line 78 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkstartx_pandora

private

Definition at line 65 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkstartx_pmtrack

private

Definition at line 73 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkstarty_pandora

private

Definition at line 66 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkstarty_pmtrack

private

Definition at line 74 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkstartz_pandora

private

Definition at line 67 of file CheckCRT_module.cc.

std::vector<double> pdsp::CheckCRT::trkstartz_pmtrack

private

Definition at line 75 of file CheckCRT_module.cc.

---

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

• protoduneana/protoduneana/singlephase/T0Ana/CheckCRT_module.cc