|
|
#include <RPCDigiValid.h>
Definition at line 17 of file RPCDigiValid.h.
◆ RPCDigiValid()
◆ ~RPCDigiValid()
RPCDigiValid::~RPCDigiValid |
( |
| ) |
|
|
override |
◆ analyze()
Reimplemented from DQMEDAnalyzer.
Definition at line 34 of file RPCDigiValid.cc.
45 PSimHitContainer::const_iterator simIt;
48 std::map<RPCDetId, std::vector<double>> allsims;
50 for (simIt =
simHit->begin(); simIt !=
simHit->end(); simIt++) {
52 const RPCRoll *soll = dynamic_cast<const RPCRoll *>(rpcGeom->
roll(Rsid));
53 int ptype = simIt->particleType();
56 std::vector<double> buff;
57 if (allsims.find(Rsid) != allsims.end()) {
61 buff.push_back(simIt->localPosition().x());
72 if (Rsid.
region() == (+1)) {
76 }
else if (Rsid.
region() == (-1)) {
85 for (detUnitIt = rpcDigis->begin(); detUnitIt != rpcDigis->end(); ++detUnitIt) {
86 const RPCDetId Rsid = (*detUnitIt).first;
87 const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->
roll(Rsid));
90 std::vector<double> sims;
91 if (allsims.find(Rsid) != allsims.end()) {
100 if (Rsid.
region() == (+1)) {
103 }
else if (Rsid.
region() == (-1)) {
109 const double digiTime = digiIt->hasTime() ? digiIt->time() : digiIt->bx() * 25;
111 if (digiIt->hasTime()) {
120 if (sims.size() == 1 && ndigi == 1) {
125 if (Rsid.
ring() == -2)
127 else if (Rsid.
ring() == -1)
129 else if (Rsid.
ring() == 0)
131 else if (Rsid.
ring() == 1)
133 else if (Rsid.
ring() == 2)
151 if (Rsid.
ring() == 2) {
153 if (Rsid.
roll() == 1)
155 if (Rsid.
roll() == 2)
157 if (Rsid.
roll() == 3)
161 if (Rsid.
roll() == 1)
163 if (Rsid.
roll() == 2)
165 if (Rsid.
roll() == 3)
169 if (Rsid.
ring() == 3) {
170 if (Rsid.
roll() == 1)
172 if (Rsid.
roll() == 2)
174 if (Rsid.
roll() == 3)
179 if (Rsid.
region() == (+1)) {
189 if (Rsid.
region() == (-1)) {
References funct::abs(), RPCRoll::centreOfStrip(), edm::EventSetup::get(), get, RPCRoll::isIRPC(), RPCDetId::layer(), AlCaHLTBitMon_ParallelJobs::p, OfflineHarvestingSequence_cosmic::ptype, FastTimerService_cff::range, RPCDetId::region(), RPCDetId::ring(), RPCGeometry::roll(), RPCDetId::roll(), rpcPointValidation_cfi::simHit, SiPixelPhase1TrackingParticleV_cfi::simHitToken, RPCDetId::station(), and GeomDet::toGlobal().
◆ bookHistograms()
Implements DQMEDAnalyzer.
Definition at line 203 of file RPCDigiValid.cc.
204 booker.setCurrentFolder(
"RPCDigisV/RPCDigis");
206 xyview = booker.book2D(
"X_Vs_Y_View",
"X_Vs_Y_View", 155, -775., 775., 155, -775., 775.);
208 xyvDplu4 = booker.book2D(
"Dplu4_XvsY",
"Dplu4_XvsY", 155, -775., 775., 155, -775., 775.);
209 xyvDmin4 = booker.book2D(
"Dmin4_XvsY",
"Dmin4_XvsY", 155, -775., 775., 155, -775., 775.);
211 rzview = booker.book2D(
"R_Vs_Z_View",
"R_Vs_Z_View", 216, -1080., 1080., 52, 260., 780.);
212 Res = booker.book1D(
"Digi_SimHit_difference",
"Digi_SimHit_difference", 300, -8, 8);
213 ResWmin2 = booker.book1D(
"W_Min2_Residuals",
"W_Min2_Residuals", 400, -8, 8);
214 ResWmin1 = booker.book1D(
"W_Min1_Residuals",
"W_Min1_Residuals", 400, -8, 8);
215 ResWzer0 = booker.book1D(
"W_Zer0_Residuals",
"W_Zer0_Residuals", 400, -8, 8);
216 ResWplu1 = booker.book1D(
"W_Plu1_Residuals",
"W_Plu1_Residuals", 400, -8, 8);
217 ResWplu2 = booker.book1D(
"W_Plu2_Residuals",
"W_Plu2_Residuals", 400, -8, 8);
219 ResLayer1_barrel = booker.book1D(
"ResLayer1_barrel",
"ResLayer1_barrel", 400, -8, 8);
220 ResLayer2_barrel = booker.book1D(
"ResLayer2_barrel",
"ResLayer2_barrel", 400, -8, 8);
221 ResLayer3_barrel = booker.book1D(
"ResLayer3_barrel",
"ResLayer3_barrel", 400, -8, 8);
222 ResLayer4_barrel = booker.book1D(
"ResLayer4_barrel",
"ResLayer4_barrel", 400, -8, 8);
223 ResLayer5_barrel = booker.book1D(
"ResLayer5_barrel",
"ResLayer5_barrel", 400, -8, 8);
224 ResLayer6_barrel = booker.book1D(
"ResLayer6_barrel",
"ResLayer6_barrel", 400, -8, 8);
226 BxDist = booker.book1D(
"Bunch_Crossing",
"Bunch_Crossing", 20, -10., 10.);
227 StripProf = booker.book1D(
"Strip_Profile",
"Strip_Profile", 100, 0, 100);
229 BxDisc_4Plus = booker.book1D(
"BxDisc_4Plus",
"BxDisc_4Plus", 20, -10., 10.);
230 BxDisc_4Min = booker.book1D(
"BxDisc_4Min",
"BxDisc_4Min", 20, -10., 10.);
233 ResDmin1 = booker.book1D(
"Disk_Min1_Residuals",
"Disk_Min1_Residuals", 400, -8, 8);
234 ResDmin2 = booker.book1D(
"Disk_Min2_Residuals",
"Disk_Min2_Residuals", 400, -8, 8);
235 ResDmin3 = booker.book1D(
"Disk_Min3_Residuals",
"Disk_Min3_Residuals", 400, -8, 8);
236 ResDplu1 = booker.book1D(
"Disk_Plu1_Residuals",
"Disk_Plu1_Residuals", 400, -8, 8);
237 ResDplu2 = booker.book1D(
"Disk_Plu2_Residuals",
"Disk_Plu2_Residuals", 400, -8, 8);
238 ResDplu3 = booker.book1D(
"Disk_Plu3_Residuals",
"Disk_Plu3_Residuals", 400, -8, 8);
240 ResDmin4 = booker.book1D(
"Disk_Min4_Residuals",
"Disk_Min4_Residuals", 400, -8, 8);
241 ResDplu4 = booker.book1D(
"Disk_Plu4_Residuals",
"Disk_Plu4_Residuals", 400, -8, 8);
243 Res_Endcap1_Ring2_A = booker.book1D(
"Res_Endcap1_Ring2_A",
"Res_Endcap1_Ring2_A", 400, -8, 8);
244 Res_Endcap1_Ring2_B = booker.book1D(
"Res_Endcap1_Ring2_B",
"Res_Endcap1_Ring2_B", 400, -8, 8);
245 Res_Endcap1_Ring2_C = booker.book1D(
"Res_Endcap1_Ring2_C",
"Res_Endcap1_Ring2_C", 400, -8, 8);
247 Res_Endcap23_Ring2_A = booker.book1D(
"Res_Endcap23_Ring2_A",
"Res_Endcap23_Ring2_A", 400, -8, 8);
248 Res_Endcap23_Ring2_B = booker.book1D(
"Res_Endcap23_Ring2_B",
"Res_Endcap23_Ring2_B", 400, -8, 8);
249 Res_Endcap23_Ring2_C = booker.book1D(
"Res_Endcap23_Ring2_C",
"Res_Endcap23_Ring2_C", 400, -8, 8);
257 booker.book1D(
"DigiTimeAll",
"Digi time including present electronics;Digi time (ns)", 100, -12.5, 12.5);
258 hDigiTime = booker.book1D(
"DigiTime",
"Digi time only with timing information;Digi time (ns)", 100, -12.5, 12.5);
259 hDigiTimeIRPC = booker.book1D(
"DigiTimeIRPC",
"IRPC Digi time;Digi time (ns)", 100, -12.5, 12.5);
260 hDigiTimeNoIRPC = booker.book1D(
"DigiTimeNoIRPC",
"non-IRPC Digi time;Digi time (ns)", 100, -12.5, 12.5);
References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), and dqm::implementation::NavigatorBase::setCurrentFolder().
◆ BxDisc_4Min
◆ BxDisc_4Plus
◆ BxDist
◆ digiLabel
std::string RPCDigiValid::digiLabel |
|
private |
◆ hDigiTime
◆ hDigiTimeAll
◆ hDigiTimeIRPC
◆ hDigiTimeNoIRPC
◆ outputFile_
std::string RPCDigiValid::outputFile_ |
|
private |
◆ Res
◆ Res_Endcap123_Ring3_A
◆ Res_Endcap123_Ring3_B
◆ Res_Endcap123_Ring3_C
◆ Res_Endcap1_Ring2_A
◆ Res_Endcap1_Ring2_B
◆ Res_Endcap1_Ring2_C
◆ Res_Endcap23_Ring2_A
◆ Res_Endcap23_Ring2_B
◆ Res_Endcap23_Ring2_C
◆ ResDmin1
◆ ResDmin2
◆ ResDmin3
◆ ResDmin4
◆ ResDplu1
◆ ResDplu2
◆ ResDplu3
◆ ResDplu4
◆ ResLayer1_barrel
◆ ResLayer2_barrel
◆ ResLayer3_barrel
◆ ResLayer4_barrel
◆ ResLayer5_barrel
◆ ResLayer6_barrel
◆ ResWmin1
◆ ResWmin2
◆ ResWplu1
◆ ResWplu2
◆ ResWzer0
◆ rpcDigiToken
◆ rzview
◆ simHitToken
◆ StripProf
◆ xyvDmin4
◆ xyvDplu4
◆ xyview
MonitorElement * Res_Endcap123_Ring3_B
MonitorElement * ResWmin2
MonitorElement * ResDmin1
MonitorElement * ResLayer4_barrel
MonitorElement * Res_Endcap123_Ring3_C
int region() const
Region id: 0 for Barrel, +/-1 For +/- Endcap.
MonitorElement * Res_Endcap23_Ring2_B
MonitorElement * ResLayer1_barrel
MonitorElement * Res_Endcap1_Ring2_C
MonitorElement * ResLayer5_barrel
MonitorElement * Res_Endcap23_Ring2_C
const RPCRoll * roll(RPCDetId id) const
Return a roll given its id.
MonitorElement * ResDplu4
MonitorElement * ResDplu3
T getUntrackedParameter(std::string const &, T const &) const
MonitorElement * ResWplu1
std::vector< RPCDigi >::const_iterator const_iterator
MonitorElement * hDigiTimeNoIRPC
MonitorElement * Res_Endcap1_Ring2_A
MonitorElement * ResDplu1
MonitorElement * xyvDplu4
MonitorElement * hDigiTime
LocalPoint centreOfStrip(int strip) const
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
MonitorElement * ResDmin2
MonitorElement * BxDisc_4Plus
MonitorElement * ResWplu2
edm::EDGetTokenT< edm::PSimHitContainer > simHitToken
MonitorElement * BxDisc_4Min
MonitorElement * StripProf
MonitorElement * ResLayer3_barrel
MonitorElement * xyvDmin4
MonitorElement * ResDplu2
MonitorElement * ResWzer0
std::pair< const_iterator, const_iterator > Range
MonitorElement * hDigiTimeIRPC
MonitorElement * Res_Endcap1_Ring2_B
MonitorElement * ResDmin3
Abs< T >::type abs(const T &t)
MonitorElement * ResWmin1
MonitorElement * ResLayer2_barrel
MonitorElement * hDigiTimeAll
edm::EDGetTokenT< RPCDigiCollection > rpcDigiToken
MonitorElement * Res_Endcap123_Ring3_A
MonitorElement * Res_Endcap23_Ring2_A
MonitorElement * ResDmin4
MonitorElement * ResLayer6_barrel