Inheritance diagram for CTPPSDiamondDQMSource:

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
	CTPPSDiamondDQMSource (const edm::ParameterSet &)

virtual	~CTPPSDiamondDQMSource ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

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

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

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

	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 &	itemsToGetFromEvent () 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

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

virtual	~EDConsumerBase () noexcept(false)

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

void	beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &)

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 &)

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

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_

int	EC_difference_45_

int	EC_difference_56_

bool	excludeMultipleHits_

GlobalPlots	globalPlot_

double	minimumStripAngleForTomography_

std::unordered_map< unsigned int, PlanePlots >	planePlots_

std::unordered_map< unsigned int, PotPlots >	potPlots_

edm::EDGetTokenT< edm::DetSetVector< CTPPSDiamondRecHit > >	tokenDiamondHit_

edm::EDGetTokenT< edm::DetSetVector< CTPPSDiamondLocalTrack > >	tokenDiamondTrack_

edm::EDGetTokenT< edm::DetSetVector< CTPPSDiamondDigi > >	tokenDigi_

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

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

edm::EDGetTokenT< edm::DetSetVector< TotemVFATStatus > >	tokenStatus_

unsigned int	verbosity_

Static Private Attributes
static const int	CHANNEL_OF_VFAT_CLOCK = 30

static const int	CTPPS_DIAMOND_NUM_OF_CHANNELS = 12

static const int	CTPPS_DIAMOND_NUM_OF_PLANES = 4

static const int	CTPPS_DIAMOND_RP_ID = 6

static const int	CTPPS_DIAMOND_STATION_ID = 1

static const int	CTPPS_FAR_RP_ID = 3

static const int	CTPPS_FED_ID_45 = 583

static const int	CTPPS_FED_ID_56 = 582

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	HPTDC_BIN_WIDTH_NS = 25./1024

static const double	SEC_PER_LUMI_SECTION = 23.31

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 37 of file CTPPSDiamondDQMSource.cc.

Constructor & Destructor Documentation

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

Definition at line 319 of file CTPPSDiamondDQMSource.cc.

                                                                         :
   tokenStatus_      ( consumes< edm::DetSetVector<TotemVFATStatus> >       ( ps.getParameter<edm::InputTag>( "tagStatus" ) ) ),
   tokenLocalTrack_  ( consumes< edm::DetSetVector<TotemRPLocalTrack> >     ( ps.getParameter<edm::InputTag>( "tagLocalTrack" ) ) ),
   tokenDigi_        ( consumes< edm::DetSetVector<CTPPSDiamondDigi> >      ( ps.getParameter<edm::InputTag>( "tagDigi" ) ) ),
   tokenDiamondHit_  ( consumes< edm::DetSetVector<CTPPSDiamondRecHit> >    ( ps.getParameter<edm::InputTag>( "tagDiamondRecHits" ) ) ),
   tokenDiamondTrack_( consumes< edm::DetSetVector<CTPPSDiamondLocalTrack> >( ps.getParameter<edm::InputTag>( "tagDiamondLocalTracks" ) ) ),
   tokenFEDInfo_     ( consumes< std::vector<TotemFEDInfo> >                ( ps.getParameter<edm::InputTag>( "tagFEDInfo" ) ) ),
   excludeMultipleHits_           ( ps.getParameter<bool>( "excludeMultipleHits" ) ),
   minimumStripAngleForTomography_( ps.getParameter<double>( "minimumStripAngleForTomography" ) ),
   verbosity_                     ( ps.getUntrackedParameter<unsigned int>( "verbosity", 0 ) ),
   EC_difference_56_( -500 ), EC_difference_45_( -500 )
 {}

CTPPSDiamondDQMSource::~CTPPSDiamondDQMSource ( )

virtual

Definition at line 334 of file CTPPSDiamondDQMSource.cc.

335 {}

Member Function Documentation

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

protected

Definition at line 384 of file CTPPSDiamondDQMSource.cc.

References CTPPSDetId::arm(), edm::EventBase::bunchCrossing(), CHANNEL_OF_VFAT_CLOCK, channelPlots_, CTPPS_FAR_RP_ID, CTPPS_FED_ID_45, CTPPS_FED_ID_56, CTPPS_NEAR_RP_ID, TauDecayModes::dec, DISPLAY_RESOLUTION_FOR_HITS_MM, EC_difference_45_, EC_difference_56_, excludeMultipleHits_, MonitorElement::Fill(), HPTDCErrorFlags::getErrorId(), globalPlot_, CTPPSDiamondDQMSource::GlobalPlots::h_trackCorr_hor, HPTDC_BIN_WIDTH_NS, mps_fire::i, edm::HandleBase::isValid(), minimumStripAngleForTomography_, planePlots_, potPlots_, TrackInfoProducer_cfi::rechits, CTPPSDiamondDetId::setChannel(), CTPPSDiamondDetId::setPlane(), mps_update::status, tokenDiamondHit_, tokenDiamondTrack_, tokenDigi_, tokenFEDInfo_, tokenLocalTrack_, tokenStatus_, HiIsolationCommonParameters_cff::track, l1t::tracks, and verbosity_.

 {
   // get event data
   edm::Handle< edm::DetSetVector<TotemVFATStatus> > diamondVFATStatus;
   event.getByToken( tokenStatus_, diamondVFATStatus );
 
   edm::Handle< edm::DetSetVector<TotemRPLocalTrack> > stripTracks;
   event.getByToken( tokenLocalTrack_, stripTracks );
 
   edm::Handle< edm::DetSetVector<CTPPSDiamondDigi> > diamondDigis;
   event.getByToken( tokenDigi_, diamondDigis );
 
   edm::Handle< std::vector<TotemFEDInfo> > fedInfo;
   event.getByToken( tokenFEDInfo_, fedInfo );
 
   edm::Handle< edm::DetSetVector<CTPPSDiamondRecHit> > diamondRecHits;
   event.getByToken( tokenDiamondHit_, diamondRecHits );
 
   edm::Handle< edm::DetSetVector<CTPPSDiamondLocalTrack> > diamondLocalTracks;
   event.getByToken( tokenDiamondTrack_, diamondLocalTracks );
 
   // check validity
   bool valid = true;
   valid &= diamondVFATStatus.isValid();
   valid &= diamondDigis.isValid();
   valid &= fedInfo.isValid();
 
   if ( !valid ) {
     if ( verbosity_ ) {
       edm::LogProblem("CTPPSDiamondDQMSource")
         << "ERROR in TotemDQMModuleRP::analyze > some of the required inputs are not valid. Skipping this event.\n"
         << "    diamondVFATStatus.isValid = " << diamondVFATStatus.isValid() << "\n"
         << "    diamondDigis.isValid = " << diamondDigis.isValid() << "\n"
         << "    fedInfo.isValid = " << fedInfo.isValid();
     }
 
     return;
   }
 
   //------------------------------
   // RP Plots
   //------------------------------  
   
   //------------------------------
   // Correlation Plots
   //------------------------------
   
   for ( const auto& ds1 : *stripTracks ) {
     for ( const auto& tr1 : ds1 ) {
       if ( ! tr1.isValid() )  continue;
 
       CTPPSDetId rpId1( ds1.detId() );
       unsigned int arm1 = rpId1.arm();
       unsigned int stNum1 = rpId1.station();
       unsigned int rpNum1 = rpId1.rp();
       if (stNum1 != 0 || ( rpNum1 != 2 && rpNum1 != 3 ) )  continue;
       unsigned int idx1 = arm1*3 + rpNum1-2;
 
       for ( const auto& ds2 : *stripTracks ) {
         for ( const auto& tr2 : ds2 ) {
           if ( ! tr2.isValid() )  continue;
 
           CTPPSDetId rpId2(ds2.detId());
           unsigned int arm2 = rpId2.arm();
           unsigned int stNum2 = rpId2.station();
           unsigned int rpNum2 = rpId2.rp();
           if (stNum2 != 0 || ( rpNum2 != 2 && rpNum2 != 3 ) )  continue;
           unsigned int idx2 = arm2*3 + rpNum2-2;
 
           if ( idx1 >= idx2 ) globalPlot_.h_trackCorr_hor->Fill( 5-idx1, idx2 ); // strips-strips
         }
       }
       for ( const auto& ds2 : *diamondLocalTracks ) {
         for ( const auto& tr2 : ds2 ) {
           if ( ! tr2.isValid() ) continue;
           if ( tr2.getOOTIndex() != 1 ) continue;
           if ( excludeMultipleHits_ && tr2.getMultipleHits() > 0 ) continue;
 
           CTPPSDetId diamId2( ds2.detId() );
           unsigned int arm2 = diamId2.arm();
           if ( idx1 >= arm2*3+2 )
             globalPlot_.h_trackCorr_hor->Fill( 5-idx1, arm2*3+2 ); // strips-diamonds
           else
             globalPlot_.h_trackCorr_hor->Fill( 5-(arm2*3+2 ),idx1 ); // strips-diamonds
         }
       }
     }
   }
 
   for ( const auto& ds1 : *diamondLocalTracks ) {
     for ( const auto& tr1 : ds1 ) {
       if ( ! tr1.isValid() ) continue;
       if ( excludeMultipleHits_ && tr1.getMultipleHits() > 0 ) continue;
       if ( tr1.getOOTIndex() != 1 ) continue;
 
       CTPPSDetId diamId1( ds1.detId() );
       unsigned int arm1 = diamId1.arm();
 
       globalPlot_.h_trackCorr_hor->Fill( 5-(arm1*3+2), arm1*3+2 ); // diamonds-diamonds
 
       for ( const auto& ds2 : *diamondLocalTracks ) {
         for ( const auto& tr2 : ds2 ) {
           if ( ! tr2.isValid() ) continue;
           if ( excludeMultipleHits_ && tr2.getMultipleHits() > 0 ) continue;
           if ( tr2.getOOTIndex() != 1 ) continue;
 
           CTPPSDetId diamId2( ds2.detId() );
           unsigned int arm2 = diamId2.arm();
           if ( arm1 > arm2 ) globalPlot_.h_trackCorr_hor->Fill( 5-(arm1*3+2), arm2*3+2 ); // diamonds-diamonds
         }
       }
     }
   }
 
 
   // Using CTPPSDiamondDigi
   for ( const auto& digis : *diamondDigis ) {
     const CTPPSDiamondDetId detId( digis.detId() );
     CTPPSDiamondDetId detId_pot( digis.detId() );
 
     for ( const auto& digi : digis ) {
       detId_pot.setPlane( 0 );
       detId_pot.setChannel( 0 );
       if ( potPlots_.find( detId_pot ) == potPlots_.end() ) continue;
       //Leading without trailing investigation
       if      ( digi.getLeadingEdge() == 0 && digi.getTrailingEdge() == 0 ) potPlots_[detId_pot].leadingWithoutTrailingCumulativePot->Fill( 1 );
       else if ( digi.getLeadingEdge() != 0 && digi.getTrailingEdge() == 0 ) potPlots_[detId_pot].leadingWithoutTrailingCumulativePot->Fill( 2 );
       else if ( digi.getLeadingEdge() == 0 && digi.getTrailingEdge() != 0 ) potPlots_[detId_pot].leadingWithoutTrailingCumulativePot->Fill( 3 );
       else if ( digi.getLeadingEdge() != 0 && digi.getTrailingEdge() != 0 ) potPlots_[detId_pot].leadingWithoutTrailingCumulativePot->Fill( 4 );
 
       if ( digi.getLeadingEdge() != 0 ) {
         // FED BX monitoring (for MINIDAQ)
         for ( const auto& fit : *fedInfo ) {
           if ( ( detId.arm() == 1 && fit.getFEDId() == CTPPS_FED_ID_56 ) || ( detId.arm() == 0 && fit.getFEDId() == CTPPS_FED_ID_45 ) ) {
             potPlots_[detId_pot].activity_per_fedbx->Fill( fit.getBX() );
             potPlots_[detId_pot].activity_per_fedbx_short->Fill( fit.getBX() );
           }
         }
       }
 
       // HPTDC Errors
       const HPTDCErrorFlags hptdcErrors = digi.getHPTDCErrorFlags();
       for ( unsigned short hptdcErrorIndex = 1; hptdcErrorIndex < 16; ++hptdcErrorIndex )
         if ( hptdcErrors.getErrorId( hptdcErrorIndex-1 ) ) potPlots_[detId_pot].HPTDCErrorFlags_cumulative->Fill( hptdcErrorIndex );
       if ( digi.getMultipleHit() ) potPlots_[detId_pot].HPTDCErrorFlags_cumulative->Fill( 16 );
     }
   }
 
   // EC Errors
   for ( const auto& vfat_status : *diamondVFATStatus ) {
     const CTPPSDiamondDetId detId( vfat_status.detId() );
     CTPPSDiamondDetId detId_pot( vfat_status.detId() );
     detId_pot.setPlane( 0 );
     detId_pot.setChannel( 0 );
     for ( const auto& status : vfat_status ) {
       if ( !status.isOK() ) continue;
       if ( potPlots_.find(detId_pot) == potPlots_.end() ) continue;
 
       // Check Event Number
       for ( const auto& optorx : *fedInfo ) {
         if ( detId.arm() == 1 && optorx.getFEDId() == CTPPS_FED_ID_56 ) {
           potPlots_[detId_pot].ECCheck->Fill((int)((optorx.getLV1()& 0xFF)-((unsigned int) status.getEC() & 0xFF)) & 0xFF);
           if ( ( static_cast<int>( ( optorx.getLV1() & 0xFF )-status.getEC() ) != EC_difference_56_ ) && ( static_cast<uint8_t>( ( optorx.getLV1() & 0xFF )-status.getEC() ) < 128 ) )
             EC_difference_56_ = static_cast<int>( optorx.getLV1() & 0xFF )-( static_cast<unsigned int>( status.getEC() ) & 0xFF );
           if ( EC_difference_56_ != 1 && EC_difference_56_ != -500 && EC_difference_56_ < 128 && EC_difference_56_ > -128 )
             if (verbosity_) edm::LogProblem("CTPPSDiamondDQMSource")  << "FED " << CTPPS_FED_ID_56 << ": ECError at EV: 0x"<< std::hex << optorx.getLV1()
                                                                       << "\t\tVFAT EC: 0x"<< static_cast<unsigned int>( status.getEC() )
                                                                       << "\twith ID: " << std::dec << detId
                                                                       << "\tdiff: " <<  EC_difference_56_;
         }
         else if ( detId.arm() == 0 && optorx.getFEDId()== CTPPS_FED_ID_45 ) {
           potPlots_[detId_pot].ECCheck->Fill((int)((optorx.getLV1()& 0xFF)-status.getEC()) & 0xFF);
           if ( ( static_cast<int>( ( optorx.getLV1() & 0xFF )-status.getEC() ) != EC_difference_45_ ) && ( static_cast<uint8_t>( ( optorx.getLV1() & 0xFF )-status.getEC() ) < 128 ) )
             EC_difference_45_ = static_cast<int>( optorx.getLV1() & 0xFF )-( static_cast<unsigned int>( status.getEC() ) & 0xFF );
           if ( EC_difference_45_ != 1 && EC_difference_45_ != -500 && EC_difference_45_ < 128 && EC_difference_45_ > -128 )
             if (verbosity_) edm::LogProblem("CTPPSDiamondDQMSource")  << "FED " << CTPPS_FED_ID_45 << ": ECError at EV: 0x"<< std::hex << optorx.getLV1()
                                                                       << "\t\tVFAT EC: 0x"<< static_cast<unsigned int>( status.getEC() )
                                                                       << "\twith ID: " << std::dec << detId
                                                                       << "\tdiff: " <<  EC_difference_45_;
         }
       }
     }
   }
 
   // Using CTPPSDiamondRecHit
   std::unordered_map<unsigned int, std::set<unsigned int> > planes;
 
   for ( const auto& rechits : *diamondRecHits ) {
     CTPPSDiamondDetId detId_pot( rechits.detId() );
     detId_pot.setPlane( 0 );
     detId_pot.setChannel( 0 );
     const CTPPSDiamondDetId detId( rechits.detId() );
 
     for ( const auto& rechit : rechits ) {
       if ( excludeMultipleHits_ && rechit.getMultipleHits() > 0 ) continue;
       planes[detId_pot].insert( detId.plane() );
 
       if ( potPlots_.find( detId_pot ) == potPlots_.end() ) continue;
 
       float UFSDShift = 0.0;
       if ( rechit.getYWidth() < 3 ) UFSDShift = 0.5;  // Display trick for UFSD that have 2 pixels with same X
         
       TH2F *hitHistoTmp = potPlots_[detId_pot].hitDistribution2d->getTH2F();
       TAxis *hitHistoTmpYAxis = hitHistoTmp->GetYaxis();
       int startBin = hitHistoTmpYAxis->FindBin( rechit.getX() - 0.5*rechit.getXWidth() );
       int numOfBins = rechit.getXWidth()/DISPLAY_RESOLUTION_FOR_HITS_MM;
       for ( int i=0; i<numOfBins; ++i) {
         hitHistoTmp->Fill( detId.plane(), hitHistoTmpYAxis->GetBinCenter(startBin+i) + UFSDShift );
       }
       
       TH2F *hitHistoOOTTmp = potPlots_[detId_pot].hitDistribution2dOOT->getTH2F();
       TAxis *hitHistoOOTTmpYAxis = hitHistoOOTTmp->GetYaxis();
       startBin = hitHistoOOTTmpYAxis->FindBin( rechit.getX() - 0.5*rechit.getXWidth() );
       numOfBins = rechit.getXWidth()/DISPLAY_RESOLUTION_FOR_HITS_MM;
       for ( int i=0; i<numOfBins; ++i) {
         hitHistoOOTTmp->Fill( detId.plane() + 0.25 * rechit.getOOTIndex(), hitHistoOOTTmpYAxis->GetBinCenter(startBin+i) );
       }
 
       potPlots_[detId_pot].leadingEdgeCumulativePot->Fill( rechit.getT() );
       potPlots_[detId_pot].timeOverThresholdCumulativePot->Fill( rechit.getToT() );
       potPlots_[detId_pot].leadingTrailingCorrelationPot->Fill( rechit.getT(), rechit.getT() + rechit.getToT() );
 
       switch ( rechit.getOOTIndex() ) {
         case 0: {
           potPlots_[detId_pot].activity_per_bx_minus1->Fill( event.bunchCrossing() );
           potPlots_[detId_pot].activity_per_bx_short_minus1->Fill( event.bunchCrossing() );
         } break;
         case 1: {
           potPlots_[detId_pot].activity_per_bx->Fill( event.bunchCrossing() );
           potPlots_[detId_pot].activity_per_bx_short->Fill( event.bunchCrossing() );
         } break;
         case 2: {
           potPlots_[detId_pot].activity_per_bx_plus1->Fill( event.bunchCrossing() );
           potPlots_[detId_pot].activity_per_bx_short_plus1->Fill( event.bunchCrossing() );
         } break;
       }      
     }
   }
   
   for ( const auto& plt : potPlots_ ) {
     plt.second.activePlanes->Fill( planes[plt.first].size() );
   }
   
   // Using CTPPSDiamondLocalTrack
   for ( const auto& tracks : *diamondLocalTracks ) {
     CTPPSDiamondDetId detId_pot( tracks.detId() );
     detId_pot.setPlane( 0 );
     detId_pot.setChannel( 0 );
     const CTPPSDiamondDetId detId( tracks.detId() );
 
     for ( const auto& track : tracks ) {
       if ( ! track.isValid() ) continue;
       if ( excludeMultipleHits_ && track.getMultipleHits() > 0 ) continue;
       if ( potPlots_.find( detId_pot ) == potPlots_.end() ) continue;
 
       TH2F *trackHistoOOTTmp = potPlots_[detId_pot].trackDistributionOOT->getTH2F();
       TAxis *trackHistoOOTTmpYAxis = trackHistoOOTTmp->GetYaxis();
       int startBin = trackHistoOOTTmpYAxis->FindBin( track.getX0() - track.getX0Sigma() );
       int numOfBins = 2*track.getX0Sigma()/DISPLAY_RESOLUTION_FOR_HITS_MM;
       for ( int i=0; i<numOfBins; ++i) {
         trackHistoOOTTmp->Fill( track.getOOTIndex(), trackHistoOOTTmpYAxis->GetBinCenter(startBin+i) );
       }
       
       if ( track.getOOTIndex() == 1 ) {
         TH1F *trackHistoInTimeTmp = potPlots_[detId_pot].trackDistribution->getTH1F();
         int startBin = trackHistoInTimeTmp->FindBin( track.getX0() - track.getX0Sigma() );
         int numOfBins = 2*track.getX0Sigma()/DISPLAY_RESOLUTION_FOR_HITS_MM;
         for ( int i=0; i<numOfBins; ++i) {
           trackHistoInTimeTmp->Fill( trackHistoInTimeTmp->GetBinCenter(startBin+i) );
         }
       }
     }
   }
   
   // Tomography of diamonds using strips
   for ( const auto& rechits : *diamondRecHits ) {
     CTPPSDiamondDetId detId_pot( rechits.detId() );
     detId_pot.setPlane( 0 );
     detId_pot.setChannel( 0 );
     const CTPPSDiamondDetId detId( rechits.detId() );
 
     for ( const auto& rechit : rechits ) {
       if ( excludeMultipleHits_ && rechit.getMultipleHits() > 0 ) continue;
       if ( !stripTracks.isValid() ) continue;
       if ( potPlots_.find( detId_pot ) == potPlots_.end() ) continue;
 
       for ( const auto& ds : *stripTracks ) {
         const CTPPSDetId stripId( ds.detId() );
         for ( const auto& striplt : ds ) {
           if ( !striplt.isValid() ) continue;
           if ( stripId.arm() != detId_pot.arm() ) continue;
           if ( striplt.getTx() > minimumStripAngleForTomography_ || striplt.getTy() > minimumStripAngleForTomography_) continue; 
           if ( stripId.rp() == CTPPS_FAR_RP_ID ) {
             switch ( rechit.getOOTIndex() ) {
               case 0: {
                 potPlots_[detId_pot].stripTomographyAllFar_minus1->Fill( striplt.getX0() + 50*detId.plane(), striplt.getY0() );
               } break;
               case 1: {
                 potPlots_[detId_pot].stripTomographyAllFar->Fill( striplt.getX0() + 50*detId.plane(), striplt.getY0() );
               } break;
               case 2: {
                 potPlots_[detId_pot].stripTomographyAllFar_plus1->Fill( striplt.getX0() + 50*detId.plane(), striplt.getY0() );
               } break;
             }
           }
           else if ( stripId.rp() == CTPPS_NEAR_RP_ID ) {
             switch ( rechit.getOOTIndex() ) {
               case 0: {
                 potPlots_[detId_pot].stripTomographyAllNear_minus1->Fill( striplt.getX0() + 50*detId.plane(), striplt.getY0() );
               } break;
               case 1: {
                 potPlots_[detId_pot].stripTomographyAllNear->Fill( striplt.getX0() + 50*detId.plane(), striplt.getY0() );
               } break;
               case 2: {
                 potPlots_[detId_pot].stripTomographyAllNear_plus1->Fill( striplt.getX0() + 50*detId.plane(), striplt.getY0() );
               } break;
             }
           }
         }
       }
     }
   }
 
   //------------------------------
   // Clock Plots
   //------------------------------
 
   for ( const auto& digis : *diamondDigis ) {
     const CTPPSDiamondDetId detId( digis.detId() );
     CTPPSDiamondDetId detId_pot( digis.detId() );
     if ( detId.channel() == CHANNEL_OF_VFAT_CLOCK ) {
       detId_pot.setPlane( 0 );
       detId_pot.setChannel( 0 );
       for ( const auto& digi : digis ) {
         if ( detId.plane() == 1 ) {
           potPlots_[detId_pot].clock_Digi1_le->Fill( HPTDC_BIN_WIDTH_NS * digi.getLeadingEdge() );
           potPlots_[detId_pot].clock_Digi1_te->Fill( HPTDC_BIN_WIDTH_NS * digi.getTrailingEdge() );
         }
         if ( detId.plane() == 3 ) {
           potPlots_[detId_pot].clock_Digi3_le->Fill( HPTDC_BIN_WIDTH_NS * digi.getLeadingEdge() );
           potPlots_[detId_pot].clock_Digi3_te->Fill( HPTDC_BIN_WIDTH_NS * digi.getTrailingEdge() );
         }
       }
     }
   }
 
   //------------------------------
   // Plane Plots
   //------------------------------
 
   // Using CTPPSDiamondDigi
   std::unordered_map<unsigned int, unsigned int> channelsPerPlane;
   for ( const auto& digis : *diamondDigis ) {
     const CTPPSDiamondDetId detId( digis.detId() );
     CTPPSDiamondDetId detId_plane( digis.detId() );
     for ( const auto& digi : digis ) {
       detId_plane.setChannel( 0 );
       if ( planePlots_.find( detId_plane ) == planePlots_.end() ) continue;
 
       planePlots_[detId_plane].threshold_voltage->Fill( detId.channel(), digi.getThresholdVoltage() );
 
       if ( digi.getLeadingEdge() != 0 ) {
         planePlots_[detId_plane].digiProfileCumulativePerPlane->Fill( detId.channel() );
         if ( channelsPerPlane.find(detId_plane) != channelsPerPlane.end() ) channelsPerPlane[detId_plane]++;
         else channelsPerPlane[detId_plane] = 0;
       }
     }
   }
 
   for ( const auto& plt : channelsPerPlane ) {
     planePlots_[plt.first].hit_multiplicity->Fill( plt.second );
   }
   
   // Using CTPPSDiamondRecHit
   for ( const auto& rechits : *diamondRecHits ) {
     CTPPSDiamondDetId detId_plane( rechits.detId() );
     detId_plane.setChannel( 0 );
     for ( const auto& rechit : rechits ) {
       if ( excludeMultipleHits_ && rechit.getMultipleHits() > 0 ) continue;
       if ( planePlots_.find( detId_plane ) != planePlots_.end() ) {
         TH1F *hitHistoTmp = planePlots_[detId_plane].hitProfile->getTH1F();
         int startBin = hitHistoTmp->FindBin( rechit.getX() - 0.5*rechit.getXWidth() );
         int numOfBins = rechit.getXWidth()/DISPLAY_RESOLUTION_FOR_HITS_MM;
         for ( int i=0; i<numOfBins; ++i) {
           if ( rechit.getOOTIndex() == 1 ) hitHistoTmp->Fill( hitHistoTmp->GetBinCenter(startBin+i) );
         }
       }
     }
   }
 
   // Tomography of diamonds using strips
   for ( const auto& rechits : *diamondRecHits ) {
     CTPPSDiamondDetId detId_plane( rechits.detId() );
     detId_plane.setChannel( 0 );
     for ( const auto& rechit : rechits ) {
       if ( excludeMultipleHits_ && rechit.getMultipleHits() > 0 ) continue;
       if ( !stripTracks.isValid() ) continue;
       if (planePlots_.find(detId_plane) == planePlots_.end()) continue;
       
       for ( const auto& ds : *stripTracks ) {
         const CTPPSDetId stripId(ds.detId());
         for ( const auto& striplt : ds ) {
           if (! striplt.isValid()) continue;
           if ( stripId.arm() != detId_plane.arm() ) continue;
           if ( striplt.getTx() > minimumStripAngleForTomography_ || striplt.getTy() > minimumStripAngleForTomography_ ) continue;
           if ( stripId.rp() == CTPPS_FAR_RP_ID ) {
             planePlots_[detId_plane].stripTomography_far->Fill( striplt.getX0(), striplt.getY0() + 50*( rechit.getOOTIndex()-1 ) );
           }
           else if ( stripId.rp() == CTPPS_NEAR_RP_ID ) {
             planePlots_[detId_plane].stripTomography_near->Fill( striplt.getX0(), striplt.getY0() + 50*( rechit.getOOTIndex()-1 ) );
           }
         }
       }
     }
   }
   //------------------------------
   // Channel Plots
   //------------------------------
 
   //Check Event Number
   for ( const auto& vfat_status : *diamondVFATStatus ) {
     const CTPPSDiamondDetId detId( vfat_status.detId() );
     for ( const auto& status : vfat_status ) {
       if ( !status.isOK() ) continue;
       if ( channelPlots_.find(detId) != channelPlots_.end() ) {
         for ( const auto& optorx : *fedInfo ) {
           if ( ( detId.arm() == 1 && optorx.getFEDId() == CTPPS_FED_ID_56 ) || ( detId.arm() == 0 && optorx.getFEDId() == CTPPS_FED_ID_45 ) ) {
             channelPlots_[detId].ECCheckPerChannel->Fill((int)((optorx.getLV1()& 0xFF)-((unsigned int) status.getEC() & 0xFF)) & 0xFF);
           }
         }
       }
     }
   }
 
   // digi profile cumulative
   for ( const auto& digis : *diamondDigis ) {
     const CTPPSDiamondDetId detId( digis.detId() );
     for ( const auto& digi : digis ) {
       if ( channelPlots_.find( detId ) != channelPlots_.end() ) {
         // HPTDC Errors
         const HPTDCErrorFlags hptdcErrors = digi.getHPTDCErrorFlags();
         for ( unsigned short hptdcErrorIndex = 1; hptdcErrorIndex < 16; ++hptdcErrorIndex )
           if ( hptdcErrors.getErrorId( hptdcErrorIndex-1 ) ) channelPlots_[detId].HPTDCErrorFlags->Fill( hptdcErrorIndex );
         if ( digi.getMultipleHit() ) channelPlots_[detId].HPTDCErrorFlags->Fill( 16 );
 
         // Check dropped trailing edges
         if      ( digi.getLeadingEdge() == 0 && digi.getTrailingEdge() == 0 ) channelPlots_[detId].leadingWithoutTrailing->Fill( 1 );
         else if ( digi.getLeadingEdge() != 0 && digi.getTrailingEdge() == 0 ) channelPlots_[detId].leadingWithoutTrailing->Fill( 2 );
         else if ( digi.getLeadingEdge() == 0 && digi.getTrailingEdge() != 0 ) channelPlots_[detId].leadingWithoutTrailing->Fill( 3 );
         else if ( digi.getLeadingEdge() != 0 && digi.getTrailingEdge() != 0 ) channelPlots_[detId].leadingWithoutTrailing->Fill( 4 );
       }
     }
   }
   
   // Using CTPPSDiamondRecHit
   for ( const auto& rechits : *diamondRecHits ) {
     CTPPSDiamondDetId detId( rechits.detId() );
     for ( const auto& rechit : rechits ) {
       if ( excludeMultipleHits_ && rechit.getMultipleHits() > 0 ) continue;
       if ( channelPlots_.find( detId ) != channelPlots_.end() ) {
         channelPlots_[detId].LeadingEdgeCumulativePerChannel->Fill( rechit.getT() );
         channelPlots_[detId].TimeOverThresholdCumulativePerChannel->Fill( rechit.getToT() );
         channelPlots_[detId].LeadingTrailingCorrelationPerChannel->Fill( rechit.getT(), rechit.getT() + rechit.getToT() );
         ++(channelPlots_[detId].hitsCounterPerLumisection);
       }
     }
   }
   
   // Tomography of diamonds using strips
   for ( const auto& rechits : *diamondRecHits ) {
     const CTPPSDiamondDetId detId( rechits.detId() );
     for ( const auto& rechit : rechits ) {
       if ( excludeMultipleHits_ && rechit.getMultipleHits() > 0 ) continue;
       if ( stripTracks.isValid() ) {
         if (channelPlots_.find(detId) == channelPlots_.end()) continue;
         for ( const auto& ds : *stripTracks ) {
           for ( const auto& striplt : ds ) {
             CTPPSDetId stripId(ds.detId());
             if ( !striplt.isValid() ) continue;
             if ( stripId.arm() != detId.arm() ) continue;
             if ( striplt.getTx() > minimumStripAngleForTomography_ || striplt.getTy() > minimumStripAngleForTomography_ ) continue;
             if ( stripId.rp() == CTPPS_FAR_RP_ID ) {
               channelPlots_[detId].stripTomography_far->Fill( striplt.getX0(), striplt.getY0() + 50*( rechit.getOOTIndex()-1 ) );
             }
             else if ( stripId.rp() == CTPPS_NEAR_RP_ID ) {
               channelPlots_[detId].stripTomography_near->Fill( striplt.getX0(), striplt.getY0() + 50*( rechit.getOOTIndex()-1 ) );
             }
           }
         }
       }
     }
   }
   
 }

void CTPPSDiamondDQMSource::beginLuminosityBlock	(	const edm::LuminosityBlock &	,
		const edm::EventSetup &
	)

protected

Definition at line 371 of file CTPPSDiamondDQMSource.cc.

References channelPlots_, plotFactory::plot, and SEC_PER_LUMI_SECTION.

 {
   for ( auto& plot : channelPlots_ ) {
     if ( plot.second.hitsCounterPerLumisection != 0 ) {
       plot.second.hit_rate->Fill( (double) plot.second.hitsCounterPerLumisection / SEC_PER_LUMI_SECTION );
     }
     plot.second.hitsCounterPerLumisection = 0;
   }
 }

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

overrideprotected

Definition at line 346 of file CTPPSDiamondDQMSource.cc.

References DQMStore::IBooker::cd(), channelPlots_, CTPPS_DIAMOND_NUM_OF_CHANNELS, CTPPS_DIAMOND_NUM_OF_PLANES, CTPPS_DIAMOND_RP_ID, CTPPS_DIAMOND_STATION_ID, CTPPS_NUM_OF_ARMS, globalPlot_, planePlots_, potPlots_, and DQMStore::IBooker::setCurrentFolder().

 {
   ibooker.cd();
   ibooker.setCurrentFolder( "CTPPS" );
   
   globalPlot_= GlobalPlots( ibooker );
 
   for ( unsigned short arm = 0; arm < CTPPS_NUM_OF_ARMS; ++arm ) {
     const CTPPSDiamondDetId rpId( arm, CTPPS_DIAMOND_STATION_ID, CTPPS_DIAMOND_RP_ID );
     potPlots_[rpId] = PotPlots( ibooker, rpId );
     for ( unsigned short pl = 0; pl < CTPPS_DIAMOND_NUM_OF_PLANES; ++pl ) {
       const CTPPSDiamondDetId plId( arm, CTPPS_DIAMOND_STATION_ID, CTPPS_DIAMOND_RP_ID, pl );
       planePlots_[plId] = PlanePlots( ibooker, plId);
       for ( unsigned short ch = 0; ch < CTPPS_DIAMOND_NUM_OF_CHANNELS; ++ch ) {
         const CTPPSDiamondDetId chId( arm, CTPPS_DIAMOND_STATION_ID, CTPPS_DIAMOND_RP_ID, pl, ch );
         channelPlots_[chId] = ChannelPlots( ibooker, chId );
       }  
     }
   }
 }

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

overrideprotected

Definition at line 340 of file CTPPSDiamondDQMSource.cc.

341 {}

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

protected

Definition at line 882 of file CTPPSDiamondDQMSource.cc.

883 {}

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

protected

Definition at line 888 of file CTPPSDiamondDQMSource.cc.

References DEFINE_FWK_MODULE.

889 {}

Member Data Documentation

const int CTPPSDiamondDQMSource::CHANNEL_OF_VFAT_CLOCK = 30

staticprivate

Definition at line 54 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 163 of file CTPPSDiamondDQMSource.cc.

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

const int CTPPSDiamondDQMSource::CTPPS_DIAMOND_NUM_OF_CHANNELS = 12

staticprivate

Definition at line 63 of file CTPPSDiamondDQMSource.cc.

Referenced by bookHistograms().

const int CTPPSDiamondDQMSource::CTPPS_DIAMOND_NUM_OF_PLANES = 4

staticprivate

Definition at line 62 of file CTPPSDiamondDQMSource.cc.

Referenced by bookHistograms().

const int CTPPSDiamondDQMSource::CTPPS_DIAMOND_RP_ID = 6

staticprivate

Definition at line 59 of file CTPPSDiamondDQMSource.cc.

Referenced by bookHistograms().

const int CTPPSDiamondDQMSource::CTPPS_DIAMOND_STATION_ID = 1

staticprivate

Definition at line 58 of file CTPPSDiamondDQMSource.cc.

Referenced by bookHistograms().

const int CTPPSDiamondDQMSource::CTPPS_FAR_RP_ID = 3

staticprivate

Definition at line 61 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

const int CTPPSDiamondDQMSource::CTPPS_FED_ID_45 = 583

staticprivate

Definition at line 64 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

const int CTPPSDiamondDQMSource::CTPPS_FED_ID_56 = 582

staticprivate

Definition at line 65 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

const int CTPPSDiamondDQMSource::CTPPS_NEAR_RP_ID = 2

staticprivate

Definition at line 60 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

const int CTPPSDiamondDQMSource::CTPPS_NUM_OF_ARMS = 2

staticprivate

Definition at line 57 of file CTPPSDiamondDQMSource.cc.

Referenced by bookHistograms().

const double CTPPSDiamondDQMSource::DISPLAY_RESOLUTION_FOR_HITS_MM = 0.1

staticprivate

Definition at line 55 of file CTPPSDiamondDQMSource.cc.

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

int CTPPSDiamondDQMSource::EC_difference_45_

private

Definition at line 125 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

int CTPPSDiamondDQMSource::EC_difference_56_

private

Definition at line 125 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

bool CTPPSDiamondDQMSource::excludeMultipleHits_

private

Definition at line 74 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

GlobalPlots CTPPSDiamondDQMSource::globalPlot_

private

Definition at line 87 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze(), and bookHistograms().

const double CTPPSDiamondDQMSource::HPTDC_BIN_WIDTH_NS = 25./1024

staticprivate

Definition at line 56 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

double CTPPSDiamondDQMSource::minimumStripAngleForTomography_

private

Definition at line 75 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 142 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 124 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze(), and bookHistograms().

const double CTPPSDiamondDQMSource::SEC_PER_LUMI_SECTION = 23.31

staticprivate

Definition at line 53 of file CTPPSDiamondDQMSource.cc.

Referenced by beginLuminosityBlock().

edm::EDGetTokenT< edm::DetSetVector<CTPPSDiamondRecHit> > CTPPSDiamondDQMSource::tokenDiamondHit_

private

Definition at line 70 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

edm::EDGetTokenT< edm::DetSetVector<CTPPSDiamondLocalTrack> > CTPPSDiamondDQMSource::tokenDiamondTrack_

private

Definition at line 71 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

edm::EDGetTokenT< edm::DetSetVector<CTPPSDiamondDigi> > CTPPSDiamondDQMSource::tokenDigi_

private

Definition at line 69 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 72 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 68 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

edm::EDGetTokenT< edm::DetSetVector<TotemVFATStatus> > CTPPSDiamondDQMSource::tokenStatus_

private

Definition at line 67 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

unsigned int CTPPSDiamondDQMSource::verbosity_

private

Definition at line 76 of file CTPPSDiamondDQMSource.cc.

Referenced by analyze().

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