Rivet::CMS_TEST_ANALYSIS Class Reference

Inheritance diagram for Rivet::CMS_TEST_ANALYSIS:

Public Member Functions
Constructors etc.
	CMS_TEST_ANALYSIS ()
	Constructor. More...
Analysis methods
void	init ()
	Book histograms and initialise projections before the run. More...
size_t	findPtBin (double ptJ)
void	analyze (const Event &event)
	Perform the per-event analysis. More...
void	finalize ()
	Normalise histograms etc., after the run. More...

Private Attributes
FastJet grooming tools (configured in constructor init list)
const fastjet::Filter	_filter
const fastjet::Filter	_trimmer
const fastjet::Pruner	_pruner

Histograms
enum Rivet::CMS_TEST_ANALYSIS::{   PT_220_300_dj =0, PT_300_450_dj, PT_450_500_dj, PT_500_600_dj,   PT_600_800_dj, PT_800_1000_dj, PT_1000_1500_dj, N_PT_BINS_dj }	BINS_dj
AIDA::IHistogram1D *	_h_ungroomedJet0pt
AIDA::IHistogram1D *	_h_ungroomedJet1pt
AIDA::IHistogram1D *	_h_ungroomedAvgJetMass_dj [N_PT_BINS_dj]
AIDA::IHistogram1D *	_h_filteredAvgJetMass_dj [N_PT_BINS_dj]
AIDA::IHistogram1D *	_h_trimmedAvgJetMass_dj [N_PT_BINS_dj]
AIDA::IHistogram1D *	_h_prunedAvgJetMass_dj [N_PT_BINS_dj]

Detailed Description

Definition at line 17 of file CMS_TEST_ANALYSIS.cc.

Member Enumeration Documentation

anonymous enum

private

Enumerator
PT_220_300_dj
PT_300_450_dj
PT_450_500_dj
PT_500_600_dj
PT_600_800_dj
PT_800_1000_dj
PT_1000_1500_dj
N_PT_BINS_dj

Definition at line 129 of file CMS_TEST_ANALYSIS.cc.

         { PT_220_300_dj=0, PT_300_450_dj, PT_450_500_dj, PT_500_600_dj,
           PT_600_800_dj, PT_800_1000_dj, PT_1000_1500_dj, N_PT_BINS_dj } BINS_dj;

Constructor & Destructor Documentation

Rivet::CMS_TEST_ANALYSIS::CMS_TEST_ANALYSIS ( )

inline

Constructor.

Definition at line 24 of file CMS_TEST_ANALYSIS.cc.

      : Analysis("CMS_TEST_ANALYSIS"),
        _filter(fastjet::Filter(fastjet::JetDefinition(fastjet::cambridge_algorithm, 0.3), fastjet::SelectorNHardest(3))),
        _trimmer(fastjet::Filter(fastjet::JetDefinition(fastjet::kt_algorithm, 0.2), fastjet::SelectorPtFractionMin(0.03))),
        _pruner(fastjet::Pruner(fastjet::cambridge_algorithm, 0.1, 0.5))
    {    }

Member Function Documentation

void Rivet::CMS_TEST_ANALYSIS::analyze ( const Event &  event )

inline

Perform the per-event analysis.

Definition at line 71 of file CMS_TEST_ANALYSIS.cc.

References _filter, _h_filteredAvgJetMass_dj, _h_prunedAvgJetMass_dj, _h_trimmedAvgJetMass_dj, _h_ungroomedAvgJetMass_dj, _pruner, _trimmer, event(), findPtBin(), N_PT_BINS_dj, and CommonMethods::weight().

                                     {
      const double weight = event.weight();

      // Look at events with >= 2 jets
      const PseudoJets& psjetsAK7 = applyProjection<FastJets>(event, "JetsAK7").pseudoJetsByPt( 50.0*GeV );
      if (psjetsAK7.size() < 2) vetoEvent;

      // Get the leading two jets and find their average pT
      const fastjet::PseudoJet& j0 = psjetsAK7[0];
      const fastjet::PseudoJet& j1 = psjetsAK7[1];
      double ptAvg = 0.5 * (j0.pt() + j1.pt());

      // Find the appropriate mean pT bin and escape if needed
      const size_t njetBin = findPtBin(ptAvg/GeV);
      if (njetBin >= N_PT_BINS_dj) vetoEvent;

      // Now run the substructure algs...
      fastjet::PseudoJet filtered0 = _filter(j0);
      fastjet::PseudoJet filtered1 = _filter(j1);
      fastjet::PseudoJet trimmed0 = _trimmer(j0);
      fastjet::PseudoJet trimmed1 = _trimmer(j1);
      fastjet::PseudoJet pruned0 = _pruner(j0);
      fastjet::PseudoJet pruned1 = _pruner(j1);

      // ... and fill the histograms
      _h_ungroomedAvgJetMass_dj[njetBin]->fill(0.5*(j0.m() + j1.m())/GeV, weight);
      _h_filteredAvgJetMass_dj[njetBin]->fill(0.5*(filtered0.m() + filtered1.m())/GeV, weight);
      _h_trimmedAvgJetMass_dj[njetBin]->fill(0.5*(trimmed0.m() + trimmed1.m())/GeV, weight);
      _h_prunedAvgJetMass_dj[njetBin]->fill(0.5*(pruned0.m() + pruned1.m())/GeV, weight);
    }

void Rivet::CMS_TEST_ANALYSIS::finalize ( void  )

inline

Normalise histograms etc., after the run.

Definition at line 104 of file CMS_TEST_ANALYSIS.cc.

References _h_filteredAvgJetMass_dj, _h_prunedAvgJetMass_dj, _h_trimmedAvgJetMass_dj, _h_ungroomedAvgJetMass_dj, i, and N_PT_BINS_dj.

                    {
      const double normalizationVal = 1000;
      for (size_t i = 0; i < N_PT_BINS_dj; ++i) {
        normalize(_h_ungroomedAvgJetMass_dj[i], normalizationVal);
        normalize(_h_filteredAvgJetMass_dj[i], normalizationVal);
        normalize(_h_trimmedAvgJetMass_dj[i], normalizationVal);
        normalize(_h_prunedAvgJetMass_dj[i], normalizationVal);
      }
    }

size_t Rivet::CMS_TEST_ANALYSIS::findPtBin ( double  ptJ )

inline

Definition at line 61 of file CMS_TEST_ANALYSIS.cc.

References N_PT_BINS_dj.

Referenced by analyze().

                                 {
      const double ptBins_dj[N_PT_BINS_dj+1] = { 220.0, 300.0, 450.0, 500.0, 600.0, 800.0, 1000.0, 1500.0};
      for (size_t ibin = 0; ibin < N_PT_BINS_dj; ++ibin) {
        if (inRange(ptJ, ptBins_dj[ibin], ptBins_dj[ibin+1])) return ibin;
      }
      return N_PT_BINS_dj;
    }

void Rivet::CMS_TEST_ANALYSIS::init ( void  )

inline

Book histograms and initialise projections before the run.

Definition at line 40 of file CMS_TEST_ANALYSIS.cc.

References _h_filteredAvgJetMass_dj, _h_prunedAvgJetMass_dj, _h_trimmedAvgJetMass_dj, _h_ungroomedAvgJetMass_dj, i, and N_PT_BINS_dj.

                {
      FinalState fs(-2.4, 2.4, 0*GeV);
      addProjection(fs, "FS");

      // Jet collections
      addProjection(FastJets(fs, FastJets::ANTIKT, 0.7), "JetsAK7");
      addProjection(FastJets(fs, FastJets::CAM, 0.8), "JetsCA8");
      addProjection(FastJets(fs, FastJets::CAM, 1.2), "JetsCA12");

      // Histograms
      for (size_t i = 0; i < N_PT_BINS_dj; ++i ) {
        _h_ungroomedAvgJetMass_dj[i] = bookHistogram1D(i+1+0*N_PT_BINS_dj, 1, 1);
        _h_filteredAvgJetMass_dj[i] = bookHistogram1D(i+1+1*N_PT_BINS_dj, 1, 1);
        _h_trimmedAvgJetMass_dj[i] = bookHistogram1D(i+1+2*N_PT_BINS_dj, 1, 1);
        _h_prunedAvgJetMass_dj[i] = bookHistogram1D(i+1+3*N_PT_BINS_dj, 1, 1);
      }
    }

Member Data Documentation

const fastjet::Filter Rivet::CMS_TEST_ANALYSIS::_filter

private

Definition at line 121 of file CMS_TEST_ANALYSIS.cc.

Referenced by analyze().

AIDA::IHistogram1D* Rivet::CMS_TEST_ANALYSIS::_h_filteredAvgJetMass_dj[N_PT_BINS_dj]

private

Definition at line 133 of file CMS_TEST_ANALYSIS.cc.

Referenced by analyze(), finalize(), and init().

AIDA::IHistogram1D* Rivet::CMS_TEST_ANALYSIS::_h_prunedAvgJetMass_dj[N_PT_BINS_dj]

private

Definition at line 135 of file CMS_TEST_ANALYSIS.cc.

Referenced by analyze(), finalize(), and init().

AIDA::IHistogram1D* Rivet::CMS_TEST_ANALYSIS::_h_trimmedAvgJetMass_dj[N_PT_BINS_dj]

private

Definition at line 134 of file CMS_TEST_ANALYSIS.cc.

Referenced by analyze(), finalize(), and init().

AIDA::IHistogram1D* Rivet::CMS_TEST_ANALYSIS::_h_ungroomedAvgJetMass_dj[N_PT_BINS_dj]

private

Definition at line 132 of file CMS_TEST_ANALYSIS.cc.

Referenced by analyze(), finalize(), and init().

AIDA::IHistogram1D* Rivet::CMS_TEST_ANALYSIS::_h_ungroomedJet0pt

private

Definition at line 131 of file CMS_TEST_ANALYSIS.cc.

AIDA::IHistogram1D * Rivet::CMS_TEST_ANALYSIS::_h_ungroomedJet1pt

private

Definition at line 131 of file CMS_TEST_ANALYSIS.cc.

const fastjet::Pruner Rivet::CMS_TEST_ANALYSIS::_pruner

private

Definition at line 123 of file CMS_TEST_ANALYSIS.cc.

Referenced by analyze().

const fastjet::Filter Rivet::CMS_TEST_ANALYSIS::_trimmer

private

Definition at line 122 of file CMS_TEST_ANALYSIS.cc.

Referenced by analyze().

enum { ... } Rivet::CMS_TEST_ANALYSIS::BINS_dj