CTPPSTrackDistributionPlotter Class Reference

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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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
bool	registeredToConsumeMany (TypeID const &, 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::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)
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 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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
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 ( 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 ( )

inlineoverride

Definition at line 38 of file CTPPSTrackDistributionPlotter.cc.

{}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 260 of file CTPPSTrackDistributionPlotter.cc.

References CTPPSDetId::arm(), armPlots, events_per_arm_, events_total_, edm::Event::getByToken(), mps_fire::i, edm::EventBase::id(), dqmiolumiharvest::j, edm::EventID::luminosityBlock(), m_mult(), rpPlots, edm::EventID::run(), CTPPSDetId::sdTrackingPixel, 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]++;
}

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

Member Data Documentation

std::map<unsigned int, ArmPlots> CTPPSTrackDistributionPlotter::armPlots

private

Definition at line 241 of file CTPPSTrackDistributionPlotter.cc.

Referenced by analyze(), CTPPSTrackDistributionPlotter(), and endJob().

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

private

Definition at line 51 of file CTPPSTrackDistributionPlotter.cc.

Referenced by analyze(), and endJob().

unsigned int CTPPSTrackDistributionPlotter::events_total_

private

Definition at line 50 of file CTPPSTrackDistributionPlotter.cc.

Referenced by analyze(), and endJob().

std::string CTPPSTrackDistributionPlotter::outputFile_

private

Definition at line 48 of file CTPPSTrackDistributionPlotter.cc.

Referenced by endJob().

std::map<unsigned int, RPPlots> CTPPSTrackDistributionPlotter::rpPlots

private

Definition at line 118 of file CTPPSTrackDistributionPlotter.cc.

Referenced by analyze(), and endJob().

edm::EDGetTokenT<CTPPSLocalTrackLiteCollection> CTPPSTrackDistributionPlotter::tracksToken_

private

Definition at line 44 of file CTPPSTrackDistributionPlotter.cc.

Referenced by analyze().

double CTPPSTrackDistributionPlotter::x_pitch_pixels_

private

Definition at line 46 of file CTPPSTrackDistributionPlotter.cc.

Referenced by analyze().