---

GridletForger.cxx

Go to the documentation of this file.

#include "FastShowerUtils/GridletForger.h"

#include "FastShowerUtils/SingleShowererSelector.h"
#include "FastShowerUtils/DoubleShowererSelector.h"
#include "FastShowerUtils/TripleShowererSelector.h"

#include "FastShowerUtils/ConfigurerBase.h"

#include "FastShowerUtils/SoftPhotonBarSelConfig.h"
#include "FastShowerUtils/SoftPhotonEcSelConfig.h"
#include "FastShowerUtils/ElectronBarSelConfig.h"
#include "FastShowerUtils/ElectronEcSelConfig.h"
#include "FastShowerUtils/HardPhotonBarSelConfig.h"
#include "FastShowerUtils/HardPhotonEcSelConfig.h"
#include "FastShowerUtils/HadBarSelConfig.h"
#include "FastShowerUtils/HadEcSelConfig.h"


#include <cmath>
#include <assert.h>
#include <iostream>
#include <algorithm>
namespace FastShower{

  const double minPhi = -M_PI;
  const double maxPhi = +M_PI;
  const double minEta = -3.2;
  const double maxEta = +3.2;
  const double phiGran = M_PI/32;  // 0.09817
  const double etaGran = 3.2/32;   // 0.1

  const double GridletForger::phiStart(-2.0);
  const double GridletForger::phiStep(0.1); // 0.098
  const double GridletForger::etaStart(-2.0);
  const double GridletForger::etaStep(0.1);

  GridletForger::GridletForger(std::vector<IShowererSelector*> v):
    IDebug(), DebugBase("GridletForger"), m_showererSelectors(v){}

  GridletForger::GridletForger(): IDebug(), DebugBase("GridletForger"){
    Moni moniMode;

    SoftPhotonBarSelConfig spbConfig(moniMode);
    m_showererSelectors.push_back(new SingleShowererSelector(spbConfig));

    SoftPhotonEcSelConfig specConfig(moniMode);
    m_showererSelectors.push_back(new SingleShowererSelector(specConfig));

    HardPhotonBarSelConfig hpbConfig(moniMode);
    m_showererSelectors.push_back(new SingleShowererSelector(hpbConfig));

    HardPhotonEcSelConfig hpecConfig(moniMode);
    m_showererSelectors.push_back(new SingleShowererSelector(hpecConfig));

    ElectronBarSelConfig ebConfig(moniMode);
    m_showererSelectors.push_back(new SingleShowererSelector(ebConfig));

    ElectronEcSelConfig eecConfig(moniMode);
    m_showererSelectors.push_back(new SingleShowererSelector(eecConfig));

    //HadBarSelConfig hbConfig(moniMode);
    //m_showererSelectors.push_back(new DoubleShowererSelector(hbConfig));

    //HadEcSelConfig hecConfig(moniMode);
    //m_showererSelectors.push_back(new DoubleShowererSelector(hecConfig));

    HadBarSelConfig hbConfig(moniMode);
    m_showererSelectors.push_back(new TripleShowererSelector(hbConfig));

    HadEcSelConfig hecConfig(moniMode);
    m_showererSelectors.push_back(new TripleShowererSelector(hecConfig));
  }


  GridletForger::GridletForger(const std::string s): IDebug(), DebugBase("GridletForger"){
    Moni moniMode;

    map<std::string, int> key_map;
    // eta scan for Bar or Ec
    key_map[ "spb"  ] = 1;
    key_map[ "spec" ] = 2;
    key_map[ "hpb"  ] = 3;
    key_map[ "hpec" ] = 4;
    key_map[ "eb"   ] = 5;
    key_map[ "eec"  ] = 6;
    key_map[ "hb"   ] = 7;
    key_map[ "hec"  ] = 8;
    // eta scan for Bar and Ec
    key_map[ "SoftPhoton"  ] = 12;
    key_map[ "HardPhoton"  ] = 34;
    key_map[ "Electron"    ] = 56;
    key_map[ "Hadron"      ] = 78;

    switch (key_map[s])
      {
      case 1:
        {
          SoftPhotonBarSelConfig spbConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(spbConfig));
        }
        break;
      case 2:
        {
          SoftPhotonEcSelConfig specConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(specConfig));
        }
        break;
      case 3:
        {
          HardPhotonBarSelConfig hpbConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(hpbConfig));
        }
        break;
      case 4:
        {
          HardPhotonEcSelConfig hpecConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(hpecConfig));
        }
        break;
      case 5:
        {
          ElectronBarSelConfig ebConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(ebConfig));
        }
        break;
      case 6:
        {
          ElectronEcSelConfig eecConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(eecConfig));
        }
        break;
      case 7:
        {
          HadBarSelConfig hbConfig(moniMode);
          m_showererSelectors.push_back(new TripleShowererSelector(hbConfig));
        }
        break;
      case 8:
        {
          HadEcSelConfig hecConfig(moniMode);
          m_showererSelectors.push_back(new TripleShowererSelector(hecConfig));
        }
        break;
      case 12:
        {
          SoftPhotonBarSelConfig spbConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(spbConfig));
          SoftPhotonEcSelConfig specConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(specConfig));
        }
        break;
      case 34:
        {
          HardPhotonBarSelConfig hpbConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(hpbConfig));
          HardPhotonEcSelConfig hpecConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(hpecConfig));
        }
        break;
      case 56:
        {
          ElectronBarSelConfig ebConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(ebConfig));
          ElectronEcSelConfig eecConfig(moniMode);
          m_showererSelectors.push_back(new SingleShowererSelector(eecConfig));
        }
        break;
      case 78:
        {
          HadBarSelConfig hbConfig(moniMode);
          m_showererSelectors.push_back(new TripleShowererSelector(hbConfig));
          HadEcSelConfig hecConfig(moniMode);
          m_showererSelectors.push_back(new TripleShowererSelector(hecConfig));
        }
        break;
      default:
        cout<<"Invalid Showerer-type request: "<<s<<endl;
        throw;
      }
    
  }


  //
  GridletForger::~GridletForger(){
    std::vector<IShowererSelector*>::iterator it = m_showererSelectors.begin();
    std::vector<IShowererSelector*>::iterator end = m_showererSelectors.end();
    for(; it!= end; ++it){delete (*it);}
  }


  Gridlet* GridletForger::makeGridlet(ParticleInfo& pi) const{

    ParticleParameters pp = makeParameters(pi);
    //    cout<<"GridletForger "<<endl;
    //    cout<<pi;
    //    cout<<pp;

    // gridlwt with center coordinate of the hit cell
    Gridlet* g = new Gridlet(pp.eta()-pp.delEta(),pp.phi()-pp.delPhi());

    if( find_if(m_showererSelectors.begin(),m_showererSelectors.end(),selectShowerer(pp,g))
        == m_showererSelectors.end() ){
      g=0;
    }
    
    return g;
  };

  //
  //  ParticleParameters 
  //  GridletForger::makeParameters(const TransportedParticleCollectionCIter& tpcCIter){
  ParticleParameters 
  GridletForger::makeParameters(const ParticleInfo& pi) const{

    // calo region?
    ParticleParameters::Region region;
    if (fabs(pi.eta())<=1.5)
      region = ParticleParameters::barrel;
    else if (fabs(pi.eta())>1.5 && fabs(pi.eta())<=maxEta)
      region = ParticleParameters::endcap;
    else if (fabs(pi.eta())>maxEta && fabs(pi.eta())<5.0)
      region=ParticleParameters::fcal;
    else
      region=ParticleParameters::beyond;

    // hit cell indices
    int    iPhi = static_cast<int>((pi.phi()-minPhi)/phiGran);
    int    iEta = static_cast<int>((pi.eta()-minEta)/etaGran);
    // center coordinate of the hit cell
    double cPhi = minPhi + (iPhi+0.5)*phiGran;
    double cEta = minEta + (iEta+0.5)*etaGran;
    // distance: impact point to cell center (assuming symmetry about Eta=0!)
    double delPhi = pi.phi() - cPhi;
    double delEta = pi.eta() - cEta;

    ParticleParameters 
      pp(pi.phi(),pi.eta(),delPhi,delEta,pi.energy(),region,pi.id());

    EnergyLimiter* energyLimiter;
    if (pp.region()==ParticleParameters::endcap && pp.isHadron()){
      energyLimiter = new EnergyLimiter(10., 50., "EnergyLimiter");
    }else{
      energyLimiter = new EnergyLimiter(10., 100., "EnergyLimiter");
    }
    pp.setEnergyLimiter(energyLimiter);
    
    return pp;
  }
  //
  void GridletForger::components(IDebug::Cpts& c) const{
    std::copy(m_showererSelectors.begin(), m_showererSelectors.end(), back_inserter(c));
  }

}//namespace

---