---

VisTrackingRecHit.cc

Go to the documentation of this file.

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

#include "VisReco/Analyzer/interface/VisTrackingRecHit.h"
#include "VisReco/Analyzer/interface/VisLocalPosition.h"
#include "VisReco/Analyzer/interface/IguanaService.h"
#include "DataFormats/TrackingRecHit/interface/TrackingRecHitFwd.h"
#include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
#include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
#include "Iguana/Framework/interface/IgCollection.h"

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

VisTrackingRecHit::VisTrackingRecHit (const edm::ParameterSet& iConfig)
    : inputTag_ (iConfig.getParameter<edm::InputTag>("visTrackingRecHitTag"))
{}

void
VisTrackingRecHit::analyze( const edm::Event& event, const edm::EventSetup& eventSetup)
{
    edm::Service<IguanaService> config;
    if (! config.isAvailable ()) 
    {
        throw cms::Exception ("Configuration")
            << "VisTrackingRecHit 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<TrackingRecHitCollection> collection;
    event.getByLabel (inputTag_, collection);
 
    edm::ESHandle<GlobalTrackingGeometry> geom;
    eventSetup.get<GlobalTrackingGeometryRecord> ().get (geom);

    if (collection.isValid () && geom.isValid ())
    {       
        IgDataStorage *storage = config->storage ();
        IgCollection &recHits = storage->getCollection("TrackingRecHits_V1");
        IgProperty POS = recHits.addProperty("pos", IgV3d());

        for (TrackingRecHitCollection::const_iterator it=collection->begin(), itEnd=collection->end(); it!=itEnd; ++it)
        {
            if ((*it).isValid () && !(*it).geographicalId ().null ())
            {
                LocalPoint point = VisLocalPosition::localPosition(&(*it), geom.product ());
                float x = geom->idToDet ((*it).geographicalId ())->surface ().toGlobal (point).x () / 100.0;
                float y = geom->idToDet ((*it).geographicalId ())->surface ().toGlobal (point).y () / 100.0;
                float z = geom->idToDet ((*it).geographicalId ())->surface ().toGlobal (point).z () / 100.0;

                IgCollectionItem irechit = recHits.create();
                irechit[POS] = IgV3d (static_cast<double>(x), static_cast<double>(y), static_cast<double>(z));
            }       
        }
    }
    else 
    {
        // friendlyName:moduleLabel:instanceName:processName
        std::string error = "### Error: TrackingRecHits "
                            + edm::TypeID (typeid (TrackingRecHitCollection)).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(VisTrackingRecHit);

---