/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/DataFormats/BTauReco/src/TauImpactParameterInfo.cc

Go to the documentation of this file.

#include "DataFormats/BTauReco/interface/TauImpactParameterInfo.h"
#include "DataFormats/TrackReco/interface/Track.h" 

using namespace edm;
using namespace reco;
using namespace std;

float reco::TauImpactParameterInfo::discriminator (double ip_min,double ip_max,double sip_min,bool use_sign, bool use3D) const {

        double discriminator = isolatedTaus->discriminator();

        const TrackRef leadingTrack = isolatedTaus->leadingSignalTrack(0.4,1.);

        if(! leadingTrack.isNull()){
          const TauImpactParameterTrackData* ipData = getTrackData(leadingTrack);
          Measurement1D ip = ipData->transverseIp;
          if(use3D) ip = ipData->ip3D;

          if( ip.value() < ip_min ||
              ip.value() > ip_max ||
              ip.significance() < sip_min ){
                discriminator = 0;
          }
        }
        return discriminator;
}
float reco::TauImpactParameterInfo::discriminator() const {
        //default discriminator: returns the value of the discriminator of the jet tag
        return isolatedTaus->discriminator();
}

const reco::TauImpactParameterTrackData* TauImpactParameterInfo::getTrackData(const reco::TrackRef & trackRef) const {

        reco::TrackTauImpactParameterAssociationCollection::const_iterator iter
            = trackDataMap.find(trackRef);

        if (iter != trackDataMap.end()) return &(iter->val);

        return 0; // if track not found return 0
}

void reco::TauImpactParameterInfo::storeTrackData(const reco::TrackRef & trackRef,
                  const reco::TauImpactParameterTrackData& trackData) {

        trackDataMap.insert(trackRef, trackData);
}

void reco::TauImpactParameterInfo::setIsolatedTauTag(const IsolatedTauTagInfoRef & isolationRef){
        isolatedTaus = isolationRef;
}

const IsolatedTauTagInfoRef& reco::TauImpactParameterInfo::getIsolatedTauTag() const {
        return isolatedTaus;
}