#include <RPCDigiValid.h>

Inheritance diagram for RPCDigiValid:

Public Member Functions
	RPCDigiValid (const edm::ParameterSet &ps)

	~RPCDigiValid () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Private Attributes
MonitorElement *	BxDisc_4Min

MonitorElement *	BxDisc_4Plus

MonitorElement *	BxDist

std::string	digiLabel

MonitorElement *	hDigiTime

MonitorElement *	hDigiTimeAll

MonitorElement *	hDigiTimeIRPC

MonitorElement *	hDigiTimeNoIRPC

std::string	outputFile_

MonitorElement *	Res

MonitorElement *	Res_Endcap123_Ring3_A

MonitorElement *	Res_Endcap123_Ring3_B

MonitorElement *	Res_Endcap123_Ring3_C

MonitorElement *	Res_Endcap1_Ring2_A

MonitorElement *	Res_Endcap1_Ring2_B

MonitorElement *	Res_Endcap1_Ring2_C

MonitorElement *	Res_Endcap23_Ring2_A

MonitorElement *	Res_Endcap23_Ring2_B

MonitorElement *	Res_Endcap23_Ring2_C

MonitorElement *	ResDmin1

MonitorElement *	ResDmin2

MonitorElement *	ResDmin3

MonitorElement *	ResDmin4

MonitorElement *	ResDplu1

MonitorElement *	ResDplu2

MonitorElement *	ResDplu3

MonitorElement *	ResDplu4

MonitorElement *	ResLayer1_barrel

MonitorElement *	ResLayer2_barrel

MonitorElement *	ResLayer3_barrel

MonitorElement *	ResLayer4_barrel

MonitorElement *	ResLayer5_barrel

MonitorElement *	ResLayer6_barrel

MonitorElement *	ResWmin1

MonitorElement *	ResWmin2

MonitorElement *	ResWplu1

MonitorElement *	ResWplu2

MonitorElement *	ResWzer0

edm::EDGetTokenT< RPCDigiCollection >	rpcDigiToken

MonitorElement *	rzview

edm::EDGetTokenT< edm::PSimHitContainer >	simHitToken

MonitorElement *	StripProf

MonitorElement *	xyvDmin4

MonitorElement *	xyvDplu4

MonitorElement *	xyview

Detailed Description

Definition at line 17 of file RPCDigiValid.h.

Constructor & Destructor Documentation

RPCDigiValid::RPCDigiValid ( const edm::ParameterSet & ps )

Definition at line 18 of file RPCDigiValid.cc.

References edm::ParameterSet::getUntrackedParameter(), and SiPixelPhase1TrackingParticleV_cfi::simHitToken.

                                                  {
   //  Init the tokens for run data retrieval - stanislav
   //  ps.getUntackedParameter<InputTag> retrieves a InputTag from the
   //  configuration. The second param is default value module, instance and
   //  process labels may be passed in a single string if separated by colon ':'
   //  (@see the edm::InputTag constructor documentation)
   simHitToken = consumes<PSimHitContainer>(
       ps.getUntrackedParameter<edm::InputTag>("simHitTag", edm::InputTag("g4SimHits:MuonRPCHits")));
   rpcDigiToken = consumes<RPCDigiCollection>(
       ps.getUntrackedParameter<edm::InputTag>("rpcDigiTag", edm::InputTag("simMuonRPCDigis")));
 
   outputFile_ = ps.getUntrackedParameter<string>("outputFile", "rpcDigiValidPlots.root");
 }

RPCDigiValid::~RPCDigiValid ( )

override

Definition at line 32 of file RPCDigiValid.cc.

32 {}

Member Function Documentation

void RPCDigiValid::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotected

Definition at line 34 of file RPCDigiValid.cc.

References funct::abs(), RPCRoll::centreOfStrip(), edm::EventSetup::get(), RPCRoll::isIRPC(), RPCDetId::layer(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), RPCDetId::region(), RPCDetId::ring(), RPCGeometry::roll(), RPCDetId::roll(), rpcPointValidation_cfi::simHit, SiPixelPhase1TrackingParticleV_cfi::simHitToken, RPCDetId::station(), GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                            {
   // Get the RPC Geometry
   edm::ESHandle<RPCGeometry> rpcGeom;
   eventSetup.get<MuonGeometryRecord>().get(rpcGeom);
 
   edm::Handle<PSimHitContainer> simHit;
   edm::Handle<RPCDigiCollection> rpcDigis;
   event.getByToken(simHitToken, simHit);
   event.getByToken(rpcDigiToken, rpcDigis);
 
   // Loop on simhits
   PSimHitContainer::const_iterator simIt;
 
   // loop over Simhit
   std::map<RPCDetId, std::vector<double>> allsims;
 
   for (simIt = simHit->begin(); simIt != simHit->end(); simIt++) {
     RPCDetId Rsid = (RPCDetId)(*simIt).detUnitId();
     const RPCRoll *soll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(Rsid));
     int ptype = simIt->particleType();
 
     if (ptype == 13 || ptype == -13) {
       std::vector<double> buff;
       if (allsims.find(Rsid) != allsims.end()) {
         buff = allsims[Rsid];
       }
 
       buff.push_back(simIt->localPosition().x());
 
       allsims[Rsid] = buff;
     }
     GlobalPoint p = soll->toGlobal(simIt->localPosition());
 
     double sim_x = p.x();
     double sim_y = p.y();
 
     xyview->Fill(sim_x, sim_y);
 
     if (Rsid.region() == (+1)) {
       if (Rsid.station() == 4) {
         xyvDplu4->Fill(sim_x, sim_y);
       }
     } else if (Rsid.region() == (-1)) {
       if (Rsid.station() == 4) {
         xyvDmin4->Fill(sim_x, sim_y);
       }
     }
     rzview->Fill(p.z(), p.perp());
   }
   // loop over Digis
   RPCDigiCollection::DigiRangeIterator detUnitIt;
   for (detUnitIt = rpcDigis->begin(); detUnitIt != rpcDigis->end(); ++detUnitIt) {
     const RPCDetId Rsid = (*detUnitIt).first;
     const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(Rsid));
 
     const RPCDigiCollection::Range &range = (*detUnitIt).second;
     std::vector<double> sims;
     if (allsims.find(Rsid) != allsims.end()) {
       sims = allsims[Rsid];
     }
 
     int ndigi = 0;
     for (RPCDigiCollection::const_iterator digiIt = range.first; digiIt != range.second; ++digiIt) {
       StripProf->Fill(digiIt->strip());
       BxDist->Fill(digiIt->bx());
       // bx for 4 endcaps
       if (Rsid.region() == (+1)) {
         if (Rsid.station() == 4)
           BxDisc_4Plus->Fill(digiIt->bx());
       } else if (Rsid.region() == (-1)) {
         if (Rsid.station() == 4)
           BxDisc_4Min->Fill(digiIt->bx());
       }
 
       // Fill timing information
       const double digiTime = digiIt->hasTime() ? digiIt->time() : digiIt->bx() * 25;
       hDigiTimeAll->Fill(digiTime);
       if (digiIt->hasTime()) {
         hDigiTime->Fill(digiTime);
         if (roll->isIRPC())
           hDigiTimeIRPC->Fill(digiTime);
         else
           hDigiTimeNoIRPC->Fill(digiTime);
       }
     }
 
     if (sims.size() == 1 && ndigi == 1) {
       double dis = roll->centreOfStrip(range.first->strip()).x() - sims[0];
       Res->Fill(dis);
 
       if (Rsid.region() == 0) {
         if (Rsid.ring() == -2)
           ResWmin2->Fill(dis);
         else if (Rsid.ring() == -1)
           ResWmin1->Fill(dis);
         else if (Rsid.ring() == 0)
           ResWzer0->Fill(dis);
         else if (Rsid.ring() == 1)
           ResWplu1->Fill(dis);
         else if (Rsid.ring() == 2)
           ResWplu2->Fill(dis);
         // barrel layers
         if (Rsid.station() == 1 && Rsid.layer() == 1)
           ResLayer1_barrel->Fill(dis);
         if (Rsid.station() == 1 && Rsid.layer() == 2)
           ResLayer2_barrel->Fill(dis);
         if (Rsid.station() == 2 && Rsid.layer() == 1)
           ResLayer3_barrel->Fill(dis);
         if (Rsid.station() == 2 && Rsid.layer() == 2)
           ResLayer4_barrel->Fill(dis);
         if (Rsid.station() == 3)
           ResLayer5_barrel->Fill(dis);
         if (Rsid.station() == 4)
           ResLayer6_barrel->Fill(dis);
       }
       // endcap layers residuals
       if (Rsid.region() != 0) {
         if (Rsid.ring() == 2) {
           if (abs(Rsid.station()) == 1) {
             if (Rsid.roll() == 1)
               Res_Endcap1_Ring2_A->Fill(dis);
             if (Rsid.roll() == 2)
               Res_Endcap1_Ring2_B->Fill(dis);
             if (Rsid.roll() == 3)
               Res_Endcap1_Ring2_C->Fill(dis);
           }
           if (abs(Rsid.station()) == 2 || abs(Rsid.station()) == 3) {
             if (Rsid.roll() == 1)
               Res_Endcap23_Ring2_A->Fill(dis);
             if (Rsid.roll() == 2)
               Res_Endcap23_Ring2_B->Fill(dis);
             if (Rsid.roll() == 3)
               Res_Endcap23_Ring2_C->Fill(dis);
           }
         }
         if (Rsid.ring() == 3) {
           if (Rsid.roll() == 1)
             Res_Endcap123_Ring3_A->Fill(dis);
           if (Rsid.roll() == 2)
             Res_Endcap123_Ring3_B->Fill(dis);
           if (Rsid.roll() == 3)
             Res_Endcap123_Ring3_C->Fill(dis);
         }
       }
 
       if (Rsid.region() == (+1)) {
         if (Rsid.station() == 1)
           ResDplu1->Fill(dis);
         else if (Rsid.station() == 2)
           ResDplu2->Fill(dis);
         else if (Rsid.station() == 3)
           ResDplu3->Fill(dis);
         else if (Rsid.station() == 4)
           ResDplu4->Fill(dis);
       }
       if (Rsid.region() == (-1)) {
         if (Rsid.station() == 1)
           ResDmin1->Fill(dis);
         else if (Rsid.station() == 2)
           ResDmin2->Fill(dis);
         else if (Rsid.station() == 3)
           ResDmin3->Fill(dis);
         else if (Rsid.station() == 4)
           ResDmin4->Fill(dis);
       }
     }
   }
 }

void RPCDigiValid::bookHistograms	(	DQMStore::IBooker &	booker,
		edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

overrideprotected

Definition at line 203 of file RPCDigiValid.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), and DQMStore::IBooker::setCurrentFolder().

                                                                                                          {
   booker.setCurrentFolder("RPCDigisV/RPCDigis");
 
   xyview = booker.book2D("X_Vs_Y_View", "X_Vs_Y_View", 155, -775., 775., 155, -775., 775.);
 
   xyvDplu4 = booker.book2D("Dplu4_XvsY", "Dplu4_XvsY", 155, -775., 775., 155, -775., 775.);
   xyvDmin4 = booker.book2D("Dmin4_XvsY", "Dmin4_XvsY", 155, -775., 775., 155, -775., 775.);
 
   rzview = booker.book2D("R_Vs_Z_View", "R_Vs_Z_View", 216, -1080., 1080., 52, 260., 780.);
   Res = booker.book1D("Digi_SimHit_difference", "Digi_SimHit_difference", 300, -8, 8);
   ResWmin2 = booker.book1D("W_Min2_Residuals", "W_Min2_Residuals", 400, -8, 8);
   ResWmin1 = booker.book1D("W_Min1_Residuals", "W_Min1_Residuals", 400, -8, 8);
   ResWzer0 = booker.book1D("W_Zer0_Residuals", "W_Zer0_Residuals", 400, -8, 8);
   ResWplu1 = booker.book1D("W_Plu1_Residuals", "W_Plu1_Residuals", 400, -8, 8);
   ResWplu2 = booker.book1D("W_Plu2_Residuals", "W_Plu2_Residuals", 400, -8, 8);
 
   ResLayer1_barrel = booker.book1D("ResLayer1_barrel", "ResLayer1_barrel", 400, -8, 8);
   ResLayer2_barrel = booker.book1D("ResLayer2_barrel", "ResLayer2_barrel", 400, -8, 8);
   ResLayer3_barrel = booker.book1D("ResLayer3_barrel", "ResLayer3_barrel", 400, -8, 8);
   ResLayer4_barrel = booker.book1D("ResLayer4_barrel", "ResLayer4_barrel", 400, -8, 8);
   ResLayer5_barrel = booker.book1D("ResLayer5_barrel", "ResLayer5_barrel", 400, -8, 8);
   ResLayer6_barrel = booker.book1D("ResLayer6_barrel", "ResLayer6_barrel", 400, -8, 8);
 
   BxDist = booker.book1D("Bunch_Crossing", "Bunch_Crossing", 20, -10., 10.);
   StripProf = booker.book1D("Strip_Profile", "Strip_Profile", 100, 0, 100);
 
   BxDisc_4Plus = booker.book1D("BxDisc_4Plus", "BxDisc_4Plus", 20, -10., 10.);
   BxDisc_4Min = booker.book1D("BxDisc_4Min", "BxDisc_4Min", 20, -10., 10.);
 
   // endcap residuals
   ResDmin1 = booker.book1D("Disk_Min1_Residuals", "Disk_Min1_Residuals", 400, -8, 8);
   ResDmin2 = booker.book1D("Disk_Min2_Residuals", "Disk_Min2_Residuals", 400, -8, 8);
   ResDmin3 = booker.book1D("Disk_Min3_Residuals", "Disk_Min3_Residuals", 400, -8, 8);
   ResDplu1 = booker.book1D("Disk_Plu1_Residuals", "Disk_Plu1_Residuals", 400, -8, 8);
   ResDplu2 = booker.book1D("Disk_Plu2_Residuals", "Disk_Plu2_Residuals", 400, -8, 8);
   ResDplu3 = booker.book1D("Disk_Plu3_Residuals", "Disk_Plu3_Residuals", 400, -8, 8);
 
   ResDmin4 = booker.book1D("Disk_Min4_Residuals", "Disk_Min4_Residuals", 400, -8, 8);
   ResDplu4 = booker.book1D("Disk_Plu4_Residuals", "Disk_Plu4_Residuals", 400, -8, 8);
 
   Res_Endcap1_Ring2_A = booker.book1D("Res_Endcap1_Ring2_A", "Res_Endcap1_Ring2_A", 400, -8, 8);
   Res_Endcap1_Ring2_B = booker.book1D("Res_Endcap1_Ring2_B", "Res_Endcap1_Ring2_B", 400, -8, 8);
   Res_Endcap1_Ring2_C = booker.book1D("Res_Endcap1_Ring2_C", "Res_Endcap1_Ring2_C", 400, -8, 8);
 
   Res_Endcap23_Ring2_A = booker.book1D("Res_Endcap23_Ring2_A", "Res_Endcap23_Ring2_A", 400, -8, 8);
   Res_Endcap23_Ring2_B = booker.book1D("Res_Endcap23_Ring2_B", "Res_Endcap23_Ring2_B", 400, -8, 8);
   Res_Endcap23_Ring2_C = booker.book1D("Res_Endcap23_Ring2_C", "Res_Endcap23_Ring2_C", 400, -8, 8);
 
   Res_Endcap123_Ring3_A = booker.book1D("Res_Endcap123_Ring3_A", "Res_Endcap123_Ring3_A", 400, -8, 8);
   Res_Endcap123_Ring3_B = booker.book1D("Res_Endcap123_Ring3_B", "Res_Endcap123_Ring3_B", 400, -8, 8);
   Res_Endcap123_Ring3_C = booker.book1D("Res_Endcap123_Ring3_C", "Res_Endcap123_Ring3_C", 400, -8, 8);
 
   // Timing informations
   hDigiTimeAll =
       booker.book1D("DigiTimeAll", "Digi time including present electronics;Digi time (ns)", 100, -12.5, 12.5);
   hDigiTime = booker.book1D("DigiTime", "Digi time only with timing information;Digi time (ns)", 100, -12.5, 12.5);
   hDigiTimeIRPC = booker.book1D("DigiTimeIRPC", "IRPC Digi time;Digi time (ns)", 100, -12.5, 12.5);
   hDigiTimeNoIRPC = booker.book1D("DigiTimeNoIRPC", "non-IRPC Digi time;Digi time (ns)", 100, -12.5, 12.5);
 }

Member Data Documentation

MonitorElement* RPCDigiValid::BxDisc_4Min

private

Definition at line 71 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::BxDisc_4Plus

private

Definition at line 70 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::BxDist

private

Definition at line 35 of file RPCDigiValid.h.

std::string RPCDigiValid::digiLabel

private

Definition at line 79 of file RPCDigiValid.h.

MonitorElement * RPCDigiValid::hDigiTime

private

Definition at line 76 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::hDigiTimeAll

private

Definition at line 76 of file RPCDigiValid.h.

MonitorElement * RPCDigiValid::hDigiTimeIRPC

private

Definition at line 76 of file RPCDigiValid.h.

MonitorElement * RPCDigiValid::hDigiTimeNoIRPC

private

Definition at line 76 of file RPCDigiValid.h.

std::string RPCDigiValid::outputFile_

private

Definition at line 78 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res

private

Definition at line 29 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res_Endcap123_Ring3_A

private

Definition at line 63 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res_Endcap123_Ring3_B

private

Definition at line 64 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res_Endcap123_Ring3_C

private

Definition at line 65 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res_Endcap1_Ring2_A

private

Definition at line 55 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res_Endcap1_Ring2_B

private

Definition at line 56 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res_Endcap1_Ring2_C

private

Definition at line 57 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res_Endcap23_Ring2_A

private

Definition at line 59 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res_Endcap23_Ring2_B

private

Definition at line 60 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::Res_Endcap23_Ring2_C

private

Definition at line 61 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResDmin1

private

Definition at line 47 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResDmin2

private

Definition at line 48 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResDmin3

private

Definition at line 49 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResDmin4

private

Definition at line 68 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResDplu1

private

Definition at line 50 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResDplu2

private

Definition at line 51 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResDplu3

private

Definition at line 52 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResDplu4

private

Definition at line 69 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResLayer1_barrel

private

Definition at line 39 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResLayer2_barrel

private

Definition at line 40 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResLayer3_barrel

private

Definition at line 41 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResLayer4_barrel

private

Definition at line 42 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResLayer5_barrel

private

Definition at line 43 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResLayer6_barrel

private

Definition at line 44 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResWmin1

private

Definition at line 31 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResWmin2

private

Definition at line 30 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResWplu1

private

Definition at line 33 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResWplu2

private

Definition at line 34 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::ResWzer0

private

Definition at line 32 of file RPCDigiValid.h.

edm::EDGetTokenT<RPCDigiCollection> RPCDigiValid::rpcDigiToken

private

Definition at line 83 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::rzview

private

Definition at line 28 of file RPCDigiValid.h.

edm::EDGetTokenT<edm::PSimHitContainer> RPCDigiValid::simHitToken

private

Definition at line 82 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::StripProf

private

Definition at line 36 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::xyvDmin4

private

Definition at line 73 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::xyvDplu4

private

Definition at line 72 of file RPCDigiValid.h.

MonitorElement* RPCDigiValid::xyview

private

Definition at line 27 of file RPCDigiValid.h.

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation