//=================================================
 //
 //
 // 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 "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) {}
   /**
    * Returns true if particle fails acceptance cuts:
    * pseudoRapidity - from CLHEP momentum
    * phi            - rotated by B field 
    */
   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 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
 
 
 
 
 
 
 

Due to the LXR bug, the updates fail sometimes to remove references to deleted files. The Saturday's full rebuilds fix these problems
This page was automatically generated by the LXR engine.