15 #include "TProfile2D.h" 33 #include "Utilities/AssociationUtil.h" 35 #include "Utilities/func/MathUtil.h" 59 #include "TMVA/Tools.h" 60 #include "TMVA/Reader.h" 61 #include "TMVA/MethodCuts.h" 67 namespace rb{
class Cluster;
class Track;}
68 namespace sim{
class Particle;}
69 namespace simb{
class MCFlux;
class MCTruth;
class MCNeutrino;}
171 _otree = tfs->
make<TTree>(
"Xbeam",
"Xbeam particle");
283 if( mclist->size()!=1 ){
330 if( nu.
CCNC() != 0 )
return;
332 if( flux->
fntype != 12 )
return;
343 bool matchSlice=
false;
345 for(
unsigned int sliceIdx = 0; sliceIdx < slices->size(); ++sliceIdx){
351 for(
size_t h = 0;
h < chits->size(); ++
h)
364 std::vector< art::Ptr<simb::MCFlux> > fluxs = fmf.at(0);
366 if(fluxs.size() == 1){
double E(const int i=0) const
std::string fCellHitLabel
label for generator information
back track the reconstruction to the simulation
double Theta() const
angle between incoming and outgoing leptons, in radians
const simb::MCNeutrino & GetNeutrino() const
double Py(const int i=0) const
void analyze(const art::Event &evt)
const rb::RecoHit MakeRecoHit(art::Ptr< rb::CellHit > const &chit, double const &w)
SubRunNumber_t subRun() const
std::vector< NeutrinoEffPur > SliceToNeutrinoInteractions(const std::vector< const rb::CellHit * > &sliceHits, const std::vector< const rb::CellHit * > &allHits, bool sortPur=false) const
Given a collection of hits (often a slice), returns vector of structures of neutrino interactions...
MaybeLogger_< ELseverityLevel::ELsev_info, false > LogInfo
const simb::MCParticle & Nu() const
double Pt() const
transverse momentum of interaction, in GeV/c
double Px(const int i=0) const
std::string fMCFluxLabel
label for slices
rb::RecoHit MakeRecoHit(rb::CellHit const &cellhit, double w)
A collection of associated CellHits.
charged current quasi-elastic
object containing MC flux information
Calibrated quantities relying on position in the orthogonal view. To generate a rb::CellHit from a rb...
int InteractionType() const
DEFINE_ART_MODULE(ROCKMRE)
const simb::MCParticle & Lepton() const
void push_back(Ptr< U > const &p)
T get(std::string const &key) const
Vertex location in position and time.
Perform a "2 point" Hough transform on a collection of hits.
Example routine for calculating the "ultimate e-m mother" of a particle in a simulated event...
This class describes a particle created in the detector Monte Carlo simulation.
double Vx(const int i=0) const
T * make(ARGS...args) const
TProfile * slEff_Etrue
label for cell hits
Data resulting from a Hough transform on the cell hit positions.
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
double Pz(const int i=0) const
double Vz(const int i=0) const
bool IsNoise() const
Is the noise flag set?
double totgoodpot
normalized by 10^12 POT
Event generator information.
void endSubRun(const art::SubRun &sr)
double Vy(const int i=0) const
Encapsulate the geometry of one entire detector (near, far, ndos)
void reconfigure(const fhicl::ParameterSet &p)