BasicHepMCHeavyIonValidation Class Reference

#include <BasicHepMCHeavyIonValidation.h>

Inheritance diagram for BasicHepMCHeavyIonValidation:

Public Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
	BasicHepMCHeavyIonValidation (const edm::ParameterSet &)
void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override
	~BasicHepMCHeavyIonValidation () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Attributes
MonitorElement *	eccentricity
MonitorElement *	event_plane_angle
edm::InputTag	hepmcCollection_
edm::EDGetTokenT< edm::HepMCProduct >	hepmcCollectionToken_
MonitorElement *	impact_parameter
MonitorElement *	N_Nwounded_collisions
MonitorElement *	Ncoll
MonitorElement *	Ncoll_hard
MonitorElement *	nEvt
	PDT table. More...
MonitorElement *	Npart_proj
MonitorElement *	Npart_targ
MonitorElement *	Nwounded_N_collisions
MonitorElement *	Nwounded_Nwounded_collisions
bool	QWdebug_
MonitorElement *	sigma_inel_NN
MonitorElement *	spectator_neutrons
MonitorElement *	spectator_protons
WeightManager	wmanager_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 32 of file BasicHepMCHeavyIonValidation.h.

Constructor & Destructor Documentation

BasicHepMCHeavyIonValidation()

BasicHepMCHeavyIonValidation::BasicHepMCHeavyIonValidation ( const edm::ParameterSet &  iPSet )

explicit

Definition at line 15 of file BasicHepMCHeavyIonValidation.cc.

References edm::ParameterSet::getUntrackedParameter(), hepmcCollection_, hepmcCollectionToken_, and QWdebug_.

    : wmanager_(iPSet, consumesCollector()), hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")) {
  QWdebug_ = iPSet.getUntrackedParameter<bool>("QWdebug", false);

  hepmcCollectionToken_ = consumes<HepMCProduct>(hepmcCollection_);
}

~BasicHepMCHeavyIonValidation()

BasicHepMCHeavyIonValidation::~BasicHepMCHeavyIonValidation ( )

override

Definition at line 22 of file BasicHepMCHeavyIonValidation.cc.

{}

Member Function Documentation

analyze()

void BasicHepMCHeavyIonValidation::analyze ( edm::Event const &  iEvent, edm::EventSetup const &  iSetup  )

overridevirtual

counters to zero for every event

Gathering the HepMCProduct information

Reimplemented from DQMEDAnalyzer.

Definition at line 81 of file BasicHepMCHeavyIonValidation.cc.

References gather_cfg::cout, eccentricity, event_plane_angle, dqm::impl::MonitorElement::Fill(), edm::HepMCProduct::GetEvent(), hepmcCollectionToken_, iEvent, impact_parameter, N_Nwounded_collisions, Ncoll, Ncoll_hard, nEvt, Npart_proj, Npart_targ, Nwounded_N_collisions, Nwounded_Nwounded_collisions, QWdebug_, sigma_inel_NN, spectator_neutrons, spectator_protons, WeightManager::weight(), mps_merge::weight, and wmanager_.

                                                                                              {

  edm::Handle<HepMCProduct> evt;
  iEvent.getByToken(hepmcCollectionToken_, evt);

  //Get EVENT
  //HepMC::GenEvent *myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));

  const HepMC::HeavyIon* ion = evt->GetEvent()->heavy_ion();

  if (!ion) {
    if (QWdebug_)
      std::cout << "!!QW!! HeavyIon == null" << std::endl;
    return;
  }

  double weight = wmanager_.weight(iEvent);
  nEvt->Fill(0.5, weight);

  Ncoll_hard->Fill(ion->Ncoll_hard(), weight);
  Npart_proj->Fill(ion->Npart_proj(), weight);
  Npart_targ->Fill(ion->Npart_targ(), weight);
  Ncoll->Fill(ion->Ncoll(), weight);
  N_Nwounded_collisions->Fill(ion->N_Nwounded_collisions(), weight);
  Nwounded_N_collisions->Fill(ion->Nwounded_N_collisions(), weight);
  Nwounded_Nwounded_collisions->Fill(ion->Nwounded_Nwounded_collisions(), weight);
  spectator_neutrons->Fill(ion->spectator_neutrons(), weight);
  spectator_protons->Fill(ion->spectator_protons(), weight);
  impact_parameter->Fill(ion->impact_parameter(), weight);
  event_plane_angle->Fill(ion->event_plane_angle(), weight);
  eccentricity->Fill(ion->eccentricity(), weight);
  sigma_inel_NN->Fill(ion->sigma_inel_NN(), weight);

  //delete myGenEvent;
}  //analyze

bookHistograms()

void BasicHepMCHeavyIonValidation::bookHistograms ( DQMStore::IBooker &  i, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Setting the DQM top directories

Booking the ME's

Implements DQMEDAnalyzer.

Definition at line 24 of file BasicHepMCHeavyIonValidation.cc.

References eccentricity, event_plane_angle, mps_fire::i, impact_parameter, N_Nwounded_collisions, Ncoll, Ncoll_hard, nEvt, Npart_proj, Npart_targ, Nwounded_N_collisions, Nwounded_Nwounded_collisions, pi, sigma_inel_NN, spectator_neutrons, and spectator_protons.

                                                                                                         {
  DQMHelper dqm(&i);
  i.setCurrentFolder("Generator/HeavyIon");

  // Number of analyzed events
  nEvt = dqm.book1dHisto("nEvt", "n analyzed Events", 1, 0., 1., "", "Number of Events");

  Ncoll_hard =
      dqm.book1dHisto("Ncoll_hard", "Ncoll_hard", 700, 0, 700, "Number of hard scatterings", "Number of Events");
  Npart_proj =
      dqm.book1dHisto("Npart_proj", "Npart_proj", 250, 0, 250, "Number of projectile participants", "Number of Events");
  Npart_targ =
      dqm.book1dHisto("Npart_targ", "Npart_targ", 250, 0, 250, "Number of target participants", "Number of Events");
  Ncoll = dqm.book1dHisto("Ncoll", "Ncoll", 700, 0, 700, "Number of N-N collisions", "Number of Events");
  N_Nwounded_collisions = dqm.book1dHisto("N_Nwounded_collisions",
                                          "N_Nwounded_collisions",
                                          250,
                                          0,
                                          250,
                                          "Number of N-N wounded collisions",
                                          "Number of Events");
  Nwounded_N_collisions = dqm.book1dHisto("Nwounded_N_collisions",
                                          "Nwounded_N_collisions",
                                          250,
                                          0,
                                          250,
                                          "Number of N wounded-N collisions",
                                          "Number of Events");
  Nwounded_Nwounded_collisions = dqm.book1dHisto("Nwounded_Nwounded_collisions",
                                                 "Nwounded_Nwounded_collisions",
                                                 250,
                                                 0,
                                                 250,
                                                 "Number of N wounded-N wounded collisions",
                                                 "Number of Events");
  spectator_neutrons = dqm.book1dHisto(
      "spectator_neutrons", "spectator_neutrons", 250, 0, 250, "Number of spectator neutrons", "Number of Events");
  spectator_protons = dqm.book1dHisto(
      "spectator_protons", "spectator_protons", 250, 0, 250, "Number of spectator protons", "Number of Events");
  impact_parameter = dqm.book1dHisto(
      "impact_parameter", "impact_parameter", 50, 0, 50, "Empact parameter of collision (fm)", "Number of Events");
  event_plane_angle = dqm.book1dHisto(
      "event_plane_angle", "event_plane_angle", 200, 0, 2 * CLHEP::pi, "#phi_{event plane} (rad)", "Number of Events");
  eccentricity = dqm.book1dHisto("eccentricity", "eccentricity", 200, 0, 1.0, "Eccentricity", "Number of Events");
  sigma_inel_NN = dqm.book1dHisto("sigma_inel_NN",
                                  "sigma_inel_NN",
                                  200,
                                  0,
                                  10.0,
                                  "#sigma{nucleon-nucleon inelastic cross-section}",
                                  "Number of Events");

  return;
}

Member Data Documentation

eccentricity

MonitorElement* BasicHepMCHeavyIonValidation::eccentricity

private

Definition at line 62 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

event_plane_angle

MonitorElement* BasicHepMCHeavyIonValidation::event_plane_angle

private

Definition at line 61 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

hepmcCollection_

edm::InputTag BasicHepMCHeavyIonValidation::hepmcCollection_

private

Definition at line 42 of file BasicHepMCHeavyIonValidation.h.

Referenced by BasicHepMCHeavyIonValidation().

hepmcCollectionToken_

edm::EDGetTokenT<edm::HepMCProduct> BasicHepMCHeavyIonValidation::hepmcCollectionToken_

private

Definition at line 68 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and BasicHepMCHeavyIonValidation().

impact_parameter

MonitorElement* BasicHepMCHeavyIonValidation::impact_parameter

private

Definition at line 60 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

N_Nwounded_collisions

MonitorElement* BasicHepMCHeavyIonValidation::N_Nwounded_collisions

private

Definition at line 55 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

Ncoll

MonitorElement* BasicHepMCHeavyIonValidation::Ncoll

private

Definition at line 54 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

Ncoll_hard

MonitorElement* BasicHepMCHeavyIonValidation::Ncoll_hard

private

Definition at line 51 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

nEvt

MonitorElement* BasicHepMCHeavyIonValidation::nEvt

private

PDT table.

Definition at line 48 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

Npart_proj

MonitorElement* BasicHepMCHeavyIonValidation::Npart_proj

private

Definition at line 52 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

Npart_targ

MonitorElement* BasicHepMCHeavyIonValidation::Npart_targ

private

Definition at line 53 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

Nwounded_N_collisions

MonitorElement* BasicHepMCHeavyIonValidation::Nwounded_N_collisions

private

Definition at line 56 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

Nwounded_Nwounded_collisions

MonitorElement* BasicHepMCHeavyIonValidation::Nwounded_Nwounded_collisions

private

Definition at line 57 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

QWdebug_

bool BasicHepMCHeavyIonValidation::QWdebug_

private

Definition at line 43 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and BasicHepMCHeavyIonValidation().

sigma_inel_NN

MonitorElement* BasicHepMCHeavyIonValidation::sigma_inel_NN

private

Definition at line 65 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

spectator_neutrons

MonitorElement* BasicHepMCHeavyIonValidation::spectator_neutrons

private

Definition at line 58 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

spectator_protons

MonitorElement* BasicHepMCHeavyIonValidation::spectator_protons

private

Definition at line 59 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

wmanager_

WeightManager BasicHepMCHeavyIonValidation::wmanager_

private

Definition at line 41 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze().