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 DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

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

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

void	beginLuminosityBlock (const edm::LuminosityBlock &, 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	endLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) override

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

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 42 of file TotemTimingDQMSource.cc.

Constructor & Destructor Documentation

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

Definition at line 409 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 431 of file TotemTimingDQMSource.cc.

431 {}

Member Function Documentation

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

overrideprotected

Definition at line 500 of file TotemTimingDQMSource.cc.

References edm::EventBase::bunchCrossing(), channelPlots_, MillePedeFileConverter_cfg::e, geometry, edm::EventSetup::get(), CTPPSGeometry::getSensor(), mps_fire::i, 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;
 
   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_));
         ++(channelPlots_[detId].hitsCounterPerLumisection);
       }
     }
   }
   // 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::beginLuminosityBlock	(	const edm::LuminosityBlock &	,
		const edm::EventSetup &
	)

overrideprotected

Definition at line 491 of file TotemTimingDQMSource.cc.

References plotFactory::plot, and potPlots_.

 {
   for ( auto& plot : potPlots_ )
     plot.second.hitDistribution2d_lumisection->Reset();
 }

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

overrideprotected

Definition at line 459 of file TotemTimingDQMSource.cc.

References DQMStore::IBooker::cd(), channelPlots_, CTPPS_NUM_OF_ARMS, globalPlot_, planePlots_, potPlots_, 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 435 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::endLuminosityBlock	(	const edm::LuminosityBlock &	,
		const edm::EventSetup &
	)

overrideprotected

Definition at line 808 of file TotemTimingDQMSource.cc.

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

                                                                      {
   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());
         plot.second.hitRate->Fill(
             chId.plane(), chId.channel(),
             (double)chPlot.second.hitsCounterPerLumisection * HIT_RATE_FACTOR);
       }
     }
   }
 
   for (auto &plot : channelPlots_) {
     if (plot.second.hitsCounterPerLumisection != 0) {
       plot.second.hitRate->Fill((double)plot.second.hitsCounterPerLumisection *
                                 HIT_RATE_FACTOR);
     }
     plot.second.hitsCounterPerLumisection = 0;
   }
 }

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

overrideprotected

Definition at line 860 of file TotemTimingDQMSource.cc.

References DEFINE_FWK_MODULE.

860 {}

Member Data Documentation

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

private

Definition at line 204 of file TotemTimingDQMSource.cc.

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

const float TotemTimingDQMSource::COS_8_DEG = 0.990268

staticprivate

Definition at line 90 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::CTPPS_FAR_RP_ID = 3

staticprivate

Definition at line 85 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::CTPPS_NEAR_RP_ID = 2

staticprivate

Definition at line 84 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::CTPPS_NUM_OF_ARMS = 2

staticprivate

Definition at line 76 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms().

const double TotemTimingDQMSource::DISPLAY_RESOLUTION_FOR_HITS_MM = 0.1

staticprivate

Definition at line 68 of file TotemTimingDQMSource.cc.

const double TotemTimingDQMSource::DQM_FRACTION_OF_EVENTS = 1.

staticprivate

Definition at line 65 of file TotemTimingDQMSource.cc.

GlobalPlots TotemTimingDQMSource::globalPlot_

private

Definition at line 116 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms(), and endLuminosityBlock().

const double TotemTimingDQMSource::HIT_RATE_FACTOR

staticprivate

Initial value:

=

DQM_FRACTION_OF_EVENTS / SEC_PER_LUMI_SECTION

Definition at line 67 of file TotemTimingDQMSource.cc.

Referenced by endLuminosityBlock().

const double TotemTimingDQMSource::INV_DISPLAY_RESOLUTION_FOR_HITS_MM

staticprivate

Initial value:

=

1. / DISPLAY_RESOLUTION_FOR_HITS_MM

Definition at line 71 of file TotemTimingDQMSource.cc.

Referenced by analyze(), TotemTimingDQMSource::PlanePlots::PlanePlots(), and TotemTimingDQMSource::PotPlots::PotPlots().

const double TotemTimingDQMSource::LHC_CLOCK_PERIOD_NS = 24.95

staticprivate

Definition at line 64 of file TotemTimingDQMSource.cc.

Referenced by analyze().

double TotemTimingDQMSource::maximumStripAngleForTomography_

private

Definition at line 100 of file TotemTimingDQMSource.cc.

Referenced by analyze().

double TotemTimingDQMSource::minimumStripAngleForTomography_

private

Definition at line 99 of file TotemTimingDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 176 of file TotemTimingDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 161 of file TotemTimingDQMSource.cc.

Referenced by analyze(), beginLuminosityBlock(), bookHistograms(), and endLuminosityBlock().

const double TotemTimingDQMSource::SAMPIC_ADC_V = 1./256

staticprivate

Definition at line 75 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const double TotemTimingDQMSource::SAMPIC_MAX_NUMBER_OF_SAMPLES = 64

staticprivate

Definition at line 74 of file TotemTimingDQMSource.cc.

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

staticprivate

Definition at line 73 of file TotemTimingDQMSource.cc.

unsigned int TotemTimingDQMSource::samplesForNoise_

private

Definition at line 101 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const double TotemTimingDQMSource::SEC_PER_LUMI_SECTION = 23.31

staticprivate

Definition at line 61 of file TotemTimingDQMSource.cc.

const float TotemTimingDQMSource::SIN_8_DEG = -0.139173

staticprivate

Definition at line 91 of file TotemTimingDQMSource.cc.

edm::TimeValue_t TotemTimingDQMSource::timeOfPreviousEvent_

private

Definition at line 103 of file TotemTimingDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 94 of file TotemTimingDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 97 of file TotemTimingDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 93 of file TotemTimingDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 95 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const double TotemTimingDQMSource::TOMOGRAPHY_RESOLUTION_MM = 1

staticprivate

Definition at line 72 of file TotemTimingDQMSource.cc.

Referenced by TotemTimingDQMSource::ChannelPlots::ChannelPlots(), and TotemTimingDQMSource::PotPlots::PotPlots().

const int TotemTimingDQMSource::TOTEM_STATION_210 = 0

staticprivate

Definition at line 78 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const int TotemTimingDQMSource::TOTEM_STATION_220 = 2

staticprivate

Definition at line 79 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const int TotemTimingDQMSource::TOTEM_STRIP_MAX_RP_ID = 5

staticprivate

Definition at line 83 of file TotemTimingDQMSource.cc.

Referenced by analyze().

const int TotemTimingDQMSource::TOTEM_STRIP_MIN_RP_ID = 4

staticprivate

Definition at line 82 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::TOTEM_TIMING_BOT_RP_ID = 1

staticprivate

Definition at line 81 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 88 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::TOTEM_TIMING_FED_ID_56

staticprivate

Initial value:

=

FEDNumbering::MINTotemRPTimingVerticalFEDID

Definition at line 89 of file TotemTimingDQMSource.cc.

const int TotemTimingDQMSource::TOTEM_TIMING_NUM_OF_CHANNELS = 12

staticprivate

Definition at line 87 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms().

const int TotemTimingDQMSource::TOTEM_TIMING_NUM_OF_PLANES = 4

staticprivate

Definition at line 86 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms().

const int TotemTimingDQMSource::TOTEM_TIMING_STATION_ID = 2

staticprivate

Definition at line 77 of file TotemTimingDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

const int TotemTimingDQMSource::TOTEM_TIMING_TOP_RP_ID = 0

staticprivate

Definition at line 80 of file TotemTimingDQMSource.cc.

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

unsigned int TotemTimingDQMSource::verbosity_

private

Definition at line 102 of file TotemTimingDQMSource.cc.

Referenced by analyze().

float TotemTimingDQMSource::verticalShiftBot_

private

Definition at line 105 of file TotemTimingDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

float TotemTimingDQMSource::verticalShiftTop_

private

Definition at line 105 of file TotemTimingDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

Classes

Public Member Functions

Protected Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation