The Evt class respresent a Monte Carlo (MC) event as well as an offline event. More...

#include <Evt.hh>

Inheritance diagram for Evt:

Public Member Functions
	Evt ()

void	print (std::ostream &out=std::cout) const
	Print event. More...

double	get_ns ()
	Get UTC time. More...

void	clear ()
	Reset event. More...

void	clear_hits ()
	Clear the hit vectors and all the references to hits in the tracks. More...

double	getMCtime (const Hit &h)
	Get the original Monte-Carlo-time of the hit. More...

Hit *	get_hit_by_id (int req_id, bool use_mc_hits=false)
	Get pointer to hit with given identifier. More...

std::vector< Trk * >	primary_trks ()
	Get all the tracks that are 'primary'. More...

Trk *	neutrino ()
	Get a pointer to the (first) neutrino from the MC track list. More...

Trk *	primary_neutrino ()
	Get a pointer to primary neutrino from the MC track list. More...

Trk *	leading_lepton ()
	Get a pointer to leading lepton from the MC track list. More...

Trk *	get_parent_of (Trk &child)
	Get the (first) parent of the track 'child'. More...

int	idx (const std::string &key) const
	Get index in user data of the item with given key. More...

bool	haveusr (const std::string &key) const
	Check availability of user data of the item with given key. More...

int	idxusr_checked (const std::string &key) const
	Get index in user data of the item with given key. More...

double	getusr (const std::string &key) const
	Get user data item with given key. More...

void	setusr (const std::string &key, double value)
	Set user data item with given key. More...

bool	delusr (const std::string &key)
	Remove (first) user data item with given key. More...

void	clearusr ()
	Clear user data. More...

void	printusr (std::ostream &out=std::cout)
	Print user data (i.e. More...

const char *	__str__ () const
	Get class name. More...

Public Attributes
int	id
	offline event identifier More...

int	det_id
	detector identifier from DAQ More...

int	mc_id
	identifier of the MC event (as found in ascii or antcc file). More...

int	run_id
	DAQ run identifier. More...

int	mc_run_id
	MC run identifier. More...

int	frame_index
	from the raw data More...

ULong64_t	trigger_mask
	trigger mask from raw data (i.e. the trigger bits) More...

ULong64_t	trigger_counter
	trigger counter More...

unsigned int	overlays
	number of overlaying triggered events More...

TTimeStamp	t
	UTC time of the start of the timeslice the event came from. More...

std::vector< Hit >	hits
	list of hits More...

std::vector< Trk >	trks
	list of reconstructed tracks (can be several because of prefits,showers, etc). More...

std::vector< double >	w
	MC: Weights w[0]=w1, w[1]=w2, w[2]]=w3 (see e.g. Tag list) More...

std::vector< double >	w2list
	MC: factors that make up w[1]=w2 (see e.g. Tag list) More...

std::vector< double >	w3list
	MC: atmospheric flux information. More...

double	mc_t
	MC: time of the mc event. More...

std::vector< Hit >	mc_hits
	MC: list of MC truth hits. More...

std::vector< Trk >	mc_trks
	MC: list of MC truth tracks. More...

TString	comment
	user can use this as he/she likes More...

int	index
	user can use this as he/she likes More...

int	flags
	user can use this as he/she likes More...

std::vector< double >	usr
	user data More...

std::vector< std::string >	usr_names
	user keys More...

TObject *	any
	Pointer to "any" user data. More...

Detailed Description

The Evt class respresent a Monte Carlo (MC) event as well as an offline event.

Some data from the online (DAQ) event are copied.

Definition at line 19 of file Evt.hh.

Constructor & Destructor Documentation

Evt::Evt ( )

inline

Definition at line 54 of file Evt.hh.

         :
     id(0), det_id(0), mc_id(0), run_id(0), mc_run_id(0), frame_index(0),
     trigger_mask(0), trigger_counter(0),
     overlays(0), t(0), mc_t(0), index(0), flags(0) {}

Member Function Documentation

void Evt::print ( std::ostream & out = std::cout ) const

inline

Print event.

Parameters

out	output stream

Definition at line 65 of file Evt.hh.

   {
     out << "Evt: id=" << id << 
       " run_id=" << run_id << 
       " #hits=" << hits.size() <<
       " #mc_hits=" << mc_hits.size() <<
       " #trks=" << trks.size() <<
       " #mc_trks=" << mc_trks.size() << std::endl;
   }

double Evt::get_ns ( )

inline

Get UTC time.

Returns: time [ns]

Definition at line 81 of file Evt.hh.

   {
     return this->t.GetSec()*1e9 + this->t.GetNanoSec(); 
   }

void Evt::clear ( )

inline

Reset event.

Definition at line 89 of file Evt.hh.

   {
     *this = Evt();
   }

void Evt::clear_hits ( )

inline

Clear the hit vectors and all the references to hits in the tracks.

Definition at line 97 of file Evt.hh.

   {
     hits.clear();
     mc_hits.clear();
     for (auto& t: trks )    t.hit_ids.clear();
     for (auto& t: mc_trks ) t.hit_ids.clear(); 
   }

double Evt::getMCtime ( const Hit & h )

inline

Get the original Monte-Carlo-time of the hit.

Apply a correction to the hit time from the DAQ, so that the time returned by this fucntion matches the original time used in the Monte Carlo event.

Parameters

h hit

Returns: time [ns]

Definition at line 116 of file Evt.hh.

   {
     return h.t+this->get_ns() - this->mc_t;
   }

Hit* Evt::get_hit_by_id	(	int	req_id,
		bool	use_mc_hits = `false`
	)

inline

Get pointer to hit with given identifier.

This method throws a run-time exception if no hit is available.

Parameters

req_id	hit identifier
use_mc_hits	option to use MC hit

Returns: pointer to hit.

Definition at line 129 of file Evt.hh.

   {
     std::vector<Hit>& col = use_mc_hits? mc_hits : hits;
 
     if ( req_id < int(col.size()) ) // try to see if we are lucky
       {
         if ( col[req_id].id == req_id ) return &(col[id]);
       }
 
     for( auto& h: col ) 
       {
         if ( h.id == req_id ) return &(col[id]);
       }
 
     THROW(Exception, "Could not find hit with id=" << req_id);
   }

std::vector<Trk*> Evt::primary_trks ( )

inline

Get all the tracks that are 'primary'.

Here, primary means the tracks have no parents.

Note that this method only works if MC parent-child relations are availabe.

Returns: list of pointers to primary tracks

Definition at line 155 of file Evt.hh.

   {
     std::vector<Trk*> r;
     for(auto& t: mc_trks )
       {
         if ( t.is_primary() ) r.push_back(&t);
       }
     return r;
   }

Trk* Evt::neutrino ( )

inline

Get a pointer to the (first) neutrino from the MC track list.

Returns: pointers to neutrino (may be NULL)

Definition at line 170 of file Evt.hh.

   {
     for(auto& t: mc_trks )
       {
         if ( t.is_neutrino() ) return &t;
       }
     return nullptr;
   }

Trk* Evt::primary_neutrino ( )

inline

Get a pointer to primary neutrino from the MC track list.

Note that this method only works if MC parent-child relations are availabe.

Returns: pointer to primary neutrino (may be NULL)

Definition at line 186 of file Evt.hh.

   {
     for( auto& t : mc_trks ) 
       {
         if ( t.is_neutrino() and t.is_primary() ) return &t;
       }
     return nullptr;
   }

Trk* Evt::leading_lepton ( )

inline

Get a pointer to leading lepton from the MC track list.

Here, leading means the lepton that has a neutrino as mother.

Note that this method only works if MC parent-child relations are availabe.

Returns: pointer to leadig lepton (may be NULL)

Definition at line 203 of file Evt.hh.

   {
     int nu_id = primary_neutrino()->id;
 
     for(auto& t: mc_trks )
       {
         if ( t.is_lepton()        &&
              t.mother_id == nu_id &&
              !t.is_orphan() ) return &t;
       }
     return nullptr;
   }

Trk* Evt::get_parent_of ( Trk & child )

inline

Get the (first) parent of the track 'child'.

This method throws a run-time exception if no parent is available.

Parameters

child child particle

Returns: pointer to parent

Definition at line 223 of file Evt.hh.

   {  
     for(auto& t: mc_trks )
       {
         if (child.mother_id == t.id ) return &t;
       }
     THROW(Exception, "Failed to find parent of track" << child.id << " mother_id="<< child.mother_id);
   }  

int AAObject::idx ( const std::string & key ) const

inlineinherited

Get index in user data of the item with given key.

Parameters

key key

Returns: index (-1 if key does not exists)

Definition at line 29 of file AAObject.hh.

   {
     auto i = std::find (usr_names.begin(), usr_names.end(), key );
     if (i == usr_names.end()) return -1;
     return i - usr_names.begin();
   }

bool AAObject::haveusr ( const std::string & key ) const

inlineinherited

Check availability of user data of the item with given key.

Parameters

key key

Returns: true if available; else false

Definition at line 42 of file AAObject.hh.

   {
     return idx( key ) >= 0;
   }

int AAObject::idxusr_checked ( const std::string & key ) const

inlineinherited

Get index in user data of the item with given key.

This method throws a run-time exception if no user data are available.

Parameters

key key

Returns: index (-1 if key does not exists)

Definition at line 54 of file AAObject.hh.

   {
     int r = idx( key );
     if (r < 0)
       {
         THROW(Exception, "No user data for key " << key << " in aanet object of type " << this -> ClassName());
       }
     return r;
   }

double AAObject::getusr ( const std::string & key ) const

inlineinherited

Get user data item with given key.

This method throws a run-time exception if no user data are available.

Parameters

key key

Returns: value

Definition at line 72 of file AAObject.hh.

   {
     const int i = idx( key );
 
     if ( i  < 0 )
       {
         THROW(Exception, "No user data for key " << key << " in aanet object of type " << this -> ClassName());
       }
 
     if ( unsigned(i) >= usr.size() )
       {
         THROW(Exception, "Warning: inconsistent user data " << i << " >= " << usr.size());
       }
 
     return usr[i];
   }

void AAObject::setusr	(	const std::string &	key,
		double	value
	)

inlineinherited

Set user data item with given key.

Parameters

key	key
value	value

Definition at line 95 of file AAObject.hh.

   {
     int i = idx( key );
     if (i < 0)
       {
         if ( usr.size() < usr_names.size() )
           {
             // this should not happen, but let's just add empty data
             usr.resize( usr_names.size() );
           }
         else
           {
             // this is possible, add empty ("") names
             usr_names.resize( usr.size() );
           }
 
         usr_names.push_back( key );
         usr.push_back( value );
       }
     else
       {
         usr[i] = value;
       }
   }  

bool AAObject::delusr ( const std::string & key )

inlineinherited

Remove (first) user data item with given key.

Parameters

key key

Returns: true if data have been removed; else false

Definition at line 126 of file AAObject.hh.

   {
     int i = idx( key );
     if ( i < 0 ) return false;
 
     usr.erase ( usr.begin() + i );
     usr_names.erase( usr_names.begin() + i );
     return true;
   }

void AAObject::clearusr ( )

inlineinherited

Clear user data.

Definition at line 139 of file AAObject.hh.

   {
     usr.resize(0);
     usr_names.resize(0);
   }

void AAObject::printusr ( std::ostream & out = std::cout )

inlineinherited

Print user data (i.e.

list of all pairs of keys and values).

Parameters

out	output stream

Definition at line 150 of file AAObject.hh.

   {
     unsigned n = std::max( usr.size(), usr_names.size() );
     
     for (unsigned i=0; i< n ;i++)
       {
         std::string name = "(unnamed)";
         if ( i < usr_names.size() && usr_names[i] != "" ) name = usr_names[i];
         out << i << " \t " << name << " : \t ";
         if ( i < usr.size() ) out << usr[i] << std::endl;
         else                  out << "(none)" << std::endl;
       }
   }

const char* AAObject::__str__ ( ) const

inlineinherited

Get class name.

Returns: class name

Definition at line 175 of file AAObject.hh.

175 { return Class_Name(); }

Member Data Documentation

int Evt::id

offline event identifier

Definition at line 21 of file Evt.hh.

int Evt::det_id

detector identifier from DAQ

Definition at line 22 of file Evt.hh.

int Evt::mc_id

identifier of the MC event (as found in ascii or antcc file).

Definition at line 23 of file Evt.hh.

int Evt::run_id

DAQ run identifier.

Definition at line 25 of file Evt.hh.

int Evt::mc_run_id

MC run identifier.

Definition at line 26 of file Evt.hh.

int Evt::frame_index

from the raw data

Definition at line 28 of file Evt.hh.

ULong64_t Evt::trigger_mask

trigger mask from raw data (i.e. the trigger bits)

Definition at line 29 of file Evt.hh.

ULong64_t Evt::trigger_counter

trigger counter

Definition at line 30 of file Evt.hh.

unsigned int Evt::overlays

number of overlaying triggered events

Definition at line 31 of file Evt.hh.

TTimeStamp Evt::t

UTC time of the start of the timeslice the event came from.

Definition at line 32 of file Evt.hh.

std::vector<Hit> Evt::hits

list of hits

Definition at line 35 of file Evt.hh.

std::vector<Trk> Evt::trks

list of reconstructed tracks (can be several because of prefits,showers, etc).

Definition at line 36 of file Evt.hh.

std::vector<double> Evt::w

MC: Weights w[0]=w1, w[1]=w2, w[2]]=w3 (see e.g. Tag list)

Definition at line 39 of file Evt.hh.

std::vector<double> Evt::w2list

MC: factors that make up w[1]=w2 (see e.g. Tag list)

Definition at line 40 of file Evt.hh.

std::vector<double> Evt::w3list

MC: atmospheric flux information.

Definition at line 41 of file Evt.hh.

double Evt::mc_t

MC: time of the mc event.

Definition at line 43 of file Evt.hh.

std::vector<Hit> Evt::mc_hits

MC: list of MC truth hits.

Definition at line 44 of file Evt.hh.

std::vector<Trk> Evt::mc_trks

MC: list of MC truth tracks.

Definition at line 45 of file Evt.hh.

TString Evt::comment

user can use this as he/she likes

Definition at line 48 of file Evt.hh.

int Evt::index

user can use this as he/she likes

Definition at line 49 of file Evt.hh.

int Evt::flags

user can use this as he/she likes

Definition at line 50 of file Evt.hh.

std::vector<double> AAObject::usr

inherited

user data

Definition at line 20 of file AAObject.hh.

std::vector<std::string> AAObject::usr_names

inherited

user keys

Definition at line 21 of file AAObject.hh.

TObject* AAObject::any

inherited

Pointer to "any" user data.

Definition at line 177 of file AAObject.hh.

The documentation for this struct was generated from the following file:

Evt.hh

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation