protoana::ProtoDUNEDataUtils Class Reference

#include <ProtoDUNEDataUtils.h>

Public Member Functions
	ProtoDUNEDataUtils (fhicl::ParameterSet const &pset)
	~ProtoDUNEDataUtils ()
void	reconfigure (fhicl::ParameterSet const &pset)
bool	IsBeamTrigger (art::Event const &evt) const
int	GetNActiveFembsForAPA (art::Event const &evt, int apa) const
	Get number of active fembs in an APA. More...
bool	CheckTimeStampConsistencyForAPAs (art::Event const &evt, std::set< int > apas, ULong64_t &timestamp, ULong64_t &timestamp2, int &apainconsist) const
	Check for consistency of timestamp values for a set of APAs. True if consistent, false if there are mismatches. More...

Private Attributes
art::InputTag	fTimingTag
art::InputTag	fRawDigitTag
art::InputTag	fRawDigitTimeStampTag

Detailed Description

Definition at line 21 of file ProtoDUNEDataUtils.h.

Constructor & Destructor Documentation

protoana::ProtoDUNEDataUtils::ProtoDUNEDataUtils

(

fhicl::ParameterSet const &

pset

)

Definition at line 12 of file ProtoDUNEDataUtils.cxx.

                                                                        {
  this->reconfigure(p);
}

protoana::ProtoDUNEDataUtils::~ProtoDUNEDataUtils

(

)

Definition at line 16 of file ProtoDUNEDataUtils.cxx.

                                               {

}

Member Function Documentation

bool protoana::ProtoDUNEDataUtils::CheckTimeStampConsistencyForAPAs	(	art::Event const &	evt,
		std::set< int >	apas,
		ULong64_t &	timestamp,
		ULong64_t &	timestamp2,
		int &	apainconsist
	)		const

Check for consistency of timestamp values for a set of APAs. True if consistent, false if there are mismatches.

Definition at line 120 of file ProtoDUNEDataUtils.cxx.

{
  art::ServiceHandle<dune::PdspChannelMapService> channelMap;
  std::vector<art::Ptr<raw::RDTimeStamp> > TSlist;
  auto RawTSListHandle = evt.getHandle< std::vector<raw::RDTimeStamp> >(fRawDigitTimeStampTag);
  art::fill_ptr_vector(TSlist, RawTSListHandle);

  timestamp = 0;
  timestamp2 = 0;
  apainconsist = 0;
  bool tsi = false;

  for(auto const & tsptr : TSlist) 
    {
      const raw::RDTimeStamp & rdts = *tsptr;
      uint16_t chan = rdts.GetFlags();
      int iapa = channelMap->APAFromOfflineChannel(chan);
      if (apas.find(iapa) != apas.end())
        {
          ULong64_t ts = rdts.GetTimeStamp();
          if (tsi)
            {
              if (ts != timestamp)
                {
                  apainconsist = iapa;
                  timestamp2 = ts;
                  return false;
                }
            }
          else
            {
              tsi = true;
              timestamp = ts;
              timestamp2 = ts;
            }
        }
    }
  //std::cout << "All timestamps equal at: " << timestamp << " for APAs: ";
  //for (auto const &iter : apas)
  // {
  //  std::cout << iter << " "; 
  //}
  //std::cout << std::endl;
  return true;
}

int protoana::ProtoDUNEDataUtils::GetNActiveFembsForAPA	(	art::Event const &	evt,
		int	apa
	)		const

Get number of active fembs in an APA.

Definition at line 49 of file ProtoDUNEDataUtils.cxx.

                                                                                         {


  // Get pd channel map
  art::ServiceHandle<dune::PdspChannelMapService> channelMap;

  // set only saves unique elements
  std::set<int> apaset;

  // Get raw digits time stamps
  std::vector<art::Ptr<raw::RDTimeStamp> > digitTSlist;
  

  // Get raw digits
  std::vector<art::Ptr<raw::RawDigit> > digitlist;

  auto RawdigitListHandle = evt.getHandle< std::vector<raw::RawDigit> >(fRawDigitTag);
  if (RawdigitListHandle){
        
    art::fill_ptr_vector(digitlist, RawdigitListHandle);  

    for(auto const & dptr : digitlist) {
      const raw::RawDigit& digit = *dptr;
      // Get the channel number for this digit
      uint32_t chan = digit.Channel();
    
      int iapa = channelMap->APAFromOfflineChannel(chan);
      if(iapa != apa) continue;
      // Get the channel FEMB and WIB
      int WIB = channelMap->WIBFromOfflineChannel(chan); // 0-4
      int FEMB = channelMap->FEMBFromOfflineChannel(chan); // 1-4
      //int FEMBchan = channelMap->FEMBChannelFromOfflineChannel(chan);
      int iFEMB = ((WIB*4)+(FEMB-1)); //index of the FEMB 0-19

      apaset.insert(iFEMB);
    }
  }

  else{ // if raw digits have been dropped use RDTimeStamps instead

    auto RawdigitTSListHandle = evt.getHandle< std::vector<raw::RDTimeStamp> >(fRawDigitTimeStampTag);
    art::fill_ptr_vector(digitTSlist, RawdigitTSListHandle);  
        
    for(auto const & dptr : digitTSlist) {
        
      const raw::RDTimeStamp & digit = *dptr;
    
      // Get the channel number for this digit
      uint16_t chan = digit.GetFlags();
    
      int iapa = channelMap->APAFromOfflineChannel(chan);
      if(iapa != apa) continue;
      // Get the channel FEMB and WIB
      int WIB = channelMap->WIBFromOfflineChannel(chan); // 0-4
      int FEMB = channelMap->FEMBFromOfflineChannel(chan); // 1-4
      //int FEMBchan = channelMap->FEMBChannelFromOfflineChannel(chan);
      int iFEMB = ((WIB*4)+(FEMB-1)); //index of the FEMB 0-19

      apaset.insert(iFEMB);
    }
  }
  return (apaset.size());
}

bool protoana::ProtoDUNEDataUtils::IsBeamTrigger

(

art::Event const &

evt

)

const

Returns true if the ProtoDUNE trigger says this is a beam trigger

Definition at line 27 of file ProtoDUNEDataUtils.cxx.

                                                                        {

  bool isBeam = false;

  // Accessing the trigger information as done in DataPrepModule
  // The information is stored in the time stamps
  auto timeStamps = evt.getHandle<std::vector<raw::RDTimeStamp>>(fTimingTag);

  // Return false if we have no time stamps
  if (!timeStamps) return isBeam;
  // We should only have one RDTimeStamp
  if (timeStamps->size() > 1) return isBeam;
  
  // Access the trigger information. Beam trigger flag = 0xc
  const raw::RDTimeStamp& timeStamp = timeStamps->at(0);
  isBeam = (timeStamp.GetFlags() == 0xc);
  
  return isBeam;
}

void protoana::ProtoDUNEDataUtils::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 20 of file ProtoDUNEDataUtils.cxx.

                                                                      {
  fTimingTag            = p.get<art::InputTag>("TimingTag");
  fRawDigitTag          = p.get<art::InputTag>("RawDigitTag");
  fRawDigitTimeStampTag = p.get<art::InputTag>("RawDigitTimeStampTag");
}

Member Data Documentation

art::InputTag protoana::ProtoDUNEDataUtils::fRawDigitTag

private

Definition at line 46 of file ProtoDUNEDataUtils.h.

art::InputTag protoana::ProtoDUNEDataUtils::fRawDigitTimeStampTag

private

Definition at line 47 of file ProtoDUNEDataUtils.h.

art::InputTag protoana::ProtoDUNEDataUtils::fTimingTag

private

Definition at line 45 of file ProtoDUNEDataUtils.h.

---

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

• duneprototypes/duneprototypes/Protodune/singlephase/DataUtils/ProtoDUNEDataUtils.h
• duneprototypes/duneprototypes/Protodune/singlephase/DataUtils/ProtoDUNEDataUtils.cxx