---

TryToShower.cxx

Go to the documentation of this file.

#include "AtlfastAlgs/Calorimeter.h" 
#include "FastShowerUtils/ParticleInfo.h" 
#include "FastShowerUtils/Gridlet.h" 
#include "FastShowerUtils/GridletForger.h" 
#include "GaudiKernel/MsgStream.h" 
namespace Atlfast{

  using FastShower::ParticleInfo;
  using FastShower::Gridlet;
  using FastShower::GridletForger;

  TryToShower::TryToShower(Calorimeter* c):m_calorimeter(c){
  }
  bool TryToShower::operator()(const TransportedParticle* tp){
   
    MsgStream& log = m_calorimeter->msgStream();

    const GridletForger* gf = m_calorimeter->gridletForger();
    if(gf == 0) return false;
    
    const HepMC::GenParticle* particle = tp->particle();
    ParticleInfo pi( tp->phi(),
                     particle->momentum().pseudoRapidity(),
                     particle->momentum().e(),
                     particle->pdg_id()
                     );
    Gridlet* gridlet = gf->makeGridlet(pi);

    if(!gridlet) {
      //      log<<MSG::DEBUG<<"Particle refused by Showerer "<<endl;
      return false;
    }
    //    else{
      //      log<<MSG::DEBUG<<"Particle accepted by Showerer : gridlet:"<<endl;
      //      log<<MSG::DEBUG<<gridlet<<endl;
    //    }

    std::vector<GridletElement*> e = gridlet->eElements();
    m_calorimeter->addEelements( e.begin(), e.end() );
    std::vector<GridletElement*> h = gridlet->hElements();
    m_calorimeter->addHelements( h.begin(), h.end() );

    return true;
  }  
}//namespace

---