Inheritance diagram for CTPPSProtonReconstructionPlotter:

Classes
struct	ArmCorrelationPlots

struct	MultiRPPlots

struct	SingleMultiCorrelationPlots

struct	SingleRPPlots

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

	~CTPPSProtonReconstructionPlotter () 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

Static Private Member Functions
static void	profileToRMSGraph (TProfile p, TGraphErrors g)

Private Attributes
std::map< unsigned int, ArmCorrelationPlots >	armCorrelationPlots_

signed int	maxNonEmptyEvents_

std::map< unsigned int, MultiRPPlots >	multiRPPlots_

signed int	n_non_empty_events_

std::string	outputFile_

std::unique_ptr< TProfile >	p_x_L_diffNF_vs_x_L_N_

std::unique_ptr< TProfile >	p_x_R_diffNF_vs_x_R_N_

std::unique_ptr< TProfile >	p_y_L_diffNF_vs_y_L_N_

std::unique_ptr< TProfile >	p_y_R_diffNF_vs_y_R_N_

unsigned int	rpId_45_F_

unsigned int	rpId_45_N_

unsigned int	rpId_56_F_

unsigned int	rpId_56_N_

std::map< unsigned int, SingleMultiCorrelationPlots >	singleMultiCorrelationPlots_

std::map< unsigned int, SingleRPPlots >	singleRPPlots_

edm::EDGetTokenT< reco::ForwardProtonCollection >	tokenRecoProtonsMultiRP_

edm::EDGetTokenT< reco::ForwardProtonCollection >	tokenRecoProtonsSingleRP_

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)

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 29 of file CTPPSProtonReconstructionPlotter.cc.

Constructor & Destructor Documentation

CTPPSProtonReconstructionPlotter::CTPPSProtonReconstructionPlotter ( const edm::ParameterSet & ps )

explicit

Definition at line 344 of file CTPPSProtonReconstructionPlotter.cc.

                                                                                             :
   tokenTracks_             (consumes<CTPPSLocalTrackLiteCollection>(ps.getParameter<edm::InputTag>("tagTracks"))),
   tokenRecoProtonsSingleRP_(consumes<reco::ForwardProtonCollection>(ps.getParameter<InputTag>("tagRecoProtonsSingleRP"))),
   tokenRecoProtonsMultiRP_ (consumes<reco::ForwardProtonCollection>(ps.getParameter<InputTag>("tagRecoProtonsMultiRP"))),
 
   rpId_45_N_(ps.getParameter<unsigned int>("rpId_45_N")),
   rpId_45_F_(ps.getParameter<unsigned int>("rpId_45_F")),
   rpId_56_N_(ps.getParameter<unsigned int>("rpId_56_N")),
   rpId_56_F_(ps.getParameter<unsigned int>("rpId_56_F")),
 
   outputFile_(ps.getParameter<string>("outputFile")),
   maxNonEmptyEvents_(ps.getUntrackedParameter<signed int>("maxNonEmptyEvents", -1)),
 
   p_x_L_diffNF_vs_x_L_N_(new TProfile("p_x_L_diffNF_vs_x_L_N", ";x_{LN};x_{LF} - x_{LN}", 100, 0., +20.)),
   p_x_R_diffNF_vs_x_R_N_(new TProfile("p_x_R_diffNF_vs_x_R_N", ";x_{RN};x_{RF} - x_{RN}", 100, 0., +20.)),
 
   p_y_L_diffNF_vs_y_L_N_(new TProfile("p_y_L_diffNF_vs_y_L_N", ";y_{LN};y_{LF} - y_{LN}", 100, -20., +20.)),
   p_y_R_diffNF_vs_y_R_N_(new TProfile("p_y_R_diffNF_vs_y_R_N", ";y_{RN};y_{RF} - y_{RN}", 100, -20., +20.)),
 
   n_non_empty_events_(0)
 {}

CTPPSProtonReconstructionPlotter::~CTPPSProtonReconstructionPlotter ( )

inlineoverride

Definition at line 33 of file CTPPSProtonReconstructionPlotter.cc.

References analyze(), and endJob().

33 {}

Member Function Documentation

void CTPPSProtonReconstructionPlotter::analyze	(	const edm::Event &	event,
		const edm::EventSetup &
	)

overrideprivate

Definition at line 368 of file CTPPSProtonReconstructionPlotter.cc.

References CTPPSDetId::arm(), armCorrelationPlots_, Exception, CTPPSLocalTrackLite::getX(), CTPPSLocalTrackLite::getY(), training_settings::idx, maxNonEmptyEvents_, multiRPPlots_, n_non_empty_events_, p_x_L_diffNF_vs_x_L_N_, p_x_R_diffNF_vs_x_R_N_, p_y_L_diffNF_vs_y_L_N_, p_y_R_diffNF_vs_y_R_N_, rpId_45_F_, rpId_45_N_, rpId_56_F_, rpId_56_N_, singleMultiCorrelationPlots_, singleRPPlots_, tokenRecoProtonsMultiRP_, tokenRecoProtonsSingleRP_, and tokenTracks_.

Referenced by ~CTPPSProtonReconstructionPlotter().

 {
   // get input
   edm::Handle<CTPPSLocalTrackLiteCollection> hTracks;
   event.getByToken(tokenTracks_, hTracks);
 
   Handle<reco::ForwardProtonCollection> hRecoProtonsSingleRP;
   event.getByToken(tokenRecoProtonsSingleRP_, hRecoProtonsSingleRP);
 
   Handle<reco::ForwardProtonCollection> hRecoProtonsMultiRP;
   event.getByToken(tokenRecoProtonsMultiRP_, hRecoProtonsMultiRP);
 
   if (!hRecoProtonsSingleRP->empty())
     n_non_empty_events_++;
 
   if (maxNonEmptyEvents_ > 0 && n_non_empty_events_ > maxNonEmptyEvents_)
     throw cms::Exception("CTPPSProtonReconstructionPlotter") << "Number of non empty events reached maximum.";
 
   // track plots
   const CTPPSLocalTrackLite *tr_L_N = nullptr;
   const CTPPSLocalTrackLite *tr_L_F = nullptr;
   const CTPPSLocalTrackLite *tr_R_N = nullptr;
   const CTPPSLocalTrackLite *tr_R_F = nullptr;
 
   for (const auto &tr : *hTracks)
   {
     CTPPSDetId rpId(tr.getRPId());
     unsigned int decRPId = rpId.arm()*100 + rpId.station()*10 + rpId.rp();
 
     if (decRPId == rpId_45_N_) tr_L_N = &tr;
     if (decRPId == rpId_45_F_) tr_L_F = &tr;
     if (decRPId == rpId_56_N_) tr_R_N = &tr;
     if (decRPId == rpId_56_F_) tr_R_F = &tr;
   }
 
   if (tr_L_N && tr_L_F)
   {
     p_x_L_diffNF_vs_x_L_N_->Fill(tr_L_N->getX(), tr_L_F->getX() - tr_L_N->getX());
     p_y_L_diffNF_vs_y_L_N_->Fill(tr_L_N->getY(), tr_L_F->getY() - tr_L_N->getY());
   }
 
   if (tr_R_N && tr_R_F)
   {
     p_x_R_diffNF_vs_x_R_N_->Fill(tr_R_N->getX(), tr_R_F->getX() - tr_R_N->getX());
     p_y_R_diffNF_vs_y_R_N_->Fill(tr_R_N->getY(), tr_R_F->getY() - tr_R_N->getY());
   }
 
   // make single-RP-reco plots
   for (const auto &proton : *hRecoProtonsSingleRP)
   {
     CTPPSDetId rpId((*proton.contributingLocalTracks().begin())->getRPId());
     unsigned int decRPId = rpId.arm()*100 + rpId.station()*10 + rpId.rp();
     singleRPPlots_[decRPId].fill(proton);
   }
 
   // make multi-RP-reco plots
   for (const auto &proton : *hRecoProtonsMultiRP)
   {
     CTPPSDetId rpId((*proton.contributingLocalTracks().begin())->getRPId());
     unsigned int armId = rpId.arm();
     multiRPPlots_[armId].fill(proton);
   }
 
   // make correlation plots
   for (const auto &proton_s : *hRecoProtonsSingleRP)
   {
     for (const auto &proton_m : *hRecoProtonsMultiRP)
     {
       // only compare object from the same arm
       CTPPSDetId rpId_s((*proton_s.contributingLocalTracks().begin())->getRPId());
       CTPPSDetId rpId_m((*proton_m.contributingLocalTracks().begin())->getRPId());
 
       if (rpId_s.arm() != rpId_m.arm())
         continue;
 
       // build index
       const unsigned int idx = rpId_s.arm()*1000 + rpId_s.station()*100 + rpId_s.rp()*10 + rpId_m.arm();
 
       // fill plots
       singleMultiCorrelationPlots_[idx].fill(proton_s, proton_m);
     }
   }
 
   // arm correlation plots
   const reco::ForwardProton *p_arm0_s_N = nullptr, *p_arm0_s_F = nullptr, *p_arm0_m = nullptr;
   const reco::ForwardProton *p_arm1_s_N = nullptr, *p_arm1_s_F = nullptr, *p_arm1_m = nullptr;
 
   for (const auto &proton : *hRecoProtonsSingleRP)
   {
     CTPPSDetId rpId((*proton.contributingLocalTracks().begin())->getRPId());
     const unsigned int rpDecId = rpId.arm()*100 + rpId.station()*10 + rpId.rp();
 
     if (rpDecId == rpId_45_N_) p_arm0_s_N = &proton;
     if (rpDecId == rpId_45_F_) p_arm0_s_F = &proton;
 
     if (rpDecId == rpId_56_N_) p_arm1_s_N = &proton;
     if (rpDecId == rpId_56_F_) p_arm1_s_F = &proton;
   }
 
   for (const auto & proton : *hRecoProtonsMultiRP)
   {
     CTPPSDetId rpId((*proton.contributingLocalTracks().begin())->getRPId());
     const unsigned int arm = rpId.arm();
 
     if (arm == 0) p_arm0_m = &proton;
     if (arm == 1) p_arm1_m = &proton;
   }
 
   if (p_arm0_s_N && p_arm0_s_F && p_arm0_m)
     armCorrelationPlots_[0].fill(*p_arm0_s_N, *p_arm0_s_F, *p_arm0_m);
 
   if (p_arm1_s_N && p_arm1_s_F && p_arm1_m)
     armCorrelationPlots_[1].fill(*p_arm1_s_N, *p_arm1_s_F, *p_arm1_m);
 }

void CTPPSProtonReconstructionPlotter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 485 of file CTPPSProtonReconstructionPlotter.cc.

References armCorrelationPlots_, DEFINE_FWK_MODULE, multiRPPlots_, n_non_empty_events_, outputFile_, p_x_L_diffNF_vs_x_L_N_, p_x_R_diffNF_vs_x_R_N_, p_y_L_diffNF_vs_y_L_N_, p_y_R_diffNF_vs_y_R_N_, singleMultiCorrelationPlots_, and singleRPPlots_.

Referenced by ~CTPPSProtonReconstructionPlotter().

 {
   LogInfo("CTPPSProtonReconstructionPlotter") << "n_non_empty_events = " << n_non_empty_events_;
 
   auto f_out = std::make_unique<TFile>(outputFile_.c_str(), "recreate");
 
   p_x_L_diffNF_vs_x_L_N_->Write();
   p_x_R_diffNF_vs_x_R_N_->Write();
 
   p_y_L_diffNF_vs_y_L_N_->Write();
   p_y_R_diffNF_vs_y_R_N_->Write();
 
   TDirectory *d_singleRPPlots = f_out->mkdir("singleRPPlots");
   for (const auto& it : singleRPPlots_) {
     gDirectory = d_singleRPPlots->mkdir(Form("rp%u", it.first));
     it.second.write();
   }
 
   TDirectory *d_multiRPPlots = f_out->mkdir("multiRPPlots");
   for (const auto& it : multiRPPlots_) {
     gDirectory = d_multiRPPlots->mkdir(Form("arm%u", it.first));
     it.second.write();
   }
 
   TDirectory *d_singleMultiCorrelationPlots = f_out->mkdir("singleMultiCorrelationPlots");
   for (const auto& it : singleMultiCorrelationPlots_) {
     unsigned int si_rp = it.first / 10;
     unsigned int mu_arm = it.first % 10;
 
     gDirectory = d_singleMultiCorrelationPlots->mkdir(Form("si_rp%u_mu_arm%u", si_rp, mu_arm));
     it.second.write();
   }
 
   TDirectory *d_armCorrelationPlots = f_out->mkdir("armCorrelationPlots");
   for (const auto& it : armCorrelationPlots_) {
     gDirectory = d_armCorrelationPlots->mkdir(Form("arm%u", it.first));
     it.second.write();
   }
 }

static void CTPPSProtonReconstructionPlotter::profileToRMSGraph	(	TProfile *	p,
		TGraphErrors *	g
	)

inlinestaticprivate

Definition at line 51 of file CTPPSProtonReconstructionPlotter.cc.

References EnergyCorrector::c, training_settings::idx, N, and mathSSE::sqrt().

Referenced by CTPPSProtonReconstructionPlotter::MultiRPPlots::write().

     {
       for (int bi = 1; bi <= p->GetNbinsX(); ++bi)
       {
         double c = p->GetBinCenter(bi);
 
         double N = p->GetBinEntries(bi);
         double Sy = p->GetBinContent(bi) * N;
         double Syy = p->GetSumw2()->At(bi);
 
         double si_sq = Syy/N - Sy*Sy/N/N;
         double si = (si_sq >= 0.) ? sqrt(si_sq) : 0.;
         double si_unc_sq = si_sq / 2. / N;  // Gaussian approximation
         double si_unc = (si_unc_sq >= 0.) ? sqrt(si_unc_sq) : 0.;
 
         int idx = g->GetN();
         g->SetPoint(idx, c, si);
         g->SetPointError(idx, 0., si_unc);
       }
     }