pdsp::CheckT0 Class Reference

Inheritance diagram for pdsp::CheckT0:

Public Member Functions
	CheckT0 (fhicl::ParameterSet const &p)
	CheckT0 (CheckT0 const &)=delete
	CheckT0 (CheckT0 &&)=delete
CheckT0 &	operator= (CheckT0 const &)=delete
CheckT0 &	operator= (CheckT0 &&)=delete
void	beginJob () override
void	analyze (art::Event const &e) 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 Attributes
TTree *	t0tree
int	run
int	event
std::vector< int >	trackid
std::vector< double >	t0crt2
std::vector< double >	t0crt1
std::vector< double >	t0pandora
std::vector< double >	t0anodep
std::vector< double >	t0truth
std::vector< double >	trackstartx
std::vector< double >	trackstarty
std::vector< double >	trackstartz
std::vector< double >	trackendx
std::vector< double >	trackendy
std::vector< double >	trackendz
std::vector< double >	trackstartx_sce
std::vector< double >	trackstarty_sce
std::vector< double >	trackstartz_sce
std::vector< double >	trackendx_sce
std::vector< double >	trackendy_sce
std::vector< double >	trackendz_sce
std::vector< double >	crt1x
std::vector< double >	crt1y
std::vector< double >	crt1z
std::vector< double >	crt2x0
std::vector< double >	crt2y0
std::vector< double >	crt2z0
std::vector< double >	crt2x1
std::vector< double >	crt2y1
std::vector< double >	crt2z1
double	rdtimestamp_evt
std::vector< double >	rdtimestamp_digits

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 44 of file CheckT0_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 99 of file CheckT0_module.cc.

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

}

pdsp::CheckT0::CheckT0 ( CheckT0 const &  )

delete

pdsp::CheckT0::CheckT0 ( CheckT0 &&  )

delete

Member Function Documentation

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

overridevirtual

Implements art::EDAnalyzer.

Definition at line 106 of file CheckT0_module.cc.

{
  // Implementation of required member function here.
  run = e.run();
  event = e.id().event();
  trackid.clear();
  t0crt2.clear();
  t0crt1.clear();
  t0pandora.clear();
  t0anodep.clear();
  t0truth.clear();
  trackstartx.clear();
  trackstarty.clear();
  trackstartz.clear();
  trackendx.clear();
  trackendy.clear();
  trackendz.clear();
  trackstartx_sce.clear();
  trackstarty_sce.clear();
  trackstartz_sce.clear();
  trackendx_sce.clear();
  trackendy_sce.clear();
  trackendz_sce.clear();
  crt1x.clear();
  crt1y.clear();
  crt1z.clear();
  crt2x0.clear();
  crt2y0.clear();
  crt2z0.clear();
  crt2x1.clear();
  crt2y1.clear();
  crt2z1.clear();
  rdtimestamp_evt = -1;
  rdtimestamp_digits.clear();

  //Services
  art::ServiceHandle<cheat::BackTrackerService> bt_serv;
  art::ServiceHandle<cheat::ParticleInventoryService> pi_serv;
  //Space charge service
  auto const* SCE = lar::providerFrom<spacecharge::SpaceChargeService>();
  
   // Reconstruciton information
  std::vector < art::Ptr < recob::Track > > trackList;
  auto trackListHandle = e.getHandle< std::vector < recob::Track > >("pandoraTrack");
  if (trackListHandle) {
    art::fill_ptr_vector(trackList, trackListHandle);
  }
  else return;

  //Get hits associated with track
  art::FindManyP < recob::Hit > hitsFromTrack(trackListHandle, e, "pandoraTrack");

  //Get PFParticles
  auto pfpListHandle = e.getHandle< std::vector<recob::PFParticle> >("pandora");
  
  //Get Track-PFParticle association
  art::FindManyP<recob::PFParticle> fmpfp(trackListHandle, e, "pandoraTrack");

  //Get Pandora T0-PFParticle association
  art::FindManyP<anab::T0> fmt0pandora(pfpListHandle, e, "pandora");

  //Get Anode Piercing T0
  art::FindManyP<anab::T0> fmt0anodepiercer(pfpListHandle, e, "anodepiercerst0");
  
  //Get 2-CRT T0
  art::FindManyP<anab::T0> fmt0crt2(trackListHandle, e, "crtreco");
  art::FindManyP<anab::CosmicTag> fmctcrt2(trackListHandle, e, "crtreco");

  //Get 1-CRT T0
  art::FindManyP<anab::T0> fmt0crt1(trackListHandle, e, "crttag");
  art::FindManyP<anab::CosmicTag> fmctcrt1(trackListHandle, e, "crttag");
  
  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService>()->DataFor(e);
  auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService>()->DataFor(e, clockData);

  for (auto const& track : trackList){
    int this_trackid = track.key();
    double this_t0crt2 = -DBL_MAX;
    double this_t0crt1 = -DBL_MAX;
    double this_t0pandora = -DBL_MAX;
    double this_t0anodep = -DBL_MAX;
    double this_t0truth = -DBL_MAX;
    double this_trackstartx = -DBL_MAX;
    double this_trackstarty = -DBL_MAX;
    double this_trackstartz = -DBL_MAX;
    double this_trackendx = -DBL_MAX;
    double this_trackendy = -DBL_MAX;
    double this_trackendz = -DBL_MAX;
    double this_trackstartx_sce = -DBL_MAX;
    double this_trackstarty_sce = -DBL_MAX;
    double this_trackstartz_sce = -DBL_MAX;
    double this_trackendx_sce = -DBL_MAX;
    double this_trackendy_sce = -DBL_MAX;
    double this_trackendz_sce = -DBL_MAX;
    double this_crt1x = -DBL_MAX;
    double this_crt1y = -DBL_MAX;
    double this_crt1z = -DBL_MAX;
    double this_crt2x0 = -DBL_MAX;
    double this_crt2y0 = -DBL_MAX;
    double this_crt2z0 = -DBL_MAX;
    double this_crt2x1 = -DBL_MAX;
    double this_crt2y1 = -DBL_MAX;
    double this_crt2z1 = -DBL_MAX;
    
    if (fmt0crt2.isValid()){
      auto const& vt0crt2 = fmt0crt2.at(track.key());
      if (!vt0crt2.empty()) this_t0crt2 = vt0crt2[0]->Time();
    }

    if (fmt0crt1.isValid()){
      auto const& vt0crt1 = fmt0crt1.at(track.key());
      if (!vt0crt1.empty()) this_t0crt1 = vt0crt1[0]->Time();
    }

    if (fmctcrt2.isValid()){
      auto const& vctcrt2 = fmctcrt2.at(track.key());
      if (!vctcrt2.empty()){
        this_crt2x0 = vctcrt2[0]->EndPoint1()[0];
        this_crt2y0 = vctcrt2[0]->EndPoint1()[1];
        this_crt2z0 = vctcrt2[0]->EndPoint1()[2];
        this_crt2x1 = vctcrt2[0]->EndPoint2()[0];
        this_crt2y1 = vctcrt2[0]->EndPoint2()[1];
        this_crt2z1 = vctcrt2[0]->EndPoint2()[2];
      }
    }

    if (fmctcrt1.isValid()){
      auto const& vctcrt1 = fmctcrt1.at(track.key());
      if (!vctcrt1.empty()){
        this_crt1x = vctcrt1[0]->EndPoint1()[0];
        this_crt1y = vctcrt1[0]->EndPoint1()[1];
        this_crt1z = vctcrt1[0]->EndPoint1()[2];
      }
    }

    if (fmpfp.isValid()){
      auto const &pfps = fmpfp.at(track.key());
      if(!pfps.empty()){
        //Find T0 for PFParticle
        if (fmt0pandora.isValid()){
          auto const &t0s = fmt0pandora.at(pfps[0].key());
          if(!t0s.empty()){
            this_t0pandora = t0s[0]->Time();
          }
        }
        if (fmt0anodepiercer.isValid()){
          auto const &t0aps = fmt0anodepiercer.at(pfps[0].key());
          if (!t0aps.empty()){
            this_t0anodep = t0aps[0]->Time();
          }
        }
      }
    }
    if (this_t0crt2 > -DBL_MAX ||
        this_t0crt1 > -DBL_MAX ||
        this_t0pandora > -DBL_MAX ||
        this_t0anodep > -DBL_MAX){

      auto const & allHits = hitsFromTrack.at(track.key());

      if (!e.isRealData()){
        // Find true particle for reconstructed track
        int TrackID = 0;
        std::map<int,double> trkide;
        for(size_t h = 0; h < allHits.size(); ++h){
          art::Ptr<recob::Hit> hit = allHits[h];
          std::vector<sim::TrackIDE> TrackIDs = bt_serv->HitToTrackIDEs(clockData, hit);
          for(size_t e = 0; e < TrackIDs.size(); ++e){
            trkide[TrackIDs[e].trackID] += TrackIDs[e].energy;
          }         
        }
        // Work out which IDE despoited the most charge in the hit if there was more than one.
        double maxe = -1;
        double tote = 0;
        for (std::map<int,double>::iterator ii = trkide.begin(); ii!=trkide.end(); ++ii){
          tote += ii->second;
          if ((ii->second)>maxe){
            maxe = ii->second;
            TrackID = ii->first;
          }
        }
        // Now have trackID, so get PdG code and T0 etc.
        const simb::MCParticle *particle = pi_serv->TrackIdToParticle_P(TrackID);
        if (particle){
          this_t0truth = particle->T();
        }
      }
      
      this_trackstartx = track->Vertex().X();
      this_trackstarty = track->Vertex().Y();
      this_trackstartz = track->Vertex().Z();
      this_trackendx = track->End().X();
      this_trackendy = track->End().Y();
      this_trackendz = track->End().Z();

      if (std::abs(this_t0pandora+DBL_MAX)<1e-10){
        //no pandora t0 found, correct for t0
        double ticksOffset = 0;
        if (this_t0crt2 > -DBL_MAX) ticksOffset = this_t0crt2/500.+detProp.GetXTicksOffset(allHits[0]->WireID());
        else if (this_t0crt1 > -DBL_MAX) ticksOffset = this_t0crt1/500.+detProp.GetXTicksOffset(allHits[0]->WireID());
        else if (this_t0anodep > -DBL_MAX) ticksOffset = this_t0anodep/500.+detProp.GetXTicksOffset(allHits[0]->WireID());
        double xOffset = detProp.ConvertTicksToX(ticksOffset,allHits[0]->WireID());
        this_trackstartx -= xOffset;
        this_trackendx -= xOffset;
      }

      auto const & posOffsets = SCE->GetCalPosOffsets(geo::Point_t(this_trackstartx, this_trackstarty, this_trackstartz), allHits[0]->WireID().TPC);
      this_trackstartx_sce = this_trackstartx - posOffsets.X();
      this_trackstarty_sce = this_trackstarty + posOffsets.Y();
      this_trackstartz_sce = this_trackstartz + posOffsets.Z();
      this_trackendx_sce = this_trackendx - posOffsets.X();
      this_trackendy_sce = this_trackendy + posOffsets.Y();
      this_trackendz_sce = this_trackendz + posOffsets.Z();

      trackid.push_back(this_trackid);
      t0crt2.push_back(this_t0crt2);
      t0crt1.push_back(this_t0crt1);
      t0pandora.push_back(this_t0pandora);
      t0anodep.push_back(this_t0anodep);
      t0truth.push_back(this_t0truth);
      trackstartx.push_back(this_trackstartx);
      trackstarty.push_back(this_trackstarty);
      trackstartz.push_back(this_trackstartz);
      trackendx.push_back(this_trackendx);
      trackendy.push_back(this_trackendy);
      trackendz.push_back(this_trackendz);
      trackstartx_sce.push_back(this_trackstartx_sce);
      trackstarty_sce.push_back(this_trackstarty_sce);
      trackstartz_sce.push_back(this_trackstartz_sce);
      trackendx_sce.push_back(this_trackendx_sce);
      trackendy_sce.push_back(this_trackendy_sce);
      trackendz_sce.push_back(this_trackendz_sce);
      crt1x.push_back(this_crt1x);
      crt1y.push_back(this_crt1y);
      crt1z.push_back(this_crt1z);
      crt2x0.push_back(this_crt2x0);
      crt2y0.push_back(this_crt2y0);
      crt2z0.push_back(this_crt2z0);
      crt2x1.push_back(this_crt2x1);
      crt2y1.push_back(this_crt2y1);
      crt2z1.push_back(this_crt2z1);
    }
  }

  //RDTimeStamps for the event
  std::vector < art::Ptr < raw::RDTimeStamp > > rdtsevtList;
  art::InputTag itag("timingrawdecoder","daq");
  auto rdts_evt = e.getHandle< std::vector < raw::RDTimeStamp > >(itag); 
  if (rdts_evt) {
    art::fill_ptr_vector(rdtsevtList, rdts_evt);
  }
  
  if (!rdtsevtList.empty()){
    rdtimestamp_evt = rdtsevtList[0]->GetTimeStamp();
  }
  
  //RDTimeStamps for each raw digit (channel)
  art::InputTag itag2("tpcrawdecoder","daq");
  std::vector < art::Ptr < raw::RDTimeStamp > > rdtsdigitList;
  auto rdts_digit = e.getHandle < std::vector < raw::RDTimeStamp > >(itag2);
  if (rdts_digit) {
    art::fill_ptr_vector(rdtsdigitList, rdts_digit);
  }
  
  for (const auto & rdts : rdtsdigitList){
    rdtimestamp_digits.push_back(rdts->GetTimeStamp());
  }

  if (!trackid.empty()) t0tree->Fill();
}

void pdsp::CheckT0::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 377 of file CheckT0_module.cc.

                           {
  art::ServiceHandle<art::TFileService> fileServiceHandle;
  t0tree = fileServiceHandle->make<TTree>("t0tree", "t0 info");
  t0tree->Branch("run", &run, "run/I");
  t0tree->Branch("event", &event, "event/I");
  t0tree->Branch("trackid", &trackid);
  t0tree->Branch("t0crt2", &t0crt2);
  t0tree->Branch("t0crt1", &t0crt1);
  t0tree->Branch("t0pandora",&t0pandora);
  t0tree->Branch("t0anodep",&t0anodep);
  t0tree->Branch("t0truth", &t0truth);
  t0tree->Branch("trackstartx",&trackstartx);
  t0tree->Branch("trackstarty",&trackstarty);
  t0tree->Branch("trackstartz",&trackstartz);
  t0tree->Branch("trackendx",&trackendx);
  t0tree->Branch("trackendy",&trackendy);
  t0tree->Branch("trackendz",&trackendz);
  t0tree->Branch("trackstartx_sce",&trackstartx_sce);
  t0tree->Branch("trackstarty_sce",&trackstarty_sce);
  t0tree->Branch("trackstartz_sce",&trackstartz_sce);
  t0tree->Branch("trackendx_sce",&trackendx_sce);
  t0tree->Branch("trackendy_sce",&trackendy_sce);
  t0tree->Branch("trackendz_sce",&trackendz_sce);
  t0tree->Branch("crt1x",&crt1x);
  t0tree->Branch("crt1y",&crt1y);
  t0tree->Branch("crt1z",&crt1z);
  t0tree->Branch("crt2x0",&crt2x0);
  t0tree->Branch("crt2y0",&crt2y0);
  t0tree->Branch("crt2z0",&crt2z0);
  t0tree->Branch("crt2x1",&crt2x1);
  t0tree->Branch("crt2y1",&crt2y1);
  t0tree->Branch("crt2z1",&crt2z1);
  t0tree->Branch("rdtimestamp_evt", &rdtimestamp_evt, "rdtimestamp_evt/D");
  t0tree->Branch("rdtimestamp_digits", &rdtimestamp_digits);
}

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

delete

CheckT0& pdsp::CheckT0::operator= ( CheckT0 &&  )

delete

Member Data Documentation

std::vector<double> pdsp::CheckT0::crt1x

private

Definition at line 84 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::crt1y

private

Definition at line 85 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::crt1z

private

Definition at line 86 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::crt2x0

private

Definition at line 87 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::crt2x1

private

Definition at line 90 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::crt2y0

private

Definition at line 88 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::crt2y1

private

Definition at line 91 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::crt2z0

private

Definition at line 89 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::crt2z1

private

Definition at line 92 of file CheckT0_module.cc.

int pdsp::CheckT0::event

private

Definition at line 65 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::rdtimestamp_digits

private

Definition at line 94 of file CheckT0_module.cc.

double pdsp::CheckT0::rdtimestamp_evt

private

Definition at line 93 of file CheckT0_module.cc.

int pdsp::CheckT0::run

private

Definition at line 64 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::t0anodep

private

Definition at line 70 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::t0crt1

private

Definition at line 68 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::t0crt2

private

Definition at line 67 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::t0pandora

private

Definition at line 69 of file CheckT0_module.cc.

TTree* pdsp::CheckT0::t0tree

private

Definition at line 63 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::t0truth

private

Definition at line 71 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackendx

private

Definition at line 75 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackendx_sce

private

Definition at line 81 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackendy

private

Definition at line 76 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackendy_sce

private

Definition at line 82 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackendz

private

Definition at line 77 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackendz_sce

private

Definition at line 83 of file CheckT0_module.cc.

std::vector<int> pdsp::CheckT0::trackid

private

Definition at line 66 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackstartx

private

Definition at line 72 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackstartx_sce

private

Definition at line 78 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackstarty

private

Definition at line 73 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackstarty_sce

private

Definition at line 79 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackstartz

private

Definition at line 74 of file CheckT0_module.cc.

std::vector<double> pdsp::CheckT0::trackstartz_sce

private

Definition at line 80 of file CheckT0_module.cc.

---

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

• protoduneana/protoduneana/singlephase/T0Ana/CheckT0_module.cc