Inheritance diagram for TotemTimingDQMSource:

Classes
struct	ChannelPlots
	plots related to one Diamond channel More...

struct	GlobalPlots
	plots related to the whole system More...

struct	PlanePlots
	plots related to one Diamond plane More...

struct	PotPlots
	plots related to one Diamond detector package More...

Public Member Functions
	TotemTimingDQMSource (const edm::ParameterSet &)

	~TotemTimingDQMSource () override

Public Member Functions inherited from one::DQMEDAnalyzer< edm::LuminosityBlockCache< totemds::Cache > >
	DQMEDAnalyzer ()=default

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

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

	~DQMEDAnalyzer () override=default

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

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

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override

void	endRun (const edm::Run &, const edm::EventSetup &) override

std::shared_ptr< totemds::Cache >	globalBeginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) const override

void	globalEndLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) override

Private Attributes
std::unordered_map< unsigned int, ChannelPlots >	channelPlots_

GlobalPlots	globalPlot_

double	maximumStripAngleForTomography_

double	minimumStripAngleForTomography_

std::unordered_map< unsigned int, PlanePlots >	planePlots_

std::unordered_map< unsigned int, PotPlots >	potPlots_

unsigned int	samplesForNoise_

edm::TimeValue_t	timeOfPreviousEvent_

edm::EDGetTokenT< edm::DetSetVector< TotemTimingDigi > >	tokenDigi_

edm::EDGetTokenT< std::vector< TotemFEDInfo > >	tokenFEDInfo_

edm::EDGetTokenT< edm::DetSetVector< TotemRPLocalTrack > >	tokenLocalTrack_

edm::EDGetTokenT< edm::DetSetVector< TotemTimingRecHit > >	tokenRecHit_

unsigned int	verbosity_

float	verticalShiftBot_

float	verticalShiftTop_

Static Private Attributes
static const float	COS_8_DEG = 0.990268

static const int	CTPPS_FAR_RP_ID = 3

static const int	CTPPS_NEAR_RP_ID = 2

static const int	CTPPS_NUM_OF_ARMS = 2

static const double	DISPLAY_RESOLUTION_FOR_HITS_MM = 0.1

static const double	DQM_FRACTION_OF_EVENTS = 1.

static const double	HIT_RATE_FACTOR

static const double	INV_DISPLAY_RESOLUTION_FOR_HITS_MM

static const double	LHC_CLOCK_PERIOD_NS = 24.95

static const double	SAMPIC_ADC_V = 1./256

static const double	SAMPIC_MAX_NUMBER_OF_SAMPLES = 64

static const double	SAMPIC_SAMPLING_PERIOD_NS = 1. / 7.8e9

static const double	SEC_PER_LUMI_SECTION = 23.31

static const float	SIN_8_DEG = -0.139173

static const double	TOMOGRAPHY_RESOLUTION_MM = 1

static const int	TOTEM_STATION_210 = 0

static const int	TOTEM_STATION_220 = 2

static const int	TOTEM_STRIP_MAX_RP_ID = 5

static const int	TOTEM_STRIP_MIN_RP_ID = 4

static const int	TOTEM_TIMING_BOT_RP_ID = 1

static const int	TOTEM_TIMING_FED_ID_45

static const int	TOTEM_TIMING_FED_ID_56

static const int	TOTEM_TIMING_NUM_OF_CHANNELS = 12

static const int	TOTEM_TIMING_NUM_OF_PLANES = 4

static const int	TOTEM_TIMING_STATION_ID = 2

static const int	TOTEM_TIMING_TOP_RP_ID = 0

Detailed Description

Definition at line 50 of file TotemTimingDQMSource.cc.

Constructor & Destructor Documentation

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

Definition at line 415 of file TotemTimingDQMSource.cc.

     : tokenLocalTrack_(consumes<edm::DetSetVector<TotemRPLocalTrack>>(
           ps.getParameter<edm::InputTag>("tagLocalTrack"))),
       tokenDigi_(consumes<edm::DetSetVector<TotemTimingDigi>>(
           ps.getParameter<edm::InputTag>("tagDigi"))),
       tokenRecHit_(consumes<edm::DetSetVector<TotemTimingRecHit>>(
           ps.getParameter<edm::InputTag>("tagRecHits"))),
       // tokenTrack_(consumes<edm::DetSetVector<TotemTimingLocalTrack>>(
       //     ps.getParameter<edm::InputTag>("tagLocalTracks"))),
       tokenFEDInfo_(consumes<std::vector<TotemFEDInfo>>(
           ps.getParameter<edm::InputTag>("tagFEDInfo"))),
       minimumStripAngleForTomography_(
           ps.getParameter<double>("minimumStripAngleForTomography")),
       maximumStripAngleForTomography_(
           ps.getParameter<double>("maximumStripAngleForTomography")),
       samplesForNoise_(
           ps.getUntrackedParameter<unsigned int>("samplesForNoise", 5)),
       verbosity_(ps.getUntrackedParameter<unsigned int>("verbosity", 0)),
       timeOfPreviousEvent_(0) {}

TotemTimingDQMSource::~TotemTimingDQMSource ( )

override

Definition at line 437 of file TotemTimingDQMSource.cc.

437 {}

Member Function Documentation

void TotemTimingDQMSource::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overrideprotected

Definition at line 510 of file TotemTimingDQMSource.cc.

References edm::EventBase::bunchCrossing(), channelPlots_, MillePedeFileConverter_cfg::e, geometry, edm::EventSetup::get(), edm::Event::getLuminosityBlock(), CTPPSGeometry::getSensor(), mps_fire::i, edm::LuminosityBlock::index(), INV_DISPLAY_RESOLUTION_FOR_HITS_MM, edm::HandleBase::isValid(), LHC_CLOCK_PERIOD_NS, maximumStripAngleForTomography_, minimumStripAngleForTomography_, TotemTimingRecHit::NO_T_AVAILABLE, planePlots_, potPlots_, TrackInfoProducer_cfi::rechits, SAMPIC_ADC_V, samplesForNoise_, TotemTimingDetId::setChannel(), TotemRPDetId::setPlane(), TotemTimingDetId::setPlane(), edm::EventBase::time(), timeOfPreviousEvent_, tokenDigi_, tokenFEDInfo_, tokenLocalTrack_, tokenRecHit_, TOTEM_STATION_210, TOTEM_STATION_220, TOTEM_STRIP_MAX_RP_ID, TOTEM_TIMING_BOT_RP_ID, TOTEM_TIMING_TOP_RP_ID, DetGeomDesc::translation(), edm::Timestamp::value(), verbosity_, verticalShiftBot_, verticalShiftTop_, x, and y.

 {
   // get event setup data
   edm::ESHandle<CTPPSGeometry> geometry;
   eventSetup.get<VeryForwardRealGeometryRecord>().get(geometry);
 
   // get event data
   edm::Handle<edm::DetSetVector<TotemRPLocalTrack>> stripTracks;
   event.getByToken(tokenLocalTrack_, stripTracks);
 
   edm::Handle<edm::DetSetVector<TotemTimingDigi>> timingDigis;
   event.getByToken(tokenDigi_, timingDigis);
 
   edm::Handle<std::vector<TotemFEDInfo>> fedInfo;
   event.getByToken(tokenFEDInfo_, fedInfo);
 
   edm::Handle<edm::DetSetVector<TotemTimingRecHit>> timingRecHits;
   event.getByToken(tokenRecHit_, timingRecHits);
 
   // edm::Handle<edm::DetSetVector<TotemTimingLocalTrack>> timingLocalTracks;
   // event.getByToken(timingLocalTracks, timingLocalTracks);
 
   // check validity
   bool valid = true;
   valid &= timingDigis.isValid();
   valid &= fedInfo.isValid();
 
   if (!valid)
   {
     if (verbosity_)
     {
       edm::LogProblem("TotemTimingDQMSource")
           << "ERROR in TotemTimingDQMSource::analyze > some of the required inputs "
              "are not valid. Skipping this event.\n"
           << "    timingDigis.isValid = " << timingDigis.isValid() << "\n"
           << "    fedInfo.isValid = " << fedInfo.isValid();
     }
 
     return;
   }
 
   // Using TotemTimingDigi
   std::set<uint8_t> boardSet;
   std::unordered_map<unsigned int, unsigned int> channelsPerPlane;
   std::unordered_map<unsigned int, unsigned int> channelsPerPlaneWithTime;
 
   auto lumiCache = luminosityBlockCache(event.getLuminosityBlock().index());
   for (const auto &digis : *timingDigis)
   {
     const TotemTimingDetId detId(digis.detId());
     TotemTimingDetId detId_pot(digis.detId());
     detId_pot.setPlane(0);
     detId_pot.setChannel(0);
     TotemTimingDetId detId_plane(digis.detId());
     detId_plane.setChannel(0);
 
     for (const auto &digi : digis)
     {
       // Pot Plots
       if (potPlots_.find(detId_pot) != potPlots_.end())
       {
         potPlots_[detId_pot].activityPerBX->Fill(event.bunchCrossing());
 
         potPlots_[detId_pot].digiDistribution->Fill(detId.plane(), detId.channel());
 
         for (auto it = digi.getSamplesBegin(); it != digi.getSamplesEnd(); ++it)
           potPlots_[detId_pot].dataSamplesRaw->Fill(*it);
 
         float boardId = digi.getEventInfo().getHardwareBoardId() +
                         0.5 * digi.getEventInfo().getHardwareSampicId();
         potPlots_[detId_pot].digiSent->Fill(boardId, digi.getHardwareChannelId());
         if (boardSet.find(digi.getEventInfo().getHardwareId()) == boardSet.end())
         {
           // This guarantees that every board is counted only once
           boardSet.insert(digi.getEventInfo().getHardwareId());
           std::bitset<16> chMap(digi.getEventInfo().getChannelMap());
           for (int i = 0; i < 16; ++i)
           {
             if (chMap.test(i))
             {
               potPlots_[detId_pot].digiAll->Fill(boardId, i);
             }
           }
         }
 
         potPlots_[detId_pot].planesWithDigisSet.insert(detId.plane());
       }
 
       // Plane Plots
       if (planePlots_.find(detId_plane) != planePlots_.end())
       {
         planePlots_[detId_plane].digiDistribution->Fill(detId.channel());
 
         if (channelsPerPlane.find(detId_plane) != channelsPerPlane.end())
           channelsPerPlane[detId_plane]++;
         else
           channelsPerPlane[detId_plane] = 0;
       }
 
       // Channel Plots
       if (channelPlots_.find(detId) != channelPlots_.end()) {
         channelPlots_[detId].activityPerBX->Fill(event.bunchCrossing());
 
         for (auto it = digi.getSamplesBegin(); it != digi.getSamplesEnd(); ++it)
           channelPlots_[detId].dataSamplesRaw->Fill(*it);
         for (unsigned short i = 0; i < samplesForNoise_; ++i)
           channelPlots_[detId].noiseSamples->Fill(SAMPIC_ADC_V *
                                                   digi.getSampleAt(i));
 
         unsigned int cellOfMax =
             std::max_element(digi.getSamplesBegin(), digi.getSamplesEnd()) -
             digi.getSamplesBegin();
         channelPlots_[detId].cellOfMax->Fill((int)cellOfMax);
 
         if (timeOfPreviousEvent_ != 0)
           channelPlots_[detId].hitTime->Fill(
               1e-3 * LHC_CLOCK_PERIOD_NS *
               (event.time().value() - timeOfPreviousEvent_));
         ++(lumiCache->hitsCounterMap[detId]);
       }
     }
   }
   // End digis
 
   for (const auto &rechits : *timingRecHits)
   {
     const TotemTimingDetId detId(rechits.detId());
     TotemTimingDetId detId_pot(rechits.detId());
     detId_pot.setPlane(0);
     detId_pot.setChannel(0);
     TotemTimingDetId detId_plane(rechits.detId());
     detId_plane.setChannel(0);
 
     for (const auto &rechit : rechits)
     {
       if (potPlots_.find(detId_pot) != potPlots_.end())
       {
         potPlots_[detId_pot].amplitude->Fill(rechit.getAmplitude());
 
 
         TH2F *hitHistoTmp = potPlots_[detId_pot].hitDistribution2d->getTH2F();
         TAxis *hitHistoTmpYAxis = hitHistoTmp->GetYaxis();
         float yCorrected = rechit.getY();
         yCorrected += (detId.rp()==TOTEM_TIMING_TOP_RP_ID)?
                       verticalShiftTop_ :
                       verticalShiftBot_;
         float x_shift = detId.plane();
         x_shift += (rechit.getX()>2)? 0.25 : 0;
         int startBin = hitHistoTmpYAxis->FindBin( yCorrected - 0.5*rechit.getYWidth() );
         int numOfBins = rechit.getYWidth()*INV_DISPLAY_RESOLUTION_FOR_HITS_MM;
         for ( int i=0; i<numOfBins; ++i)
         {
           potPlots_[detId_pot].hitDistribution2d->Fill(
                     detId.plane() + 0.25 * (rechit.getX()>2),
                     hitHistoTmpYAxis->GetBinCenter(startBin+i) );
           potPlots_[detId_pot].hitDistribution2d_lumisection->Fill(
                     x_shift,
                     hitHistoTmpYAxis->GetBinCenter(startBin+i) );
         }
 
         //All plots with Time
         if ( rechit.getT() != TotemTimingRecHit::NO_T_AVAILABLE )
         {
           for ( int i=0; i<numOfBins; ++i)
             potPlots_[detId_pot].hitDistribution2dWithTime->Fill(
                      detId.plane() + 0.25 * (rechit.getX()>2),
                      hitHistoTmpYAxis->GetBinCenter(startBin+i) );
 
           potPlots_[detId_pot].recHitTime->Fill(rechit.getT());
           potPlots_[detId_pot].planesWithTimeSet.insert(detId.plane());
 
           // Plane Plots
           if (planePlots_.find(detId_plane) != planePlots_.end())
           {
             // Visualization tricks
             float x_shift = (rechit.getX()>2)? 15 : 0;
             TH1F *hitProfileHistoTmp = planePlots_[detId_plane].hitProfile->getTH1F();
             int numOfBins = rechit.getYWidth()*INV_DISPLAY_RESOLUTION_FOR_HITS_MM;
             if (detId.rp()==TOTEM_TIMING_TOP_RP_ID)
             {
               float yCorrected = rechit.getY()+verticalShiftTop_-0.5*rechit.getYWidth()+x_shift;
               int startBin = hitProfileHistoTmp->FindBin(yCorrected);
               for ( int i=0; i<numOfBins; ++i)
                 hitProfileHistoTmp->Fill( hitProfileHistoTmp->GetBinCenter(startBin+i) );
             }
             else
             {
               float yCorrected = rechit.getY()+verticalShiftBot_+0.5*rechit.getYWidth()+(15-x_shift);
               int startBin = hitProfileHistoTmp->FindBin(yCorrected);
               int totBins = hitProfileHistoTmp->GetNbinsX();
               for ( int i=0; i<numOfBins; ++i)
                 hitProfileHistoTmp->Fill( hitProfileHistoTmp->GetBinCenter(totBins-startBin+i) );
             }
 
             if (channelsPerPlaneWithTime.find(detId_plane) != channelsPerPlaneWithTime.end())
               channelsPerPlaneWithTime[detId_plane]++;
             else
               channelsPerPlaneWithTime[detId_plane] = 0;
           }
 
           if (channelPlots_.find(detId) != channelPlots_.end())
           {
             potPlots_[detId_pot].tirggerCellTime->Fill(rechit.getSampicThresholdTime());
             channelPlots_[detId].tirggerCellTime->Fill(rechit.getSampicThresholdTime());
             channelPlots_[detId].recHitTime->Fill(rechit.getT());
             channelPlots_[detId].amplitude->Fill(rechit.getAmplitude());
           }
         }
       }
     }
   }
   // End RecHits
 
   // Tomography of timing using strips
   for (const auto &rechits : *timingRecHits)
   {
     const TotemTimingDetId detId(rechits.detId());
     TotemTimingDetId detId_pot(rechits.detId());
     detId_pot.setPlane(0);
     detId_pot.setChannel(0);
     TotemTimingDetId detId_plane(rechits.detId());
     detId_plane.setChannel(0);
 
     float y_shift=(detId.rp()==TOTEM_TIMING_TOP_RP_ID)?20:5;
 
     for (const auto &rechit : rechits)
     {
       if (rechit.getT() != TotemTimingRecHit::NO_T_AVAILABLE        &&
                 potPlots_.find(detId_pot) != potPlots_.end()        &&
                 planePlots_.find(detId_plane) != planePlots_.end()  &&
                 channelPlots_.find(detId) != channelPlots_.end()      )
       {
         if ( stripTracks.isValid() )
         {
           for (const auto &ds : *stripTracks) {
             const CTPPSDetId stripId(ds.detId());
             // mean position of U and V planes
             TotemRPDetId plId_V(stripId);
             plId_V.setPlane(0);
             TotemRPDetId plId_U(stripId);
             plId_U.setPlane(1);
 
             double rp_x = 0;
             double rp_y = 0;
             try
             {
               rp_x = (geometry->getSensor(plId_V)->translation().x() +
                       geometry->getSensor(plId_U)->translation().x())/2;
               rp_y = (geometry->getSensor(plId_V)->translation().y() +
                       geometry->getSensor(plId_U)->translation().y())/2;
             }
             catch(const cms::Exception &)
             {
               continue;
             }
 
             for (const auto &striplt : ds) {
               if (striplt.isValid() && stripId.arm() == detId.arm()) {
                 if (striplt.getTx() > maximumStripAngleForTomography_ ||
                     striplt.getTy() > maximumStripAngleForTomography_)
                   continue;
                 if (striplt.getTx() < minimumStripAngleForTomography_ ||
                     striplt.getTy() < minimumStripAngleForTomography_)
                   continue;
                 if (stripId.rp() - detId.rp() ==
                     (TOTEM_STRIP_MAX_RP_ID - TOTEM_TIMING_BOT_RP_ID))
                 {
                   double x = striplt.getX0() - rp_x;
                   double y = striplt.getY0() - rp_y;
                   if (stripId.station() == TOTEM_STATION_210)
                   {
                     potPlots_[detId_pot].stripTomography210->Fill(x + detId.plane()*50, y + y_shift);
                     channelPlots_[detId].stripTomography210->Fill(x, y + y_shift);
                   }
                   else if (stripId.station() == TOTEM_STATION_220)
                   {
                     potPlots_[detId_pot].stripTomography220->Fill(x+detId.plane()*50, y + y_shift);
                     channelPlots_[detId].stripTomography220->Fill(x, y + y_shift);
                   }
                 }
               }
             }
           }
         }
       }
     }
   }
 
 
   for (auto &plt : potPlots_) {
     plt.second.planesWithDigis->Fill(plt.second.planesWithDigisSet.size());
     plt.second.planesWithDigisSet.clear();
     plt.second.planesWithTime->Fill(plt.second.planesWithTimeSet.size());
     plt.second.planesWithTimeSet.clear();
   }
 
   for (const auto &plt : channelsPerPlane) {
     planePlots_[plt.first].hitMultiplicity->Fill(plt.second);
   }
   for (const auto &plt : channelsPerPlaneWithTime) {
     planePlots_[plt.first].hitMultiplicityWithTime->Fill(plt.second);
   }
 
   timeOfPreviousEvent_ = event.time().value();
 }

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

overrideprotected

Definition at line 465 of file TotemTimingDQMSource.cc.

References DQMStore::IBooker::cd(), channelPlots_, CTPPS_NUM_OF_ARMS, globalPlot_, planePlots_, potPlots_, year_2016_postTS2_cff::rpId, DQMStore::IBooker::setCurrentFolder(), TOTEM_TIMING_BOT_RP_ID, TOTEM_TIMING_NUM_OF_CHANNELS, TOTEM_TIMING_NUM_OF_PLANES, TOTEM_TIMING_STATION_ID, and TOTEM_TIMING_TOP_RP_ID.

 {
   ibooker.cd();
   ibooker.setCurrentFolder("CTPPS");
 
   globalPlot_ = GlobalPlots(ibooker);
 
   for (unsigned short arm = 0; arm < CTPPS_NUM_OF_ARMS; ++arm)
   {
     for (unsigned short rp = TOTEM_TIMING_TOP_RP_ID;
          rp <= TOTEM_TIMING_BOT_RP_ID; ++rp)
       {
       const TotemTimingDetId rpId(arm, TOTEM_TIMING_STATION_ID, rp);
       potPlots_[rpId] = PotPlots(ibooker, rpId);
       for (unsigned short pl = 0; pl < TOTEM_TIMING_NUM_OF_PLANES; ++pl)
       {
         const TotemTimingDetId plId(arm, TOTEM_TIMING_STATION_ID, rp, pl);
         planePlots_[plId] = PlanePlots(ibooker, plId);
         for (unsigned short ch = 0; ch < TOTEM_TIMING_NUM_OF_CHANNELS; ++ch)
         {
           const TotemTimingDetId chId(arm, TOTEM_TIMING_STATION_ID, rp, pl, ch);
           channelPlots_[chId] = ChannelPlots(ibooker, chId);
         }
       }
     }
   }
 }

void TotemTimingDQMSource::dqmBeginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overrideprotected

Definition at line 441 of file TotemTimingDQMSource.cc.

References relativeConstraints::geom, edm::EventSetup::get(), CTPPSGeometry::getSensorNoThrow(), DetGeomDesc::params(), edm::ESHandle< T >::product(), TOTEM_TIMING_BOT_RP_ID, TOTEM_TIMING_STATION_ID, TOTEM_TIMING_TOP_RP_ID, DetGeomDesc::translation(), verticalShiftBot_, and verticalShiftTop_.

 {
   // Get detector shifts from the geometry (if present)
   edm::ESHandle<CTPPSGeometry> geometry_;
   iSetup.get<VeryForwardRealGeometryRecord>().get( geometry_ );
   const CTPPSGeometry *geom = geometry_.product();
   const TotemTimingDetId detid_top(0, TOTEM_TIMING_STATION_ID, TOTEM_TIMING_BOT_RP_ID, 0, 0 );
   const TotemTimingDetId detid_bot(0, TOTEM_TIMING_STATION_ID, TOTEM_TIMING_TOP_RP_ID, 0, 7 );
   verticalShiftTop_ = 0;
   verticalShiftBot_ = 0;
   {
     const DetGeomDesc* det_top = geom->getSensorNoThrow( detid_top );
     if(det_top) {
       verticalShiftTop_ = det_top->translation().y() + det_top->params().at( 1 );
     }
     const DetGeomDesc* det_bot = geom->getSensorNoThrow( detid_bot );
     if(det_bot)
       verticalShiftBot_ = det_bot->translation().y() + det_bot->params().at( 1 );
   }
 }

void TotemTimingDQMSource::endRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

overrideprotected

Definition at line 877 of file TotemTimingDQMSource.cc.

References DEFINE_FWK_MODULE.

877 {}

std::shared_ptr< totemds::Cache > TotemTimingDQMSource::globalBeginLuminosityBlock	(	const edm::LuminosityBlock &	,
		const edm::EventSetup &
	)		const

overrideprotected

Definition at line 498 of file TotemTimingDQMSource.cc.

References edmIntegrityCheck::d, plotFactory::plot, and potPlots_.

 {
   auto d = std::make_shared<totemds::Cache>();
   d->hitDistribution2dMap.reserve(potPlots_.size());
   for ( auto& plot : potPlots_ )
     d->hitDistribution2dMap[plot.first] = std::unique_ptr<TH2F>(static_cast<TH2F*>(plot.second.hitDistribution2d_lumisection->getTH2F()->Clone()));
   return d;
 }

void TotemTimingDQMSource::globalEndLuminosityBlock	(	const edm::LuminosityBlock &	iLumi,
		const edm::EventSetup &
	)

overrideprotected

Definition at line 819 of file TotemTimingDQMSource.cc.

References CTPPSDetId::arm(), TotemTimingDetId::channel(), channelPlots_, TotemTimingDQMSource::GlobalPlots::digiSentPercentage, MonitorElement::getTH2F(), globalPlot_, HIT_RATE_FACTOR, edm::LuminosityBlock::index(), TotemTimingDetId::plane(), plotFactory::plot, potPlots_, MonitorElement::Reset(), CTPPSDetId::rp(), and year_2016_postTS2_cff::rpId.

                                                                            {
   auto lumiCache = luminosityBlockCache(iLumi.index());
   for ( auto& plot : potPlots_ ) {
     *(plot.second.hitDistribution2d_lumisection->getTH2F())=*(lumiCache->hitDistribution2dMap[plot.first]);
   }
 
   globalPlot_.digiSentPercentage->Reset();
   TH2F *hitHistoGlobalTmp = globalPlot_.digiSentPercentage->getTH2F();
   for (auto &plot : potPlots_) {
     TH2F *hitHistoTmp = plot.second.digiSentPercentage->getTH2F();
     TH2F *histoSent = plot.second.digiSent->getTH2F();
     TH2F *histoAll = plot.second.digiAll->getTH2F();
 
     hitHistoTmp->Divide(histoSent, histoAll);
     hitHistoTmp->Scale(100);
     hitHistoGlobalTmp->Add(hitHistoTmp,1);
 
     plot.second.baseline->Reset();
     plot.second.noiseRMS->Reset();
     plot.second.meanAmplitude->Reset();
     plot.second.cellOfMax->Reset();
     plot.second.hitRate->Reset();
     TotemTimingDetId rpId(plot.first);
     for (auto &chPlot : channelPlots_) {
       TotemTimingDetId chId(chPlot.first);
       if (chId.arm() == rpId.arm() && chId.rp() == rpId.rp()) {
         plot.second.baseline->Fill(
             chId.plane(), chId.channel(),
             chPlot.second.noiseSamples->getTH1F()->GetMean());
         plot.second.noiseRMS->Fill(
             chId.plane(), chId.channel(),
             chPlot.second.noiseSamples->getTH1F()->GetRMS());
         plot.second.meanAmplitude->Fill(
             chId.plane(), chId.channel(),
             chPlot.second.amplitude->getTH1F()->GetMean());
         plot.second.cellOfMax->Fill(
             chId.plane(), chId.channel(),
             chPlot.second.cellOfMax->getTH1F()->GetMean());
         auto hitsCounterPerLumisection = lumiCache->hitsCounterMap[chPlot.first];
         plot.second.hitRate->Fill(
             chId.plane(), chId.channel(),
             (double)hitsCounterPerLumisection * HIT_RATE_FACTOR);
       }
     }
   }
 
   for (auto &plot : channelPlots_) {
         auto hitsCounterPerLumisection = lumiCache->hitsCounterMap[plot.first];
     if (hitsCounterPerLumisection != 0) {
       plot.second.hitRate->Fill((double)hitsCounterPerLumisection *
                                 HIT_RATE_FACTOR);
     }
   }
 }

Member Data Documentation

std::unordered_map<unsigned int, ChannelPlots> TotemTimingDQMSource::channelPlots_

private

Definition at line 211 of file TotemTimingDQMSource.cc.

Referenced by analyze(), bookHistograms(), and globalEndLuminosityBlock().

const float TotemTimingDQMSource::COS_8_DEG = 0.990268

staticprivate

Definition at line 98 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::CTPPS_FAR_RP_ID = 3

staticprivate

Definition at line 93 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::CTPPS_NEAR_RP_ID = 2

staticprivate

Definition at line 92 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::CTPPS_NUM_OF_ARMS = 2

staticprivate

Definition at line 84 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms().

const double TotemTimingDQMSource::DISPLAY_RESOLUTION_FOR_HITS_MM = 0.1

staticprivate

Definition at line 76 of file TotemTimingDQMSource.cc.

const double TotemTimingDQMSource::DQM_FRACTION_OF_EVENTS = 1.

staticprivate

Definition at line 73 of file TotemTimingDQMSource.cc.

GlobalPlots TotemTimingDQMSource::globalPlot_

private

Definition at line 124 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms(), and globalEndLuminosityBlock().

const double TotemTimingDQMSource::HIT_RATE_FACTOR

staticprivate

Initial value:

=

DQM_FRACTION_OF_EVENTS / SEC_PER_LUMI_SECTION

Definition at line 75 of file TotemTimingDQMSource.cc.

Referenced by globalEndLuminosityBlock().

const double TotemTimingDQMSource::INV_DISPLAY_RESOLUTION_FOR_HITS_MM

staticprivate

Initial value:

=

1. / DISPLAY_RESOLUTION_FOR_HITS_MM

Definition at line 79 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const double TotemTimingDQMSource::LHC_CLOCK_PERIOD_NS = 24.95

staticprivate

Definition at line 72 of file TotemTimingDQMSource.cc.

Referenced by analyze().

double TotemTimingDQMSource::maximumStripAngleForTomography_

private

Definition at line 108 of file TotemTimingDQMSource.cc.

Referenced by analyze().

double TotemTimingDQMSource::minimumStripAngleForTomography_

private

Definition at line 107 of file TotemTimingDQMSource.cc.

Referenced by analyze().

std::unordered_map<unsigned int, PlanePlots> TotemTimingDQMSource::planePlots_

private

Definition at line 184 of file TotemTimingDQMSource.cc.

Referenced by analyze(), and bookHistograms().

std::unordered_map<unsigned int, PotPlots> TotemTimingDQMSource::potPlots_

private

Definition at line 169 of file TotemTimingDQMSource.cc.

Referenced by analyze(), bookHistograms(), globalBeginLuminosityBlock(), and globalEndLuminosityBlock().

const double TotemTimingDQMSource::SAMPIC_ADC_V = 1./256

staticprivate

Definition at line 83 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const double TotemTimingDQMSource::SAMPIC_MAX_NUMBER_OF_SAMPLES = 64

staticprivate

Definition at line 82 of file TotemTimingDQMSource.cc.

const double TotemTimingDQMSource::SAMPIC_SAMPLING_PERIOD_NS = 1. / 7.8e9

staticprivate

Definition at line 81 of file TotemTimingDQMSource.cc.

unsigned int TotemTimingDQMSource::samplesForNoise_

private

Definition at line 109 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const double TotemTimingDQMSource::SEC_PER_LUMI_SECTION = 23.31

staticprivate

Definition at line 69 of file TotemTimingDQMSource.cc.

const float TotemTimingDQMSource::SIN_8_DEG = -0.139173

staticprivate

Definition at line 99 of file TotemTimingDQMSource.cc.

edm::TimeValue_t TotemTimingDQMSource::timeOfPreviousEvent_

private

Definition at line 111 of file TotemTimingDQMSource.cc.

Referenced by analyze().

edm::EDGetTokenT<edm::DetSetVector<TotemTimingDigi> > TotemTimingDQMSource::tokenDigi_

private

Definition at line 102 of file TotemTimingDQMSource.cc.

Referenced by analyze().

edm::EDGetTokenT<std::vector<TotemFEDInfo> > TotemTimingDQMSource::tokenFEDInfo_

private

Definition at line 105 of file TotemTimingDQMSource.cc.

Referenced by analyze().

edm::EDGetTokenT<edm::DetSetVector<TotemRPLocalTrack> > TotemTimingDQMSource::tokenLocalTrack_

private

Definition at line 101 of file TotemTimingDQMSource.cc.

Referenced by analyze().

edm::EDGetTokenT<edm::DetSetVector<TotemTimingRecHit> > TotemTimingDQMSource::tokenRecHit_

private

Definition at line 103 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const double TotemTimingDQMSource::TOMOGRAPHY_RESOLUTION_MM = 1

staticprivate

Definition at line 80 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::TOTEM_STATION_210 = 0

staticprivate

Definition at line 86 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const int TotemTimingDQMSource::TOTEM_STATION_220 = 2

staticprivate

Definition at line 87 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const int TotemTimingDQMSource::TOTEM_STRIP_MAX_RP_ID = 5

staticprivate

Definition at line 91 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const int TotemTimingDQMSource::TOTEM_STRIP_MIN_RP_ID = 4

staticprivate

Definition at line 90 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::TOTEM_TIMING_BOT_RP_ID = 1

staticprivate

Definition at line 89 of file TotemTimingDQMSource.cc.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

const int TotemTimingDQMSource::TOTEM_TIMING_FED_ID_45

staticprivate

Initial value:

=

FEDNumbering::MAXTotemRPTimingVerticalFEDID

Definition at line 96 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::TOTEM_TIMING_FED_ID_56

staticprivate

Initial value:

=

FEDNumbering::MINTotemRPTimingVerticalFEDID

Definition at line 97 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::TOTEM_TIMING_NUM_OF_CHANNELS = 12

staticprivate

Definition at line 95 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms().

const int TotemTimingDQMSource::TOTEM_TIMING_NUM_OF_PLANES = 4

staticprivate

Definition at line 94 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms().

const int TotemTimingDQMSource::TOTEM_TIMING_STATION_ID = 2

staticprivate

Definition at line 85 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

const int TotemTimingDQMSource::TOTEM_TIMING_TOP_RP_ID = 0

staticprivate

Definition at line 88 of file TotemTimingDQMSource.cc.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

unsigned int TotemTimingDQMSource::verbosity_

private

Definition at line 110 of file TotemTimingDQMSource.cc.

Referenced by analyze().

float TotemTimingDQMSource::verticalShiftBot_

private

Definition at line 113 of file TotemTimingDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

float TotemTimingDQMSource::verticalShiftTop_

private

Definition at line 113 of file TotemTimingDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

Classes

Public Member Functions

Protected Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation