StopMuFinderFar.cxx

Go to the documentation of this file.

#include <list>
#include <cmath>

#include "UgliGeometry/UgliGeomHandle.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "CandData/CandRecord.h"
#include "CandData/CandHeader.h"
#include "RecoBase/CandFitTrackListHandle.h"
#include "RecoBase/CandFitTrackHandle.h"
#include "DataUtil/GetCandidate.h"
#include "DataUtil/GetCandHeader.h"
#include "DataUtil/CDL2STL.h"

#include "StopMuFinderFactory.h"
#include "StopMuFinderFar.h"

CVSID("$Id: StopMuFinderFar.cxx,v 1.3 2007/02/04 05:49:29 rhatcher Exp $");

using std::list;
using std::string;
using DataUtil::GetCandidate;
using DataUtil::GetCandHeader;
using DataUtil::CDL2STLlist;

static const string STOP_MU_FINDER_FAR = "StopMuFinderFar";

namespace {

    StopMuFinder* CreateSMF() { return new StopMuFinderFar; }

    bool registered = StopMuFinderFactory::Instance().
            RegisterStopMuFinder(STOP_MU_FINDER_FAR, CreateSMF);

} // namespace 


StopMuFinderFar::StopMuFinderFar()
{
    MSG("StopMuFilter",Msg::kDebug) 
        << "StopMuFinderFar::StopMuFinderFar()" 
        << endl;

    DefaultConfig();
} // StopMuFinderFar::StopMuFinderFar() 



void StopMuFinderFar::DefaultConfig(void)
{
    MSG("StopMuFilter",Msg::kDebug)
        << "StopMuFinderFar::DefaultConfig()"
        << endl;

    fMinZDistance   =       0.5;
    fMinRadius      =       0.5;
    fMaxRadius      =       3.5;
    fMinNPlanes     =       10;
}



void StopMuFinderFar::Configure(const Registry & registry)
{
    MSG("StopMuFilter",Msg::kDebug)
        << "StopMuFinderFar::Configure(const Registry &)"
        << endl;

    if ( registry.KeyExists("CandTrackList") ) {
        fCandTrackListName = std::string(registry.GetCharString("CandTrackList"));
    }
    if ( registry.KeyExists("MinZDistance") &&
            registry.GetType("MinZDistance") == typeid(Double_t) ) {
        fMinZDistance = registry.GetDouble("MinZDistance");
    }
    if ( registry.KeyExists("MinRadius") &&
            registry.GetType("MinRadius") == typeid(Double_t) ) {
        fMinRadius = registry.GetDouble("MinRadius");
    }
    if ( registry.KeyExists("MaxRadius") &&
            registry.GetType("MaxRadius") == typeid(Double_t) ) {
        fMaxRadius = registry.GetDouble("MaxRadius");
    }
    if ( registry.KeyExists("MinNPlanes") &&
            registry.GetType("MinNPlanes") == typeid(Int_t) ) {
        fMinNPlanes = registry.GetInt("MinNPlanes");
    }
}



Bool_t StopMuFinderFar::IsStopping(const MomNavigator* mom)
{
    MSG("StopMuFilter",Msg::kDebug)
        << "StopMuFinderFar::IsStopping(const MomNavigator*)"
        << endl;

    // get CandFitTrackList
    const CandFitTrackListHandle* candFitTrackList =
        GetCandidate<CandFitTrackListHandle> (mom, fCandTrackListName.c_str());

    // no candidate list found
    if ( !candFitTrackList ) return kFALSE;

    // set Z extent
    // doing it every event makes filter slow, so
    // instead do it once and save in static variables
    static Bool_t   zExtentSet  = kFALSE;
    static Double_t zMinSM1     = 0.;
    static Double_t zMaxSM1     = 0.;
    static Double_t zMinSM2     = 0.;
    static Double_t zMaxSM2     = 0.;

    static const Int_t Sm1 = 0;    // due to some bizarro convention id of SM1 is "0"
    static const Int_t Sm2 = 1;    // and id of SM2 is "1"

    if ( ! zExtentSet ) {
        // Get CandRecord from mom
        const CandHeader* candHeader = GetCandHeader ( mom );
        if ( !candHeader ) {
            MSG("StopMuFinderFar",Msg::kDebug)
                << "No CandHeader!!" << endl;
            return kFALSE;
        }

        const VldContext vldc = candHeader->GetVldContext();

        UgliGeomHandle ugh(vldc);

        ugh.GetZExtent(zMinSM1, zMaxSM1, Sm1);
        ugh.GetZExtent(zMinSM2, zMaxSM2, Sm2);

        zExtentSet = kTRUE;
    }

    typedef list<const CandFitTrackHandle*> FitTrackList;
    typedef list<const CandFitTrackHandle*>::const_iterator FitTrackListItor;

    // Convert daughter list to STL
    FitTrackList candFitTrackListSTL =
        CDL2STLlist<CandFitTrackHandle>(*candFitTrackList);

    // for now only look at events with only one track
    if ( candFitTrackListSTL.empty() ) return kFALSE;

    const CandFitTrackHandle* pcfth = *candFitTrackListSTL.begin();

    // check number of planes
    if ( pcfth->GetNPlane() < fMinNPlanes ) return kFALSE;

    // get track end coord's
    Double_t endu = pcfth->GetEndU();
    Double_t endv = pcfth->GetEndV();
    Double_t endz = pcfth->GetEndZ();

    // check endtrack position along Z axis
    if ( endz<(zMinSM1+fMinZDistance) ||
        ( endz>(zMaxSM1-fMinZDistance) && endz<(zMinSM2+fMinZDistance) ) ||
        endz>(zMaxSM2-fMinZDistance) ) return kFALSE;

    // check against radius
    Double_t radius = sqrt(endu*endu + endv*endv);
    if ( radius<fMinRadius || radius>fMaxRadius ) return kFALSE;

    return kTRUE;
} // void StopMuFinderFar::IsStopping(const MomNavigator* mom)



StopMuFinderFar::~StopMuFinderFar()
{
    MSG("StopMuFilter",Msg::kDebug)
        << "StopMuFinderFar::~StopMuFinderFar()"
        << endl;
} // StopMuFinderFar::~StopMuFinderFar()

---