eic-opticks/sysrap/SFastSim_Debug.hh

0001 #pragma once
0002 /**
0003 SFastSim_Debug.h
0004 ===========================
0005
0006 HMM: this has been mostly replaced by ModelTrigger_Debug
0007
0008 **/
0009
0010 #include "plog/Severity.h"
0011 #include <vector>
0012 #include "SYSRAP_API_EXPORT.hh"
0013
0014 struct SYSRAP_API SFastSim_Debug
0015 {
0016     static const plog::Severity LEVEL ;
0017     static std::vector<SFastSim_Debug> record ;
0018     static constexpr const unsigned NUM_QUAD = 4u ;
0019     static constexpr const char* NAME = "SFastSim_Debug.npy" ;
0020     static constexpr int LIMIT = 100000 ;
0021     static void Save(const char* dir);
0022     void add();
0023     void fill(double value);
0024
0025     double posx ;  // pos = fs[:,0,:3]
0026     double posy ;
0027     double posz ;
0028     double time ;  // tim = fs[:,0,3]
0029
0030     double dirx ;  // mom = fs[:,1,:3]
0031     double diry ;
0032     double dirz ;
0033     double dist1 ; // ds1 = fs[:,1,3]
0034
0035     double polx ;  // pol = fs[:,2,:3]
0036     double poly ;
0037     double polz ;
0038     double dist2 ; // ds2 = fs[:,2,3]
0039
0040     // TODO: make the layout closer to sphoton.h, stop wasting bits
0041
0042     double ModelTrigger ; // trg = fs[:,3,0].astype(np.int64)   ## wasting 63 bits
0043     double whereAmI ;     // wai = fs[:,3,1].astype(np.int64)   ## wasting 62-63 bits
0044     double c ;            // c   = fs[:,3,2]
0045     double PhotonId ;     // pid = fs[:,3,3].astype(np.int64)   ## again wasting bits
0046
0047     // NB for python parsing check line terminations with set list
0048     // TODO: improve py parsing to cope with comments
0049 };
0050