TrackClusterAssociationMVAAlgorithm.h

Go to the documentation of this file.

/*
 *
 * TrackClusterAssociationMVAAlgorithm.h header template automatically generated by a class generator
 * Creation date : jeu. avr. 9 2015
 *
 * This file is part of APRILContent libraries.
 * 
 * APRILContent is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 * based upon these libraries are permitted. Any copy of these libraries
 * must include this copyright notice.
 * 
 * APRILContent is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with APRILContent.  If not, see <http://www.gnu.org/licenses/>.
 * 
 * @author Remi Ete
 * @copyright CNRS , IPNL
 */


#ifndef TRACKCLUSTERASSOCIATIONMVAALGORITHM_H
#define TRACKCLUSTERASSOCIATIONMVAALGORITHM_H

#include "Pandora/Algorithm.h"

#include "TMVA/Reader.h"
#include "TMVA/MethodCuts.h"

namespace april_content
{

class TrackClusterAssociationMVAAlgorithm : public pandora::Algorithm
{
 public:
    class Factory : public pandora::AlgorithmFactory
    {
    public:
        pandora::Algorithm *CreateAlgorithm() const;
    };

private:
    pandora::StatusCode Run();
    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);

    struct AssociationConstraints
    {
        const pandora::Cluster *m_pCluster;
        float                   m_chi2;
        float                   m_trackClusterAngle;
        float                   m_distanceToHelix;
        float                   m_compatibility;
    };

    typedef std::vector<AssociationConstraints> AssociationConstraintsVector;
    typedef std::map<const pandora::Track*, AssociationConstraintsVector> AssociationConstraintsMap;
    typedef std::map<const pandora::Track*, const pandora::Cluster*> AssociationMap;
    typedef std::map<const pandora::Cluster*, const pandora::Track*> MultiAssociationMap;

    pandora::StatusCode ExtractAssociationLists(pandora::ClusterList &clusterList, pandora::TrackList &trackList);

    pandora::StatusCode BuildPossibleAssociations(const pandora::ClusterList &clusterList, const pandora::TrackList &trackList,
            AssociationMap &associationMap);

    pandora::StatusCode PerformAssociations(const AssociationMap &associationMap);

    pandora::StatusCode ComputeInnerCentroid(const pandora::Cluster *const pCluster, const pandora::Track *const pTrack, pandora::CartesianVector &innerCentroid) const;

    bool PassesInitialCuts(const pandora::Cluster *const pCluster, const pandora::Track *const pTrack) const;

    float MVADeteriminAssociation(const pandora::Cluster* const pCluster, const pandora::Track* const pTrack) const;

    pandora::StatusCode CalculateAssociationConstraints(const pandora::Cluster *const pCluster, const pandora::Track *const pTrack, AssociationConstraints &constraints);

    bool FitConstraints(const AssociationConstraints &constraints) const;

    static bool SortByBestCompatibility(const AssociationConstraints &lhs, const AssociationConstraints &rhs);

private:
    unsigned int                   m_maxClusterInnerPseudoLayer;   
    unsigned int                   m_nFirstClusterPseudoLayer;     

    float                          m_maxDistanceToHelix;           
    float                          m_maxTrackClusterAngle;         

    bool                           m_allowMultiAssociations;       
    bool                           m_useEnergyCompatibility;       

    bool                           m_usePhotonClusters;

    static TMVA::Reader*           m_reader;
};

//------------------------------------------------------------------------------------------------------------------------------------------

inline pandora::Algorithm *TrackClusterAssociationMVAAlgorithm::Factory::CreateAlgorithm() const
{
    return new TrackClusterAssociationMVAAlgorithm();
}

} 

#endif  //  TRACKCLUSTERASSOCIATIONMVAALGORITHM_H