UserPi0Cut Struct Reference

#include <Pi0Cut.h>

Inheritance diagram for UserPi0Cut:

Public Member Functions
	UserPi0Cut (double MEBP=0.025, double MEEP=0.045, double MA=20., double TL=25., double MinM=0.12, double MaxM=0.15, double c=1e+10)
bool	operator() (CDPi0 &_this)
bool	isGoodPhoton (EvtRecTrack *gTrk)
Public Member Functions inherited from DCSelectionFunction< CDPi0 >
	DCSelectionFunction ()
virtual	~DCSelectionFunction ()

Static Public Member Functions
static void	SetForTrack (EvtRecEvent *_recEvt, EvtRecTrackCol *_recTrkCol)

Public Attributes
double	MinEnergyOfBarrelPhoton
double	MinEnergyOfEndcapPhoton
double	MinAngle
double	TimeLimit
double	MinMass
double	MaxMass
double	Chisq

Static Public Attributes
static EvtRecEvent *	recEvt = 0
static EvtRecTrackCol *	recTrkCol = 0

Detailed Description

Definition at line 10 of file Pi0Cut.h.

Constructor & Destructor Documentation

UserPi0Cut()

UserPi0Cut::UserPi0Cut ( double MEBP = 0.025, double MEEP = 0.045, double MA = 20., double TL = 25., double MinM = 0.12, double MaxM = 0.15, double c = 1e+10 )

inline

Definition at line 11 of file Pi0Cut.h.

      : MinEnergyOfBarrelPhoton( MEBP )
      , MinEnergyOfEndcapPhoton( MEEP )
      , MinAngle( MA )
      , TimeLimit( TL )
      , Chisq( c )
      , MinMass( MinM )
      , MaxMass( MaxM ) {}

Member Function Documentation

isGoodPhoton()

bool UserPi0Cut::isGoodPhoton ( EvtRecTrack * gTrk )

inline

Definition at line 40 of file Pi0Cut.h.

                                         {
    if ( !gTrk->isEmcShowerValid() ) return false;
 
    RecEmcShower* emcTrk = gTrk->emcShower();
    // check time
    if ( emcTrk->time() > TimeLimit ) return false;
    // check energy
    int flag = ( emcTrk->cellId() & 0x000F0000 ) >> 16;
    if ( ( flag == 1 ) && emcTrk->energy() < MinEnergyOfBarrelPhoton ) return false;
    if ( ( flag == 0 || flag == 2 ) && emcTrk->energy() < MinEnergyOfEndcapPhoton )
      return false;
    // check angle
    Hep3Vector emcpos( emcTrk->x(), emcTrk->y(), emcTrk->z() );
    if ( recEvt->totalCharged() > 0 )
    {
      double dthe  = 200.;
      double dphi  = 200.;
      double dang1 = 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 );
        double     thed   = extpos.theta() - emcpos.theta();
        double     phid   = extpos.deltaPhi( emcpos );
        thed =
            fmod( thed + CLHEP::twopi + CLHEP::twopi + CLHEP::pi, CLHEP::twopi ) - CLHEP::pi;
        phid =
            fmod( phid + CLHEP::twopi + CLHEP::twopi + CLHEP::pi, CLHEP::twopi ) - CLHEP::pi;
 
        if ( fabs( thed ) < fabs( dthe ) ) dthe = thed;
        if ( fabs( phid ) < fabs( dphi ) ) dphi = phid;
        if ( angd1 < dang1 ) dang1 = angd1;
      }
      if ( dang1 < 200 )
      {
        dthe  = dthe * 180 / ( CLHEP::pi );
        dphi  = dphi * 180 / ( CLHEP::pi );
        dang1 = dang1 * 180 / ( CLHEP::pi );
        if ( dang1 < MinAngle ) return false;
 
      } // End of "dang1 < 200" IF
 
    } // End of "recEv
    return true;
  }

Referenced by operator()().

operator()()

bool UserPi0Cut::operator() ( CDPi0 & _this )

inlinevirtual

Implements DCSelectionFunction< CDPi0 >.

Definition at line 29 of file Pi0Cut.h.

                                  {
    EvtRecTrack* g1 = const_cast<EvtRecTrack*>( _this.navPi0()->hiEnGamma() );
    EvtRecTrack* g2 = const_cast<EvtRecTrack*>( _this.navPi0()->loEnGamma() );
    if ( !isGoodPhoton( g1 ) ) return false;
    if ( !isGoodPhoton( g2 ) ) return false;
    if ( _this.navPi0()->chisq() > Chisq ) return false;
    if ( _this.navPi0()->unconMass() < MinMass || _this.navPi0()->unconMass() > MaxMass )
      return false;
    return true;
  }

SetForTrack()

void UserPi0Cut::SetForTrack ( EvtRecEvent * _recEvt, EvtRecTrackCol * _recTrkCol )

inlinestatic

Definition at line 91 of file Pi0Cut.h.

                                                                              {
    recEvt    = _recEvt;
    recTrkCol = _recTrkCol;
  }

Referenced by Pi0RecAlg::execute().

Member Data Documentation

Chisq

double UserPi0Cut::Chisq

Definition at line 27 of file Pi0Cut.h.

Referenced by operator()(), and UserPi0Cut().

MaxMass

double UserPi0Cut::MaxMass

Definition at line 26 of file Pi0Cut.h.

Referenced by operator()(), and UserPi0Cut().

MinAngle

double UserPi0Cut::MinAngle

Definition at line 23 of file Pi0Cut.h.

Referenced by isGoodPhoton(), and UserPi0Cut().

MinEnergyOfBarrelPhoton

double UserPi0Cut::MinEnergyOfBarrelPhoton

Definition at line 21 of file Pi0Cut.h.

Referenced by isGoodPhoton(), and UserPi0Cut().

MinEnergyOfEndcapPhoton

double UserPi0Cut::MinEnergyOfEndcapPhoton

Definition at line 22 of file Pi0Cut.h.

Referenced by isGoodPhoton(), and UserPi0Cut().

MinMass

double UserPi0Cut::MinMass

Definition at line 25 of file Pi0Cut.h.

Referenced by operator()(), and UserPi0Cut().

recEvt

EvtRecEvent * UserPi0Cut::recEvt = 0

static

Definition at line 95 of file Pi0Cut.h.

Referenced by isGoodPhoton(), Pi0::make_gamma_list(), and SetForTrack().

recTrkCol

EvtRecTrackCol * UserPi0Cut::recTrkCol = 0

static

Definition at line 96 of file Pi0Cut.h.

Referenced by isGoodPhoton(), Pi0::make_gamma_list(), and SetForTrack().

TimeLimit

double UserPi0Cut::TimeLimit

Definition at line 24 of file Pi0Cut.h.

Referenced by isGoodPhoton(), and UserPi0Cut().

---

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

• Pi0Cut.h
• Pi0Cut.cxx