Inheritance diagram for DiMuonValidation:

Public Types
enum	Variable { CosThetaCS = 0, DeltaEta = 1, EtaMinus = 2, EtaPlus = 3, PhiCS = 4, PhiMinus = 5, PhiPlus = 6, Pt = 7, VarNumber = 8 }

Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	DiMuonValidation (const edm::ParameterSet &pset)

	~DiMuonValidation () override=default

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () 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

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)

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

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Static Public Attributes
static constexpr double	mu_mass2_ = 0.105658 * 0.105658

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

Private Attributes
int	compressionSettings_

float	eBeam_

DiLeptonHelp::PlotsVsDiLeptonRegion	InvMassInEtaBins = DiLeptonHelp::PlotsVsDiLeptonRegion(1.5)

DiLeptonHelp::PlotsVsDiLeptonRegion	InvMassVsCosThetaCSInEtaBins = DiLeptonHelp::PlotsVsDiLeptonRegion(1.5)

DiLeptonHelp::PlotsVsDiLeptonRegion	InvMassVsPhiMinusInEtaBins = DiLeptonHelp::PlotsVsDiLeptonRegion(1.5)

DiLeptonHelp::PlotsVsDiLeptonRegion	InvMassVsPhiPlusInEtaBins = DiLeptonHelp::PlotsVsDiLeptonRegion(1.5)

double	pair_etamaxneg_

double	pair_etamaxpos_

double	pair_etaminneg_

double	pair_etaminpos_

double	pair_mass_max_

double	pair_mass_min_

int	pair_mass_nbins_

TH1F *	th1f_mass

TH2D *	th2d_mass_variables_ [Variable::VarNumber]

TH3D *	th3d_mass_vs_eta_phi_minus_

TH3D *	th3d_mass_vs_eta_phi_plus_

edm::EDGetTokenT< reco::TrackCollection >	theTrackCollectionToken_

std::string	TkTag_

int	variable_CosThetaCS_nbins_

double	variable_CosThetaCS_xmax_

double	variable_CosThetaCS_xmin_

int	variable_DeltaEta_nbins_

double	variable_DeltaEta_xmax_

double	variable_DeltaEta_xmin_

int	variable_EtaMinus_nbins_

double	variable_EtaMinus_xmax_

double	variable_EtaMinus_xmin_

int	variable_EtaPlus_nbins_

double	variable_EtaPlus_xmax_

double	variable_EtaPlus_xmin_

int	variable_PairPt_nbins_

double	variable_PairPt_xmax_

double	variable_PairPt_xmin_

int	variable_PhiCS_nbins_

double	variable_PhiCS_xmax_

double	variable_PhiCS_xmin_

int	variable_PhiMinus_nbins_

double	variable_PhiMinus_xmax_

double	variable_PhiMinus_xmin_

int	variable_PhiPlus_nbins_

double	variable_PhiPlus_xmax_

double	variable_PhiPlus_xmin_

int	variables_bins_number_ [Variable::VarNumber]

double	variables_max_ [Variable::VarNumber]

double	variables_min_ [Variable::VarNumber]

std::string	variables_name_ [Variable::VarNumber]

std::string	variables_title_ [Variable::VarNumber]

Additional Inherited Members
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 35 of file DiMuonValidation.cc.

Member Enumeration Documentation

◆ Variable

enum DiMuonValidation::Variable

Enumerator
CosThetaCS
DeltaEta
EtaMinus
EtaPlus
PhiCS
PhiMinus
PhiPlus
Pt
VarNumber

Definition at line 108 of file DiMuonValidation.cc.

                 {
     CosThetaCS = 0,
     DeltaEta = 1,
     EtaMinus = 2,
     EtaPlus = 3,
     PhiCS = 4,
     PhiMinus = 5,
     PhiPlus = 6,
     Pt = 7,
     VarNumber = 8
   };

Constructor & Destructor Documentation

◆ DiMuonValidation()

DiMuonValidation::DiMuonValidation ( const edm::ParameterSet & pset )

inlineexplicit

Definition at line 37 of file DiMuonValidation.cc.

References HLT_2024v12_cff::DeltaEta, TFileService::kSharedResource, theTrackCollectionToken_, TkTag_, variable_CosThetaCS_nbins_, variable_CosThetaCS_xmax_, variable_CosThetaCS_xmin_, variable_DeltaEta_nbins_, variable_DeltaEta_xmax_, variable_DeltaEta_xmin_, variable_EtaMinus_nbins_, variable_EtaMinus_xmax_, variable_EtaMinus_xmin_, variable_EtaPlus_nbins_, variable_EtaPlus_xmax_, variable_EtaPlus_xmin_, variable_PairPt_nbins_, variable_PairPt_xmax_, variable_PairPt_xmin_, variable_PhiCS_nbins_, variable_PhiCS_xmax_, variable_PhiCS_xmin_, variable_PhiMinus_nbins_, variable_PhiMinus_xmax_, variable_PhiMinus_xmin_, variable_PhiPlus_nbins_, variable_PhiPlus_xmax_, variable_PhiPlus_xmin_, variables_bins_number_, variables_max_, and variables_min_.

       : eBeam_(pset.getParameter<double>("eBeam")),
         compressionSettings_(pset.getUntrackedParameter<int>("compressionSettings", -1)),
         TkTag_(pset.getParameter<std::string>("TkTag")),
         pair_mass_min_(pset.getParameter<double>("Pair_mass_min")),
         pair_mass_max_(pset.getParameter<double>("Pair_mass_max")),
         pair_mass_nbins_(pset.getParameter<int>("Pair_mass_nbins")),
         pair_etaminpos_(pset.getParameter<double>("Pair_etaminpos")),
         pair_etamaxpos_(pset.getParameter<double>("Pair_etamaxpos")),
         pair_etaminneg_(pset.getParameter<double>("Pair_etaminneg")),
         pair_etamaxneg_(pset.getParameter<double>("Pair_etamaxneg")),
         variable_CosThetaCS_xmin_(pset.getParameter<double>("Variable_CosThetaCS_xmin")),
         variable_CosThetaCS_xmax_(pset.getParameter<double>("Variable_CosThetaCS_xmax")),
         variable_CosThetaCS_nbins_(pset.getParameter<int>("Variable_CosThetaCS_nbins")),
         variable_DeltaEta_xmin_(pset.getParameter<double>("Variable_DeltaEta_xmin")),
         variable_DeltaEta_xmax_(pset.getParameter<double>("Variable_DeltaEta_xmax")),
         variable_DeltaEta_nbins_(pset.getParameter<int>("Variable_DeltaEta_nbins")),
         variable_EtaMinus_xmin_(pset.getParameter<double>("Variable_EtaMinus_xmin")),
         variable_EtaMinus_xmax_(pset.getParameter<double>("Variable_EtaMinus_xmax")),
         variable_EtaMinus_nbins_(pset.getParameter<int>("Variable_EtaMinus_nbins")),
         variable_EtaPlus_xmin_(pset.getParameter<double>("Variable_EtaPlus_xmin")),
         variable_EtaPlus_xmax_(pset.getParameter<double>("Variable_EtaPlus_xmax")),
         variable_EtaPlus_nbins_(pset.getParameter<int>("Variable_EtaPlus_nbins")),
         variable_PhiCS_xmin_(pset.getParameter<double>("Variable_PhiCS_xmin")),
         variable_PhiCS_xmax_(pset.getParameter<double>("Variable_PhiCS_xmax")),
         variable_PhiCS_nbins_(pset.getParameter<int>("Variable_PhiCS_nbins")),
         variable_PhiMinus_xmin_(pset.getParameter<double>("Variable_PhiMinus_xmin")),
         variable_PhiMinus_xmax_(pset.getParameter<double>("Variable_PhiMinus_xmax")),
         variable_PhiMinus_nbins_(pset.getParameter<int>("Variable_PhiMinus_nbins")),
         variable_PhiPlus_xmin_(pset.getParameter<double>("Variable_PhiPlus_xmin")),
         variable_PhiPlus_xmax_(pset.getParameter<double>("Variable_PhiPlus_xmax")),
         variable_PhiPlus_nbins_(pset.getParameter<int>("Variable_PhiPlus_nbins")),
         variable_PairPt_xmin_(pset.getParameter<double>("Variable_PairPt_xmin")),
         variable_PairPt_xmax_(pset.getParameter<double>("Variable_PairPt_xmax")),
         variable_PairPt_nbins_(pset.getParameter<int>("Variable_PairPt_nbins")) {
     usesResource(TFileService::kSharedResource);
     theTrackCollectionToken_ = consumes<reco::TrackCollection>(TkTag_);
 
     variables_min_[Variable::CosThetaCS] = variable_CosThetaCS_xmin_;
     variables_min_[Variable::DeltaEta] = variable_DeltaEta_xmin_;
     variables_min_[Variable::EtaMinus] = variable_EtaMinus_xmin_;
     variables_min_[Variable::EtaPlus] = variable_EtaPlus_xmin_;
     variables_min_[Variable::PhiCS] = variable_PhiCS_xmin_;
     variables_min_[Variable::PhiMinus] = variable_PhiMinus_xmin_;
     variables_min_[Variable::PhiPlus] = variable_PhiPlus_xmin_;
     variables_min_[Variable::Pt] = variable_PairPt_xmin_;
 
     variables_max_[Variable::CosThetaCS] = variable_CosThetaCS_xmax_;
     variables_max_[Variable::DeltaEta] = variable_DeltaEta_xmax_;
     variables_max_[Variable::EtaMinus] = variable_EtaMinus_xmax_;
     variables_max_[Variable::EtaPlus] = variable_EtaPlus_xmax_;
     variables_max_[Variable::PhiCS] = variable_PhiCS_xmax_;
     variables_max_[Variable::PhiMinus] = variable_PhiMinus_xmax_;
     variables_max_[Variable::PhiPlus] = variable_PhiPlus_xmax_;
     variables_max_[Variable::Pt] = variable_PairPt_xmax_;
 
     variables_bins_number_[Variable::CosThetaCS] = variable_CosThetaCS_nbins_;
     variables_bins_number_[Variable::DeltaEta] = variable_DeltaEta_nbins_;
     variables_bins_number_[Variable::EtaMinus] = variable_EtaMinus_nbins_;
     variables_bins_number_[Variable::EtaPlus] = variable_EtaPlus_nbins_;
     variables_bins_number_[Variable::PhiCS] = variable_PhiCS_nbins_;
     variables_bins_number_[Variable::PhiMinus] = variable_PhiMinus_nbins_;
     variables_bins_number_[Variable::PhiPlus] = variable_PhiPlus_nbins_;
     variables_bins_number_[Variable::Pt] = variable_PairPt_nbins_;
   }

◆ ~DiMuonValidation()

DiMuonValidation::~DiMuonValidation ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 205 of file DiMuonValidation.cc.

References HLT_2024v12_cff::delta_eta, HLT_2024v12_cff::DeltaEta, eBeam_, DiLeptonHelp::PlotsVsDiLeptonRegion::fillTH1Plots(), DiLeptonHelp::PlotsVsDiLeptonRegion::fillTH2Plots(), iEvent, InvMassInEtaBins, InvMassVsCosThetaCSInEtaBins, InvMassVsPhiMinusInEtaBins, InvMassVsPhiPlusInEtaBins, mu_mass2_, BPhysicsValidation_cfi::muminus, BPhysicsValidation_cfi::muplus, pair_etamaxneg_, pair_etamaxpos_, pair_etaminneg_, pair_mass_max_, pair_mass_min_, funct::pow(), dttmaxenums::R, mathSSE::sqrt(), edm::swap(), th1f_mass, th2d_mass_variables_, th3d_mass_vs_eta_phi_minus_, th3d_mass_vs_eta_phi_plus_, and theTrackCollectionToken_.

                                                                                   {
   using namespace edm;
 
   const reco::TrackCollection& tC = iEvent.get(theTrackCollectionToken_);
 
   DiMuonValid::LV LV_mother(0., 0., 0., 0.);
   for (reco::TrackCollection::const_iterator track1 = tC.begin(); track1 != tC.end(); track1++) {
     DiMuonValid::LV LV_track1(track1->px(),
                               track1->py(),
                               track1->pz(),
                               sqrt((track1->p() * track1->p()) + mu_mass2_));  //old 106
 
     for (reco::TrackCollection::const_iterator track2 = track1 + 1; track2 != tC.end(); track2++) {
       if (track1->charge() == track2->charge()) {
         continue;
       }  // only reconstruct opposite charge pair
 
       DiMuonValid::LV LV_track2(
           track2->px(), track2->py(), track2->pz(), sqrt((track2->p() * track2->p()) + mu_mass2_));
 
       LV_mother = LV_track1 + LV_track2;
       double mother_mass = LV_mother.M();
       th1f_mass->Fill(mother_mass);
       double mother_pt = LV_mother.Pt();
 
       int charge1 = track1->charge();
       double etaMu1 = track1->eta();
       double phiMu1 = track1->phi();
       double ptMu1 = track1->pt();
 
       int charge2 = track2->charge();
       double etaMu2 = track2->eta();
       double phiMu2 = track2->phi();
       double ptMu2 = track2->pt();
 
       if (charge1 < 0) {  // use Mu+ for charge1, Mu- for charge2
         std::swap(charge1, charge2);
         std::swap(etaMu1, etaMu2);
         std::swap(phiMu1, phiMu2);
         std::swap(ptMu1, ptMu2);
       }
 
       const auto& tplus = track1->charge() > 0 ? track1 : track2;
       const auto& tminus = track1->charge() < 0 ? track1 : track2;
       TLorentzVector p4_tplus(tplus->px(), tplus->py(), tplus->pz(), sqrt((tplus->p() * tplus->p()) + mu_mass2_));
       TLorentzVector p4_tminus(tminus->px(), tminus->py(), tminus->pz(), sqrt((tminus->p() * tminus->p()) + mu_mass2_));
       std::pair<TLorentzVector, TLorentzVector> tktk_p4 = std::make_pair(p4_tplus, p4_tminus);
       InvMassInEtaBins.fillTH1Plots(mother_mass, tktk_p4);
       InvMassVsPhiPlusInEtaBins.fillTH2Plots(mother_mass, phiMu1, tktk_p4);
       InvMassVsPhiMinusInEtaBins.fillTH2Plots(mother_mass, phiMu2, tktk_p4);
 
       //eta cut
       if (etaMu1 < pair_etaminpos_ || etaMu1 > pair_etamaxpos_ || etaMu2 < pair_etaminneg_ ||
           etaMu2 > pair_etamaxneg_) {
         continue;
       }
 
       double delta_eta = etaMu1 - etaMu2;
 
       double muplus = 1.0 / sqrt(2.0) * (LV_track1.E() + LV_track1.Z());
       double muminus = 1.0 / sqrt(2.0) * (LV_track1.E() - LV_track1.Z());
       double mubarplus = 1.0 / sqrt(2.0) * (LV_track2.E() + LV_track2.Z());
       double mubarminus = 1.0 / sqrt(2.0) * (LV_track2.E() - LV_track2.Z());
       //double costheta = 2.0 / Q.mag() / sqrt(pow(Q.mag(), 2) + pow(Q.Pt(), 2)) * (muplus * mubarminus - muminus * mubarplus);
       double costhetaCS = 2.0 / LV_mother.mag() / sqrt(pow(LV_mother.mag(), 2) + pow(LV_mother.Pt(), 2)) *
                           (muplus * mubarminus - muminus * mubarplus);
 
       InvMassVsCosThetaCSInEtaBins.fillTH2Plots(mother_mass, costhetaCS, tktk_p4);
 
       DiMuonValid::LV Pbeam(0., 0., eBeam_, eBeam_);
       auto R = Pbeam.Vect().Cross(LV_mother.Vect());
       auto Runit = R.Unit();
       auto Qt = LV_mother.Vect();
       Qt.SetZ(0);
       auto Qtunit = Qt.Unit();
       DiMuonValid::LV D(LV_track1 - LV_track2);
       auto Dt = D.Vect();
       Dt.SetZ(0);
       double tanphi =
           sqrt(pow(LV_mother.mag(), 2) + pow(LV_mother.Pt(), 2)) / LV_mother.mag() * Dt.Dot(Runit) / Dt.Dot(Qtunit);
       double phiCS = atan(tanphi);
 
       if (mother_mass > pair_mass_min_ && mother_mass < pair_mass_max_) {
         th2d_mass_variables_[Variable::CosThetaCS]->Fill(mother_mass, costhetaCS, 1);
         th2d_mass_variables_[Variable::DeltaEta]->Fill(mother_mass, delta_eta, 1);
         th2d_mass_variables_[Variable::EtaMinus]->Fill(mother_mass, etaMu2, 1);
         th2d_mass_variables_[Variable::EtaPlus]->Fill(mother_mass, etaMu1, 1);
         th2d_mass_variables_[Variable::PhiCS]->Fill(mother_mass, phiCS, 1);
         th2d_mass_variables_[Variable::PhiMinus]->Fill(mother_mass, phiMu2, 1);
         th2d_mass_variables_[Variable::PhiPlus]->Fill(mother_mass, phiMu1, 1);
         th2d_mass_variables_[Variable::Pt]->Fill(mother_mass, mother_pt, 1);
 
         th3d_mass_vs_eta_phi_plus_->Fill(mother_mass, etaMu1, phiMu1);
         th3d_mass_vs_eta_phi_minus_->Fill(mother_mass, etaMu2, phiMu2);
       }
     }
   }
 }

◆ beginJob()

void DiMuonValidation::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 305 of file DiMuonValidation.cc.

References DiLeptonHelp::PlotsVsDiLeptonRegion::bookSet(), compressionSettings_, dqm-mbProfile::format, compareTotals::fs, mps_fire::i, InvMassInEtaBins, InvMassVsCosThetaCSInEtaBins, InvMassVsPhiMinusInEtaBins, InvMassVsPhiPlusInEtaBins, pair_mass_max_, pair_mass_min_, pair_mass_nbins_, AlCaHLTBitMon_QueryRunRegistry::string, th1f_mass, th2d_mass_variables_, th3d_mass_vs_eta_phi_minus_, th3d_mass_vs_eta_phi_plus_, variables_bins_number_, variables_max_, variables_min_, variables_name_, and variables_title_.

                                 {
   edm::Service<TFileService> fs;
   if (compressionSettings_ > 0) {
     fs->file().SetCompressionSettings(compressionSettings_);
   }
 
   th1f_mass = fs->make<TH1F>("hMass", "mass;m_{#mu#mu} [GeV];events", 200, 0., 200.);
 
   for (int i = 0; i < Variable::VarNumber; i++) {
     std::string th2d_name = fmt::sprintf("th2d_mass_%s", variables_name_[i].c_str());
     th2d_mass_variables_[i] =
         fs->make<TH2D>(th2d_name.c_str(),
                        fmt::format("{};M_{{#mu^{{-}}#mu^{{+}}}} [GeV];{}", th2d_name, variables_title_[i]).c_str(),
                        pair_mass_nbins_,
                        pair_mass_min_,
                        pair_mass_max_,
                        variables_bins_number_[i],
                        variables_min_[i],
                        variables_max_[i]);
   }
 
   // 3D histogram for eta/phi map (mu+)
   th3d_mass_vs_eta_phi_plus_ =
       fs->make<TH3D>("th3d_mass_vs_eta_phi_plus",
                      "th3d_mass_vs_eta_phi_plus;M_{#mu^{-}#mu^{+}} [GeV];#mu^{+} #eta;#mu^{+} #phi [rad]",
                      pair_mass_nbins_,
                      pair_mass_min_,
                      pair_mass_max_,
                      variables_bins_number_[Variable::EtaPlus],
                      variables_min_[Variable::EtaPlus],
                      variables_max_[Variable::EtaPlus],
                      variables_bins_number_[Variable::PhiPlus],
                      variables_min_[Variable::PhiPlus],
                      variables_max_[Variable::PhiPlus]);
 
   // 3D histogram for eta/phi map (mu+)
   th3d_mass_vs_eta_phi_minus_ =
       fs->make<TH3D>("th3d_mass_vs_eta_phi_minus",
                      "th3d_mass_vs_eta_phi_minus;M_{#mu^{-}#mu^{+}} [GeV];#mu^{-} #eta;#mu^{-} #phi [rad]",
                      pair_mass_nbins_,
                      pair_mass_min_,
                      pair_mass_max_,
                      variables_bins_number_[Variable::EtaMinus],
                      variables_min_[Variable::EtaMinus],
                      variables_max_[Variable::EtaMinus],
                      variables_bins_number_[Variable::PhiMinus],
                      variables_min_[Variable::PhiMinus],
                      variables_max_[Variable::PhiMinus]);
 
   // Z-> mm mass in eta bins
   TFileDirectory dirResMassEta = fs->mkdir("TkTkMassInEtaBins");
   InvMassInEtaBins.bookSet(dirResMassEta, th1f_mass);
 
   // Z->mm mass vs cos phi Collins-Soper in eta bins
   TFileDirectory dirResMassVsCosThetaCSInEta = fs->mkdir("TkTkMassVsCosThetaCSInEtaBins");
   InvMassVsCosThetaCSInEtaBins.bookSet(dirResMassVsCosThetaCSInEta, th2d_mass_variables_[Variable::CosThetaCS]);
 
   // Z->mm mass vs phi minus in eta bins
   TFileDirectory dirResMassVsPhiMinusInEta = fs->mkdir("TkTkMassVsPhiMinusInEtaBins");
   InvMassVsPhiMinusInEtaBins.bookSet(dirResMassVsPhiMinusInEta, th2d_mass_variables_[Variable::PhiMinus]);
 
   // Z->mm mass vs phi plus in eta bins
   TFileDirectory dirResMassVsPhiPlusInEta = fs->mkdir("TkTkMassVsPhiPlusInEtaBins");
   InvMassVsPhiPlusInEtaBins.bookSet(dirResMassVsPhiPlusInEta, th2d_mass_variables_[Variable::PhiPlus]);
 }

◆ fillDescriptions()

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

static

Definition at line 372 of file DiMuonValidation.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, M_PI, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                   {
   edm::ParameterSetDescription desc;
   desc.setComment("Validates alignment payloads by evaluating bias in Z->mm mass distributions");
   desc.addUntracked<int>("compressionSettings", -1);
 
   desc.add<double>("eBeam", 3500.)->setComment("beam energy in GeV");
   desc.add<std::string>("TkTag", "ALCARECOTkAlZMuMu");
 
   desc.add<double>("Pair_mass_min", 60.);
   desc.add<double>("Pair_mass_max", 120.);
   desc.add<int>("Pair_mass_nbins", 120);
   desc.add<double>("Pair_etaminpos", -2.4);
   desc.add<double>("Pair_etamaxpos", 2.4);
   desc.add<double>("Pair_etaminneg", -2.4);
   desc.add<double>("Pair_etamaxneg", 2.4);
 
   desc.add<double>("Variable_CosThetaCS_xmin", -1.);
   desc.add<double>("Variable_CosThetaCS_xmax", 1.);
   desc.add<int>("Variable_CosThetaCS_nbins", 20);
 
   desc.add<double>("Variable_DeltaEta_xmin", -4.8);
   desc.add<double>("Variable_DeltaEta_xmax", 4.8);
   desc.add<int>("Variable_DeltaEta_nbins", 20);
 
   desc.add<double>("Variable_EtaMinus_xmin", -2.4);
   desc.add<double>("Variable_EtaMinus_xmax", 2.4);
   desc.add<int>("Variable_EtaMinus_nbins", 12);
 
   desc.add<double>("Variable_EtaPlus_xmin", -2.4);
   desc.add<double>("Variable_EtaPlus_xmax", 2.4);
   desc.add<int>("Variable_EtaPlus_nbins", 12);
 
   desc.add<double>("Variable_PhiCS_xmin", -M_PI / 2);
   desc.add<double>("Variable_PhiCS_xmax", M_PI / 2);
   desc.add<int>("Variable_PhiCS_nbins", 20);
 
   desc.add<double>("Variable_PhiMinus_xmin", -M_PI);
   desc.add<double>("Variable_PhiMinus_xmax", M_PI);
   desc.add<int>("Variable_PhiMinus_nbins", 16);
 
   desc.add<double>("Variable_PhiPlus_xmin", -M_PI);
   desc.add<double>("Variable_PhiPlus_xmax", M_PI);
   desc.add<int>("Variable_PhiPlus_nbins", 16);
 
   desc.add<double>("Variable_PairPt_xmin", 0.);
   desc.add<double>("Variable_PairPt_xmax", 100.);
   desc.add<int>("Variable_PairPt_nbins", 100);
 
   descriptions.addWithDefaultLabel(desc);
 }

Member Data Documentation

◆ compressionSettings_

int DiMuonValidation::compressionSettings_

private

Definition at line 128 of file DiMuonValidation.cc.

Referenced by beginJob().

◆ eBeam_

float DiMuonValidation::eBeam_

private

Definition at line 127 of file DiMuonValidation.cc.

Referenced by analyze().

◆ InvMassInEtaBins

DiLeptonHelp::PlotsVsDiLeptonRegion DiMuonValidation::InvMassInEtaBins = DiLeptonHelp::PlotsVsDiLeptonRegion(1.5)

private

Definition at line 194 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ InvMassVsCosThetaCSInEtaBins

DiLeptonHelp::PlotsVsDiLeptonRegion DiMuonValidation::InvMassVsCosThetaCSInEtaBins = DiLeptonHelp::PlotsVsDiLeptonRegion(1.5)

private

Definition at line 197 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ InvMassVsPhiMinusInEtaBins

DiLeptonHelp::PlotsVsDiLeptonRegion DiMuonValidation::InvMassVsPhiMinusInEtaBins = DiLeptonHelp::PlotsVsDiLeptonRegion(1.5)

private

Definition at line 196 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ InvMassVsPhiPlusInEtaBins

DiLeptonHelp::PlotsVsDiLeptonRegion DiMuonValidation::InvMassVsPhiPlusInEtaBins = DiLeptonHelp::PlotsVsDiLeptonRegion(1.5)

private

Definition at line 195 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ mu_mass2_

constexpr double DiMuonValidation::mu_mass2_ = 0.105658 * 0.105658

static

Definition at line 120 of file DiMuonValidation.cc.

Referenced by analyze().

◆ pair_etamaxneg_

double DiMuonValidation::pair_etamaxneg_

private

Definition at line 137 of file DiMuonValidation.cc.

Referenced by analyze().

◆ pair_etamaxpos_

double DiMuonValidation::pair_etamaxpos_

private

Definition at line 135 of file DiMuonValidation.cc.

Referenced by analyze().

◆ pair_etaminneg_

double DiMuonValidation::pair_etaminneg_

private

Definition at line 136 of file DiMuonValidation.cc.

Referenced by analyze().

◆ pair_etaminpos_

double DiMuonValidation::pair_etaminpos_

private

Definition at line 134 of file DiMuonValidation.cc.

◆ pair_mass_max_

double DiMuonValidation::pair_mass_max_

private

Definition at line 132 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ pair_mass_min_

double DiMuonValidation::pair_mass_min_

private

Definition at line 131 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ pair_mass_nbins_

int DiMuonValidation::pair_mass_nbins_

private

Definition at line 133 of file DiMuonValidation.cc.

Referenced by beginJob().

◆ th1f_mass

TH1F* DiMuonValidation::th1f_mass

private

Definition at line 173 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ th2d_mass_variables_

TH2D* DiMuonValidation::th2d_mass_variables_[Variable::VarNumber]

private

Definition at line 174 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ th3d_mass_vs_eta_phi_minus_

TH3D* DiMuonValidation::th3d_mass_vs_eta_phi_minus_

private

Definition at line 176 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ th3d_mass_vs_eta_phi_plus_

TH3D* DiMuonValidation::th3d_mass_vs_eta_phi_plus_

private

Definition at line 175 of file DiMuonValidation.cc.

Referenced by analyze(), and beginJob().

◆ theTrackCollectionToken_

edm::EDGetTokenT<reco::TrackCollection> DiMuonValidation::theTrackCollectionToken_

private

Definition at line 171 of file DiMuonValidation.cc.

Referenced by analyze(), and DiMuonValidation().

◆ TkTag_

std::string DiMuonValidation::TkTag_

private

Definition at line 129 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_CosThetaCS_nbins_

int DiMuonValidation::variable_CosThetaCS_nbins_

private

Definition at line 141 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_CosThetaCS_xmax_

double DiMuonValidation::variable_CosThetaCS_xmax_

private

Definition at line 140 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_CosThetaCS_xmin_

double DiMuonValidation::variable_CosThetaCS_xmin_

private

Definition at line 139 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_DeltaEta_nbins_

int DiMuonValidation::variable_DeltaEta_nbins_

private

Definition at line 145 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_DeltaEta_xmax_

double DiMuonValidation::variable_DeltaEta_xmax_

private

Definition at line 144 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_DeltaEta_xmin_

double DiMuonValidation::variable_DeltaEta_xmin_

private

Definition at line 143 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_EtaMinus_nbins_

int DiMuonValidation::variable_EtaMinus_nbins_

private

Definition at line 149 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_EtaMinus_xmax_

double DiMuonValidation::variable_EtaMinus_xmax_

private

Definition at line 148 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_EtaMinus_xmin_

double DiMuonValidation::variable_EtaMinus_xmin_

private

Definition at line 147 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_EtaPlus_nbins_

int DiMuonValidation::variable_EtaPlus_nbins_

private

Definition at line 153 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_EtaPlus_xmax_

double DiMuonValidation::variable_EtaPlus_xmax_

private

Definition at line 152 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_EtaPlus_xmin_

double DiMuonValidation::variable_EtaPlus_xmin_

private

Definition at line 151 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PairPt_nbins_

int DiMuonValidation::variable_PairPt_nbins_

private

Definition at line 169 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PairPt_xmax_

double DiMuonValidation::variable_PairPt_xmax_

private

Definition at line 168 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PairPt_xmin_

double DiMuonValidation::variable_PairPt_xmin_

private

Definition at line 167 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PhiCS_nbins_

int DiMuonValidation::variable_PhiCS_nbins_

private

Definition at line 157 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PhiCS_xmax_

double DiMuonValidation::variable_PhiCS_xmax_

private

Definition at line 156 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PhiCS_xmin_

double DiMuonValidation::variable_PhiCS_xmin_

private

Definition at line 155 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PhiMinus_nbins_

int DiMuonValidation::variable_PhiMinus_nbins_

private

Definition at line 161 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PhiMinus_xmax_

double DiMuonValidation::variable_PhiMinus_xmax_

private

Definition at line 160 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PhiMinus_xmin_

double DiMuonValidation::variable_PhiMinus_xmin_

private

Definition at line 159 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PhiPlus_nbins_

int DiMuonValidation::variable_PhiPlus_nbins_

private

Definition at line 165 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PhiPlus_xmax_

double DiMuonValidation::variable_PhiPlus_xmax_

private

Definition at line 164 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variable_PhiPlus_xmin_

double DiMuonValidation::variable_PhiPlus_xmin_

private

Definition at line 163 of file DiMuonValidation.cc.

Referenced by DiMuonValidation().

◆ variables_bins_number_

int DiMuonValidation::variables_bins_number_[Variable::VarNumber]

private

Definition at line 190 of file DiMuonValidation.cc.

Referenced by beginJob(), and DiMuonValidation().

◆ variables_max_

double DiMuonValidation::variables_max_[Variable::VarNumber]

private

Definition at line 192 of file DiMuonValidation.cc.

Referenced by beginJob(), and DiMuonValidation().

◆ variables_min_

double DiMuonValidation::variables_min_[Variable::VarNumber]

private

Definition at line 191 of file DiMuonValidation.cc.

Referenced by beginJob(), and DiMuonValidation().

◆ variables_name_

std::string DiMuonValidation::variables_name_[Variable::VarNumber]

private

Initial value:

= {

"CosThetaCS", "DeltaEta", "EtaMinus", "EtaPlus", "PhiCS", "PhiMinus", "PhiPlus", "Pt"}

Definition at line 178 of file DiMuonValidation.cc.

Referenced by beginJob().

◆ variables_title_

std::string DiMuonValidation::variables_title_[Variable::VarNumber]

private

Initial value:

= {"Cos#theta_{CS}",
                                                       "#Delta#eta(#mu^{-},#mu^{+})",
                                                       "#mu^{-} #eta",
                                                       "#mu^{+} #eta",
                                                       "#phi_{CS} [rad]",
                                                       "#mu^{-} #phi [rad]",
                                                       "#mu^{+} #phi [rad]",
                                                       "p_{T} [GeV]"}

Definition at line 181 of file DiMuonValidation.cc.

Referenced by beginJob().

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ Variable

Constructor & Destructor Documentation

◆ DiMuonValidation()

◆ ~DiMuonValidation()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ fillDescriptions()

Member Data Documentation

◆ compressionSettings_

◆ eBeam_

◆ InvMassInEtaBins

◆ InvMassVsCosThetaCSInEtaBins

◆ InvMassVsPhiMinusInEtaBins

◆ InvMassVsPhiPlusInEtaBins

◆ mu_mass2_

◆ pair_etamaxneg_

◆ pair_etamaxpos_

◆ pair_etaminneg_

◆ pair_etaminpos_

◆ pair_mass_max_

◆ pair_mass_min_

◆ pair_mass_nbins_

◆ th1f_mass

◆ th2d_mass_variables_

◆ th3d_mass_vs_eta_phi_minus_

◆ th3d_mass_vs_eta_phi_plus_

◆ theTrackCollectionToken_

◆ TkTag_

◆ variable_CosThetaCS_nbins_

◆ variable_CosThetaCS_xmax_

◆ variable_CosThetaCS_xmin_

◆ variable_DeltaEta_nbins_

◆ variable_DeltaEta_xmax_

◆ variable_DeltaEta_xmin_

◆ variable_EtaMinus_nbins_

◆ variable_EtaMinus_xmax_

◆ variable_EtaMinus_xmin_

◆ variable_EtaPlus_nbins_

◆ variable_EtaPlus_xmax_

◆ variable_EtaPlus_xmin_

◆ variable_PairPt_nbins_

◆ variable_PairPt_xmax_

◆ variable_PairPt_xmin_

◆ variable_PhiCS_nbins_

◆ variable_PhiCS_xmax_

◆ variable_PhiCS_xmin_

◆ variable_PhiMinus_nbins_

◆ variable_PhiMinus_xmax_

◆ variable_PhiMinus_xmin_

◆ variable_PhiPlus_nbins_

◆ variable_PhiPlus_xmax_

◆ variable_PhiPlus_xmin_

◆ variables_bins_number_

◆ variables_max_

◆ variables_min_

◆ variables_name_

◆ variables_title_