---

VisDTRecHit Class Reference

#include <VisReco/Analyzer/interface/VisDTRecHit.h>

Inheritance diagram for VisDTRecHit:

List of all members.

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
	VisDTRecHit (const edm::ParameterSet &)
virtual	~VisDTRecHit (void)
Private Attributes
edm::InputTag	inputTag_

---

Detailed Description

Definition at line 8 of file VisDTRecHit.h.

---

Constructor & Destructor Documentation

VisDTRecHit::VisDTRecHit

(

const edm::ParameterSet &

iConfig

)

[explicit]

Definition at line 22 of file VisDTRecHit.cc.

    : inputTag_(iConfig.getParameter<edm::InputTag>("visDTRecHitTag"))
{}

virtual VisDTRecHit::~VisDTRecHit

(

void

)

[inline, virtual]

Definition at line 12 of file VisDTRecHit.h.

        {}

---

Member Function Documentation

void VisDTRecHit::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)			[virtual]

Implements edm::EDAnalyzer.

Definition at line 26 of file VisDTRecHit.cc.

References IgCollection::addProperty(), DTTopology::cellLenght(), collection, edmplugin::standard::config(), IgCollection::create(), error, Exception, edm::EventSetup::get(), IgDataStorage::getCollection(), inputTag_, edm::InputTag::instance(), edm::Service< T >::isAvailable(), edm::Handle< T >::isValid(), edm::ESHandle< T >::isValid(), edm::InputTag::label(), DTRecHit1D::localPosition(), edm::InputTag::process(), DTLayer::specificTopology(), GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{
    edm::Service<IguanaService> config;
    if (! config.isAvailable ()) 
    {
        throw cms::Exception ("Configuration")
            << "VisDTRecHit 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::ESHandle<DTGeometry> geom;
    eventSetup.get<MuonGeometryRecord>().get(geom);

    edm::Handle<DTRecHitCollection> collection;
    event.getByLabel (inputTag_, collection);

    if ( collection.isValid () && geom.isValid ())
    {
        IgDataStorage* storage = config->storage();
        IgCollection& recHits = storage->getCollection("DTRecHits_V1");

        IgProperty WIRE_ID = recHits.addProperty("wireId", int(0));
        IgProperty LAYER_ID = recHits.addProperty("layerId", int(0));
        IgProperty SUPERLAYER_ID = recHits.addProperty("superLayerId", int(0));
        IgProperty SECTOR_ID = recHits.addProperty("sectorId", int(0));
        IgProperty STATION_ID = recHits.addProperty("stationId", int(0));
        IgProperty WHEEL_ID = recHits.addProperty("wheelId", int(0));
        IgProperty DIGITIME = recHits.addProperty("digitime", 0.0);

        IgProperty LPLUS_GLOBALPOS = recHits.addProperty("lPlusGlobalPos", IgV3d());
        IgProperty LMINUS_GLOBALPOS = recHits.addProperty("lMinusGlobalPos", IgV3d());
        IgProperty RPLUS_GLOBALPOS = recHits.addProperty("rPlusGlobalPos", IgV3d());
        IgProperty RMINUS_GLOBALPOS = recHits.addProperty("rMinusGlobalPos", IgV3d());
        IgProperty LGLOBALPOS = recHits.addProperty("lGlobalPos", IgV3d());
        IgProperty RGLOBALPOS = recHits.addProperty("rGlobalPos", IgV3d());

        for ( DTRecHitCollection::const_iterator dit = collection->begin();
              dit != collection->end(); ++dit )
        {
            IgCollectionItem recHit = recHits.create();

            const DTRecHit1D* lrechit = (*dit).componentRecHit(Left);
            const DTRecHit1D* rrechit = (*dit).componentRecHit(Right);
            
            int wireId = (*dit).wireId().wire();
            recHit[WIRE_ID] = static_cast<int>(wireId);
            
            int layerId = (*dit).wireId().layerId().layer();
            recHit[LAYER_ID] = static_cast<int>(layerId);
             
            int superLayerId = (*dit).wireId().layerId().superlayerId().superLayer();
            recHit[SUPERLAYER_ID] = static_cast<int>(superLayerId);
             
            int sectorId = (*dit).wireId().layerId().superlayerId().chamberId().sector();
            recHit[SECTOR_ID] = static_cast<int>(sectorId);
             
            int stationId = (*dit).wireId().layerId().superlayerId().chamberId().station();
            recHit[STATION_ID] = static_cast<int>(stationId);
             
            int wheelId = (*dit).wireId().layerId().superlayerId().chamberId().wheel();
            recHit[WHEEL_ID] = static_cast<int>(wheelId);
            
            double digitime = (*dit).digiTime();
            recHit[DIGITIME] = static_cast<double>(digitime);

            const DTLayer* layer = geom->layer((*dit).wireId());  
            const DTTopology& topo = layer->specificTopology();
               
            LocalPoint lLocalPos = lrechit->localPosition();
            LocalPoint rLocalPos = rrechit->localPosition();
               
            LocalPoint yPlusLeft(lLocalPos.x(),topo.cellLenght()/2.,0.);
            LocalPoint yMinusLeft(lLocalPos.x(),-topo.cellLenght()/2.,0.);
            LocalPoint yPlusRight(rLocalPos.x(),topo.cellLenght()/2.,0.);
            LocalPoint yMinusRight(rLocalPos.x(),-topo.cellLenght()/2.,0.);
               
            GlobalPoint lPlusGlobalPos = layer->toGlobal(yPlusLeft);
            GlobalPoint lMinusGlobalPos = layer->toGlobal(yMinusLeft);
            GlobalPoint rPlusGlobalPos = layer->toGlobal(yPlusRight);
            GlobalPoint rMinusGlobalPos = layer->toGlobal(yMinusRight);
               
            GlobalPoint lGlobalPos = layer->toGlobal(lLocalPos);
            GlobalPoint rGlobalPos = layer->toGlobal(rLocalPos);

            recHit[LPLUS_GLOBALPOS] = IgV3d(static_cast<double>(lPlusGlobalPos.x()/100.0),
                                            static_cast<double>(lPlusGlobalPos.y()/100.0),
                                            static_cast<double>(lPlusGlobalPos.z()/100.0));

            recHit[LMINUS_GLOBALPOS] = IgV3d(static_cast<double>(lMinusGlobalPos.x()/100.0),
                                              static_cast<double>(lMinusGlobalPos.y()/100.0),
                                              static_cast<double>(lMinusGlobalPos.z()/100.0));
             
            recHit[RPLUS_GLOBALPOS] = IgV3d(static_cast<double>(rPlusGlobalPos.x()/100.0),
                                             static_cast<double>(rPlusGlobalPos.y()/100.0),
                                             static_cast<double>(rPlusGlobalPos.z()/100.0));

            recHit[RMINUS_GLOBALPOS] = IgV3d(static_cast<double>(rMinusGlobalPos.x()/100.0),
                                             static_cast<double>(rMinusGlobalPos.y()/100.0),
                                             static_cast<double>(rMinusGlobalPos.z()/100.0));

            recHit[LGLOBALPOS] = IgV3d(static_cast<double>(lGlobalPos.x()/100.0),
                                       static_cast<double>(lGlobalPos.y()/100.0),
                                       static_cast<double>(lGlobalPos.z()/100.0));

            recHit[RGLOBALPOS] = IgV3d(static_cast<double>(rGlobalPos.x()/100.0),
                                        static_cast<double>(rGlobalPos.y()/100.0),
                                        static_cast<double>(rGlobalPos.z()/100.0));
        }
        
    }
    
    else 
    {
        // friendlyName:moduleLabel:instanceName:processName
        std::string error = "### Error: DTRecHits "
                            + edm::TypeID (typeid (DTRecHitCollection)).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;
    }
}

---

Member Data Documentation

edm::InputTag VisDTRecHit::inputTag_ [private]

Definition at line 17 of file VisDTRecHit.h.

Referenced by analyze().

---

The documentation for this class was generated from the following files:

• VisReco/Analyzer/interface/VisDTRecHit.h
• VisReco/Analyzer/src/VisDTRecHit.cc

---