EventGenerator/interface/BasicHepMCHeavyIonValidation.h

0001 #ifndef BASICHEPMCHEAVYIONVALIDATION_H
0002 #define BASICHEPMCHEAVYIONVALIDATION_H
0003
0004 /*class BasicHepMCHeavyIonValidation
0005  *
0006  *  Class to fill Event Generator dqm monitor elements; works on HepMCProduct HepMC::HeavyIon
0007  *  Quan Wang - 04/2013
0008  *
0009  */
0010
0011 // framework & common header files
0012 #include "FWCore/Framework/interface/Event.h"
0013 #include "FWCore/Framework/interface/EventSetup.h"
0014 #include "FWCore/Framework/interface/Run.h"
0015
0016 #include "DataFormats/Common/interface/Handle.h"
0017 #include "FWCore/Framework/interface/ESHandle.h"
0018 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0019 #include "FWCore/Utilities/interface/InputTag.h"
0020
0021 //DQM services
0022 #include "DQMServices/Core/interface/DQMStore.h"
0023 #include "FWCore/ServiceRegistry/interface/Service.h"
0024 #include "DQMServices/Core/interface/DQMEDAnalyzer.h"
0025
0026 #include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
0027
0028 #include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"
0029
0030 #include "Validation/EventGenerator/interface/WeightManager.h"
0031
0032 class BasicHepMCHeavyIonValidation : public DQMEDAnalyzer {
0033 public:
0034   explicit BasicHepMCHeavyIonValidation(const edm::ParameterSet&);
0035   ~BasicHepMCHeavyIonValidation() override;
0036
0037   void bookHistograms(DQMStore::IBooker& i, edm::Run const&, edm::EventSetup const&) override;
0038   void analyze(edm::Event const&, edm::EventSetup const&) override;
0039
0040 private:
0041   WeightManager wmanager_;
0042   edm::InputTag hepmcCollection_;
0043   bool QWdebug_;
0044
0045   /// PDT table
0046   //edm::ESHandle<HepPDT::ParticleDataTable> fPDGTable ;
0047
0048   MonitorElement* nEvt;
0049
0050   // Additional information stored in HeavyIon structure
0051   MonitorElement* Ncoll_hard;                    // Number of hard scatterings
0052   MonitorElement* Npart_proj;                    // Number of projectile participants
0053   MonitorElement* Npart_targ;                    // Number of target participants
0054   MonitorElement* Ncoll;                         // Number of NN (nucleon-nucleon) collisions
0055   MonitorElement* N_Nwounded_collisions;         // Number of N-Nwounded collisions
0056   MonitorElement* Nwounded_N_collisions;         // Number of Nwounded-N collisons
0057   MonitorElement* Nwounded_Nwounded_collisions;  // Number of Nwounded-Nwounded collisions
0058   MonitorElement* spectator_neutrons;            // Number of spectator neutrons
0059   MonitorElement* spectator_protons;             // Number of spectator protons
0060   MonitorElement* impact_parameter;              // Impact Parameter(fm) of collision
0061   MonitorElement* event_plane_angle;             // Azimuthal angle of event plane
0062   MonitorElement* eccentricity;                  // eccentricity of participating nucleons
0063                                                  // in the transverse plane
0064                                                  // (as in phobos nucl-ex/0510031)
0065   MonitorElement* sigma_inel_NN;                 // nucleon-nucleon inelastic
0066                                                  // (including diffractive) cross-section
0067
0068   edm::EDGetTokenT<edm::HepMCProduct> hepmcCollectionToken_;
0069 };
0070
0071 #endif