icarus::crt::CRTEventProducer Class Reference

Inheritance diagram for icarus::crt::CRTEventProducer:

Public Member Functions
	CRTEventProducer (fhicl::ParameterSet const &p)
	CRTEventProducer (CRTEventProducer const &)=delete
	CRTEventProducer (CRTEventProducer &&)=delete
CRTEventProducer &	operator= (CRTEventProducer const &)=delete
CRTEventProducer &	operator= (CRTEventProducer &&)=delete
void	produce (art::Event &e) override
void	beginJob () override
void	endJob () override
void	reconfigure (fhicl::ParameterSet const &p)

Private Attributes
std::string	fInFile
float	fTimeWindow
icarus::crt::CRTPreProcessTree *	fAnaTree
std::map< size_t, size_t >	fEventEntryStart
std::map< size_t, size_t >	fEventEntryEnd

Detailed Description

Definition at line 47 of file CRTEventProducer_module.cc.

Constructor & Destructor Documentation

icarus::crt::CRTEventProducer::CRTEventProducer ( fhicl::ParameterSet const &  p )

explicit

Definition at line 82 of file CRTEventProducer_module.cc.

  : EDProducer{p}
  // Initialize member data here, if know don't want to reconfigure on the fly
  {
 
   // Call appropriate produces<>() functions here.
    produces< std::vector<icarus::crt::CRTData> >();
    
    reconfigure(p);

    /*std::cout << "opening file " << fInFile << std::endl;
    TFile fin(fInFile.c_str());
    TTree* tree=(TTree*)fin.FindObjectAny("anaTree");
    if(tree==nullptr) std::cout << "anaTree not found!" << std::endl;
    fAnaTree = new icarus::crt::CRTPreProcessTree(tree);
    std::cout << "constructed CRTPreProcessTree with " << fAnaTree->GetNEntries() << " entries" << std::endl;
    std::cout << "earliest timestamp: " << fAnaTree->GetAbsTime(0) << std::endl;
*/
  } // CRTEventProducer()

icarus::crt::CRTEventProducer::CRTEventProducer ( CRTEventProducer const &  )

delete

icarus::crt::CRTEventProducer::CRTEventProducer ( CRTEventProducer &&  )

delete

Member Function Documentation

void icarus::crt::CRTEventProducer::beginJob ( )

override

Definition at line 110 of file CRTEventProducer_module.cc.

  {
        std::cout << "opening file " << fInFile << std::endl;
        TFile fin(fInFile.c_str());
        TTree* tree=(TTree*)fin.FindObjectAny("anaTree");
        if(tree==nullptr) std::cout << "anaTree not found!" << std::endl;
        fAnaTree = new icarus::crt::CRTPreProcessTree(tree);

        std::cout << "running 'beginJob'" << std::endl;
        size_t event=0;
        uint64_t ti = fAnaTree->GetAbsTime(0);
        fEventEntryStart[0] = 0;

        std::cout << "time slicing tree with " << fAnaTree->GetNEntries() << " entries" << std::endl;

        for(size_t ientry =0; ientry<fAnaTree->GetNEntries(); ientry++){
                if(fAnaTree->GetAbsTime(ientry)>=ti &&
                   fAnaTree->GetAbsTime(ientry)<ti+fTimeWindow*1.0e9)   
                        continue;

                else if(fAnaTree->GetAbsTime(ientry)>=ti+fTimeWindow*1.0e9 ||
                        ientry == fAnaTree->GetNEntries()-1){
                        ti += fTimeWindow*1.0e9;
                        fEventEntryEnd[event] = ientry-1;
                        event++;
                        if (ientry < fAnaTree->GetNEntries()-1)
                                fEventEntryStart[event] = ientry;
                        
                }
        }

        std::cout << "found " << event << " time slices" << std::endl;
        if(fEventEntryEnd.size()!=fEventEntryStart.size())
                std::cout << "event map start/end size mismatch!" << std::endl;

        fin.Close();
        delete fAnaTree;

  } // CRTEventProducer::beginJob()

void icarus::crt::CRTEventProducer::endJob ( )

override

Definition at line 189 of file CRTEventProducer_module.cc.

  {

  } // CRTEventProducer::endJob()

CRTEventProducer& icarus::crt::CRTEventProducer::operator= ( CRTEventProducer const &  )

delete

CRTEventProducer& icarus::crt::CRTEventProducer::operator= ( CRTEventProducer &&  )

delete

void icarus::crt::CRTEventProducer::produce ( art::Event &  e )

override

Definition at line 150 of file CRTEventProducer_module.cc.

  {
    TFile fin(fInFile.c_str());
    TTree* tree=(TTree*)fin.FindObjectAny("anaTree");
    if(tree==nullptr) std::cout << "anaTree not found!" << std::endl;
    fAnaTree = new icarus::crt::CRTPreProcessTree(tree);

    std::unique_ptr< std::vector<icarus::crt::CRTData> > crtdata( new std::vector<icarus::crt::CRTData>);

    size_t eveId = (size_t) event.event();
    std::cout << "processing event " << eveId << " with " << fEventEntryEnd[eveId]+1-fEventEntryStart[eveId] << " triggers" << std::endl;

    int entry=0;
    for(size_t ientry=fEventEntryStart[eveId]; ientry<=fEventEntryEnd[eveId]; ientry++){
        fAnaTree->Load(ientry);
        CRTData data;
        data.fMac5 = fAnaTree->Mac5();
        data.fTs0 = fAnaTree->GetAbsTime();
        data.fTs1 = data.fTs0;
        for(size_t ichan=0; ichan<32; ichan++){
                if(fAnaTree->Above(ichan)&&fAnaTree->Active(ichan)) {
                        data.fAdc[ichan]  = fAnaTree->PE(ichan);
                }//if channel active and above threshold
                else
                        data.fAdc[ichan] = 0;
        }//for channels

        crtdata->push_back(data);
        entry++;

    }//for entries in time slice

    delete fAnaTree;
    fin.Close();

    event.put(std::move(crtdata));

  } // CRTEventProducer::produce()

void icarus::crt::CRTEventProducer::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 102 of file CRTEventProducer_module.cc.

  {
    //fCrtModuleLabel = (p.get<art::InputTag> ("CrtModuleLabel"));
    fInFile = p.get<string>("infile");
    fTimeWindow = p.get<float>("time_window");

  } // CRTEventProducer::reconfigure()

Member Data Documentation

icarus::crt::CRTPreProcessTree* icarus::crt::CRTEventProducer::fAnaTree

private

Definition at line 75 of file CRTEventProducer_module.cc.

std::map<size_t,size_t> icarus::crt::CRTEventProducer::fEventEntryEnd

private

Definition at line 78 of file CRTEventProducer_module.cc.

std::map<size_t,size_t> icarus::crt::CRTEventProducer::fEventEntryStart

private

Definition at line 77 of file CRTEventProducer_module.cc.

std::string icarus::crt::CRTEventProducer::fInFile

private

Definition at line 73 of file CRTEventProducer_module.cc.

float icarus::crt::CRTEventProducer::fTimeWindow

private

Definition at line 74 of file CRTEventProducer_module.cc.

---

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

• CRTEventProducer_module.cc