SiPixelRecHitsValid.h

Go to the documentation of this file.

#ifndef SiPixelRecHitsValid_h
#define SiPixelRecHitsValid_h

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/EDProduct.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

//DWM histogram services
#include "DQMServices/Core/interface/DQMStore.h"

#include "FWCore/ServiceRegistry/interface/Service.h"

//Simhit stuff
#include "SimDataFormats/TrackingHit/interface/PSimHit.h"
#include "SimTracker/TrackerHitAssociation/interface/TrackerHitAssociator.h"

#include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"
#include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHitCollection.h"
#include "DataFormats/DetId/interface/DetId.h"

#include "Geometry/CommonTopologies/interface/PixelTopology.h"
#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetType.h"
#include "Geometry/CommonDetUnit/interface/GeomDetType.h"
#include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"

#include "FWCore/Utilities/interface/InputTag.h"

#include <string>
#include "DQMServices/Core/interface/MonitorElement.h"


class SiPixelRecHitsValid : public edm::EDAnalyzer {

   public:
    //Constructor
    SiPixelRecHitsValid(const edm::ParameterSet& conf);

    //Destructor
    ~SiPixelRecHitsValid();

   protected:

    virtual void analyze(const edm::Event& e, const edm::EventSetup& c);
    void beginJob();
    void endJob();

   private:
    DQMStore* dbe_;
    std::string outputFile_;

    edm::ParameterSet conf_;

    void fillBarrel(const SiPixelRecHit &,const PSimHit &, DetId, const PixelGeomDetUnit *);    
    void fillForward(const SiPixelRecHit &, const PSimHit &, DetId, const PixelGeomDetUnit *);

    //Clusters BPIX
    MonitorElement* clustYSizeModule[8];
    MonitorElement* clustXSizeLayer[3];
    MonitorElement* clustChargeLayer1Modules[8];
    MonitorElement* clustChargeLayer2Modules[8];
    MonitorElement* clustChargeLayer3Modules[8];

    //Cluster FPIX
    MonitorElement* clustXSizeDisk1Plaquettes[7];
    MonitorElement* clustXSizeDisk2Plaquettes[7];
    MonitorElement* clustYSizeDisk1Plaquettes[7];
    MonitorElement* clustYSizeDisk2Plaquettes[7];
    MonitorElement* clustChargeDisk1Plaquettes[7];
    MonitorElement* clustChargeDisk2Plaquettes[7];

    //RecHits BPIX
    MonitorElement* recHitXResAllB;
    MonitorElement* recHitYResAllB;
    MonitorElement* recHitXFullModules;
    MonitorElement* recHitXHalfModules;
    MonitorElement* recHitYAllModules;
    MonitorElement* recHitXResFlippedLadderLayers[3];
    MonitorElement* recHitXResNonFlippedLadderLayers[3];
    MonitorElement* recHitYResLayer1Modules[8];
    MonitorElement* recHitYResLayer2Modules[8];
    MonitorElement* recHitYResLayer3Modules[8];

    //RecHits FPIX
    MonitorElement* recHitXResAllF;
    MonitorElement* recHitYResAllF;
    MonitorElement* recHitXPlaquetteSize1;
    MonitorElement* recHitXPlaquetteSize2;
    MonitorElement* recHitYPlaquetteSize2;
    MonitorElement* recHitYPlaquetteSize3;
    MonitorElement* recHitYPlaquetteSize4;
    MonitorElement* recHitYPlaquetteSize5;
    MonitorElement* recHitXResDisk1Plaquettes[7];
    MonitorElement* recHitXResDisk2Plaquettes[7];
    MonitorElement* recHitYResDisk1Plaquettes[7];
    MonitorElement* recHitYResDisk2Plaquettes[7];

    // Pull distributions
    //RecHits BPIX
    MonitorElement* recHitXPullAllB;
    MonitorElement* recHitYPullAllB;

    MonitorElement* recHitXPullFlippedLadderLayers[3];
    MonitorElement* recHitXPullNonFlippedLadderLayers[3];
    MonitorElement* recHitYPullLayer1Modules[8];
    MonitorElement* recHitYPullLayer2Modules[8];
    MonitorElement* recHitYPullLayer3Modules[8];

    //RecHits FPIX
    MonitorElement* recHitXPullAllF;
    MonitorElement* recHitYPullAllF;

    MonitorElement* recHitXPullDisk1Plaquettes[7];
    MonitorElement* recHitXPullDisk2Plaquettes[7];
    MonitorElement* recHitYPullDisk1Plaquettes[7];
    MonitorElement* recHitYPullDisk2Plaquettes[7];

        edm::InputTag src_;
};

#endif