CMS 3D CMS Logo
MultiplicityCorrelatorHistogramMaker.cc
Go to the documentation of this file.
1 #include "DPGAnalysis/SiStripTools/interface/MultiplicityCorrelatorHistogramMaker.h"
2 #include "FWCore/Framework/interface/Event.h"
3 #include "FWCore/Framework/interface/Run.h"
4 #include "FWCore/MessageLogger/interface/MessageLogger.h"
5 #include "FWCore/ParameterSet/interface/ParameterSet.h"
6 #include "FWCore/ServiceRegistry/interface/Service.h"
7 #include "CommonTools/UtilAlgos/interface/TFileService.h"
8 #include "TH2F.h"
9 #include "TH1F.h"
10 #include <cmath>
11 
12 MultiplicityCorrelatorHistogramMaker::MultiplicityCorrelatorHistogramMaker(edm::ConsumesCollector&& iC):
13  m_rhm(iC, false), m_fhm(iC, true), m_runHisto(false), m_runHistoBXProfile(false), m_runHistoBX(false), m_runHisto2D(false), m_runHistoProfileBX(false),
14  m_scfact(1.), m_yvsxmult(nullptr),
15  m_atanyoverx(nullptr), m_atanyoverxrun(nullptr), m_atanyoverxvsbxrun(nullptr), m_atanyoverxvsbxrun2D(nullptr),
16  m_yvsxmultrun(nullptr), m_yvsxmultprofvsbxrun(nullptr), m_xvsymultprofvsbxrun(nullptr)
17 {}
18 
19 MultiplicityCorrelatorHistogramMaker::MultiplicityCorrelatorHistogramMaker(const edm::ParameterSet& iConfig, edm::ConsumesCollector&& iC):
20  m_rhm(iC, false), m_fhm(iC, true),
21  m_runHisto(iConfig.getParameter<bool>("runHisto")),
22  m_runHistoBXProfile(iConfig.getParameter<bool>("runHistoBXProfile")),
23  m_runHistoBX(iConfig.getParameter<bool>("runHistoBX")),
24  m_runHisto2D(iConfig.getParameter<bool>("runHisto2D")),
25  m_runHistoProfileBX(iConfig.getUntrackedParameter<bool>("runHistoProfileBX",false)),
26  m_scfact(iConfig.getUntrackedParameter<double>("scaleFactor",1.)),
27  m_atanyoverxrun(nullptr), m_atanyoverxvsbxrun(nullptr), m_atanyoverxvsbxrun2D(nullptr), m_yvsxmultrun(nullptr),
28  m_yvsxmultprofvsbxrun(nullptr), m_xvsymultprofvsbxrun(nullptr)
29  {
30 
31  edm::Service<TFileService> tfserv;
32 
33  char hname[300];
34  sprintf(hname,"%sVs%s",
35  iConfig.getParameter<std::string>("yDetLabel").c_str(),
36  iConfig.getParameter<std::string>("xDetLabel").c_str());
37  char htitle[300];
38  sprintf(htitle,"%s Vs %s multiplicity",
39  iConfig.getParameter<std::string>("yDetLabel").c_str(),
40  iConfig.getParameter<std::string>("xDetLabel").c_str());
41 
42  m_yvsxmult = tfserv->make<TH2F>(hname,htitle,
43  iConfig.getParameter<unsigned int>("xBins"),0.,iConfig.getParameter<double>("xMax"),
44  iConfig.getParameter<unsigned int>("yBins"),0.,iConfig.getParameter<double>("yMax"));
45 
46  if(m_runHisto && m_runHisto2D) {
47  m_yvsxmultrun = m_rhm.makeTH2F(hname,htitle,
48  iConfig.getParameter<unsigned int>("xBins"),0.,iConfig.getParameter<double>("xMax"),
49  iConfig.getParameter<unsigned int>("yBins"),0.,iConfig.getParameter<double>("yMax"));
50  }
51 
52  if(m_runHisto && m_runHistoProfileBX) {
53  sprintf(hname,"%sVs%sprofvsbx",
54  iConfig.getParameter<std::string>("yDetLabel").c_str(),
55  iConfig.getParameter<std::string>("xDetLabel").c_str());
56  sprintf(htitle,"%s Vs %s multiplicity vs BX",
57  iConfig.getParameter<std::string>("yDetLabel").c_str(),
58  iConfig.getParameter<std::string>("xDetLabel").c_str());
59  m_yvsxmultprofvsbxrun = m_fhm.makeTProfile2D(hname,htitle,
60  3564,-0.5,3564-0.5,
61  iConfig.getParameter<unsigned int>("xBins"),0.,iConfig.getParameter<double>("xMax"));
62  sprintf(hname,"%sVs%sprofvsbx",
63  iConfig.getParameter<std::string>("xDetLabel").c_str(),
64  iConfig.getParameter<std::string>("yDetLabel").c_str());
65  sprintf(htitle,"%s Vs %s multiplicity vs BX",
66  iConfig.getParameter<std::string>("xDetLabel").c_str(),
67  iConfig.getParameter<std::string>("yDetLabel").c_str());
68  m_xvsymultprofvsbxrun = m_fhm.makeTProfile2D(hname,htitle,
69  3564,-0.5,3564-0.5,
70  iConfig.getParameter<unsigned int>("yBins"),0.,iConfig.getParameter<double>("yMax"));
71  }
72 
73  sprintf(hname,"%sOver%s",
74  iConfig.getParameter<std::string>("yDetLabel").c_str(),
75  iConfig.getParameter<std::string>("xDetLabel").c_str());
76  sprintf(htitle,"atan (%4.2f*%s / %s multiplicity ratio)",
77  m_scfact,
78  iConfig.getParameter<std::string>("yDetLabel").c_str(),
79  iConfig.getParameter<std::string>("xDetLabel").c_str()
80  );
81 
82  m_atanyoverx = tfserv->make<TH1F>(hname,htitle,
83  iConfig.getParameter<unsigned int>("rBins"),0.,1.6);
84 
85  if(m_runHisto) {
86  sprintf(hname,"%sOver%srun",
87  iConfig.getParameter<std::string>("yDetLabel").c_str(),
88  iConfig.getParameter<std::string>("xDetLabel").c_str());
89  m_atanyoverxrun = m_rhm.makeTH1F(hname,htitle,
90  iConfig.getParameter<unsigned int>("rBins"),0.,1.6);
91  if(m_runHistoBX) {
92 
93  sprintf(hname,"%sOver%svsbx2D",
94  iConfig.getParameter<std::string>("yDetLabel").c_str(),
95  iConfig.getParameter<std::string>("xDetLabel").c_str());
96  sprintf(htitle,"atan (%4.2f*%s / %s multiplicity ratio)",
97  m_scfact,
98  iConfig.getParameter<std::string>("yDetLabel").c_str(),
99  iConfig.getParameter<std::string>("xDetLabel").c_str()
100  );
101  m_atanyoverxvsbxrun2D = m_fhm.makeTH2F(hname,htitle,3564,-0.5,3564-0.5,
102  iConfig.getParameter<unsigned int>("rBins"),0.,1.6);
103  }
104  if(m_runHistoBXProfile) {
105  sprintf(hname,"%sOver%svsbx",
106  iConfig.getParameter<std::string>("yDetLabel").c_str(),
107  iConfig.getParameter<std::string>("xDetLabel").c_str());
108  sprintf(htitle,"atan (%4.2f*%s / %s multiplicity ratio)",
109  m_scfact,
110  iConfig.getParameter<std::string>("yDetLabel").c_str(),
111  iConfig.getParameter<std::string>("xDetLabel").c_str()
112  );
113  m_atanyoverxvsbxrun = m_fhm.makeTProfile(hname,htitle,3564,-0.5,3564-0.5);
114  }
115  }
116 
117 }
118 
119 
120 MultiplicityCorrelatorHistogramMaker::~MultiplicityCorrelatorHistogramMaker() { }
121 
122 void MultiplicityCorrelatorHistogramMaker::beginRun(const edm::Run& iRun) {
123 
124  m_rhm.beginRun(iRun);
125  m_fhm.beginRun(iRun);
126 
127 }
128 
129 void MultiplicityCorrelatorHistogramMaker::fill(const edm::Event& iEvent, const int xmult, const int ymult) {
130 
131 
132 
133  const int bx = iEvent.bunchCrossing();
134 
135  if(m_yvsxmult) m_yvsxmult->Fill(xmult,ymult);
136  if(m_atanyoverx) m_atanyoverx->Fill(atan2(ymult*m_scfact,xmult));
137 
138  if(m_yvsxmultrun && *m_yvsxmultrun) (*m_yvsxmultrun)->Fill(xmult,ymult);
139  if(m_atanyoverxrun && *m_atanyoverxrun) (*m_atanyoverxrun)->Fill(atan2(ymult*m_scfact,xmult));
140  if(m_atanyoverxvsbxrun && *m_atanyoverxvsbxrun) (*m_atanyoverxvsbxrun)->Fill(bx%3564,atan2(ymult*m_scfact,xmult));
141  if(m_atanyoverxvsbxrun2D && *m_atanyoverxvsbxrun2D) (*m_atanyoverxvsbxrun2D)->Fill(bx%3564,atan2(ymult*m_scfact,xmult));
142 
143  if(m_yvsxmultprofvsbxrun && *m_yvsxmultprofvsbxrun) (*m_yvsxmultprofvsbxrun)->Fill(bx%3564,xmult,ymult);
144  if(m_xvsymultprofvsbxrun && *m_xvsymultprofvsbxrun) (*m_xvsymultprofvsbxrun)->Fill(bx%3564,ymult,xmult);
145 }
146 
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
MultiplicityCorrelatorHistogramMaker::fill
void fill(const edm::Event &iEvent, const int xmult, const int ymult)
Definition: MultiplicityCorrelatorHistogramMaker.cc:129
nullptr
#define nullptr
Definition: GCC11Compatibility.h:37
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:64
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
RunHistogramManager::makeTH1F
TH1F ** makeTH1F(const char *name, const char *title, const unsigned int nbinx, const double xmin, const double xmax)
Definition: RunHistogramManager.cc:35
iEvent
int iEvent
Definition: GenABIO.cc:224
RunHistogramManager::makeTProfile
TProfile ** makeTProfile(const char *name, const char *title, const unsigned int nbinx, const double xmin, const double xmax)
Definition: RunHistogramManager.cc:59
electrons_cff.bool
bool
Definition: electrons_cff.py:359
RunHistogramManager::makeTProfile2D
TProfile2D ** makeTProfile2D(const char *name, const char *title, const unsigned int nbinx, const double xmin, const double xmax, const unsigned int nbiny, const double ymin, const double ymax)
Definition: RunHistogramManager.cc:84
RunHistogramManager::beginRun
void beginRun(const edm::Run &iRun)
Definition: RunHistogramManager.cc:96
RunHistogramManager::makeTH2F
TH2F ** makeTH2F(const char *name, const char *title, const unsigned int nbinx, const double xmin, const double xmax, const unsigned int nbiny, const double ymin, const double ymax)
Definition: RunHistogramManager.cc:72
edm::Run
Definition: Run.h:45
Generated for CMSSW Reference Manual by doxygen 1.8.11