DAQToOffline::BadTimings Class Reference

Inheritance diagram for DAQToOffline::BadTimings:

Public Member Functions
	BadTimings (const fhicl::ParameterSet &pset)
void	analyze (const art::Event &evt)
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 ()
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 &)
std::string const &	processName () const
bool	wantAllEvents () const
bool	wantEvent (Event const &e)
fhicl::ParameterSetID	selectorConfig () const
Handle< TriggerResults >	getTriggerResults (Event const &e) 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
std::string	fTPCLabel
std::string	fTPCInstance
std::string	fPTBLabel
std::string	fPTBInstance
lbne::TpcNanoSlice::Header::nova_timestamp_t	fRCETriggerTimestamp
lbne::PennMilliSlice::Header::timestamp_t	fPTBTriggerTimestamp
art::ServiceHandle< geo::Geometry >	fGeometry
art::ServiceHandle< art::TFileService >	fTFS
std::string	fPTBMapFile
std::map< int, int >	fPTBMap
std::vector< bool >	fWithinTrigger
TTree *	fTree
int	fRun
int	fEvent
int	fRCE
int	fDiffTriggerTimestamps
int	fMicrosliceWithTrigger
TH1I *	hDiffTriggerTimestamps
TH1I *	hTriggerStart

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &paths, fhicl::ParameterSet const &config)
detail::ProcessAndEventSelectors &	processAndEventSelectors ()
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 48 of file BadTimings_module.cc.

Constructor & Destructor Documentation

DAQToOffline::BadTimings::BadTimings ( const fhicl::ParameterSet &  pset )

explicit

Definition at line 84 of file BadTimings_module.cc.

                                                                : EDAnalyzer(pset) {
  fTPCLabel = pset.get<std::string>("TPCLabel");
  fTPCInstance = pset.get<std::string>("TPCInstance");
  fPTBLabel = pset.get<std::string>("PTBLabel");
  fPTBInstance = pset.get<std::string>("PTBInstance");
  fPTBMapFile = pset.get<std::string>("PTBMapFile");

  fWithinTrigger = std::vector<bool>(16,false);

  hDiffTriggerTimestamps = fTFS->make<TH1I>("PTBRCEDiffTimestamps",";RCE Trigger Timestamp - PTB Trigger Timestamp (NOvA ticks);",100,0,2000);
  hTriggerStart = fTFS->make<TH1I>("TriggerStart",";Microslice number containing trigger;",15,0,15);
  hTriggerStart->GetXaxis()->SetNdivisions(15);
  hTriggerStart->GetXaxis()->CenterLabels();

  fTree = fTFS->make<TTree>("BadTimings","BadTimings");
  fTree->Branch("Run",                  &fRun);
  fTree->Branch("Event",                &fEvent);
  fTree->Branch("RCE",                  &fRCE);
  // fTree->Branch("RCETriggerTimestamp",  &fRCETriggerTimestamp);
  // fTree->Branch("PTBTriggerTimestamp",  &fPTBTriggerTimestamp);
  fTree->Branch("DiffTriggerTimestamps",&fDiffTriggerTimestamps);
  fTree->Branch("MicrosliceWithTrigger",&fMicrosliceWithTrigger);

  DAQToOffline::BuildPTBChannelMap("", fPTBMapFile, fPTBMap);
}

Member Function Documentation

void DAQToOffline::BadTimings::analyze

(

const art::Event &

evt

)

Definition at line 110 of file BadTimings_module.cc.

                                                        {

  fRun = evt.run();
  fEvent = evt.event();

  // Get the TPC data out of the event
  art::Handle<artdaq::Fragments> RCERawFragments;
  evt.getByLabel(fTPCLabel, fTPCInstance, RCERawFragments);

  // Get the counter data out of the event
  art::Handle<artdaq::Fragments> PTBRawFragments;
  evt.getByLabel(fPTBLabel, fPTBInstance, PTBRawFragments);

  // Look at PTB
  for (size_t fragIndex = 0; fragIndex < PTBRawFragments->size(); ++fragIndex) {
    const artdaq::Fragment &singleFragment = (*PTBRawFragments)[fragIndex];
    lbne::PennMilliSliceFragment msf(singleFragment);

    lbne::PennMicroSlice::Payload_Header *word_header = nullptr;
    //lbne::PennMicroSlice::Payload_Counter *word_p_counter = nullptr;
    lbne::PennMicroSlice::Payload_Trigger *word_p_trigger = nullptr;
    lbne::PennMicroSlice::Payload_Timestamp *previous_timestamp = nullptr;
    lbne::PennMicroSlice::Payload_Header *future_timestamp_header = nullptr;
    lbne::PennMicroSlice::Payload_Timestamp *future_timestamp = nullptr;
    uint8_t* payload_data = nullptr;
    uint32_t payload_index = 0;

    uint16_t counter, trigger, timestamp, payloadCount;
    payloadCount = msf.payloadCount(counter, trigger, timestamp);

    if (trigger != 1) {
      //std::cout << "Event has no trigger in PTB data" << std::endl;
      return;
    }

    while (payload_index < uint32_t(payloadCount-1)) {

      payload_data = msf.get_next_payload(payload_index,word_header);
      if (payload_data == nullptr)
        continue;

      if (word_header->data_packet_type == lbne::PennMicroSlice::DataTypeTimestamp)
        previous_timestamp = reinterpret_cast<lbne::PennMicroSlice::Payload_Timestamp*>(payload_data);
      if (word_header->data_packet_type != lbne::PennMicroSlice::DataTypeTrigger)
        continue;

      word_p_trigger = reinterpret_cast<lbne::PennMicroSlice::Payload_Trigger*>(payload_data);
      if (!word_p_trigger->has_muon_trigger()) {
        std::cout << "Trigger is not from muon" << std::endl;
        return;
      }

      if (previous_timestamp != nullptr)
        fPTBTriggerTimestamp = word_header->get_full_timestamp_post(previous_timestamp->nova_timestamp);
      else {
        if (future_timestamp == nullptr) {
          future_timestamp = reinterpret_cast<lbne::PennMicroSlice::Payload_Timestamp*>(msf.get_next_timestamp(future_timestamp_header));
          if (future_timestamp == nullptr) {
            std::cerr << "CAN'T FIND PTB TIMESTAMP WORDS IN MILLISLICE FRAGMENT!!! Logic will fail." << std::endl;
            return;
          }
        }
        fPTBTriggerTimestamp = word_header->get_full_timestamp_pre(future_timestamp->nova_timestamp);
      }

    } // loop over payloads

  } // fragments

  // Look at RCE
  // Create a map containing (fragmentID, fragIndex) for the event, will be used to check if each channel is present
  std::map<unsigned int,unsigned int> mapFragID;
  for (size_t fragIndex = 0; fragIndex < RCERawFragments->size(); ++fragIndex) {
    const artdaq::Fragment &singleFragment = (*RCERawFragments)[fragIndex];
    unsigned int fragmentID = singleFragment.fragmentID();
    mapFragID.insert(std::pair<unsigned int, unsigned int>(fragmentID,fragIndex));
  }

  int maxRCEs = 16;

  // Loop over the RCEs
  for (std::map<unsigned int,unsigned int>::iterator rceIt = mapFragID.begin(); rceIt != mapFragID.end() and std::distance(mapFragID.begin(),rceIt) < maxRCEs; ++rceIt) {

    fRCE = rceIt->first - 100;

    // Get millislice
    const artdaq::Fragment &singleFragment = (*RCERawFragments)[rceIt->second];
    lbne::TpcMilliSliceFragment millisliceFragment(singleFragment);
    //uint64_t millisliceTimestamp = millisliceFragment.microSlice(0)->software_message() & (((long int)0x1 << 56) - 1);

    // Get microslices
    auto numMicroSlices = millisliceFragment.microSliceCount();
    for (unsigned int i_micro = 0; i_micro < numMicroSlices; i_micro++) {

      std::unique_ptr<const lbne::TpcMicroSlice> microSlice = millisliceFragment.microSlice(i_micro);
      uint64_t rceTriggerTimestamp = microSlice->software_message() & (((long int)0x1 << 56) - 1);

      // Get nanoslices
      auto numNanoSlices = microSlice->nanoSliceCount();

      if (!fWithinTrigger[fRCE] and numNanoSlices) {

        // New trigger
        fWithinTrigger[fRCE] = true;
        std::cout << std::endl << "New trigger (RCE " << fRCE << ") on microslice " << i_micro << " (event " << evt.event() << ")" << std::endl;
        //uint64_t nanoSliceTimestamp = microSlice->nanoSlice(0)->nova_timestamp();

        fRCETriggerTimestamp = rceTriggerTimestamp;
        fMicrosliceWithTrigger = (rceTriggerTimestamp - microSlice->nanoSlice(0)->nova_timestamp())/*time*/ * (15.625/500)/*nova->tpc tick*/ * 1/1000/* nano->micro*/;
        fDiffTriggerTimestamps = fRCETriggerTimestamp - fPTBTriggerTimestamp;

        hDiffTriggerTimestamps->Fill(fDiffTriggerTimestamps);
        hTriggerStart->Fill(fMicrosliceWithTrigger);

        fTree->Fill();

      }
      else if (fWithinTrigger[fRCE] and !numNanoSlices) {

        // End of trigger
        fWithinTrigger[fRCE] = false;
        std::cout << "End of trigger (RCE " << fRCE << ") on microslice " << i_micro << " (event " << evt.event() << ")" << std::endl;

      }

    } // microslices

  } // rces

  return;

}

Member Data Documentation

int DAQToOffline::BadTimings::fDiffTriggerTimestamps

private

Definition at line 76 of file BadTimings_module.cc.

int DAQToOffline::BadTimings::fEvent

private

Definition at line 75 of file BadTimings_module.cc.

art::ServiceHandle<geo::Geometry> DAQToOffline::BadTimings::fGeometry

private

Definition at line 66 of file BadTimings_module.cc.

int DAQToOffline::BadTimings::fMicrosliceWithTrigger

private

Definition at line 77 of file BadTimings_module.cc.

std::string DAQToOffline::BadTimings::fPTBInstance

private

Definition at line 59 of file BadTimings_module.cc.

std::string DAQToOffline::BadTimings::fPTBLabel

private

Definition at line 59 of file BadTimings_module.cc.

std::map<int,int> DAQToOffline::BadTimings::fPTBMap

private

Definition at line 70 of file BadTimings_module.cc.

std::string DAQToOffline::BadTimings::fPTBMapFile

private

Definition at line 69 of file BadTimings_module.cc.

lbne::PennMilliSlice::Header::timestamp_t DAQToOffline::BadTimings::fPTBTriggerTimestamp

private

Definition at line 62 of file BadTimings_module.cc.

int DAQToOffline::BadTimings::fRCE

private

Definition at line 75 of file BadTimings_module.cc.

lbne::TpcNanoSlice::Header::nova_timestamp_t DAQToOffline::BadTimings::fRCETriggerTimestamp

private

Definition at line 61 of file BadTimings_module.cc.

int DAQToOffline::BadTimings::fRun

private

Definition at line 75 of file BadTimings_module.cc.

art::ServiceHandle<art::TFileService> DAQToOffline::BadTimings::fTFS

private

Definition at line 67 of file BadTimings_module.cc.

std::string DAQToOffline::BadTimings::fTPCInstance

private

Definition at line 58 of file BadTimings_module.cc.

std::string DAQToOffline::BadTimings::fTPCLabel

private

Definition at line 58 of file BadTimings_module.cc.

TTree* DAQToOffline::BadTimings::fTree

private

Definition at line 74 of file BadTimings_module.cc.

std::vector<bool> DAQToOffline::BadTimings::fWithinTrigger

private

Definition at line 72 of file BadTimings_module.cc.

TH1I* DAQToOffline::BadTimings::hDiffTriggerTimestamps

private

Definition at line 79 of file BadTimings_module.cc.

TH1I* DAQToOffline::BadTimings::hTriggerStart

private

Definition at line 80 of file BadTimings_module.cc.

---

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

• dunetpc/dune/daqinput35t/BadTimings_module.cc