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 ()
bool	operator() (CDPi0 &iArg) const

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 12 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 43 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 32 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()

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

inlinestatic

Definition at line 90 of file Pi0Cut.h.

{ recEvt =_recEvt; recTrkCol=_recTrkCol;}

Referenced by Pi0RecAlg::execute().

Member Data Documentation

Chisq

double UserPi0Cut::Chisq

Definition at line 30 of file Pi0Cut.h.

Referenced by operator()().

MaxMass

double UserPi0Cut::MaxMass

Definition at line 29 of file Pi0Cut.h.

Referenced by Pi0RecAlg::execute(), and operator()().

MinAngle

double UserPi0Cut::MinAngle

Definition at line 26 of file Pi0Cut.h.

Referenced by isGoodPhoton().

MinEnergyOfBarrelPhoton

double UserPi0Cut::MinEnergyOfBarrelPhoton

Definition at line 24 of file Pi0Cut.h.

Referenced by isGoodPhoton().

MinEnergyOfEndcapPhoton

double UserPi0Cut::MinEnergyOfEndcapPhoton

Definition at line 25 of file Pi0Cut.h.

Referenced by isGoodPhoton().

MinMass

double UserPi0Cut::MinMass

Definition at line 28 of file Pi0Cut.h.

Referenced by Pi0RecAlg::execute(), and operator()().

recEvt

EvtRecEvent * UserPi0Cut::recEvt = 0

static

Definition at line 91 of file Pi0Cut.h.

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

recTrkCol

EvtRecTrackCol * UserPi0Cut::recTrkCol = 0

static

Definition at line 92 of file Pi0Cut.h.

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

TimeLimit

double UserPi0Cut::TimeLimit

Definition at line 27 of file Pi0Cut.h.

Referenced by isGoodPhoton().

---

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

• Pi0Cut.h
• Pi0Cut.cxx