12 if (strRunType ==
"online") {
14 }
else if (strRunType ==
"offline") {
16 }
else if (strRunType ==
"relval") {
18 }
else if (strRunType ==
"allplots") {
41 const auto&& superchambers =
station->superChambers();
47 const auto&
chambers = superchambers.front()->chambers();
53 return chambers.front()->nEtaPartitions();
63 for (
auto sch :
station->superChambers()) {
64 for (
auto pchamber : sch->chambers()) {
67 for (
auto iEta : pchamber->etaPartitions()) {
78 const int region_number =
region->region();
81 const auto&& superchambers =
station->superChambers();
83 const int station_number =
station->station();
84 const int num_superchambers = (station_number == 1 ? 36 : 18);
88 const int num_vfat = num_etas * max_vfat;
89 const int strip1st = (station_number == 2 ? 1 : 0);
94 Int_t nMinIdxChamber = 1048576;
95 Int_t nMaxIdxChamber = -1048576;
96 for (
auto sch : superchambers) {
97 auto nIdxChamber = sch->chambers().front()->id().chamber();
98 if (nMinIdxChamber > nIdxChamber)
99 nMinIdxChamber = nIdxChamber;
100 if (nMaxIdxChamber < nIdxChamber)
101 nMaxIdxChamber = nIdxChamber;
104 const auto&
chambers = superchambers.front()->chambers();
107 int layer_number = pchamber->id().layer();
108 ME3IdsKey key3(region_number, station_number, layer_number);
131 Int_t re1 = std::get<0>(key1), st1 = std::get<1>(key1), la1 = std::get<2>(key1);
132 Int_t re2 = std::get<0>(key2), st2 = std::get<1>(key2), la2 = std::get<2>(key2);
133 if (re1 < 0 && re2 > 0)
135 if (re1 > 0 && re2 < 0)
137 Bool_t bRes = (re1 < 0);
138 Int_t sum1 = 256 *
std::abs(re1) + 16 * st1 + 1 * la1;
139 Int_t sum2 = 256 *
std::abs(re2) + 16 * st2 + 1 * la2;
149 std::vector<ME3IdsKey> listLayers;
152 listLayers.push_back(
key);
154 for (Int_t
i = 0;
i < (Int_t)listLayers.size();
i++)
160 h2Res->setYTitle(
"Layer");
162 if (h2Res ==
nullptr)
166 h2Res->setBinLabel(
i, Form(
"%i",
i), 1);
167 for (Int_t
i = 1;
i <= (Int_t)listLayers.size();
i++) {
168 auto key = listLayers[
i - 1];
173 h2Res->setBinLabel(
i,
label, 2);
175 h2Res->setBinContent(0,
i, nNumCh);
189 ME2IdsKey key2{gid.region(), gid.station()};
190 ME3IdsKey key3{gid.region(), gid.station(), gid.layer()};
191 ME4IdsKey key3WithChamber{gid.region(), gid.station(), gid.layer(), gid.chamber()};
207 Int_t nCh = gid.chamber();
208 Int_t nLa = gid.layer();
209 char cLS = (nCh % 2 == 0 ?
'L' :
'S');
218 ME4IdsKey key4{gid.region(), gid.station(), gid.layer(), eId.
ieta()};
219 ME3IdsKey key2WithEta{gid.region(), gid.station(), eId.
ieta()};
236 auto strSuffixName = Form(
"_GE%d1-E%02i", gid.station(), eId.
ieta());
237 auto strSuffixTitle = Form(
" GE%d1-E%02i", gid.station(), eId.
ieta());
248 auto listSuperChambers =
station->superChambers();
250 Bool_t bDoneEven =
false, bDoneOdd =
false;
253 for (
auto superchamber : listSuperChambers) {
254 Int_t chamberNo = superchamber->id().chamber();
255 if (chamberNo % 2 == 0 && bDoneEven)
257 if (chamberNo % 2 != 0 && bDoneOdd)
260 auto& etaPartitions = superchamber->chambers().front()->etaPartitions();
268 std::vector<Float_t> listRadiusLower, listRadiusUpper;
269 for (
auto iEta : etaPartitions) {
272 LocalPoint lp1(0.0, -0.5 * fHeight), lp2(0.0, 0.5 * fHeight);
273 auto& surface =
iEta->surface();
274 GlobalPoint gp1 = surface.toGlobal(lp1), gp2 = surface.toGlobal(lp2);
275 Float_t fR1 = gp1.
perp(), fR2 = gp2.perp();
277 listRadiusLower.push_back(fRL);
278 listRadiusUpper.push_back(fRH);
284 std::sort(listRadiusLower.begin(), listRadiusLower.end());
285 std::sort(listRadiusUpper.begin(), listRadiusUpper.end());
287 std::vector<Float_t>& listR =
290 listR.push_back(listRadiusLower.front());
291 for (
int i = 1;
i < (
int)listRadiusLower.size();
i++) {
292 listR.push_back(0.5 * (listRadiusLower[
i] + listRadiusUpper[
i - 1]));
294 listR.push_back(listRadiusUpper.back());
296 if (chamberNo % 2 == 0)
298 if (chamberNo % 2 != 0)
301 if (bDoneEven && bDoneOdd)
309 auto listSuperChambers =
station->superChambers();
312 std::vector<std::pair<Int_t, std::pair<std::pair<Float_t, Float_t>, Bool_t>>> listDivPhi;
315 for (
auto superchamber : listSuperChambers) {
316 auto iEta = superchamber->chambers().front()->etaPartitions().front();
319 Float_t fWidthStrip =
std::abs(
iEta->centreOfStrip((Int_t)1).x() -
iEta->centreOfStrip((Int_t)0).x());
321 Int_t nStripMid = (Int_t)
iEta->strip(lpRef);
322 Int_t nFirstStrip = 1 - ((nNumStripEta / 2) - nStripMid);
323 Int_t nLastStrip = nFirstStrip + nNumStripEta - 1;
325 auto& surface =
iEta->surface();
326 LocalPoint lpF =
iEta->centreOfStrip((Float_t)(nFirstStrip - 0.5));
327 LocalPoint lpL =
iEta->centreOfStrip((Float_t)(nLastStrip + 0.5));
331 Float_t fPhiF = gpF.
phi();
332 Float_t fPhiL = gpL.
phi();
333 if (fPhiF * fPhiL < 0 &&
std::abs(fPhiF) > 0.5 * 3.14159265359) {
335 fPhiF += 2 * 3.14159265359;
337 fPhiL += 2 * 3.14159265359;
339 Bool_t bFlipped = fPhiF > fPhiL;
340 Float_t fPhiMin =
std::min(fPhiF, fPhiL);
341 Float_t fPhiMax =
std::max(fPhiF, fPhiL);
343 listDivPhi.emplace_back();
344 listDivPhi.back().first =
iEta->id().chamber();
345 listDivPhi.back().second.first.first = fPhiMin;
346 listDivPhi.back().second.first.second = fPhiMax;
347 listDivPhi.back().second.second = bFlipped;
351 for (
auto p : listDivPhi) {
353 stationInfo.
fMinPhi_ =
p.second.first.first;
float stripLength() const override
bool checkRefs(const std::vector< T *> &)
dqm::impl::MonitorElement * CreateSummaryHist(DQMStore::IBooker &ibooker, TString strName)
virtual int ProcessWithMEMap4(BookingHelper &bh, ME4IdsKey key)
TString getSuffixName(Int_t region_id)
std::vector< GEMDetId > listChamberId_
virtual int ProcessWithMEMap2AbsReWithEta(BookingHelper &bh, ME3IdsKey key)
Geom::Phi< T > phi() const
int initGeometry(edm::EventSetup const &iSetup)
std::tuple< Int_t, Int_t, Int_t, Int_t > ME4IdsKey
int readGeometryRadiusInfoChamber(const GEMStation *station, MEStationInfo &stationInfo)
int getNumModule(const int)
std::map< ME3IdsKey, MEStationInfo > mapStationInfo_
int keyToRegion(ME2IdsKey key)
int readGeometryPhiInfoChamber(const GEMStation *station, MEStationInfo &stationInfo)
std::map< ME3IdsKey, int > mapStationToIdx_
static const int maxChan_
int keyToLayer(ME3IdsKey key)
Log< level::Error, false > LogError
TString getSuffixTitle(Int_t region_id)
std::vector< Float_t > listRadiusOddChamber_
std::string log_category_
std::map< ME3IdsKey, bool > MEMap2AbsReWithEtaCheck_
std::map< ME2IdsKey, bool > MEMap2Check_
virtual int ProcessWithMEMap3(BookingHelper &bh, ME3IdsKey key)
virtual void setXTitle(std::string const &title)
int GenerateMEPerChamber(DQMStore::IBooker &ibooker)
std::vector< Float_t > listRadiusEvenChamber_
std::map< ME3IdsKey, bool > MEMap3Check_
int getNumEtaPartitions(const GEMStation *)
virtual int ProcessWithMEMap3WithChamber(BookingHelper &bh, ME4IdsKey key)
edm::ESGetToken< GEMGeometry, MuonGeometryRecord > geomToken_
Abs< T >::type abs(const T &t)
int getMaxVFAT(const int)
int keyToStation(ME2IdsKey key)
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
std::map< ME4IdsKey, bool > MEMap4Check_
constexpr int ieta() const
const GEMGeometry * GEMGeometry_
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
virtual int ProcessWithMEMap2WithEta(BookingHelper &bh, ME3IdsKey key)
int SortingLayers(std::vector< ME3IdsKey > &listLayers)
GEMDQMBase(const edm::ParameterSet &cfg)
std::map< ME3IdsKey, bool > MEMap2WithEtaCheck_
std::map< GEMDetId, std::vector< const GEMEtaPartition * > > mapEtaPartition_
const std::vector< const GEMRegion * > & regions() const
Return a vector of all GEM regions.
std::tuple< Int_t, Int_t > ME2IdsKey
std::tuple< Int_t, Int_t, Int_t > ME3IdsKey
static char chambers[264][20]
std::map< ME4IdsKey, bool > MEMap3WithChCheck_
virtual int ProcessWithMEMap2(BookingHelper &bh, ME2IdsKey key)