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 final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () const

	~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, pwdgSkimBPark_cfi::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_, 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();
   }
 }