FillTruth.h

Go to the documentation of this file.

#ifndef CAF_FILLTRUTH_H
#define CAF_FILLTRUTH_H

#include "StandardRecord/StandardRecord.h"

#include "NovaDAQConventions/DAQConventions.h"

#include "canvas/Persistency/Common/PtrVector.h"

#include <vector>
#include <iterator>
#include <nusimdata/SimulationBase/MCGeneratorInfo.h>
#include "nusimdata/SimulationBase/MCTruth.h"

namespace rb{class Cluster; class CellHit;}
namespace sim{class Particle; class TrueEnergy;}
namespace cheat{struct NeutrinoEffPur; struct TrackIDE;}
namespace rwgt{class GENIEReweightTable;}
namespace g4rwgt{class G4WeightTable;}
namespace fxwgt{class FluxWeights;}
namespace g4rwgt{class G4WeightTable;}

namespace caf
{
  /// Function to fill the particle truth for a set of hits
  SRParticleTruth FillParticleTruth(const std::vector<rb::Cluster>& sliceList,
                                    const art::PtrVector<rb::CellHit>& hits,
                                    std::vector<cheat::TrackIDE>& allTracks,
                                    std::vector<cheat::TrackIDE>& sliceTracks,
                                    std::vector<cheat::TrackIDE> &allTracksBirks,
                                    std::vector<cheat::TrackIDE> &sliceTracksBirks,
                                    int sliceIdx);
   
  void AddParticleToVec(const sim::Particle& part,
                        std::vector<cheat::TrackIDE>& allTracks,
                        std::vector<cheat::TrackIDE>& sliceTracks,
                        std::vector<cheat::TrackIDE> &allTracksBirks,
                        std::vector<cheat::TrackIDE> &sliceTracksBirks,
                        std::vector<SRTrueParticle>* vec, 
                        const std::vector<sim::TrueEnergy>& TrueEnergies);

  float FindDaughterVisE(const sim::Particle& part,
                        std::vector<cheat::TrackIDE>& allTracks);

  float FindDaughterVisENonRecur(const sim::Particle& part,
                         std::vector<cheat::TrackIDE>& allTracks);

  void AddCosmicTruthToVec(const cheat::NeutrinoEffPur* effPur,
                           std::vector<cheat::TrackIDE>& allTracks,
                           std::vector<cheat::TrackIDE>& sliceTracks,
                           std::vector<cheat::TrackIDE> &allTracksBirks,
                           std::vector<cheat::TrackIDE> &sliceTracksBirks,
                           std::vector<SRCosmic>* vec);

  void FindAndAddMichels(std::vector<const sim::Particle*> particles,
                         std::vector<cheat::TrackIDE>& allTracks,
                         std::vector<SRTrueMichelE>* michelVec);

  void AddTrueMichelEToVec(const sim::Particle& michel,
                           const sim::Particle& motherMuon,
                           std::vector<cheat::TrackIDE>& allTracks,
                           std::vector<SRTrueMichelE>* vec);

  void AddPreFSI( const art::Ptr<simb::MCTruth>& truth, SRNeutrino& nu );

  double TrueNeutrinoDistance(novadaq::cnv::DetId det, const SRNeutrino& nu);

  std::vector<SRGenieWeights> GenieReweightTable(const rwgt::GENIEReweightTable& table);

  SRFluxWeights FluxReweights(const fxwgt::FluxWeights& flxwgts);

  SRGeant4Weights Geant4Reweights(const g4rwgt::G4WeightTable& g4wgts);

  caf::generator_ CAFGeneratorEnum(simb::Generator_t simbGeneratorEnum);

  std::vector<unsigned int> DecodeGeneratorVersion(const std::string& versionString);
}

#endif