Rivet Analyses Reference

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/LeadingParticlesFinalState.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
#include "Rivet/Projections/FastJets.hh"

namespace Rivet {


  /// @brief Measurement of gamma + jets + X triple differential cross-sections
  ///
  /// @author David Grellscheid
  class CMS_2014_I1266056 : public Analysis {
  public:

    // Constructor
    CMS_2014_I1266056()
      : Analysis("CMS_2014_I1266056")
    {    }


    // Book histograms and initialise projections before the run
    void init() {
      // Final state
      FinalState fs((Cuts::etaIn(-3, 3)));
      declare(fs, "FS");

      // Leading photon
      LeadingParticlesFinalState photonfs(FinalState((Cuts::etaIn(-2.5, 2.5) && Cuts::pT >=  40.0*GeV)));
      photonfs.addParticleId(PID::PHOTON);
      declare(photonfs, "LeadingPhoton");

      // FS excluding the leading photon
      VetoedFinalState vfs(fs);
      vfs.addVetoOnThisFinalState(photonfs);
      declare(vfs, "JetFS");

      // Jets
      FastJets jetpro(vfs, FastJets::ANTIKT, 0.5);
      //jetpro.useInvisibles();
      declare(jetpro, "Jets");

      book(_h_phverycentral_jetcentral, 1, 1, 1);
      book(_h_phcentral_jetcentral    , 2, 1, 1);
      book(_h_phforward_jetcentral    , 3, 1, 1);
      book(_h_phveryforward_jetcentral, 4, 1, 1);

      book(_h_phverycentral_jetforward, 1, 1, 2);
      book(_h_phcentral_jetforward    , 2, 1, 2);
      book(_h_phforward_jetforward    , 3, 1, 2);
      book(_h_phveryforward_jetforward, 4, 1, 2);

    }

    // Perform the per-event analysis
    void analyze(const Event& event) {

      // Get the photon
      const FinalState& photonfs = applyProjection<FinalState>(event, "LeadingPhoton");
      if (photonfs.particles().empty()) vetoEvent;
      const FourMomentum photon = photonfs.particles().front().momentum();

      // Get the jet
      Jets jets = applyProjection<FastJets>(event, "Jets").jetsByPt(30.0*GeV);
      if (jets.empty()) vetoEvent;
      FourMomentum leadingJet;
      for ( const Jet & j : jets ) {
        leadingJet = j.momentum();
        // keep the first separated jet
        if (deltaR(photon, leadingJet) > 0.5)
          break;
      }
      if (deltaR(photon, leadingJet) < 0.5)
          vetoEvent;

      // Veto if leading jet is outside plotted rapidity regions
      if (leadingJet.abseta() > 2.5) vetoEvent;

      // TODO: photon isolation 'IsoGamma' needed?

      // Fill histos
      const double abs_jet_eta = leadingJet.abseta();
      const double photon_pt = photon.pT()/GeV;
      const double abs_photon_eta = photon.abseta();

      if (abs_jet_eta < 1.5) {
        if      (abs_photon_eta < 0.9)  _h_phverycentral_jetcentral->fill(photon_pt);
        else if (abs_photon_eta < 1.44) _h_phcentral_jetcentral->fill(    photon_pt);
        else if (abs_photon_eta < 1.57) {}
        else if (abs_photon_eta < 2.1)  _h_phforward_jetcentral->fill(    photon_pt);
        else if (abs_photon_eta < 2.5)  _h_phveryforward_jetcentral->fill(photon_pt);
      }
      else if (abs_jet_eta < 2.5) {
        if      (abs_photon_eta < 0.9)  _h_phverycentral_jetforward->fill(photon_pt);
        else if (abs_photon_eta < 1.44) _h_phcentral_jetforward->fill(    photon_pt);
        else if (abs_photon_eta < 1.57) {}
        else if (abs_photon_eta < 2.1)  _h_phforward_jetforward->fill(    photon_pt);
        else if (abs_photon_eta < 2.5)  _h_phveryforward_jetforward->fill(photon_pt);
      }
    }
    


    /// Normalise histograms etc., after the run
    void finalize() {
      const double scale_jetcentral = crossSection()/sumOfWeights(); // *3 (jet eta < 1.5)
      scale(_h_phverycentral_jetcentral, scale_jetcentral); // * 1.8 (photon eta < 0.9)
      scale(_h_phcentral_jetcentral    , scale_jetcentral); // * 1.08 (0.9 .. 1.44)
      scale(_h_phforward_jetcentral    , scale_jetcentral); // * 1.06 (1.57 .. 2.1)
      scale(_h_phveryforward_jetcentral, scale_jetcentral); // * 0.8  (2.1 .. 2.5)

      const double scale_jetforward = crossSection()/sumOfWeights(); // *2 (1.5 < eta < 2.5)
      scale(_h_phverycentral_jetforward, scale_jetforward); // .. as above ..
      scale(_h_phcentral_jetforward    , scale_jetforward); // .. as above ..
      scale(_h_phforward_jetforward    , scale_jetforward); // .. as above ..
      scale(_h_phveryforward_jetforward, scale_jetforward); // .. as above ..

    }


  private:

    Histo1DPtr _h_phverycentral_jetcentral;
    Histo1DPtr _h_phcentral_jetcentral    ;
    Histo1DPtr _h_phforward_jetcentral    ;
    Histo1DPtr _h_phveryforward_jetcentral;

    Histo1DPtr _h_phverycentral_jetforward;
    Histo1DPtr _h_phcentral_jetforward    ;
    Histo1DPtr _h_phforward_jetforward    ;
    Histo1DPtr _h_phveryforward_jetforward;

  };


  // The hook for the plugin system
  RIVET_DECLARE_PLUGIN(CMS_2014_I1266056);


}