Feature/adding blip to caf

sbnobj/SBND/Blip/BlipDataTypes.cc

@@ -0,0 +1 @@
+#include "sbnobj/SBND/Blip/BlipDataTypes.h"

sbnobj/SBND/Blip/BlipDataTypes.h

@@ -0,0 +1,174 @@
+#ifndef BLIPDATATYPE
+#define BLIPDATATYPE
+#include "lardataobj/RecoBase/Hit.h"
+#include "nusimdata/SimulationBase/MCParticle.h"
+#include <vector>
+#include <map>
+#include <set>
+typedef std::vector<int>        vint_t;
+typedef std::vector<bool>       vbool_t;
+typedef std::vector<float>      vfloat_t;
+typedef std::set<int>           si_t;
+typedef std::map<int,float>     mif_t;
+
+namespace blip {
+  const int kNplanes    = 3;
+  const int kNTPCs      = 2;
+
+
+  //###################################################
+  //  Data structures
+  //###################################################
+
+  struct ParticleInfo {
+    simb::MCParticle particle;
+    int   trackId           = -9;
+    int   index             = -9;
+    int   isPrimary         = -9;
+    int   numTrajPts          = -9;
+    double depEnergy         = -9;
+    int   depElectrons        = -9;
+    double numElectrons      = -9;
+    double mass              = -9;
+    double E                 = -9;
+    double endE              = -9;
+    double KE                = -9;
+    double endKE             = -9;
+    double P                 = -9; 
+    double Px                = -9; 
+    double Py                = -9; 
+    double Pz                = -9; 
+    double pathLength        = -9;
+    double time              = -9;
+    double endtime           = -9;
+    TVector3 startPoint;      
+    TVector3 endPoint;
+    TVector3 position;
+  };
+
+  // True energy depositions
+  struct TrueBlip {
+    int       ID            = -9;     // unique blip ID
+    int       Cryostat      = -9;     // Cryostat ID
+    int       TPC           = -9;     // TPC ID
+    float     Time          = -9;     // time of particle interaction
+    int       TimeTick      = -9;     // time tick
+    float     DriftTime     = -9;     // drift time [us]
+    float     Energy        = 0;      // energy dep [MeV]
+    int       DepElectrons  = 0;      // deposited electrons
+    int       NumElectrons  = 0;      // electrons reaching wires
+    int       LeadG4ID      = -9;     // lead G4 track ID
+    int       LeadG4Index   = -9;     // lead G4 track index
+    int       LeadG4PDG     = -9;     // lead G4 PDG
+    float     LeadCharge    = -9;     // lead G4 charge dep
+    TVector3  Position;               // XYZ position
+    mif_t     G4ChargeMap;          
+    mif_t     G4PDGMap;
+  };
+
+  struct HitInfo {
+    int   hitid         = -9;
+    int   cryo          = -9;
+    int   tpc           = -9;
+    int   plane         = -9;
+    int   wire          = -9;
+    int   chan          = -9;
+    float amp           = -9;
+    float rms           = -9;
+    int   trkid         = -9;
+    int   shwrid        = -9;
+    int   clustid       = -9;
+    int   blipid        = -9;
+    bool  ismatch       = false;
+    float integralADC    = -999;     // [ADCs] from integral
+    float sigmaintegral = -999;
+    float sumADC        = -999;     // [ADCs] from sum 
+    float charge        = -999;     // [e-]
+    float peakTime      = -999999;
+    float driftTime     = -999999;  // [tick]
+    float gof           = -9;
+    int   g4trkid       = -9;
+    int   g4pdg         = -999;
+    int   g4charge      = -999;     // [e-]
+    float g4frac        = -99;      
+    float g4energy      = -999;     // [MeV]
+  };
+
+  struct HitClust {
+    int     ID              = -9;
+    bool    isValid         = false;
+    int     CenterChan      = -999;
+    int     CenterWire      = -999;
+    bool    isTruthMatched  = false;
+    bool    isMerged        = false;
+    bool    isMatched       = false;
+    int     DeadWireSep     = 99;
+    int     Cryostat        = -9;
+    int     TPC             = -9;
+    int     Plane           = -9;
+    int     NHits           = -9;
+    int     NWires          = -9;
+    float   ADCs            = -999;
+    float   Amplitude       = -999;
+    float   Charge          = -999;
+    float   SigmaCharge     = -999;
+    float   TimeTick        = -999;
+    float   Time            = -999;
+    float   StartHitTime    = -999;
+    float   EndHitTime      = -999;
+    float   StartTime       = -999;
+    float   EndTime         = -999;
+    float   Timespan        = -999;
+    float   RMS             = -999;
+    int     StartWire       = -999;
+    int     EndWire         = -999;
+    int     NPulseTrainHits = -9;
+    float   GoodnessOfFit   = -999;
+    int     BlipID          = -9;
+    int     EdepID          = -9;
+    si_t    HitIDs;
+    si_t    Wires;
+    si_t    Chans;
+    si_t    G4IDs;
+
+    std::map<int,TVector3> IntersectLocations;
+  };
+
+  struct Blip {
+
+    int       ID              = -9;         // Blip ID / index
+    bool      isValid         = false;      // Blip passes basic checks
+    int       Cryostat        = -9;         // Cryostat
+    int       TPC             = -9;         // TPC
+    int       NPlanes         = -9;         // Num. matched planes
+    int       MaxWireSpan     = -9;         // Maximum span of wires on any plane cluster
+    float     TimeTick        = -999;       // Readout time [ticks]
+    float     Time            = -999;       // Drift time [us]
+    float     Charge          = -9;         // Charge on calorimetry plane
+    float     Energy          = -999;       // Energy (const dE/dx, fcl-configurable)
+    float     EnergyESTAR     = -999;       // Energy (ESTAR method from ArgoNeuT)
+    float     EnergyPSTAR     = -999;       // Energy (PSTAR method similar with ESTAR method from ArgoNeuT)
+    float     ProxTrkDist     = -9;         // Distance to cloest track
+    int       ProxTrkID       = -9;         // ID of closest track
+    bool      inCylinder      = false;      // Is it in a cone/cylinder region? 
+
+    TVector3  Position;                     // 3D position TVector3
+    float     SigmaYZ         = -9.;        // Uncertainty in YZ intersect [cm]
+    float     dX              = -9;         // Equivalent length along drift direction [cm] 
+    float     dYZ             = -9;         // Approximate length scale in YZ space [cm]
+
+    // Plane/cluster-specific information
+    blip::HitClust clusters[kNplanes];
+
+    // Truth-matched energy deposition
+    blip::TrueBlip truth;
+
+    // Prototype getter functions
+    double X() { return Position.X(); }
+    double Y() { return Position.Y(); }
+    double Z() { return Position.Z(); }
+
+  };
+
+}
+#endif

sbnobj/SBND/Blip/CMakeLists.txt

@@ -0,0 +1,11 @@
+cet_make_library( LIBRARY_NAME sbndobj_BlipDataTypes
+  SOURCE
+    BlipDataTypes.cc
+  LIBRARIES
+    lardataobj::RecoBase
+    nusimdata::SimulationBase
+)
+art_dictionary(DICTIONARY_LIBRARIES sbndobj_BlipDataTypes) 
+install_headers()
+install_source()
+install_fhicl()

sbnobj/SBND/Blip/classes.h

@@ -0,0 +1,36 @@
+//
+// Build a dictionary.
+//
+// $Id: classes.h,v 1.8 2010/04/12 18:12:28  Exp $
+// $Author:  $
+// $Date: 2010/04/12 18:12:28 $
+// 
+// Original author Rob Kutschke, modified by wes
+//
+
+#include "canvas/Persistency/Common/Wrapper.h"
+
+// data-products
+// lardataobj
+//#include "lardata/Utilities/AssociationUtil.h"
+#include "canvas/Persistency/Common/Assns.h"
+#include "lardataobj/RecoBase/PFParticle.h"
+#include "lardataobj/RecoBase/Hit.h"
+#include "nusimdata/SimulationBase/MCTruth.h"
+#include "sbnobj/SBND/Blip/BlipDataTypes.h"
+#include "lardataobj/RecoBase/SpacePoint.h"
+
+//
+// Only include objects that we would like to be able to put into the event.
+// Do not include the objects they contain internally.
+//
+
+template class art::Assns<recob::Hit,blip::Blip,void>;
+template class art::Wrapper<art::Assns<recob::Hit,blip::Blip,void> >;
+template class std::vector<blip::Blip>;
+template class art::Wrapper<std::vector<blip::Blip> >;
+template class std::map<int, TVector3>;
+template class art::Assns<blip::Blip,recob::Hit,void>;
+template class art::Wrapper<art::Assns<blip::Blip,recob::Hit,void> >;
+template class art::Assns<blip::Blip,recob::SpacePoint,void>;
+template class art::Wrapper<art::Assns<blip::Blip,recob::SpacePoint,void> >;

sbnobj/SBND/Blip/classes_def.xml

@@ -0,0 +1,14 @@
+<lcgdict>
+<class name="art::Wrapper<std::vector<blip::Blip> >"/>
+<class name="blip::Blip"/>
+<class name="std::vector<blip::Blip>"/>
+<class name="blip::HitClust"/>
+<class name="blip::TrueBlip"/>
+<class name="std::map<int,TVector3>"/>
+<class name="art::Assns<blip::Blip,recob::Hit,void>"/>
+<class name="art::Assns<recob::Hit, blip::Blip,void>"/>
+<class name="art::Assns<blip::Blip,recob::SpacePoint,void>"/>
+<class name="art::Wrapper<art::Assns<blip::Blip,recob::Hit,void> >"/>
+<class name="art::Wrapper<art::Assns<recob::Hit,blip::Blip,void> >"/>
+<class name="art::Wrapper<art::Assns<blip::Blip,recob::SpacePoint,void> >"/>
+</lcgdict>

sbnobj/SBND/CMakeLists.txt

@@ -3,3 +3,4 @@ add_subdirectory(Commissioning)
 add_subdirectory(Trigger)
 add_subdirectory(ToF)
 add_subdirectory(Timing)
+add_subdirectory(Blip)