[Rivet-svn] rivet: ALICE_2012_I1127497 new centrality calculation method + f...

[Rivet Mercurial]

details:   https://rivet.hepforge.org/hg/rivet/rev/26ed1f040772
branches:  heavyion
changeset: 6012:26ed1f040772
user:      Przemyslaw Karczmarczyk <przemyslaw.karczmarczyk at cern.ch>
date:      Mon Sep 11 16:56:07 2017 +0200
description:
ALICE_2012_I1127497 new centrality calculation method + fix in ALICE_2010_I880049 analysis

diffs (truncated from 86 to 50 lines):

--- a/analyses/pluginALICE/ALICE_2010_I880049.cc	Mon Aug 21 15:50:57 2017 +0200
+++ b/analyses/pluginALICE/ALICE_2010_I880049.cc	Mon Sep 11 16:56:07 2017 +0200
@@ -22,7 +22,7 @@
       HeavyIonAnalysis::init();
       
       // Set centrality method
-      addCentralityMethod(HeavyIonAnalysis::ImpactParameter, 10000, "method1");
+      addCentralityMethod(HeavyIonAnalysis::ImpactParameter, 10000, "ImpactParameterMethod");
       
       // Charged final states with |eta| < 0.5 and pT > 50 MeV
       const Cut& cut = Cuts::abseta < 0.5 && Cuts::pT > 50*MeV;
--- a/analyses/pluginALICE/ALICE_2012_I1127497.cc	Mon Aug 21 15:50:57 2017 +0200
+++ b/analyses/pluginALICE/ALICE_2012_I1127497.cc	Mon Sep 11 16:56:07 2017 +0200
@@ -23,11 +23,16 @@
 
 
     void init() {
-
-      // charged final states
-      const ChargedFinalState cfs(-0.8, 0.8, 150*MeV);
+      HeavyIonAnalysis::init();
+      
+      // Set centrality method
+      addCentralityMethod(HeavyIonAnalysis::ImpactParameter, 10000, "method1");
+      
+      // Charged final states with |eta| < 0.5 and pT > 150 MeV
+      const Cut& cut = Cuts::abseta < 0.5 && Cuts::pT > 150*MeV;
+      const ChargedFinalState cfs(cut);
       addProjection(cfs,"CFS");
-
+      
       // @debug Member "log" is declared below.
       //log.open("/data/mcplots/mcplots/scripts/mcprod/testarea/logs/log.log");
 
@@ -71,18 +76,25 @@
       Particles chargedParticles = charged.particlesByPt();
 
       // @note Choose this convention here since JEWEL return a centrality value of -1 for a vacuum run.
-      float centr = -1;
+      //float centr = -1;
       
       // @note It has to be checked explicitly, if the pointer is sane and is no null_ptr. The standard
       // constructor of HepMC does not initialize the pointer and returns a null_ptr in case no heavy-ion
       // was written to the HepMC file. This is e.g. the case for the current pp generators. Trying to
       // access this naivly may therefore lead to a segmentation fault.
-      if( event.genEvent()->heavy_ion() ) {
+      //if( event.genEvent()->heavy_ion() ) {
         // JEWEL saves the centrality as the impact parameter. This should be changed later in JEWEL directly.
         // Note, that changes of HepMC may be required as well. However, given the current implementation,
         // the centrality can be accessed via