|
[Rivet-svn] r4064 - in branches/2012-06-aidarivet: . data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgFri Dec 7 09:13:54 GMT 2012
Author: hoeth Date: Fri Dec 7 09:13:54 2012 New Revision: 4064 Log: Added ALICE_2012_I1181770 analysis Added: branches/2012-06-aidarivet/data/anainfo/ALICE_2012_I1181770.info branches/2012-06-aidarivet/data/plotinfo/ALICE_2012_I1181770.plot branches/2012-06-aidarivet/data/refdata/ALICE_2012_I1181770.aida branches/2012-06-aidarivet/src/Analyses/ALICE_2012_I1181770.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 6 21:38:55 2012 (r4063) +++ branches/2012-06-aidarivet/ChangeLog Fri Dec 7 09:13:54 2012 (r4064) @@ -1,3 +1,7 @@ +2012-12-07 Hendrik Hoeth <hendrik.hoeth at cern.ch> + + * Added ALICE_2012_I1181770 analysis + 2012-12-06 Hendrik Hoeth <hendrik.hoeth at cern.ch> * Added ATLAS_2012_I1188891 analysis Added: branches/2012-06-aidarivet/data/anainfo/ALICE_2012_I1181770.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/anainfo/ALICE_2012_I1181770.info Fri Dec 7 09:13:54 2012 (r4064) @@ -0,0 +1,39 @@ +Name: ALICE_2012_I1181770 +Year: 2012 +Summary: Measurement of inelastic, single- and double-diffraction cross sections in proton--proton collisions at the LHC with ALICE +Experiment: ALICE +Collider: LHC +SpireID: 9732586 +InspireID: 1181770 +Status: VALIDATED +Authors: + - Martin Poghosyan <Martin.Poghosyan at cern.ch> + - Sercan Sen <Sercan.Sen at cern.ch> + - Burak Bilki <bbilki at gmail.com> +References: + - arXiv:1208.4968 [hep-ex] +RunInfo: + Inelastic events (non-diffractive and inelastic diffractive). +NumEvents: 100K +Beams: [p+, p+] +Energies: [900, 2760, 7000] +Description: + 'Measurements of cross sections of inelastic and diffractive processes in proton-proton collisions at $\sqrt{s}$=900, 2760 and 7000 GeV. + The fractions of diffractive processes in inelastic collisions were determined from a study of gaps in charged particle pseudorapidity distributions. + Single-diffractive events are selected with $M_{X} < 200$ GeV $/c^{2}$ and double-diffractive events defined as NSD events with $\Delta\eta >3$. + To measure the inelastic cross section, beam properties were determined with van der Meer scans using a simulation of diffraction adjusted to data.' +BibKey: :2012sja +BibTeX: '@article{:2012sja, + author = "Abelev, Betty and others", + title = "{Measurement of inelastic, single- and double-diffraction + cross sections in proton--proton collisions at the LHC + with ALICE}", + collaboration = "ALICE Collaboration", + journal = "Eur. Phys. J. C", + year = "2012", + eprint = "1208.4968", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "CERN-PH-EP-2012-138", + SLACcitation = "%%CITATION = ARXIV:1208.4968;%%", +}' Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Dec 6 21:38:55 2012 (r4063) +++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Fri Dec 7 09:13:54 2012 (r4064) @@ -11,6 +11,7 @@ ALICE_2010_S8706239.info \ ALICE_2011_S8909580.info \ ALICE_2011_S8945144.info \ + ALICE_2012_I1181770.info \ ARGUS_1993_S2653028.info \ ARGUS_1993_S2669951.info \ ARGUS_1993_S2789213.info \ Added: branches/2012-06-aidarivet/data/plotinfo/ALICE_2012_I1181770.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/plotinfo/ALICE_2012_I1181770.plot Fri Dec 7 09:13:54 2012 (r4064) @@ -0,0 +1,28 @@ +# BEGIN PLOT /ALICE_2012_I1181770/d +XLabel=$\sqrt{s}$ [GeV] +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1181770/d01-x01-y0 +Title=Production ratios of SD with $M_{X} < 200$ GeV$/c^{2}$ to INEL +YLabel=$\sigma_\text{SD} / \sigma_\text{inel}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1181770/d02-x01-y0 +Title=Production ratios of DD with $\Delta\eta >3$ to INEL +YLabel=$\sigma_\text{DD} / \sigma_\text{inel}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1181770/d03-x01-y0 +Title=Single diffraction cross-section for $M_{X} < 200$ GeV$/c^{2}$ +YLabel=$\sigma_\text{SD}$ [mb] +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1181770/d04-x01-y0 +Title=Double diffraction cross-section for $\Delta\eta >3$ +YLabel=$\sigma_\text{DD}$ [mb] +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1181770/d05-x01-y0 +Title=Inleastic cross-section +YLabel=$\sigma_\text{inel}$ [mb] +# END PLOT Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Dec 6 21:38:55 2012 (r4063) +++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Fri Dec 7 09:13:54 2012 (r4064) @@ -11,6 +11,7 @@ ALICE_2010_S8706239.plot \ ALICE_2011_S8909580.plot \ ALICE_2011_S8945144.plot \ + ALICE_2012_I1181770.plot \ ARGUS_1993_S2653028.plot \ ARGUS_1993_S2669951.plot \ ARGUS_1993_S2789213.plot \ Added: branches/2012-06-aidarivet/data/refdata/ALICE_2012_I1181770.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/refdata/ALICE_2012_I1181770.aida Fri Dec 7 09:13:54 2012 (r4064) @@ -0,0 +1,140 @@ +<?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.0" package="LWH"/> + <dataPointSet name="d01-x01-y01" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="900.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.21" errorPlus="0.03" errorMinus="0.03"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y02" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="2760.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.20" errorPlus="0.07" errorMinus="0.08"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y03" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="7000.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.20" errorPlus="0.04" errorMinus="0.07"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="900.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.11" errorPlus="0.03" errorMinus="0.03"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y02" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="2760.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.12" errorPlus="0.05" errorMinus="0.05"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y03" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="7000.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.12" errorPlus="0.05" errorMinus="0.04"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="900.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="11.2" errorPlus="1.6" errorMinus="2.1"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y02" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="2760.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="12.2" errorPlus="3.9" errorMinus="5.3"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y03" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="7000.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="14.9" errorPlus="3.4" errorMinus="5.9"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y01" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="900.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="5.6" errorPlus="2.0" errorMinus="2.0"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y02" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="2760.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="7.8" errorPlus="3.2" errorMinus="3.2"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y03" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="7000.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="9.0" errorPlus="2.6" errorMinus="2.6"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y01" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="900.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="52.5" errorPlus="2.0" errorMinus="3.3"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y02" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="2760.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="62.8" errorPlus="2.6832815729997477" errorMinus="4.17612260356422"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y03" + title="" path="/REF/ALICE_2012_I1181770" dimension="2"> + <dimension dim="0" title="" /> + <dimension dim="1" title="" /> + <dataPoint> + <measurement value="7000.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="73.2" errorPlus="3.2802438933713454" errorMinus="5.2839379254491625"/> + </dataPoint> + </dataPointSet> +</aida> Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/refdata/Makefile.am Thu Dec 6 21:38:55 2012 (r4063) +++ branches/2012-06-aidarivet/data/refdata/Makefile.am Fri Dec 7 09:13:54 2012 (r4064) @@ -19,6 +19,7 @@ ALICE_2010_S8706239.aida \ ALICE_2011_S8909580.aida \ ALICE_2011_S8945144.aida \ + ALICE_2012_I1181770.aida \ ARGUS_1993_S2653028.aida \ ARGUS_1993_S2669951.aida \ ARGUS_1993_S2789213.aida \ Added: branches/2012-06-aidarivet/src/Analyses/ALICE_2012_I1181770.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/src/Analyses/ALICE_2012_I1181770.cc Fri Dec 7 09:13:54 2012 (r4064) @@ -0,0 +1,182 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/ChargedFinalState.hh" + +namespace Rivet { + + class ALICE_2012_I1181770 : public Analysis { + public: + + ALICE_2012_I1181770() + : Analysis("ALICE_2012_I1181770") + { } + + public: + + void init() { + addProjection(ChargedFinalState(),"CFS"); + + if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) { + _h_xsec_sd = bookHistogram1D(3, 1, 1); + _h_xsec_dd = bookHistogram1D(4, 1, 1); + _h_xsec_inel = bookHistogram1D(5, 1, 1); + } else if (fuzzyEquals(sqrtS()/GeV, 2760, 1E-3)) { + _h_xsec_sd = bookHistogram1D(3, 1, 2); + _h_xsec_dd = bookHistogram1D(4, 1, 2); + _h_xsec_inel = bookHistogram1D(5, 1, 2); + } else if (fuzzyEquals(sqrtS()/GeV, 7000, 1E-3)) { + _h_xsec_sd = bookHistogram1D(3, 1, 3); + _h_xsec_dd = bookHistogram1D(4, 1, 3); + _h_xsec_inel = bookHistogram1D(5, 1, 3); + } + } + + void analyze(const Event& event) { + const double weight = event.weight(); + + const ChargedFinalState& cfs = applyProjection<ChargedFinalState>(event, "CFS"); + + // fill INEL plots for each event + if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) { + _h_xsec_inel->fill( 900/GeV, weight ); + } else if (fuzzyEquals(sqrtS()/GeV, 2760, 1E-3)) { + _h_xsec_inel->fill( 2760/GeV, weight ); + } else if (fuzzyEquals(sqrtS()/GeV, 7000, 1E-3)) { + _h_xsec_inel->fill( 7000/GeV, weight ); + } + + double Eslowest(0.0), pslowest(0.0), yslowest(999.); + double Efastest(0.0), pfastest(0.0), yfastest(-999.); + int pidslowest(0), pidfastest(0); + + double LRG(0.0), etapre(0.0), gapbwd(0.0), gapfwd(0.0); + unsigned int num_p(0); + + FourMomentum leadP(0.,0.,0.,0.); + double Elead(0.0), plead(0.0); + + foreach(const Particle& p, cfs.particlesByEta()) { //sorted from minus to plus + const PdgId pid = p.pdgId(); + double y = p.momentum().rapidity(); + double eta = p.momentum().eta(); + //SD case + if (y < yslowest) { + Eslowest = p.momentum().E(); + pslowest = p.momentum().vector3().mod(); + yslowest = p.momentum().rapidity(); + pidslowest = pid; + } + if (y > yfastest) { + Efastest = p.momentum().E(); + pfastest = p.momentum().vector3().mod(); + yfastest = p.momentum().rapidity(); + pidfastest = pid; + } + + num_p += 1; + // DD case + if (num_p==1) { + etapre = p.momentum().eta(); + } else if (num_p > 1) { + if (num_p==2) gapbwd = fabs(eta-etapre); + + double gap = fabs(eta-etapre); + LRG = (gap > LRG ? gap : LRG); // largest gap + + if (num_p==cfs.size()) gapfwd = fabs(eta-etapre); + etapre = eta; + } + } + + // Mx calculation + if (pidslowest==2212 && pidfastest==2212) { + if (fabs(yslowest) > fabs(yfastest)) { + Elead = Eslowest; + plead = pslowest; + } else if (fabs(yslowest) < fabs(yfastest)) { + Elead = Efastest; + plead = pfastest; + } else { + Elead = Eslowest; + plead = pslowest; + if ( (double)rand() / (double)RAND_MAX > 0.5) { // generate random number in [0.,1.] range and make decision randomly + Elead = Efastest; + plead = pfastest; + } + } + } else if (pidslowest==2212) { + Elead = Eslowest; + plead = pslowest; + } else if (pidfastest==2212) { + Elead = Efastest; + plead = pfastest; + } + + double Mx = sqrt((sqrtS()/GeV-Elead-plead)*(sqrtS()/GeV-Elead+plead)); + bool singleDiff = false; + + // Fill SD + if (Mx < 200.) { + singleDiff = true; + if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) { + _h_xsec_sd->fill( 900/GeV, weight); + } else if (fuzzyEquals(sqrtS()/GeV, 2760, 1E-3)) { + _h_xsec_sd->fill( 2760/GeV, weight); + } else if (fuzzyEquals(sqrtS()/GeV, 7000, 1E-3)) { + _h_xsec_sd->fill( 7000/GeV, weight); + } + } + + if ( singleDiff ) vetoEvent; // DD events are defined as NSD with large gap. + + // also remove SD-like events in NSD events + if ( std::abs(gapbwd-LRG) < std::numeric_limits<double>::epsilon() || std::abs(gapfwd-LRG) < std::numeric_limits<double>::epsilon() ) vetoEvent; + + // Fill DD plots + if (LRG > 3.) { + if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) { + _h_xsec_dd->fill( 900/GeV, weight); + } else if (fuzzyEquals(sqrtS()/GeV, 2760, 1E-3)) { + _h_xsec_dd->fill( 2760/GeV, weight); + } else if (fuzzyEquals(sqrtS()/GeV, 7000, 1E-3)) { + _h_xsec_dd->fill( 7000/GeV, weight); + } + } + } + + void finalize() { + + // get the ratio plots: SD/inel, DD/inel + const string dir = histoDir(); + + if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) { + histogramFactory().divide( dir + "/d01-x01-y01", *_h_xsec_sd , *_h_xsec_inel ); + histogramFactory().divide( dir + "/d02-x01-y01", *_h_xsec_dd , *_h_xsec_inel ); + } else if (fuzzyEquals(sqrtS()/GeV, 2760, 1E-3)) { + histogramFactory().divide( dir + "/d01-x01-y02", *_h_xsec_sd , *_h_xsec_inel ); + histogramFactory().divide( dir + "/d02-x01-y02", *_h_xsec_dd , *_h_xsec_inel ); + } else if (fuzzyEquals(sqrtS()/GeV, 7000, 1E-3)) { + histogramFactory().divide( dir + "/d01-x01-y03", *_h_xsec_sd , *_h_xsec_inel ); + histogramFactory().divide( dir + "/d02-x01-y03", *_h_xsec_dd , *_h_xsec_inel ); + } + + scale(_h_xsec_sd, crossSection()/millibarn/sumOfWeights()); + scale(_h_xsec_dd, crossSection()/millibarn/sumOfWeights()); + scale(_h_xsec_inel, crossSection()/millibarn/sumOfWeights()); + + } + + private: + + AIDA::IHistogram1D *_h_xsec_sd; + AIDA::IHistogram1D *_h_xsec_dd; + AIDA::IHistogram1D *_h_xsec_inel; + + }; + + // Hook for the plugin system + DECLARE_RIVET_PLUGIN(ALICE_2012_I1181770); + +} Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Dec 6 21:38:55 2012 (r4063) +++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Fri Dec 7 09:13:54 2012 (r4064) @@ -40,7 +40,8 @@ ALICE_2010_S8625980.cc \ ALICE_2010_S8706239.cc \ ALICE_2011_S8909580.cc \ - ALICE_2011_S8945144.cc + ALICE_2011_S8945144.cc \ + ALICE_2012_I1181770.cc lib_LTLIBRARIES += RivetATLASAnalyses.la
More information about the Rivet-svn mailing list |