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 >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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...
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 50 of file MuonSimHitsValidAnalyzer.h.

Member Typedef Documentation

FloatVector

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

Definition at line 52 of file MuonSimHitsValidAnalyzer.h.

IntegerVector

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

Definition at line 53 of file MuonSimHitsValidAnalyzer.h.

LongIntegerVector

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

Definition at line 54 of file MuonSimHitsValidAnalyzer.h.

UnsigIntegerVector

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

Definition at line 55 of file MuonSimHitsValidAnalyzer.h.

Constructor & Destructor Documentation

MuonSimHitsValidAnalyzer()

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.

    : 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;
}

References costeta, DTHitsToken_, fName, 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.

~MuonSimHitsValidAnalyzer()

MuonSimHitsValidAnalyzer::~MuonSimHitsValidAnalyzer ( )

override

Definition at line 73 of file MuonSimHitsValidAnalyzer.cc.

{}

Member Function Documentation

analyze()

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 254 of file MuonSimHitsValidAnalyzer.cc.

                                                                                          {
  ++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;
}

References count, fillDT(), getAllProvenances, mps_fire::i, iEvent, nevt, printProvenanceInfo, and verbosity.

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 75 of file MuonSimHitsValidAnalyzer.cc.

                                                                      {
  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);
}

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

fillDT()

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

private

iterator to access containers

access the DT access the DT geometry

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 303 of file MuonSimHitsValidAnalyzer.cc.

                                                                                         {
  TString eventout;
  if (verbosity > 0)
    eventout = "\nGathering DT info:";
 
  edm::PSimHitContainer::const_iterator itHit;
 
  edm::ESHandle<DTGeometry> theDTGeometry;
  iSetup.get<MuonGeometryRecord>().get(theDTGeometry);
  if (!theDTGeometry.isValid()) {
    edm::LogWarning("MuonSimHitsValidAnalyzer::fillDT")
        << "Unable to find MuonGeometryRecord for the DTGeometry in event!";
    return;
  }
  const DTGeometry& theDTMuon(*theDTGeometry);
 
  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;
}

References funct::abs(), funct::cos(), costeta, hgcalTestNeighbor_cfi::detector, dMuon, DTHitsToken_, PVValHelper::eta, JetChargeProducer_cfi::exp, dqm::impl::MonitorElement::Fill(), edm::EventSetup::get(), get, globposeta, globposphi, globposx, globposy, globposz, mps_fire::i, iden, DTGeometry::idToDetUnit(), iEvent, 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().

Member Data Documentation

costeta

Float_t MuonSimHitsValidAnalyzer::costeta

private

Definition at line 84 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

count

unsigned int MuonSimHitsValidAnalyzer::count

private

private statistics information

Definition at line 128 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze().

DTHitsToken_

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

private

Input tags.

Definition at line 89 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

fName

std::string MuonSimHitsValidAnalyzer::fName

private

parameter information

Definition at line 67 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

getAllProvenances

bool MuonSimHitsValidAnalyzer::getAllProvenances

private

Definition at line 70 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze(), and MuonSimHitsValidAnalyzer().

globposeta

Float_t MuonSimHitsValidAnalyzer::globposeta

private

Definition at line 86 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

globposphi

Float_t MuonSimHitsValidAnalyzer::globposphi

private

Definition at line 86 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

globposx

Float_t MuonSimHitsValidAnalyzer::globposx

private

Definition at line 85 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

globposy

Float_t MuonSimHitsValidAnalyzer::globposy

private

Definition at line 85 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

globposz

Float_t MuonSimHitsValidAnalyzer::globposz

private

Definition at line 85 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

iden

unsigned int MuonSimHitsValidAnalyzer::iden

private

Definition at line 76 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

label

std::string MuonSimHitsValidAnalyzer::label

private

Definition at line 69 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

layer

Int_t MuonSimHitsValidAnalyzer::layer

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

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

meAllDTHits

MonitorElement* MuonSimHitsValidAnalyzer::meAllDTHits

private

Definition at line 93 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meChamberOccup

MonitorElement* MuonSimHitsValidAnalyzer::meChamberOccup

private

Definition at line 121 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meCosTheta

MonitorElement* MuonSimHitsValidAnalyzer::meCosTheta

private

Definition at line 123 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meEnergyLoss

MonitorElement* MuonSimHitsValidAnalyzer::meEnergyLoss

private

Definition at line 96 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalEta

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalEta

private

Definition at line 124 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalPhi

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalPhi

private

Definition at line 125 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsY

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsY

private

Definition at line 103 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsYW0

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYW0

private

Definition at line 111 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsYWm1

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWm1

private

Definition at line 110 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsYWm2

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWm2

private

Definition at line 109 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsYWp1

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWp1

private

Definition at line 112 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsYWp2

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWp2

private

Definition at line 113 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsZ

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZ

private

Definition at line 102 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsZW0

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZW0

private

Definition at line 106 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsZWm1

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWm1

private

Definition at line 105 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsZWm2

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWm2

private

Definition at line 104 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsZWp1

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWp1

private

Definition at line 107 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meGlobalXvsZWp2

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWp2

private

Definition at line 108 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meHitRadius

MonitorElement* MuonSimHitsValidAnalyzer::meHitRadius

private

Definition at line 122 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meLayerOccup

MonitorElement* MuonSimHitsValidAnalyzer::meLayerOccup

private

Definition at line 118 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meLocalXvsY

MonitorElement* MuonSimHitsValidAnalyzer::meLocalXvsY

private

Definition at line 101 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meLocalXvsZ

MonitorElement* MuonSimHitsValidAnalyzer::meLocalXvsZ

private

Definition at line 100 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meLossMomIron

MonitorElement* MuonSimHitsValidAnalyzer::meLossMomIron

private

Definition at line 99 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meMomentumMB1

MonitorElement* MuonSimHitsValidAnalyzer::meMomentumMB1

private

Definition at line 97 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meMomentumMB4

MonitorElement* MuonSimHitsValidAnalyzer::meMomentumMB4

private

Definition at line 98 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meMuDTHits

MonitorElement* MuonSimHitsValidAnalyzer::meMuDTHits

private

Definition at line 94 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

mePathMuon

MonitorElement* MuonSimHitsValidAnalyzer::mePathMuon

private

Definition at line 120 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meSectorOccup

MonitorElement* MuonSimHitsValidAnalyzer::meSectorOccup

private

Definition at line 116 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meStationOccup

MonitorElement* MuonSimHitsValidAnalyzer::meStationOccup

private

Definition at line 115 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meSuperLOccup

MonitorElement* MuonSimHitsValidAnalyzer::meSuperLOccup

private

Definition at line 117 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meToF

MonitorElement* MuonSimHitsValidAnalyzer::meToF

private

Definition at line 95 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meWheelOccup

MonitorElement* MuonSimHitsValidAnalyzer::meWheelOccup

private

Definition at line 114 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

meWireOccup

MonitorElement* MuonSimHitsValidAnalyzer::meWireOccup

private

Definition at line 119 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

mom1

Float_t MuonSimHitsValidAnalyzer::mom1

private

Definition at line 83 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

mom4

Float_t MuonSimHitsValidAnalyzer::mom4

private

Definition at line 83 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

mome1

Float_t MuonSimHitsValidAnalyzer::mome1

private

Definition at line 83 of file MuonSimHitsValidAnalyzer.h.

mome4

Float_t MuonSimHitsValidAnalyzer::mome4

private

Definition at line 83 of file MuonSimHitsValidAnalyzer.h.

nRawGenPart

int MuonSimHitsValidAnalyzer::nRawGenPart

private

G4MC info.

Definition at line 74 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

nummu_CSC

Int_t MuonSimHitsValidAnalyzer::nummu_CSC

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

nummu_DT

Int_t MuonSimHitsValidAnalyzer::nummu_DT

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

nummu_RPC

Int_t MuonSimHitsValidAnalyzer::nummu_RPC

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

path

Int_t MuonSimHitsValidAnalyzer::path

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

pathchamber

Int_t MuonSimHitsValidAnalyzer::pathchamber

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

pow6

Float_t MuonSimHitsValidAnalyzer::pow6

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

printProvenanceInfo

bool MuonSimHitsValidAnalyzer::printProvenanceInfo

private

Definition at line 71 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze(), and MuonSimHitsValidAnalyzer().

radius

Float_t MuonSimHitsValidAnalyzer::radius

private

Definition at line 84 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

region

Int_t MuonSimHitsValidAnalyzer::region

private

Definition at line 78 of file MuonSimHitsValidAnalyzer.h.

ring

Int_t MuonSimHitsValidAnalyzer::ring

private

Definition at line 78 of file MuonSimHitsValidAnalyzer.h.

Referenced by geometryXMLparser.CSCAlignable::index().

roll

Int_t MuonSimHitsValidAnalyzer::roll

private

Definition at line 78 of file MuonSimHitsValidAnalyzer.h.

sector

Int_t MuonSimHitsValidAnalyzer::sector

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

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

sinteta

Float_t MuonSimHitsValidAnalyzer::sinteta

private

Definition at line 84 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

station

Int_t MuonSimHitsValidAnalyzer::station

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

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

subsector

Int_t MuonSimHitsValidAnalyzer::subsector

private

Definition at line 78 of file MuonSimHitsValidAnalyzer.h.

superlayer

Int_t MuonSimHitsValidAnalyzer::superlayer

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

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

touch1

Int_t MuonSimHitsValidAnalyzer::touch1

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

touch4

Int_t MuonSimHitsValidAnalyzer::touch4

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

touche1

Int_t MuonSimHitsValidAnalyzer::touche1

private

Definition at line 81 of file MuonSimHitsValidAnalyzer.h.

touche4

Int_t MuonSimHitsValidAnalyzer::touche4

private

Definition at line 81 of file MuonSimHitsValidAnalyzer.h.

verbosity

int MuonSimHitsValidAnalyzer::verbosity

private

Definition at line 68 of file MuonSimHitsValidAnalyzer.h.

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

wheel

Int_t MuonSimHitsValidAnalyzer::wheel

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

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

wire

Int_t MuonSimHitsValidAnalyzer::wire

private

Definition at line 77 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().