Hit.hh

Go to the documentation of this file.

#ifndef HIT_HH_INCLUDED
#define HIT_HH_INCLUDED

#include "TObject.h"
#include "TString.h"
#include "Vec.hh"

struct Hit :
  public TObject 
{  
  int id;                               
  
  // straight from the data
  int dom_id;                          ///< module identifier from the data (unique in the detector).
  unsigned int channel_id;             ///< PMT channel id {0,1, .., 31} local to moduke
  unsigned int tdc;                    ///< hit tdc (=time in ns)
  unsigned int tot;                    ///< tot value as stored in raw data (int for pyroot)
  int trig;                            ///< non-zero if the hit is a trigger hit.

  int pmt_id;                          ///< global PMT identifier as found in evt files

  // values after calibration
  double t;                            ///< hit time (from calibration or MC truth)
  double a;                            ///< hit amplitude (in p.e.)
  Vec pos;                             ///< hit position 
  Vec dir;                             ///< hit direction; i.e. direction of the PMT

  double pure_t;                       ///< photon time before pmt simultion (MC only)
  double pure_a;                       ///< amptitude before pmt simution (MC only)
  int type;                            ///< particle type or parametrisation used for hit (mc only)
  int origin;                          ///< track id of the track that created this hit
  
  unsigned pattern_flags;              ///< some number that you can use to flag the hit
  
  /**
   * Default constructor.
   */
  Hit(): id(0), dom_id(0), channel_id(0), tdc(0), tot(0), trig(0), pmt_id(0), 
         t(0), a(0), pure_t(0), pure_a(0), type(0), origin(0), pattern_flags(0) {} 
  //virtual ~Hit() {}
  

  /**
   * Read hit (useful in python).
   *
   * \param  h             hit
   */
  void read(const Hit& h) { *this = h;}

  /**
   * Write hit (useful in python).
   *
   * \param  h             hit
   */
  void write(Hit& h) const { h = *this;}

  /**
   * Print hit.
   *
   * \param  out           output stream
   */
  void print( std::ostream& out = std::cout ) const
  { 
    out << "Hit: id=" << id << " dom="<<dom_id << " channel="<<channel_id;
    out << " pmt="<<pmt_id << " t=" << t << " tot="<<tot;
    out << " pos="; pos.print(out);
    out << " dir="; dir.print(out);
  } 
 
  ClassDefNV(Hit,105) // reserve <100 for antcc class of the same name
};

#endif