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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () 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

ESProxyIndex const *	esGetTokenIndices (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::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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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)

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

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

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)

Detailed Description

Definition at line 32 of file CTPPSTrackDistributionPlotter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 179 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 ( )

inlineoverride

Definition at line 36 of file CTPPSTrackDistributionPlotter.cc.

References analyze(), and endJob().

36 {}

Member Function Documentation

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

overrideprivate

Definition at line 193 of file CTPPSTrackDistributionPlotter.cc.

References protons_cff::arm, armPlots, events_per_arm_, events_total_, edm::Event::getByToken(), edm::EventBase::id(), edm::EventID::luminosityBlock(), year_2016_postTS2_cff::rpId, rpPlots, edm::EventID::run(), CTPPSDetId::sdTrackingPixel, l1t::tracks, tracksToken_, and x_pitch_pixels_.

Referenced by ~CTPPSTrackDistributionPlotter().

                                                                                         {
   // 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.getRPId());
     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.fill(trk.getX(), trk.getY(), trk.getTime());
 
     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.getRPId());
 
     for (const auto& t2 : *tracks) {
       const CTPPSDetId rpId2(t2.getRPId());
 
       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.getX(), t1.getY(), t2.getX(), t2.getY());
     }
   }
 
   // update counters
   events_total_++;
 
   if (m_mult[0] > 0)
     events_per_arm_[0]++;
   if (m_mult[1] > 0)
     events_per_arm_[1]++;
 }

void CTPPSTrackDistributionPlotter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 251 of file CTPPSTrackDistributionPlotter.cc.

References armPlots, DEFINE_FWK_MODULE, events_per_arm_, events_total_, outputFile_, and rpPlots.

Referenced by ~CTPPSTrackDistributionPlotter().

                                            {
   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();
   }
 }