MITConversionUtilities.cc

Go to the documentation of this file.

#include "DQM/PhysicsHWW/interface/MITConversionUtilities.h"
#include "DataFormats/TrackReco/interface/TrackBase.h"
#include "TMath.h"

namespace HWWFunctions {

  bool isMITConversion(HWW& hww, unsigned int elidx, 
           int nWrongHitsMax, 
           float probMin,
           float dlMin,
           bool matchCTF,
           bool requireArbitratedMerged) {

    unsigned int nconvs = hww.convs_isConverted().size();
    if(nconvs == 0) 
      return false;
    bool isGoodConversion = false;

    for(unsigned int iconv = 0; iconv < nconvs; iconv++) {
      
      bool conversionMatchFound = false;
      for(unsigned int itk = 0; itk < hww.convs_tkidx().at(iconv).size(); itk++) {

        if(hww.convs_tkalgo().at(iconv)[itk] == reco::TrackBase::gsf && hww.convs_tkidx().at(iconv)[itk] == hww.els_gsftrkidx().at(elidx))
    conversionMatchFound = true;
        if(matchCTF) {
          switch(hww.convs_tkalgo().at(iconv)[itk]) {
          case reco::TrackBase::initialStep:
          case reco::TrackBase::lowPtTripletStep:
          case reco::TrackBase::pixelPairStep:
          case reco::TrackBase::detachedTripletStep:
          case reco::TrackBase::mixedTripletStep:
          case reco::TrackBase::pixelLessStep:
          case reco::TrackBase::tobTecStep:
          case reco::TrackBase::jetCoreRegionalStep:
            if(hww.convs_tkidx().at(iconv)[itk] == hww.els_trkidx().at(elidx))
              conversionMatchFound = true;
            break;
          default:
            break;
          }
        }
      
        if(conversionMatchFound)
    break;
      }
      
      
      if(conversionMatchFound==false)
        continue;
      
      if( TMath::Prob( hww.convs_chi2().at(iconv), (Int_t)hww.convs_ndof().at(iconv) )  > probMin && hww.convs_dl().at(iconv) > dlMin ) isGoodConversion = true;
      if(requireArbitratedMerged) {
        if(hww.convs_quality().at(iconv) & 4)
    isGoodConversion = true;
        else 
    isGoodConversion = false;
      }

      for(unsigned int j = 0; j < hww.convs_nHitsBeforeVtx().at(iconv).size(); j++) {
        if(hww.convs_nHitsBeforeVtx().at(iconv)[j] > nWrongHitsMax)
    isGoodConversion = false;
      }
        
      if(isGoodConversion)
        break;
        
        
    }//loop over convserions


    return isGoodConversion;
  }
}