---

VisJet.cc

Go to the documentation of this file.

//<<<<<< INCLUDES                                                       >>>>>>

#include "VisReco/Analyzer/interface/VisJet.h"
#include "VisReco/Analyzer/interface/IguanaService.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "Iguana/Framework/interface/IgCollection.h"
#include <iostream>
#include <sstream>

//<<<<<< PRIVATE DEFINES                                                >>>>>>
//<<<<<< PRIVATE CONSTANTS                                              >>>>>>
//<<<<<< PRIVATE TYPES                                                  >>>>>>
//<<<<<< PRIVATE VARIABLE DEFINITIONS                                   >>>>>>
//<<<<<< PUBLIC VARIABLE DEFINITIONS                                    >>>>>>
//<<<<<< CLASS STRUCTURE INITIALIZATION                                 >>>>>>
//<<<<<< PRIVATE FUNCTION DEFINITIONS                                   >>>>>>
//<<<<<< PUBLIC FUNCTION DEFINITIONS                                    >>>>>>
//<<<<<< MEMBER FUNCTION DEFINITIONS                                    >>>>>>

using namespace edm::service;

VisJet::VisJet (const edm::ParameterSet& iConfig)
    : inputTag_ (iConfig.getParameter<edm::InputTag>("visJetTag")),
      energyCut_(iConfig.getUntrackedParameter<double>("energyCut", 0.1))
{}

void
VisJet::analyze( const edm::Event& event, const edm::EventSetup& eventSetup)
{
    edm::Service<IguanaService> config;
    if (! config.isAvailable ()) 
    {
        throw cms::Exception ("Configuration")
            << "VisJet requires the IguanaService\n"
            "which is not present in the configuration file.\n"
            "You must add the service in the configuration file\n"
            "or remove the module that requires it";
    }
    
    edm::Handle<reco::CaloJetCollection> collection;
    event.getByLabel (inputTag_, collection);
    
    std::vector<edm::Handle<CaloTowerCollection> > towerCollections;
    event.getManyByType (towerCollections);

    if (collection.isValid ())
    {   
        IgDataStorage *storage = config->storage ();
        IgCollection &jets = storage->getCollection ("Jets_V1");

        IgProperty ET = jets.addProperty ("et", 0.0); 
        IgProperty ETA = jets.addProperty ("eta", 0.0);
        IgProperty THETA = jets.addProperty ("theta", 0.0);
        IgProperty PHI = jets.addProperty ("phi", 0.0);

        for (reco::CaloJetCollection::const_iterator it = collection->begin (), itEnd = collection->end (); it != itEnd; ++it)
        {
            IgCollectionItem ijet = jets.create ();
            ijet[ET]    = static_cast<double>((*it).et ());
            ijet[ETA]   = static_cast<double>((*it).eta ());
            ijet[THETA] = static_cast<double>(2.0 * atan (exp (-(*it).eta ())));
            ijet[PHI]   = static_cast<double>((*it).phi());
        }
    }    
    else 
    {
        // friendlyName:moduleLabel:instanceName:processName
        std::string error = "### Error: Jets "
                            + edm::TypeID (typeid (reco::CaloJetCollection)).friendlyClassName () + ":" 
                            + inputTag_.label() + ":"
                            + inputTag_.instance() + ":" 
                            + inputTag_.process() + " are not found.";

        IgDataStorage *storage = config->storage ();
        IgCollection &collection = storage->getCollection ("Errors_V1");
        IgProperty ERROR_MSG = collection.addProperty ("Error", std::string ());
        IgCollectionItem item = collection.create ();
        item [ERROR_MSG] = error;
    }
}

DEFINE_FWK_MODULE(VisJet);

---