---

VisDTDigi Class Reference

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

Inheritance diagram for VisDTDigi:

List of all members.

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

---

Detailed Description

Definition at line 8 of file VisDTDigi.h.

---

Constructor & Destructor Documentation

VisDTDigi::VisDTDigi

(

const edm::ParameterSet &

iConfig

)

[explicit]

Definition at line 21 of file VisDTDigi.cc.

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

virtual VisDTDigi::~VisDTDigi

(

void

)

[inline, virtual]

Definition at line 12 of file VisDTDigi.h.

        {}

---

Member Function Documentation

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

Implements edm::EDAnalyzer.

Definition at line 25 of file VisDTDigi.cc.

References IgCollection::addProperty(), DTTopology::cellHeight(), DTTopology::cellLenght(), DTTopology::cellWidth(), DTLayer::chamber(), DTSuperLayerId::chamberId(), collection, edmplugin::standard::config(), IgCollection::create(), error, Exception, edm::EventSetup::get(), IgDataStorage::getCollection(), DTChamber::id(), inputTag_, edm::InputTag::instance(), edm::Service< T >::isAvailable(), edm::Handle< T >::isValid(), edm::ESHandle< T >::isValid(), edm::InputTag::label(), DTLayerId::layer(), edm::InputTag::process(), range, DTChamberId::sector(), DTLayer::specificTopology(), DTChamberId::station(), DTSuperLayerId::superLayer(), DTLayerId::superlayerId(), GeomDet::surface(), Surface::toGlobal(), DTChamberId::wheel(), DTTopology::wirePosition(), 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")
            << "VisDTDigi 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<DTDigiCollection> collection;
    event.getByLabel (inputTag_, collection);

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

        IgProperty WIREN = digis.addProperty("wireNumber", int(0));
        IgProperty LAYER_ID = digis.addProperty("layerId", int(0));
        IgProperty SUPERLAYER_ID = digis.addProperty("superLayerId", int(0));
        IgProperty SECTOR_ID = digis.addProperty("sectorId", int(0));
        IgProperty STATION_ID = digis.addProperty("stationId", int(0));
        IgProperty WHEEL_ID = digis.addProperty("wheelId", int(0));

        IgProperty POS = digis.addProperty("pos", IgV3d());
        //IgProperty AXIS = digis.addProperty("axis", IgV3d());
        //IgProperty ANGLE = digis.addProperty("angle", 0.0);
        IgProperty COUNT = digis.addProperty("countsTDC", int(0));
        IgProperty NUMBER = digis.addProperty("number", int(0));

        IgProperty CELL_WIDTH = digis.addProperty("cellWidth", 0.0);
        IgProperty CELL_LENGTH = digis.addProperty("cellLength", 0.0);
        IgProperty CELL_HEIGHT = digis.addProperty("cellHeight", 0.0);
        
        for ( DTDigiCollection::DigiRangeIterator dri = collection->begin();
              dri != collection->end(); ++dri )
        {
            const DTLayerId& dtlayerId = (*dri).first;
            const DTDigiCollection::Range& range = (*dri).second;

            const DTLayer* layer = geom->layer(dtlayerId);

            const DTChamber* chamber = layer->chamber();
            const DTChamberId chamberId = chamber->id();

            //const GeomDetUnit* det = geom->idToDetUnit((*dri).first);

            int layerId = dtlayerId.layer();
            int superLayerId = dtlayerId.superlayerId().superLayer();
            int sectorId = dtlayerId.superlayerId().chamberId().sector();                 
            int stationId = dtlayerId.superlayerId().chamberId().station();
            int wheelId = dtlayerId.superlayerId().chamberId().wheel();

            for ( DTDigiCollection::const_iterator dit = range.first;
                  dit != range.second; ++dit )
            {
                IgCollectionItem digi = digis.create();

                digi[LAYER_ID] = static_cast<int>(layerId);
                digi[SUPERLAYER_ID] = static_cast<int>(superLayerId);
                digi[SECTOR_ID] = static_cast<int>(sectorId);
                digi[STATION_ID] = static_cast<int>(stationId);
                digi[WHEEL_ID] = static_cast<int>(wheelId);
           
                int wireNumber = (*dit).wire();
                digi[WIREN] = static_cast<int>(wireNumber);

                double localXPos = layer->specificTopology().wirePosition(wireNumber);
                LocalPoint localPos(localXPos, 0.0, 0.0);
                GlobalPoint pos = layer->surface().toGlobal(localPos);
            
                digi[POS] = IgV3d(static_cast<double>(pos.x()/100.0),
                                  static_cast<double>(pos.y()/100.0),
                                  static_cast<double>(pos.z()/100.0));

                //NOTE: May need to hand-code transformation by hand.
                // Leave this for next version.
                /*
                float angle;
                SbVec3f axis;
                VisTrackingGeometry::createRotation(det, axis, angle);
                */

                int countsTDC = (*dit).countsTDC();
                digi[COUNT] = static_cast<int>(countsTDC);
                int number = (*dit).number();
                digi[NUMBER] = static_cast<int>(number);
                
                digi[CELL_WIDTH] =      
                    static_cast<double>(layer->specificTopology().cellWidth()/100.0);
                digi[CELL_LENGTH] = 
                    static_cast<double>(layer->specificTopology().cellLenght()/100.0);  
                digi[CELL_HEIGHT] = 
                    static_cast<double>(layer->specificTopology().cellHeight()/100.0);
            }
        }
        
    }
    
    else 
    {
        // friendlyName:moduleLabel:instanceName:processName
        std::string error = "### Error: DTDigis "
                            + edm::TypeID (typeid (DTDigiCollection)).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 VisDTDigi::inputTag_ [private]

Definition at line 17 of file VisDTDigi.h.

Referenced by analyze().

---

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

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

---