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;
}