AlignmentMonitorMuonSystemMap1D Class Reference

#include <AlignmentMonitorMuonSystemMap1D.h>

Inheritance diagram for AlignmentMonitorMuonSystemMap1D:

Public Member Functions
void	afterAlignment (const edm::EventSetup &iSetup)
	AlignmentMonitorMuonSystemMap1D (const edm::ParameterSet &cfg)
void	book ()
	Book or retrieve histograms; MUST be reimplemented. More...
void	event (const edm::Event &iEvent, const edm::EventSetup &iSetup, const ConstTrajTrackPairCollection &iTrajTracks)
	Called for each event (by "run()"): may be reimplemented. More...
	~AlignmentMonitorMuonSystemMap1D ()
Public Member Functions inherited from AlignmentMonitorBase
	AlignmentMonitorBase (const edm::ParameterSet &cfg, std::string name)
	Constructor. More...
void	beginOfJob (AlignableTracker *pTracker, AlignableMuon *pMuon, AlignmentParameterStore *pStore)
	Called at beginning of job: don't reimplement. More...
void	duringLoop (const edm::Event &iEvent, const edm::EventSetup &iSetup, const ConstTrajTrackPairCollection &iTrajTracks)
	Called for each event: don't reimplement. More...
void	endOfJob ()
	Called at end of processing: don't implement. More...
void	endOfLoop (const edm::EventSetup &iSetup)
	Called at end of loop: don't reimplement. More...
void	startingNewLoop ()
	Called at beginning of loop: don't reimplement. More...
virtual	~AlignmentMonitorBase ()
	Destructor. More...

Private Member Functions
std::string	num02d (int num)

Private Attributes
bool	m_allowTIDTEC
long	m_counter_13numhits
long	m_counter_2numhits
long	m_counter_csc
long	m_counter_cscnumhits
long	m_counter_dt
long	m_counter_event
long	m_counter_track
long	m_counter_trackokay
long	m_counter_trackpt
MuonSystemMapPlot1D *	m_CSCvsphi_me11 [2]
MuonSystemMapPlot1D *	m_CSCvsphi_me12 [2]
MuonSystemMapPlot1D *	m_CSCvsphi_me13 [2]
MuonSystemMapPlot1D *	m_CSCvsphi_me14 [2]
MuonSystemMapPlot1D *	m_CSCvsphi_me21 [2]
MuonSystemMapPlot1D *	m_CSCvsphi_me22 [2]
MuonSystemMapPlot1D *	m_CSCvsphi_me31 [2]
MuonSystemMapPlot1D *	m_CSCvsphi_me32 [2]
MuonSystemMapPlot1D *	m_CSCvsphi_me41 [2]
MuonSystemMapPlot1D *	m_CSCvsphi_me42 [2]
MuonSystemMapPlot1D *	m_CSCvsr_me1 [2][36]
MuonSystemMapPlot1D *	m_CSCvsr_me2 [2][36]
MuonSystemMapPlot1D *	m_CSCvsr_me3 [2][36]
MuonSystemMapPlot1D *	m_CSCvsr_me4 [2][36]
MuonSystemMapPlot1D *	m_DTvsphi_station1 [5]
MuonSystemMapPlot1D *	m_DTvsphi_station2 [5]
MuonSystemMapPlot1D *	m_DTvsphi_station3 [5]
MuonSystemMapPlot1D *	m_DTvsphi_station4 [5]
MuonSystemMapPlot1D *	m_DTvsz_station1 [12]
MuonSystemMapPlot1D *	m_DTvsz_station2 [12]
MuonSystemMapPlot1D *	m_DTvsz_station3 [12]
MuonSystemMapPlot1D *	m_DTvsz_station4 [14]
double	m_maxTrackerRedChi2
double	m_maxTrackPt
int	m_minCSCHits
int	m_minDT13Hits
int	m_minDT2Hits
int	m_minTrackerHits
double	m_minTrackPt
std::vector < MuonSystemMapPlot1D * >	m_plots

Friends
class	MuonSystemMapPlot1D

Additional Inherited Members
Public Types inherited from AlignmentMonitorBase
typedef std::pair< const Trajectory *, const reco::Track * >	ConstTrajTrackPair
typedef std::vector < ConstTrajTrackPair >	ConstTrajTrackPairCollection
Protected Member Functions inherited from AlignmentMonitorBase
TH1F *	book1D (std::string dir, std::string name, std::string title, int nchX, double lowX, double highX)
TH2F *	book2D (std::string dir, std::string name, std::string title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
TProfile *	bookProfile (std::string dir, std::string name, std::string title, int nchX, double lowX, double highX, int nchY=1, double lowY=0., double highY=0., const char *option="s")
TFileDirectory *	directory (std::string dir)
int	iteration ()
AlignableMuon *	pMuon ()
AlignableNavigator *	pNavigator ()
AlignmentParameterStore *	pStore ()
AlignableTracker *	pTracker ()

Detailed Description

Date:

2009/08/29 18:18:07

Revision:

1.2

Author

J. Pivarski - Texas A&M University pivar.nosp@m.ski@.nosp@m.physi.nosp@m.cs.t.nosp@m.amu.e.nosp@m.du

Definition at line 28 of file AlignmentMonitorMuonSystemMap1D.h.

Constructor & Destructor Documentation

AlignmentMonitorMuonSystemMap1D::AlignmentMonitorMuonSystemMap1D

(

const edm::ParameterSet &

cfg

)

Definition at line 27 of file AlignmentMonitorMuonSystemMap1D.cc.

   : AlignmentMonitorBase(cfg, "AlignmentMonitorMuonSystemMap1D")
   , m_minTrackPt(cfg.getParameter<double>("minTrackPt"))
   , m_maxTrackPt(cfg.getParameter<double>("maxTrackPt"))
   , m_minTrackerHits(cfg.getParameter<int>("minTrackerHits"))
   , m_maxTrackerRedChi2(cfg.getParameter<double>("maxTrackerRedChi2"))
   , m_allowTIDTEC(cfg.getParameter<bool>("allowTIDTEC"))
   , m_minDT13Hits(cfg.getParameter<int>("minDT13Hits"))
   , m_minDT2Hits(cfg.getParameter<int>("minDT2Hits"))
   , m_minCSCHits(cfg.getParameter<int>("minCSCHits"))
{}

AlignmentMonitorMuonSystemMap1D::~AlignmentMonitorMuonSystemMap1D ( )

inline

Definition at line 31 of file AlignmentMonitorMuonSystemMap1D.h.

{};

Member Function Documentation

void AlignmentMonitorMuonSystemMap1D::afterAlignment ( const edm::EventSetup &  iSetup )

virtual

Called after updating AlignableTracker and AlignableMuon (by "endOfLoop()"): may be reimplemented

Reimplemented from AlignmentMonitorBase.

Definition at line 296 of file AlignmentMonitorMuonSystemMap1D.cc.

References gather_cfg::cout, m_counter_13numhits, m_counter_2numhits, m_counter_csc, m_counter_cscnumhits, m_counter_dt, m_counter_event, m_counter_track, m_counter_trackokay, and m_counter_trackpt.

                                                                                {
   std::cout << "monitor m_counter_event = " << m_counter_event << std::endl;
   std::cout << "monitor m_counter_track = " << m_counter_track << std::endl;
   std::cout << "monitor m_counter_trackpt = " << m_counter_trackpt << std::endl;
   std::cout << "monitor m_counter_trackokay = " << m_counter_trackokay << std::endl;
   std::cout << "monitor m_counter_dt = " << m_counter_dt << std::endl;
   std::cout << "monitor m_counter_13numhits = " << m_counter_13numhits << std::endl;
   std::cout << "monitor m_counter_2numhits = " << m_counter_2numhits << std::endl;
   std::cout << "monitor m_counter_csc = " << m_counter_csc << std::endl;
   std::cout << "monitor m_counter_cscnumhits = " << m_counter_cscnumhits << std::endl;
}

void AlignmentMonitorMuonSystemMap1D::book ( )

virtual

Book or retrieve histograms; MUST be reimplemented.

Implements AlignmentMonitorBase.

Definition at line 69 of file AlignmentMonitorMuonSystemMap1D.cc.

References Reference_intrackfit_cff::endcap, m_counter_13numhits, m_counter_2numhits, m_counter_csc, m_counter_cscnumhits, m_counter_dt, m_counter_event, m_counter_track, m_counter_trackokay, m_counter_trackpt, m_CSCvsphi_me11, m_CSCvsphi_me12, m_CSCvsphi_me13, m_CSCvsphi_me14, m_CSCvsphi_me21, m_CSCvsphi_me22, m_CSCvsphi_me31, m_CSCvsphi_me32, m_CSCvsphi_me41, m_CSCvsphi_me42, m_CSCvsr_me1, m_CSCvsr_me2, m_CSCvsr_me3, m_CSCvsr_me4, m_DTvsphi_station1, m_DTvsphi_station2, m_DTvsphi_station3, m_DTvsphi_station4, m_DTvsz_station1, m_DTvsz_station2, m_DTvsz_station3, m_DTvsz_station4, M_PI, m_plots, MuonSystemMapPlot1D, and num02d().

                                           {
  for (int sector = 1;  sector <= 14;  sector++) {
    if (sector <= 12) {
      m_DTvsz_station1[sector-1] = new MuonSystemMapPlot1D(std::string("DTvsz_st1sec") + num02d(sector), this, 60, -660., 660., true);  m_plots.push_back(m_DTvsz_station1[sector-1]);
      m_DTvsz_station2[sector-1] = new MuonSystemMapPlot1D(std::string("DTvsz_st2sec") + num02d(sector), this, 60, -660., 660., true);  m_plots.push_back(m_DTvsz_station2[sector-1]);
      m_DTvsz_station3[sector-1] = new MuonSystemMapPlot1D(std::string("DTvsz_st3sec") + num02d(sector), this, 60, -660., 660., true);  m_plots.push_back(m_DTvsz_station3[sector-1]);
    }
    m_DTvsz_station4[sector-1] = new MuonSystemMapPlot1D(std::string("DTvsz_st4sec") + num02d(sector), this, 60, -660., 660., false);  m_plots.push_back(m_DTvsz_station4[sector-1]);
  }

  for (int endcap = 1;  endcap <= 2;  endcap++) {
    for (int chamber = 1;  chamber <= 36;  chamber++) {
      m_CSCvsr_me1[endcap-1][chamber-1] = new MuonSystemMapPlot1D(std::string("CSCvsr_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("1ch") + num02d(chamber), this, 60, 100., 700., false);  m_plots.push_back(m_CSCvsr_me1[endcap-1][chamber-1]);
      m_CSCvsr_me2[endcap-1][chamber-1] = new MuonSystemMapPlot1D(std::string("CSCvsr_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("2ch") + num02d(chamber), this, 60, 100., 700., false);  m_plots.push_back(m_CSCvsr_me2[endcap-1][chamber-1]);
      m_CSCvsr_me3[endcap-1][chamber-1] = new MuonSystemMapPlot1D(std::string("CSCvsr_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("3ch") + num02d(chamber), this, 60, 100., 700., false);  m_plots.push_back(m_CSCvsr_me3[endcap-1][chamber-1]);
      m_CSCvsr_me4[endcap-1][chamber-1] = new MuonSystemMapPlot1D(std::string("CSCvsr_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("4ch") + num02d(chamber), this, 60, 100., 700., false);  m_plots.push_back(m_CSCvsr_me4[endcap-1][chamber-1]);
    }
  }

  for (int wheel = -2;  wheel <= 2;  wheel++) {
    std::string s_wheel;
    if (wheel == -2)      s_wheel = std::string("A");
    else if (wheel == -1) s_wheel = std::string("B");
    else if (wheel == 0)  s_wheel = std::string("C");
    else if (wheel == +1) s_wheel = std::string("D");
    else if (wheel == +2) s_wheel = std::string("E");

    m_DTvsphi_station1[wheel+2] = new MuonSystemMapPlot1D(std::string("DTvsphi_st1wh") + s_wheel, this, 180, -M_PI, M_PI, true);  m_plots.push_back(m_DTvsphi_station1[wheel+2]);
    m_DTvsphi_station2[wheel+2] = new MuonSystemMapPlot1D(std::string("DTvsphi_st2wh") + s_wheel, this, 180, -M_PI, M_PI, true);  m_plots.push_back(m_DTvsphi_station2[wheel+2]);
    m_DTvsphi_station3[wheel+2] = new MuonSystemMapPlot1D(std::string("DTvsphi_st3wh") + s_wheel, this, 180, -M_PI, M_PI, true);  m_plots.push_back(m_DTvsphi_station3[wheel+2]);
    m_DTvsphi_station4[wheel+2] = new MuonSystemMapPlot1D(std::string("DTvsphi_st4wh") + s_wheel, this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_DTvsphi_station4[wheel+2]);
  }

  for (int endcap = 1;  endcap <= 2;  endcap++) {
    m_CSCvsphi_me11[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("11"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me11[endcap-1]);
    m_CSCvsphi_me12[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("12"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me12[endcap-1]);
    m_CSCvsphi_me13[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("13"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me13[endcap-1]);
    m_CSCvsphi_me14[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("14"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me14[endcap-1]);
    m_CSCvsphi_me21[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("21"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me21[endcap-1]);
    m_CSCvsphi_me22[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("22"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me22[endcap-1]);
    m_CSCvsphi_me31[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("31"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me31[endcap-1]);
    m_CSCvsphi_me32[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("32"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me32[endcap-1]);
    m_CSCvsphi_me41[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("41"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me41[endcap-1]);
    m_CSCvsphi_me42[endcap-1] = new MuonSystemMapPlot1D(std::string("CSCvsphi_me") + (endcap == 1 ? std::string("p") : std::string("m")) + std::string("42"), this, 180, -M_PI, M_PI, false);  m_plots.push_back(m_CSCvsphi_me42[endcap-1]);
  }

  m_counter_event = 0;
  m_counter_track = 0;
  m_counter_trackpt = 0;
  m_counter_trackokay = 0;
  m_counter_dt = 0;
  m_counter_13numhits = 0;
  m_counter_2numhits = 0;
  m_counter_csc = 0;
  m_counter_cscnumhits = 0;
}

void AlignmentMonitorMuonSystemMap1D::event ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup, const ConstTrajTrackPairCollection &  iTrajTracks  )

virtual

Called for each event (by "run()"): may be reimplemented.

Reimplemented from AlignmentMonitorBase.

Definition at line 126 of file AlignmentMonitorMuonSystemMap1D.cc.

References MuonResidualsFromTrack::chamberIds(), MuonResidualsFromTrack::chamberResidual(), DeDxDiscriminatorTools::charge(), reco::TrackBase::charge(), MuonChamberResidual::chi2(), ExpressReco_HICollisions_FallBack::chi2, MuonResidualsFromTrack::contains_TIDTEC(), MuonSubdetId::CSC, MuonSubdetId::DT, Reference_intrackfit_cff::endcap, MuonSystemMapPlot1D::fill_dxdz(), MuonSystemMapPlot1D::fill_dydz(), MuonSystemMapPlot1D::fill_x(), MuonSystemMapPlot1D::fill_y(), edm::EventSetup::get(), MuonChamberResidual::global_residual(), MuonChamberResidual::global_resslope(), MuonChamberResidual::global_trackpos(), ExpressReco_HICollisions_FallBack::id, MuonChamberResidual::kCSC, MuonChamberResidual::kDT13, MuonChamberResidual::kDT2, m_allowTIDTEC, m_counter_13numhits, m_counter_2numhits, m_counter_csc, m_counter_cscnumhits, m_counter_dt, m_counter_event, m_counter_track, m_counter_trackokay, m_counter_trackpt, m_CSCvsphi_me11, m_CSCvsphi_me12, m_CSCvsphi_me13, m_CSCvsphi_me14, m_CSCvsphi_me21, m_CSCvsphi_me22, m_CSCvsphi_me31, m_CSCvsphi_me32, m_CSCvsphi_me41, m_CSCvsphi_me42, m_CSCvsr_me1, m_CSCvsr_me2, m_CSCvsr_me3, m_CSCvsr_me4, m_DTvsphi_station1, m_DTvsphi_station2, m_DTvsphi_station3, m_DTvsphi_station4, m_DTvsz_station1, m_DTvsz_station2, m_DTvsz_station3, m_DTvsz_station4, m_maxTrackerRedChi2, m_maxTrackPt, m_minCSCHits, m_minDT13Hits, m_minDT2Hits, m_minTrackerHits, DetId::Muon, MuonChamberResidual::ndof(), NULL, MuonChamberResidual::numHits(), phi, AlignmentMonitorBase::pNavigator(), funct::pow(), reco::TrackBase::pt(), ExpressReco_HICollisions_FallBack::pt, dttmaxenums::R, relativeConstraints::ring, mathSSE::sqrt(), relativeConstraints::station, ExpressReco_HICollisions_FallBack::track, MuonResidualsFromTrack::trackerNumHits(), MuonResidualsFromTrack::trackerRedChi2(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                                               {
   m_counter_event++;

  edm::ESHandle<GlobalTrackingGeometry> globalGeometry;
  iSetup.get<GlobalTrackingGeometryRecord>().get(globalGeometry);

  for (ConstTrajTrackPairCollection::const_iterator trajtrack = trajtracks.begin();  trajtrack != trajtracks.end();  ++trajtrack) {
    const Trajectory* traj = (*trajtrack).first;
    const reco::Track* track = (*trajtrack).second;

    m_counter_track++;

    if (m_minTrackPt < track->pt()  &&  track->pt() < m_maxTrackPt) {
      char charge = (track->charge() > 0 ? 1 : -1);
      // double qoverpt = track->charge() / track->pt();
      // double qoverpz = track->charge() / track->pz();
      MuonResidualsFromTrack muonResidualsFromTrack(globalGeometry, traj, pNavigator(), 1000.);

      m_counter_trackpt++;

      if (muonResidualsFromTrack.trackerNumHits() >= m_minTrackerHits  &&  muonResidualsFromTrack.trackerRedChi2() < m_maxTrackerRedChi2  &&  (m_allowTIDTEC  ||  !muonResidualsFromTrack.contains_TIDTEC())) {
    std::vector<DetId> chamberIds = muonResidualsFromTrack.chamberIds();

    m_counter_trackokay++;

    for (std::vector<DetId>::const_iterator chamberId = chamberIds.begin();  chamberId != chamberIds.end();  ++chamberId) {

      if (chamberId->det() == DetId::Muon  &&  chamberId->subdetId() == MuonSubdetId::DT) {
        MuonChamberResidual *dt13 = muonResidualsFromTrack.chamberResidual(*chamberId, MuonChamberResidual::kDT13);
        MuonChamberResidual *dt2 = muonResidualsFromTrack.chamberResidual(*chamberId, MuonChamberResidual::kDT2);
        DTChamberId id(chamberId->rawId());

        m_counter_dt++;

        if (dt13 != NULL  &&  dt13->numHits() >= m_minDT13Hits) {
           m_counter_13numhits++;

          double residual = dt13->global_residual();
          double resslope = dt13->global_resslope();
          double chi2 = dt13->chi2();
          int dof = dt13->ndof();

          GlobalPoint trackpos = dt13->global_trackpos();
          double phi = atan2(trackpos.y(), trackpos.x());
          double z = trackpos.z();

          assert(1 <= id.sector()  &&  id.sector() <= 14);

          if (id.station() == 1) m_DTvsz_station1[id.sector()-1]->fill_x(charge, z, residual, chi2, dof);
          if (id.station() == 2) m_DTvsz_station2[id.sector()-1]->fill_x(charge, z, residual, chi2, dof);
          if (id.station() == 3) m_DTvsz_station3[id.sector()-1]->fill_x(charge, z, residual, chi2, dof);
          if (id.station() == 4) m_DTvsz_station4[id.sector()-1]->fill_x(charge, z, residual, chi2, dof);

          if (id.station() == 1) m_DTvsz_station1[id.sector()-1]->fill_dxdz(charge, z, resslope, chi2, dof);
          if (id.station() == 2) m_DTvsz_station2[id.sector()-1]->fill_dxdz(charge, z, resslope, chi2, dof);
          if (id.station() == 3) m_DTvsz_station3[id.sector()-1]->fill_dxdz(charge, z, resslope, chi2, dof);
          if (id.station() == 4) m_DTvsz_station4[id.sector()-1]->fill_dxdz(charge, z, resslope, chi2, dof);

          if (id.station() == 1) m_DTvsphi_station1[id.wheel()+2]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 2) m_DTvsphi_station2[id.wheel()+2]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 3) m_DTvsphi_station3[id.wheel()+2]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 4) m_DTvsphi_station4[id.wheel()+2]->fill_x(charge, phi, residual, chi2, dof);

          if (id.station() == 1) m_DTvsphi_station1[id.wheel()+2]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 2) m_DTvsphi_station2[id.wheel()+2]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 3) m_DTvsphi_station3[id.wheel()+2]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 4) m_DTvsphi_station4[id.wheel()+2]->fill_dxdz(charge, phi, resslope, chi2, dof);
        }

        if (dt2 != NULL  &&  dt2->numHits() >= m_minDT2Hits) {
           m_counter_2numhits++;

          double residual = dt2->global_residual();
          double resslope = dt2->global_resslope();
          double chi2 = dt2->chi2();
          int dof = dt2->ndof();

          GlobalPoint trackpos = dt2->global_trackpos();
          double phi = atan2(trackpos.y(), trackpos.x());
          double z = trackpos.z();

          assert(1 <= id.sector()  &&  id.sector() <= 14);

          if (id.station() == 1) m_DTvsz_station1[id.sector()-1]->fill_y(charge, z, residual, chi2, dof);
          if (id.station() == 2) m_DTvsz_station2[id.sector()-1]->fill_y(charge, z, residual, chi2, dof);
          if (id.station() == 3) m_DTvsz_station3[id.sector()-1]->fill_y(charge, z, residual, chi2, dof);
          if (id.station() == 4) m_DTvsz_station4[id.sector()-1]->fill_y(charge, z, residual, chi2, dof);

          if (id.station() == 1) m_DTvsz_station1[id.sector()-1]->fill_dydz(charge, z, resslope, chi2, dof);
          if (id.station() == 2) m_DTvsz_station2[id.sector()-1]->fill_dydz(charge, z, resslope, chi2, dof);
          if (id.station() == 3) m_DTvsz_station3[id.sector()-1]->fill_dydz(charge, z, resslope, chi2, dof);
          if (id.station() == 4) m_DTvsz_station4[id.sector()-1]->fill_dydz(charge, z, resslope, chi2, dof);

          if (id.station() == 1) m_DTvsphi_station1[id.wheel()+2]->fill_y(charge, phi, residual, chi2, dof);
          if (id.station() == 2) m_DTvsphi_station2[id.wheel()+2]->fill_y(charge, phi, residual, chi2, dof);
          if (id.station() == 3) m_DTvsphi_station3[id.wheel()+2]->fill_y(charge, phi, residual, chi2, dof);
          if (id.station() == 4) m_DTvsphi_station4[id.wheel()+2]->fill_y(charge, phi, residual, chi2, dof);

          if (id.station() == 1) m_DTvsphi_station1[id.wheel()+2]->fill_dydz(charge, phi, resslope, chi2, dof);
          if (id.station() == 2) m_DTvsphi_station2[id.wheel()+2]->fill_dydz(charge, phi, resslope, chi2, dof);
          if (id.station() == 3) m_DTvsphi_station3[id.wheel()+2]->fill_dydz(charge, phi, resslope, chi2, dof);
          if (id.station() == 4) m_DTvsphi_station4[id.wheel()+2]->fill_dydz(charge, phi, resslope, chi2, dof);
        }
      }

      else if (chamberId->det() == DetId::Muon  &&  chamberId->subdetId() == MuonSubdetId::CSC) {
        MuonChamberResidual *csc = muonResidualsFromTrack.chamberResidual(*chamberId, MuonChamberResidual::kCSC);
        CSCDetId id(chamberId->rawId());

        m_counter_csc++;

        if (csc != NULL  &&  csc->numHits() >= m_minCSCHits) {
           m_counter_cscnumhits++;

          double residual = csc->global_residual();
          double resslope = csc->global_resslope();
          double chi2 = csc->chi2();
          int dof = csc->ndof();

          GlobalPoint trackpos = csc->global_trackpos();
          double phi = atan2(trackpos.y(), trackpos.x());
          double R = sqrt(pow(trackpos.x(), 2) + pow(trackpos.y(), 2));

          int chamber = id.chamber() - 1;
          if (id.station() > 1  &&  id.ring() == 1) chamber *= 2;

          assert(1 <= id.endcap()  &&  id.endcap() <= 2  &&  0 <= chamber  &&  chamber <= 35);

          if (id.station() == 1) m_CSCvsr_me1[id.endcap()-1][chamber]->fill_x(charge, R, residual, chi2, dof);
          if (id.station() == 2) m_CSCvsr_me2[id.endcap()-1][chamber]->fill_x(charge, R, residual, chi2, dof);
          if (id.station() == 3) m_CSCvsr_me3[id.endcap()-1][chamber]->fill_x(charge, R, residual, chi2, dof);
          if (id.station() == 4) m_CSCvsr_me4[id.endcap()-1][chamber]->fill_x(charge, R, residual, chi2, dof);

          if (id.station() == 1) m_CSCvsr_me1[id.endcap()-1][chamber]->fill_dxdz(charge, R, resslope, chi2, dof);
          if (id.station() == 2) m_CSCvsr_me2[id.endcap()-1][chamber]->fill_dxdz(charge, R, resslope, chi2, dof);
          if (id.station() == 3) m_CSCvsr_me3[id.endcap()-1][chamber]->fill_dxdz(charge, R, resslope, chi2, dof);
          if (id.station() == 4) m_CSCvsr_me4[id.endcap()-1][chamber]->fill_dxdz(charge, R, resslope, chi2, dof);

          if (id.station() == 1  &&  id.ring() == 1) m_CSCvsphi_me11[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 1  &&  id.ring() == 2) m_CSCvsphi_me12[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 1  &&  id.ring() == 3) m_CSCvsphi_me13[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 1  &&  id.ring() == 4) m_CSCvsphi_me14[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 2  &&  id.ring() == 1) m_CSCvsphi_me21[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 2  &&  id.ring() == 2) m_CSCvsphi_me22[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 3  &&  id.ring() == 1) m_CSCvsphi_me31[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 3  &&  id.ring() == 2) m_CSCvsphi_me32[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 4  &&  id.ring() == 1) m_CSCvsphi_me41[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);
          if (id.station() == 4  &&  id.ring() == 2) m_CSCvsphi_me42[id.endcap()-1]->fill_x(charge, phi, residual, chi2, dof);

          if (id.station() == 1  &&  id.ring() == 1) m_CSCvsphi_me11[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 1  &&  id.ring() == 2) m_CSCvsphi_me12[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 1  &&  id.ring() == 3) m_CSCvsphi_me13[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 1  &&  id.ring() == 4) m_CSCvsphi_me14[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 2  &&  id.ring() == 1) m_CSCvsphi_me21[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 2  &&  id.ring() == 2) m_CSCvsphi_me22[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 3  &&  id.ring() == 1) m_CSCvsphi_me31[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 3  &&  id.ring() == 2) m_CSCvsphi_me32[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 4  &&  id.ring() == 1) m_CSCvsphi_me41[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
          if (id.station() == 4  &&  id.ring() == 2) m_CSCvsphi_me42[id.endcap()-1]->fill_dxdz(charge, phi, resslope, chi2, dof);
        }
      }

      else { assert(false); }

    } // end loop over chambers
      } // end if track has enough tracker hits
    } // end if track has acceptable momentum
  } // end loop over tracks
}

std::string AlignmentMonitorMuonSystemMap1D::num02d ( int  num )

private

Definition at line 39 of file AlignmentMonitorMuonSystemMap1D.cc.

Referenced by book().

                                                         {
  int tens = num / 10;
  int ones = num % 10;
  
  std::string s_tens, s_ones;
  if (tens == 0) s_tens = std::string("0");
  if (tens == 1) s_tens = std::string("1");
  if (tens == 2) s_tens = std::string("2");
  if (tens == 3) s_tens = std::string("3");
  if (tens == 4) s_tens = std::string("4");
  if (tens == 5) s_tens = std::string("5");
  if (tens == 6) s_tens = std::string("6");
  if (tens == 7) s_tens = std::string("7");
  if (tens == 8) s_tens = std::string("8");
  if (tens == 9) s_tens = std::string("9");

  if (ones == 0) s_ones = std::string("0");
  if (ones == 1) s_ones = std::string("1");
  if (ones == 2) s_ones = std::string("2");
  if (ones == 3) s_ones = std::string("3");
  if (ones == 4) s_ones = std::string("4");
  if (ones == 5) s_ones = std::string("5");
  if (ones == 6) s_ones = std::string("6");
  if (ones == 7) s_ones = std::string("7");
  if (ones == 8) s_ones = std::string("8");
  if (ones == 9) s_ones = std::string("9");

  return s_tens + s_ones;
}

Friends And Related Function Documentation

friend class MuonSystemMapPlot1D

friend

Definition at line 38 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book().

Member Data Documentation

bool AlignmentMonitorMuonSystemMap1D::m_allowTIDTEC

private

Definition at line 46 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by event().

long AlignmentMonitorMuonSystemMap1D::m_counter_13numhits

private

Definition at line 83 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by afterAlignment(), book(), and event().

long AlignmentMonitorMuonSystemMap1D::m_counter_2numhits

private

Definition at line 84 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by afterAlignment(), book(), and event().

long AlignmentMonitorMuonSystemMap1D::m_counter_csc

private

Definition at line 85 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by afterAlignment(), book(), and event().

long AlignmentMonitorMuonSystemMap1D::m_counter_cscnumhits

private

Definition at line 86 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by afterAlignment(), book(), and event().

long AlignmentMonitorMuonSystemMap1D::m_counter_dt

private

Definition at line 82 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by afterAlignment(), book(), and event().

long AlignmentMonitorMuonSystemMap1D::m_counter_event

private

Definition at line 78 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by afterAlignment(), book(), and event().

long AlignmentMonitorMuonSystemMap1D::m_counter_track

private

Definition at line 79 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by afterAlignment(), book(), and event().

long AlignmentMonitorMuonSystemMap1D::m_counter_trackokay

private

Definition at line 81 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by afterAlignment(), book(), and event().

long AlignmentMonitorMuonSystemMap1D::m_counter_trackpt

private

Definition at line 80 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by afterAlignment(), book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me11[2]

private

Definition at line 65 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me12[2]

private

Definition at line 66 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me13[2]

private

Definition at line 67 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me14[2]

private

Definition at line 68 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me21[2]

private

Definition at line 69 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me22[2]

private

Definition at line 70 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me31[2]

private

Definition at line 71 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me32[2]

private

Definition at line 72 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me41[2]

private

Definition at line 73 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsphi_me42[2]

private

Definition at line 74 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsr_me1[2][36]

private

Definition at line 55 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsr_me2[2][36]

private

Definition at line 56 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsr_me3[2][36]

private

Definition at line 57 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_CSCvsr_me4[2][36]

private

Definition at line 58 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_DTvsphi_station1[5]

private

Definition at line 60 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_DTvsphi_station2[5]

private

Definition at line 61 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_DTvsphi_station3[5]

private

Definition at line 62 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_DTvsphi_station4[5]

private

Definition at line 63 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_DTvsz_station1[12]

private

Definition at line 51 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_DTvsz_station2[12]

private

Definition at line 52 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_DTvsz_station3[12]

private

Definition at line 53 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

MuonSystemMapPlot1D* AlignmentMonitorMuonSystemMap1D::m_DTvsz_station4[14]

private

Definition at line 54 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book(), and event().

double AlignmentMonitorMuonSystemMap1D::m_maxTrackerRedChi2

private

Definition at line 45 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by event().

double AlignmentMonitorMuonSystemMap1D::m_maxTrackPt

private

Definition at line 43 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by event().

int AlignmentMonitorMuonSystemMap1D::m_minCSCHits

private

Definition at line 49 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by event().

int AlignmentMonitorMuonSystemMap1D::m_minDT13Hits

private

Definition at line 47 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by event().

int AlignmentMonitorMuonSystemMap1D::m_minDT2Hits

private

Definition at line 48 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by event().

int AlignmentMonitorMuonSystemMap1D::m_minTrackerHits

private

Definition at line 44 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by event().

double AlignmentMonitorMuonSystemMap1D::m_minTrackPt

private

Definition at line 42 of file AlignmentMonitorMuonSystemMap1D.h.

std::vector<MuonSystemMapPlot1D*> AlignmentMonitorMuonSystemMap1D::m_plots

private

Definition at line 76 of file AlignmentMonitorMuonSystemMap1D.h.

Referenced by book().