← Back to branch summary

~ubuntu-branches/ubuntu/saucy/rivet/saucy-proposed

« back to all changes in this revision

Viewing changes to src/Analyses/MC_WWJETS.cc

  • Committer: Package Import Robot
  • Author(s): Lifeng Sun
  • Date: 2013-05-07 09:24:27 UTC
  • mfrom: (1.2.1) (2.1.2 experimental)
  • Revision ID: package-import@ubuntu.com-20130507092427-rpfijl4mn3y87ek7
Tags: 1.8.3-1
New upstream release.

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/Analyses/MC_WWJETS.cc
38
38
      jetinput
39
39
          .addVetoOnThisFinalState(wenufinder)
40
40
          .addVetoOnThisFinalState(wmnufinder);
41
 
      FastJets jetpro(jetinput, FastJets::KT, 0.7);
 
41
      FastJets jetpro(jetinput, FastJets::ANTIKT, 0.4);
42
42
      addProjection(jetpro, "Jets");
43
43
 
44
 
      // properties of the pair momentum
45
 
      _h_WW_pT = bookHistogram1D("WW_pT", logBinEdges(100, 1.0, 0.5*sqrtS()));
46
 
      _h_WW_pT_peak = bookHistogram1D("WW_pT_peak", 25, 0.0, 25.0);
47
 
      _h_WW_eta = bookHistogram1D("WW_eta", 40, -7.0, 7.0);
48
 
      _h_WW_phi = bookHistogram1D("WW_phi", 25, 0.0, TWOPI);
49
 
      _h_WW_m = bookHistogram1D("WW_m", logBinEdges(100, 150.0, 180.0+0.25*sqrtS()));
50
 
 
51
 
      // correlations between the WW
52
 
      _h_WW_dphi = bookHistogram1D("WW_dphi", 25, 0.0, PI);  /// @todo non-linear?
53
 
      _h_WW_deta = bookHistogram1D("WW_deta", 25, -7.0, 7.0);
54
 
      _h_WW_dR = bookHistogram1D("WW_dR", 25, 0.5, 7.0);
55
 
      _h_WW_dpT = bookHistogram1D("WW_dpT", logBinEdges(100, 1.0, 0.5*sqrtS()));
56
 
      _h_WW_costheta_planes = bookHistogram1D("WW_costheta_planes", 25, -1.0, 1.0);
57
 
 
58
 
      /// @todo fuer WW: missing ET
59
 
 
60
 
      // properties of the W bosons
61
 
      _h_W_pT = bookHistogram1D("W_pT", logBinEdges(100, 10.0, 0.25*sqrtS()));
62
 
      _h_W_eta = bookHistogram1D("W_eta", 70, -7.0, 7.0);
63
 
 
64
 
      // properties of the leptons
65
 
      _h_Wl_pT = bookHistogram1D("Wl_pT", logBinEdges(100, 30.0, 0.1
66
 
                                                      *sqrtS()));
67
 
      _h_Wl_eta = bookHistogram1D("Wl_eta", 40, -3.5, 3.5);
68
 
 
69
 
      // correlations between the opposite charge leptons
70
 
      _h_WeWm_dphi = bookHistogram1D("WeWm_dphi", 25, 0.0, PI);
71
 
      _h_WeWm_deta = bookHistogram1D("WeWm_deta", 25, -5.0, 5.0);
72
 
      _h_WeWm_dR = bookHistogram1D("WeWm_dR", 25, 0.5, 5.0);
73
 
      _h_WeWm_m = bookHistogram1D("WeWm_m", 100, 0.0, 300.0);
74
 
 
75
44
      // correlations with jets
76
45
      _h_WW_jet1_deta = bookHistogram1D("WW_jet1_deta", 70, -7.0, 7.0);
77
46
      _h_WW_jet1_dR = bookHistogram1D("WW_jet1_dR", 25, 1.5, 7.0);
109
78
      FourMomentum mm=wmnufinder.constituentLeptons()[0].momentum();
110
79
      FourMomentum mnu=wmnufinder.constituentNeutrinos()[0].momentum();
111
80
 
112
 
      _h_WW_pT->fill(ww.pT(),weight);
113
 
      _h_WW_pT_peak->fill(ww.pT(),weight);
114
 
      _h_WW_eta->fill(ww.eta(),weight);
115
 
      _h_WW_phi->fill(ww.azimuthalAngle(),weight);
116
 
      double mww2=ww.mass2();
117
 
      if (mww2>0.0) _h_WW_m->fill(sqrt(mww2), weight);
118
 
 
119
 
      _h_WW_dphi->fill(mapAngle0ToPi(wenu.phi()-wmnu.phi()), weight);
120
 
      _h_WW_deta->fill(wenu.eta()-wmnu.eta(), weight);
121
 
      _h_WW_dR->fill(deltaR(wenu,wmnu), weight);
122
 
      _h_WW_dpT->fill(fabs(wenu.pT()-wmnu.pT()), weight);
123
 
 
124
 
      Vector3 crossWenu = ep.vector3().cross(enu.vector3());
125
 
      Vector3 crossWmnu = mm.vector3().cross(mnu.vector3());
126
 
      double costheta = crossWenu.dot(crossWmnu)/crossWenu.mod()/crossWmnu.mod();
127
 
      _h_WW_costheta_planes->fill(costheta, weight);
128
 
 
129
 
      _h_W_pT->fill(wenu.pT(),weight);
130
 
      _h_W_pT->fill(wmnu.pT(),weight);
131
 
      _h_W_eta->fill(wenu.eta(),weight);
132
 
      _h_W_eta->fill(wmnu.eta(),weight);
133
 
 
134
 
      _h_Wl_pT->fill(ep.pT(), weight);
135
 
      _h_Wl_pT->fill(mm.pT(), weight);
136
 
      _h_Wl_eta->fill(ep.eta(), weight);
137
 
      _h_Wl_eta->fill(mm.eta(), weight);
138
 
 
139
 
      _h_WeWm_dphi->fill(mapAngle0ToPi(ep.phi()-mm.phi()), weight);
140
 
      _h_WeWm_deta->fill(ep.eta()-mm.eta(), weight);
141
 
      _h_WeWm_dR->fill(deltaR(ep,mm), weight);
142
 
      double m2=FourMomentum(ep+mm).mass2();
143
 
      if (m2 < 0) m2 = 0.0;
144
 
      _h_WeWm_m->fill(sqrt(m2), weight);
145
 
 
146
 
      const FastJets& jetpro = applyProjection<FastJets>(e, "Jets");
147
 
      const Jets& jets = jetpro.jetsByPt(20.0*GeV);
 
81
      const Jets& jets = applyProjection<FastJets>(e, "Jets").jetsByPt(m_jetptcut);
148
82
      if (jets.size() > 0) {
149
83
        _h_WW_jet1_deta->fill(ww.eta()-jets[0].momentum().eta(), weight);
150
84
        _h_WW_jet1_dR->fill(deltaR(ww, jets[0].momentum()), weight);
170
104
    /// Finalize
171
105
    void finalize() {
172
106
      double norm=crossSection()/sumOfWeights();
173
 
      scale(_h_WW_pT, norm);
174
 
      scale(_h_WW_pT_peak, norm);
175
 
      scale(_h_WW_eta, norm);
176
 
      scale(_h_WW_phi, norm);
177
 
      scale(_h_WW_m, norm);
178
 
      scale(_h_WW_dphi, norm);
179
 
      scale(_h_WW_deta, norm);
180
 
      scale(_h_WW_dR, norm);
181
 
      scale(_h_WW_dpT, norm);
182
 
      scale(_h_WW_costheta_planes, norm);
183
 
      scale(_h_W_pT, norm);
184
 
      scale(_h_W_eta, norm);
185
 
      scale(_h_Wl_pT, norm);
186
 
      scale(_h_Wl_eta, norm);
187
 
      scale(_h_WeWm_dphi, norm);
188
 
      scale(_h_WeWm_deta, norm);
189
 
      scale(_h_WeWm_dR, norm);
190
 
      scale(_h_WeWm_m, norm);
191
107
      scale(_h_WW_jet1_deta, norm);
192
108
      scale(_h_WW_jet1_dR, norm);
193
109
      scale(_h_We_jet1_dR, norm);
204
120
 
205
121
    /// @name Histograms
206
122
    //@{
207
 
    AIDA::IHistogram1D * _h_WW_pT;
208
 
    AIDA::IHistogram1D * _h_WW_pT_peak;
209
 
    AIDA::IHistogram1D * _h_WW_eta;
210
 
    AIDA::IHistogram1D * _h_WW_phi;
211
 
    AIDA::IHistogram1D * _h_WW_m;
212
 
    AIDA::IHistogram1D * _h_WW_dphi;
213
 
    AIDA::IHistogram1D * _h_WW_deta;
214
 
    AIDA::IHistogram1D * _h_WW_dR;
215
 
    AIDA::IHistogram1D * _h_WW_dpT;
216
 
    AIDA::IHistogram1D * _h_WW_costheta_planes;
217
 
    AIDA::IHistogram1D * _h_W_pT;
218
 
    AIDA::IHistogram1D * _h_W_eta;
219
 
    AIDA::IHistogram1D * _h_Wl_pT;
220
 
    AIDA::IHistogram1D * _h_Wl_eta;
221
 
    AIDA::IHistogram1D * _h_WeWm_dphi;
222
 
    AIDA::IHistogram1D * _h_WeWm_deta;
223
 
    AIDA::IHistogram1D * _h_WeWm_dR;
224
 
    AIDA::IHistogram1D * _h_WeWm_m;
225
123
    AIDA::IHistogram1D * _h_WW_jet1_deta;
226
124
    AIDA::IHistogram1D * _h_WW_jet1_dR;
227
125
    AIDA::IHistogram1D * _h_We_jet1_dR;