ConversionVertexFinder Class Reference

#include <ConversionVertexFinder.h>

Public Member Functions
	ConversionVertexFinder (const edm::ParameterSet &config)
TransientVertex	run (const std::vector< reco::TransientTrack > &pair)
bool	run (const std::vector< reco::TransientTrack > &pair, reco::Vertex &the_vertex)
	~ConversionVertexFinder ()

Private Attributes
edm::ParameterSet	conf_
KinematicConstrainedVertexFitter *	kcvFitter_
double	maxDelta_
int	maxNbrOfIterations_
double	maxReducedChiSq_
double	minChiSqImprovement_

Detailed Description

Author

N. Marinelli - Univ. of Notre Dame

Version

Definition at line 29 of file ConversionVertexFinder.h.

Constructor & Destructor Documentation

ConversionVertexFinder()

ConversionVertexFinder::ConversionVertexFinder

(

const edm::ParameterSet &

config

)

Definition at line 19 of file ConversionVertexFinder.cc.

References conf_, edm::ParameterSet::getParameter(), kcvFitter_, LogDebug, maxDelta_, maxNbrOfIterations_, maxReducedChiSq_, minChiSqImprovement_, and KinematicConstrainedVertexFitter::setParameters().

                                                                            : conf_(config) {
  LogDebug("ConversionVertexFinder") << "ConversionVertexFinder CTOR  "
                                     << "\n";
  maxDelta_ = conf_.getParameter<double>("maxDelta");
  maxReducedChiSq_ = conf_.getParameter<double>("maxReducedChiSq");
  minChiSqImprovement_ = conf_.getParameter<double>("minChiSqImprovement");
  maxNbrOfIterations_ = conf_.getParameter<int>("maxNbrOfIterations");
  kcvFitter_ = new KinematicConstrainedVertexFitter();
  kcvFitter_->setParameters(conf_);
}

~ConversionVertexFinder()

ConversionVertexFinder::~ConversionVertexFinder

(

)

Definition at line 30 of file ConversionVertexFinder.cc.

References kcvFitter_, and LogDebug.

                                                {
  LogDebug("ConversionVertexFinder") << "ConversionVertexFinder DTOR "
                                     << "\n";
  delete kcvFitter_;
}

Member Function Documentation

run() [1/2]

TransientVertex ConversionVertexFinder::run

(

const std::vector< reco::TransientTrack > &

pair

)

Definition at line 135 of file ConversionVertexFinder.cc.

References MillePedeFileConverter_cfg::e, LogDebug, metname, AlCaHLTBitMon_QueryRunRegistry::string, and KalmanVertexFitter::vertex().

Referenced by ConvertedPhotonProducer::buildCollections(), and ConversionProducer::checkVertex().

                                                                                     {
  LogDebug("ConversionVertexFinder") << "ConversionVertexFinder run pair size " << pair.size() << "\n";

  //for ( std::vector<reco::TransientTrack>::const_iterator iTk=pair.begin(); iTk!=pair.end(); ++iTk) {
  // LogDebug("ConversionVertexFinder") << "  ConversionVertexFinder  Tracks in the pair  charge " << iTk->charge() << " Num of RecHits " << iTk->recHitsSize() << " inner momentum " << iTk->track().innerMomentum() << "\n";
  //}

  reco::Vertex theVertex;
  KalmanVertexFitter fitter(true);
  TransientVertex transientVtx;

  const std::string metname = "ConversionVertexFinder| ConversionVertexFinder";
  try {
    transientVtx = fitter.vertex(pair);

  } catch (cms::Exception& e) {
    edm::LogWarning(metname) << "cms::Exception caught in ConversionVertexFinder::run\n" << e.explainSelf();
  }

  return transientVtx;
}

run() [2/2]

bool ConversionVertexFinder::run	(	const std::vector< reco::TransientTrack > &	pair,
		reco::Vertex &	the_vertex
	)

Definition at line 36 of file ConversionVertexFinder.cc.

References ChiSquaredProbability(), conf_, gather_cfg::cout, KinematicConstrainedVertexFitter::fit(), KinematicConstrainedVertexFitterT< nTrk, nConstraint >::fit(), newFWLiteAna::found, KinematicConstrainedVertexFitter::getNit(), KinematicConstrainedVertexFitterT< nTrk, nConstraint >::getNit(), kcvFitter_, GetRecoTauVFromDQM_MC_cff::kk, EgHLTOffHistBins_cfi::mass, KinematicParticleFactoryFromTransientTrack::particle(), ecalTrigSettings_cff::particles, ColinearityKinematicConstraint::PhiTheta, and KinematicConstrainedVertexFitterT< nTrk, nConstraint >::setParameters().

                                                                                                   {
  std::vector<reco::TransientTrack> pair = _pair;
  bool found = false;

  if (pair.size() < 2)
    return found;

  float sigma = 0.00000001;
  float chi = 0.;
  float ndf = 0.;
  float mass = 0.000511;

  KinematicParticleFactoryFromTransientTrack pFactory;

  std::vector<RefCountedKinematicParticle> particles;

  particles.push_back(
      pFactory.particle(pair[0], mass, chi, ndf, sigma, *pair[0].innermostMeasurementState().freeState()));
  particles.push_back(
      pFactory.particle(pair[1], mass, chi, ndf, sigma, *pair[1].innermostMeasurementState().freeState()));

#ifdef OldKineFit
  ColinearityKinematicConstraint constr(ColinearityKinematicConstraint::PhiTheta);

  RefCountedKinematicTree myTree = kcvFitter_->fit(particles, &constr);
#else

  // bizzare way to the get the field...
  const MagneticField* mf = pair[0].field();

  ColinearityKinematicConstraintT<colinearityKinematic::PhiTheta> constr;
  KinematicConstrainedVertexFitterT<2, 2> kcvFitter(mf);
  kcvFitter.setParameters(conf_);
  RefCountedKinematicTree myTree = kcvFitter.fit(particles, &constr);

#ifdef KineFitDebug

  ColinearityKinematicConstraint oldconstr(ColinearityKinematicConstraint::PhiTheta);

  RefCountedKinematicTree oldTree = kcvFitter_->fit(particles, &oldconstr);

  if (oldTree->isValid()) {
    std::cout << "old " << kcvFitter_->getNit() << std::endl;
    RefCountedKinematicVertex gamma_dec_vertex = oldTree->currentDecayVertex();
    if (gamma_dec_vertex->vertexIsValid())
      std::cout << gamma_dec_vertex->chiSquared() << " " << gamma_dec_vertex->degreesOfFreedom() << std::endl;
    std::cout << oldTree->currentParticle()->currentState().globalMomentum()
              << oldTree->currentParticle()->currentState().globalPosition() << std::endl;
    std::vector<RefCountedKinematicParticle> fStates = oldTree->finalStateParticles();
    for (unsigned int kk = 0; kk < fStates.size(); kk++) {
      std::cout << fStates[kk]->currentState().globalMomentum() << fStates[kk]->currentState().globalPosition()
                << std::endl;
      std::cout << fStates[kk]->currentState().kinematicParametersError().matrix() << std::endl;
    }
  } else
    std::cout << "old invalid " << kcvFitter_->getNit() << std::endl;

  if (myTree->isValid()) {
    std::cout << "new " << kcvFitter.getNit() << std::endl;
    RefCountedKinematicVertex gamma_dec_vertex = myTree->currentDecayVertex();
    if (gamma_dec_vertex->vertexIsValid())
      std::cout << gamma_dec_vertex->chiSquared() << " " << gamma_dec_vertex->degreesOfFreedom() << std::endl;
    std::cout << myTree->currentParticle()->currentState().globalMomentum()
              << myTree->currentParticle()->currentState().globalPosition() << std::endl;
    std::vector<RefCountedKinematicParticle> fStates = myTree->finalStateParticles();
    for (unsigned int kk = 0; kk < fStates.size(); kk++) {
      std::cout << fStates[kk]->currentState().globalMomentum() << fStates[kk]->currentState().globalPosition()
                << std::endl;
      std::cout << fStates[kk]->currentState().kinematicParametersError().matrix() << std::endl;
    }
  } else
    std::cout << "new invalid " << kcvFitter.getNit() << std::endl;

#endif  // TemplateKineFitDebug

#endif

  if (myTree->isValid()) {
    myTree->movePointerToTheTop();
    RefCountedKinematicParticle the_photon = myTree->currentParticle();
    if (the_photon->currentState().isValid()) {
      //const ParticleMass photon_mass = the_photon->currentState().mass();
      RefCountedKinematicVertex gamma_dec_vertex;
      gamma_dec_vertex = myTree->currentDecayVertex();
      if (gamma_dec_vertex->vertexIsValid()) {
        const float chi2Prob =
            ChiSquaredProbability(gamma_dec_vertex->chiSquared(), gamma_dec_vertex->degreesOfFreedom());
        if (chi2Prob > 0.) {  // no longer cut here, only ask positive probability here
          //const math::XYZPoint vtxPos(gamma_dec_vertex->position());
          the_vertex = *gamma_dec_vertex;
          found = true;
        }
      }
    }
  }

  return found;
}

Member Data Documentation

conf_

edm::ParameterSet ConversionVertexFinder::conf_

private

Definition at line 40 of file ConversionVertexFinder.h.

Referenced by ConversionVertexFinder(), and run().

kcvFitter_

KinematicConstrainedVertexFitter* ConversionVertexFinder::kcvFitter_

private

Definition at line 45 of file ConversionVertexFinder.h.

Referenced by ConversionVertexFinder(), run(), and ~ConversionVertexFinder().

maxDelta_

double ConversionVertexFinder::maxDelta_

private

Definition at line 41 of file ConversionVertexFinder.h.

Referenced by ConversionVertexFinder().

maxNbrOfIterations_

int ConversionVertexFinder::maxNbrOfIterations_

private

Definition at line 44 of file ConversionVertexFinder.h.

Referenced by ConversionVertexFinder().

maxReducedChiSq_

double ConversionVertexFinder::maxReducedChiSq_

private

Definition at line 42 of file ConversionVertexFinder.h.

Referenced by ConversionVertexFinder().

minChiSqImprovement_

double ConversionVertexFinder::minChiSqImprovement_

private

Definition at line 43 of file ConversionVertexFinder.h.

Referenced by ConversionVertexFinder().