MyMuoncountertimes::Muoncountertimes Class Reference

Inheritance diagram for MyMuoncountertimes::Muoncountertimes:

Public Member Functions
	Muoncountertimes (fhicl::ParameterSet const &p)
virtual	~Muoncountertimes ()
void	beginJob ()
void	endJob ()
void	reconfigure (fhicl::ParameterSet const &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 Member Functions
void	ResetVars ()

Private Attributes
int	run
int	subrun
int	event
int	event_Count = 0
std::vector< TH1D * >	fHist_Vec
int	fLower_Block1
int	fUpper_Block1
int	fLower_Block2
int	fUpper_Block2
int	fCoincidence_Window
double	fCombinedTimeDelay
std::string	fCounterModuleLabel
std::string	fFragType
std::string	fRawDataLabel
std::string	fPTBMapFile
std::string	fPTBMapDir
std::map< int, int >	fPTBMap

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 52 of file MuonCounterTimes_module.cc.

Constructor & Destructor Documentation

MyMuoncountertimes::Muoncountertimes::Muoncountertimes ( fhicl::ParameterSet const &  p )

explicit

Definition at line 104 of file MuonCounterTimes_module.cc.

  : EDAnalyzer(pset)
{
   // Read in the parameters from the .fcl file.
    this->reconfigure(pset);
}

MyMuoncountertimes::Muoncountertimes::~Muoncountertimes ( )

virtual

Definition at line 133 of file MuonCounterTimes_module.cc.

{
  // Clean up dynamic memory and other resources here.
}

Member Function Documentation

void MyMuoncountertimes::Muoncountertimes::analyze

(

const art::Event &

evt

)

Definition at line 138 of file MuonCounterTimes_module.cc.

{
  // Implementation of required member function here.
  ResetVars();

  art::ServiceHandle<geo::Geometry> geom;

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

  // New way to get counter hits
  /*
  art::Handle< std::vector< raw::ExternalTrigger> > externalTriggerListHandle;
  evt.getByLabel(fCounterModuleLabel, externalTriggerListHandle);
  std::vector< art::Ptr< raw::ExternalTrigger> > trigs;
  art::fill_ptr_vector(trigs,externalTriggerListHandle);
 */

  art::Handle< artdaq::Fragments > fragment;

  evt.getByLabel(fRawDataLabel, fFragType, fragment);

  // Check if there is PTB data in this event
  // Don't crash code if not present, just don't save anything
  try { fragment->size(); }
  catch(std::exception const&) {
    mf::LogWarning("MuonCounterTimes") << "WARNING: Raw PTB data not found in event " << evt.event();
    return;
  }

 // Check that the data is valid
  if(!fragment.isValid()){
    mf::LogError("MuonCounterTimes") << "Run: " << evt.run()
                                 << ", SubRun: " << evt.subRun()
                                 << ", Event: " << evt.event()
                                 << " is NOT VALID";
    throw cet::exception("raw NOT VALID");
    return;
  }

  std::vector<std::vector<int>>       coincidences(2);
  std::vector<std::vector<int>>       trig_ID_Vec(2);
  std::vector<std::vector<long long>> trig_Times_Vec(2); 
  lbne::PennMicroSlice::Payload_Timestamp *FirstPTBTimestamp = nullptr;
  auto trigs = DAQToOffline::PennFragmentToExternalTrigger(*fragment, fPTBMap, FirstPTBTimestamp);
  
  unsigned int total_Hits = trigs.size();

  for(unsigned int i = 0; i < total_Hits; i++)
  {
    int       auxdetid  = trigs.at(i).GetTrigID();
    long long trig_Time = trigs.at(i).GetTrigTime();   

    if(auxdetid>=fLower_Block1 && auxdetid<=fUpper_Block1)
    {
      trig_ID_Vec.at(0).push_back(auxdetid); 
      trig_Times_Vec.at(0).push_back(trig_Time);
    }
    else if(auxdetid>=fLower_Block2 && auxdetid<=fUpper_Block2)
    {
      trig_ID_Vec.at(1).push_back(auxdetid); 
      trig_Times_Vec.at(1).push_back(trig_Time);
    }
  }

  for(unsigned int i = 0; i < trig_ID_Vec.at(0).size(); i++)
  {
    for(unsigned int j = 0; j < trig_ID_Vec.at(1).size(); j++)
    {
      if(j<=i)
      {
        if(std::abs(trig_Times_Vec.at(0).at(i)-trig_Times_Vec.at(1).at(j))<fCoincidence_Window)
        {
          coincidences.at(0).push_back(trig_ID_Vec.at(0).at(i));
          coincidences.at(1).push_back(trig_ID_Vec.at(1).at(j));
        }
      }
    }
  }

  for(unsigned int i = 0; i < coincidences.at(0).size(); i++)
  {
    fHist_Vec.at(coincidences.at(0).at(i)-fLower_Block1)->Fill(coincidences.at(1).at(i));
  }
  event_Count+=1;
}

void MyMuoncountertimes::Muoncountertimes::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 226 of file MuonCounterTimes_module.cc.

{
  // Implementation of optional member function here.

  fHist_Vec.resize(fUpper_Block1-fLower_Block1+1);

  art::ServiceHandle<art::TFileService> tfs;

  for(unsigned int i = 0; i < fHist_Vec.size(); i++)
  {
    TString key     = Form("h%i", i);
    TH1D *tempHist  = tfs->make<TH1D>(key, key, fUpper_Block2-fLower_Block2+1, fLower_Block2, fUpper_Block2+1); 
    fHist_Vec.at(i) = tempHist; 
  }
}

void MyMuoncountertimes::Muoncountertimes::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 242 of file MuonCounterTimes_module.cc.

{
  for(unsigned int i = 0; i < fHist_Vec.size(); i++)
  {
    TString fHist_Title1 = Form("Coincidence Plots, run %i, ", run); 
    TString fHist_Title2 = Form("counter %i.", i+fLower_Block1);
    TString fHist_Title  = fHist_Title1 + fHist_Title2;

    TString fHist_Name1 = Form("Run%i", run);
    TString fHist_Name2 = Form("Counter%i", i+fLower_Block1);
    TString fHist_Name  = fHist_Name1 + fHist_Name2;
 
    fHist_Vec.at(i)->GetXaxis()->SetTitle("Coincidence Partner");
    fHist_Vec.at(i)->GetYaxis()->SetTitle("Frequency");
    fHist_Vec.at(i)->GetYaxis()->SetTitleOffset(1.);

    fHist_Vec.at(i)->SetTitle(fHist_Title);
    fHist_Vec.at(i)->SetName(fHist_Name);

    fHist_Vec.at(i)->Sumw2();
  }
}

void MyMuoncountertimes::Muoncountertimes::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 111 of file MuonCounterTimes_module.cc.

{   // Read parameters from the .fcl file. The names in the arguments
  // to p.get<TYPE> must match names in the .fcl file.

  fLower_Block1       = p.get<int>("LowerBlock1");
  fUpper_Block1       = p.get<int>("UpperBlock1");
  fLower_Block2       = p.get<int>("LowerBlock2");
  fUpper_Block2       = p.get<int>("UpperBlock2");
  fCoincidence_Window = p.get<int>("CoincidenceWindow");
  fCounterModuleLabel = p.get< std::string >("CounterModuleLabel");
  //    fLArG4ModuleLabel(p.get< std::string >("LArGeantModuleLabel", "largeant"));
  fCombinedTimeDelay = p.get< double >("CombinedTimeDelay");
  fFragType = p.get<std::string>("FragType");
  fRawDataLabel = p.get<std::string>("DataLabel");
  fPTBMapFile   = p.get<std::string>("PTBMapFile");
  fPTBMapDir    = p.get<std::string>("PTBMapDir");

  DAQToOffline::BuildPTBChannelMap(fPTBMapDir, fPTBMapFile, fPTBMap);

  return;
}

void MyMuoncountertimes::Muoncountertimes::ResetVars ( void  )

private

Definition at line 265 of file MuonCounterTimes_module.cc.

                                {

  run = -99999;
  subrun = -99999;
  event = -99999;
}

Member Data Documentation

int MyMuoncountertimes::Muoncountertimes::event

private

Definition at line 83 of file MuonCounterTimes_module.cc.

int MyMuoncountertimes::Muoncountertimes::event_Count = 0

private

Definition at line 84 of file MuonCounterTimes_module.cc.

int MyMuoncountertimes::Muoncountertimes::fCoincidence_Window

private

Definition at line 93 of file MuonCounterTimes_module.cc.

double MyMuoncountertimes::Muoncountertimes::fCombinedTimeDelay

private

Definition at line 94 of file MuonCounterTimes_module.cc.

std::string MyMuoncountertimes::Muoncountertimes::fCounterModuleLabel

private

Definition at line 95 of file MuonCounterTimes_module.cc.

std::string MyMuoncountertimes::Muoncountertimes::fFragType

private

Definition at line 96 of file MuonCounterTimes_module.cc.

std::vector<TH1D*> MyMuoncountertimes::Muoncountertimes::fHist_Vec

private

Definition at line 87 of file MuonCounterTimes_module.cc.

int MyMuoncountertimes::Muoncountertimes::fLower_Block1

private

Definition at line 89 of file MuonCounterTimes_module.cc.

int MyMuoncountertimes::Muoncountertimes::fLower_Block2

private

Definition at line 91 of file MuonCounterTimes_module.cc.

std::map<int,int> MyMuoncountertimes::Muoncountertimes::fPTBMap

private

Definition at line 101 of file MuonCounterTimes_module.cc.

std::string MyMuoncountertimes::Muoncountertimes::fPTBMapDir

private

Definition at line 99 of file MuonCounterTimes_module.cc.

std::string MyMuoncountertimes::Muoncountertimes::fPTBMapFile

private

Definition at line 98 of file MuonCounterTimes_module.cc.

std::string MyMuoncountertimes::Muoncountertimes::fRawDataLabel

private

Definition at line 97 of file MuonCounterTimes_module.cc.

int MyMuoncountertimes::Muoncountertimes::fUpper_Block1

private

Definition at line 90 of file MuonCounterTimes_module.cc.

int MyMuoncountertimes::Muoncountertimes::fUpper_Block2

private

Definition at line 92 of file MuonCounterTimes_module.cc.

int MyMuoncountertimes::Muoncountertimes::run

private

Definition at line 81 of file MuonCounterTimes_module.cc.

int MyMuoncountertimes::Muoncountertimes::subrun

private

Definition at line 82 of file MuonCounterTimes_module.cc.

---

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

• dunetpc/dune/MuonCounterTimes/MuonCounterTimes_module.cc