CaloJet.cc

Go to the documentation of this file.

// CaloJet.cc
// Fedor Ratnikov UMd
// $Id: CaloJet.cc,v 1.22 2009/04/16 20:04:20 srappocc Exp $
#include <sstream>

#include "FWCore/Utilities/interface/Exception.h"
#include "DataFormats/RecoCandidate/interface/RecoCaloTowerCandidate.h"

//Own header file
#include "DataFormats/JetReco/interface/CaloJet.h"

using namespace reco;

 CaloJet::CaloJet (const LorentzVector& fP4, const Point& fVertex, 
           const Specific& fSpecific)
   : Jet (fP4, fVertex),
     m_specific (fSpecific)
 {}

 CaloJet::CaloJet (const LorentzVector& fP4, const Point& fVertex, 
          const Specific& fSpecific, 
          const Jet::Constituents& fConstituents)
   : Jet (fP4, fVertex, fConstituents),
     m_specific (fSpecific)
 {}

CaloJet::CaloJet (const LorentzVector& fP4, 
          const Specific& fSpecific, 
          const Jet::Constituents& fConstituents)
  : Jet (fP4, Point(0,0,0), fConstituents),
    m_specific (fSpecific)
{}


//float CaloJet::physicsEtaQuick (float fZVertex) const {
//  return Jet::physicsEta (fZVertex, eta());
//}

//float CaloJet::physicsEtaDetailed (float fZVertex) const {
//  Jet::LorentzVector correctedMomentum;
//  std::vector<const Candidate*> towers = getJetConstituentsQuick ();
//  for (unsigned i = 0; i < towers.size(); ++i) {
//    const Candidate* c = towers[i];
//    double etaRef = Jet::physicsEta (fZVertex, c->eta());
//    math::PtEtaPhiMLorentzVectorD p4 (c->p()/cosh(etaRef), etaRef, c->phi(), c->mass());
//    correctedMomentum += p4;
//  }
//  return correctedMomentum.eta();
//}

//deprecated, with 3d vertex correction not clear anymore!
//CaloJet::LorentzVector CaloJet::physicsP4 (float fZVertex) const {
//  double physicsEta = Jet::physicsEta (fZVertex, eta());
//  math::PtEtaPhiMLorentzVectorD p4 (p()/cosh(physicsEta), physicsEta, phi(), mass());
//  return CaloJet::LorentzVector (p4);
//}

CaloJet::LorentzVector CaloJet::physicsP4 (const Particle::Point &vertex) const {
  return Jet::physicsP4(vertex,*this,this->vertex());
}

CaloJet::LorentzVector CaloJet::detectorP4 () const {
  return Jet::detectorP4(this->vertex(),*this);
}


CaloTowerPtr CaloJet::getCaloConstituent (unsigned fIndex) const {
   Constituent dau = daughterPtr (fIndex);

   if ( dau.isNonnull() && dau.isAvailable() ) {

   const CaloTower* towerCandidate = dynamic_cast <const CaloTower*> (dau.get());

    if (towerCandidate) {
//      return towerCandidate;
// 086     Ptr(ProductID const& productID, T const* item, key_type item_key) :
      return edm::Ptr<CaloTower> (dau.id(), towerCandidate, dau.key() );
    }
    else {
      throw cms::Exception("Invalid Constituent") << "CaloJet constituent is not of CaloTowere type";
    }

   }

   else {
     return CaloTowerPtr();
   }
}


std::vector <CaloTowerPtr > CaloJet::getCaloConstituents () const {
  std::vector <CaloTowerPtr> result;
  for (unsigned i = 0;  i <  numberOfDaughters (); i++) result.push_back (getCaloConstituent (i));
  return result;
}


CaloJet* CaloJet::clone () const {
  return new CaloJet (*this);
}

bool CaloJet::overlap( const Candidate & ) const {
  return false;
}

std::string CaloJet::print () const {
  std::ostringstream out;
  out << Jet::print () // generic jet info
      << "    CaloJet specific:" << std::endl
      << "      energy fractions em/had: " << emEnergyFraction () << '/' << energyFractionHadronic () << std::endl
      << "      em energy in EB/EE/HF: " << emEnergyInEB() << '/' << emEnergyInEE() << '/' << emEnergyInHF() << std::endl
      << "      had energy in HB/HO/HE/HF: " << hadEnergyInHB() << '/' << hadEnergyInHO() << '/' << hadEnergyInHE() << '/' << hadEnergyInHF() << std::endl
      << "      constituent towers area: " << towersArea() << std::endl;
  out << "      Towers:" << std::endl;
  std::vector <CaloTowerPtr > towers = getCaloConstituents ();
  for (unsigned i = 0; i < towers.size (); i++) {
    if (towers[i].get ()) {
      out << "      #" << i << " " << *(towers[i]) << std::endl;
    }
    else {
      out << "      #" << i << " tower is not available in the event"  << std::endl;
    }
  }
  return out.str ();
}

//----------------------------------------------------------
// here are methods extracting information from constituents
//----------------------------------------------------------

std::vector<CaloTowerDetId> CaloJet::getTowerIndices() const {
  std::vector<CaloTowerDetId> result;
  std::vector <CaloTowerPtr> towers = getCaloConstituents ();
  for (unsigned i = 0; i < towers.size(); ++i) {
    result.push_back (towers[i]->id());
  }
  return result;
}