---

BremRBin.cxx

Go to the documentation of this file.

#include "AtlfastAlgs/BremRBin.h"

namespace Atlfast{

  BremRBin::BremRBin(int randSeed) {
    m_randomEngine = new HepJamesRandom(randSeed);
  }

  BremRBin::BremRBin(const BremRBin& other): m_randomEngine(NULL) {
    *this = other;
  }
  
  BremRBin& BremRBin::operator=(const BremRBin& other) {
    if (this != &other) {
      if (m_randomEngine) delete m_randomEngine;
      m_randomEngine = new HepJamesRandom(*other.m_randomEngine);
      m_ptBins = other.m_ptBins;
    }
    return *this;
  }
  BremRBin::~BremRBin() {
    delete m_randomEngine;
  }



  void BremRBin::addBin(TestValue value, BremPTScaleBin bin) {
    m_ptBins[value] = bin;
  }
  double BremRBin::calculatePScale() const {
    double random = m_randomEngine->flat();
    map<TestValue, BremPTScaleBin>::const_iterator iter = m_ptBins.begin();
    map<TestValue, BremPTScaleBin>::const_iterator end = m_ptBins.end();
 
     for (; iter != end; ++iter) {
        if ( random < (iter->first).value() ) {
          return iter->second.calculatePScale(random);
        }else{
        }
      }
      return 1.0;
  }
  const map<TestValue, BremPTScaleBin>* BremRBin::bins() const {return &m_ptBins;} 
  
  MsgStream& operator<<(MsgStream& s, BremRBin& bin) {
    map<TestValue, BremPTScaleBin>::const_iterator iter = bin.bins()->begin();
    map<TestValue, BremPTScaleBin>::const_iterator end = bin.bins()->end();
    int a=0;
    s << setprecision(4);
    s << "RBin with " << bin.bins()->size() << " pT Bins" << endreq;
    for (; iter!=end; ++iter) {
      if (a==12) { a=0; s << endreq; }
      s << setw(6) << iter->first.value();
      ++a;
    }
    s<< endreq;
    return s;
  }
  
}//namepace

  

---