ATLAS_2012_I1204447 Class Reference
Inheritance diagram for ATLAS_2012_I1204447:
![]()
Collaboration diagram for ATLAS_2012_I1204447:
![]()
Detailed DescriptionDefinition at line 14 of file ATLAS_2012_I1204447.cc. Constructor & Destructor Documentation
Constructor. Definition at line 18 of file ATLAS_2012_I1204447.cc. : Analysis("ATLAS_2012_I1204447") { } Member Function Documentation
Untemplated function to do the work... Definition at line 33 of file ProjectionApplier.cc. References ProjectionApplier::_allowProjReg, ProjectionApplier::getProjHandler(), ProjectionApplier::name(), Projection::name(), and ProjectionHandler::registerProjection(). Referenced by ProjectionApplier::addProjection(). { if (!_allowProjReg) { cerr << "Trying to register projection '" << proj.name() << "' before init phase in '" << this->name() << "'." << endl; exit(2); } const Projection& reg = getProjHandler().registerProjection(*this, proj, name); return reg; }
Register a data object in the histogram system.
Definition at line 721 of file Analysis.cc. References Analysis::_analysisobjects. Referenced by Analysis::bookHisto1D(), Analysis::bookHisto2D(), Analysis::bookProfile1D(), Analysis::bookProfile2D(), Analysis::bookScatter2D(), and H1_1995_S3167097::finalize(). { _analysisobjects.push_back(ao); }
Register a contained projection. The type of the argument is used to instantiate a new projection internally: this new object is applied to events rather than the argument object. Hence you are advised to only use locally-scoped Projection objects in your Projection and Analysis constructors, and to avoid polymorphism (e.g. handling Definition at line 116 of file ProjectionApplier.hh. References ProjectionApplier::_addProjection(). Referenced by FastJets::_init1(), FastJets::_init2(), FastJets::_init3(), VetoedFinalState::addVetoOnThisFinalState(), BeamThrust::BeamThrust(), CDF_2009_S8057893::CDF_2009_S8057893::init(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), DressedLeptons::DressedLeptons(), FinalState::FinalState(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), HeavyHadrons::HeavyHadrons(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), CMS_2010_S8547297::init(), CMS_2010_S8656010::init(), ATLAS_2010_S8894728::init(), ALICE_2012_I1181770::init(), ATLAS_2011_S8994773::init(), CMS_2011_S8950903::init(), CMS_2012_PAS_QCD_11_010::init(), CMS_2011_S8968497::init(), CMS_2012_I1193338::init(), CMS_2012_I1184941::init(), CMSTOTEM_2014_I1294140::init(), ATLAS_2014_I1298811::init(), ATLAS_2011_I894867::init(), ALICE_2011_S8909580::init(), LHCB_2013_I1208105::init(), CMS_2011_S8941262::init(), CMS_2011_S8973270::init(), CMS_2012_I1090423::init(), ATLAS_2010_CONF_2010_049::init(), CDF_2007_S7057202::init(), CMS_2011_S9086218::init(), TOTEM_2012_I1115294::init(), CMS_2011_S9215166::init(), CMS_2012_I1087342::init(), ATLAS_2014_I1282441::init(), ALICE_2011_S8945144::init(), ATLAS_2012_I1091481::init(), MC_MUONS::init(), MC_TAUS::init(), ATLAS_2010_S8591806::init(), ATLAS_2011_S9002537::init(), CMS_2011_I954992::init(), CMS_2011_S8957746::init(), CMS_2011_S8978280::init(), TOTEM_2012_002::init(), CMS_2011_S9088458::init(), CMS_2013_I1265659::init(), D0_2011_I895662::init(), ATLAS_2011_I930220::init(), LHCF_2012_I1115479::init(), MC_ELECTRONS::init(), CMS_2013_I1273574::init(), D0_2010_S8570965::init(), ATLAS_2011_I925932::init(), MC_JETTAGS::init(), CDF_2012_NOTE10874::init(), ATLAS_2012_I1124167::init(), STAR_2006_S6500200::init(), STAR_2008_S7993412::init(), CMS_2013_I1256943::init(), CMS_QCD_10_024::init(), UA5_1987_S1640666::init(), CDF_1997_S3541940::init(), MC_HINC::init(), MC_JETS::init(), MC_KTSPLITTINGS::init(), CDF_2005_S6080774::init(), CDF_2006_S6450792::init(), SFM_1984_S1178091::init(), ARGUS_1993_S2789213::init(), ATLAS_2011_S9128077::init(), CMS_2013_I1261026::init(), CMS_2013_I1258128::init(), CMS_2014_I1298810::init(), D0_2000_S4480767::init(), BABAR_2007_S7266081::init(), BELLE_2008_I786560::init(), CDF_1993_S2742446::init(), CDF_2000_S4155203::init(), MC_WWINC::init(), CDF_2008_S7782535::init(), MC_ZJETS::init(), CDF_2008_S8093652::init(), MC_ZZINC::init(), STAR_2006_S6870392::init(), UA5_1982_S875503::init(), UA5_1989_S1926373::init(), ATLAS_2013_I1190187::init(), CMS_2013_I1209721::init(), D0_2008_S6879055::init(), ATLAS_2012_I1204447::init(), ARGUS_1993_S2669951::init(), LHCB_2012_I1208102::init(), ATLAS_2011_I944826::init(), MC_ZINC::init(), ALEPH_1991_S2435284::init(), MC_WINC::init(), MC_WJETS::init(), MC_WWKTSPLITTINGS::init(), MC_ZZJETS::init(), ATLAS_2012_I1183818::init(), ATLAS_2011_S9035664::init(), ATLAS_2012_I1204784::init(), CMS_2012_I1102908::init(), ATLAS_2013_I1219109::init(), CMS_2013_I1272853::init(), ATLAS_2013_I1243871::init(), ATLAS_2014_I1268975::init(), D0_2000_I499943::init(), ATLAS_2014_I1325553::init(), E735_1998_S3905616::init(), MC_ZZKTSPLITTINGS::init(), MC_DIJET::init(), MC_DIPHOTON::init(), ALEPH_2002_S4823664::init(), DELPHI_1999_S3960137::init(), ALEPH_1996_S3196992::init(), UA5_1986_S1583476::init(), EXAMPLE_CUTS::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), MC_ZKTSPLITTINGS::init(), MC_HJETS::init(), MC_HKTSPLITTINGS::init(), CDF_2005_S6217184::init(), MC_PHOTONINC::init(), MC_WKTSPLITTINGS::init(), ATLAS_2011_S8924791::init(), OPAL_1995_S3198391::init(), OPAL_1996_S3257789::init(), OPAL_1997_S3608263::init(), OPAL_1998_S3702294::init(), OPAL_1998_S3749908::init(), OPAL_2000_S4418603::init(), ATLAS_2010_S8817804::init(), PDG_TAUS::init(), D0_1996_S3324664::init(), D0_2007_S7075677::init(), D0_2009_S8202443::init(), D0_2010_S8821313::init(), ATLAS_2011_I945498::init(), MC_PHOTONKTSPLITTINGS::init(), MC_PHOTONJETS::init(), MC_WWJETS::init(), CMS_2012_I1107658::init(), D0_2008_S7554427::init(), D0_2008_S7863608::init(), D0_2010_S8671338::init(), MC_IDENTIFIED::init(), H1_1995_S3167097::init(), CDF_1994_S2952106::init(), ATLAS_2011_I954993::init(), ATLAS_2011_S9212353::init(), ATLAS_2014_I1306615::init(), ATLAS_2014_I1307756::init(), ATLAS_2012_CONF_2012_104::init(), ATLAS_2012_CONF_2012_105::init(), UA1_1990_S2044935::init(), ZEUS_2001_S4815815::init(), MC_LEADJETUE::init(), LHCB_2011_I919315::init(), MC_TTBAR::init(), ATLAS_2012_I1119557::init(), CDF_2008_S7828950::init(), MC_GENERIC::init(), ATLAS_2012_I1180197::init(), MC_PHOTONJETUE::init(), ATLAS_2012_I1188891::init(), CDF_2008_S7540469::init(), D0_2001_S4674421::init(), ATLAS_2014_I1300647::init(), D0_2008_S7837160::init(), ATLAS_2011_S9131140::init(), CDF_2008_S8095620::init(), DELPHI_1995_S3137023::init(), ATLAS_2012_CONF_2012_109::init(), ATLAS_2012_I1126136::init(), ATLAS_2011_S8983313::init(), ATLAS_2013_I1217863_Z::init(), STAR_2009_UE_HELEN::init(), OPAL_1997_S3396100::init(), STAR_2006_S6860818::init(), ATLAS_2011_CONF_2011_098::init(), JADE_1998_S3612880::init(), CMS_2012_I941555::init(), ATLAS_2012_I1095236::init(), JADE_OPAL_2000_S4300807::init(), D0_1996_S3214044::init(), ATLAS_2012_CONF_2012_103::init(), CDF_1996_S3108457::init(), EXAMPLE::init(), ATLAS_2011_I926145::init(), ALICE_2010_S8624100::init(), MC_QCD_PARTONS::init(), ATLAS_2012_I1186556::init(), ATLAS_2013_I1217863_W::init(), CMS_2014_I1303894::init(), ALEPH_2004_S5765862::init(), ALICE_2010_S8706239::init(), MC_HFJETS::init(), ATLAS_2014_I1288706::init(), ATLAS_2012_I1112263::init(), ATLAS_2011_S9212183::init(), ATLAS_2012_I1190891::init(), ATLAS_2012_CONF_2012_001::init(), CMS_2013_I1218372::init(), ATLAS_2012_I1125575::init(), ATLAS_2012_I1125961::init(), LHCB_2013_I1218996::init(), ALICE_2010_S8625980::init(), CDF_2009_NOTE_9936::init(), CDF_1996_S3349578::init(), MC_SUSY::init(), ATLAS_2011_S8971293::init(), CDF_2009_S8383952::init(), CDF_2009_S8436959::init(), CDF_1998_S3618439::init(), ATLAS_2013_I1263495::init(), D0_2009_S8320160::init(), CDF_2000_S4266730::init(), D0_2009_S8349509::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), D0_2006_S6438750::init(), ALEPH_2001_S4656318::init(), SLD_2002_S4869273::init(), ATLAS_2010_S8919674::init(), ATLAS_2012_I1082936::init(), D0_2010_S8566488::init(), DELPHI_2002_069_CONF_603::init(), LHCB_2014_I1281685::init(), LHCB_2011_I917009::init(), UA5_1988_S1867512::init(), SLD_1996_S3398250::init(), MC_PHOTONS::init(), ATLAS_2011_I921594::init(), CDF_1996_S3418421::init(), LHCB_2012_I1119400::init(), ATLAS_2011_S9108483::init(), MC_WPOL::init(), CDF_2006_S6653332::init(), CDF_2008_S7541902::init(), ATLAS_2014_I1306294::init(), ATLAS_2014_I1304688::init(), D0_2008_S7662670::init(), ATLAS_2013_I1230812::init(), ATLAS_2012_I1199269::init(), ATLAS_2014_I1307243::init(), DELPHI_2000_S4328825::init(), ATLAS_2012_I943401::init(), OPAL_2002_S5361494::init(), TASSO_1990_S2148048::init(), ATLAS_2012_I946427::init(), ATLAS_2012_I1117704::init(), ATLAS_2011_S9019561::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2011_S9225137::init(), ATLAS_2012_I1082009::init(), ATLAS_2012_CONF_2012_153::init(), ATLAS_2012_I1084540::init(), CMS_2013_I1224539_WJET::init(), CMS_2013_I1224539_DIJET::init(), ATLAS_2012_I1083318::init(), CDF_2009_S8233977::init(), ATLAS_2013_I1217867::init(), CMS_2013_I1224539_ZJET::init(), CDF_2010_S8591881_DY::init(), CDF_2010_S8591881_QCD::init(), D0_2004_S5992206::init(), ATLAS_2010_S8914702::init(), ATLAS_2014_I1312627::init(), ATLAS_2010_S8918562::init(), ATLAS_2011_S9120807::init(), CDF_2001_S4751469::init(), D0_2008_S7719523::init(), ATLAS_2011_S9041966::init(), ATLAS_2012_I1094568::init(), DELPHI_1996_S3430090::init(), ALEPH_1996_S3486095::init(), ALEPH_1999_S4193598::init(), ATLAS_2012_I1093738::init(), OPAL_2004_S6132243::init(), OPAL_1994_S2927284::init(), BELLE_2013_I1216515::init(), STAR_2008_S7869363::init(), ATLAS_2011_S9126244::init(), BABAR_2007_S6895344::init(), MC_VH2BB::init(), BABAR_2005_S6181155::init(), BELLE_2001_S4598261::init(), ATLAS_2011_I919017::init(), OPAL_2001_S4553896::init(), BABAR_2013_I1238276::init(), DELPHI_2003_WUD_03_11::init(), ATLAS_2012_I1203852::init(), CDF_2004_S5839831::init(), BABAR_2003_I593379::init(), ATLAS_2012_I1094061::init(), ARGUS_1993_S2653028::init(), CLEO_2004_S5809304::init(), OPAL_1998_S3780481::init(), ATLAS_2012_I1093734::init(), OPAL_1993_S2692198::init(), ATLAS_2014_I1279489::init(), ATLAS_2012_I1094564::init(), H1_1994_S2919893::init(), H1_2000_S4129130::init(), SLD_2004_S5693039::init(), SLD_1999_S3743934::init(), PDG_HADRON_MULTIPLICITIES::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), MissingMomentum::MissingMomentum(), NeutralFinalState::NeutralFinalState(), NonHadronicFinalState::NonHadronicFinalState(), ParisiTensor::ParisiTensor(), PrimaryHadrons::PrimaryHadrons(), PromptFinalState::PromptFinalState(), Sphericity::Sphericity(), Spherocity::Spherocity(), TauFinder::TauFinder(), Thrust::Thrust(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), VetoedFinalState::VetoedFinalState(), VisibleFinalState::VisibleFinalState(), WFinder::WFinder(), and ZFinder::ZFinder(). { const Projection& reg = _addProjection(proj, name); const PROJ& rtn = dynamic_cast<const PROJ&>(reg); return rtn; }
List of registered analysis data objects. Definition at line 694 of file Analysis.hh. References Analysis::_analysisobjects. Referenced by Analysis::getAnalysisObject(). { return _analysisobjects; } Perform the per-event analysis.
Implements Analysis. Definition at line 88 of file ATLAS_2012_I1204447.cc. References ATLAS_2012_I1204447::_h_e_n, ATLAS_2012_I1204447::_h_HTjets_all, ATLAS_2012_I1204447::_h_HTlep_all, ATLAS_2012_I1204447::_h_Meff_all, ATLAS_2012_I1204447::_h_MET_all, ATLAS_2012_I1204447::_h_mu_n, ATLAS_2012_I1204447::_h_pt_1_2ltau, ATLAS_2012_I1204447::_h_pt_1_3l, ATLAS_2012_I1204447::_h_pt_2_2ltau, ATLAS_2012_I1204447::_h_pt_2_3l, ATLAS_2012_I1204447::_h_pt_3_2ltau, ATLAS_2012_I1204447::_h_pt_3_3l, ATLAS_2012_I1204447::_h_tau_n, ATLAS_2012_I1204447::_use_fiducial_lepton_efficiency, ParticleBase::abseta(), FourVector::abseta(), Particle::abspid(), ATLAS_2012_I1204447::apply_reco_eff(), ParticleBase::azimuthalAngle(), Rivet::deltaR(), Rivet::Cuts::eta, ATLAS_2012_I1204447::fillEventCountsPerSR(), Particle::genParticle(), ATLAS_2012_I1204447::get_prong_number(), ATLAS_2012_I1204447::get_tau_neutrino_mom(), Rivet::GeV, Particle::hasAncestor(), Rivet::inRange(), ATLAS_2012_I1204447::isonZ(), Rivet::MINUSPI_PLUSPI, Particle::momentum(), Jet::momentum(), MSG_DEBUG, Rivet::PID::MUON, Rivet::particles(), ParticleBase::perp(), Rivet::Cuts::phi, Rivet::Cuts::pT, ParticleBase::pT(), FourMomentum::pT(), ParticleBase::rapidity(), Rivet::sortByPt(), Rivet::PID::TAU, and vetoEvent. { // Muons Particles muon_candidates; const Particles charged_tracks = applyProjection<ChargedFinalState>(event, "CFS").particles(); const Particles visible_particles = applyProjection<VisibleFinalState>(event, "VFS").particles(); foreach (const Particle& mu, applyProjection<IdentifiedFinalState>(event, "muons").particlesByPt()) { // Calculate pTCone30 variable (pT of all tracks within dR<0.3 - pT of muon itself) double pTinCone = -mu.pT(); foreach (const Particle& track, charged_tracks) { if (deltaR(mu.momentum(), track.momentum()) < 0.3) pTinCone += track.pT(); } // Calculate eTCone30 variable (pT of all visible particles within dR<0.3) double eTinCone = 0.; foreach (const Particle& visible_particle, visible_particles) { if (visible_particle.abspid() != PID::MUON && inRange(deltaR(mu.momentum(), visible_particle.momentum()), 0.1, 0.3)) eTinCone += visible_particle.pT(); } // Apply reconstruction efficiency and simulate reco int muon_id = 13; if ( mu.hasAncestor(15) || mu.hasAncestor(-15)) muon_id = 14; const double eff = (_use_fiducial_lepton_efficiency) ? apply_reco_eff(muon_id, mu) : 1.0; const bool keep_muon = rand()/static_cast<double>(RAND_MAX) <= eff; // Keep muon if pTCone30/pT < 0.15 and eTCone30/pT < 0.2 and reconstructed if (keep_muon && pTinCone/mu.pT() <= 0.15 && eTinCone/mu.pT() < 0.2) muon_candidates.push_back(mu); } // Electrons Particles electron_candidates; foreach (const Particle& e, applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt()) { // Neglect electrons in crack regions if (inRange(e.abseta(), 1.37, 1.52)) continue; // Calculate pTCone30 variable (pT of all tracks within dR<0.3 - pT of electron itself) double pTinCone = -e.pT(); foreach (const Particle& track, charged_tracks) { if (deltaR(e.momentum(), track.momentum()) < 0.3) pTinCone += track.pT(); } // Calculate eTCone30 variable (pT of all visible particles (except muons) within dR<0.3) double eTinCone = 0.; foreach (const Particle& visible_particle, visible_particles) { if (visible_particle.abspid() != PID::MUON && inRange(deltaR(e.momentum(), visible_particle.momentum()), 0.1, 0.3)) eTinCone += visible_particle.pT(); } // Apply reconstruction efficiency and simulate reco int elec_id = 11; if (e.hasAncestor(15) || e.hasAncestor(-15)) elec_id = 12; const double eff = (_use_fiducial_lepton_efficiency) ? apply_reco_eff(elec_id, e) : 1.0; const bool keep_elec = rand()/static_cast<double>(RAND_MAX) <= eff; // Keep electron if pTCone30/pT < 0.13 and eTCone30/pT < 0.2 and reconstructed if (keep_elec && pTinCone/e.pT() <= 0.13 && eTinCone/e.pT() < 0.2) electron_candidates.push_back(e); } // Taus /// @todo This could benefit from a tau finder projection Particles tau_candidates; foreach (const Particle& tau, applyProjection<UnstableFinalState>(event, "UFS").particlesByPt()) { // Only pick taus out of all unstable particles if (tau.abspid() != PID::TAU) continue; // Check that tau has decayed into daughter particles /// @todo Huh? Unstable taus with no decay vtx? Can use Particle.isStable()? But why in this situation? if (tau.genParticle()->end_vertex() == 0) continue; // Calculate visible tau pT from pT of tau neutrino in tau decay for pT and |eta| cuts const FourMomentum tau_neutrino_mom = get_tau_neutrino_mom(tau); const FourMomentum visible_tau_mom = tau.momentum() - tau_neutrino_mom; if (visible_tau_mom.pT() < 15*GeV || visible_tau_mom.abseta() > 2.5) continue; // Get prong number (number of tracks) in tau decay and check if tau decays leptonically unsigned int nprong = 0; bool lep_decaying_tau = false; get_prong_number(tau.genParticle(), nprong, lep_decaying_tau); // Apply reconstruction efficiency int tau_id = 15; if (nprong == 1) tau_id = 15; else if (nprong == 3) tau_id = 16; // Get fiducial lepton efficiency simulate reco efficiency const double eff = (_use_fiducial_lepton_efficiency) ? apply_reco_eff(tau_id, tau) : 1.0; const bool keep_tau = rand()/static_cast<double>(RAND_MAX) <= eff; // Keep tau if nprong = 1, it decays hadronically, and it's reconstructed by the detector if ( !lep_decaying_tau && nprong == 1 && keep_tau) tau_candidates.push_back(tau); } // Jets (all anti-kt R=0.4 jets with pT > 25 GeV and eta < 4.9) Jets jet_candidates; foreach (const Jet& jet, applyProjection<FastJets>(event, "AntiKtJets04").jetsByPt(25*GeV)) { if (jet.abseta() < 4.9) jet_candidates.push_back(jet); } // ETmiss Particles vfs_particles = applyProjection<VisibleFinalState>(event, "VFS").particles(); FourMomentum pTmiss; foreach (const Particle& p, vfs_particles) pTmiss -= p.momentum(); double eTmiss = pTmiss.pT()/GeV; //------------------ // Overlap removal // electron - electron Particles electron_candidates_2; for (size_t ie = 0; ie < electron_candidates.size(); ++ie) { const Particle & e = electron_candidates[ie]; bool away = true; // If electron pair within dR < 0.1: remove electron with lower pT for (size_t ie2=0; ie2 < electron_candidates_2.size(); ++ie2) { if ( deltaR( e.momentum(), electron_candidates_2[ie2].momentum()) < 0.1 ) { away = false; break; } } // If isolated keep it if ( away ) electron_candidates_2.push_back( e ); } // jet - electron Jets recon_jets; foreach (const Jet& jet, jet_candidates) { bool away = true; // if jet within dR < 0.2 of electron: remove jet foreach (const Particle& e, electron_candidates_2) { if (deltaR(e.momentum(), jet.momentum()) < 0.2) { away = false; break; } } // jet - tau if (away) { // If jet within dR < 0.2 of tau: remove jet foreach (const Particle& tau, tau_candidates) { if (deltaR(tau.momentum(), jet.momentum()) < 0.2) { away = false; break; } } } // If isolated keep it if ( away ) recon_jets.push_back( jet ); } // electron - jet Particles recon_leptons, recon_e; for (size_t ie = 0; ie < electron_candidates_2.size(); ++ie) { const Particle& e = electron_candidates_2[ie]; // If electron within 0.2 < dR < 0.4 from any jets: remove electron bool away = true; foreach (const Jet& jet, recon_jets) { if (deltaR(e.momentum(), jet.momentum()) < 0.4) { away = false; break; } } // electron - muon // if electron within dR < 0.1 of a muon: remove electron if (away) { foreach (const Particle& mu, muon_candidates) { if (deltaR(mu.momentum(), e.momentum()) < 0.1) { away = false; break; } } } // If isolated keep it if (away) { recon_e += e; recon_leptons += e; } } // tau - electron Particles recon_tau; foreach ( const Particle& tau, tau_candidates ) { bool away = true; // If tau within dR < 0.2 of an electron: remove tau foreach ( const Particle& e, recon_e ) { if (deltaR( tau.momentum(), e.momentum()) < 0.2) { away = false; break; } } // tau - muon // If tau within dR < 0.2 of a muon: remove tau if (away) { foreach (const Particle& mu, muon_candidates) { if (deltaR(tau.momentum(), mu.momentum()) < 0.2) { away = false; break; } } } // If isolated keep it if (away) recon_tau.push_back( tau ); } // Muon - jet isolation Particles recon_mu, trigger_mu; // If muon within dR < 0.4 of a jet, remove muon foreach (const Particle& mu, muon_candidates) { bool away = true; foreach (const Jet& jet, recon_jets) { if ( deltaR( mu.momentum(), jet.momentum()) < 0.4 ) { away = false; break; } } if (away) { recon_mu.push_back( mu ); recon_leptons.push_back( mu ); if (mu.abseta() < 2.4) trigger_mu.push_back( mu ); } } // End overlap removal //------------------ // Jet cleaning if (rand()/static_cast<double>(RAND_MAX) <= 0.42) { foreach (const Jet& jet, recon_jets) { const double eta = jet.rapidity(); const double phi = jet.azimuthalAngle(MINUSPI_PLUSPI); if (jet.pT() > 25*GeV && inRange(eta, -0.1, 1.5) && inRange(phi, -0.9, -0.5)) vetoEvent; } } // Post-isolation event cuts // Require at least 3 charged tracks in event if (charged_tracks.size() < 3) vetoEvent; // And at least one e/mu passing trigger if (!( !recon_e .empty() && recon_e[0] .pT() > 25*GeV) && !( !trigger_mu.empty() && trigger_mu[0].pT() > 25*GeV) ) { MSG_DEBUG("Hardest lepton fails trigger"); vetoEvent; } // And only accept events with at least 2 electrons and muons and at least 3 leptons in total if (recon_mu.size() + recon_e.size() + recon_tau.size() < 3 || recon_leptons.size() < 2) vetoEvent; // Now it's worth getting the event weight const double weight = event.weight(); // Sort leptons by decreasing pT sortByPt(recon_leptons); sortByPt(recon_tau); // Calculate HTlep, fill lepton pT histograms & store chosen combination of 3 leptons double HTlep = 0.; Particles chosen_leptons; if ( recon_leptons.size() > 2 ) { _h_pt_1_3l->fill(recon_leptons[0].perp()/GeV, weight); _h_pt_2_3l->fill(recon_leptons[1].perp()/GeV, weight); _h_pt_3_3l->fill(recon_leptons[2].perp()/GeV, weight); HTlep = (recon_leptons[0].pT() + recon_leptons[1].pT() + recon_leptons[2].pT())/GeV; chosen_leptons.push_back( recon_leptons[0] ); chosen_leptons.push_back( recon_leptons[1] ); chosen_leptons.push_back( recon_leptons[2] ); } else { _h_pt_1_2ltau->fill(recon_leptons[0].perp()/GeV, weight); _h_pt_2_2ltau->fill(recon_leptons[1].perp()/GeV, weight); _h_pt_3_2ltau->fill(recon_tau[0].perp()/GeV, weight); HTlep = (recon_leptons[0].pT() + recon_leptons[1].pT() + recon_tau[0].pT())/GeV ; chosen_leptons.push_back( recon_leptons[0] ); chosen_leptons.push_back( recon_leptons[1] ); chosen_leptons.push_back( recon_tau[0] ); } // Number of prompt e/mu and had taus _h_e_n ->fill(recon_e.size() , weight); _h_mu_n ->fill(recon_mu.size() , weight); _h_tau_n->fill(recon_tau.size(), weight); // Calculate HTjets double HTjets = 0.; foreach ( const Jet & jet, recon_jets ) HTjets += jet.perp()/GeV; // Calculate meff double meff = eTmiss + HTjets; Particles all_leptons; foreach ( const Particle & e , recon_e ) { meff += e.perp()/GeV; all_leptons.push_back( e ); } foreach ( const Particle & mu, recon_mu ) { meff += mu.perp()/GeV; all_leptons.push_back( mu ); } foreach ( const Particle & tau, recon_tau ) { meff += tau.perp()/GeV; all_leptons.push_back( tau ); } // Fill histogram of kinematic variables _h_HTlep_all ->fill(HTlep , weight); _h_HTjets_all->fill(HTjets, weight); _h_MET_all ->fill(eTmiss, weight); _h_Meff_all ->fill(meff , weight); // Determine signal region (3l/2ltau, onZ/offZ) string basic_signal_region; if ( recon_mu.size() + recon_e.size() > 2 ) basic_signal_region += "3l_"; else if ( (recon_mu.size() + recon_e.size() == 2) && (recon_tau.size() > 0)) basic_signal_region += "2ltau_"; // Is there an OSSF pair or a three lepton combination with an invariant mass close to the Z mass int onZ = isonZ(chosen_leptons); if (onZ == 1) basic_signal_region += "onZ"; else if (onZ == 0) basic_signal_region += "offZ"; // Check in which signal regions this event falls and adjust event counters fillEventCountsPerSR(basic_signal_region, onZ, HTlep, eTmiss, HTjets, meff, weight); }
Function giving fiducial lepton efficiency. Definition at line 659 of file ATLAS_2012_I1204447.cc. References Rivet::Cuts::eta, ParticleBase::eta(), Rivet::GeV, Rivet::Cuts::pt, and ParticleBase::pT(). Referenced by ATLAS_2012_I1204447::analyze(). { float pt = p.pT()/GeV; float eta = p.eta(); double eff = 0.; //double err = 0.; if (flavor == 11) { // weight prompt electron -- now including data/MC ID SF in eff. //float rho = 0.820; float p0 = 7.34; float p1 = 0.8977; //float ep0= 0.5 ; float ep1= 0.0087; eff = p1 - p0/pt; //double err0 = ep0/pt; // d(eff)/dp0 //double err1 = ep1; // d(eff)/dp1 //err = sqrt(err0*err0 + err1*err1 - 2*rho*err0*err1); double avgrate = 0.6867; float wz_ele_eta[] = {0.588717,0.603674,0.666135,0.747493,0.762202,0.675051,0.751606,0.745569,0.665333,0.610432,0.592693,}; //float ewz_ele_eta[] ={0.00292902,0.002476,0.00241209,0.00182319,0.00194339,0.00299785,0.00197339,0.00182004,0.00241793,0.00245997,0.00290394,}; int ibin = 3; if (eta >= -2.5 && eta < -2.0) ibin = 0; if (eta >= -2.0 && eta < -1.5) ibin = 1; if (eta >= -1.5 && eta < -1.0) ibin = 2; if (eta >= -1.0 && eta < -0.5) ibin = 3; if (eta >= -0.5 && eta < -0.1) ibin = 4; if (eta >= -0.1 && eta < 0.1) ibin = 5; if (eta >= 0.1 && eta < 0.5) ibin = 6; if (eta >= 0.5 && eta < 1.0) ibin = 7; if (eta >= 1.0 && eta < 1.5) ibin = 8; if (eta >= 1.5 && eta < 2.0) ibin = 9; if (eta >= 2.0 && eta < 2.5) ibin = 10; double eff_eta = wz_ele_eta[ibin]; //double err_eta = ewz_ele_eta[ibin]; eff = (eff*eff_eta)/avgrate; } if (flavor == 12) { // weight electron from tau //float rho = 0.884; float p0 = 6.799; float p1 = 0.842; //float ep0= 0.664; float ep1= 0.016; eff = p1 - p0/pt; //double err0 = ep0/pt; // d(eff)/dp0 //double err1 = ep1; // d(eff)/dp1 //err = sqrt(err0*err0 + err1*err1 - 2*rho*err0*err1); double avgrate = 0.5319; float wz_elet_eta[] = {0.468945,0.465953,0.489545,0.58709,0.59669,0.515829,0.59284,0.575828,0.498181,0.463536,0.481738,}; //float ewz_elet_eta[] ={0.00933795,0.00780868,0.00792679,0.00642083,0.00692652,0.0101568,0.00698452,0.00643524,0.0080002,0.00776238,0.0094699,}; int ibin = 3; if (eta >= -2.5 && eta < -2.0) ibin = 0; if (eta >= -2.0 && eta < -1.5) ibin = 1; if (eta >= -1.5 && eta < -1.0) ibin = 2; if (eta >= -1.0 && eta < -0.5) ibin = 3; if (eta >= -0.5 && eta < -0.1) ibin = 4; if (eta >= -0.1 && eta < 0.1) ibin = 5; if (eta >= 0.1 && eta < 0.5) ibin = 6; if (eta >= 0.5 && eta < 1.0) ibin = 7; if (eta >= 1.0 && eta < 1.5) ibin = 8; if (eta >= 1.5 && eta < 2.0) ibin = 9; if (eta >= 2.0 && eta < 2.5) ibin = 10; double eff_eta = wz_elet_eta[ibin]; //double err_eta = ewz_elet_eta[ibin]; eff = (eff*eff_eta)/avgrate; } if (flavor == 13) {// weight prompt muon //if eta>0.1 float p0 = -18.21; float p1 = 14.83; float p2 = 0.9312; //float ep0= 5.06; float ep1= 1.9; float ep2=0.00069; if ( fabs(eta) < 0.1) { p0 = 7.459; p1 = 2.615; p2 = 0.5138; //ep0 = 10.4; ep1 = 4.934; ep2 = 0.0034; } double arg = ( pt-p0 )/( 2.*p1 ) ; eff = 0.5 * p2 * (1.+erf(arg)); //err = 0.1*eff; } if (flavor == 14) {// weight muon from tau if (fabs(eta) < 0.1) { float p0 = -1.756; float p1 = 12.38; float p2 = 0.4441; //float ep0= 10.39; float ep1= 7.9; float ep2=0.022; double arg = ( pt-p0 )/( 2.*p1 ) ; eff = 0.5 * p2 * (1.+erf(arg)); //err = 0.1*eff; } else { float p0 = 2.102; float p1 = 0.8293; //float ep0= 0.271; float ep1= 0.0083; eff = p1 - p0/pt; //double err0 = ep0/pt; // d(eff)/dp0 //double err1 = ep1; // d(eff)/dp1 //err = sqrt(err0*err0 + err1*err1 - 2*rho*err0*err1); } } if (flavor == 15) {// weight hadronic tau 1p float wz_tau1p[] = {0.0249278,0.146978,0.225049,0.229212,0.21519,0.206152,0.201559,0.197917,0.209249,0.228336,0.193548,}; //float ewz_tau1p[] ={0.00178577,0.00425252,0.00535052,0.00592126,0.00484684,0.00612941,0.00792099,0.0083006,0.0138307,0.015568,0.0501751,}; int ibin = 0; if (pt > 15) ibin = 1; if (pt > 20) ibin = 2; if (pt > 25) ibin = 3; if (pt > 30) ibin = 4; if (pt > 40) ibin = 5; if (pt > 50) ibin = 6; if (pt > 60) ibin = 7; if (pt > 80) ibin = 8; if (pt > 100) ibin = 9; if (pt > 200) ibin = 10; eff = wz_tau1p[ibin]; //err = ewz_tau1p[ibin]; double avgrate = 0.1718; float wz_tau1p_eta[] = {0.162132,0.176393,0.139619,0.178813,0.185144,0.210027,0.203937,0.178688,0.137034,0.164216,0.163713,}; //float ewz_tau1p_eta[] ={0.00706705,0.00617989,0.00506798,0.00525172,0.00581865,0.00865675,0.00599245,0.00529877,0.00506368,0.00617025,0.00726219,}; ibin = 3; if (eta >= -2.5 && eta < -2.0) ibin = 0; if (eta >= -2.0 && eta < -1.5) ibin = 1; if (eta >= -1.5 && eta < -1.0) ibin = 2; if (eta >= -1.0 && eta < -0.5) ibin = 3; if (eta >= -0.5 && eta < -0.1) ibin = 4; if (eta >= -0.1 && eta < 0.1) ibin = 5; if (eta >= 0.1 && eta < 0.5) ibin = 6; if (eta >= 0.5 && eta < 1.0) ibin = 7; if (eta >= 1.0 && eta < 1.5) ibin = 8; if (eta >= 1.5 && eta < 2.0) ibin = 9; if (eta >= 2.0 && eta < 2.5) ibin = 10; double eff_eta = wz_tau1p_eta[ibin]; //double err_eta = ewz_tau1p_eta[ibin]; eff = (eff*eff_eta)/avgrate; } if (flavor == 16) { //weight hadronic tau 3p float wz_tau3p[] = {0.000587199,0.00247181,0.0013031,0.00280112,}; //float ewz_tau3p[] ={0.000415091,0.000617187,0.000582385,0.00197792,}; int ibin = 0; if (pt > 15) ibin = 1; if (pt > 20) ibin = 2; if (pt > 40) ibin = 3; if (pt > 80) ibin = 4; eff = wz_tau3p[ibin]; //err = ewz_tau3p[ibin]; } return eff; }
Apply the supplied projection on event. Definition at line 70 of file ProjectionApplier.hh. References ProjectionApplier::_applyProjection(). Referenced by DISFinalState::project(). { return pcast<PROJ>(_applyProjection(evt, proj)); }
Apply the supplied projection on event. Definition at line 77 of file ProjectionApplier.hh. References ProjectionApplier::_applyProjection(). { return pcast<PROJ>(_applyProjection(evt, proj)); }
Apply the named projection on event. Definition at line 84 of file ProjectionApplier.hh. References ProjectionApplier::_applyProjection(). { return pcast<PROJ>(_applyProjection(evt, name)); }
Helper for histogram asymmetry calculation.
Definition at line 639 of file Analysis.cc. { const string path = s->path(); *s = YODA::asymm(*h1, *h2); s->setPath(path); }
Helper for histogram asymmetry calculation.
Names & emails of paper/analysis authors. Names and email of authors in 'NAME <EMAIL>' format. The first name in the list should be the primary contact person. Definition at line 132 of file Analysis.hh. References AnalysisInfo::authors(), and Analysis::info(). Incoming beam IDs for this run. Definition at line 35 of file Analysis.cc. References AnalysisHandler::beamIds(), and Analysis::handler(). Referenced by UA5_1982_S875503::finalize(), and UA5_1982_S875503::init().
Incoming beams for this run. Definition at line 31 of file Analysis.cc. References AnalysisHandler::beams(), and Analysis::handler(). Referenced by BABAR_2005_S6181155::analyze(), BABAR_2007_S6895344::analyze(), CLEO_2004_S5809304::analyze(), BABAR_2013_I1238276::analyze(), BELLE_2013_I1216515::analyze(), OPAL_1994_S2927284::analyze(), SLD_2004_S5693039::analyze(), OPAL_1998_S3780481::analyze(), SLD_1999_S3743934::analyze(), OPAL_1995_S3198391::analyze(), OPAL_1997_S3608263::analyze(), OPAL_2000_S4418603::analyze(), ALEPH_2002_S4823664::analyze(), DELPHI_1999_S3960137::analyze(), OPAL_1996_S3257789::analyze(), OPAL_1998_S3702294::analyze(), DELPHI_1995_S3137023::analyze(), SLD_2002_S4869273::analyze(), ALEPH_2001_S4656318::analyze(), DELPHI_2002_069_CONF_603::analyze(), OPAL_1998_S3749908::analyze(), OPAL_1997_S3396100::analyze(), ALEPH_1999_S4193598::analyze(), MC_WPOL::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_2004_S5765862::analyze(), DELPHI_1996_S3430090::analyze(), MC_WPOL::init(), and Analysis::isCompatible().
BibTeX citation key for this article. Definition at line 185 of file Analysis.hh. References AnalysisInfo::bibKey(), and Analysis::info().
BibTeX citation entry for this article. Definition at line 190 of file Analysis.hh. References AnalysisInfo::bibTeX(), and Analysis::info().
Book a 1D histogram with nbins uniformly distributed across the range lower - upper . Definition at line 200 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). Referenced by ATLAS_2014_I1304688::_histogram(), Analysis::bookHisto1D(), ATLAS_2012_I1118269::init(), CMS_2010_S8547297::init(), CMS_2010_S8656010::init(), LHCB_2010_I867355::init(), ALICE_2012_I1181770::init(), CMS_2011_S8950903::init(), ATLAS_2010_CONF_2010_049::init(), CMS_2011_S8941262::init(), CMS_2011_S8968497::init(), CMS_2011_S8973270::init(), CMS_2012_I1090423::init(), CMS_2012_I1184941::init(), CMS_2012_I1193338::init(), CMSTOTEM_2014_I1294140::init(), ATLAS_2014_I1298811::init(), ATLAS_2011_I894867::init(), ALICE_2011_S8909580::init(), TOTEM_2012_I1115294::init(), CMS_2011_S9086218::init(), CMS_2012_I1087342::init(), CMS_2011_S9215166::init(), ATLAS_2014_I1282441::init(), ALICE_2011_S8945144::init(), ATLAS_2012_I1091481::init(), CDF_2007_S7057202::init(), ATLAS_2010_S8591806::init(), CMS_2011_I954992::init(), CMS_2011_S8957746::init(), CMS_2011_S8978280::init(), CMS_2013_I1265659::init(), D0_2011_I895662::init(), LHCF_2012_I1115479::init(), ATLAS_2011_I930220::init(), TOTEM_2012_002::init(), CMS_2013_I1273574::init(), D0_2010_S8570965::init(), MC_JETTAGS::init(), ATLAS_2011_I925932::init(), STAR_2006_S6500200::init(), ATLAS_2012_I1124167::init(), UA5_1987_S1640666::init(), CMS_2013_I1256943::init(), CMS_QCD_10_024::init(), CDF_1997_S3541940::init(), SFM_1984_S1178091::init(), ARGUS_1993_S2789213::init(), CMS_2013_I1261026::init(), CMS_2013_I1258128::init(), CMS_2014_I1298810::init(), D0_2000_S4480767::init(), ATLAS_2011_S9128077::init(), BABAR_2007_S7266081::init(), BELLE_2008_I786560::init(), MC_HINC::init(), CDF_1993_S2742446::init(), CDF_2000_S4155203::init(), CDF_2006_S6450792::init(), CDF_2005_S6080774::init(), LHCB_2012_I1208102::init(), MC_JetSplittings::init(), MC_PDFS::init(), MC_WWINC::init(), MC_ZJETS::init(), MC_ZZINC::init(), CDF_2008_S8093652::init(), STAR_2006_S6870392::init(), UA5_1982_S875503::init(), UA5_1989_S1926373::init(), ARGUS_1993_S2669951::init(), ATLAS_2012_I1204447::init(), CMS_2013_I1209721::init(), ATLAS_2013_I1190187::init(), D0_2008_S6879055::init(), ATLAS_2011_I944826::init(), ALEPH_1991_S2435284::init(), ATLAS_2012_I1183818::init(), CMS_2012_I1102908::init(), ATLAS_2011_S9035664::init(), ATLAS_2012_I1204784::init(), CMS_2013_I1272853::init(), ATLAS_2013_I1219109::init(), D0_2000_I499943::init(), ATLAS_2014_I1268975::init(), ATLAS_2014_I1325553::init(), E735_1998_S3905616::init(), MC_DIJET::init(), MC_DIPHOTON::init(), MC_WINC::init(), MC_WJETS::init(), MC_ZINC::init(), MC_ZZJETS::init(), MC_ParticleAnalysis::init(), PDG_TAUS::init(), UA5_1986_S1583476::init(), ALEPH_1996_S3196992::init(), ATLAS_2010_S8817804::init(), D0_1996_S3324664::init(), D0_2007_S7075677::init(), D0_2009_S8202443::init(), D0_2010_S8821313::init(), DELPHI_1999_S3960137::init(), EXAMPLE_CUTS::init(), ALEPH_2002_S4823664::init(), MC_HJETS::init(), CDF_1988_S1865951::init(), MC_JetAnalysis::init(), CDF_1990_S2089246::init(), MC_PHOTONINC::init(), OPAL_1995_S3198391::init(), OPAL_1996_S3257789::init(), OPAL_1997_S3608263::init(), OPAL_1998_S3702294::init(), OPAL_1998_S3749908::init(), OPAL_2000_S4418603::init(), ATLAS_2011_I945498::init(), ATLAS_2011_I954993::init(), CMS_2012_I1107658::init(), D0_2008_S7554427::init(), D0_2008_S7863608::init(), D0_2010_S8671338::init(), H1_1995_S3167097::init(), MC_IDENTIFIED::init(), CDF_1994_S2952106::init(), MC_PHOTONJETS::init(), MC_WWJETS::init(), CDF_2008_S7828950::init(), ATLAS_2012_I1119557::init(), ZEUS_2001_S4815815::init(), ATLAS_2012_I1180197::init(), ATLAS_2012_I1188891::init(), MC_PHOTONJETUE::init(), D0_2001_S4674421::init(), ATLAS_2011_S9131140::init(), D0_2008_S7837160::init(), ATLAS_2011_S9212353::init(), ATLAS_2014_I1306615::init(), CDF_2008_S7540469::init(), ATLAS_2012_CONF_2012_105::init(), ATLAS_2012_CONF_2012_104::init(), LHCB_2011_I919315::init(), ATLAS_2014_I1300647::init(), UA1_1990_S2044935::init(), MC_GENERIC::init(), MC_TTBAR::init(), OPAL_1997_S3396100::init(), ATLAS_2012_I1095236::init(), ALICE_2010_S8624100::init(), ATLAS_2012_I1126136::init(), ATLAS_2012_CONF_2012_103::init(), JADE_1998_S3612880::init(), CDF_1996_S3108457::init(), ATLAS_2011_S8983313::init(), ATLAS_2011_CONF_2011_098::init(), ATLAS_2012_I1186556::init(), JADE_OPAL_2000_S4300807::init(), STAR_2006_S6860818::init(), D0_1996_S3214044::init(), CMS_2012_I941555::init(), CDF_2008_S8095620::init(), EXAMPLE::init(), ATLAS_2013_I1217863_W::init(), ALEPH_2004_S5765862::init(), ATLAS_2011_I926145::init(), CMS_2014_I1303894::init(), ATLAS_2012_CONF_2012_109::init(), MC_XS::init(), ATLAS_2013_I1217863_Z::init(), DELPHI_1995_S3137023::init(), ATLAS_2012_I1190891::init(), ATLAS_2012_I1112263::init(), ATLAS_2012_I1125961::init(), ALICE_2010_S8625980::init(), ATLAS_2012_CONF_2012_001::init(), CMS_2013_I1218372::init(), ALICE_2010_S8706239::init(), ATLAS_2014_I1288706::init(), ATLAS_2011_S9212183::init(), LHCB_2013_I1218996::init(), MC_HFJETS::init(), CDF_1996_S3349578::init(), CDF_2009_NOTE_9936::init(), ATLAS_2013_I1263495::init(), CDF_2009_S8436959::init(), CDF_2000_S4266730::init(), MC_SUSY::init(), D0_2009_S8320160::init(), ATLAS_2011_S8971293::init(), CDF_2009_S8383952::init(), D0_2009_S8349509::init(), CDF_2001_S4563131::init(), CDF_1998_S3618439::init(), CDF_2001_S4517016::init(), ATLAS_2012_I1082936::init(), ALEPH_2001_S4656318::init(), DELPHI_2002_069_CONF_603::init(), D0_2006_S6438750::init(), D0_2010_S8566488::init(), ATLAS_2010_S8919674::init(), SLD_2002_S4869273::init(), SLD_1996_S3398250::init(), MC_PHOTONS::init(), CDF_1996_S3418421::init(), ATLAS_2011_I921594::init(), LHCB_2012_I1119400::init(), LHCB_2014_I1281685::init(), CDF_2008_S7541902::init(), CDF_2006_S6653332::init(), ATLAS_2011_S9108483::init(), D0_2008_S7662670::init(), MC_WPOL::init(), ATLAS_2014_I1306294::init(), ATLAS_2013_I1230812::init(), ATLAS_2014_I1307243::init(), ATLAS_2012_I943401::init(), TASSO_1990_S2148048::init(), ATLAS_2012_I946427::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2011_S9225137::init(), ATLAS_2012_I1117704::init(), ATLAS_2012_I1199269::init(), ATLAS_2011_S9019561::init(), ATLAS_2012_I1082009::init(), ATLAS_2012_I1084540::init(), CMS_2013_I1224539_DIJET::init(), ATLAS_2012_CONF_2012_153::init(), ATLAS_2012_I1083318::init(), CMS_2013_I1224539_WJET::init(), CDF_2009_S8233977::init(), ATLAS_2013_I1217867::init(), CMS_2013_I1224539_ZJET::init(), D0_2004_S5992206::init(), ATLAS_2010_S8914702::init(), ATLAS_2014_I1312627::init(), CDF_2001_S4751469::init(), ATLAS_2010_S8918562::init(), ATLAS_2011_S9120807::init(), D0_2008_S7719523::init(), ATLAS_2011_S9041966::init(), DELPHI_1996_S3430090::init(), ALEPH_1996_S3486095::init(), ALEPH_1999_S4193598::init(), ATLAS_2012_I1093738::init(), OPAL_2004_S6132243::init(), OPAL_1994_S2927284::init(), BELLE_2013_I1216515::init(), STAR_2008_S7869363::init(), BABAR_2007_S6895344::init(), MC_VH2BB::init(), BABAR_2005_S6181155::init(), BELLE_2001_S4598261::init(), ATLAS_2011_I919017::init(), OPAL_2001_S4553896::init(), BABAR_2013_I1238276::init(), DELPHI_2003_WUD_03_11::init(), ATLAS_2012_I1203852::init(), CDF_2004_S5839831::init(), BABAR_2003_I593379::init(), ARGUS_1993_S2653028::init(), CLEO_2004_S5809304::init(), OPAL_1998_S3780481::init(), OPAL_1993_S2692198::init(), ATLAS_2012_I1094564::init(), H1_1994_S2919893::init(), H1_2000_S4129130::init(), SLD_2004_S5693039::init(), SLD_1999_S3743934::init(), PDG_HADRON_MULTIPLICITIES::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), ATLAS_2014_I1279489::initialisePlots(), ATLAS_2012_I1125575::initializeHistograms(), ATLAS_2012_I1094568::initializePlots(), and ATLAS_2011_S9126244::initializePlots(). { const string path = histoPath(hname); Histo1DPtr hist( new Histo1D(nbins, lower, upper, path, title) ); addAnalysisObject(hist); MSG_TRACE("Made histogram " << hname << " for " << name()); hist->setAnnotation("XLabel", xtitle); hist->setAnnotation("YLabel", ytitle); return hist; }
Book a 1D histogram with non-uniform bins defined by the vector of bin edges binedges . Definition at line 215 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). { const string path = histoPath(hname); Histo1DPtr hist( new Histo1D(binedges, path, title) ); addAnalysisObject(hist); MSG_TRACE("Made histogram " << hname << " for " << name()); hist->setAnnotation("XLabel", xtitle); hist->setAnnotation("YLabel", ytitle); return hist; }
Book a 1D histogram with binning from a reference scatter. Definition at line 230 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). { const string path = histoPath(hname); Histo1DPtr hist( new Histo1D(refscatter, path) ); addAnalysisObject(hist); MSG_TRACE("Made histogram " << hname << " for " << name()); hist->setTitle(title); hist->setAnnotation("XLabel", xtitle); hist->setAnnotation("YLabel", ytitle); return hist; }
Book a 1D histogram, using the binnings in the reference data histogram. Definition at line 246 of file Analysis.cc. References Analysis::bookHisto1D(), and Analysis::refData(). { const Scatter2D& refdata = refData(hname); return bookHisto1D(hname, refdata, title, xtitle, ytitle); }
Book a 1D histogram, using the binnings in the reference data histogram. The paper, dataset and x/y-axis IDs will be used to build the histo name in the HepData standard way. Definition at line 255 of file Analysis.cc. References Analysis::bookHisto1D(), and Analysis::makeAxisCode(). { const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId); return bookHisto1D(axisCode, title, xtitle, ytitle); }
Book a 2D histogram with nxbins and nybins uniformly distributed across the ranges xlower - xupper and ylower - yupper respectively along the x- and y-axis.
Definition at line 270 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). Referenced by MC_PDFS::init(). { const string path = histoPath(hname); Histo2DPtr hist( new Histo2D(nxbins, xlower, xupper, nybins, ylower, yupper, path, title) ); addAnalysisObject(hist); MSG_TRACE("Made 2D histogram " << hname << " for " << name()); hist->setAnnotation("XLabel", xtitle); hist->setAnnotation("YLabel", ytitle); hist->setAnnotation("ZLabel", ztitle); return hist; }
Book a 2D histogram with non-uniform bins defined by the vectorx of bin edges xbinedges and ybinedges. Definition at line 289 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). { const string path = histoPath(hname); Histo2DPtr hist( new Histo2D(xbinedges, ybinedges, path, title) ); addAnalysisObject(hist); MSG_TRACE("Made 2D histogram " << hname << " for " << name()); hist->setAnnotation("XLabel", xtitle); hist->setAnnotation("YLabel", ytitle); hist->setAnnotation("ZLabel", ztitle); return hist; }
Book a 1D profile histogram with nbins uniformly distributed across the range lower - upper . Definition at line 349 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). Referenced by Analysis::bookProfile1D(), ATLAS_2010_S8894728::init(), ATLAS_2011_S8994773::init(), CMS_2012_PAS_QCD_11_010::init(), ATLAS_2014_I1298811::init(), ALICE_2011_S8945144::init(), ATLAS_2010_S8591806::init(), CDF_2012_NOTE10874::init(), STAR_2008_S7993412::init(), ATLAS_2012_I1124167::init(), CMS_2013_I1261026::init(), CDF_2008_S7782535::init(), ATLAS_2013_I1243871::init(), D0_1996_S3324664::init(), ATLAS_2011_S8924791::init(), CDF_2005_S6217184::init(), CMS_2012_I1107658::init(), MC_GENERIC::init(), MC_LEADJETUE::init(), UA1_1990_S2044935::init(), MC_PHOTONJETUE::init(), STAR_2009_UE_HELEN::init(), STAR_2006_S6860818::init(), ALICE_2010_S8706239::init(), ALEPH_2001_S4656318::init(), DELPHI_2002_069_CONF_603::init(), MC_PHOTONS::init(), MC_WPOL::init(), ATLAS_2014_I1307243::init(), CDF_2009_S8233977::init(), CDF_2010_S8591881_DY::init(), CDF_2010_S8591881_QCD::init(), CDF_2001_S4751469::init(), ATLAS_2010_S8918562::init(), DELPHI_1996_S3430090::init(), CDF_2004_S5839831::init(), H1_1994_S2919893::init(), H1_2000_S4129130::init(), ATLAS_2014_I1279489::initialisePlots(), ATLAS_2011_S9126244::initializePlots(), and ATLAS_2012_I1125575::initializeProfiles(). { const string path = histoPath(hname); Profile1DPtr prof( new Profile1D(nbins, lower, upper, path, title) ); addAnalysisObject(prof); MSG_TRACE("Made profile histogram " << hname << " for " << name()); prof->setAnnotation("XLabel", xtitle); prof->setAnnotation("YLabel", ytitle); return prof; }
Book a 1D profile histogram with non-uniform bins defined by the vector of bin edges binedges . Definition at line 364 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). { const string path = histoPath(hname); Profile1DPtr prof( new Profile1D(binedges, path, title) ); addAnalysisObject(prof); MSG_TRACE("Made profile histogram " << hname << " for " << name()); prof->setAnnotation("XLabel", xtitle); prof->setAnnotation("YLabel", ytitle); return prof; }
Book a 1D profile histogram with binning from a reference scatter. Definition at line 379 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). { const string path = histoPath(hname); Profile1DPtr prof( new Profile1D(refscatter, path) ); addAnalysisObject(prof); MSG_TRACE("Made profile histogram " << hname << " for " << name()); prof->setTitle(title); prof->setAnnotation("XLabel", xtitle); prof->setAnnotation("YLabel", ytitle); return prof; }
Book a 1D profile histogram, using the binnings in the reference data histogram. Definition at line 395 of file Analysis.cc. References Analysis::bookProfile1D(), and Analysis::refData(). { const Scatter2D& refdata = refData(hname); return bookProfile1D(hname, refdata, title, xtitle, ytitle); }
Book a 1D profile histogram, using the binnings in the reference data histogram. The paper, dataset and x/y-axis IDs will be used to build the histo name in the HepData standard way. Definition at line 404 of file Analysis.cc. References Analysis::bookProfile1D(), and Analysis::makeAxisCode(). { const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId); return bookProfile1D(axisCode, title, xtitle, ytitle); }
Book a 2D profile histogram with nxbins and nybins uniformly distributed across the ranges xlower - xupper and ylower - yupper respectively along the x- and y-axis. Definition at line 417 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). { const string path = histoPath(hname); Profile2DPtr prof( new Profile2D(nxbins, xlower, xupper, nybins, ylower, yupper, path, title) ); addAnalysisObject(prof); MSG_TRACE("Made 2D profile histogram " << hname << " for " << name()); prof->setAnnotation("XLabel", xtitle); prof->setAnnotation("YLabel", ytitle); prof->setAnnotation("ZLabel", ztitle); return prof; }
Book a 2D profile histogram with non-uniform bins defined by the vectorx of bin edges xbinedges and ybinedges. Definition at line 436 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name(). { const string path = histoPath(hname); Profile2DPtr prof( new Profile2D(xbinedges, ybinedges, path, title) ); addAnalysisObject(prof); MSG_TRACE("Made 2D profile histogram " << hname << " for " << name()); prof->setAnnotation("XLabel", xtitle); prof->setAnnotation("YLabel", ytitle); prof->setAnnotation("ZLabel", ztitle); return prof; }
Book a 2-dimensional data point set with the given name.
Definition at line 506 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, Analysis::name(), Analysis::refData(), and Rivet::s. Referenced by Analysis::bookScatter2D(), MC_JetAnalysis::finalize(), MC_ParticleAnalysis::finalize(), DELPHI_2000_S4328825::finalize(), OPAL_2002_S5361494::finalize(), BABAR_2007_S7266081::finalize(), LHCB_2011_I917009::finalize(), ARGUS_1993_S2669951::finalize(), ATLAS_2014_I1307243::finalize(), ALEPH_2004_S5765862::finalize(), SLD_2004_S5693039::finalize(), ATLAS_2011_S9126244::finalize(), SLD_1999_S3743934::finalize(), ATLAS_2014_I1279489::finalizeEfficiencies(), ATLAS_2014_I1312627::hInit(), ATLAS_2010_S8894728::init(), ALICE_2012_I1181770::init(), ALICE_2011_S8909580::init(), LHCB_2013_I1208105::init(), ATLAS_2011_S9002537::init(), CMS_2011_S8978280::init(), CMS_2011_S9088458::init(), STAR_2006_S6500200::init(), CMS_2014_I1298810::init(), ATLAS_2011_S9128077::init(), MC_JetSplittings::init(), CDF_2008_S7782535::init(), MC_WINC::init(), ATLAS_2011_I944826::init(), CMS_2012_I1102908::init(), MC_JetAnalysis::init(), CDF_2005_S6217184::init(), MC_ParticleAnalysis::init(), ATLAS_2011_I945498::init(), CDF_1994_S2952106::init(), ATLAS_2012_I1188891::init(), D0_2001_S4674421::init(), D0_2008_S7837160::init(), MC_GENERIC::init(), STAR_2006_S6860818::init(), MC_XS::init(), CMS_2013_I1218372::init(), LHCB_2012_I1119400::init(), CDF_1996_S3418421::init(), UA5_1988_S1867512::init(), CDF_2008_S7541902::init(), ATLAS_2013_I1230812::init(), ATLAS_2014_I1307243::init(), ATLAS_2012_I1083318::init(), D0_2008_S7719523::init(), ATLAS_2012_I1093734::init(), SLD_2004_S5693039::init(), SLD_1999_S3743934::init(), ATLAS_2012_I1094568::initializePlots(), ATLAS_2011_S9126244::initializePlots(), and SLD_1996_S3398250::multiplicity_subtract(). { Scatter2DPtr s; const string path = histoPath(hname); if (copy_pts) { const Scatter2D& refdata = refData(hname); s.reset( new Scatter2D(refdata, path) ); foreach (Point2D& p, s->points()) p.setY(0, 0); } else { s.reset( new Scatter2D(path) ); } addAnalysisObject(s); MSG_TRACE("Made scatter " << hname << " for " << name()); s->setTitle(title); s->setAnnotation("XLabel", xtitle); s->setAnnotation("YLabel", ytitle); return s; }
Book a 2-dimensional data point set, using the binnings in the reference data histogram. The paper, dataset and x/y-axis IDs will be used to build the histo name in the HepData standard way.
Definition at line 496 of file Analysis.cc. References Analysis::bookScatter2D(), and Analysis::makeAxisCode(). { const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId); return bookScatter2D(axisCode, copy_pts, title, xtitle, ytitle); }
Book a 2-dimensional data point set with equally spaced x-points in a range. The y values and errors will be set to 0. Definition at line 529 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, Analysis::name(), Rivet::Cuts::pt, and Rivet::s. { const string path = histoPath(hname); Scatter2DPtr s( new Scatter2D(path) ); const double binwidth = (upper-lower)/npts; for (size_t pt = 0; pt < npts; ++pt) { const double bincentre = lower + (pt + 0.5) * binwidth; s->addPoint(bincentre, 0, binwidth/2.0, 0); } addAnalysisObject(s); MSG_TRACE("Made scatter " << hname << " for " << name()); s->setTitle(title); s->setAnnotation("XLabel", xtitle); s->setAnnotation("YLabel", ytitle); return s; }
Book a 2-dimensional data point set based on provided contiguous "bin edges". The y values and errors will be set to 0. Definition at line 550 of file Analysis.cc. References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, Analysis::name(), Rivet::Cuts::pt, and Rivet::s. { const string path = histoPath(hname); Scatter2DPtr s( new Scatter2D(path) ); for (size_t pt = 0; pt < binedges.size()-1; ++pt) { const double bincentre = (binedges[pt] + binedges[pt+1]) / 2.0; const double binwidth = binedges[pt+1] - binedges[pt]; s->addPoint(bincentre, 0, binwidth/2.0, 0); } addAnalysisObject(s); MSG_TRACE("Made scatter " << hname << " for " << name()); s->setTitle(title); s->setAnnotation("XLabel", xtitle); s->setAnnotation("YLabel", ytitle); return s; }
Collider on which the experiment ran. Definition at line 170 of file Analysis.hh. References AnalysisInfo::collider(), and Analysis::info().
Get the process cross-section in pb. Throws if this hasn't been set. Definition at line 155 of file Analysis.cc. References Analysis::_crossSection, Analysis::_gotCrossSection, and Analysis::name(). Referenced by MC_JetSplittings::finalize(), MC_JetAnalysis::finalize(), MC_ParticleAnalysis::finalize(), CDF_2006_S6450792::finalize(), D0_2000_S4480767::finalize(), TOTEM_2012_002::finalize(), CMS_2011_S9086218::finalize(), CDF_2001_S4563131::finalize(), CDF_2008_S7828950::finalize(), D0_2010_S8671338::finalize(), MC_XS::finalize(), CDF_2008_S8093652::finalize(), BELLE_2013_I1216515::finalize(), CDF_2000_S4155203::finalize(), CDF_2000_S4266730::finalize(), CMS_2011_S8941262::finalize(), CMS_2012_I1184941::finalize(), CDF_1998_S3618439::finalize(), CMS_2012_I1087342::finalize(), LHCB_2012_I1208102::finalize(), MC_HINC::finalize(), MC_ZJETS::finalize(), CMS_2011_I954992::finalize(), MC_HJETS::finalize(), MC_WJETS::finalize(), ATLAS_2012_I1118269::finalize(), ATLAS_2011_I894867::finalize(), BABAR_2005_S6181155::finalize(), CDF_2009_S8383952::finalize(), CDF_2001_S4517016::finalize(), D0_2010_S8566488::finalize(), STAR_2006_S6870392::finalize(), MC_ZINC::finalize(), LHCB_2010_I867355::finalize(), CDF_2009_S8436959::finalize(), CMS_2012_I1193338::finalize(), D0_2011_I895662::finalize(), CDF_2005_S6080774::finalize(), ATLAS_2014_I1288706::finalize(), ATLAS_2014_I1325553::finalize(), ATLAS_2014_I1282441::finalize(), D0_2006_S6438750::finalize(), CMS_2013_I1273574::finalize(), MC_DIPHOTON::finalize(), D0_2000_I499943::finalize(), ALICE_2012_I1181770::finalize(), CLEO_2004_S5809304::finalize(), D0_2008_S7662670::finalize(), D0_2008_S7863608::finalize(), MC_ZZJETS::finalize(), MC_WINC::finalize(), ATLAS_2011_S9035664::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ATLAS_2011_I954993::finalize(), MC_WWJETS::finalize(), CDF_2009_S8233977::finalize(), CMS_2014_I1298810::finalize(), LHCB_2013_I1218996::finalize(), ATLAS_2013_I1219109::finalize(), D0_2010_S8570965::finalize(), CMS_2011_S8973270::finalize(), ATLAS_2012_I946427::finalize(), CDF_2008_S7541902::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2011_I921594::finalize(), ATLAS_2013_I1263495::finalize(), MC_ZZINC::finalize(), MC_WWINC::finalize(), ATLAS_2014_I1312627::finalize(), CDF_2006_S6653332::finalize(), ATLAS_2013_I1217863_Z::finalize(), D0_2008_S7719523::finalize(), ATLAS_2013_I1217863_W::finalize(), ATLAS_2014_I1307756::finalize(), ATLAS_2011_I926145::finalize(), CMS_2013_I1256943::finalize(), ATLAS_2014_I1306294::finalize(), ATLAS_2011_S9108483::finalize(), CDF_2008_S7540469::finalize(), ATLAS_2012_I1199269::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2014_I1307243::finalize(), ATLAS_2012_I1083318::finalize(), ATLAS_2014_I1304688::finalize(), ATLAS_2012_CONF_2012_105::finalize(), CMS_2014_I1303894::finalize(), ATLAS_2012_CONF_2012_104::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2012_I1186556::finalize(), ATLAS_2012_I1084540::finalize(), MC_VH2BB::finalize(), ATLAS_2012_CONF_2012_103::finalize(), ATLAS_2012_I1117704::finalize(), ATLAS_2013_I1190187::finalize(), ATLAS_2012_I1093738::finalize(), ATLAS_2012_I1125961::finalize(), ATLAS_2011_CONF_2011_098::finalize(), ATLAS_2011_S8983313::finalize(), ATLAS_2012_I1126136::finalize(), ATLAS_2012_I1095236::finalize(), ATLAS_2012_I1190891::finalize(), ATLAS_2011_S9019561::finalize(), ATLAS_2012_I1112263::finalize(), ATLAS_2011_S9212183::finalize(), ATLAS_2012_CONF_2012_109::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2011_S9225137::finalize(), ATLAS_2012_CONF_2012_001::finalize(), ATLAS_2012_I1203852::finalize(), ATLAS_2012_I943401::finalize(), ATLAS_2012_I1180197::finalize(), ATLAS_2012_CONF_2012_153::finalize(), ATLAS_2012_I1204447::finalize(), ATLAS_2011_S9212353::finalize(), ATLAS_2011_S9041966::finalize(), and ATLAS_2011_I919017::finalize(). { if (!_gotCrossSection || std::isnan(_crossSection)) { string errMsg = "You did not set the cross section for the analysis " + name(); throw Error(errMsg); } return _crossSection; }
Get the process cross-section per generated event in pb. Throws if this hasn't been set. Definition at line 163 of file Analysis.cc. References Analysis::_crossSection, and Analysis::sumOfWeights(). Referenced by CDF_2007_S7057202::finalize(), CDF_1988_S1865951::finalize(), LHCB_2011_I919315::finalize(), MC_PHOTONINC::finalize(), ATLAS_2014_I1268975::finalize(), ATLAS_2010_S8817804::finalize(), MC_PHOTONJETS::finalize(), ATLAS_2012_I1082936::finalize(), ATLAS_2011_I930220::finalize(), D0_2009_S8349509::finalize(), UA1_1990_S2044935::finalize(), D0_2001_S4674421::finalize(), MC_WPOL::finalize(), ATLAS_2011_S9128077::finalize(), ATLAS_2013_I1230812::finalize(), ATLAS_2011_I945498::finalize(), and ATLAS_2014_I1306615::finalize(). { const double sumW = sumOfWeights(); assert(sumW != 0.0); return _crossSection / sumW; }
Get a full description of the analysis. Full textual description of this analysis, what it is useful for, what experimental techniques are applied, etc. Should be treated as a chunk of restructuredText (http://docutils.sourceforge.net/rst.html), with equations to be rendered as LaTeX with amsmath operators. Definition at line 151 of file Analysis.hh. References AnalysisInfo::description(), and Analysis::info(). { return info().description(); }
Helper for histogram division.
Definition at line 574 of file Analysis.cc. Referenced by D0_2008_S7837160::calc_asymm(), MC_JetAnalysis::finalize(), MC_ParticleAnalysis::finalize(), CMS_2011_S9088458::finalize(), ATLAS_2012_I1094061::HistoPair::finalize(), CDF_1996_S3418421::finalize(), ALICE_2011_S8909580::finalize(), STAR_2006_S6500200::finalize(), CMS_2011_S8978280::finalize(), ALICE_2012_I1181770::finalize(), MC_WINC::finalize(), CMS_2014_I1298810::finalize(), LHCB_2011_I917009::finalize(), MC_GENERIC::finalize(), ATLAS_2014_I1312627::finalize(), ATLAS_2012_I1188891::finalize(), D0_2008_S7719523::finalize(), STAR_2006_S6860818::finalize(), LHCB_2012_I1119400::finalize(), CMS_2013_I1218372::finalize(), ATLAS_2011_S9128077::finalize(), ATLAS_2014_I1307243::finalize(), ATLAS_2011_I944826::finalize(), SLD_2004_S5693039::finalize(), ATLAS_2011_S9126244::finalize(), SLD_1999_S3743934::finalize(), and ATLAS_2014_I1279489::finalizeEfficiencies().
Helper for histogram division with raw YODA objects.
Helper for profile histogram division.
Helper for profile histogram division with raw YODA objects.
Helper for 2D histogram division.
Helper for 2D histogram division with raw YODA objects.
Helper for 2D profile histogram division.
Helper for 2D profile histogram division with raw YODA objects
Helper for histogram efficiency calculation.
Definition at line 626 of file Analysis.cc. Referenced by CMS_2012_I1102908::finalize(), and ATLAS_2014_I1307243::finalize(). { const string path = s->path(); *s = YODA::efficiency(*h1, *h2); s->setPath(path); }
Helper for histogram efficiency calculation.
Experiment which performed and published this analysis. Definition at line 165 of file Analysis.hh. References AnalysisInfo::experiment(), and Analysis::info(). { return info().experiment(); }
function fills map EventCountsPerSR by looping over all signal regions and looking if the event falls into this signal region Definition at line 579 of file ATLAS_2012_I1204447.cc. References ATLAS_2012_I1204447::_eventCountsPerSR, ATLAS_2012_I1204447::getCutsPerSignalRegion(), and Rivet::toString(). Referenced by ATLAS_2012_I1204447::analyze(). { // Get cut values for HTlep, loop over them and add event if cut is passed vector<int> cut_values = getCutsPerSignalRegion("HTlep", onZ); for (size_t i = 0; i < cut_values.size(); i++) { if (HTlep > cut_values[i]) _eventCountsPerSR[("HTlep_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight; } // Get cut values for METStrong, loop over them and add event if cut is passed cut_values = getCutsPerSignalRegion("METStrong", onZ); for (size_t i = 0; i < cut_values.size(); i++) { if (eTmiss > cut_values[i] && HTjets > 100.) _eventCountsPerSR[("METStrong_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight; } // Get cut values for METWeak, loop over them and add event if cut is passed cut_values = getCutsPerSignalRegion("METWeak", onZ); for (size_t i = 0; i < cut_values.size(); i++) { if (eTmiss > cut_values[i] && HTjets <= 100.) _eventCountsPerSR[("METWeak_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight; } // Get cut values for Meff, loop over them and add event if cut is passed cut_values = getCutsPerSignalRegion("Meff", onZ); for (size_t i = 0; i < cut_values.size(); i++) { if (meff > cut_values[i]) _eventCountsPerSR[("Meff_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight; } // Get cut values for MeffStrong, loop over them and add event if cut is passed cut_values = getCutsPerSignalRegion("MeffStrong", onZ); for (size_t i = 0; i < cut_values.size(); i++) { if (meff > cut_values[i] && eTmiss > 75.) _eventCountsPerSR[("MeffStrong_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight; } }
Normalise histograms etc., after the run. Reimplemented from Analysis. Definition at line 425 of file ATLAS_2012_I1204447.cc. References ATLAS_2012_I1204447::_eventCountsPerSR, ATLAS_2012_I1204447::_h_e_n, ATLAS_2012_I1204447::_h_excluded, ATLAS_2012_I1204447::_h_HTjets_all, ATLAS_2012_I1204447::_h_HTlep_all, ATLAS_2012_I1204447::_h_Meff_all, ATLAS_2012_I1204447::_h_MET_all, ATLAS_2012_I1204447::_h_mu_n, ATLAS_2012_I1204447::_h_pt_1_2ltau, ATLAS_2012_I1204447::_h_pt_1_3l, ATLAS_2012_I1204447::_h_pt_2_2ltau, ATLAS_2012_I1204447::_h_pt_2_3l, ATLAS_2012_I1204447::_h_pt_3_2ltau, ATLAS_2012_I1204447::_h_pt_3_3l, ATLAS_2012_I1204447::_h_tau_n, ATLAS_2012_I1204447::_signal_regions, Analysis::crossSection(), Rivet::femtobarn, ATLAS_2012_I1204447::getUpperLimit(), Analysis::scale(), and Analysis::sumOfWeights(). { // Normalize to an integrated luminosity of 1 fb-1 double norm = crossSection()/femtobarn/sumOfWeights(); string best_signal_region = ""; double ratio_best_SR = 0.; // Loop over all signal regions and find signal region with best sensitivity (ratio signal events/visible cross-section) for (size_t i = 0; i < _signal_regions.size(); i++) { double signal_events = _eventCountsPerSR[_signal_regions[i]] * norm; // Use expected upper limits to find best signal region double UL95 = getUpperLimit(_signal_regions[i], false); double ratio = signal_events / UL95; if (ratio > ratio_best_SR) { best_signal_region = _signal_regions[i]; ratio_best_SR = ratio; } } double signal_events_best_SR = _eventCountsPerSR[best_signal_region] * norm; double exp_UL_best_SR = getUpperLimit(best_signal_region, false); double obs_UL_best_SR = getUpperLimit(best_signal_region, true); // Print out result cout << "----------------------------------------------------------------------------------------" << endl; cout << "Best signal region: " << best_signal_region << endl; cout << "Normalized number of signal events in this best signal region (per fb-1): " << signal_events_best_SR << endl; cout << "Efficiency*Acceptance: " << _eventCountsPerSR[best_signal_region]/sumOfWeights() << endl; cout << "Cross-section [fb]: " << crossSection()/femtobarn << endl; cout << "Expected visible cross-section (per fb-1): " << exp_UL_best_SR << endl; cout << "Ratio (signal events / expected visible cross-section): " << ratio_best_SR << endl; cout << "Observed visible cross-section (per fb-1): " << obs_UL_best_SR << endl; cout << "Ratio (signal events / observed visible cross-section): " << signal_events_best_SR/obs_UL_best_SR << endl; cout << "----------------------------------------------------------------------------------------" << endl; cout << "Using the EXPECTED limits (visible cross-section) of the analysis: " << endl; if (signal_events_best_SR > exp_UL_best_SR) { cout << "Since the number of signal events > the visible cross-section, this model/grid point is EXCLUDED with 95% CL." << endl; _h_excluded->fill(1); } else { cout << "Since the number of signal events < the visible cross-section, this model/grid point is NOT EXCLUDED." << endl; _h_excluded->fill(0); } cout << "----------------------------------------------------------------------------------------" << endl; cout << "Using the OBSERVED limits (visible cross-section) of the analysis: " << endl; if (signal_events_best_SR > obs_UL_best_SR) { cout << "Since the number of signal events > the visible cross-section, this model/grid point is EXCLUDED with 95% CL." << endl; _h_excluded->fill(1); } else { cout << "Since the number of signal events < the visible cross-section, this model/grid point is NOT EXCLUDED." << endl; _h_excluded->fill(0); } cout << "----------------------------------------------------------------------------------------" << endl; // Normalize to cross section if (norm != 0) { scale(_h_HTlep_all, norm); scale(_h_HTjets_all, norm); scale(_h_MET_all, norm); scale(_h_Meff_all, norm); scale(_h_pt_1_3l, norm); scale(_h_pt_2_3l, norm); scale(_h_pt_3_3l, norm); scale(_h_pt_1_2ltau, norm); scale(_h_pt_2_2ltau, norm); scale(_h_pt_3_2ltau, norm); scale(_h_e_n, norm); scale(_h_mu_n, norm); scale(_h_tau_n, norm); scale(_h_excluded, signal_events_best_SR); } }
Function calculating the prong number of taus Definition at line 636 of file ATLAS_2012_I1204447.cc. References Rivet::PID::charge(), Rivet::PID::ELECTRON, Rivet::PID::MUON, and Rivet::PID::TAU. Referenced by ATLAS_2012_I1204447::analyze(). { assert(p != NULL); //const int tau_barcode = p->barcode(); const GenVertex* dv = p->end_vertex(); assert(dv != NULL); for (GenVertex::particles_out_const_iterator pp = dv->particles_out_const_begin(); pp != dv->particles_out_const_end(); ++pp) { // If they have status 1 and are charged they will produce a track and the prong number is +1 if ((*pp)->status() == 1 ) { const int id = (*pp)->pdg_id(); if (Rivet::PID::charge(id) != 0 ) ++nprong; // Check if tau decays leptonically // @todo Can a tau decay include a tau in its decay daughters?! if ((abs(id) == PID::ELECTRON || abs(id) == PID::MUON || abs(id) == PID::TAU) && abs(p->pdg_id()) == PID::TAU) lep_decaying_tau = true; } // If the status of the daughter particle is 2 it is unstable and the further decays are checked else if ((*pp)->status() == 2 ) { get_prong_number(*pp, nprong, lep_decaying_tau); } } }
Function returning 4-vector of daughter-particle if it is a tau neutrino Definition at line 623 of file ATLAS_2012_I1204447.cc. References Particle::abspid(), Particle::genParticle(), Rivet::PID::NU_TAU, and Rivet::PID::TAU. Referenced by ATLAS_2012_I1204447::analyze(). { assert(p.abspid() == PID::TAU); const GenVertex* dv = p.genParticle()->end_vertex(); assert(dv != NULL); for (GenVertex::particles_out_const_iterator pp = dv->particles_out_const_begin(); pp != dv->particles_out_const_end(); ++pp) { if ((*pp)->pdg_id() == PID::NU_TAU) return FourMomentum((*pp)->momentum()); } return FourMomentum(); }
Get a data object from the histogram system
Definition at line 711 of file Analysis.hh. References Analysis::analysisObjects(), and Analysis::histoPath(). { foreach (const AnalysisObjectPtr& ao, analysisObjects()) { if (ao->path() == histoPath(name)) return dynamic_pointer_cast<AO>(ao); } throw Exception("Data object " + histoPath(name) + " not found"); }
Get a data object from the histogram system (non-const)
Definition at line 722 of file Analysis.hh. References Analysis::analysisObjects(), and Analysis::histoPath(). { foreach (const AnalysisObjectPtr& ao, analysisObjects()) { if (ao->path() == histoPath(name)) return dynamic_pointer_cast<AO>(ao); } throw Exception("Data object " + histoPath(name) + " not found"); }
Function giving all cut vales per kinematic variable (taking onZ for MET into account) Definition at line 547 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::fillEventCountsPerSR(), and ATLAS_2012_I1204447::getSignalRegions(). { vector<int> cutValues; // Cut values for HTlep if (signal_region.compare("HTlep") == 0) { cutValues.push_back(0); cutValues.push_back(100); cutValues.push_back(150); cutValues.push_back(200); cutValues.push_back(300); } // Cut values for METStrong (HTjets > 100 GeV) and METWeak (HTjets < 100 GeV) else if (signal_region.compare("METStrong") == 0 || signal_region.compare("METWeak") == 0) { if (onZ == 0) cutValues.push_back(0); else if (onZ == 1) cutValues.push_back(20); cutValues.push_back(50); cutValues.push_back(75); } // Cut values for Meff and MeffStrong (MET > 75 GeV) if (signal_region.compare("Meff") == 0 || signal_region.compare("MeffStrong") == 0) { cutValues.push_back(0); cutValues.push_back(150); cutValues.push_back(300); cutValues.push_back(500); } return cutValues; }
Get a named Histo1D object from the histogram system. Definition at line 737 of file Analysis.hh. References Analysis::name(). { return getAnalysisObject<Histo1D>(name); }
Get a named Histo1D object from the histogram system (non-const) Definition at line 742 of file Analysis.hh. References Analysis::name(). { return getAnalysisObject<Histo1D>(name); }
Get a Histo1D object from the histogram system by axis ID codes (non-const) Definition at line 747 of file Analysis.hh. References Analysis::makeAxisCode(). { return getAnalysisObject<Histo1D>(makeAxisCode(datasetId, xAxisId, yAxisId)); }
Get a Histo1D object from the histogram system by axis ID codes (non-const) Definition at line 752 of file Analysis.hh. References Analysis::makeAxisCode(). { return getAnalysisObject<Histo1D>(makeAxisCode(datasetId, xAxisId, yAxisId)); } Get a Log object based on the name() property of the calling analysis object. Reimplemented from ProjectionApplier. Definition at line 82 of file Analysis.cc. References Analysis::name(). { string logname = "Rivet.Analysis." + name(); return Log::getLog(logname); }
Get a named Profile1D object from the histogram system. Definition at line 779 of file Analysis.hh. References Analysis::name(). { return getAnalysisObject<Profile1D>(name); }
Get a named Profile1D object from the histogram system (non-const) Definition at line 784 of file Analysis.hh. References Analysis::name(). { return getAnalysisObject<Profile1D>(name); }
Get a Profile1D object from the histogram system by axis ID codes (non-const) Definition at line 789 of file Analysis.hh. References Analysis::makeAxisCode(). { return getAnalysisObject<Profile1D>(makeAxisCode(datasetId, xAxisId, yAxisId)); }
Get a Profile1D object from the histogram system by axis ID codes (non-const) Definition at line 794 of file Analysis.hh. References Analysis::makeAxisCode(). { return getAnalysisObject<Profile1D>(makeAxisCode(datasetId, xAxisId, yAxisId)); }
Get the named projection, specifying return type via a template argument. Definition at line 52 of file ProjectionApplier.hh. References ProjectionHandler::getProjection(), and ProjectionApplier::getProjHandler(). Referenced by ProjectionApplier::_applyProjection(), Rivet::pcmp(), and Hemispheres::project(). { const Projection& p = getProjHandler().getProjection(*this, name); return pcast<PROJ>(p); }
Get the named projection (non-templated, so returns as a reference to a Projection base class). Definition at line 60 of file ProjectionApplier.hh. References ProjectionHandler::getProjection(), and ProjectionApplier::getProjHandler(). { return getProjHandler().getProjection(*this, name); }
Get the contained projections, including recursion. Definition at line 45 of file ProjectionApplier.hh. References ProjectionHandler::DEEP, ProjectionHandler::getChildProjections(), and ProjectionApplier::getProjHandler(). Referenced by Projection::beamPairs(). { return getProjHandler().getChildProjections(*this, ProjectionHandler::DEEP); }
Get a reference to the ProjectionHandler for this thread. Definition at line 97 of file ProjectionApplier.hh. References ProjectionApplier::_projhandler. Referenced by ProjectionApplier::_addProjection(), ProjectionApplier::getProjection(), ProjectionApplier::getProjections(), and ProjectionApplier::~ProjectionApplier(). { return _projhandler; }
Get a named Scatter2D object from the histogram system. Definition at line 821 of file Analysis.hh. References Analysis::name(). Referenced by ATLAS_2014_I1279489::finalizeEfficiencies(). { return getAnalysisObject<Scatter2D>(name); }
Get a named Scatter2D object from the histogram system (non-const) Definition at line 826 of file Analysis.hh. References Analysis::name(). { return getAnalysisObject<Scatter2D>(name); }
Get a Scatter2D object from the histogram system by axis ID codes (non-const) Definition at line 831 of file Analysis.hh. References Analysis::makeAxisCode(). { return getAnalysisObject<Scatter2D>(makeAxisCode(datasetId, xAxisId, yAxisId)); }
Get a Scatter2D object from the histogram system by axis ID codes (non-const) Definition at line 836 of file Analysis.hh. References Analysis::makeAxisCode(). { return getAnalysisObject<Scatter2D>(makeAxisCode(datasetId, xAxisId, yAxisId)); }
Helper functions. Function giving a list of all signal regions Definition at line 510 of file ATLAS_2012_I1204447.cc. References ATLAS_2012_I1204447::getCutsPerSignalRegion(), and Rivet::toString(). Referenced by ATLAS_2012_I1204447::init(). { // List of basic signal regions vector<string> basic_signal_regions; basic_signal_regions.push_back("3l_offZ"); basic_signal_regions.push_back("3l_onZ"); basic_signal_regions.push_back("2ltau_offZ"); basic_signal_regions.push_back("2ltau_onZ"); // List of kinematic variables vector<string> kinematic_variables; kinematic_variables.push_back("HTlep"); kinematic_variables.push_back("METStrong"); kinematic_variables.push_back("METWeak"); kinematic_variables.push_back("Meff"); kinematic_variables.push_back("MeffStrong"); vector<string> signal_regions; // Loop over all kinematic variables and basic signal regions for (size_t i0 = 0; i0 < kinematic_variables.size(); i0++) { for (size_t i1 = 0; i1 < basic_signal_regions.size(); i1++) { // Is signal region onZ? int onZ = (basic_signal_regions[i1].find("onZ") != string::npos) ? 1 : 0; // Get cut values for this kinematic variable vector<int> cut_values = getCutsPerSignalRegion(kinematic_variables[i0], onZ); // Loop over all cut values for (size_t i2 = 0; i2 < cut_values.size(); i2++) { // push signal region into vector signal_regions.push_back( (kinematic_variables[i0] + "_" + basic_signal_regions[i1] + "_cut_" + toString(i2)) ); } } } return signal_regions; }
Function giving observed upper limit (visible cross-section) Definition at line 831 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::finalize(). { map<string,double> upperLimitsObserved; upperLimitsObserved["HTlep_3l_offZ_cut_0"] = 11.; upperLimitsObserved["HTlep_3l_offZ_cut_100"] = 8.7; upperLimitsObserved["HTlep_3l_offZ_cut_150"] = 4.0; upperLimitsObserved["HTlep_3l_offZ_cut_200"] = 4.4; upperLimitsObserved["HTlep_3l_offZ_cut_300"] = 1.6; upperLimitsObserved["HTlep_2ltau_offZ_cut_0"] = 25.; upperLimitsObserved["HTlep_2ltau_offZ_cut_100"] = 14.; upperLimitsObserved["HTlep_2ltau_offZ_cut_150"] = 6.1; upperLimitsObserved["HTlep_2ltau_offZ_cut_200"] = 3.3; upperLimitsObserved["HTlep_2ltau_offZ_cut_300"] = 1.2; upperLimitsObserved["HTlep_3l_onZ_cut_0"] = 48.; upperLimitsObserved["HTlep_3l_onZ_cut_100"] = 38.; upperLimitsObserved["HTlep_3l_onZ_cut_150"] = 14.; upperLimitsObserved["HTlep_3l_onZ_cut_200"] = 7.2; upperLimitsObserved["HTlep_3l_onZ_cut_300"] = 4.5; upperLimitsObserved["HTlep_2ltau_onZ_cut_0"] = 85.; upperLimitsObserved["HTlep_2ltau_onZ_cut_100"] = 53.; upperLimitsObserved["HTlep_2ltau_onZ_cut_150"] = 11.0; upperLimitsObserved["HTlep_2ltau_onZ_cut_200"] = 5.2; upperLimitsObserved["HTlep_2ltau_onZ_cut_300"] = 3.0; upperLimitsObserved["METStrong_3l_offZ_cut_0"] = 2.6; upperLimitsObserved["METStrong_3l_offZ_cut_50"] = 2.1; upperLimitsObserved["METStrong_3l_offZ_cut_75"] = 2.1; upperLimitsObserved["METStrong_2ltau_offZ_cut_0"] = 4.2; upperLimitsObserved["METStrong_2ltau_offZ_cut_50"] = 3.1; upperLimitsObserved["METStrong_2ltau_offZ_cut_75"] = 2.6; upperLimitsObserved["METStrong_3l_onZ_cut_20"] = 11.0; upperLimitsObserved["METStrong_3l_onZ_cut_50"] = 6.4; upperLimitsObserved["METStrong_3l_onZ_cut_75"] = 5.1; upperLimitsObserved["METStrong_2ltau_onZ_cut_20"] = 5.9; upperLimitsObserved["METStrong_2ltau_onZ_cut_50"] = 3.4; upperLimitsObserved["METStrong_2ltau_onZ_cut_75"] = 1.2; upperLimitsObserved["METWeak_3l_offZ_cut_0"] = 11.; upperLimitsObserved["METWeak_3l_offZ_cut_50"] = 5.3; upperLimitsObserved["METWeak_3l_offZ_cut_75"] = 3.1; upperLimitsObserved["METWeak_2ltau_offZ_cut_0"] = 23.; upperLimitsObserved["METWeak_2ltau_offZ_cut_50"] = 4.3; upperLimitsObserved["METWeak_2ltau_offZ_cut_75"] = 3.1; upperLimitsObserved["METWeak_3l_onZ_cut_20"] = 41.; upperLimitsObserved["METWeak_3l_onZ_cut_50"] = 16.; upperLimitsObserved["METWeak_3l_onZ_cut_75"] = 8.0; upperLimitsObserved["METWeak_2ltau_onZ_cut_20"] = 80.; upperLimitsObserved["METWeak_2ltau_onZ_cut_50"] = 4.4; upperLimitsObserved["METWeak_2ltau_onZ_cut_75"] = 1.8; upperLimitsObserved["Meff_3l_offZ_cut_0"] = 11.; upperLimitsObserved["Meff_3l_offZ_cut_150"] = 8.1; upperLimitsObserved["Meff_3l_offZ_cut_300"] = 3.1; upperLimitsObserved["Meff_3l_offZ_cut_500"] = 2.1; upperLimitsObserved["Meff_2ltau_offZ_cut_0"] = 25.; upperLimitsObserved["Meff_2ltau_offZ_cut_150"] = 12.; upperLimitsObserved["Meff_2ltau_offZ_cut_300"] = 3.9; upperLimitsObserved["Meff_2ltau_offZ_cut_500"] = 2.2; upperLimitsObserved["Meff_3l_onZ_cut_0"] = 48.; upperLimitsObserved["Meff_3l_onZ_cut_150"] = 37.; upperLimitsObserved["Meff_3l_onZ_cut_300"] = 11.; upperLimitsObserved["Meff_3l_onZ_cut_500"] = 4.8; upperLimitsObserved["Meff_2ltau_onZ_cut_0"] = 85.; upperLimitsObserved["Meff_2ltau_onZ_cut_150"] = 28.; upperLimitsObserved["Meff_2ltau_onZ_cut_300"] = 5.9; upperLimitsObserved["Meff_2ltau_onZ_cut_500"] = 1.9; upperLimitsObserved["MeffStrong_3l_offZ_cut_0"] = 3.8; upperLimitsObserved["MeffStrong_3l_offZ_cut_150"] = 3.8; upperLimitsObserved["MeffStrong_3l_offZ_cut_300"] = 2.8; upperLimitsObserved["MeffStrong_3l_offZ_cut_500"] = 2.1; upperLimitsObserved["MeffStrong_2ltau_offZ_cut_0"] = 3.9; upperLimitsObserved["MeffStrong_2ltau_offZ_cut_150"] = 4.0; upperLimitsObserved["MeffStrong_2ltau_offZ_cut_300"] = 2.9; upperLimitsObserved["MeffStrong_2ltau_offZ_cut_500"] = 1.5; upperLimitsObserved["MeffStrong_3l_onZ_cut_0"] = 10.0; upperLimitsObserved["MeffStrong_3l_onZ_cut_150"] = 10.0; upperLimitsObserved["MeffStrong_3l_onZ_cut_300"] = 6.8; upperLimitsObserved["MeffStrong_3l_onZ_cut_500"] = 3.9; upperLimitsObserved["MeffStrong_2ltau_onZ_cut_0"] = 1.6; upperLimitsObserved["MeffStrong_2ltau_onZ_cut_150"] = 1.4; upperLimitsObserved["MeffStrong_2ltau_onZ_cut_300"] = 1.5; upperLimitsObserved["MeffStrong_2ltau_onZ_cut_500"] = 0.9; // Expected upper limits are also given but not used in this analysis map<string,double> upperLimitsExpected; upperLimitsExpected["HTlep_3l_offZ_cut_0"] = 11.; upperLimitsExpected["HTlep_3l_offZ_cut_100"] = 8.5; upperLimitsExpected["HTlep_3l_offZ_cut_150"] = 4.6; upperLimitsExpected["HTlep_3l_offZ_cut_200"] = 3.6; upperLimitsExpected["HTlep_3l_offZ_cut_300"] = 1.9; upperLimitsExpected["HTlep_2ltau_offZ_cut_0"] = 23.; upperLimitsExpected["HTlep_2ltau_offZ_cut_100"] = 14.; upperLimitsExpected["HTlep_2ltau_offZ_cut_150"] = 6.4; upperLimitsExpected["HTlep_2ltau_offZ_cut_200"] = 3.6; upperLimitsExpected["HTlep_2ltau_offZ_cut_300"] = 1.5; upperLimitsExpected["HTlep_3l_onZ_cut_0"] = 33.; upperLimitsExpected["HTlep_3l_onZ_cut_100"] = 25.; upperLimitsExpected["HTlep_3l_onZ_cut_150"] = 12.; upperLimitsExpected["HTlep_3l_onZ_cut_200"] = 6.5; upperLimitsExpected["HTlep_3l_onZ_cut_300"] = 3.1; upperLimitsExpected["HTlep_2ltau_onZ_cut_0"] = 94.; upperLimitsExpected["HTlep_2ltau_onZ_cut_100"] = 61.; upperLimitsExpected["HTlep_2ltau_onZ_cut_150"] = 9.9; upperLimitsExpected["HTlep_2ltau_onZ_cut_200"] = 4.5; upperLimitsExpected["HTlep_2ltau_onZ_cut_300"] = 1.9; upperLimitsExpected["METStrong_3l_offZ_cut_0"] = 3.1; upperLimitsExpected["METStrong_3l_offZ_cut_50"] = 2.4; upperLimitsExpected["METStrong_3l_offZ_cut_75"] = 2.3; upperLimitsExpected["METStrong_2ltau_offZ_cut_0"] = 4.8; upperLimitsExpected["METStrong_2ltau_offZ_cut_50"] = 3.3; upperLimitsExpected["METStrong_2ltau_offZ_cut_75"] = 2.1; upperLimitsExpected["METStrong_3l_onZ_cut_20"] = 8.7; upperLimitsExpected["METStrong_3l_onZ_cut_50"] = 4.9; upperLimitsExpected["METStrong_3l_onZ_cut_75"] = 3.8; upperLimitsExpected["METStrong_2ltau_onZ_cut_20"] = 7.3; upperLimitsExpected["METStrong_2ltau_onZ_cut_50"] = 2.8; upperLimitsExpected["METStrong_2ltau_onZ_cut_75"] = 1.5; upperLimitsExpected["METWeak_3l_offZ_cut_0"] = 10.; upperLimitsExpected["METWeak_3l_offZ_cut_50"] = 4.7; upperLimitsExpected["METWeak_3l_offZ_cut_75"] = 3.0; upperLimitsExpected["METWeak_2ltau_offZ_cut_0"] = 21.; upperLimitsExpected["METWeak_2ltau_offZ_cut_50"] = 4.0; upperLimitsExpected["METWeak_2ltau_offZ_cut_75"] = 2.6; upperLimitsExpected["METWeak_3l_onZ_cut_20"] = 30.; upperLimitsExpected["METWeak_3l_onZ_cut_50"] = 10.; upperLimitsExpected["METWeak_3l_onZ_cut_75"] = 5.4; upperLimitsExpected["METWeak_2ltau_onZ_cut_20"] = 88.; upperLimitsExpected["METWeak_2ltau_onZ_cut_50"] = 5.5; upperLimitsExpected["METWeak_2ltau_onZ_cut_75"] = 2.2; upperLimitsExpected["Meff_3l_offZ_cut_0"] = 11.; upperLimitsExpected["Meff_3l_offZ_cut_150"] = 8.8; upperLimitsExpected["Meff_3l_offZ_cut_300"] = 3.7; upperLimitsExpected["Meff_3l_offZ_cut_500"] = 2.1; upperLimitsExpected["Meff_2ltau_offZ_cut_0"] = 23.; upperLimitsExpected["Meff_2ltau_offZ_cut_150"] = 13.; upperLimitsExpected["Meff_2ltau_offZ_cut_300"] = 4.9; upperLimitsExpected["Meff_2ltau_offZ_cut_500"] = 2.4; upperLimitsExpected["Meff_3l_onZ_cut_0"] = 33.; upperLimitsExpected["Meff_3l_onZ_cut_150"] = 25.; upperLimitsExpected["Meff_3l_onZ_cut_300"] = 9.; upperLimitsExpected["Meff_3l_onZ_cut_500"] = 3.9; upperLimitsExpected["Meff_2ltau_onZ_cut_0"] = 94.; upperLimitsExpected["Meff_2ltau_onZ_cut_150"] = 35.; upperLimitsExpected["Meff_2ltau_onZ_cut_300"] = 6.8; upperLimitsExpected["Meff_2ltau_onZ_cut_500"] = 2.5; upperLimitsExpected["MeffStrong_3l_offZ_cut_0"] = 3.9; upperLimitsExpected["MeffStrong_3l_offZ_cut_150"] = 3.9; upperLimitsExpected["MeffStrong_3l_offZ_cut_300"] = 3.0; upperLimitsExpected["MeffStrong_3l_offZ_cut_500"] = 2.0; upperLimitsExpected["MeffStrong_2ltau_offZ_cut_0"] = 3.8; upperLimitsExpected["MeffStrong_2ltau_offZ_cut_150"] = 3.9; upperLimitsExpected["MeffStrong_2ltau_offZ_cut_300"] = 3.1; upperLimitsExpected["MeffStrong_2ltau_offZ_cut_500"] = 1.6; upperLimitsExpected["MeffStrong_3l_onZ_cut_0"] = 6.9; upperLimitsExpected["MeffStrong_3l_onZ_cut_150"] = 7.1; upperLimitsExpected["MeffStrong_3l_onZ_cut_300"] = 4.9; upperLimitsExpected["MeffStrong_3l_onZ_cut_500"] = 3.0; upperLimitsExpected["MeffStrong_2ltau_onZ_cut_0"] = 2.4; upperLimitsExpected["MeffStrong_2ltau_onZ_cut_150"] = 2.5; upperLimitsExpected["MeffStrong_2ltau_onZ_cut_300"] = 2.0; upperLimitsExpected["MeffStrong_2ltau_onZ_cut_500"] = 1.1; if (observed) return upperLimitsObserved[signal_region]; else return upperLimitsExpected[signal_region]; }
Access the controlling AnalysisHandler object. Definition at line 288 of file Analysis.hh. References Analysis::_analysishandler. Referenced by Analysis::beamIds(), Analysis::beams(), Analysis::histoDir(), Analysis::numEvents(), Analysis::sqrtS(), and Analysis::sumOfWeights(). { return *_analysishandler; }
Get the canonical histogram "directory" path for this analysis.
Definition at line 40 of file Analysis.cc. References Analysis::handler(), Analysis::name(), and AnalysisHandler::runName(). Referenced by Analysis::histoPath().
Get the canonical histogram path for the named histogram in this analysis. Definition at line 56 of file Analysis.cc. References Analysis::histoDir(). Referenced by Analysis::bookHisto1D(), Analysis::bookHisto2D(), Analysis::bookProfile1D(), Analysis::bookProfile2D(), Analysis::bookScatter2D(), ATLAS_2012_I1094061::HistoPair::finalize(), H1_1995_S3167097::finalize(), and Analysis::getAnalysisObject(). { const string path = histoDir() + "/" + hname; return path; }
Get the canonical histogram path for the numbered histogram in this analysis. Definition at line 62 of file Analysis.cc. References Analysis::histoDir(), and Analysis::makeAxisCode(). { return histoDir() + "/" + makeAxisCode(datasetId, xAxisId, yAxisId); }
Get the actual AnalysisInfo object in which all this metadata is stored. Definition at line 104 of file Analysis.hh. References Analysis::_info. Referenced by Analysis::authors(), Analysis::bibKey(), Analysis::bibTeX(), Analysis::collider(), Analysis::description(), Analysis::experiment(), Analysis::inspireId(), Analysis::name(), Analysis::needsCrossSection(), Analysis::references(), Analysis::requiredBeams(), Analysis::requiredEnergies(), Analysis::runInfo(), Analysis::setNeedsCrossSection(), Analysis::setRequiredBeams(), Analysis::setRequiredEnergies(), Analysis::spiresId(), Analysis::status(), Analysis::summary(), Analysis::todos(), and Analysis::year().
Get the actual AnalysisInfo object in which all this metadata is stored (non-const). Definition at line 246 of file Analysis.hh. References Analysis::_info.
Book histograms and initialise projections before the run. Reimplemented from Analysis. Definition at line 23 of file ATLAS_2012_I1204447.cc. References ATLAS_2012_I1204447::_eventCountsPerSR, ATLAS_2012_I1204447::_h_e_n, ATLAS_2012_I1204447::_h_excluded, ATLAS_2012_I1204447::_h_HTjets_all, ATLAS_2012_I1204447::_h_HTlep_all, ATLAS_2012_I1204447::_h_Meff_all, ATLAS_2012_I1204447::_h_MET_all, ATLAS_2012_I1204447::_h_mu_n, ATLAS_2012_I1204447::_h_pt_1_2ltau, ATLAS_2012_I1204447::_h_pt_1_3l, ATLAS_2012_I1204447::_h_pt_2_2ltau, ATLAS_2012_I1204447::_h_pt_2_3l, ATLAS_2012_I1204447::_h_pt_3_2ltau, ATLAS_2012_I1204447::_h_pt_3_3l, ATLAS_2012_I1204447::_h_tau_n, ATLAS_2012_I1204447::_signal_regions, ATLAS_2012_I1204447::_use_fiducial_lepton_efficiency, Rivet::Cuts::abseta, IdentifiedFinalState::acceptIdPair(), ProjectionApplier::addProjection(), VetoedFinalState::addVetoPairId(), FastJets::ANTIKT, Analysis::bookHisto1D(), Rivet::PID::ELECTRON, ATLAS_2012_I1204447::getSignalRegions(), Rivet::GeV, and Rivet::PID::MUON. { // To calculate the acceptance without having the fiducial lepton efficiencies included, this part can be turned off _use_fiducial_lepton_efficiency = true; // Random numbers for simulation of ATLAS detector reconstruction efficiency srand(160385); // Read in all signal regions _signal_regions = getSignalRegions(); // Set number of events per signal region to 0 for (size_t i = 0; i < _signal_regions.size(); i++) _eventCountsPerSR[_signal_regions[i]] = 0.0; // Final state including all charged and neutral particles const FinalState fs(-5.0, 5.0, 1*GeV); addProjection(fs, "FS"); // Final state including all charged particles addProjection(ChargedFinalState(Cuts::abseta < 2.5 && Cuts::pT > 1*GeV), "CFS"); // Final state including all visible particles (to calculate MET, Jets etc.) addProjection(VisibleFinalState(Cuts::abseta < 5.0), "VFS"); // Final state including all AntiKt 04 Jets VetoedFinalState vfs; vfs.addVetoPairId(PID::MUON); addProjection(FastJets(vfs, FastJets::ANTIKT, 0.4), "AntiKtJets04"); // Final state including all unstable particles (including taus) addProjection(UnstableFinalState(Cuts::abseta < 5.0 && Cuts::pT > 5*GeV), "UFS"); // Final state including all electrons IdentifiedFinalState elecs(Cuts::abseta < 2.47 && Cuts::pT > 10*GeV); elecs.acceptIdPair(PID::ELECTRON); addProjection(elecs, "elecs"); // Final state including all muons IdentifiedFinalState muons(Cuts::abseta < 2.5 && Cuts::pT > 10*GeV); muons.acceptIdPair(PID::MUON); addProjection(muons, "muons"); // Book histograms _h_HTlep_all = bookHisto1D("HTlep_all" , 30, 0, 1500); _h_HTjets_all = bookHisto1D("HTjets_all", 30, 0, 1500); _h_MET_all = bookHisto1D("MET_all" , 20, 0, 1000); _h_Meff_all = bookHisto1D("Meff_all" , 30, 0, 3000); _h_e_n = bookHisto1D("e_n" , 10, -0.5, 9.5); _h_mu_n = bookHisto1D("mu_n" , 10, -0.5, 9.5); _h_tau_n = bookHisto1D("tau_n", 10, -0.5, 9.5); _h_pt_1_3l = bookHisto1D("pt_1_3l", 100, 0, 2000); _h_pt_2_3l = bookHisto1D("pt_2_3l", 100, 0, 2000); _h_pt_3_3l = bookHisto1D("pt_3_3l", 100, 0, 2000); _h_pt_1_2ltau = bookHisto1D("pt_1_2ltau", 100, 0, 2000); _h_pt_2_2ltau = bookHisto1D("pt_2_2ltau", 100, 0, 2000); _h_pt_3_2ltau = bookHisto1D("pt_3_2ltau", 100, 0, 2000); _h_excluded = bookHisto1D("excluded", 2, -0.5, 1.5); }
Get the Inspire ID code for this analysis. Definition at line 119 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::inspireId().
Helper for converting a differential histo to an integral one.
Definition at line 700 of file Analysis.cc.
Helper for converting a differential histo to an integral one.
Definition at line 707 of file Analysis.cc.
Check if analysis is compatible with the provided beam particle IDs and energies. Definition at line 101 of file Analysis.cc. Referenced by Analysis::isCompatible(). { return isCompatible(beams.first.pid(), beams.second.pid(), beams.first.energy(), beams.second.energy()); }
Check if analysis is compatible with the provided beam particle IDs and energies. Definition at line 107 of file Analysis.cc. References Analysis::beams(), and Analysis::isCompatible(). { PdgIdPair beams(beam1, beam2); pair<double,double> energies(e1, e2); return isCompatible(beams, energies); }
Check if analysis is compatible with the provided beam particle IDs and energies. Function checking if there is an OSSF lepton pair or a combination of 3 leptons with an invariant mass close to the Z mass
Definition at line 996 of file ATLAS_2012_I1204447.cc. References Rivet::GeV, Rivet::Cuts::mass, Particle::momentum(), and Particle::pid(). Referenced by ATLAS_2012_I1204447::analyze(). { int onZ = 0; double best_mass_2 = 999.; double best_mass_3 = 999.; // Loop over all 2 particle combinations to find invariant mass of OSSF pair closest to Z mass foreach ( const Particle& p1, particles ) { foreach ( const Particle& p2, particles ) { double mass_difference_2_old = fabs(91.0 - best_mass_2); double mass_difference_2_new = fabs(91.0 - (p1.momentum() + p2.momentum()).mass()/GeV); // If particle combination is OSSF pair calculate mass difference to Z mass if ( (p1.pid()*p2.pid() == -121 || p1.pid()*p2.pid() == -169) ) { // Get invariant mass closest to Z mass if (mass_difference_2_new < mass_difference_2_old) best_mass_2 = (p1.momentum() + p2.momentum()).mass()/GeV; // In case there is an OSSF pair take also 3rd lepton into account (e.g. from FSR and photon to electron conversion) foreach ( const Particle & p3 , particles ) { double mass_difference_3_old = fabs(91.0 - best_mass_3); double mass_difference_3_new = fabs(91.0 - (p1.momentum() + p2.momentum() + p3.momentum()).mass()/GeV); if (mass_difference_3_new < mass_difference_3_old) best_mass_3 = (p1.momentum() + p2.momentum() + p3.momentum()).mass()/GeV; } } } } // Pick the minimum invariant mass of the best OSSF pair combination and the best 3 lepton combination // If this mass is in a 20 GeV window around the Z mass, the event is classified as onZ double best_mass = min(best_mass_2, best_mass_3); if (fabs(91.0 - best_mass) < 20) onZ = 1; return onZ; }
Get the internal histogram name for given d, x and y (cf. HepData) Definition at line 67 of file Analysis.cc. Referenced by Analysis::bookHisto1D(), Analysis::bookProfile1D(), Analysis::bookScatter2D(), Analysis::getHisto1D(), Analysis::getProfile1D(), Analysis::getScatter2D(), Analysis::histoPath(), and Analysis::refData(). { stringstream axisCode; axisCode << "d"; if (datasetId < 10) axisCode << 0; axisCode << datasetId; axisCode << "-x"; if (xAxisId < 10) axisCode << 0; axisCode << xAxisId; axisCode << "-y"; if (yAxisId < 10) axisCode << 0; axisCode << yAxisId; return axisCode.str(); }
Get the name of the analysis. By default this is computed by combining the results of the experiment, year and Spires ID metadata methods and you should only override it if there's a good reason why those won't work. Implements ProjectionApplier. Definition at line 114 of file Analysis.hh. References Analysis::_defaultname, Analysis::info(), and AnalysisInfo::name(). Referenced by Analysis::_cacheRefData(), Analysis::Analysis(), Analysis::bookHisto1D(), Analysis::bookHisto2D(), Analysis::bookProfile1D(), Analysis::bookProfile2D(), Analysis::bookScatter2D(), Analysis::crossSection(), Analysis::getHisto1D(), Analysis::getLog(), Analysis::getProfile1D(), Analysis::getScatter2D(), Analysis::histoDir(), Analysis::normalize(), Analysis::refData(), and Analysis::scale(). { return (info().name().empty()) ? _defaultname : info().name(); }
Return true if this analysis needs to know the process cross-section.
Definition at line 229 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::needsCrossSection(). { return info().needsCrossSection(); }
Normalize the given histogram, histo, to area = norm.
Definition at line 652 of file Analysis.cc. References MSG_ERROR, MSG_TRACE, MSG_WARNING, and Analysis::name(). Referenced by CMS_2011_S8950903::finalize(), TOTEM_2012_002::finalize(), CMS_2011_S8968497::finalize(), CMS_2012_I1090423::finalize(), D0_2010_S8671338::finalize(), D0_2007_S7075677::finalize(), MC_JETTAGS::finalize(), CDF_2009_NOTE_9936::finalize(), D0_2008_S7554427::finalize(), MC_HINC::finalize(), MC_ZJETS::finalize(), EXAMPLE_CUTS::finalize(), MC_HJETS::finalize(), BABAR_2005_S6181155::finalize(), CMS_2013_I1265659::finalize(), MC_ZINC::finalize(), ALICE_2010_S8624100::finalize(), ATLAS_2011_S8971293::finalize(), ALICE_2010_S8625980::finalize(), CDF_1996_S3418421::finalize(), D0_2009_S8320160::finalize(), CDF_1993_S2742446::finalize(), CMS_2011_S8957746::finalize(), D0_1996_S3324664::finalize(), ATLAS_2012_I1119557::finalize(), CMS_2013_I1273574::finalize(), DELPHI_1995_S3137023::finalize(), D0_2010_S8821313::finalize(), SLD_2002_S4869273::finalize(), MC_ZZJETS::finalize(), CMS_2013_I1272853::finalize(), ATLAS_2011_I925932::finalize(), D0_1996_S3214044::finalize(), D0_2008_S6879055::finalize(), CDF_1996_S3108457::finalize(), ATLAS_2014_I1300647::finalize(), ATLAS_2012_I1204784::finalize(), EXAMPLE::finalize(), CDF_1994_S2952106::finalize(), CMS_2013_I1224539_DIJET::finalize(), JADE_1998_S3612880::finalize(), ALEPH_2001_S4656318::finalize(), DELPHI_2002_069_CONF_603::finalize(), D0_2004_S5992206::finalize(), OPAL_1997_S3396100::finalize(), MC_PHOTONS::finalize(), MC_HFJETS::finalize(), MC_GENERIC::finalize(), D0_2001_S4674421::finalize(), MC_ZZINC::finalize(), OPAL_2001_S4553896::finalize(), ATLAS_2013_I1217867::finalize(), CMS_2012_I1107658::finalize(), TASSO_1990_S2148048::finalize(), CMS_2013_I1209721::finalize(), ATLAS_2013_I1217863_Z::finalize(), STAR_2008_S7869363::finalize(), ATLAS_2013_I1217863_W::finalize(), OPAL_1998_S3780481::finalize(), CDF_1997_S3541940::finalize(), CMS_2013_I1224539_WJET::finalize(), CMS_2013_I1256943::finalize(), CMS_2013_I1224539_ZJET::finalize(), DELPHI_2003_WUD_03_11::finalize(), JADE_OPAL_2000_S4300807::finalize(), ATLAS_2014_I1298811::finalize(), ATLAS_2012_I1124167::finalize(), CMS_2014_I1303894::finalize(), H1_1994_S2919893::finalize(), CDF_2001_S4751469::finalize(), MC_TTBAR::finalize(), ALEPH_2004_S5765862::finalize(), ATLAS_2013_I1190187::finalize(), CDF_2004_S5839831::finalize(), CDF_1996_S3349578::finalize(), ATLAS_2012_I1094564::finalize(), ATLAS_2012_I1203852::finalize(), ALEPH_1996_S3486095::finalize(), DELPHI_1996_S3430090::finalize(), ATLAS_2014_I1279489::finalizePlots(), CMS_2013_I1258128::normalizeByContents(), and ATLAS_2011_I919017::safeinvscale(). { if (!histo) { MSG_ERROR("Failed to normalize histo=NULL in analysis " << name() << " (norm=" << norm << ")"); return; } MSG_TRACE("Normalizing histo " << histo->path() << " to " << norm); try { histo->normalize(norm, includeoverflows); } catch (YODA::Exception& we) { MSG_WARNING("Could not normalize histo " << histo->path()); return; } }
Normalize the given histogram, histo, to area = norm.
Get the number of events seen (via the analysis handler). Use in the finalize phase only. Definition at line 88 of file Analysis.cc. References Analysis::handler(), and AnalysisHandler::numEvents(). Referenced by STAR_2008_S7993412::analyze(), STAR_2006_S6870392::analyze(), CDF_2008_S7540469::analyze(), D0_2008_S7863608::analyze(), D0_2009_S8349509::analyze(), ATLAS_2014_I1325553::finalize(), ATLAS_2013_I1219109::finalize(), ATLAS_2013_I1217863_Z::finalize(), ATLAS_2013_I1217863_W::finalize(), ATLAS_2014_I1307756::finalize(), and ATLAS_2014_I1307243::finalize().
Get reference data for a named histo. Definition at line 183 of file Analysis.cc. References Analysis::_cacheRefData(), Analysis::_refdata, MSG_ERROR, MSG_TRACE, and Analysis::name(). Referenced by ATLAS_2011_S8994773::analyze(), CDF_2001_S4751469::analyze(), ATLAS_2010_S8894728::analyze(), Analysis::bookHisto1D(), Analysis::bookProfile1D(), Analysis::bookScatter2D(), DELPHI_2000_S4328825::finalize(), OPAL_2002_S5361494::finalize(), ATLAS_2011_S9128077::finalize(), ALEPH_2004_S5765862::finalize(), ATLAS_2010_S8894728::init(), ALICE_2011_S8909580::init(), LHCB_2013_I1208105::init(), ATLAS_2011_S9002537::init(), CMS_2011_S9088458::init(), STAR_2006_S6500200::init(), ATLAS_2011_S9128077::init(), CMS_2012_I1102908::init(), CDF_1994_S2952106::init(), ATLAS_2012_I1188891::init(), D0_2008_S7837160::init(), CMS_2013_I1218372::init(), LHCB_2011_I917009::init(), CDF_1996_S3418421::init(), LHCB_2012_I1119400::init(), ATLAS_2012_I1093734::init(), SLD_1999_S3743934::init(), ATLAS_2014_I1279489::initialisePlots(), ATLAS_2012_I1125575::initializeHistograms(), ATLAS_2012_I1094568::initializePlots(), ATLAS_2011_S9126244::initializePlots(), and Analysis::refData().
Get reference data for a numbered histo. Definition at line 194 of file Analysis.cc. References Analysis::makeAxisCode(), and Analysis::refData(). { const string hname = makeAxisCode(datasetId, xAxisId, yAxisId); return refData(hname); }
Journal, and preprint references. Definition at line 180 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::references(). { return info().references(); }
Unregister a data object from the histogram system (by name) Definition at line 725 of file Analysis.cc. References Analysis::_analysisobjects. Referenced by ATLAS_2014_I1279489::finalizeEfficiencies(). { for (vector<AnalysisObjectPtr>::iterator it = _analysisobjects.begin(); it != _analysisobjects.end(); ++it) { if ((*it)->path() == path) { _analysisobjects.erase(it); break; } } }
Unregister a data object from the histogram system (by pointer) Definition at line 734 of file Analysis.cc. References Analysis::_analysisobjects. { for (vector<AnalysisObjectPtr>::iterator it = _analysisobjects.begin(); it != _analysisobjects.end(); ++it) { if (*it == ao) { _analysisobjects.erase(it); break; } } }
Return the allowed pairs of incoming beams required by this analysis. Definition at line 206 of file Analysis.hh. References AnalysisInfo::beams(), and Analysis::info().
Sets of valid beam energy pairs, in GeV. Definition at line 217 of file Analysis.hh. References AnalysisInfo::energies(), and Analysis::info(). Referenced by Analysis::setRequiredEnergies().
Information about the events needed as input for this analysis. Event types, energies, kinematic cuts, particles to be considered stable, etc. etc. Should be treated as a restructuredText bullet list (http://docutils.sourceforge.net/rst.html) Definition at line 160 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::runInfo().
Multiplicatively scale the given histogram, histo, by factor scale.
Definition at line 667 of file Analysis.cc. References MSG_ERROR, MSG_TRACE, MSG_WARNING, and Analysis::name(). Referenced by MC_JetSplittings::finalize(), MC_JetAnalysis::finalize(), MC_ParticleAnalysis::finalize(), CDF_2006_S6450792::finalize(), ALEPH_1991_S2435284::finalize(), D0_2000_S4480767::finalize(), TOTEM_2012_I1115294::finalize(), E735_1998_S3905616::finalize(), D0_2010_S8671338::finalize(), ATLAS_2010_S8591806::finalize(), CDF_2001_S4563131::finalize(), MC_XS::finalize(), UA5_1987_S1640666::finalize(), CDF_2007_S7057202::finalize(), CDF_2000_S4155203::finalize(), CDF_2008_S8093652::finalize(), BELLE_2013_I1216515::finalize(), BABAR_2007_S6895344::finalize(), CDF_2000_S4266730::finalize(), CMS_2011_S8941262::finalize(), BELLE_2008_I786560::finalize(), CDF_1998_S3618439::finalize(), CMS_2012_I1184941::finalize(), LHCB_2012_I1208102::finalize(), CDF_1990_S2089246::finalize(), CMS_2012_I1087342::finalize(), CMS_2010_S8656010::finalize(), CMS_2011_I954992::finalize(), CDF_1988_S1865951::finalize(), ATLAS_2012_I1118269::finalize(), MC_WJETS::finalize(), UA5_1982_S875503::finalize(), BABAR_2005_S6181155::finalize(), ATLAS_2011_I894867::finalize(), CMSTOTEM_2014_I1294140::finalize(), CDF_2009_S8383952::finalize(), STAR_2006_S6870392::finalize(), CMS_QCD_10_024::finalize(), LHCB_2010_I867355::finalize(), BELLE_2001_S4598261::finalize(), OPAL_1995_S3198391::finalize(), CMS_2010_S8547297::finalize(), OPAL_1994_S2927284::finalize(), OPAL_1997_S3608263::finalize(), OPAL_2000_S4418603::finalize(), CDF_2009_S8436959::finalize(), CMS_2012_I1193338::finalize(), LHCB_2011_I919315::finalize(), ALICE_2010_S8625980::finalize(), D0_2011_I895662::finalize(), MC_PDFS::finalize(), ALEPH_2002_S4823664::finalize(), ALICE_2010_S8706239::finalize(), OPAL_1996_S3257789::finalize(), UA5_1989_S1926373::finalize(), ATLAS_2014_I1288706::finalize(), CDF_2005_S6080774::finalize(), DELPHI_1999_S3960137::finalize(), ATLAS_2014_I1282441::finalize(), ALICE_2011_S8945144::finalize(), MC_IDENTIFIED::finalize(), OPAL_1998_S3702294::finalize(), BABAR_2013_I1238276::finalize(), D0_2006_S6438750::finalize(), MC_DIPHOTON::finalize(), CMS_2013_I1273574::finalize(), ALICE_2011_S8909580::finalize(), STAR_2006_S6500200::finalize(), UA5_1986_S1583476::finalize(), MC_PHOTONINC::finalize(), ATLAS_2011_S9131140::finalize(), CMS_2011_S8978280::finalize(), D0_2000_I499943::finalize(), H1_1995_S3167097::finalize(), BABAR_2003_I593379::finalize(), PDG_TAUS::finalize(), ALICE_2012_I1181770::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7662670::finalize(), MC_WINC::finalize(), ARGUS_1993_S2653028::finalize(), CLEO_2004_S5809304::finalize(), ATLAS_2011_S9035664::finalize(), D0_2009_S8202443::finalize(), D0_2008_S6879055::finalize(), CMS_2011_S9215166::finalize(), MC_PHOTONJETS::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ATLAS_2011_I954993::finalize(), MC_WWJETS::finalize(), ALEPH_1999_S4193598::finalize(), BABAR_2007_S7266081::finalize(), CDF_2009_S8233977::finalize(), CMS_2014_I1298810::finalize(), SLD_1996_S3398250::finalize(), CMS_2013_I1261026::finalize(), ATLAS_2011_I930220::finalize(), ATLAS_2013_I1219109::finalize(), SFM_1984_S1178091::finalize(), CMS_2012_I941555::finalize(), OPAL_1998_S3749908::finalize(), D0_2009_S8349509::finalize(), D0_2010_S8570965::finalize(), CMS_2011_S8973270::finalize(), ATLAS_2012_I946427::finalize(), ATLAS_2012_I1082009::finalize(), CDF_2008_S7541902::finalize(), ARGUS_1993_S2669951::finalize(), ATLAS_2010_S8919674::finalize(), UA1_1990_S2044935::finalize(), ATLAS_2011_I921594::finalize(), ATLAS_2013_I1263495::finalize(), MC_WPOL::finalize(), MC_WWINC::finalize(), CMS_2012_I1107658::finalize(), ATLAS_2014_I1312627::finalize(), CDF_2006_S6653332::finalize(), ATLAS_2013_I1217863_Z::finalize(), ALEPH_1996_S3196992::finalize(), OPAL_1993_S2692198::finalize(), D0_2008_S7719523::finalize(), ATLAS_2013_I1217863_W::finalize(), CDF_2008_S8095620::finalize(), OPAL_1998_S3780481::finalize(), STAR_2006_S6860818::finalize(), CMS_2013_I1218372::finalize(), ATLAS_2012_I1091481::finalize(), ATLAS_2011_S9128077::finalize(), ATLAS_2011_I926145::finalize(), LHCB_2014_I1281685::finalize(), ATLAS_2010_S8918562::finalize(), JADE_OPAL_2000_S4300807::finalize(), ATLAS_2014_I1306294::finalize(), ATLAS_2011_S9108483::finalize(), CDF_2008_S7540469::finalize(), ATLAS_2012_I1199269::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2012_I1083318::finalize(), ATLAS_2014_I1304688::finalize(), ATLAS_2012_I1183818::finalize(), ATLAS_2012_CONF_2012_105::finalize(), OPAL_2004_S6132243::finalize(), ARGUS_1993_S2789213::finalize(), ATLAS_2012_CONF_2012_104::finalize(), H1_1994_S2919893::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2012_I1186556::finalize(), ATLAS_2012_I1084540::finalize(), MC_VH2BB::finalize(), ATLAS_2012_CONF_2012_103::finalize(), ATLAS_2011_I944826::finalize(), H1_2000_S4129130::finalize(), ATLAS_2013_I1230812::finalize(), ATLAS_2012_I1117704::finalize(), ALEPH_2004_S5765862::finalize(), ATLAS_2011_I945498::finalize(), ATLAS_2013_I1190187::finalize(), ATLAS_2012_I1093738::finalize(), SLD_2004_S5693039::finalize(), ATLAS_2012_I1125961::finalize(), ATLAS_2011_CONF_2011_098::finalize(), ATLAS_2011_S8983313::finalize(), ATLAS_2012_I1126136::finalize(), ATLAS_2012_I1095236::finalize(), ATLAS_2012_I1190891::finalize(), ATLAS_2011_S9019561::finalize(), ATLAS_2012_I1112263::finalize(), ATLAS_2011_S9212183::finalize(), ATLAS_2012_CONF_2012_109::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2011_S9225137::finalize(), ATLAS_2014_I1306615::finalize(), ATLAS_2012_CONF_2012_001::finalize(), ATLAS_2012_I1203852::finalize(), ATLAS_2012_I943401::finalize(), ATLAS_2012_I1180197::finalize(), ALEPH_1996_S3486095::finalize(), ATLAS_2012_CONF_2012_153::finalize(), DELPHI_1996_S3430090::finalize(), SLD_1999_S3743934::finalize(), ATLAS_2012_I1204447::finalize(), ATLAS_2011_S9212353::finalize(), ATLAS_2011_S9041966::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), ATLAS_2012_I1125575::finalizeHistograms(), ATLAS_2011_I919017::safeinvscale(), and BinnedHistogram< T >::scale(). { if (!histo) { MSG_ERROR("Failed to scale histo=NULL in analysis " << name() << " (scale=" << scale << ")"); return; } if (std::isnan(scale) || std::isinf(scale)) { MSG_ERROR("Failed to scale histo=" << histo->path() << " in analysis: " << name() << " (invalid scale factor = " << scale << ")"); scale = 0; } MSG_TRACE("Scaling histo " << histo->path() << " by factor " << scale); try { histo->scaleW(scale); } catch (YODA::Exception& we) { MSG_WARNING("Could not scale histo " << histo->path()); return; } // // Transforming the histo into a scatter after scaling // vector<double> x, y, ex, ey; // for (size_t i = 0, N = histo->numBins(); i < N; ++i) { // x.push_back( histo->bin(i).midpoint() ); // ex.push_back(histo->bin(i).width()*0.5); // y.push_back(histo->bin(i).height()*scale); // ey.push_back(histo->bin(i).heightErr()*scale); // } // string title = histo->title(); // Scatter2DPtr dps( new Scatter2D(x, y, ex, ey, hpath, title) ); // addAnalysisObject(dps); }
Multiplicatively scale the given histogram, histo, by factor scale.
Set the cross section from the generator. Definition at line 149 of file Analysis.cc. References Analysis::_crossSection, and Analysis::_gotCrossSection. { _crossSection = xs; _gotCrossSection = true; return *this; }
Declare whether this analysis needs to know the process cross-section from the generator.
Definition at line 234 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::setNeedsCrossSection(). Referenced by ATLAS_2011_I954993::ATLAS_2011_I954993(), ATLAS_2014_I1306294::ATLAS_2014_I1306294(), ATLAS_2014_I1312627::ATLAS_2014_I1312627(), MC_JetAnalysis::MC_JetAnalysis(), MC_JetSplittings::MC_JetSplittings(), and MC_ParticleAnalysis::MC_ParticleAnalysis(). { info().setNeedsCrossSection(needed); return *this; }
Declare the allowed pairs of incoming beams required by this analysis. Definition at line 210 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::setBeams().
Declare the list of valid beam energy pairs, in GeV. Definition at line 221 of file Analysis.hh. References Analysis::info(), Analysis::requiredEnergies(), and AnalysisInfo::setEnergies(). { info().setEnergies(requiredEnergies); return *this; }
Get the SPIRES ID code for this analysis (~deprecated). Definition at line 124 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::spiresId().
Centre of mass energy for this run. Definition at line 27 of file Analysis.cc. References Analysis::handler(), and AnalysisHandler::sqrtS(). Referenced by ATLAS_2011_I894867::analyze(), PDG_HADRON_MULTIPLICITIES::analyze(), CMS_2012_I1193338::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), TOTEM_2012_002::analyze(), CMS_2012_I1184941::analyze(), SLD_2004_S5693039::analyze(), CMS_2011_I954992::analyze(), ALICE_2012_I1181770::analyze(), CMS_2011_S9215166::analyze(), SFM_1984_S1178091::analyze(), ALICE_2010_S8625980::analyze(), UA1_1990_S2044935::analyze(), CMS_2013_I1218372::analyze(), JADE_1998_S3612880::analyze(), TASSO_1990_S2148048::analyze(), ATLAS_2013_I1190187::analyze(), MC_TTBAR::analyze(), ATLAS_2010_S8894728::analyze(), ATLAS_2010_S8918562::analyze(), ALEPH_2004_S5765862::analyze(), CDF_2004_S5839831::analyze(), ATLAS_2012_I1203852::analyze(), ALICE_2010_S8625980::finalize(), DELPHI_2000_S4328825::finalize(), OPAL_2002_S5361494::finalize(), JADE_1998_S3612880::finalize(), UA1_1990_S2044935::finalize(), ATLAS_2010_S8918562::finalize(), ATLAS_2012_I1093734::finalize(), ALEPH_2004_S5765862::finalize(), CDF_2004_S5839831::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), ATLAS_2010_S8894728::init(), CMS_2010_S8547297::init(), ALICE_2012_I1181770::init(), ATLAS_2011_S8994773::init(), CMSTOTEM_2014_I1294140::init(), CMS_2011_S9215166::init(), ATLAS_2012_I1091481::init(), CMS_2011_S8978280::init(), CDF_2012_NOTE10874::init(), CMS_QCD_10_024::init(), SFM_1984_S1178091::init(), MC_HINC::init(), UA5_1989_S1926373::init(), MC_JetSplittings::init(), MC_WWINC::init(), MC_ZZINC::init(), ATLAS_2011_I944826::init(), MC_WINC::init(), MC_DIPHOTON::init(), MC_ZINC::init(), MC_ZZJETS::init(), MC_PHOTONINC::init(), CDF_1988_S1865951::init(), UA5_1986_S1583476::init(), MC_ParticleAnalysis::init(), CDF_1990_S2089246::init(), MC_JetAnalysis::init(), MC_WWJETS::init(), UA1_1990_S2044935::init(), ALEPH_2004_S5765862::init(), JADE_1998_S3612880::init(), JADE_OPAL_2000_S4300807::init(), ALICE_2010_S8624100::init(), ALICE_2010_S8625980::init(), CMS_2013_I1218372::init(), LHCB_2012_I1119400::init(), UA5_1988_S1867512::init(), LHCB_2011_I917009::init(), MC_PHOTONS::init(), MC_WPOL::init(), TASSO_1990_S2148048::init(), ATLAS_2010_S8918562::init(), OPAL_2004_S6132243::init(), CDF_2004_S5839831::init(), ATLAS_2012_I1094061::init(), ATLAS_2012_I1093734::init(), PDG_HADRON_MULTIPLICITIES::init(), and PDG_HADRON_MULTIPLICITIES_RATIOS::init().
Whether this analysis is trusted (in any way!) Definition at line 195 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::status().
Get a short description of the analysis. Short (one sentence) description used as an index entry. Use description() to provide full descriptive paragraphs of analysis details. Definition at line 141 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::summary().
Get the sum of event weights seen (via the analysis handler). Use in the finalize phase only. Definition at line 93 of file Analysis.cc. References Analysis::handler(), and AnalysisHandler::sumOfWeights(). Referenced by Analysis::crossSectionPerEvent(), MC_JetSplittings::finalize(), MC_JetAnalysis::finalize(), MC_ParticleAnalysis::finalize(), CDF_2006_S6450792::finalize(), ALEPH_1991_S2435284::finalize(), D0_2000_S4480767::finalize(), CMS_2011_S9086218::finalize(), CDF_2001_S4563131::finalize(), CDF_2008_S7828950::finalize(), D0_2010_S8671338::finalize(), LHCF_2012_I1115479::finalize(), MC_XS::finalize(), CDF_2008_S8093652::finalize(), CDF_2000_S4155203::finalize(), BELLE_2013_I1216515::finalize(), BABAR_2007_S6895344::finalize(), CDF_2000_S4266730::finalize(), CMS_2011_S8941262::finalize(), CDF_1998_S3618439::finalize(), CMS_2012_I1184941::finalize(), CMS_2012_I1087342::finalize(), LHCB_2012_I1208102::finalize(), CMS_2010_S8656010::finalize(), CMS_2011_I954992::finalize(), ATLAS_2012_I1118269::finalize(), MC_WJETS::finalize(), BABAR_2005_S6181155::finalize(), ATLAS_2011_I894867::finalize(), CDF_2001_S4517016::finalize(), CDF_2009_S8383952::finalize(), STAR_2006_S6870392::finalize(), D0_2010_S8566488::finalize(), LHCB_2010_I867355::finalize(), OPAL_1995_S3198391::finalize(), CMS_2010_S8547297::finalize(), OPAL_1994_S2927284::finalize(), OPAL_1997_S3608263::finalize(), OPAL_2000_S4418603::finalize(), CDF_2009_S8436959::finalize(), CMS_2012_I1193338::finalize(), D0_2011_I895662::finalize(), MC_PDFS::finalize(), ALEPH_2002_S4823664::finalize(), OPAL_1996_S3257789::finalize(), CDF_2005_S6080774::finalize(), ATLAS_2014_I1288706::finalize(), ATLAS_2014_I1325553::finalize(), BABAR_2013_I1238276::finalize(), OPAL_1998_S3702294::finalize(), DELPHI_1999_S3960137::finalize(), ATLAS_2014_I1282441::finalize(), ALICE_2011_S8945144::finalize(), MC_IDENTIFIED::finalize(), D0_2006_S6438750::finalize(), MC_DIPHOTON::finalize(), CMS_2013_I1273574::finalize(), DELPHI_1995_S3137023::finalize(), ALICE_2011_S8909580::finalize(), CMS_2011_S8978280::finalize(), D0_2000_I499943::finalize(), ALICE_2012_I1181770::finalize(), CLEO_2004_S5809304::finalize(), D0_2008_S7662670::finalize(), D0_2008_S7863608::finalize(), MC_WINC::finalize(), ATLAS_2011_S9035664::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ATLAS_2011_I954993::finalize(), MC_WWJETS::finalize(), CDF_2009_S8233977::finalize(), CMS_2014_I1298810::finalize(), LHCB_2013_I1218996::finalize(), ATLAS_2013_I1219109::finalize(), OPAL_1998_S3749908::finalize(), OPAL_1997_S3396100::finalize(), D0_2010_S8570965::finalize(), CMS_2011_S8973270::finalize(), ATLAS_2012_I946427::finalize(), CDF_2008_S7541902::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2011_I921594::finalize(), ATLAS_2013_I1263495::finalize(), MC_WWINC::finalize(), ATLAS_2014_I1312627::finalize(), CDF_2006_S6653332::finalize(), ALEPH_1996_S3196992::finalize(), ATLAS_2013_I1217863_Z::finalize(), OPAL_1993_S2692198::finalize(), ATLAS_2013_I1217863_W::finalize(), D0_2008_S7719523::finalize(), OPAL_1998_S3780481::finalize(), STAR_2006_S6860818::finalize(), ATLAS_2014_I1307756::finalize(), ATLAS_2011_I926145::finalize(), CMS_2013_I1256943::finalize(), JADE_OPAL_2000_S4300807::finalize(), ATLAS_2014_I1306294::finalize(), ATLAS_2011_S9108483::finalize(), CDF_2008_S7540469::finalize(), ATLAS_2012_I1199269::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2014_I1307243::finalize(), ATLAS_2012_I1083318::finalize(), ATLAS_2014_I1304688::finalize(), ATLAS_2012_CONF_2012_105::finalize(), CMS_2014_I1303894::finalize(), ATLAS_2012_CONF_2012_104::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2012_I1186556::finalize(), ATLAS_2012_I1084540::finalize(), MC_VH2BB::finalize(), ATLAS_2012_CONF_2012_103::finalize(), ATLAS_2012_I1117704::finalize(), ALEPH_2004_S5765862::finalize(), ATLAS_2013_I1190187::finalize(), ATLAS_2012_I1093738::finalize(), SLD_2004_S5693039::finalize(), ATLAS_2012_I1125961::finalize(), ATLAS_2011_CONF_2011_098::finalize(), ATLAS_2011_S8983313::finalize(), ATLAS_2012_I1126136::finalize(), ATLAS_2012_I1095236::finalize(), ATLAS_2012_I1190891::finalize(), ATLAS_2011_S9019561::finalize(), ATLAS_2012_I1112263::finalize(), ATLAS_2011_S9212183::finalize(), ATLAS_2012_CONF_2012_109::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2011_S9225137::finalize(), ATLAS_2012_CONF_2012_001::finalize(), ATLAS_2012_I1203852::finalize(), ATLAS_2012_I943401::finalize(), ATLAS_2012_I1180197::finalize(), ALEPH_1996_S3486095::finalize(), ATLAS_2012_CONF_2012_153::finalize(), SLD_1999_S3743934::finalize(), ATLAS_2012_I1204447::finalize(), ATLAS_2011_S9212353::finalize(), ATLAS_2011_S9041966::finalize(), and PDG_HADRON_MULTIPLICITIES::finalize(). { return handler().sumOfWeights(); }
Any work to be done on this analysis. Definition at line 200 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::todos().
When the original experimental analysis was published. Definition at line 175 of file Analysis.hh. References Analysis::info(), and AnalysisInfo::year(). Member Data Documentation
Flag to forbid projection registration in analyses until the init phase. Definition at line 143 of file ProjectionApplier.hh. Referenced by ProjectionApplier::_addProjection(), and Analysis::Analysis().
Definition at line 1050 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::fillEventCountsPerSR(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1041 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1042 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1039 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Histograms. Definition at line 1039 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1039 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1039 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1041 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1040 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1040 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1040 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1040 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1040 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1040 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Definition at line 1041 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
List of signal regions and event counts per signal region. Definition at line 1049 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::finalize(), and ATLAS_2012_I1204447::init().
Fiducial efficiencies to model the effects of the ATLAS detector. Definition at line 1046 of file ATLAS_2012_I1204447.cc. Referenced by ATLAS_2012_I1204447::analyze(), and ATLAS_2012_I1204447::init(). The documentation for this class was generated from the following file: Generated on Tue Mar 24 2015 17:36:12 for The Rivet MC analysis system by ![]() |