---

RTestParameterGenerator.cxx

Go to the documentation of this file.

#include "FastShowerUtils/RTestParameterGenerator.h"

#include <cmath>
namespace FastShower{
  
  RTestParameterGenerator::RTestParameterGenerator(int nEvents):
    m_nEvents(nEvents), m_curEvnt(nEvents), m_nParticles(0){
    
   cout<<"RTestParameterGenerator begin"<<endl;

   double delEnergy = 0.0;   
   double energy    = 0.0;
   // generate energies: 0,1,3,6,10,15,21,28,36,45,55,66,78,91,105,120
   while(energy < 120.0){
     energy += delEnergy;
     m_energies.push_back(energy);
     ++delEnergy;
   }
   m_enIter = m_energies.begin();
   
   m_ids.push_back(11);
   m_ids.push_back(22);
   m_ids.push_back(211);
   m_pidIter = m_ids.begin();
   
   double eta;
   double cellEtaWidth = 3.2/32;
   double delEta  = cellEtaWidth/10;
   eta = 0.200;
   while(eta<(0.200+cellEtaWidth+0.0005)){
    m_etas.push_back(eta);
    eta += delEta;
   }
   eta = 1.900;
   while(eta<(1.900+cellEtaWidth+0.0005)){
    m_etas.push_back(eta);
    eta += delEta;
   }
   m_etaIter = m_etas.begin();

   double cellPhiWidth = M_PI/32;
   double delPhi  = cellPhiWidth/10;
   double phi = 0.0;
   while(phi<(cellPhiWidth+0.0005)){
    m_phis.push_back(phi);
    phi += delPhi;
   }
   m_phiIter = m_phis.begin();
   

   std::ostream_iterator<double>                     dOut(cout, "\n");
   std::ostream_iterator<int>                        iOut(cout, "\n");
   std::ostream_iterator<double>                     eOut(cout, "\n");
   std::ostream_iterator<double>                     pOut(cout, "\n");
   
   cout<<"RTestParameterGenerator vectors:"<<endl;
   
   cout<<"Energies"<<endl;
   std::copy(m_energies.begin(), m_energies.end(), dOut);
   
   cout<<"Ids"<<endl;
   std::copy(m_ids.begin(), m_ids.end(), iOut);
   
   cout<<"Etas"<<endl;
   std::copy(m_etas.begin(), m_etas.end(), eOut);
   
   cout<<"Phis"<<endl;
   std::copy(m_phis.begin(), m_phis.end(), pOut);
   
   cout<<"RTestParameterGenerator Constructor end"<<endl;
  }
  ParticleInfo* RTestParameterGenerator::next(){
    if(m_curEvnt == 0){
      m_curEvnt = m_nEvents;
      m_nParticles = 0;
      return 0;
    }
    
    ParticleInfo* p = new ParticleInfo(*m_phiIter, *m_etaIter, *m_enIter, *m_pidIter);
    
    if(++m_enIter  == m_energies.end()){
      m_enIter  = m_energies.begin();
      ++m_etaIter;
      if(m_etaIter == m_etas.end()){
        m_etaIter = m_etas.begin();
        ++m_phiIter;
        if(m_phiIter == m_phis.end()){
          m_phiIter = m_phis.begin();
          ++m_pidIter;
          if(m_pidIter == m_ids.end()){
            m_pidIter = m_ids.begin();
            --m_curEvnt;
          }
        }
      }
    }
    
    ++m_nParticles;
    //    cout<<m_nParticles<<" Generated in "<<m_nEvents<<" Events"<<endl;
    
    return p;
  }
  void RTestParameterGenerator::report(){
    cout<<m_nParticles<<" Generated in "<<m_nEvents<<" Events"<<endl;
  }
    
}//namespace
    
    
  
  

---