---

ShowerDemoDumper.cxx

Go to the documentation of this file.

#include "FastShowerUtils/ShowerDemoDumper.h"
#include "FastShowerUtils/Gridlet.h"
#include "FastShowerUtils/GridletForger.h"
#include "FastShowerUtils/IConfigurer.h"

#include "FastShowerUtils/ParticleInfo.h"
#include "FastShowerUtils/Pinger/ComplexPinger.h"
#include "FastShowerUtils/Pinger/ComplexArmer.h"
#include "FastShowerUtils/Pinger/ComplexStatsCollector.h"
#include "FastShowerUtils/ProcessedFlat.h"
#include "FastShowerUtils/LinearProcessor.h"

#include<cmath>
namespace FastShower{
  ShowerDemoDumper::ShowerDemoDumper(int nEvents, 
                                     std::string key, int pdgId,
                                     double minEta, double maxEta,
                                     double energy):
    m_nEvents(nEvents), m_pdgId(pdgId), 
    m_minEta(minEta), m_maxEta(maxEta), m_energy(energy){
    cout<<"ShowerDemoDumper(int nEvents, string key, int pdgId, "
        <<"double minEta, double maxEta, double energy): "<<m_nEvents<<" "
        <<m_key<<" "<<m_pdgId<<" "<<m_minEta<<" "<<m_maxEta<<" "<<m_energy<<endl;
  }


  void ShowerDemoDumper::initialise(const std::string s){
    cout<<"ShowerDemoDumper - initialising"<<endl;
    m_gridletForger = new GridletForger(s);
    cout<<"*****************************************"<<endl;
    cout<<"*       GridletForger Structure         *"<<endl;
    cout<<"*****************************************"<<endl;
    ComplexPinger p;
    IDebug* gfDbg = gridletForger();
    gfDbg->ping(p);
  } 


  void ShowerDemoDumper::execute(){
    std::string eFile("Test.txt");
    std::fstream ofile;
    ofile.open(eFile.c_str(), ios::out);
    cout<<"*****************************************"<<endl;
    cout<<"*       ShowerDemoDumper Execute        *"<<endl;
    cout<<"*****************************************"<<endl;

    for(int ievt = 0; ievt!=m_nEvents; ++ievt){
      cout<<"DEMO  Event: "<<ievt<<endl<<endl;
      ofile<<"DEMO  Event: "<<ievt<<endl<<endl;
      //
      ParticleInfo pi = makePI(m_pdgId, m_energy);
      ofile<<pi;
      //
      Gridlet* g = gridletForger()->makeGridlet(pi);
      cout<<*g;
      ofile<<*g;
      delete g;
      //reset what needs to be reset before using the showerer;
      //ComplexArmer ca;
      //IDebug* gfDbg = gridletForger();
      //gfDbg->ping( ca );
      //
      //      ComplexStatsCollector sc;
      //      gfDbg->ping( sc );
      //
    }
  }


  void ShowerDemoDumper::finalise() const{
    cout<<"ShowerDemoDumper: Finalise "<<endl;
    ComplexStatsCollector sc;
    IDebug* gfDbg = gridletForger();
    gfDbg->ping( sc );
  }


  GridletForger* ShowerDemoDumper::gridletForger() const{
    return m_gridletForger;
  }


  ParticleInfo ShowerDemoDumper::makePI(int pdgId, double energy) const {
    ProcessedFlat pf;
    LinearProcessor etaLp( (m_maxEta-m_minEta), m_minEta );
    LinearProcessor phiLp( 2*M_PI, -M_PI );
    double eta = pf.sample(&etaLp);
    double phi = pf.sample(&phiLp);
    ParticleInfo pi(phi,eta,energy,pdgId);
    return pi;
  }
}//namespace

---