Inheritance diagram for CTPPSTrackDistributionPlotter:

Classes
struct	ArmPlots

struct	RPPlots

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

	~CTPPSTrackDistributionPlotter () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () 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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	endJob () override

Private Attributes
std::map< unsigned int, ArmPlots >	armPlots

std::map< unsigned int, unsigned int >	events_per_arm_

unsigned int	events_total_

std::string	outputFile_

std::map< unsigned int, RPPlots >	rpPlots

edm::EDGetTokenT< CTPPSLocalTrackLiteCollection >	tracksToken_

double	x_pitch_pixels_

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 34 of file CTPPSTrackDistributionPlotter.cc.

Constructor & Destructor Documentation

◆ CTPPSTrackDistributionPlotter()

CTPPSTrackDistributionPlotter::CTPPSTrackDistributionPlotter ( const edm::ParameterSet & iConfig )

explicit

Definition at line 246 of file CTPPSTrackDistributionPlotter.cc.

References armPlots, and edm::ParameterSet::getParameter().

     : tracksToken_(consumes<CTPPSLocalTrackLiteCollection>(iConfig.getParameter<edm::InputTag>("tagTracks"))),
       x_pitch_pixels_(iConfig.getUntrackedParameter<double>("x_pitch_pixels", 150E-3)),
       outputFile_(iConfig.getParameter<std::string>("outputFile")),
       events_total_(0) {
   armPlots[0].rpId_N = iConfig.getParameter<unsigned int>("rpId_45_N");
   armPlots[0].rpId_F = iConfig.getParameter<unsigned int>("rpId_45_F");
 
   armPlots[1].rpId_N = iConfig.getParameter<unsigned int>("rpId_56_N");
   armPlots[1].rpId_F = iConfig.getParameter<unsigned int>("rpId_56_F");
 }

◆ ~CTPPSTrackDistributionPlotter()

CTPPSTrackDistributionPlotter::~CTPPSTrackDistributionPlotter ( )

inlineoverride

Definition at line 38 of file CTPPSTrackDistributionPlotter.cc.

38 {}

Member Function Documentation

◆ analyze()

void CTPPSTrackDistributionPlotter::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 260 of file CTPPSTrackDistributionPlotter.cc.

References protons_cff::arm, CTPPSDetId::arm(), armPlots, events_per_arm_, events_total_, mps_fire::i, iEvent, dqmiolumiharvest::j, rpPlots, CTPPSDetId::sdTrackingPixel, RandomServiceHelper::t1, RandomServiceHelper::t2, DiMuonV_cfg::tracks, tracksToken_, and x_pitch_pixels_.

                                                                                         {
   // get input
   edm::Handle<CTPPSLocalTrackLiteCollection> tracks;
   iEvent.getByToken(tracksToken_, tracks);
 
   // process tracks
   std::map<unsigned int, unsigned int> m_mult;
 
   for (const auto& trk : *tracks) {
     CTPPSDetId rpId(trk.rpId());
     unsigned int rpDecId = rpId.arm() * 100 + rpId.station() * 10 + rpId.rp();
     bool rpPixel = (rpId.subdetId() == CTPPSDetId::sdTrackingPixel);
 
     auto& pl = rpPlots[rpDecId];
     if (!pl.initialized)
       pl.init(rpPixel, x_pitch_pixels_);
 
     pl.fillOneTrack(trk.x(), trk.y(), trk.time());
 
     m_mult[rpId.arm()]++;
   }
 
   for (unsigned int arm = 0; arm < 2; ++arm) {
     auto& st = armPlots[arm].m_stat[iEvent.id().run()][iEvent.id().luminosityBlock()];
     st.s1++;
     st.sN += m_mult[arm];
   }
 
   for (const auto& t1 : *tracks) {
     const CTPPSDetId rpId1(t1.rpId());
 
     for (const auto& t2 : *tracks) {
       const CTPPSDetId rpId2(t2.rpId());
 
       if (rpId1.arm() != rpId2.arm())
         continue;
 
       auto& ap = armPlots[rpId1.arm()];
 
       const unsigned int rpDecId1 = rpId1.arm() * 100 + rpId1.station() * 10 + rpId1.rp();
       const unsigned int rpDecId2 = rpId2.arm() * 100 + rpId2.station() * 10 + rpId2.rp();
 
       if (rpDecId1 == ap.rpId_N && rpDecId2 == ap.rpId_F)
         ap.fill(t1.x(), t1.y(), t2.x(), t2.y());
     }
   }
 
   for (unsigned int i = 0; i < tracks->size(); ++i) {
     for (unsigned int j = 0; j < tracks->size(); ++j) {
       if (i == j)
         continue;
 
       const auto& tr_i = tracks->at(i);
       const auto& tr_j = tracks->at(j);
 
       if (tr_i.rpId() != tr_j.rpId())
         continue;
 
       CTPPSDetId rpId(tr_i.rpId());
       unsigned int rpDecId = rpId.arm() * 100 + rpId.station() * 10 + rpId.rp();
 
       auto& pl = rpPlots[rpDecId];
 
       pl.h2_de_x_vs_x->Fill(tr_i.x(), tr_j.x() - tr_i.x());
       pl.h2_de_x_vs_y->Fill(tr_i.y(), tr_j.x() - tr_i.x());
       pl.h2_de_y_vs_x->Fill(tr_i.x(), tr_j.y() - tr_i.y());
       pl.h2_de_y_vs_y->Fill(tr_i.y(), tr_j.y() - tr_i.y());
     }
   }
 
   // update counters
   events_total_++;
 
   if (m_mult[0] > 0)
     events_per_arm_[0]++;
   if (m_mult[1] > 0)
     events_per_arm_[1]++;
 }

◆ endJob()

void CTPPSTrackDistributionPlotter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 341 of file CTPPSTrackDistributionPlotter.cc.

References armPlots, events_per_arm_, events_total_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, outputFile_, and rpPlots.

                                            {
   edm::LogInfo("CTPPSTrackDistributionPlotter")
       << "    events processed: " << events_total_ << " (" << std::scientific << double(events_total_) << ")\n"
       << "    events with tracks in sector 45: " << events_per_arm_[0] << " (" << double(events_per_arm_[0]) << ")\n"
       << "    events with tracks in sector 56: " << events_per_arm_[1] << " (" << double(events_per_arm_[1]) << ")";
 
   auto f_out = std::make_unique<TFile>(outputFile_.c_str(), "recreate");
 
   for (const auto& it : rpPlots) {
     gDirectory = f_out->mkdir(Form("RP %u", it.first));
     it.second.write();
   }
 
   for (const auto& it : armPlots) {
     gDirectory = f_out->mkdir(Form("arm %u", it.first));
     it.second.write();
   }
 }