DTSegment4DQuality.h

Go to the documentation of this file.

#ifndef Validation_DTSegment4D_H
#define Validation_DTSegment4D_H

#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "Histograms.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"

#include <vector>
#include <map>
#include <string>

namespace edm {
  class ParameterSet;
  class Event;
  class EventSetup;
}

class TFile;

class DTSegment4DQuality : public edm::EDAnalyzer {
public:
  DTSegment4DQuality(const edm::ParameterSet& pset);

  virtual ~DTSegment4DQuality();

  // Operations

  void analyze(const edm::Event & event, const edm::EventSetup& eventSetup);

  virtual void beginRun(const edm::Run& iRun, const edm::EventSetup &setup);
  
  // Write the histos to file
  void endJob();
  void endLuminosityBlock(edm::LuminosityBlock const& lumiSeg,
              edm::EventSetup const& c);

protected:

private: 

  // The file which will store the histos
  //TFile *theFile;
  // Switch for debug output
  bool debug;
  // Root file name
  std::string rootFileName;
  //Labels to read from event
  edm::InputTag simHitLabel;
  edm::InputTag segment4DLabel;
  edm::EDGetTokenT<edm::PSimHitContainer> simHitToken_;
  edm::EDGetTokenT<DTRecSegment4DCollection> segment4DToken_;
  //Sigma resolution on position
  double sigmaResX;
  double sigmaResY;
  //Sigma resolution on angle
  double sigmaResAlpha;
  double sigmaResBeta;

  HRes4DHit *h4DHit;
  HRes4DHit *h4DHit_W0;
  HRes4DHit *h4DHit_W1;
  HRes4DHit *h4DHit_W2;
  HRes4DHit *h4DHitWS[3][4];

  HEff4DHit *hEff_All;
  HEff4DHit *hEff_W0;
  HEff4DHit *hEff_W1;
  HEff4DHit *hEff_W2;
  HEff4DHit *hEffWS[3][4];

  DQMStore* dbe_;
  bool doall;
  bool local;
};

#endif