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 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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

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

Detailed Description

Definition at line 48 of file MuonSimHitsValidAnalyzer.h.

Member Typedef Documentation

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

Definition at line 53 of file MuonSimHitsValidAnalyzer.h.

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

Definition at line 54 of file MuonSimHitsValidAnalyzer.h.

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

Definition at line 55 of file MuonSimHitsValidAnalyzer.h.

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

Definition at line 56 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 15 of file MuonSimHitsValidAnalyzer.cc.

References costeta, DTHitsToken_, fName, getAllProvenances, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), globposx, globposy, checklumidiff::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;
}

MuonSimHitsValidAnalyzer::~MuonSimHitsValidAnalyzer ( )

override

Definition at line 72 of file MuonSimHitsValidAnalyzer.cc.

{
}

Member Function Documentation

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

override

keep track of number of events processed

get event id information

look at information available in the event

Definition at line 254 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 &    )

override

Definition at line 76 of file MuonSimHitsValidAnalyzer.cc.

References DQMStore::IBooker::book1D(), DQMStore::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, nullptr, and DQMStore::IBooker::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 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 309 of file MuonSimHitsValidAnalyzer.cc.

References funct::abs(), funct::cos(), costeta, gamEcalExtractorBlocks_cff::detector, dMuon, DTHitsToken_, PVValHelper::eta, JetChargeProducer_cfi::exp, MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), globposeta, globposphi, globposx, globposy, globposz, mps_fire::i, iden, DTGeometry::idToDetUnit(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), 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> 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;
}

Member Data Documentation

Float_t MuonSimHitsValidAnalyzer::costeta

private

Definition at line 86 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

unsigned int MuonSimHitsValidAnalyzer::count

private

private statistics information

Definition at line 130 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze().

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

private

Input tags.

Definition at line 91 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

std::string MuonSimHitsValidAnalyzer::fName

private

parameter information

Definition at line 69 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

bool MuonSimHitsValidAnalyzer::getAllProvenances

private

Definition at line 72 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::globposeta

private

Definition at line 88 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Float_t MuonSimHitsValidAnalyzer::globposphi

private

Definition at line 88 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Float_t MuonSimHitsValidAnalyzer::globposx

private

Definition at line 87 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::globposy

private

Definition at line 87 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::globposz

private

Definition at line 87 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

unsigned int MuonSimHitsValidAnalyzer::iden

private

Definition at line 78 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

std::string MuonSimHitsValidAnalyzer::label

private

Definition at line 71 of file MuonSimHitsValidAnalyzer.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::_sort_list(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::foundIn(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::fullFilename(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputEventContent(), MuonSimHitsValidAnalyzer(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), Vispa.Plugins.ConfigEditor.ToolDataAccessor.ToolDataAccessor::properties(), Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::properties(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::readConnections(), Vispa.Plugins.ConfigEditor.ToolDataAccessor.ToolDataAccessor::updateProcess(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::usedBy().

Int_t MuonSimHitsValidAnalyzer::layer

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

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

MonitorElement* MuonSimHitsValidAnalyzer::meAllDTHits

private

Definition at line 95 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meChamberOccup

private

Definition at line 123 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meCosTheta

private

Definition at line 125 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meEnergyLoss

private

Definition at line 98 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalEta

private

Definition at line 126 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalPhi

private

Definition at line 127 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsY

private

Definition at line 105 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYW0

private

Definition at line 113 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWm1

private

Definition at line 112 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWm2

private

Definition at line 111 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWp1

private

Definition at line 114 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsYWp2

private

Definition at line 115 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZ

private

Definition at line 104 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZW0

private

Definition at line 108 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWm1

private

Definition at line 107 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWm2

private

Definition at line 106 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWp1

private

Definition at line 109 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meGlobalXvsZWp2

private

Definition at line 110 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meHitRadius

private

Definition at line 124 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meLayerOccup

private

Definition at line 120 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meLocalXvsY

private

Definition at line 103 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meLocalXvsZ

private

Definition at line 102 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meLossMomIron

private

Definition at line 101 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meMomentumMB1

private

Definition at line 99 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meMomentumMB4

private

Definition at line 100 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meMuDTHits

private

Definition at line 96 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::mePathMuon

private

Definition at line 122 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meSectorOccup

private

Definition at line 118 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meStationOccup

private

Definition at line 117 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meSuperLOccup

private

Definition at line 119 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meToF

private

Definition at line 97 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meWheelOccup

private

Definition at line 116 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

MonitorElement* MuonSimHitsValidAnalyzer::meWireOccup

private

Definition at line 121 of file MuonSimHitsValidAnalyzer.h.

Referenced by bookHistograms(), and fillDT().

Float_t MuonSimHitsValidAnalyzer::mom1

private

Definition at line 85 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::mom4

private

Definition at line 85 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::mome1

private

Definition at line 85 of file MuonSimHitsValidAnalyzer.h.

Float_t MuonSimHitsValidAnalyzer::mome4

private

Definition at line 85 of file MuonSimHitsValidAnalyzer.h.

int MuonSimHitsValidAnalyzer::nRawGenPart

private

G4MC info.

Definition at line 76 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::nummu_CSC

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::nummu_DT

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::nummu_RPC

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

Referenced by MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::path

private

Definition at line 81 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Int_t MuonSimHitsValidAnalyzer::pathchamber

private

Definition at line 81 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Float_t MuonSimHitsValidAnalyzer::pow6

private

Definition at line 84 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

bool MuonSimHitsValidAnalyzer::printProvenanceInfo

private

Definition at line 73 of file MuonSimHitsValidAnalyzer.h.

Referenced by analyze(), and MuonSimHitsValidAnalyzer().

Float_t MuonSimHitsValidAnalyzer::radius

private

Definition at line 86 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::region

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Int_t MuonSimHitsValidAnalyzer::ring

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Referenced by geometryXMLparser.CSCAlignable::index().

Int_t MuonSimHitsValidAnalyzer::roll

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Int_t MuonSimHitsValidAnalyzer::sector

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

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

Float_t MuonSimHitsValidAnalyzer::sinteta

private

Definition at line 86 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT(), and MuonSimHitsValidAnalyzer().

Int_t MuonSimHitsValidAnalyzer::station

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

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

Int_t MuonSimHitsValidAnalyzer::subsector

private

Definition at line 80 of file MuonSimHitsValidAnalyzer.h.

Int_t MuonSimHitsValidAnalyzer::superlayer

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

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

Int_t MuonSimHitsValidAnalyzer::touch1

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Int_t MuonSimHitsValidAnalyzer::touch4

private

Definition at line 82 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().

Int_t MuonSimHitsValidAnalyzer::touche1

private

Definition at line 83 of file MuonSimHitsValidAnalyzer.h.

Int_t MuonSimHitsValidAnalyzer::touche4

private

Definition at line 83 of file MuonSimHitsValidAnalyzer.h.

int MuonSimHitsValidAnalyzer::verbosity

private

Definition at line 70 of file MuonSimHitsValidAnalyzer.h.

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

Int_t MuonSimHitsValidAnalyzer::wheel

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

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

Int_t MuonSimHitsValidAnalyzer::wire

private

Definition at line 79 of file MuonSimHitsValidAnalyzer.h.

Referenced by fillDT().