BasicHepMCHeavyIonValidation.h

Go to the documentation of this file.

#ifndef BASICHEPMCHEAVYIONVALIDATION_H
#define BASICHEPMCHEAVYIONVALIDATION_H

/*class BasicHepMCHeavyIonValidation
 *  
 *  Class to fill Event Generator dqm monitor elements; works on HepMCProduct HepMC::HeavyIon
 *  Quan Wang - 04/2013
 *
 */

// framework & common header files
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Run.h"

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"

//DQM services
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/MonitorElement.h"

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"

#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"

#include "Validation/EventGenerator/interface/WeightManager.h"

class BasicHepMCHeavyIonValidation : public edm::EDAnalyzer
{
    public:
        explicit BasicHepMCHeavyIonValidation(const edm::ParameterSet&);
        virtual ~BasicHepMCHeavyIonValidation();
        virtual void beginJob();
        virtual void endJob();  
        virtual void analyze(const edm::Event&, const edm::EventSetup&);
        virtual void beginRun(const edm::Run&, const edm::EventSetup&);
        virtual void endRun(const edm::Run&, const edm::EventSetup&);

    private:
        WeightManager wmanager_;
        edm::InputTag hepmcCollection_;
        bool QWdebug_;

        //edm::ESHandle<HepPDT::ParticleDataTable> fPDGTable ;

        DQMStore *dbe;

        MonitorElement* nEvt;

        // Additional information stored in HeavyIon structure
        MonitorElement* Ncoll_hard;     // Number of hard scatterings
        MonitorElement* Npart_proj;     // Number of projectile participants
        MonitorElement* Npart_targ;     // Number of target participants
        MonitorElement* Ncoll;          // Number of NN (nucleon-nucleon) collisions
        MonitorElement* N_Nwounded_collisions;      // Number of N-Nwounded collisions
        MonitorElement* Nwounded_N_collisions;      // Number of Nwounded-N collisons
        MonitorElement* Nwounded_Nwounded_collisions;   // Number of Nwounded-Nwounded collisions
        MonitorElement* spectator_neutrons;     // Number of spectator neutrons
        MonitorElement* spectator_protons;      // Number of spectator protons
        MonitorElement* impact_parameter;       // Impact Parameter(fm) of collision
        MonitorElement* event_plane_angle;      // Azimuthal angle of event plane
        MonitorElement* eccentricity;       // eccentricity of participating nucleons
                            // in the transverse plane 
                            // (as in phobos nucl-ex/0510031) 
        MonitorElement* sigma_inel_NN;      // nucleon-nucleon inelastic 
                            // (including diffractive) cross-section

        edm::EDGetTokenT<edm::HepMCProduct> hepmcCollectionToken_;
};

#endif