NueSkimmer.h

Go to the documentation of this file.

//
//  NueSkimmer.h
//
//  Created by Brian Rebel on 2/18/16.
//

#ifndef NueSkimmer_h
#define NueSkimmer_h

#include "art/Framework/Core/EDProducer.h"
#include "art/Framework/Principal/Event.h"
#include "canvas/Utilities/InputTag.h"
#include "fhiclcpp/ParameterSet.h"

#include "RecoBase/Cluster.h"
#include "RecoBase/Shower.h"
#include "ShowerLID/ShowerLID.h"
#include "AnalysisSkimmer/skimmers/BaseSkimmer.h"

namespace skim {
  class NueSkimmer;
}

namespace me {
  class SlcME;
  class TrkME;
}

// this module is a producer because it is assumed to be after a path that
// contains the CellHitSkimmer, which has already removed events without
// interesting slices
class skim::NueSkimmer : public skim::BaseSkimmer {
  
public:
  
  explicit NueSkimmer(fhicl::ParameterSet const & p);

  void reconfigure(fhicl::ParameterSet const & p);

  void Skim(art::Event           & e,
            art::EDProducer const& prod);
  
private:
  
  void CopyShowerLID(art::Event                                     & e,
                     art::EDProducer                           const& prod,
                     std::vector< art::Ptr<rb::Shower> >             const& showerIn,
                     std::vector<slid::ShowerLID>                   & slidOut,
                     std::vector<rb::Cluster>                               & sliceOut,
                     std::vector<rb::Shower>                                & showerOut,
                     art::Assns<rb::Cluster, slid::ShowerLID>       & slcSLDAssn,
                     art::Assns<rb::Shower, slid::ShowerLID>            & shwSLDAssn);

  void CopyProngCVN(art::Event                                   & e,
                    art::EDProducer                         const& prod,
                    std::vector< art::Ptr<rb::Prong> >      const& prongIn,
                    std::vector<rb::PID>                         & pidOut,
                    std::vector<rb::Cluster>                     & sliceOut,
                    std::vector<rb::Prong>                       & prongOut,
                    art::Assns<rb::Cluster, rb::PID>             & slcPIDAssn,
                    art::Assns<rb::Prong,   rb::PID>             & pngPIDAssn);

  void CopyShowerPngAssn(art::Event                                   & e,
                         art::EDProducer                         const& prod,
                         std::vector< art::Ptr<rb::Shower> >     const& showerIn,
                         std::vector< art::Ptr<rb::Prong>  >     const& prongIn,
                         std::vector<rb::Shower>                      & showerOut,
                         std::vector<rb::Prong>                       & prongOut,
                         art::Assns<rb::Shower,  rb::Prong>           & shwPngAssn);

  void CopyMichelSlice(art::Event                             & e,
                       art::EDProducer                   const& prod,
                       std::vector<art::Ptr<me::SlcME> > const& objIn,
                       std::vector<me::SlcME>                 & objOut,
                       std::vector<rb::Cluster>               & sliceOut,
                       art::Assns<rb::Cluster, me::SlcME>     & slcMEAssn);

  void CopyMichelTrack(art::Event                             & e,
                       art::EDProducer                   const& prod,
                       std::vector<art::Ptr<me::TrkME> > const& objIn,
                       std::vector<me::TrkME>                 & objOut,
                       std::vector<rb::Cluster>               & sliceOut,
                       art::Assns<rb::Cluster, me::TrkME>     & slcMEAssn);

  std::string   fVertexLabel;    ///< label of module creating vertex
  std::string   fShowerLabel;    ///< label of module creating showers
  art::InputTag fProngTag;       ///< input tag of module creating showers
  std::string   fLIDLabel;       ///< label of module creating lids
  std::string   fLEMLabel;       ///< label of module creating lems
  std::string   fCosRejLabel;    ///< label of module creating cosrejs
  std::string   fCAFLabel;       ///< label of module creating StandardRecords
  std::string   fGeneratorLabel; ///< label of module generating events
  std::string   fFluxWeightLabel;///< label of module creating flux weights
  std::string   fCVNLabel;       ///< label of module creating cvns
  std::string   fNueVetoLabel;   ///< label of module creating nuevetos
  std::string   fMichelLabel;    ///< label of module creating slice Michel e's
  std::string   fCVNProngLabel;  ///< label of module creating prongs cvn
};


#endif /* NueSkimmer_h */