SRTrueInteraction.h

Go to the documentation of this file.

 
 
#ifndef DUNEANAOBJ_SRTRUEINTERACTION_H
#define DUNEANAOBJ_SRTRUEINTERACTION_H
 
#include <limits>
 
#include "duneanaobj/StandardRecord/SREnums.h"
#include "duneanaobj/StandardRecord/SRTrueParticle.h"
#include "duneanaobj/StandardRecord/SRVector3D.h"
 
// // This is pretty ugly. The proper fix is to get a new version of castxml that
// // is built with at least clang v9.
// // See https://github.com/CastXML/CastXML/issues/178
// #ifndef __castxml__
// #include <string>
// #else
// namespace std{class string{};}
// #endif
 
// a less ugly hack, works with castxml v0_4_2 and e20, c7 and e26 but not c14
// not needed for castxml v0_5_1.  Now e20, c7, c14 work but not e26
 
//#ifdef __GNUC__
//#undef _GLIBCXX_HAVE_BUILTIN_IS_CONSTANT_EVALUATED
//#endif
 
#include <string>
 
#include <vector>
 
namespace caf
{
 
   class SRDialWeights
   {
   public:
     std::vector<float> weights;
   };
 
   class SRTrueInteraction
   {
     private:
       // just to keep the typing under control below
       static constexpr float NaN = std::numeric_limits<float>::signaling_NaN();
 
     public:
       long int  id       = -1;
       
       long int  genieIdx = -1;
 
       int   pdg     = 0;         
       int   pdgorig = 0;         
 
       bool           iscc      = false;                          
       ScatteringMode mode      = ScatteringMode::kUnknownMode;   
       int            targetPDG = 0;                              
 
       int   hitnuc   = 0;
       float removalE = NaN;        
 
       float       E         = NaN; 
       SRVector3D  vtx;             
       SRVector3D  momentum;        
       bool  isvtxcont = false;     
 
       float time         = NaN;    
       float bjorkenX     = NaN;    
       float inelasticity = NaN;    
       float Q2           = NaN;    
       float q0           = NaN;    
       float modq         = NaN;    
       float W            = NaN;    
       float t            = NaN;    
 
       // GENIE truth stuff
       bool  ischarm    = false;    
       bool  isseaquark = false;    
       int   resnum     = 0;        
       float xsec       = NaN;      
 
       float genweight = NaN;       
 
       float baseline             = NaN; 
       SRVector3D prod_vtx;              
       SRVector3D parent_dcy_mom;        
       int        parent_dcy_mode = -1;  
       int        parent_pdg      = 0;   
       float      parent_dcy_E    = NaN; 
       float      imp_weight      = NaN; 
 
       Generator                   generator = Generator::kUnknownGenerator;  
       std::vector<unsigned int>   genVersion;                                
 
       // fixme: For unknown reasons, deserializing SRTrueInteractions with this branch sometimes (but not always) causes `free: invalid pointer` segfaults.
       //        For the moment, we just disable it, but we need to track down what's going wrong at some point.
//       std::string                 genConfigString;                           ///< String associated with generator configuration. (For GENIE 3+, this is the "Comprehensive Model Configuration".)
 
       // shortcuts to avoid iterating through the particle stack for these simple quantities
       int        nproton  = 0;    
       int        nneutron = 0;    
       int        npip     = 0;    
       int        npim     = 0;    
       int        npi0     = 0;    
 
       int                         nprim = 0;        
       std::vector<SRTrueParticle> prim;             
       int                         nprefsi = 0;      
       std::vector<SRTrueParticle> prefsi;           
       int                         nsec = 0;         
       std::vector<SRTrueParticle> sec;              
 
       float                       xsec_cvwgt = NaN; 
 
       std::vector<SRDialWeights> syst_dials; 
 
   };
 
} // caf
 
#endif //DUNEANAOBJ_SRTRUEINTERACTION_H