CalSectionReject.cxx

Go to the documentation of this file.

//=================================================
//
//
// Selector function object used by CalSection to
// make acceptance cuts on particles.
//
// Cuts on:
//
// pseudoRapidity - from CLHEP momentum
// phi            - given separately (eg rotated by B field) 
//
// Author: Peter Sherwood
//
//
//===================================================
#include "AtlfastAlgs/CalSectionReject.h"
#include "AtlfastEvent/ITransportedParticle.h"
#include "AtlfastEvent/EPileupDeposit.h"
#include "FastShowerUtils/Gridlet.h"
#include "GaudiKernel/MsgStream.h"
#include <iomanip>
#include <iostream>
//#include <ios>
//#include <iomanip>
#include <fstream>

namespace Atlfast{
  using std::ios;
  using std::setprecision;
  using std::setw;
  using std::ostream;

  using FastShower::GridletElement;

  CalSectionReject::CalSectionReject(double etaMin, double etaMax,
                                     double phiMin, double phiMax):
    m_etaMin(etaMin), m_etaMax(etaMax),
    m_phiMin(phiMin), m_phiMax(phiMax) {}
  
  CalSectionReject::CalSectionReject(const CalSectionReject& other):
    m_etaMin(other.m_etaMin), m_etaMax(other.m_etaMax),
    m_phiMin(other.m_phiMin), m_phiMax(other.m_phiMax) {}
  bool CalSectionReject::operator()(const ITransportedParticle* p) {
    double eta = p->eta();
    //double eta = ((p->particle())->momentum()).pseudoRapidity();
    ::Phi phi = p->phi();
    if(eta<=m_etaMin || eta > m_etaMax) {
      return true;
    }
    if(phi<=m_phiMin || phi > m_phiMax){
      return true;
    }
    return false; 
  }
  
  bool CalSectionReject::operator()(const EPileupDeposit* p) {
    double eta = p->eta();
    ::Phi phi = p->phi();
    if(eta<=m_etaMin || eta > m_etaMax) return true;
    if(phi<=m_phiMin || phi > m_phiMax) return true;
    return false; 
  }
  
  bool CalSectionReject::operator()(const GridletElement* p) {
    double eta = p->eta();
    ::Phi phi = p->phi();
    if(eta<=m_etaMin || eta > m_etaMax) {return true;}
    if(phi<=m_phiMin || phi > m_phiMax) {return true;}
    return false; 
  }

  bool CalSectionReject::operator()(const Gridlet* g) {
    double eta = g->eta0();
    ::Phi phi = g->phi0();
    if(eta<=m_etaMin || eta > m_etaMax) return true;
    if(phi<=m_phiMin || phi > m_phiMax) return true;
    return false; 
  }
  //debug
  MsgStream& CalSectionReject::writeout(MsgStream& stream) const{
    stream<<setiosflags(ios::fixed);
    stream<<setprecision(3);
    stream<<"Lower limits (phi, eta):"<<endreq;
    stream<< m_phiMin <<" "<<m_etaMin<<endreq;
    stream<<"Upper limit (phi, eta):"<<endreq;
    stream<< m_phiMax <<" "<<m_etaMax<<endreq;
    return stream;
  }
  //debug
  //  ostream& CalSectionReject::writeout(ostream& stream) const{
  //  stream<<setiosflags(ios::fixed);
  //  stream<<setprecision(3);
  //  stream<<"Lower limits (phi, eta):"<<endreq;
  //  stream<< m_phiMin <<" "<<m_etaMin<<endreq;
  //  stream<<"Upper limit (phi, eta):"<<endreq;
  //  stream<< m_phiMax <<" "<<m_etaMax<<endreq;
  //  return stream;
  // }


  MsgStream& operator <<( MsgStream& stream, const CalSectionReject& c) {
    return c.writeout(stream);
  }
  MsgStream& operator <<( MsgStream& stream, const CalSectionReject* const c) {
    return (*c).writeout(stream);
  }

  //  ostream& operator <<( ostream& stream, const CalSectionReject& c) {
  //    return c.writeout(stream);
  //  }
  //  ostream& operator <<( ostream& stream, const CalSectionReject* const c) {
  //    return (*c).writeout(stream);
  //  }



}//namespace

---