PPFXFluxReader.h

Go to the documentation of this file.

#include "art/Framework/Core/FileBlock.h"
#include "art/Framework/Core/ProductRegistryHelper.h"
#include "art/Framework/IO/Sources/SourceHelper.h"
#include "art/Framework/Principal/fwd.h"
#include "fhiclcpp/ParameterSet.h"

#include <fstream>
#include <map>
#include <vector>

#include "FluxInterface.h"

class TH1D;
class TTree;
class TFile;

namespace fluxr {
  class PPFXFluxReader {
  public:
    // Required constructor
    PPFXFluxReader(fhicl::ParameterSet const& pset,
                   art::ProductRegistryHelper& helper,
                   art::SourceHelper const& pm);

    // Required by FileReaderSource:
    void closeCurrentFile();
    void readFile(std::string const& name, art::FileBlock*& fb);
    bool readNext(art::RunPrincipal* const& inR,
                  art::SubRunPrincipal* const& inSR,
                  art::RunPrincipal*& outR,
                  art::SubRunPrincipal*& outSR,
                  art::EventPrincipal*& outE);

  private:
    art::SourceHelper const& fSourceHelper;
    art::SubRunID fSubRunID;

    uint32_t fEventCounter;
    uint32_t fEntry;
    int fMaxEvents;         //fhicl parameter.  Maximum number of events.
    uint32_t fSkipEvents;   // fhicl parameter.  Number of events to skip.
    std::string fInputType; //fhicl parameter.  Maximum number of events.
    float fPOT;

    FluxInterface* fFluxDriver;
    TFile* fFluxInputFile;
    //int                           fNuPdgCode[4];
    TH1D* fHFlux[4];
    TH1D* fHFluxParent[4][4];
    TH1D* fHFluxSec[4][5];

    fhicl::ParameterSet fConfigPS;

    art::TypeLabel fTLmctruth;
    art::TypeLabel fTLmcflux;
    art::TypeLabel fTLdk2nu;
    art::TypeLabel fTLnuchoice;
  }; // PPFXFluxReader
}