DTVDriftCalibration.h

Go to the documentation of this file.

#ifndef DTVDriftCalibration_H
#define DTVDriftCalibration_H

#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"

#include "CalibMuon/DTCalibration/interface/vDriftHistos.h"
#include "CalibMuon/DTCalibration/interface/DTTMax.h"
#include "CalibMuon/DTCalibration/interface/DTSegmentSelector.h"

#include "DTCalibrationMap.h"

#include <string>
#include <vector>

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

class TFile;
class DTMeanTimerFitter;
class DTGeometry;
class MuonGeometryRecord;

class DTVDriftCalibration : public edm::one::EDAnalyzer<> {
public:
  DTVDriftCalibration(const edm::ParameterSet& pset);

  ~DTVDriftCalibration() override;

  // Operations

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

  void endJob() override;

protected:
private:
  std::unique_ptr<DTSegmentSelector> select_;

  // The class containing TMax information
  typedef DTTMax::TMax TMax;

  dtcalibration::Histograms histograms_;

  // class to create/manage histos for each partition (SL)
  class cellInfo {
  public:
    cellInfo(const TString& name) { histos = new hTMaxCell(name); }

    ~cellInfo() { delete histos; }

    void add(const std::vector<const TMax*>& tMaxes);
    void update() { addedCells.clear(); }
    hTMaxCell* getHists() { return histos; }

  private:
    cellInfo() {}
    cellInfo(const cellInfo&) {}

    std::vector<dttmaxenums::TMaxCells> addedCells;
    hTMaxCell* histos;
  };

  TH1F* hChi2;
  h2DSegm* h2DSegmRZ;
  h2DSegm* h2DSegmRPhi;
  h4DSegm* h4DSegmAllCh;

  // Divide cellInfo by given granularity (to be implemented)
  // DTVDriftCalibration::cellInfo* partition(const DTWireId& wireId);

  // Specify the granularity for the TMax histograms
  enum TMaxGranularity { byChamber, bySL, byPartition };
  TMaxGranularity theGranularity;

  // The label used to retrieve 4D segments from the event
  edm::EDGetTokenT<DTRecSegment4DCollection> theRecHits4DToken;

  const edm::ESGetToken<DTGeometry, MuonGeometryRecord> theDTGeomToken;

  // Debug flag
  bool debug;

  // The label used to retrieve digis from the event
  std::string digiLabel;

  // The file which will contain the tMax histograms
  TFile* theFile;

  // The fitter
  std::unique_ptr<DTMeanTimerFitter> theFitter;

  // Perform the vDrift and t0 evaluation or just fill the
  //  tMaxHists (if you read the dataset in different jobs)
  bool findVDriftAndT0;

  // The name of the output text file
  std::string theVDriftOutputFile;

  // Map of wires and cellInfo with coarse granularity
  std::map<DTWireId, cellInfo*> theWireIdAndCellMap;

  // Switch for checking of noisy channels
  //bool checkNoisyChannels;

  // The module for t0 subtraction
  std::unique_ptr<DTTTrigBaseSync> theSync;  //FIXME: should be const

  // parameter set for DTCalibrationMap constructor
  edm::ParameterSet theCalibFilePar;

  // Maximum value for the 4D Segment chi2
  //double theMaxChi2;

  // Maximum incident angle for Phi Seg
  //double theMaxPhiAngle;

  // Maximum incident angle for Theta Seg
  //double theMaxZAngle;

  // Choose the chamber you want to calibrate
  std::string theCalibChamber;

  // which format to be created
  bool writeLegacyVDriftDB;
};
#endif