---

BremEtaBin.cxx

Go to the documentation of this file.

#include "AtlfastAlgs/BremEtaBin.h"

namespace Atlfast{

  BremEtaBin::BremEtaBin(int randSeed) {
      m_randomEngine = new HepJamesRandom(randSeed);
    }
  BremEtaBin::BremEtaBin(const BremEtaBin& other): m_randomEngine(NULL){
    *this = other;
  }

  BremEtaBin& BremEtaBin::operator=(const BremEtaBin& other) {
    if (this != &other) {
      if (m_randomEngine) delete m_randomEngine;
      m_randomEngine = new HepJamesRandom(*other.m_randomEngine);
    map<TestValue, BremRBin*>::const_iterator iter = other.m_rBins.begin();
    map<TestValue, BremRBin*>::const_iterator end = other.m_rBins.end();
    for (;iter != end; ++iter) {
      m_rBins[iter->first] = new BremRBin( *(iter->second) );
    }
    }
    return *this;
  }
  BremEtaBin::~BremEtaBin() {
    map<TestValue, BremRBin*>::iterator iter = m_rBins.begin();
    map<TestValue, BremRBin*>::iterator end = m_rBins.end();
    for (;iter !=end; ++iter) {
      delete iter->second;
    }
    delete m_randomEngine;
  }

  void BremEtaBin::addBin(TestValue bremR, BremRBin* bin) {
      m_rBins[bremR] = bin;
    }

  double BremEtaBin::calculatePScale(double random) const {
      map<TestValue, BremRBin*>::const_iterator iter = m_rBins.begin();
      map<TestValue, BremRBin*>::const_iterator end = m_rBins.end();
      for (;iter != end; ++iter) {
        if (random < (iter->first).value() ) {
          return iter->second->calculatePScale();
        }else{
        }
      }
      return m_rBins.rbegin()->second->calculatePScale();
    }
  

  const map<TestValue, BremRBin*>* BremEtaBin::bins() const {return &m_rBins;}

  MsgStream& operator<<(MsgStream& s, BremEtaBin& bin) {
    map<TestValue, BremRBin*>::const_iterator iter = bin.bins()->begin();
    map<TestValue, BremRBin*>::const_iterator end = bin.bins()->end();
    s << "----------------- New Eta Bin  with " << bin.bins()->size() << " R Bins ----------------" << endreq;
    int i=1;
    for (; iter!=end; ++iter) {
      s << "Bremstrahlung Radius " << i << " = " << (iter->first).value() << endreq;
      s << *(iter->second) << endreq;
      ++i;
    }
    return s;
  }

}//namepace

  

---