source: trunk/examples/novice/gemc/src/identifier.cc

// %%%%%%%%%%
// G4 headers
// %%%%%%%%%%
#include "G4UnitsTable.hh"

// %%%%%%%%%%%%
// gemc headers
// %%%%%%%%%%%%
#include "detector.h"
#include "identifier.h"

// %%%%%%%%%%%
// C++ headers
// %%%%%%%%%%%
#include <iostream>
using namespace std;


bool identifier::operator == (const identifier& I) const
{
 if(I.name == this->name && I.rule == this->rule && I.id == this->id && fabs(I.time - this->time) <= this->TimeWindow)
    return true;
 else
  return false;
}


bool identifier::operator < (const identifier& I) const
{
 if(this->name == I.name)
    if(this->id < I.id)     return true;
    if(this->time < I.time) return true;

 return false;
}

bool identifier::operator > (const identifier& I) const
{
 if(this->name == I.name)
    if(this->id > I.id)     return true;
    if(this->time > I.time) return true;

 return false;
}


ostream &operator<<(ostream &stream, vector<identifier> Iden)
{
 for(int i=0; i<Iden.size(); i++)
 {
    cout << "      id " << i+1 ;
    cout << "           " ;
    cout.width(10);
    cout << Iden[i].name << " " << Iden[i].id << endl;
 }
 cout << "      identifier time:  " <<  Iden[0].time/ns << " ns - TimeWindow: " << Iden[0].TimeWindow/ns << " ns." << endl;

 return stream;
}


#include <cmath>
vector<int> set_IDshifts(vector<int> maxs)
{
 vector<int> y;
 int  max;
 for(int i=0; i<maxs.size(); i++)
 {
   max = 0;
   while( ( 1 << max )  <  maxs[i] ) max++;
   y.push_back( max );
 }
 return y;
}

// Sets the ncopy ID accordingly to Geant4 Volumes copy number
vector<identifier> SetId(vector<identifier> Iden, G4VTouchable* TH, double time, double TimeWindow)
{
  vector<identifier> identity = Iden;

 // Look for "ncopy" flag, set to volume copy number
 for(int i=0; i<identity.size(); i++)
 {
    if(identity[i].rule.find("ncopy") != string::npos)
    {
       // h=1 don't need to check volume itself
       for(int h=0; h<TH-> GetHistoryDepth(); h++)
       {
          string pname = TH->GetVolume(h)->GetName();
          int    pcopy = TH->GetVolume(h)->GetCopyNo();
          if(pname.find(identity[i].name) != string::npos) identity[i].id = pcopy;
       }
    }
 
    // Make sure id is not still zero
    if(identity[i].id == 0) 
    {
      cout << " Something is wrong. Identity not completely set." << endl;
      cout << identity;
      cout << " Exiting. " << endl;
      exit(0);
    }
    identity[i].time       = time;
    identity[i].TimeWindow = TimeWindow;
 }

 return identity;
}