|
[Rivet-svn] r2917 - in trunk: data/anainfo src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgMon Feb 7 13:40:59 GMT 2011
Author: buckley Date: Mon Feb 7 13:40:58 2011 New Revision: 2917 Log: Finishing ALICE analysis metadata Modified: trunk/data/anainfo/ALICE_2010_S8706239.info trunk/src/Analyses/ALICE_2010_S8706239.cc Modified: trunk/data/anainfo/ALICE_2010_S8706239.info ============================================================================== --- trunk/data/anainfo/ALICE_2010_S8706239.info Sat Feb 5 23:01:06 2011 (r2916) +++ trunk/data/anainfo/ALICE_2010_S8706239.info Mon Feb 7 13:40:58 2011 (r2917) @@ -1,6 +1,6 @@ Name: ALICE_2010_S8706239 Year: 2010 -Summary: Charged +Summary: Charged particle $\langle p_\perp \rangle$ vs. $N_\text{ch}$ in pp collisions at 900 GeV Experiment: ALICE Collider: LHC SpiresID: 8706239 Modified: trunk/src/Analyses/ALICE_2010_S8706239.cc ============================================================================== --- trunk/src/Analyses/ALICE_2010_S8706239.cc Sat Feb 5 23:01:06 2011 (r2916) +++ trunk/src/Analyses/ALICE_2010_S8706239.cc Mon Feb 7 13:40:58 2011 (r2917) @@ -18,7 +18,6 @@ : Analysis("ALICE_2010_S8706239"), _Nevt_after_cuts(0.0) { - /// @todo Set whether your finalize method needs the generator cross section setNeedsCrossSection(false); setBeams(PROTON, PROTON); } @@ -36,9 +35,9 @@ ChargedFinalState cfs(-0.8, 0.8, 0.15); addProjection(cfs, "CFS"); - + _h_pT = bookHistogram1D(4, 1, 1); - + _h_pT_Nch_015 = bookProfile1D(11, 1, 1); _h_pT_Nch_05 = bookProfile1D(12, 1, 1); @@ -48,9 +47,9 @@ /// Perform the per-event analysis void analyze(const Event& event) { const double weight = event.weight(); - + const ChargedFinalState& charged = applyProjection<ChargedFinalState>(event, "CFS"); - + _Nevt_after_cuts += weight; // Get number of particles that fulfill certain pT requirements @@ -58,8 +57,8 @@ int Nch_05 = 0; foreach (const Particle& p, charged.particles()) { double pT = p.momentum().pT()/GeV; - if (pT < 4.0) Nch_015++; - if (pT > 0.5 && pT < 4.0) Nch_05++; + if (pT < 4.0) Nch_015++; + if (pT > 0.5 && pT < 4.0) Nch_05++; } // Now we can fill histograms @@ -77,20 +76,19 @@ /// Normalise histograms etc., after the run void finalize() { - scale(_h_pT, 1.0/_Nevt_after_cuts); - } //@} + private: /// @name Histograms //@{ AIDA::IHistogram1D* _h_pT; - + AIDA::IProfile1D* _h_pT_Nch_015 ; AIDA::IProfile1D* _h_pT_Nch_05 ;
More information about the Rivet-svn mailing list |