---

VisPFRecHit Class Reference

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

Inheritance diagram for VisPFRecHit:

List of all members.

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

---

Detailed Description

Definition at line 8 of file VisPFRecHit.h.

---

Constructor & Destructor Documentation

VisPFRecHit::VisPFRecHit

(

const edm::ParameterSet &

iConfig

)

Definition at line 22 of file VisPFRecHit.cc.

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

virtual VisPFRecHit::~VisPFRecHit

(

)

[inline, virtual]

Definition at line 12 of file VisPFRecHit.h.

        {}

---

Member Function Documentation

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

Implements edm::EDAnalyzer.

Definition at line 26 of file VisPFRecHit.cc.

References IgCollection::addProperty(), IgAssociationSet::associate(), collection, edmplugin::standard::config(), IgCollection::create(), PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, relval_parameters_module::energy, error, ETA, eta, Exception, edm::EventSetup::get(), IgDataStorage::getAssociationSet(), IgDataStorage::getCollection(), i, inputTag_, edm::InputTag::instance(), edm::Service< T >::isAvailable(), edm::Handle< T >::isValid(), edm::ESHandle< T >::isValid(), edm::InputTag::label(), sistrip::LAYER, PHI, phi, edm::InputTag::process(), and EZArrayFL< T >::size().

{
    edm::Service<IguanaService> config;
    
    if ( ! config.isAvailable() ) 
    {
        throw cms::Exception("Configuration")
            << "VisPFRecHit 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::PFRecHitCollection> collection;
    event.getByLabel(inputTag_, collection);

    edm::ESHandle<CaloGeometry> caloGeometry;
    eventSetup.get<CaloGeometryRecord>().get(caloGeometry);

    if ( collection.isValid() && caloGeometry.isValid() )
    {
        IgDataStorage *storage = config->storage();
        
        IgCollection& rechits = storage->getCollection("PFRecHits_V1");

        IgProperty ID = rechits.addProperty("id", 0);
        IgProperty NAME = rechits.addProperty("name", std::string());

        IgProperty E = rechits.addProperty("energy", 0.0); 
        IgProperty ETA = rechits.addProperty("eta", 0.0);
        IgProperty PHI = rechits.addProperty("phi", 0.0);

        IgProperty LAYER = rechits.addProperty("layer", std::string());

        IgCollection& corners = storage->getCollection("Corners_V1");

        IgProperty F1 = corners.addProperty("front1", IgV3d());
        IgProperty F2 = corners.addProperty("front2", IgV3d());
        IgProperty F3 = corners.addProperty("front3", IgV3d());
        IgProperty F4 = corners.addProperty("front4", IgV3d());
        
        IgProperty B1 = corners.addProperty("back1", IgV3d());
        IgProperty B2 = corners.addProperty("back2", IgV3d());
        IgProperty B3 = corners.addProperty("back3", IgV3d());
        IgProperty B4 = corners.addProperty("back4", IgV3d());

        IgAssociationSet& recHitCorners = storage->getAssociationSet("PFRecHitCorners_V1");

        unsigned int i = 0;

        for ( std::vector<reco::PFRecHit>::const_iterator rechit = collection->begin();
                 rechit != collection->end(); ++rechit, ++i )
            {
                IgCollectionItem rh = rechits.create();

                rh[ID] = i;
                rh["name"] = std::string("PFRecHit");

                if ( (*rechit).layer() == PFLayer::ECAL_BARREL )
                    rh[LAYER] = std::string("EB");
                
                else if ( (*rechit).layer() == PFLayer::ECAL_ENDCAP )
                    rh[LAYER] = std::string("EE");

                // NOTE: Add HCAL support. Is there still the detector ID problem in PF with HCAL?

                else
                    rh[LAYER] = std::string("UNDEFINED");
                   
                double energy = (*rechit).energy();
                rh["energy"] = energy;

                double eta = (*rechit).position().Eta();
                double phi = (*rechit).position().Phi();

                rh[ETA] = eta;
                rh[PHI] = phi;
                
                // NOTE: Which corners are which depends on whether it's ECAL or HCAL.
                // Fix this later.

                double f1x, f1y, f1z;
                double f2x, f2y, f2z;
                double f3x, f3y, f3z;
                double f4x, f4y, f4z;
                
                double b1x, b1y, b1z;
                double b2x, b2y, b2z;
                double b3x, b3y, b3z;
                double b4x, b4y, b4z;
                
                const CaloCellGeometry::CornersVec& cornerPoints 
                    = (*caloGeometry).getGeometry((*rechit).detId())->getCorners();
                
                assert(cornerPoints.size() == 8);
                
                IgCollectionItem cr = corners.create();

                f1x = cornerPoints[3].x()/100.0;
                f1y = cornerPoints[3].y()/100.0;
                f1z = cornerPoints[3].z()/100.0;

                cr[F1] = IgV3d(f1x, f1y, f1z);
                
                f2x = cornerPoints[2].x()/100.0;
                f2y = cornerPoints[2].y()/100.0;
                f2z = cornerPoints[2].z()/100.0;

                cr[F2] = IgV3d(f2x, f2y, f2z);

                f3x = cornerPoints[1].x()/100.0;
                f3y = cornerPoints[1].y()/100.0;
                f3z = cornerPoints[1].z()/100.0;

                cr[F3] = IgV3d(f3x, f3y, f3z);

                f4x = cornerPoints[0].x()/100.0;
                f4y = cornerPoints[0].y()/100.0;
                f4z = cornerPoints[0].z()/100.0;

                cr[F4] = IgV3d(f4x, f4y, f4z);

                b1x = cornerPoints[7].x()/100.0;
                b1y = cornerPoints[7].y()/100.0;
                b1z = cornerPoints[7].z()/100.0;

                cr["back1"] = IgV3d(b1x, b1y, b1z);

                b2x = cornerPoints[6].x()/100.0;
                b2y = cornerPoints[6].y()/100.0;
                b2z = cornerPoints[6].z()/100.0;

                cr["back2"] = IgV3d(b2x, b2y, b2z);

                b3x = cornerPoints[5].x()/100.0;
                b3y = cornerPoints[5].y()/100.0;
                b3z = cornerPoints[5].z()/100.0;

                cr["back3"] = IgV3d(b3x, b3y, b3z);

                b4x = cornerPoints[4].x()/100.0;
                b4y = cornerPoints[4].y()/100.0;
                b4z = cornerPoints[4].z()/100.0;

                cr["back4"] = IgV3d(b4x, b4y, b4z);

                recHitCorners.associate(rh, cr);
            }    
        }
    
    else 
    {
        //std::cout<<"Collection not valid. No PFRecHits?!"<<std::endl;
        
        // friendlyName:moduleLabel:instanceName:processName
        std::string error = "### Error: PFRecHits "
                            + edm::TypeID (typeid (reco::PFRecHitCollection)).friendlyClassName() + ":" 
                            + inputTag_.label() + ":"
                            + inputTag_.instance() + ":" 
                            + inputTag_.process() + " are not found.";

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

---

Member Data Documentation

edm::InputTag VisPFRecHit::inputTag_ [private]

Definition at line 18 of file VisPFRecHit.h.

Referenced by analyze().

---

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

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

---