MuonSimHitsValidAnalyzer Class Reference

#include <MuonSimHitsValidAnalyzer.h>

Inheritance diagram for MuonSimHitsValidAnalyzer:

Public Types
typedef std::vector< float >	FloatVector
typedef std::vector< int >	IntegerVector
typedef std::vector< long int >	LongIntegerVector
typedef std::vector< unsigned int >	UnsigIntegerVector
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

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	MuonSimHitsValidAnalyzer (const edm::ParameterSet &)
	~MuonSimHitsValidAnalyzer () 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 Member Functions
void	fillDT (const edm::Event &, const edm::EventSetup &)

Private Attributes
Float_t	costeta
unsigned int	count
	private statistics information More...
edm::EDGetTokenT < edm::PSimHitContainer >	DTHitsToken_
	Input tags. More...
std::string	fName
	parameter information More...
edm::ESGetToken< DTGeometry, MuonGeometryRecord >	geomToken_
bool	getAllProvenances
Float_t	globposeta
Float_t	globposphi
Float_t	globposx
Float_t	globposy
Float_t	globposz
unsigned int	iden
std::string	label
Int_t	layer
MonitorElement *	meAllDTHits
MonitorElement *	meChamberOccup
MonitorElement *	meCosTheta
MonitorElement *	meEnergyLoss
MonitorElement *	meGlobalEta
MonitorElement *	meGlobalPhi
MonitorElement *	meGlobalXvsY
MonitorElement *	meGlobalXvsYW0
MonitorElement *	meGlobalXvsYWm1
MonitorElement *	meGlobalXvsYWm2
MonitorElement *	meGlobalXvsYWp1
MonitorElement *	meGlobalXvsYWp2
MonitorElement *	meGlobalXvsZ
MonitorElement *	meGlobalXvsZW0
MonitorElement *	meGlobalXvsZWm1
MonitorElement *	meGlobalXvsZWm2
MonitorElement *	meGlobalXvsZWp1
MonitorElement *	meGlobalXvsZWp2
MonitorElement *	meHitRadius
MonitorElement *	meLayerOccup
MonitorElement *	meLocalXvsY
MonitorElement *	meLocalXvsZ
MonitorElement *	meLossMomIron
MonitorElement *	meMomentumMB1
MonitorElement *	meMomentumMB4
MonitorElement *	meMuDTHits
MonitorElement *	mePathMuon
MonitorElement *	meSectorOccup
MonitorElement *	meStationOccup
MonitorElement *	meSuperLOccup
MonitorElement *	meToF
MonitorElement *	meWheelOccup
MonitorElement *	meWireOccup
Float_t	mom1
Float_t	mom4
Float_t	mome1
Float_t	mome4
int	nRawGenPart
	G4MC info. More...
Int_t	nummu_CSC
Int_t	nummu_DT
Int_t	nummu_RPC
Int_t	path
Int_t	pathchamber
Float_t	pow6
bool	printProvenanceInfo
Float_t	radius
Int_t	region
Int_t	ring
Int_t	roll
Int_t	sector
Float_t	sinteta
Int_t	station
Int_t	subsector
Int_t	superlayer
Int_t	touch1
Int_t	touch4
Int_t	touche1
Int_t	touche4
int	verbosity
Int_t	wheel
Int_t	wire

Additional Inherited Members
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 49 of file MuonSimHitsValidAnalyzer.h.

Member Typedef Documentation

typedef std::vector<float> MuonSimHitsValidAnalyzer::FloatVector

Definition at line 51 of file MuonSimHitsValidAnalyzer.h.

typedef std::vector<int> MuonSimHitsValidAnalyzer::IntegerVector

Definition at line 52 of file MuonSimHitsValidAnalyzer.h.

typedef std::vector<long int> MuonSimHitsValidAnalyzer::LongIntegerVector

Definition at line 53 of file MuonSimHitsValidAnalyzer.h.

typedef std::vector<unsigned int> MuonSimHitsValidAnalyzer::UnsigIntegerVector

Definition at line 54 of file MuonSimHitsValidAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

get information from parameter set

get labels for input tags

print out Parameter Set information being used

Definition at line 14 of file MuonSimHitsValidAnalyzer.cc.

References costeta, DTHitsToken_, fName, geomToken_, getAllProvenances, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), globposx, globposy, cmsLHEtoEOSManager::l, label, mom1, mom4, nRawGenPart, nummu_CSC, nummu_DT, nummu_RPC, pow6, printProvenanceInfo, radius, sinteta, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

    : fName(""),
      verbosity(0),
      label(""),
      getAllProvenances(false),
      printProvenanceInfo(false),
      nRawGenPart(0),
      count(0)

{
  fName = iPSet.getUntrackedParameter<std::string>("Name");
  verbosity = iPSet.getUntrackedParameter<int>("Verbosity");
  label = iPSet.getParameter<std::string>("Label");
  edm::ParameterSet m_Prov = iPSet.getParameter<edm::ParameterSet>("ProvenanceLookup");
  getAllProvenances = m_Prov.getUntrackedParameter<bool>("GetAllProvenances");
  printProvenanceInfo = m_Prov.getUntrackedParameter<bool>("PrintProvenanceInfo");

  nRawGenPart = 0;

  DTHitsToken_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("DTHitsSrc"));


  if (verbosity) {
    Labels l;
    labelsForToken(DTHitsToken_, l);
    edm::LogInfo("MuonSimHitsValidAnalyzer::MuonSimHitsValidAnalyzer")
        << "\n===============================\n"
        << "Initialized as EDAnalyzer with parameter values:\n"
        << "    Name      = " << fName << "\n"
        << "    Verbosity = " << verbosity << "\n"
        << "    Label     = " << label << "\n"
        << "    GetProv   = " << getAllProvenances << "\n"
        << "    PrintProv = " << printProvenanceInfo
        << "\n"
        //    << "    CSCHitsSrc=  " <<CSCHitsSrc_.label()
        //    << ":" << CSCHitsSrc_.instance() << "\n"
        << "    DTHitsSrc =  " << l.module << ":" << l.productInstance
        << "\n"
        //     << "    RPCHitsSrc=  " <<RPCHitsSrc_.label()
        //     << ":" << RPCHitsSrc_.instance() << "\n"
        << "===============================\n";
  }

  pow6 = 1000000.0;
  mom4 = 0.;
  mom1 = 0;
  costeta = 0.;
  radius = 0;
  sinteta = 0.;
  globposx = 0.;
  globposy = 0;
  nummu_DT = 0;
  nummu_CSC = 0;
  nummu_RPC = 0;
  geomToken_ = esConsumes<DTGeometry, MuonGeometryRecord>();
}

MuonSimHitsValidAnalyzer::~MuonSimHitsValidAnalyzer ( )

override

Definition at line 74 of file MuonSimHitsValidAnalyzer.cc.

{}

Member Function Documentation

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

overridevirtual

keep track of number of events processed

get event id information

look at information available in the event

Reimplemented from DQMEDAnalyzer.

Definition at line 255 of file MuonSimHitsValidAnalyzer.cc.

References count, edm::EventID::event(), fillDT(), getAllProvenances, edm::Event::getAllStableProvenance(), mps_fire::i, edm::EventBase::id(), nevt, printProvenanceInfo, edm::EventID::run(), and verbosity.

                                                                                          {
  ++count;

  edm::RunNumber_t nrun = iEvent.id().run();
  edm::EventNumber_t nevt = iEvent.id().event();

  if (verbosity > 0) {
    edm::LogInfo("MuonSimHitsValidAnalyzer::analyze") << "Processing run " << nrun << ", event " << nevt;
  }

  if (getAllProvenances) {
    std::vector<const edm::StableProvenance*> AllProv;
    iEvent.getAllStableProvenance(AllProv);

    if (verbosity > 0)
      edm::LogInfo("MuonSimHitsValidAnalyzer::analyze") << "Number of Provenances = " << AllProv.size();

    if (printProvenanceInfo && (verbosity > 0)) {
      TString eventout("\nProvenance info:\n");

      for (unsigned int i = 0; i < AllProv.size(); ++i) {
        eventout += "\n       ******************************";
        eventout += "\n       Module       : ";
        eventout += AllProv[i]->moduleLabel();
        eventout += "\n       ProductID    : ";
        eventout += AllProv[i]->productID().id();
        eventout += "\n       ClassName    : ";
        eventout += AllProv[i]->className();
        eventout += "\n       InstanceName : ";
        eventout += AllProv[i]->productInstanceName();
        eventout += "\n       BranchName   : ";
        eventout += AllProv[i]->branchName();
      }
      eventout += "       ******************************\n";
      edm::LogInfo("MuonSimHitsValidAnalyzer::analyze") << eventout << "\n";
    }
  }

  fillDT(iEvent, iSetup);

  if (verbosity > 0)
    edm::LogInfo("MuonSimHitsValidAnalyzer::analyze") << "Done gathering data from event.";

  return;
}

void MuonSimHitsValidAnalyzer::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  iRun, edm::EventSetup const &    )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 76 of file MuonSimHitsValidAnalyzer.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), meAllDTHits, meChamberOccup, meCosTheta, meEnergyLoss, meGlobalEta, meGlobalPhi, meGlobalXvsY, meGlobalXvsYW0, meGlobalXvsYWm1, meGlobalXvsYWm2, meGlobalXvsYWp1, meGlobalXvsYWp2, meGlobalXvsZ, meGlobalXvsZW0, meGlobalXvsZWm1, meGlobalXvsZWm2, meGlobalXvsZWp1, meGlobalXvsZWp2, meHitRadius, meLayerOccup, meLocalXvsY, meLocalXvsZ, meLossMomIron, meMomentumMB1, meMomentumMB4, meMuDTHits, mePathMuon, meSectorOccup, meStationOccup, meSuperLOccup, meToF, meWheelOccup, meWireOccup, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                      {
  meAllDTHits = nullptr;
  meMuDTHits = nullptr;
  meToF = nullptr;
  meEnergyLoss = nullptr;
  meMomentumMB1 = nullptr;
  meMomentumMB4 = nullptr;
  meLossMomIron = nullptr;
  meLocalXvsZ = nullptr;
  meLocalXvsY = nullptr;
  meGlobalXvsZ = nullptr;
  meGlobalXvsY = nullptr;
  meGlobalXvsZWm2 = nullptr;
  meGlobalXvsZWm1 = nullptr;
  meGlobalXvsZW0 = nullptr;
  meGlobalXvsZWp1 = nullptr;
  meGlobalXvsZWp2 = nullptr;
  meGlobalXvsYWm2 = nullptr;
  meGlobalXvsYWm1 = nullptr;
  meGlobalXvsYW0 = nullptr;
  meGlobalXvsYWp1 = nullptr;
  meGlobalXvsYWp2 = nullptr;
  meWheelOccup = nullptr;
  meStationOccup = nullptr;
  meSectorOccup = nullptr;
  meSuperLOccup = nullptr;
  meLayerOccup = nullptr;
  meWireOccup = nullptr;
  mePathMuon = nullptr;
  meChamberOccup = nullptr;
  meHitRadius = nullptr;
  meCosTheta = nullptr;
  meGlobalEta = nullptr;
  meGlobalPhi = nullptr;

  Char_t histo_n[100];
  Char_t histo_t[100];

  iBooker.setCurrentFolder("MuonDTHitsV/DTHitsValidationTask");

  sprintf(histo_n, "Number_of_all_DT_hits");
  sprintf(histo_t, "Number_of_all_DT_hits");
  meAllDTHits = iBooker.book1D(histo_n, histo_t, 200, 1.0, 201.0);

  sprintf(histo_n, "Number_of_muon_DT_hits");
  sprintf(histo_t, "Number_of_muon_DT_hits");
  meMuDTHits = iBooker.book1D(histo_n, histo_t, 150, 1.0, 151.0);

  sprintf(histo_n, "Tof_of_hits ");
  sprintf(histo_t, "Tof_of_hits ");
  meToF = iBooker.book1D(histo_n, histo_t, 100, -0.5, 50.);

  sprintf(histo_n, "DT_energy_loss_keV");
  sprintf(histo_t, "DT_energy_loss_keV");
  meEnergyLoss = iBooker.book1D(histo_n, histo_t, 100, 0.0, 10.0);

  sprintf(histo_n, "Momentum_at_MB1");
  sprintf(histo_t, "Momentum_at_MB1");
  meMomentumMB1 = iBooker.book1D(histo_n, histo_t, 100, 10.0, 200.0);

  sprintf(histo_n, "Momentum_at_MB4");
  sprintf(histo_t, "Momentum_at_MB4");
  meMomentumMB4 = iBooker.book1D(histo_n, histo_t, 100, 10.0, 200.0);

  sprintf(histo_n, "Loss_of_muon_Momentum_in_Iron");
  sprintf(histo_t, "Loss_of_muon_Momentum_in_Iron");
  meLossMomIron = iBooker.book1D(histo_n, histo_t, 80, 0.0, 40.0);

  sprintf(histo_n, "Local_x-coord_vs_local_z-coord_of_muon_hit");
  sprintf(histo_t, "Local_x-coord_vs_local_z-coord_of_muon_hit");
  meLocalXvsZ = iBooker.book2D(histo_n, histo_t, 100, -150., 150., 100, -0.8, 0.8);

  sprintf(histo_n, "local_x-coord_vs_local_y-coord_of_muon_hit");
  sprintf(histo_t, "local_x-coord_vs_local_y-coord_of_muon_hit");
  meLocalXvsY = iBooker.book2D(histo_n, histo_t, 100, -150., 150., 100, -150., 150.);

  sprintf(histo_n, "Global_x-coord_vs_global_z-coord_of_muon_hit");
  sprintf(histo_t, "Global_x-coord_vs_global_z-coord_of_muon_hit");
  meGlobalXvsZ = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_y-coord_of_muon_hit");
  sprintf(histo_t, "Global_x-coord_vs_global_y-coord_of_muon_hit");
  meGlobalXvsY = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  // New histos

  sprintf(histo_n, "Global_x-coord_vs_global_z-coord_of_muon_hit_w-2");
  sprintf(histo_t, "Global_x-coord_vs_global_z-coord_of_muon_hit_w-2");
  meGlobalXvsZWm2 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_y-coord_of_muon_hit_w-2");
  sprintf(histo_t, "Global_x-coord_vs_global_y-coord_of_muon_hit_w-2");
  meGlobalXvsYWm2 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_z-coord_of_muon_hit_w-1");
  sprintf(histo_t, "Global_x-coord_vs_global_z-coord_of_muon_hit_w-1");
  meGlobalXvsZWm1 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_y-coord_of_muon_hit_w-1");
  sprintf(histo_t, "Global_x-coord_vs_global_y-coord_of_muon_hit_w-1");
  meGlobalXvsYWm1 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_z-coord_of_muon_hit_w0");
  sprintf(histo_t, "Global_x-coord_vs_global_z-coord_of_muon_hit_w0");
  meGlobalXvsZW0 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_y-coord_of_muon_hit_w0");
  sprintf(histo_t, "Global_x-coord_vs_global_y-coord_of_muon_hit_w0");
  meGlobalXvsYW0 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_z-coord_of_muon_hit_w1");
  sprintf(histo_t, "Global_x-coord_vs_global_z-coord_of_muon_hit_w1");
  meGlobalXvsZWp1 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_y-coord_of_muon_hit_w1");
  sprintf(histo_t, "Global_x-coord_vs_global_y-coord_of_muon_hit_w1");
  meGlobalXvsYWp1 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_z-coord_of_muon_hit_w2");
  sprintf(histo_t, "Global_x-coord_vs_global_z-coord_of_muon_hit_w2");
  meGlobalXvsZWp2 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  sprintf(histo_n, "Global_x-coord_vs_global_y-coord_of_muon_hit_w2");
  sprintf(histo_t, "Global_x-coord_vs_global_y-coord_of_muon_hit_w2");
  meGlobalXvsYWp2 = iBooker.book2D(histo_n, histo_t, 100, -800., 800., 100, -800., 800.);

  //

  sprintf(histo_n, "Wheel_occupancy");
  sprintf(histo_t, "Wheel_occupancy");
  meWheelOccup = iBooker.book1D(histo_n, histo_t, 10, -5.0, 5.0);

  sprintf(histo_n, "Station_occupancy");
  sprintf(histo_t, "Station_occupancy");
  meStationOccup = iBooker.book1D(histo_n, histo_t, 6, 0., 6.0);

  sprintf(histo_n, "Sector_occupancy");
  sprintf(histo_t, "Sector_occupancy");
  meSectorOccup = iBooker.book1D(histo_n, histo_t, 20, 0., 20.);

  sprintf(histo_n, "SuperLayer_occupancy");
  sprintf(histo_t, "SuperLayer_occupancy");
  meSuperLOccup = iBooker.book1D(histo_n, histo_t, 5, 0., 5.);

  sprintf(histo_n, "Layer_occupancy");
  sprintf(histo_t, "Layer_occupancy");
  meLayerOccup = iBooker.book1D(histo_n, histo_t, 6, 0., 6.);

  sprintf(histo_n, "Wire_occupancy");
  sprintf(histo_t, "Wire_occupancy");
  meWireOccup = iBooker.book1D(histo_n, histo_t, 100, 0., 100.);

  sprintf(histo_n, "path_followed_by_muon");
  sprintf(histo_t, "path_followed_by_muon");
  mePathMuon = iBooker.book1D(histo_n, histo_t, 160, 0., 160.);

  sprintf(histo_n, "chamber_occupancy");
  sprintf(histo_t, "chamber_occupancy");
  meChamberOccup = iBooker.book1D(histo_n, histo_t, 251, 0., 251.);

  sprintf(histo_n, "radius_of_hit");
  sprintf(histo_t, "radius_of_hit");
  meHitRadius = iBooker.book1D(histo_n, histo_t, 100, 0., 1200.);

  sprintf(histo_n, "costheta_of_hit");
  sprintf(histo_t, "costheta_of_hit");
  meCosTheta = iBooker.book1D(histo_n, histo_t, 100, -1., 1.);

  sprintf(histo_n, "global_eta_of_hit");
  sprintf(histo_t, "global_eta_of_hit");
  meGlobalEta = iBooker.book1D(histo_n, histo_t, 60, -2.7, 2.7);

  sprintf(histo_n, "global_phi_of_hit");
  sprintf(histo_t, "global_phi_of_hit");
  meGlobalPhi = iBooker.book1D(histo_n, histo_t, 60, -3.14, 3.14);
}

void MuonSimHitsValidAnalyzer::fillDT ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

iterator to access containers

access the DT

get DT information

cycle through container

create a DetId from the detUnitId

check that expected detector is returned

get the GeomDetUnit from the geometry using theDetUnitID

get the Surface of the hit (knows how to go from local <-> global)

gather necessary information

Muon Momentum at MB1

Muon Momentum at MB4 & Loss of Muon Momentum in Iron (between MB1 and MB4)

X-Local Coordinate vs Z-Local Coordinate

X-Local Coordinate vs Y-Local Coordinate

Global Coordinates

Z-Global Coordinate vs X-Global Coordinate Y-Global Coordinate vs X-Global Coordinate

Definition at line 304 of file MuonSimHitsValidAnalyzer.cc.

References funct::abs(), funct::cos(), costeta, DetId::det(), dMuon, DTHitsToken_, PVValHelper::eta, funct::exp(), dqm::impl::MonitorElement::Fill(), geomToken_, edm::Event::getByToken(), edm::EventSetup::getHandle(), globposeta, globposphi, globposx, globposy, globposz, mps_fire::i, iden, DTGeometry::idToDetUnit(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), dqmiolumiharvest::j, layer, meAllDTHits, meChamberOccup, meCosTheta, meEnergyLoss, meGlobalEta, meGlobalPhi, meGlobalXvsY, meGlobalXvsYW0, meGlobalXvsYWm1, meGlobalXvsYWm2, meGlobalXvsYWp1, meGlobalXvsYWp2, meGlobalXvsZ, meGlobalXvsZW0, meGlobalXvsZWm1, meGlobalXvsZWm2, meGlobalXvsZWp1, meGlobalXvsZWp2, meHitRadius, meLayerOccup, meLocalXvsY, meLocalXvsZ, meLossMomIron, meMomentumMB1, meMomentumMB4, meMuDTHits, mePathMuon, meSectorOccup, meStationOccup, meSuperLOccup, meToF, meWheelOccup, meWireOccup, mom1, mom4, nummu_DT, path, pathchamber, phi, pow6, radius, sdMuonDT, sector, funct::sin(), sinteta, station, subdetector, superlayer, GeomDet::surface(), touch1, touch4, verbosity, wheel, wire, and z.

Referenced by analyze().

                                                                                         {
  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering DT info:";

  edm::PSimHitContainer::const_iterator itHit;

  edm::ESHandle<DTGeometry> hGeom = iSetup.getHandle(geomToken_);
  if (!hGeom.isValid()) {
    edm::LogWarning("MuonSimHitsValidAnalyzer::fillDT")
        << "Unable to find MuonGeometryRecord for the DTGeometry in event!";
    return;
  }
  const DTGeometry& theDTMuon(*hGeom);

  edm::Handle<edm::PSimHitContainer> MuonDTContainer;
  iEvent.getByToken(DTHitsToken_, MuonDTContainer);
  if (!MuonDTContainer.isValid()) {
    edm::LogWarning("MuonSimHitsValidAnalyzer::fillDT") << "Unable to find MuonDTHits in event!";
    return;
  }

  touch1 = 0;
  touch4 = 0;
  nummu_DT = 0;

  meAllDTHits->Fill(MuonDTContainer->size());

  int i = 0, j = 0;
  for (itHit = MuonDTContainer->begin(); itHit != MuonDTContainer->end(); ++itHit) {
    ++i;

    DetId theDetUnitId(itHit->detUnitId());
    int detector = theDetUnitId.det();
    int subdetector = theDetUnitId.subdetId();

    if ((detector == dMuon) && (subdetector == sdMuonDT)) {
      const GeomDetUnit* theDet = theDTMuon.idToDetUnit(theDetUnitId);

      if (!theDet) {
        edm::LogWarning("MuonSimHitsValidAnalyzer::fillDT") << "Unable to get GeomDetUnit from theDTMuon for hit " << i;
        continue;
      }

      ++j;

      const BoundPlane& bsurf = theDet->surface();


      if (abs(itHit->particleType()) == 13) {
        nummu_DT++;
        meToF->Fill(itHit->tof());
        meEnergyLoss->Fill(itHit->energyLoss() * pow6);

        iden = itHit->detUnitId();

        wheel = ((iden >> 15) & 0x7) - 3;
        station = ((iden >> 22) & 0x7);
        sector = ((iden >> 18) & 0xf);
        superlayer = ((iden >> 13) & 0x3);
        layer = ((iden >> 10) & 0x7);
        wire = ((iden >> 3) & 0x7f);

        meWheelOccup->Fill((float)wheel);
        meStationOccup->Fill((float)station);
        meSectorOccup->Fill((float)sector);
        meSuperLOccup->Fill((float)superlayer);
        meLayerOccup->Fill((float)layer);
        meWireOccup->Fill((float)wire);

        // Define a quantity to take into account station, splayer and layer being hit.
        path = (station - 1) * 40 + superlayer * 10 + layer;
        mePathMuon->Fill((float)path);

        // Define a quantity to take into chamber being hit.
        pathchamber = (wheel + 2) * 50 + (station - 1) * 12 + sector;
        meChamberOccup->Fill((float)pathchamber);

        if (station == 1) {
          if (touch1 == 0) {
            mom1 = itHit->pabs();
            meMomentumMB1->Fill(mom1);
            touch1 = 1;
          }
        }

        if (station == 4) {
          if (touch4 == 0) {
            mom4 = itHit->pabs();
            touch4 = 1;
            meMomentumMB4->Fill(mom4);
            if (touch1 == 1) {
              meLossMomIron->Fill(mom1 - mom4);
            }
          }
        }

        meLocalXvsZ->Fill(itHit->localPosition().x(), itHit->localPosition().z());

        meLocalXvsY->Fill(itHit->localPosition().x(), itHit->localPosition().y());


        globposz = bsurf.toGlobal(itHit->localPosition()).z();
        globposeta = bsurf.toGlobal(itHit->localPosition()).eta();
        globposphi = bsurf.toGlobal(itHit->localPosition()).phi();

        radius = globposz * (1. + exp(-2. * globposeta)) / (1. - exp(-2. * globposeta));

        costeta = (1. - exp(-2. * globposeta)) / (1. + exp(-2. * globposeta));
        sinteta = 2. * exp(-globposeta) / (1. + exp(-2. * globposeta));

        globposx = radius * sinteta * cos(globposphi);
        globposy = radius * sinteta * sin(globposphi);

        meGlobalXvsZ->Fill(globposz, globposx);
        meGlobalXvsY->Fill(globposx, globposy);

        //  New Histos
        if (wheel == -2) {
          meGlobalXvsZWm2->Fill(globposz, globposx);
          meGlobalXvsYWm2->Fill(globposx, globposy);
        }
        if (wheel == -1) {
          meGlobalXvsZWm1->Fill(globposz, globposx);
          meGlobalXvsYWm1->Fill(globposx, globposy);
        }
        if (wheel == 0) {
          meGlobalXvsZW0->Fill(globposz, globposx);
          meGlobalXvsYW0->Fill(globposx, globposy);
        }
        if (wheel == 1) {
          meGlobalXvsZWp1->Fill(globposz, globposx);
          meGlobalXvsYWp1->Fill(globposx, globposy);
        }
        if (wheel == 2) {
          meGlobalXvsZWp2->Fill(globposz, globposx);
          meGlobalXvsYWp2->Fill(globposx, globposy);
        }
        //
        meHitRadius->Fill(radius);
        meCosTheta->Fill(costeta);
        meGlobalEta->Fill(globposeta);
        meGlobalPhi->Fill(globposphi);
      }
    } else {
      edm::LogWarning("MuonSimHitsValidAnalyzer::fillDT")
          << "MuonDT PSimHit " << i << " is expected to be (det,subdet) = (" << dMuon << "," << sdMuonDT
          << "); value returned is: (" << detector << "," << subdetector << ")";
      continue;
    }
  }

  if (verbosity > 1) {
    eventout += "\n          Number of DT muon Hits collected:......... ";
    eventout += j;
  }
  meMuDTHits->Fill((float)nummu_DT);

  if (verbosity > 0)
    edm::LogInfo("MuonSimHitsValidAnalyzer::fillDT") << eventout << "\n";
  return;
}

Member Data Documentation

Float_t MuonSimHitsValidAnalyzer::costeta

private

Definition at line 83 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

unsigned int MuonSimHitsValidAnalyzer::count

private

private statistics information

Definition at line 128 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> MuonSimHitsValidAnalyzer::DTHitsToken_

private

Input tags.

Definition at line 88 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

std::string MuonSimHitsValidAnalyzer::fName

private

parameter information

Definition at line 66 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

edm::ESGetToken<DTGeometry, MuonGeometryRecord> MuonSimHitsValidAnalyzer::geomToken_

private

Definition at line 89 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

bool MuonSimHitsValidAnalyzer::getAllProvenances

private

Definition at line 69 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::globposeta

private

Definition at line 85 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Float_t MuonSimHitsValidAnalyzer::globposphi

private

Definition at line 85 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Float_t MuonSimHitsValidAnalyzer::globposx

private

Definition at line 84 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::globposy

private

Definition at line 84 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::globposz

private

Definition at line 84 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

unsigned int MuonSimHitsValidAnalyzer::iden

private

Definition at line 75 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

std::string MuonSimHitsValidAnalyzer::label

private

Definition at line 68 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::layer

private

Definition at line 76 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), geometryXMLparser.DTAlignable::index(), and geometryXMLparser.CSCAlignable::index().

MonitorElement* MuonSimHitsValidAnalyzer::meAllDTHits

private

Definition at line 93 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meChamberOccup

private

Definition at line 121 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meCosTheta

private

Definition at line 123 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meEnergyLoss

private

Definition at line 96 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalEta

private

Definition at line 124 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalPhi

private

Definition at line 125 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsY

private

Definition at line 103 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYW0

private

Definition at line 111 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWm1

private

Definition at line 110 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWm2

private

Definition at line 109 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWp1

private

Definition at line 112 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWp2

private

Definition at line 113 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZ

private

Definition at line 102 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZW0

private

Definition at line 106 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWm1

private

Definition at line 105 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWm2

private

Definition at line 104 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWp1

private

Definition at line 107 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWp2

private

Definition at line 108 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meHitRadius

private

Definition at line 122 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meLayerOccup

private

Definition at line 118 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meLocalXvsY

private

Definition at line 101 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meLocalXvsZ

private

Definition at line 100 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meLossMomIron

private

Definition at line 99 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meMomentumMB1

private

Definition at line 97 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meMomentumMB4

private

Definition at line 98 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meMuDTHits

private

Definition at line 94 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::mePathMuon

private

Definition at line 120 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meSectorOccup

private

Definition at line 116 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meStationOccup

private

Definition at line 115 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meSuperLOccup

private

Definition at line 117 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meToF

private

Definition at line 95 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meWheelOccup

private

Definition at line 114 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meWireOccup

private

Definition at line 119 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

Float_t MuonSimHitsValidAnalyzer::mom1

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::mom4

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::mome1

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

Float_t MuonSimHitsValidAnalyzer::mome4

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

int MuonSimHitsValidAnalyzer::nRawGenPart

private

G4MC info.

Definition at line 73 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::nummu_CSC

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::nummu_DT

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::nummu_RPC

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::path

private

Definition at line 78 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Int_t MuonSimHitsValidAnalyzer::pathchamber

private

Definition at line 78 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Float_t MuonSimHitsValidAnalyzer::pow6

private

Definition at line 81 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

bool MuonSimHitsValidAnalyzer::printProvenanceInfo

private

Definition at line 70 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::radius

private

Definition at line 83 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::region

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

Int_t MuonSimHitsValidAnalyzer::ring

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

Referenced by geometryXMLparser.CSCAlignable::index().

Int_t MuonSimHitsValidAnalyzer::roll

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

Int_t MuonSimHitsValidAnalyzer::sector

private

Definition at line 76 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and geometryXMLparser.DTAlignable::index().

Float_t MuonSimHitsValidAnalyzer::sinteta

private

Definition at line 83 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::station

private

Definition at line 76 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), geometryXMLparser.DTAlignable::index(), and geometryXMLparser.CSCAlignable::index().

Int_t MuonSimHitsValidAnalyzer::subsector

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

Int_t MuonSimHitsValidAnalyzer::superlayer

private

Definition at line 76 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and geometryXMLparser.DTAlignable::index().

Int_t MuonSimHitsValidAnalyzer::touch1

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Int_t MuonSimHitsValidAnalyzer::touch4

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Int_t MuonSimHitsValidAnalyzer::touche1

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Int_t MuonSimHitsValidAnalyzer::touche4

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

int MuonSimHitsValidAnalyzer::verbosity

private

Definition at line 67 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze(), fillDT(), and MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::wheel

private

Definition at line 76 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and geometryXMLparser.DTAlignable::index().

Int_t MuonSimHitsValidAnalyzer::wire

private

Definition at line 76 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().