Inheritance diagram for CTPPSHepMCDistributionPlotter:

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

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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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 &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

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::unique_ptr< TH1D >	h_th_x_

std::unique_ptr< TH1D >	h_th_y_

std::unique_ptr< TH1D >	h_vtx_x_

std::unique_ptr< TH1D >	h_vtx_y_

std::unique_ptr< TH1D >	h_vtx_z_

std::unique_ptr< TH1D >	h_xi_

std::string	lhcInfoLabel_

std::string	outputFile_

edm::EDGetTokenT< edm::HepMCProduct >	tokenHepMC_

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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 27 of file CTPPSHepMCDistributionPlotter.cc.

Constructor & Destructor Documentation

◆ CTPPSHepMCDistributionPlotter()

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

explicit

Definition at line 51 of file CTPPSHepMCDistributionPlotter.cc.

     : tokenHepMC_(consumes<edm::HepMCProduct>(iConfig.getParameter<edm::InputTag>("tagHepMC"))),
       lhcInfoLabel_(iConfig.getParameter<std::string>("lhcInfoLabel")),
       outputFile_(iConfig.getParameter<string>("outputFile")),
  
       h_vtx_x_(new TH1D("h_vtx_x", ";vtx_x   (mm)", 100, 0., 0.)),
       h_vtx_y_(new TH1D("h_vtx_y", ";vtx_y   (mm)", 100, 0., 0.)),
       h_vtx_z_(new TH1D("h_vtx_z", ";vtx_y   (mm)", 100, 0., 0.)),
  
       h_xi_(new TH1D("h_xi", ";#xi", 100, 0., 0.30)),
       h_th_x_(new TH1D("h_th_x", ";#theta^{*}_{x}", 100, -300E-6, +300E-6)),

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 66 of file CTPPSHepMCDistributionPlotter.cc.

                                                                                                {
   // get conditions
   edm::ESHandle<LHCInfo> hLHCInfo;
   iSetup.get<LHCInfoRcd>().get(lhcInfoLabel_, hLHCInfo);
  
   // get input
   edm::Handle<edm::HepMCProduct> hHepMC;
   iEvent.getByToken(tokenHepMC_, hHepMC);
   HepMC::GenEvent *hepMCEvent = (HepMC::GenEvent *)hHepMC->GetEvent();
  
   // plot vertices
   for (HepMC::GenEvent::vertex_iterator vit = hepMCEvent->vertices_begin(); vit != hepMCEvent->vertices_end(); ++vit) {
     const auto pos = (*vit)->position();
     h_vtx_x_->Fill(pos.x());
     h_vtx_y_->Fill(pos.y());
     h_vtx_z_->Fill(pos.z());
   }
  
   // extract protons
   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();
     const double p_nom = hLHCInfo->energy();
  
     if (mom.rho() / p_nom < 0.7)
       continue;
  
     const double xi_simu = (p_nom - mom.e()) / p_nom;
     const double th_x_simu = mom.x() / mom.rho();
     const double th_y_simu = mom.y() / mom.rho();
  
     h_xi_->Fill(xi_simu);
     h_th_x_->Fill(th_x_simu);
     h_th_y_->Fill(th_y_simu);
   }

References LHCInfo::energy(), edm::EventSetup::get(), get, edm::HepMCProduct::GetEvent(), h_th_x_, h_th_y_, h_vtx_x_, h_vtx_y_, h_vtx_z_, h_xi_, iEvent, lhcInfoLabel_, and tokenHepMC_.

◆ endJob()

void CTPPSHepMCDistributionPlotter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 116 of file CTPPSHepMCDistributionPlotter.cc.

                                            {
   auto f_out = std::make_unique<TFile>(outputFile_.c_str(), "recreate");
  
   h_vtx_x_->Write();
   h_vtx_y_->Write();
   h_vtx_z_->Write();
  
   h_xi_->Write();
   h_th_x_->Write();
   h_th_y_->Write();