SiStripLAProfileBooker.h

Go to the documentation of this file.

#ifndef CalibTracker_SiStripLorentzAngle_SiStripLAProfileBooker_h
#define CalibTracker_SiStripLorentzAngle_SiStripLAProfileBooker_h
 
#include <map>
 
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
 
#include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
#include "MagneticField/Engine/interface/MagneticField.h"
 
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"
 
#include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
#include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
 
#include "DQMServices/Core/interface/DQMStore.h"
 
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
 
#include <TTree.h>
#include <TFile.h>
#include <TH1D.h>
#include <TDirectory.h>
 
class TrackerTopology;
 
class SiStripLAProfileBooker : public edm::EDAnalyzer {
public:
  typedef dqm::legacy::MonitorElement MonitorElement;
  typedef dqm::legacy::DQMStore DQMStore;
  explicit SiStripLAProfileBooker(const edm::ParameterSet& conf);
 
  ~SiStripLAProfileBooker() override;
 
  void beginRun(edm::Run const&, const edm::EventSetup& c) override;
 
  void endJob() override;
 
  void analyze(const edm::Event& e, const edm::EventSetup& c) override;
 
  void getlayer(const DetId& detid, const TrackerTopology* tTopo, std::string& name, unsigned int& layerid);
 
private:
  typedef struct {
    float thickness;
    float pitch;
    LocalVector magfield;
  } detparameters;
  typedef std::map<unsigned int, detparameters*> detparmap;
  typedef std::map<unsigned int, MonitorElement*> histomap;
 
  int trackcollsize, trajsize, RunNumber, EventNumber, ClSize, HitCharge, Type, Layer, Wheel, bw_fw, Ext_Int,
      MonoStereo, ParticleCharge;
  float sumx, hit_std_dev, barycenter, TanTrackAngle, TanTrackAngleParallel, SignCorrection, MagField, XGlobal, YGlobal,
      ZGlobal, Momentum, pt, chi2norm, EtaTrack, PhiTrack;
  int nstrip, eventcounter, size, HitNr, hitcounter, hitcounter_2ndloop, worse_double_hit, better_double_hit,
      HitPerTrack;
  int id_detector, TrackCounter, EventCounter;
  float thick_detector, pitch_detector;
  uint8_t Amplitudes[100];
 
  TTree *HitsTree, *TrackTree, *EventTree;
  TFile* hFile;
 
  TDirectory *Hit_Tree, *Track_Tree, *Event_Tree;
 
  histomap histos;
  histomap summaryhisto;
 
  DQMStore* dbe_;
 
  detparmap detmap;
  detparmap summarydetmap;
  edm::ParameterSet conf_;
  std::string treename_;
 
  const TrackerGeometry* tkGeom_ = nullptr;
  edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> tTopoToken_;
  edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkGeomToken_;
  edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> magFieldToken_;
  edm::ESGetToken<SiStripDetCabling, SiStripDetCablingRcd> detCablingToken_;
};
 
#endif