CTPPSAcceptancePlotter Class Reference

Inheritance diagram for CTPPSAcceptancePlotter:

Classes
struct	DoubleArmPlots
struct	SingleArmPlots

Public Member Functions
	CTPPSAcceptancePlotter (const edm::ParameterSet &)
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	beginJob () override
void	endJob () override

Private Attributes
std::vector< std::set< unsigned int > >	doubleArmConfigurations
std::map< std::set< unsigned int >, DoubleArmPlots >	doubleArmPlots
std::string	outputFile_
unsigned int	rpId_45_F_
unsigned int	rpId_45_N_
unsigned int	rpId_56_F_
unsigned int	rpId_56_N_
std::vector< std::set< unsigned int > >	singleArmConfigurations
std::map< std::set< unsigned int >, SingleArmPlots >	singleArmPlots
edm::EDGetTokenT< edm::HepMCProduct >	tokenHepMC_
edm::EDGetTokenT< CTPPSLocalTrackLiteCollection >	tokenTracks_

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 28 of file CTPPSAcceptancePlotter.cc.

Constructor & Destructor Documentation

CTPPSAcceptancePlotter()

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

explicit

Definition at line 132 of file CTPPSAcceptancePlotter.cc.

References doubleArmConfigurations, rpId_45_F_, rpId_45_N_, rpId_56_F_, rpId_56_N_, and singleArmConfigurations.

    : tokenHepMC_(consumes<edm::HepMCProduct>(iConfig.getParameter<edm::InputTag>("tagHepMC"))),
      tokenTracks_(consumes<CTPPSLocalTrackLiteCollection>(iConfig.getParameter<edm::InputTag>("tagTracks"))),
      rpId_45_N_(iConfig.getParameter<unsigned int>("rpId_45_N")),
      rpId_45_F_(iConfig.getParameter<unsigned int>("rpId_45_F")),
      rpId_56_N_(iConfig.getParameter<unsigned int>("rpId_56_N")),
      rpId_56_F_(iConfig.getParameter<unsigned int>("rpId_56_F")),
      outputFile_(iConfig.getParameter<string>("outputFile")) {
  singleArmConfigurations = {
      {rpId_45_N_},
      {rpId_45_F_},
      {rpId_56_N_},
      {rpId_56_F_},
      {rpId_45_N_, rpId_45_F_},
      {rpId_56_N_, rpId_56_F_},
  };

  doubleArmConfigurations = {
      {rpId_45_N_, rpId_56_N_},
      {rpId_45_F_, rpId_56_F_},
      {rpId_45_N_, rpId_45_F_, rpId_56_N_, rpId_56_F_},
  };
}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 158 of file CTPPSAcceptancePlotter.cc.

References protons_cff::arm, CTPPSDetId::arm(), doubleArmConfigurations, doubleArmPlots, edm::HepMCProduct::GetEvent(), iEvent, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, simPPS2016_cfi::rpIds, singleArmConfigurations, singleArmPlots, tokenHepMC_, tokenTracks_, and protons_cff::xi.

                                                                                  {
  // get input
  edm::Handle<edm::HepMCProduct> hHepMC;
  iEvent.getByToken(tokenHepMC_, hHepMC);
  HepMC::GenEvent *hepMCEvent = (HepMC::GenEvent *)hHepMC->GetEvent();

  edm::Handle<CTPPSLocalTrackLiteCollection> hTracks;
  iEvent.getByToken(tokenTracks_, hTracks);

  // extract protons
  bool proton_45_set = false;
  bool proton_56_set = false;
  FourVector mom_45, mom_56;

  for (auto it = hepMCEvent->particles_begin(); it != hepMCEvent->particles_end(); ++it) {
    const auto &part = *it;

    // accept only stable non-beam protons
    if (part->pdg_id() != 2212)
      continue;

    if (part->status() != 1)
      continue;

    if (part->is_beam())
      continue;

    const auto &mom = part->momentum();

    if (mom.e() < 4500.)
      continue;

    if (mom.z() > 0) {
      // 45
      if (proton_45_set) {
        LogError("CTPPSAcceptancePlotter") << "Multiple protons found in sector 45.";
        return;
      }

      proton_45_set = true;
      mom_45 = mom;
    } else {
      // 56
      if (proton_56_set) {
        LogError("CTPPSAcceptancePlotter") << "Multiple protons found in sector 56.";
        return;
      }

      proton_56_set = true;
      mom_56 = mom;
    }
  }

  // stop if protons missing
  if (!proton_45_set || !proton_56_set)
    return;

  // calculate xi's
  const double p_nom = 6500.;
  const double xi_45 = (p_nom - mom_45.e()) / p_nom;
  const double xi_56 = (p_nom - mom_56.e()) / p_nom;

  // process tracks
  map<unsigned int, bool> trackPresent;
  for (const auto &trk : *hTracks) {
    CTPPSDetId rpId(trk.rpId());
    unsigned int rpDecId = rpId.arm() * 100 + rpId.station() * 10 + rpId.rp();
    trackPresent[rpDecId] = true;
  }

  // update plots
  for (const auto &rpIds : singleArmConfigurations) {
    bool acc = true;
    signed int arm = -1;
    for (const auto rpId : rpIds) {
      acc &= trackPresent[rpId];
      arm = rpId / 100;
    }

    if (arm < 0)
      continue;

    const double xi = (arm == 0) ? xi_45 : xi_56;

    singleArmPlots[rpIds].fill(xi, acc);
  }

  for (const auto &rpIds : doubleArmConfigurations) {
    bool acc = true;
    for (const auto rpId : rpIds)
      acc &= trackPresent[rpId];

    doubleArmPlots[rpIds].fill(xi_45, xi_56, acc);
  }
}

beginJob()

void CTPPSAcceptancePlotter::beginJob ( void  )

inlineoverrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 33 of file CTPPSAcceptancePlotter.cc.

{}

endJob()

void CTPPSAcceptancePlotter::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 256 of file CTPPSAcceptancePlotter.cc.

References TrackerOfflineValidation_Dqm_cff::dirName, doubleArmPlots, outputFile_, AlCaHLTBitMon_ParallelJobs::p, and singleArmPlots.

                                    {
  auto f_out = std::make_unique<TFile>(outputFile_.c_str(), "recreate");

  for (const auto &p : singleArmPlots) {
    string dirName;
    for (const auto &rpId : p.first) {
      if (!dirName.empty())
        dirName += ",";
      dirName += Form("%u", rpId);
    }

    gDirectory = f_out->mkdir(dirName.c_str());
    p.second.write();
  }

  for (const auto &p : doubleArmPlots) {
    string dirName;
    for (const auto &rpId : p.first) {
      if (!dirName.empty())
        dirName += ",";
      dirName += Form("%u", rpId);
    }

    gDirectory = f_out->mkdir(dirName.c_str());
    p.second.write();
  }
}

Member Data Documentation

doubleArmConfigurations

std::vector<std::set<unsigned int> > CTPPSAcceptancePlotter::doubleArmConfigurations

private

Definition at line 120 of file CTPPSAcceptancePlotter.cc.

Referenced by analyze(), and CTPPSAcceptancePlotter().

doubleArmPlots

std::map<std::set<unsigned int>, DoubleArmPlots> CTPPSAcceptancePlotter::doubleArmPlots

private

Definition at line 121 of file CTPPSAcceptancePlotter.cc.

Referenced by analyze(), and endJob().

outputFile_

std::string CTPPSAcceptancePlotter::outputFile_

private

Definition at line 42 of file CTPPSAcceptancePlotter.cc.

Referenced by endJob().

rpId_45_F_

unsigned int CTPPSAcceptancePlotter::rpId_45_F_

private

Definition at line 40 of file CTPPSAcceptancePlotter.cc.

Referenced by CTPPSAcceptancePlotter().

rpId_45_N_

unsigned int CTPPSAcceptancePlotter::rpId_45_N_

private

Definition at line 40 of file CTPPSAcceptancePlotter.cc.

Referenced by CTPPSAcceptancePlotter().

rpId_56_F_

unsigned int CTPPSAcceptancePlotter::rpId_56_F_

private

Definition at line 40 of file CTPPSAcceptancePlotter.cc.

Referenced by CTPPSAcceptancePlotter().

rpId_56_N_

unsigned int CTPPSAcceptancePlotter::rpId_56_N_

private

Definition at line 40 of file CTPPSAcceptancePlotter.cc.

Referenced by CTPPSAcceptancePlotter().

singleArmConfigurations

std::vector<std::set<unsigned int> > CTPPSAcceptancePlotter::singleArmConfigurations

private

Definition at line 64 of file CTPPSAcceptancePlotter.cc.

Referenced by analyze(), and CTPPSAcceptancePlotter().

singleArmPlots

std::map<std::set<unsigned int>, SingleArmPlots> CTPPSAcceptancePlotter::singleArmPlots

private

Definition at line 65 of file CTPPSAcceptancePlotter.cc.

Referenced by analyze(), and endJob().

tokenHepMC_

edm::EDGetTokenT<edm::HepMCProduct> CTPPSAcceptancePlotter::tokenHepMC_

private

Definition at line 37 of file CTPPSAcceptancePlotter.cc.

Referenced by analyze().

tokenTracks_

edm::EDGetTokenT<CTPPSLocalTrackLiteCollection> CTPPSAcceptancePlotter::tokenTracks_

private

Definition at line 38 of file CTPPSAcceptancePlotter.cc.

Referenced by analyze().