LocalPhotonSelector Class Reference

#include <LocalPhotonSelector.h>

Inheritance diagram for LocalPhotonSelector:

Public Member Functions
	LocalPhotonSelector (const std::string &type, const std::string &name, const IInterface *parent)
bool	operator() (CDPhoton &aPhoton) override

Detailed Description

Definition at line 7 of file LocalPhotonSelector.h.

Constructor & Destructor Documentation

LocalPhotonSelector()

LocalPhotonSelector::LocalPhotonSelector	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Definition at line 11 of file LocalPhotonSelector.cxx.

    : base_class( type, name, parent ) {
  // Declare the properties
  declareProperty( "MinEnergy", m_minEnergy = 0.025 );
 
  declareProperty( "InBarrelOrEndcap", m_useBarrelEndcap = true );
  declareProperty( "MaxCosThetaBarrel", m_maxCosThetaBarrel = 0.8 );
  declareProperty( "MinCosThetaEndcap", m_minCosThetaEndcap = 0.86 );
  declareProperty( "MaxCosThetaEndcap", m_maxCosThetaEndcap = 0.92 );
  declareProperty( "MinEndcapEnergy", m_minEndcapEnergy = 0.050 );
 
  declareProperty( "ApplyTimeCut", m_applyTimeCut = true );
  declareProperty( "MinTime", m_minTime = 0. );
  declareProperty( "MaxTime", m_maxTime = 14. );
  declareProperty( "PhotonDeltaTime", m_deltaTime = 10. );
 
  declareProperty( "ApplyDangCut", m_applyDangCut = false );
  declareProperty( "MinDang", m_minDang = 20.0 );
}

Referenced by LocalPhotonSelector().

Member Function Documentation

operator()()

bool LocalPhotonSelector::operator() ( CDPhoton & aPhoton )

override

Minimum energy

Barrel/Endcap

Time, only apply timing cuts if "recEvt->totalCharged() > 0"

Dang

Definition at line 32 of file LocalPhotonSelector.cxx.

                                                        {
 
  aPhoton.setUserTag( 1 );
  SmartDataPtr<EvtRecEvent>    recEvt( evtSvc(), EventModel::EvtRec::EvtRecEvent );
  SmartDataPtr<EvtRecTrackCol> recTrkCol( evtSvc(), EventModel::EvtRec::EvtRecTrackCol );
 
  EvtRecTrack* recTrk = const_cast<EvtRecTrack*>( aPhoton.photon() );
 
  if ( !recTrk->isEmcShowerValid() ) return false;
 
  RecEmcShower*    emcTrk = recTrk->emcShower();
  double           eraw   = emcTrk->energy();
  double           phi    = emcTrk->phi();
  double           the    = emcTrk->theta();
  HepLorentzVector shP4( eraw * sin( the ) * cos( phi ), eraw * sin( the ) * sin( phi ),
                         eraw * cos( the ), eraw );
 
  double cosThetaSh = shP4.vect().cosTheta();
 
  /// Minimum energy
  if ( shP4.e() <= m_minEnergy ) return false;
 
  /// Barrel/Endcap
  if ( m_useBarrelEndcap )
  {
    bool inBarrelEndcap = false;
 
    if ( fabs( cosThetaSh ) < m_maxCosThetaBarrel ) inBarrelEndcap = true;
 
    if ( ( fabs( cosThetaSh ) > m_minCosThetaEndcap ) &&
         ( fabs( cosThetaSh ) < m_maxCosThetaEndcap ) && ( shP4.e() > m_minEndcapEnergy ) )
      inBarrelEndcap = true;
 
    if ( !inBarrelEndcap ) return false;
  }
 
  /// Time, only apply timing cuts if "recEvt->totalCharged() > 0"
  if ( m_applyTimeCut )
  {
    double time = emcTrk->time();
    if ( recEvt->totalCharged() > 0 )
    {
      if ( time < m_minTime || time > m_maxTime ) return false;
    }
    else
    {
      RecEmcShower* firstG    = ( *( recTrkCol->begin() ) )->emcShower();
      double        deltaTime = fabs( time - firstG->time() );
      if ( deltaTime > 10 ) return false;
    }
  }
 
  /// Dang
  if ( m_applyDangCut && recEvt->totalCharged() > 0 )
  {
    Hep3Vector emcpos( emcTrk->x(), emcTrk->y(), emcTrk->z() );
    double     dang = 200.;
 
    for ( int j = 0; j < recEvt->totalCharged(); j++ )
    {
      EvtRecTrackIterator jtTrk = recTrkCol->begin() + j;
      if ( !( *jtTrk )->isExtTrackValid() ) continue;
      RecExtTrack* extTrk = ( *jtTrk )->extTrack();
      if ( extTrk->emcVolumeNumber() == -1 ) continue;
      Hep3Vector extpos = extTrk->emcPosition();
      double     angd1  = extpos.angle( emcpos );
      if ( angd1 < dang ) dang = angd1;
    }
 
    if ( dang < 200 )
    {
      dang = dang * 180 / ( CLHEP::pi );
      if ( dang <= m_minDang ) return ( false );
    }
  } // End of "recEvt->totalCharged() > 0" IF
 
  return true;
}

---

The documentation for this class was generated from the following files:

• LocalPhotonSelector.h
• LocalPhotonSelector.cxx