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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

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_

Detailed Description

Definition at line 34 of file BasicHepMCHeavyIonValidation.h.

Constructor & Destructor Documentation

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 ( )

override

Definition at line 24 of file BasicHepMCHeavyIonValidation.cc.

{}

Member Function Documentation

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

override

counters to zero for every event

Gathering the HepMCProduct information

Definition at line 52 of file BasicHepMCHeavyIonValidation.cc.

References gather_cfg::cout, eccentricity, event_plane_angle, MonitorElement::Fill(), edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), hepmcCollectionToken_, 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

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

override

Setting the DQM top directories

Booking the ME's

Definition at line 26 of file BasicHepMCHeavyIonValidation.cc.

References DQMHelper::book1dHisto(), eccentricity, event_plane_angle, impact_parameter, N_Nwounded_collisions, Ncoll, Ncoll_hard, nEvt, Npart_proj, Npart_targ, Nwounded_N_collisions, Nwounded_Nwounded_collisions, pi, DQMStore::IBooker::setCurrentFolder(), 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

MonitorElement* BasicHepMCHeavyIonValidation::eccentricity

private

Definition at line 64 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::event_plane_angle

private

Definition at line 63 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag BasicHepMCHeavyIonValidation::hepmcCollection_

private

Definition at line 44 of file BasicHepMCHeavyIonValidation.h.

Referenced by BasicHepMCHeavyIonValidation().

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

private

Definition at line 70 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and BasicHepMCHeavyIonValidation().

MonitorElement* BasicHepMCHeavyIonValidation::impact_parameter

private

Definition at line 62 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::N_Nwounded_collisions

private

Definition at line 57 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::Ncoll

private

Definition at line 56 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::Ncoll_hard

private

Definition at line 53 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::nEvt

private

PDT table.

Definition at line 50 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::Npart_proj

private

Definition at line 54 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::Npart_targ

private

Definition at line 55 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::Nwounded_N_collisions

private

Definition at line 58 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::Nwounded_Nwounded_collisions

private

Definition at line 59 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

bool BasicHepMCHeavyIonValidation::QWdebug_

private

Definition at line 45 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and BasicHepMCHeavyIonValidation().

MonitorElement* BasicHepMCHeavyIonValidation::sigma_inel_NN

private

Definition at line 67 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::spectator_neutrons

private

Definition at line 60 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* BasicHepMCHeavyIonValidation::spectator_protons

private

Definition at line 61 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze(), and bookHistograms().

WeightManager BasicHepMCHeavyIonValidation::wmanager_

private

Definition at line 43 of file BasicHepMCHeavyIonValidation.h.

Referenced by analyze().