Inheritance diagram for CTPPSBeamSmearingValidator:

Classes
struct	SectorPlots

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

	~CTPPSBeamSmearingValidator () 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::unique_ptr< TH1D >	h_de_vtx_x_

std::unique_ptr< TH1D >	h_de_vtx_y_

std::unique_ptr< TH1D >	h_de_vtx_z_

std::string	outputFile_

std::map< unsigned int, SectorPlots >	sectorPlots_

edm::EDGetTokenT< edm::HepMCProduct >	tokenAfterSmearing_

edm::EDGetTokenT< edm::HepMCProduct >	tokenBeforeSmearing_

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 22 of file CTPPSBeamSmearingValidator.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 69 of file CTPPSBeamSmearingValidator.cc.

                                                                                      :
   tokenBeforeSmearing_( consumes<edm::HepMCProduct>(iConfig.getParameter<edm::InputTag>("tagBeforeSmearing")) ),
   tokenAfterSmearing_( consumes<edm::HepMCProduct>(iConfig.getParameter<edm::InputTag>("tagAfterSmearing")) ),
   outputFile_(iConfig.getParameter<string>("outputFile")),
   h_de_vtx_x_(new TH1D("h_de_vtx_x", ";#Delta vtx_{x}   (mm)", 100, 0., 0.)),
   h_de_vtx_y_(new TH1D("h_de_vtx_y", ";#Delta vtx_{y}   (mm)", 100, 0., 0.)),
   h_de_vtx_z_(new TH1D("h_de_vtx_z", ";#Delta vtx_{z}   (mm)", 100, 0., 0.))
 {}

CTPPSBeamSmearingValidator::~CTPPSBeamSmearingValidator ( )

inlineoverride

Definition at line 27 of file CTPPSBeamSmearingValidator.cc.

References analyze(), and endJob().

27 {}

Member Function Documentation

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

overrideprivate

Definition at line 80 of file CTPPSBeamSmearingValidator.cc.

References edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), h_de_vtx_x_, h_de_vtx_y_, h_de_vtx_z_, training_settings::idx, M_PI, gen::n, connectstrParser::o, sectorPlots_, tokenAfterSmearing_, and tokenBeforeSmearing_.

Referenced by ~CTPPSBeamSmearingValidator().

 {
   // get input
   edm::Handle<edm::HepMCProduct> hBeforeSmearing;
   iEvent.getByToken(tokenBeforeSmearing_, hBeforeSmearing);
   HepMC::GenEvent *orig = (HepMC::GenEvent *) hBeforeSmearing->GetEvent();
 
   edm::Handle<edm::HepMCProduct> hAfterSmearing;
   iEvent.getByToken(tokenAfterSmearing_, hAfterSmearing);
   HepMC::GenEvent *smear = (HepMC::GenEvent *) hAfterSmearing->GetEvent();
 
   // vertices
   GenEvent::vertex_const_iterator vold, vnew;
   for (vold = orig->vertices_begin(), vnew = smear->vertices_begin();
       vold != orig->vertices_end() && vnew != smear->vertices_end(); ++vold, ++vnew)
   {
     const FourVector &vo = (*vold)->position();
     const FourVector &vn = (*vnew)->position();
 
     // HepMC gives vertex in mm
     h_de_vtx_x_->Fill(vn.x() - vo.x());
     h_de_vtx_y_->Fill(vn.y() - vo.y());
     h_de_vtx_z_->Fill(vn.z() - vo.z());
   }
 
   // particles
   GenEvent::particle_const_iterator pold, pnew;
   for (pold = orig->particles_begin(), pnew = smear->particles_begin();
       pold != orig->particles_end() && pnew != smear->particles_end(); ++pold, ++pnew)
   {
     FourVector o = (*pold)->momentum(), n = (*pnew)->momentum();
 
     // determine direction region
     signed int idx = -1;
     const double thetaLim = 0.01; // rad
     double th = o.theta();
 
     if (th < thetaLim)
       idx = 0;
     if (th > (M_PI - thetaLim))
       idx = 1;
 
     if (idx < 0)
       continue;
 
     /*
         cout << "particle\n\told: [" << o.x() << ", " << o.y() << ", " << o.z() << ", " << o.t()
         << "]\n\tnew: [" << n.x() << ", " << n.y() << ", " << n.z() << ", " << n.t()
         << "]\n\tregion: " << idx << endl;
     */
 
     // fill histograms
     auto &sp = sectorPlots_[idx];
 
     double othx = o.x() / o.rho(), othy = o.y() / o.rho();
     double nthx = n.x() / n.rho(), nthy = n.y() / n.rho();
 
     sp.h_de_p->Fill(n.rho() - o.rho());
 
     sp.h_de_th_x->Fill(nthx - othx);
     sp.h_de_th_y->Fill(nthy - othy);
   }
 }

void CTPPSBeamSmearingValidator::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 146 of file CTPPSBeamSmearingValidator.cc.

References DEFINE_FWK_MODULE, h_de_vtx_x_, h_de_vtx_y_, h_de_vtx_z_, outputFile_, and sectorPlots_.

Referenced by ~CTPPSBeamSmearingValidator().

 {
   auto f_out = std::make_unique<TFile>(outputFile_.c_str(), "recreate");
 
   h_de_vtx_x_->Write();
   h_de_vtx_y_->Write();
   h_de_vtx_z_->Write();
 
   gDirectory = f_out->mkdir("sector 45");
   sectorPlots_[0].write();
 
   gDirectory = f_out->mkdir("sector 56");
   sectorPlots_[1].write();
 }