ME0HitsValidation Class Reference

#include <ME0HitsValidation.h>

Inheritance diagram for ME0HitsValidation:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	ME0HitsValidation (const edm::ParameterSet &)
	~ME0HitsValidation () override
Public Member Functions inherited from ME0BaseValidation
MonitorElement *	BookHistXY (DQMStore::IBooker &, const char *name, const char *label, unsigned int region_num, unsigned int layer_num=99)
MonitorElement *	BookHistZR (DQMStore::IBooker &, const char *name, const char *label, unsigned int region_num, unsigned int layer_num=99)
	ME0BaseValidation (const edm::ParameterSet &ps)
	~ME0BaseValidation () 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 *	me0_sh_eloss [2][6]
MonitorElement *	me0_sh_elossMu [2][6]
MonitorElement *	me0_sh_tof [2][6]
MonitorElement *	me0_sh_tofMu [2][6]
MonitorElement *	me0_sh_tot_zr [2]
MonitorElement *	me0_sh_xy [2][6]
MonitorElement *	me0_sh_zr [2][6]
Int_t	npart

Additional Inherited Members
Protected Attributes inherited from ME0BaseValidation
edm::EDGetToken	InputTagToken_
std::vector< std::string >	layerLabel
int	nBinXY_
std::vector< double >	nBinZR_
std::vector< double >	RangeZR_
std::vector< std::string >	regionLabel

Detailed Description

Definition at line 6 of file ME0HitsValidation.h.

Constructor & Destructor Documentation

ME0HitsValidation::ME0HitsValidation ( const edm::ParameterSet &  cfg )

explicit

Definition at line 5 of file ME0HitsValidation.cc.

References edm::ParameterSet::getParameter(), and ME0BaseValidation::InputTagToken_.

                                                               : ME0BaseValidation(cfg) {
  InputTagToken_ = consumes<edm::PSimHitContainer>(cfg.getParameter<edm::InputTag>("simInputLabel"));
}

ME0HitsValidation::~ME0HitsValidation ( )

override

Definition at line 58 of file ME0HitsValidation.cc.

{}

Member Function Documentation

void ME0HitsValidation::analyze ( const edm::Event &  e, const edm::EventSetup &  iSetup  )

overridevirtual

Implements ME0BaseValidation.

Definition at line 60 of file ME0HitsValidation.cc.

References funct::abs(), gather_cfg::cout, MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), hfClusterShapes_cfi::hits, triggerObjects_cff::id, ME0Geometry::idToDet(), ME0BaseValidation::InputTagToken_, createfilelist::int, edm::HandleBase::isValid(), me0_sh_eloss, me0_sh_elossMu, me0_sh_tof, me0_sh_tofMu, me0_sh_tot_zr, me0_sh_xy, me0_sh_zr, PV3DBase< T, PVType, FrameType >::perp(), spr::timeOfFlight(), and GeomDet::toGlobal().

                                                                              {
  edm::ESHandle<ME0Geometry> hGeom;
  iSetup.get<MuonGeometryRecord>().get(hGeom);
  const ME0Geometry *ME0Geometry_ = (&*hGeom);

  edm::Handle<edm::PSimHitContainer> ME0Hits;
  e.getByToken(InputTagToken_, ME0Hits);
  if (!ME0Hits.isValid()) {
    edm::LogError("ME0HitsValidation") << "Cannot get ME0Hits by Token simInputTagToken";
    return;
  }

  Float_t timeOfFlightMuon = 0.;
  Float_t energyLossMuon = 0;

  for (auto hits = ME0Hits->begin(); hits != ME0Hits->end(); hits++) {
    const ME0DetId id(hits->detUnitId());
    Int_t region = id.region();
    Int_t layer = id.layer();

    // Int_t even_odd = id.chamber()%2;
    if (ME0Geometry_->idToDet(hits->detUnitId()) == nullptr) {
      std::cout << "simHit did not matched with GEMGeometry." << std::endl;
      continue;
    }

    const LocalPoint hitLP(hits->localPosition());

    const GlobalPoint hitGP(ME0Geometry_->idToDet(hits->detUnitId())->surface().toGlobal(hitLP));
    Float_t g_r = hitGP.perp();
    Float_t g_x = hitGP.x();
    Float_t g_y = hitGP.y();
    Float_t g_z = hitGP.z();
    Float_t energyLoss = hits->energyLoss();
    Float_t timeOfFlight = hits->timeOfFlight();

    if (abs(hits->particleType()) == 13) {
      timeOfFlightMuon = hits->timeOfFlight();
      energyLossMuon = hits->energyLoss();
      // fill histos for Muons only
      me0_sh_tofMu[(int)(region / 2. + 0.5)][layer - 1]->Fill(timeOfFlightMuon);
      me0_sh_elossMu[(int)(region / 2. + 0.5)][layer - 1]->Fill(energyLossMuon * 1.e9);
    }

    me0_sh_zr[(int)(region / 2. + 0.5)][layer - 1]->Fill(g_z, g_r);
    me0_sh_tot_zr[(int)(region / 2. + 0.5)]->Fill(g_z, g_r);
    me0_sh_xy[(int)(region / 2. + 0.5)][layer - 1]->Fill(g_x, g_y);
    me0_sh_tof[(int)(region / 2. + 0.5)][layer - 1]->Fill(timeOfFlight);
    me0_sh_eloss[(int)(region / 2. + 0.5)][layer - 1]->Fill(energyLoss * 1.e9);
  }
}

void ME0HitsValidation::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  Run, edm::EventSetup const &  iSetup  )

override

Definition at line 9 of file ME0HitsValidation.cc.

References DQMStore::IBooker::book1D(), ME0BaseValidation::BookHistXY(), ME0BaseValidation::BookHistZR(), ME0BaseValidation::layerLabel, LogDebug, me0_sh_eloss, me0_sh_elossMu, me0_sh_tof, me0_sh_tofMu, me0_sh_tot_zr, me0_sh_xy, me0_sh_zr, ME0BaseValidation::regionLabel, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                               {
  LogDebug("MuonME0HitsValidation") << "Info : Loading Geometry information\n";
  ibooker.setCurrentFolder("MuonME0HitsV/ME0HitsTask");

  unsigned int nregion = 2;

  edm::LogInfo("MuonME0HitsValidation") << "+++ Info : # of region : " << nregion << std::endl;

  LogDebug("MuonME0HitsValidation") << "+++ Info : finish to get geometry information from ES.\n";

  for (unsigned int region_num = 0; region_num < nregion; region_num++) {
    me0_sh_tot_zr[region_num] = BookHistZR(ibooker, "me0_sh", "SimHit", region_num);
    for (unsigned int layer_num = 0; layer_num < 6; layer_num++) {
      me0_sh_zr[region_num][layer_num] = BookHistZR(ibooker, "me0_sh", "SimHit", region_num, layer_num);
      me0_sh_xy[region_num][layer_num] = BookHistXY(ibooker, "me0_sh", "SimHit", region_num, layer_num);
      std::string hist_name_for_tof =
          std::string("me0_sh_tof_r") + regionLabel[region_num] + "_l" + layerLabel[layer_num];
      std::string hist_name_for_tofMu =
          std::string("me0_sh_tofMuon_r") + regionLabel[region_num] + "_l" + layerLabel[layer_num];
      std::string hist_name_for_eloss =
          std::string("me0_sh_energyloss_r") + regionLabel[region_num] + "_l" + layerLabel[layer_num];
      std::string hist_name_for_elossMu =
          std::string("me0_sh_energylossMuon_r") + regionLabel[region_num] + "_l" + layerLabel[layer_num];
      std::string hist_label_for_xy = "SimHit occupancy : region" + regionLabel[region_num] + " layer " +
                                      layerLabel[layer_num] + " ; globalX [cm]; globalY[cm]";
      std::string hist_label_for_tof = "SimHit TOF : region" + regionLabel[region_num] + " layer " +
                                       layerLabel[layer_num] + " " + " ; Time of flight [ns] ; entries";
      std::string hist_label_for_tofMu = "SimHit TOF(Muon only) : region" + regionLabel[region_num] + " layer " +
                                         layerLabel[layer_num] + " " + " ; Time of flight [ns] ; entries";
      std::string hist_label_for_eloss = "SimHit energy loss : region" + regionLabel[region_num] + " layer " +
                                         layerLabel[layer_num] + " " + " ; Energy loss [eV] ; entries";
      std::string hist_label_for_elossMu = "SimHit energy loss(Muon only) : region" + regionLabel[region_num] +
                                           " layer " + layerLabel[layer_num] + " " + " ; Energy loss [eV] ; entries";

      double tof_min, tof_max;
      tof_min = 10;
      tof_max = 30;
      me0_sh_tof[region_num][layer_num] =
          ibooker.book1D(hist_name_for_tof.c_str(), hist_label_for_tof.c_str(), 40, tof_min, tof_max);
      me0_sh_tofMu[region_num][layer_num] =
          ibooker.book1D(hist_name_for_tofMu.c_str(), hist_label_for_tofMu.c_str(), 40, tof_min, tof_max);
      me0_sh_eloss[region_num][layer_num] =
          ibooker.book1D(hist_name_for_eloss.c_str(), hist_label_for_eloss.c_str(), 60, 0., 6000.);
      me0_sh_elossMu[region_num][layer_num] =
          ibooker.book1D(hist_name_for_elossMu.c_str(), hist_label_for_elossMu.c_str(), 60, 0., 6000.);
    }
  }
}

Member Data Documentation

MonitorElement* ME0HitsValidation::me0_sh_eloss[2][6]

private

Definition at line 20 of file ME0HitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* ME0HitsValidation::me0_sh_elossMu[2][6]

private

Definition at line 21 of file ME0HitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* ME0HitsValidation::me0_sh_tof[2][6]

private

Definition at line 18 of file ME0HitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* ME0HitsValidation::me0_sh_tofMu[2][6]

private

Definition at line 19 of file ME0HitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* ME0HitsValidation::me0_sh_tot_zr[2]

private

Definition at line 16 of file ME0HitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* ME0HitsValidation::me0_sh_xy[2][6]

private

Definition at line 14 of file ME0HitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* ME0HitsValidation::me0_sh_zr[2][6]

private

Definition at line 15 of file ME0HitsValidation.h.

Referenced by analyze(), and bookHistograms().

Int_t ME0HitsValidation::npart

private

Definition at line 23 of file ME0HitsValidation.h.