|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
|
#ifndef CalibTracker_SiStripLorentzAngle_SiStripLAProfileBooker_h
#define CalibTracker_SiStripLorentzAngle_SiStripLAProfileBooker_h
#include <map>
#include "FWCore/Framework/interface/one/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::one::EDAnalyzer<edm::one::WatchRuns> {
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 endRun(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
|