---

CoreDepositor.cxx

Go to the documentation of this file.

#include "FastShowerUtils/CoreDepositor.h"
#include "FastShowerUtils/ISampler.h"
#include "FastShowerUtils/Deposits.h"
#include "FastShowerUtils/ParticleParameters.h"
#include "FastShowerUtils/PolyArgs.h"

#include "FastShowerUtils/Normalisations.h"

#include <boost/array.hpp>
#include <algorithm>

#include <pair.h>


#include <iostream>
#include <iomanip>
#include <numeric>
#include <assert.h>

namespace FastShower{
  using boost::array;
  using std::pair;
  //
  //                                               //
  //Function objects used by STL algorithms        //
  //                                               //
  class ProxToGeom{
    public:

    ProxToGeom(CoreDepositor::ProxDepIter pdIter, IDeposits& deposits):
      m_proxIter(pdIter), m_deposits(deposits){}
    
    void operator()(std::pair<double, double> geom){
      m_deposits.accept(geom.first, geom.second, *m_proxIter++);
    }
    
  private:
    CoreDepositor::ProxDepIter m_proxIter;
    IDeposits& m_deposits;
  };
  
  
  class ApplyNorm{
  public:
    ApplyNorm(double norm):m_norm(norm){}
    void operator()(double& proxDep){ 
      proxDep *= m_norm;
    }
  private:
    double m_norm;
  };
  //                                               //
  //   Function objects used by STL algorithms     //
  //                                               //


  CoreDepositor::CoreDepositor(ISampler* s0,ISampler* c0,
                               ISampler* csn,ISampler* can):
    IDepositor(), DebugBase("CoreDepositor"){
    unsigned int it = 0;
    unsigned int end = m_samplers.size();

    m_samplers[it++]=s0;
    m_samplers[it++]=c0;
    m_samplers[it++]=csn;
    m_samplers[it++]=can;
    assert(it=end);
    
    m_proximityDeps.assign(0.); //eta and phi offsets wrt to cell 0 center
    array<pair<int, int>, 6> q0, q1, q2, q3;
    // 
    q0[sf]=pair<int, int>(  0, -1);
    q0[sc]=pair<int, int>(  0,  0);
    q0[sn]=pair<int, int>(  0,  1);
    q0[af]=pair<int, int>(  1, -1);
    q0[ac]=pair<int, int>(  1,  0);
    q0[an]=pair<int, int>(  1,  1);
    m_mappers[0]=q0;
    //
    q1[sf]=pair<int, int>(  0, -1);
    q1[sc]=pair<int, int>(  0,  0 );
    q1[sn]=pair<int, int>(  0,  1);
    q1[af]=pair<int, int>( -1, -1);
    q1[ac]=pair<int, int>( -1,  0 );
    q1[an]=pair<int, int>( -1,  1);
    m_mappers[1]=q1;
    //
    q2[sf]=pair<int, int>(  0,  1);
    q2[sc]=pair<int, int>(  0,  0 );
    q2[sn]=pair<int, int>(  0, -1);
    q2[af]=pair<int, int>( -1,  1);
    q2[ac]=pair<int, int>( -1,  0);
    q2[an]=pair<int, int>( -1, -1);
    m_mappers[2]=q2;
    //
    q3[sf]=pair<int, int>(  0,  1);
    q3[sc]=pair<int, int>(  0,  0);
    q3[sn]=pair<int, int>(  0, -1);
    q3[af]=pair<int, int>(  1,  1);
    q3[ac]=pair<int, int>(  1,  0);
    q3[an]=pair<int, int>(  1, -1);
    m_mappers[3]=q3;

  }

  //
  void CoreDepositor::deposit(const ParticleParameters& pp,
                              const Normalisations& ns,
                              IDeposits& geomDeposits){
    unsigned int ind = 0;
    unsigned int end = m_samplers.size();
    PolyArgs pa(pp, m_cs);
    
    // fill core samples
    for(; ind!=end; ++ind){
      m_samplers[ind]->sample(pa, m_cs);
    }

    // fill proximity map
    this->fillProximityDeps();

    //normalise core deposits
    std::for_each(m_proximityDeps.begin(),
          m_proximityDeps.end(),
          ApplyNorm(ns.core())
          );

    // proximity to geometry map
    std::for_each(
                  m_mappers[pp.quadrant()].begin(),
                  m_mappers[pp.quadrant()].end(),
                  ProxToGeom( m_proximityDeps.begin(), geomDeposits )
                  );

  }

  //
  void CoreDepositor::fillProximityDeps(){
    double sfF, acF, afF;
    sfF = 1.0 - m_cs.cell0() - m_cs.cellSN();
    assert(sfF>=0.0 && sfF<=1.0);
    if (sfF<=0.0) {
      afF = 0.0;
    } else {
      afF = (1.0 - m_cs.cellAN())/(1.0 + m_cs.cell0()/sfF);
    }
    assert(afF>=0.0 && afF<=1.0);
    if (m_cs.cell0()<=0.0) {
      acF = 0.0;
    } else {
      acF = (1.0 - m_cs.cellAN())/(1.0 + sfF/m_cs.cell0());
    }
    assert(acF>=0.0 && acF<=1.0);

    m_proximityDeps[sc] = m_cs.cell0()*m_cs.slice0();
    m_proximityDeps[sn] = m_cs.cellSN()*m_cs.slice0();
    m_proximityDeps[sf] = sfF*m_cs.slice0();
    m_proximityDeps[an] = m_cs.cellAN()*(1.0-m_cs.slice0());
    m_proximityDeps[af] = afF*(1.0-m_cs.slice0());
    m_proximityDeps[ac] = acF*(1.0-m_cs.slice0());
  }

  //
  IDepositor* CoreDepositor::clone() const {
    IDepositor* d = new CoreDepositor(*this);
    return d;
  }

  //
  void CoreDepositor::components(std::vector<const IDebug*>& v) const{
    v.clear();
    std::copy(m_samplers.begin(), m_samplers.end(), back_inserter(v));
  }
  
}//namespace

---