HitProfile.cxx

Go to the documentation of this file.

#include "Simulation/HitProfile.h"

using namespace sim;
ClassImp(HitProfile);

HitProfile::HitProfile() :
  fTrackId(0), fFraction(0), kNoiseId(65535)
{
}

//__________________________________________________________

HitProfile::~HitProfile() {
}

//__________________________________________________________

unsigned int HitProfile::NComponents() const {

  char msg[256];
  sprintf(msg,"HitProfile::NComponents() inconsistent size!");
  assert_jobc(fTrackId.size()==fFraction.size(),msg);
  
  return fTrackId.size();
}

//__________________________________________________________

std::pair<unsigned short, float> HitProfile::Component(unsigned int i) const
{
  char msg[256];
  sprintf(msg,"HitProfile::Component(%d) out of range!",i);
  assert_jobc(i<fTrackId.size(),msg);

  std::pair<unsigned short, float> component;
  component.first  = fTrackId[i];
  component.second = fFraction[i];

  return component;
}

//__________________________________________________________

std::pair<unsigned short, float> HitProfile::Component(unsigned int i, float fact) const
{
  char msg[256];
  sprintf(msg,"HitProfile::Component(%d) out of range!",i);
  assert_jobc(i<fTrackId.size(),msg);

  std::pair<unsigned short, float> component;
  component.first  = fTrackId[i];
  component.second = fFraction[i] * fact;

  return component;
}

//__________________________________________________________

unsigned short HitProfile::NoiseId()   const
{
  return kNoiseId;

}

//__________________________________________________________

unsigned short HitProfile::TrackId(unsigned int i)   const
{
  char msg[256];
  sprintf(msg,"HitProfile::TrackId(%d) out of range!",i);
  assert_jobc(i<fTrackId.size(),msg);

  return fTrackId[i];

}

//__________________________________________________________

float HitProfile::Fraction(unsigned int i)  const
{
  char msg[256];
  sprintf(msg,"HitProfile::Fraction(%d) out of range!",i);
  assert_jobc(i<fFraction.size(),msg);

  return fFraction[i];
}

//__________________________________________________________

void HitProfile::Push(std::pair<unsigned short, float> component)
{
  fTrackId.push_back(component.first);
  fFraction.push_back(component.second);
  return;
}

//__________________________________________________________

void HitProfile::Push(unsigned short trackid, float fraction)
{
  fTrackId.push_back(trackid);
  fFraction.push_back(fraction);
  return;
}

//__________________________________________________________

void HitProfile::SetFraction(unsigned int i, float fraction)
{
  char msg[256];
  sprintf(msg,"HitProfile::SetFraction(%d) out of range!",i);
  assert_jobc(i<fFraction.size(),msg);

  fFraction[i] = fraction;

  return;

}

//__________________________________________________________

void HitProfile::NormalizeFraction()
{
  if(this->NComponents()==0) return;

  float total = 0.;
  for(unsigned int j=0; j<this->NComponents(); j++) 
    total += fFraction[j];

  if(total==0) return;

  for(unsigned int j=0; j<this->NComponents(); j++) 
    fFraction[j] /= total;

  return;
}
    
//__________________________________________________________

void HitProfile::Combine()
{
  std::vector<unsigned short> copyTrackId = fTrackId;
  std::vector<float>          copyFraction = fFraction;
  
  fTrackId.clear();
  fFraction.clear();

  for(unsigned int i=0; i<copyTrackId.size(); i++) {

    bool duplicated = false;

    for(unsigned int j=0; j<fTrackId.size(); j++) {

      if(copyTrackId[i]==fTrackId[j]) {
        fFraction[j] += copyFraction[i];
        duplicated = true;
        continue;
      }
    }

    if(!duplicated) {
      fTrackId.push_back(copyTrackId[i]);
      fFraction.push_back(copyFraction[i]);
    }

  }

  return;
}

//__________________________________________________________

void HitProfile::Clear()
{
  fTrackId.clear();
  fFraction.clear();
  return;
}

---