Module to apply Proton POG quality criteria. More...

Inheritance diagram for PPSFilteredProtonProducer:

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

	~PPSFilteredProtonProducer () override=default

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	endStream () override

bool	IsTrackOK (const CTPPSLocalTrackLite &tr, unsigned int idx, std::ostringstream &log)
	check one track More...

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
unsigned int	n_protons_multi_rp_all

unsigned int	n_protons_multi_rp_kept

unsigned int	n_protons_single_rp_all
	counters More...

unsigned int	n_protons_single_rp_kept

bool	protons_multi_rp_check_valid_fit_

double	protons_multi_rp_chi_sq_max_

bool	protons_multi_rp_include_

edm::EDGetTokenT< reco::ForwardProtonCollection >	protons_multi_rp_input_token_

double	protons_multi_rp_normalised_chi_sq_max_

std::string	protons_multi_rp_output_label_

bool	protons_single_rp_include_

edm::EDGetTokenT< reco::ForwardProtonCollection >	protons_single_rp_input_token_

std::string	protons_single_rp_output_label_

double	tracks_all_local_angle_x_max_

double	tracks_all_local_angle_y_max_

std::vector< unsigned int >	tracks_pixel_forbidden_RecoInfo_values_

double	tracks_pixel_normalised_chi_sq_max_

unsigned int	tracks_pixel_number_of_hits_min_

edm::EDGetTokenT< CTPPSLocalTrackLiteCollection >	tracksToken_

bool	verbosity_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Module to apply Proton POG quality criteria.

Definition at line 21 of file PPSFilteredProtonProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 64 of file PPSFilteredProtonProducer.cc.

References edm::ParameterSet::getParameterSet(), protons_multi_rp_check_valid_fit_, protons_multi_rp_chi_sq_max_, protons_multi_rp_include_, protons_multi_rp_input_token_, protons_multi_rp_normalised_chi_sq_max_, protons_multi_rp_output_label_, protons_cff::protons_single_rp, protons_single_rp_include_, protons_single_rp_input_token_, protons_single_rp_output_label_, AlCaHLTBitMon_QueryRunRegistry::string, tracks_all_local_angle_x_max_, tracks_all_local_angle_y_max_, tracks_pixel_forbidden_RecoInfo_values_, tracks_pixel_normalised_chi_sq_max_, and tracks_pixel_number_of_hits_min_.

     : verbosity_(iConfig.getUntrackedParameter<bool>("verbosity")),
       n_protons_single_rp_all(0),
       n_protons_single_rp_kept(0),
       n_protons_multi_rp_all(0),
       n_protons_multi_rp_kept(0) {
   const auto &tracks_all = iConfig.getParameterSet("tracks_all");
   tracks_all_local_angle_x_max_ = tracks_all.getParameter<double>("local_angle_x_max");
   tracks_all_local_angle_y_max_ = tracks_all.getParameter<double>("local_angle_y_max");
 
   const auto &tracks_pixel = iConfig.getParameterSet("tracks_pixel");
   tracks_pixel_forbidden_RecoInfo_values_ =
       tracks_pixel.getParameter<std::vector<unsigned int>>("forbidden_RecoInfo_values");
   tracks_pixel_number_of_hits_min_ = tracks_pixel.getParameter<unsigned int>("number_of_hits_min");
   tracks_pixel_normalised_chi_sq_max_ = tracks_pixel.getParameter<double>("normalised_chi_sq_max");
 
   const auto &protons_single_rp = iConfig.getParameterSet("protons_single_rp");
   protons_single_rp_include_ = protons_single_rp.getParameter<bool>("include");
   protons_single_rp_input_token_ =
       consumes<reco::ForwardProtonCollection>(protons_single_rp.getParameter<edm::InputTag>("input_tag"));
   protons_single_rp_output_label_ = protons_single_rp.getParameter<std::string>("output_label");
 
   const auto &protons_multi_rp = iConfig.getParameterSet("protons_multi_rp");
   protons_multi_rp_include_ = protons_multi_rp.getParameter<bool>("include");
   protons_multi_rp_input_token_ =
       consumes<reco::ForwardProtonCollection>(protons_multi_rp.getParameter<edm::InputTag>("input_tag"));
   protons_multi_rp_output_label_ = protons_multi_rp.getParameter<std::string>("output_label");
   protons_multi_rp_check_valid_fit_ = protons_multi_rp.getParameter<bool>("check_valid_fit");
   protons_multi_rp_chi_sq_max_ = protons_multi_rp.getParameter<double>("chi_sq_max");
   protons_multi_rp_normalised_chi_sq_max_ = protons_multi_rp.getParameter<double>("normalised_chi_sq_max");
 
   if (protons_single_rp_include_)
     produces<reco::ForwardProtonCollection>(protons_single_rp_output_label_);
 
   if (protons_multi_rp_include_)
     produces<reco::ForwardProtonCollection>(protons_multi_rp_output_label_);
 }

PPSFilteredProtonProducer::~PPSFilteredProtonProducer ( )

overridedefault

Member Function Documentation

void PPSFilteredProtonProducer::endStream ( )

overrideprivate

Definition at line 255 of file PPSFilteredProtonProducer.cc.

References DEFINE_FWK_MODULE, n_protons_multi_rp_all, n_protons_multi_rp_kept, n_protons_single_rp_all, and n_protons_single_rp_kept.

                                           {
   edm::LogInfo("PPS")
       << "single-RP protons: total=" << n_protons_single_rp_all << ", kept=" << n_protons_single_rp_kept
       << " --> keep rate="
       << ((n_protons_single_rp_all > 0) ? double(n_protons_single_rp_kept) / n_protons_single_rp_all * 100. : 0.)
       << "%\n"
       << "multi-RP protons: total=" << n_protons_multi_rp_all << ", kept=" << n_protons_multi_rp_kept
       << " --> keep rate="
       << ((n_protons_multi_rp_all > 0) ? double(n_protons_multi_rp_kept) / n_protons_multi_rp_all * 100. : 0.) << "%";
 }

void PPSFilteredProtonProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 104 of file PPSFilteredProtonProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), allShiftedPlanes, createfilelist::int, mixedPlanes, protons_cff::protons_single_rp, edm::ParameterDescriptionNode::setComment(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                            {
   edm::ParameterSetDescription desc;
 
   desc.addUntracked<bool>("verbosity", false)->setComment("verbosity");
 
   edm::ParameterSetDescription tracks_all;
   tracks_all.add<double>("local_angle_x_max", 0.020)
       ->setComment("maximum absolute value of local horizontal angle, in rad");
   tracks_all.add<double>("local_angle_y_max", 0.020)
       ->setComment("maximum absolute value of local horizontal angle, in rad");
   desc.add<edm::ParameterSetDescription>("tracks_all", tracks_all)->setComment("settings for all tracks");
 
   edm::ParameterSetDescription tracks_pixel;
   const std::vector<unsigned int> def_for_RecoInfo_vals = {
       (unsigned int)CTPPSpixelLocalTrackReconstructionInfo::allShiftedPlanes,
       (unsigned int)CTPPSpixelLocalTrackReconstructionInfo::mixedPlanes};
   tracks_pixel.add<std::vector<unsigned int>>("forbidden_RecoInfo_values", def_for_RecoInfo_vals)
       ->setComment("list of forbidden RecoInfo values");
   tracks_pixel.add<unsigned int>("number_of_hits_min", 0)->setComment("minimum required number of hits");
   tracks_pixel.add<double>("normalised_chi_sq_max", 1E100)->setComment("maximum tolerated chi square / ndof");
   desc.add<edm::ParameterSetDescription>("tracks_pixel", tracks_pixel)
       ->setComment("specific settings for pixel-RP tracks");
 
   edm::ParameterSetDescription protons_single_rp;
   protons_single_rp.add<bool>("include", true)->setComment("flag whether single-RP protons should be processed");
   protons_single_rp.add<edm::InputTag>("input_tag", edm::InputTag("ctppsProtons", "singleRP"))->setComment("input tag");
   protons_single_rp.add<std::string>("output_label", "singleRP")->setComment("output label");
   desc.add<edm::ParameterSetDescription>("protons_single_rp", protons_single_rp)
       ->setComment("settings for single-RP protons");
 
   edm::ParameterSetDescription protons_multi_rp;
   protons_multi_rp.add<bool>("include", true)->setComment("flag whether multi-RP protons should be processed");
   protons_multi_rp.add<edm::InputTag>("input_tag", edm::InputTag("ctppsProtons", "multiRP"))->setComment("input tag");
   protons_multi_rp.add<std::string>("output_label", "multiRP")->setComment("output label");
   protons_multi_rp.add<bool>("check_valid_fit", true)->setComment("flag whether validFit should be checked");
   protons_multi_rp.add<double>("chi_sq_max", 1E-4)->setComment("maximum tolerated value of chi square");
   protons_multi_rp.add<double>("normalised_chi_sq_max", 1E100)
       ->setComment("maximum tolerated value of chi square / ndof, applied only if ndof > 0");
   desc.add<edm::ParameterSetDescription>("protons_multi_rp", protons_multi_rp)
       ->setComment("settings for multi-RP protons");
 
   descriptions.add("ppsFilteredProtonProducer", desc);
 }

bool PPSFilteredProtonProducer::IsTrackOK	(	const CTPPSLocalTrackLite &	tr,
		unsigned int	idx,
		std::ostringstream &	log
	)

private

check one track

Definition at line 150 of file PPSFilteredProtonProducer.cc.

References funct::abs(), spr::find(), CTPPSLocalTrackLite::getChiSquaredOverNDF(), CTPPSLocalTrackLite::getNumberOfPointsUsedForFit(), CTPPSLocalTrackLite::getPixelTrackRecoInfo(), CTPPSLocalTrackLite::getRPId(), CTPPSLocalTrackLite::getTx(), CTPPSLocalTrackLite::getTy(), createfilelist::int, convertSQLiteXML::ok, year_2016_postTS2_cff::rpId, CTPPSDetId::sdTrackingPixel, tracks_all_local_angle_x_max_, tracks_all_local_angle_y_max_, tracks_pixel_forbidden_RecoInfo_values_, tracks_pixel_normalised_chi_sq_max_, tracks_pixel_number_of_hits_min_, and verbosity_.

Referenced by produce().

                                                                                                                 {
   bool ok = true;
 
   // checks for all tracks
   ok &= (std::abs(tr.getTx()) < tracks_all_local_angle_x_max_);
   ok &= (std::abs(tr.getTy()) < tracks_all_local_angle_y_max_);
 
   // pixel checks
   const CTPPSDetId rpId(tr.getRPId());
   if (rpId.subdetId() == CTPPSDetId::sdTrackingPixel) {
     ok &= (find(tracks_pixel_forbidden_RecoInfo_values_.begin(),
                 tracks_pixel_forbidden_RecoInfo_values_.end(),
                 (unsigned int)tr.getPixelTrackRecoInfo()) == tracks_pixel_forbidden_RecoInfo_values_.end());
     ok &= (tr.getNumberOfPointsUsedForFit() >= tracks_pixel_number_of_hits_min_);
     ok &= (tr.getChiSquaredOverNDF() <= tracks_pixel_normalised_chi_sq_max_);
   }
 
   if (!ok && verbosity_)
     log << "track idx=" << idx << " does not fulfil criteria." << std::endl;
 
   return ok;
 }

void PPSFilteredProtonProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 175 of file PPSFilteredProtonProducer.cc.

References edm::Event::get(), IsTrackOK(), keep, eostools::move(), n_protons_multi_rp_all, n_protons_multi_rp_kept, n_protons_single_rp_all, n_protons_single_rp_kept, protons_multi_rp_check_valid_fit_, protons_multi_rp_chi_sq_max_, protons_multi_rp_include_, protons_multi_rp_input_token_, protons_multi_rp_normalised_chi_sq_max_, protons_multi_rp_output_label_, protons_single_rp_include_, protons_single_rp_input_token_, protons_single_rp_output_label_, edm::Event::put(), and verbosity_.

                                                                                      {
   std::ostringstream ssLog;
 
   // process single-RP protons
   if (protons_single_rp_include_) {
     reco::ForwardProtonCollection const &hInputProtons = iEvent.get(protons_single_rp_input_token_);
     std::unique_ptr<reco::ForwardProtonCollection> pOutputProtons(new reco::ForwardProtonCollection);
 
     for (const auto &proton : hInputProtons) {
       bool keep = true;
 
       // no specific checks for single-RP protons
 
       // test contributing tracks
       for (const auto &tr_ref : proton.contributingLocalTracks()) {
         if (!keep)
           break;
 
         keep &= IsTrackOK(*tr_ref, tr_ref.key(), ssLog);
       }
 
       n_protons_single_rp_all++;
 
       if (keep) {
         n_protons_single_rp_kept++;
         pOutputProtons->push_back(proton);
       } else {
         if (verbosity_)
           ssLog << "single-RP proton idx=" << n_protons_single_rp_all - 1 << " excluded." << std::endl;
       }
     }
 
     iEvent.put(std::move(pOutputProtons), protons_single_rp_output_label_);
   }
 
   // process multi-RP protons
   if (protons_multi_rp_include_) {
     reco::ForwardProtonCollection const &hInputProtons = iEvent.get(protons_multi_rp_input_token_);
     std::unique_ptr<reco::ForwardProtonCollection> pOutputProtons(new reco::ForwardProtonCollection);
 
     for (const auto &proton : hInputProtons) {
       bool keep = true;
 
       // multi-RP proton checks
       if (protons_multi_rp_check_valid_fit_)
         keep &= proton.validFit();
 
       keep &= (proton.chi2() <= protons_multi_rp_chi_sq_max_);
 
       if (proton.ndof() > 0)
         keep &= (proton.normalizedChi2() <= protons_multi_rp_normalised_chi_sq_max_);
 
       // test contributing tracks
       for (const auto &tr_ref : proton.contributingLocalTracks()) {
         if (!keep)
           break;
 
         keep &= IsTrackOK(*tr_ref, tr_ref.key(), ssLog);
       }
 
       n_protons_multi_rp_all++;
 
       if (keep) {
         n_protons_multi_rp_kept++;
         pOutputProtons->push_back(proton);
       } else {
         if (verbosity_)
           ssLog << "multi-RP proton idx=" << n_protons_multi_rp_all - 1 << " excluded." << std::endl;
       }
     }
 
     iEvent.put(std::move(pOutputProtons), protons_multi_rp_output_label_);
   }
 
   if (verbosity_ && !ssLog.str().empty())
     edm::LogInfo("PPS") << ssLog.str();
 }