|
[Rivet-svn] r4083 - in branches/2012-06-aidarivet: . data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgFri Dec 14 13:38:43 GMT 2012
Author: hoeth Date: Fri Dec 14 13:38:42 2012 New Revision: 4083 Log: Added CMS_2012_PAS_FWD_11_003 analysis Added: branches/2012-06-aidarivet/data/anainfo/CMS_2012_PAS_FWD_11_003.info branches/2012-06-aidarivet/data/plotinfo/CMS_2012_PAS_FWD_11_003.plot branches/2012-06-aidarivet/data/refdata/CMS_2012_PAS_FWD_11_003.aida branches/2012-06-aidarivet/src/Analyses/CMS_2012_PAS_FWD_11_003.cc Modified: branches/2012-06-aidarivet/ChangeLog branches/2012-06-aidarivet/data/anainfo/Makefile.am branches/2012-06-aidarivet/data/plotinfo/Makefile.am branches/2012-06-aidarivet/data/refdata/Makefile.am branches/2012-06-aidarivet/src/Analyses/Makefile.am Modified: branches/2012-06-aidarivet/ChangeLog ============================================================================== --- branches/2012-06-aidarivet/ChangeLog Thu Dec 13 15:29:18 2012 (r4082) +++ branches/2012-06-aidarivet/ChangeLog Fri Dec 14 13:38:42 2012 (r4083) @@ -1,3 +1,7 @@ +2012-12-14 Hendrik Hoeth <hendrik.hoeth at cern.ch> + + * Added CMS_2012_PAS_FWD_11_003 analysis + 2012-12-12 Andy Buckley <andy.buckley at cern.ch> * Correction to jet acceptance in CMS_2011_S9120041, from Sercan Sen: thanks! Added: branches/2012-06-aidarivet/data/anainfo/CMS_2012_PAS_FWD_11_003.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/anainfo/CMS_2012_PAS_FWD_11_003.info Fri Dec 14 13:38:42 2012 (r4083) @@ -0,0 +1,22 @@ +Name: CMS_2012_PAS_FWD_11_003 +Year: 2011 +Summary: Study of the Underlying Event at Forward Rapidity in Proton-Proton Collisions at the LHC +Experiment: CMS +Collider: LHC +SpiresID: +Status: VALIDATED +Authors: + - Samantha Dooling <samantha.dooling at desy.de> +References: + - CMS-PAS-FWD-11-003 + - http://cdsweb.cern.ch/record/1434458 +RunInfo: + Inelastic events (non-diffractive and diffractive) at $\sqrt{s}$ = 0.9, 2.76 and 7 TeV. +NumEvents: 1 Mio +Beams: [p+, p+] +Energies: [900, 2760, 7000] +Description: + 'Ratio of the energy deposited in the pseudorapidity range $-6.6 < \eta < -5.2$ + for events with a charged particle jet with $|\eta|<2$ with respect to the + energy in inclusive events, as a function of charged particle jet transverse + momentum for $\sqrt{s}=$0.9, 2.76 and 7 TeV.' Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Dec 13 15:29:18 2012 (r4082) +++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Fri Dec 14 13:38:42 2012 (r4083) @@ -131,6 +131,7 @@ CMS_2012_I1102908.info \ CMS_2012_I1107658.info \ CMS_2012_I1184941.info \ + CMS_2012_PAS_FWD_11_003.info \ CMS_2012_PAS_QCD_11_010.info \ CMS_QCD_10_024.info \ D0_1996_S3214044.info \ Added: branches/2012-06-aidarivet/data/plotinfo/CMS_2012_PAS_FWD_11_003.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/plotinfo/CMS_2012_PAS_FWD_11_003.plot Fri Dec 14 13:38:42 2012 (r4083) @@ -0,0 +1,21 @@ +# BEGIN PLOT /CMS_2012_PAS_FWD_11_003/d0 +XLabel=Leading charged jet $p_\perp$ [GeV] +YLabel=$(\mathrm{d}E^\text{hard} / \mathrm{d}\eta) / (\mathrm{d} E^\text{incl} / \mathrm{d}\eta)$ +LegendXPos=0.05 +LogY=0 +# END PLOT + + +# BEGIN PLOT /CMS_2012_PAS_FWD_11_003/d01-x01-y01 +Title=Ratio of energy deposited in $-6.6 < \eta < -5.2$ for $\sqrt{s}=0.9$ TeV +# END PLOT + +# BEGIN PLOT /CMS_2012_PAS_FWD_11_003/d02-x01-y01 +Title=Ratio of energy deposited in $-6.6 < \eta < -5.2$ for $\sqrt{s}=2.76$ TeV +# END PLOT + +# BEGIN PLOT /CMS_2012_PAS_FWD_11_003/d03-x01-y01 +Title=Ratio of energy deposited in $-6.6 < \eta < -5.2$ for $\sqrt{s}=7$ TeV +LegendXPos=0.45 +LegendYPos=0.4 +# END PLOT Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Dec 13 15:29:18 2012 (r4082) +++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Fri Dec 14 13:38:42 2012 (r4083) @@ -128,6 +128,7 @@ CMS_2012_I1102908.plot \ CMS_2012_I1107658.plot \ CMS_2012_I1184941.plot \ + CMS_2012_PAS_FWD_11_003.plot \ CMS_2012_PAS_QCD_11_010.plot \ CMS_QCD_10_024.plot \ D0_1996_S3214044.plot \ Added: branches/2012-06-aidarivet/data/refdata/CMS_2012_PAS_FWD_11_003.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/refdata/CMS_2012_PAS_FWD_11_003.aida Fri Dec 14 13:38:42 2012 (r4083) @@ -0,0 +1,99 @@ +<?xml version="1.0" encoding="ISO-8859-1" ?> +<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd"> +<aida version="3.3"> + <implementation version="1.1" package="FreeHEP"/> + <dataPointSet name="d01-x01-y01" dimension="2" + path="/REF/CMS_2012_PAS_FWD_11_003" title="energy flow ratio to pt, 0.9 TeV"> + <dataPoint> + <measurement errorPlus="0.5" value="1.5" errorMinus="0.5"/> + <measurement errorPlus="0.013615" value="1.03319" errorMinus="0.013615"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="0.5" value="2.5" errorMinus="0.5"/> + <measurement errorPlus="0.0143562" value="0.934369" errorMinus="0.0143562"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.0" value="4.0" errorMinus="1.0"/> + <measurement errorPlus="0.0169606" value="0.852281" errorMinus="0.0169606"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.25" value="6.25" errorMinus="1.25"/> + <measurement errorPlus="0.0255735" value="0.788704" errorMinus="0.0255735"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.25" value="8.75" errorMinus="1.25"/> + <measurement errorPlus="0.0305182" value="0.783938" errorMinus="0.0305182"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.5" value="12.5" errorMinus="2.5"/> + <measurement errorPlus="0.0306537" value="0.783507" errorMinus="0.0306537"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.0" value="20.0" errorMinus="5.0"/> + <measurement errorPlus="0.0389878" value="0.743184" errorMinus="0.0389878"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" + path="/REF/CMS_2012_PAS_FWD_11_003" title="energy flow ratio to pt, 2.76 TeV"> + <dataPoint> + <measurement errorPlus="0.5" value="1.5" errorMinus="0.5"/> + <measurement errorPlus="0.0142233" value="0.9948763" errorMinus="0.0142233"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="0.5" value="2.5" errorMinus="0.5"/> + <measurement errorPlus="0.0166567" value="1.05403" errorMinus="0.0166567"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.0" value="4.0" errorMinus="1.0"/> + <measurement errorPlus="0.0330071" value="1.06285" errorMinus="0.0330071"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.25" value="6.25" errorMinus="1.25"/> + <measurement errorPlus="0.0418842" value="1.05258" errorMinus="0.0418842"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.25" value="8.75" errorMinus="1.25"/> + <measurement errorPlus="0.0448775" value="1.04008" errorMinus="0.0448775"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.5" value="12.5" errorMinus="2.5"/> + <measurement errorPlus="0.0425066" value="1.01836" errorMinus="0.0425066"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.0" value="20.0" errorMinus="5.0"/> + <measurement errorPlus="0.0432732" value="0.985929" errorMinus="0.0432732"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" + path="/REF/CMS_2012_PAS_FWD_11_003" title="energy flow ratio to pt, 7 TeV"> + <dataPoint> + <measurement errorPlus="0.5" value="1.5" errorMinus="0.5"/> + <measurement errorPlus="0.0297284" value="0.8715" errorMinus="0.0297284"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="0.5" value="2.5" errorMinus="0.5"/> + <measurement errorPlus="0.0184714" value="1.04212" errorMinus="0.0184714"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.0" value="4.0" errorMinus="1.0"/> + <measurement errorPlus="0.0222113" value="1.19448" errorMinus="0.0222113"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.25" value="6.25" errorMinus="1.25"/> + <measurement errorPlus="0.0370915" value="1.30589" errorMinus="0.0370915"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.25" value="8.75" errorMinus="1.25"/> + <measurement errorPlus="0.0461365" value="1.34462" errorMinus="0.0461365"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="2.5" value="12.5" errorMinus="2.5"/> + <measurement errorPlus="0.0484821" value="1.34462" errorMinus="0.0484821"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="5.0" value="20.0" errorMinus="5.0"/> + <measurement errorPlus="0.0404986" value="1.30869" errorMinus="0.0404986"/> + </dataPoint> + </dataPointSet> +</aida> + Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Dec 13 15:29:18 2012 (r4082) +++ branches/2012-06-aidarivet/data/refdata/Makefile.am Fri Dec 14 13:38:42 2012 (r4083) @@ -85,6 +85,7 @@ CMS_2012_I1102908.aida \ CMS_2012_I1107658.aida \ CMS_2012_I1184941.aida \ + CMS_2012_PAS_FWD_11_003.aida \ CMS_2012_PAS_QCD_11_010.aida \ CMS_QCD_10_024.aida \ DELPHI_1994_S3021912.aida \ Added: branches/2012-06-aidarivet/src/Analyses/CMS_2012_PAS_FWD_11_003.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/src/Analyses/CMS_2012_PAS_FWD_11_003.cc Fri Dec 14 13:38:42 2012 (r4083) @@ -0,0 +1,172 @@ +// Samantha Dooling DESY +// February 2012 +// +// -*- C++ -*- +// ============================= +// +// Ratio of the energy deposited in the pseudorapditiy range +// -6.6 < eta < -5.2 for events with a charged particle jet +// +// ============================= +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/ChargedFinalState.hh" +#include "Rivet/Projections/FastJets.hh" +#include "Rivet/Projections/Beam.hh" +#include "Rivet/Projections/VetoedFinalState.hh" +#include "LWH/Histogram1D.h" + +namespace Rivet { + + + class CMS_2012_PAS_FWD_11_003 : public Analysis { + public: + + /// Constructor + CMS_2012_PAS_FWD_11_003() + : Analysis("CMS_2012_PAS_FWD_11_003") + { } + + void init() { + + // gives the range of eta and min pT for the final state from which I get the jets + FastJets jetpro (ChargedFinalState(-2.5, 2.5, 0.3*GeV), FastJets::ANTIKT, 0.5); + addProjection(jetpro, "Jets"); + + // skip Neutrinos and Muons + VetoedFinalState fsv(FinalState(-7.0, -4.0, 0.*GeV)); + fsv.vetoNeutrinos(); + fsv.addVetoPairId(MUON); + addProjection(fsv, "fsv"); + + // for the hadron level selection + VetoedFinalState sfsv(FinalState(-MAXRAPIDITY, MAXRAPIDITY, 0.*GeV)); + sfsv.vetoNeutrinos(); + sfsv.addVetoPairId(MUON); + addProjection(sfsv, "sfsv"); + + //counters + passedSumOfWeights = 0.; + inclEflow = 0.; + + // Temporary histograms to fill the energy flow for leading jet events. + // Ratios are calculated in finalyze(). + int id = 0; + if (fuzzyEquals(sqrtS()/GeV, 900, 1e-3)) id=1; + if (fuzzyEquals(sqrtS()/GeV, 2760, 1e-3)) id=2; + if (fuzzyEquals(sqrtS()/GeV, 7000, 1e-3)) id=3; + _tmp_jet = bookHistogram1D("eflow_jet", binEdges(id, 1, 1)); // Leading jet energy flow in pt + _tmp_njet = bookHistogram1D("number_jet", binEdges(id, 1, 1)); // Number of events in pt + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + + const double weight = event.weight(); + + // Skip if the event is empty + const FinalState& fsv = applyProjection<FinalState>(event, "fsv"); + if (fsv.empty()) vetoEvent; + + // ====================== Minimum Bias selection + + const FinalState& sfsv = applyProjection<FinalState>(event, "sfsv"); + ParticleVector parts = sfsv.particlesByRapidity(); + if (parts.empty()) vetoEvent; + + // find dymax + double dymax = 0; + int gap_pos = -1; + for (size_t i=0; i < parts.size()-1; ++i) { + double dy = parts[i+1].momentum().rapidity() - parts[i].momentum().rapidity(); + if (dy > dymax) { + dymax = dy; + gap_pos = i; + } + } + + // calculate mx2 and my2 + FourMomentum xmom; + for (int i=0; i<=gap_pos; ++i) { + xmom += parts[i].momentum(); + } + double mx2 = xmom.mass2(); + if (mx2<0) vetoEvent; + + FourMomentum ymom; + for (size_t i=gap_pos+1; i<parts.size(); ++i) { + ymom += parts[i].momentum(); + } + double my2 = ymom.mass2(); + if (my2<0) vetoEvent; + + // calculate xix and xiy and xidd + double xix = mx2 / sqr(sqrtS()); + double xiy = my2 / sqr(sqrtS()); + double xidd = mx2*my2 / sqr(sqrtS()*0.938*GeV); + + // combine the selection: xi cuts + bool passedHadronCuts = false; + if (fuzzyEquals(sqrtS()/GeV, 900, 1e-3) && (xix > 0.1 || xiy > 0.4 || xidd > 0.5)) passedHadronCuts = true; + if (fuzzyEquals(sqrtS()/GeV, 2760, 1e-3) && (xix > 0.07 || xiy > 0.2 || xidd > 0.5)) passedHadronCuts = true; + if (fuzzyEquals(sqrtS()/GeV, 7000, 1e-3) && (xix > 0.04 || xiy > 0.1 || xidd > 0.5)) passedHadronCuts = true; + if (!passedHadronCuts) vetoEvent; + + // ============================== MINIMUM BIAS EVENTS + + // loop over particles to calculate the energy + passedSumOfWeights += weight; + + foreach (const Particle& p, fsv.particles()) { + if (-5.2 > p.momentum().eta() && p.momentum().eta() > -6.6) inclEflow += weight*p.momentum().E()/GeV; + } + + // ============================== JET EVENTS + + const FastJets& jetpro = applyProjection<FastJets>(event, "Jets"); + const Jets& jets = jetpro.jetsByPt(1.0*GeV); + if (jets.size()<1) vetoEvent; + + if (fabs(jets[0].momentum().eta()) < 2.0) { + _tmp_njet->fill(jets[0].momentum().pT()/GeV, weight); + + // energy flow + foreach (const Particle& p, fsv.particles()) { + if (p.momentum().eta() > -6.6 && p.momentum().eta() < -5.2) { // ask for the CASTOR region + _tmp_jet->fill(jets[0].momentum().pT()/GeV, weight * p.momentum().E()/GeV); + } + } + } + + }// analysis + + void finalize() { + _tmp_jet->scale(passedSumOfWeights/inclEflow); + + AIDA::IHistogramFactory& hf = histogramFactory(); + if (fuzzyEquals(sqrtS()/GeV, 900, 1e-3)) hf.divide(histoDir() + "/d01-x01-y01", *_tmp_jet, *_tmp_njet); + if (fuzzyEquals(sqrtS()/GeV, 2760, 1e-3)) hf.divide(histoDir() + "/d02-x01-y01", *_tmp_jet, *_tmp_njet); + if (fuzzyEquals(sqrtS()/GeV, 7000, 1e-3)) hf.divide(histoDir() + "/d03-x01-y01", *_tmp_jet, *_tmp_njet); + + hf.destroy(_tmp_jet); + hf.destroy(_tmp_njet); + } + + private: + // counters + double passedSumOfWeights; + double inclEflow; + + // histograms + AIDA::IHistogram1D* _tmp_jet; + AIDA::IHistogram1D* _tmp_njet; + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(CMS_2012_PAS_FWD_11_003); + +} Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Dec 13 15:29:18 2012 (r4082) +++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Fri Dec 14 13:38:42 2012 (r4083) @@ -133,6 +133,7 @@ if ENABLE_PRELIMINARY RivetCMSAnalyses_la_SOURCES += \ CMS_QCD_10_024.cc \ + CMS_2012_PAS_FWD_11_003.cc \ CMS_2012_PAS_QCD_11_010.cc endif
More information about the Rivet-svn mailing list |