|
|
|
Public Member Functions | |
| void | afterAlignment (const edm::EventSetup &iSetup) |
| AlignmentMonitorSegmentDifferences (const edm::ParameterSet &cfg) | |
| void | book () |
| Book or retrieve histograms; MUST be reimplemented. | |
| void | event (const edm::Event &iEvent, const edm::EventSetup &iSetup, const ConstTrajTrackPairCollection &iTrajTracks) |
| Called for each event (by "run()"): may be reimplemented. | |
| ~AlignmentMonitorSegmentDifferences () | |
Private Attributes | |
| bool | m_allowTIDTEC |
| TProfile * | m_cscinner_resid [2][18][3] |
| TProfile * | m_cscinner_slope [2][18][3] |
| TProfile * | m_cscouter_resid [2][36][2] |
| TProfile * | m_cscouter_slope [2][36][2] |
| TProfile * | m_dt13_resid [5][12][3] |
| TProfile * | m_dt13_slope [5][12][3] |
| TProfile * | m_dt2_resid [5][12][2] |
| TProfile * | m_dt2_slope [5][12][2] |
| double | m_maxTrackerRedChi2 |
| int | m_minCSCHits |
| int | m_minDT13Hits |
| int | m_minDT2Hits |
| int | m_minTrackerHits |
| double | m_minTrackPt |
| TH1F * | m_negcscinner_resid [2][18][3] |
| TH1F * | m_negcscinner_slope [2][18][3] |
| TH1F * | m_negcscouter_resid [2][36][2] |
| TH1F * | m_negcscouter_slope [2][36][2] |
| TH1F * | m_negdt13_resid [5][12][3] |
| TH1F * | m_negdt13_slope [5][12][3] |
| TH1F * | m_negdt2_resid [5][12][2] |
| TH1F * | m_negdt2_slope [5][12][2] |
| TH1F * | m_poscscinner_resid [2][18][3] |
| TH1F * | m_poscscinner_slope [2][18][3] |
| TH1F * | m_poscscouter_resid [2][36][2] |
| TH1F * | m_poscscouter_slope [2][36][2] |
| TH1F * | m_posdt13_resid [5][12][3] |
| TH1F * | m_posdt13_slope [5][12][3] |
| TH1F * | m_posdt2_resid [5][12][2] |
| TH1F * | m_posdt2_slope [5][12][2] |
Definition at line 35 of file AlignmentMonitorSegmentDifferences.cc.
| AlignmentMonitorSegmentDifferences::AlignmentMonitorSegmentDifferences | ( | const edm::ParameterSet & | cfg | ) |
Definition at line 94 of file AlignmentMonitorSegmentDifferences.cc.
: AlignmentMonitorBase(cfg, "AlignmentMonitorSegmentDifferences") , m_minTrackPt(cfg.getParameter<double>("minTrackPt")) , 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")) { }
| AlignmentMonitorSegmentDifferences::~AlignmentMonitorSegmentDifferences | ( | ) | [inline] |
Definition at line 38 of file AlignmentMonitorSegmentDifferences.cc.
{};
| void AlignmentMonitorSegmentDifferences::afterAlignment | ( | const edm::EventSetup & | iSetup | ) | [virtual] |
Called after updating AlignableTracker and AlignableMuon (by "endOfLoop()"): may be reimplemented
Reimplemented from AlignmentMonitorBase.
Definition at line 531 of file AlignmentMonitorSegmentDifferences.cc.
{
}
| void AlignmentMonitorSegmentDifferences::book | ( | ) | [virtual] |
Book or retrieve histograms; MUST be reimplemented.
Implements AlignmentMonitorBase.
Definition at line 106 of file AlignmentMonitorSegmentDifferences.cc.
References AlignmentMonitorBase::book1D(), AlignmentMonitorBase::bookProfile(), Reference_intrackfit_cff::endcap, m_cscinner_resid, m_cscinner_slope, m_cscouter_resid, m_cscouter_slope, m_dt13_resid, m_dt13_slope, m_dt2_resid, m_dt2_slope, m_minTrackPt, m_negcscinner_resid, m_negcscinner_slope, m_negcscouter_resid, m_negcscouter_slope, m_negdt13_resid, m_negdt13_slope, m_negdt2_resid, m_negdt2_slope, m_poscscinner_resid, m_poscscinner_slope, m_poscscouter_resid, m_poscscouter_slope, m_posdt13_resid, m_posdt13_slope, m_posdt2_resid, m_posdt2_slope, AlCaRecoCosmics_cfg::name, neg, and pos.
{
for (int wheel = -2; wheel <= +2; wheel++) {
for (int sector = 1; sector <= 12; sector++) {
char num[3];
num[0] = ('0' + (sector / 10));
num[1] = ('0' + (sector % 10));
num[2] = 0;
std::string wheelletter;
if (wheel == -2) wheelletter = "A";
else if (wheel == -1) wheelletter = "B";
else if (wheel == 0) wheelletter = "C";
else if (wheel == +1) wheelletter = "D";
else if (wheel == +2) wheelletter = "E";
std::string name, pos, neg;
name = (std::string("dt13_resid_") + wheelletter + std::string("_") + std::string(num) + std::string("_12"));
m_dt13_resid[wheel+2][sector-1][0] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_posdt13_resid[wheel+2][sector-1][0] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negdt13_resid[wheel+2][sector-1][0] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("dt13_resid_") + wheelletter + std::string("_") + std::string(num) + std::string("_23"));
m_dt13_resid[wheel+2][sector-1][1] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_posdt13_resid[wheel+2][sector-1][1] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negdt13_resid[wheel+2][sector-1][1] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("dt13_resid_") + wheelletter + std::string("_") + std::string(num) + std::string("_34"));
m_dt13_resid[wheel+2][sector-1][2] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_posdt13_resid[wheel+2][sector-1][2] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negdt13_resid[wheel+2][sector-1][2] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("dt2_resid_") + wheelletter + std::string("_") + std::string(num) + std::string("_12"));
m_dt2_resid[wheel+2][sector-1][0] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -200., 200., " ");
pos = std::string("pos") + name;
m_posdt2_resid[wheel+2][sector-1][0] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -20., 20.);
neg = std::string("neg") + name;
m_negdt2_resid[wheel+2][sector-1][0] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -20., 20.);
name = (std::string("dt2_resid_") + wheelletter + std::string("_") + std::string(num) + std::string("_23"));
m_dt2_resid[wheel+2][sector-1][1] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -200., 200., " ");
pos = std::string("pos") + name;
m_posdt2_resid[wheel+2][sector-1][1] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -20., 20.);
neg = std::string("neg") + name;
m_negdt2_resid[wheel+2][sector-1][1] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -20., 20.);
name = (std::string("dt13_slope_") + wheelletter + std::string("_") + std::string(num) + std::string("_12"));
m_dt13_slope[wheel+2][sector-1][0] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_posdt13_slope[wheel+2][sector-1][0] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negdt13_slope[wheel+2][sector-1][0] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("dt13_slope_") + wheelletter + std::string("_") + std::string(num) + std::string("_23"));
m_dt13_slope[wheel+2][sector-1][1] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_posdt13_slope[wheel+2][sector-1][1] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negdt13_slope[wheel+2][sector-1][1] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("dt13_slope_") + wheelletter + std::string("_") + std::string(num) + std::string("_34"));
m_dt13_slope[wheel+2][sector-1][2] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_posdt13_slope[wheel+2][sector-1][2] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negdt13_slope[wheel+2][sector-1][2] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("dt2_slope_") + wheelletter + std::string("_") + std::string(num) + std::string("_12"));
m_dt2_slope[wheel+2][sector-1][0] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -1000., 1000., " ");
pos = std::string("pos") + name;
m_posdt2_slope[wheel+2][sector-1][0] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -100., 100.);
neg = std::string("neg") + name;
m_negdt2_slope[wheel+2][sector-1][0] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -100., 100.);
name = (std::string("dt2_slope_") + wheelletter + std::string("_") + std::string(num) + std::string("_23"));
m_dt2_slope[wheel+2][sector-1][1] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -1000., 1000., " ");
pos = std::string("pos") + name;
m_posdt2_slope[wheel+2][sector-1][1] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -100., 100.);
neg = std::string("neg") + name;
m_negdt2_slope[wheel+2][sector-1][1] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -100., 100.);
}
}
for (int endcap = 1; endcap <= 2; endcap++) {
for (int chamber = 1; chamber <= 36; chamber++) {
char num[3];
num[0] = ('0' + (chamber / 10));
num[1] = ('0' + (chamber % 10));
num[2] = 0;
std::string endcapletter;
if (endcap == 1) endcapletter = "p";
else if (endcap == 2) endcapletter = "m";
std::string name, pos, neg;
name = (std::string("cscouter_resid_") + endcapletter + std::string("_") + std::string(num) + std::string("_12"));
m_cscouter_resid[endcap-1][chamber-1][0] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscouter_resid[endcap-1][chamber-1][0] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscouter_resid[endcap-1][chamber-1][0] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("cscouter_resid_") + endcapletter + std::string("_") + std::string(num) + std::string("_23"));
m_cscouter_resid[endcap-1][chamber-1][1] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscouter_resid[endcap-1][chamber-1][1] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscouter_resid[endcap-1][chamber-1][1] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("cscouter_slope_") + endcapletter + std::string("_") + std::string(num) + std::string("_12"));
m_cscouter_slope[endcap-1][chamber-1][0] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscouter_slope[endcap-1][chamber-1][0] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscouter_slope[endcap-1][chamber-1][0] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("cscouter_slope_") + endcapletter + std::string("_") + std::string(num) + std::string("_23"));
m_cscouter_slope[endcap-1][chamber-1][1] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscouter_slope[endcap-1][chamber-1][1] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscouter_slope[endcap-1][chamber-1][1] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
}
for (int chamber = 1; chamber <= 18; chamber++) {
char num[3];
num[0] = ('0' + (chamber / 10));
num[1] = ('0' + (chamber % 10));
num[2] = 0;
std::string endcapletter;
if (endcap == 1) endcapletter = "p";
else if (endcap == 2) endcapletter = "m";
std::string name, pos, neg;
name = (std::string("cscinner_resid_") + endcapletter + std::string("_") + std::string(num) + std::string("_12"));
m_cscinner_resid[endcap-1][chamber-1][0] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscinner_resid[endcap-1][chamber-1][0] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscinner_resid[endcap-1][chamber-1][0] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("cscinner_resid_") + endcapletter + std::string("_") + std::string(num) + std::string("_23"));
m_cscinner_resid[endcap-1][chamber-1][1] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscinner_resid[endcap-1][chamber-1][1] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscinner_resid[endcap-1][chamber-1][1] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("cscinner_resid_") + endcapletter + std::string("_") + std::string(num) + std::string("_34"));
m_cscinner_resid[endcap-1][chamber-1][2] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscinner_resid[endcap-1][chamber-1][2] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscinner_resid[endcap-1][chamber-1][2] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("cscinner_slope_") + endcapletter + std::string("_") + std::string(num) + std::string("_12"));
m_cscinner_slope[endcap-1][chamber-1][0] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscinner_slope[endcap-1][chamber-1][0] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscinner_slope[endcap-1][chamber-1][0] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("cscinner_slope_") + endcapletter + std::string("_") + std::string(num) + std::string("_23"));
m_cscinner_slope[endcap-1][chamber-1][1] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscinner_slope[endcap-1][chamber-1][1] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscinner_slope[endcap-1][chamber-1][1] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
name = (std::string("cscinner_slope_") + endcapletter + std::string("_") + std::string(num) + std::string("_34"));
m_cscinner_slope[endcap-1][chamber-1][2] = bookProfile("/iterN/", name.c_str(), name.c_str(), 20, -1./m_minTrackPt, 1./m_minTrackPt, 1, -100., 100., " ");
pos = std::string("pos") + name;
m_poscscinner_slope[endcap-1][chamber-1][2] = book1D("/iterN/", pos.c_str(), pos.c_str(), 100, -10., 10.);
neg = std::string("neg") + name;
m_negcscinner_slope[endcap-1][chamber-1][2] = book1D("/iterN/", neg.c_str(), neg.c_str(), 100, -10., 10.);
}
}
}
| void AlignmentMonitorSegmentDifferences::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 294 of file AlignmentMonitorSegmentDifferences.cc.
References MuonChamberResidual::chamberAlignable(), MuonResidualsFromTrack::chamberIds(), MuonResidualsFromTrack::chamberResidual(), reco::TrackBase::charge(), MuonResidualsFromTrack::contains_TIDTEC(), CSC(), GeomDetEnumerators::DT, edm::EventSetup::get(), MuonChamberResidual::global_residual(), MuonChamberResidual::global_resslope(), MuonChamberResidual::kCSC, MuonChamberResidual::kDT13, MuonChamberResidual::kDT2, m_allowTIDTEC, m_cscinner_resid, m_cscinner_slope, m_cscouter_resid, m_cscouter_slope, m_dt13_resid, m_dt13_slope, m_dt2_resid, m_dt2_slope, m_maxTrackerRedChi2, m_minCSCHits, m_minDT13Hits, m_minDT2Hits, m_minTrackerHits, m_minTrackPt, m_negcscinner_resid, m_negcscinner_slope, m_negcscouter_resid, m_negcscouter_slope, m_negdt13_resid, m_negdt13_slope, m_negdt2_resid, m_negdt2_slope, m_poscscinner_resid, m_poscscinner_slope, m_poscscouter_resid, m_poscscouter_slope, m_posdt13_resid, m_posdt13_slope, m_posdt2_resid, m_posdt2_slope, DetId::Muon, NULL, MuonChamberResidual::numHits(), AlignmentMonitorBase::pNavigator(), reco::TrackBase::pt(), ExpressReco_HICollisions_FallBack::track, MuonResidualsFromTrack::trackerNumHits(), and MuonResidualsFromTrack::trackerRedChi2().
{
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;
if (track->pt() > m_minTrackPt) {
double qoverpt = (track->charge() > 0 ? 1. : -1.) / track->pt();
MuonResidualsFromTrack muonResidualsFromTrack(globalGeometry, traj, pNavigator(), 1000.);
if (muonResidualsFromTrack.trackerNumHits() >= m_minTrackerHits && muonResidualsFromTrack.trackerRedChi2() < m_maxTrackerRedChi2 && (m_allowTIDTEC || !muonResidualsFromTrack.contains_TIDTEC())) {
std::vector<DetId> chamberIds = muonResidualsFromTrack.chamberIds();
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);
if (dt13 != NULL && dt13->numHits() >= m_minDT13Hits) {
DTChamberId thisid(chamberId->rawId());
for (std::vector<DetId>::const_iterator otherId = chamberIds.begin(); otherId != chamberIds.end(); ++otherId) {
if (otherId->det() == DetId::Muon && otherId->subdetId() == MuonSubdetId::DT) {
DTChamberId thatid(otherId->rawId());
if (thisid.rawId() != thatid.rawId() && thisid.wheel() == thatid.wheel() && thisid.sector() == thatid.sector()) {
MuonChamberResidual *dt13other = muonResidualsFromTrack.chamberResidual(*otherId, MuonChamberResidual::kDT13);
if (dt13other != NULL && dt13other->numHits() >= m_minDT13Hits) {
double slopediff = 1000. * (dt13->global_resslope() - dt13other->global_resslope());
double length = dt13->chamberAlignable()->surface().toGlobal(LocalPoint(0,0,0)).perp() - dt13other->chamberAlignable()->surface().toGlobal(LocalPoint(0,0,0)).perp();
double residdiff = 10. * (dt13->global_residual() + length*dt13->global_resslope() - dt13other->global_residual());
if (thisid.station() == 1 && thatid.station() == 2) {
m_dt13_resid[thisid.wheel()+2][thisid.sector()-1][0]->Fill(qoverpt, residdiff);
m_dt13_slope[thisid.wheel()+2][thisid.sector()-1][0]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_posdt13_resid[thisid.wheel()+2][thisid.sector()-1][0]->Fill(residdiff);
m_posdt13_slope[thisid.wheel()+2][thisid.sector()-1][0]->Fill(slopediff);
}
else {
m_negdt13_resid[thisid.wheel()+2][thisid.sector()-1][0]->Fill(residdiff);
m_negdt13_slope[thisid.wheel()+2][thisid.sector()-1][0]->Fill(slopediff);
}
}
else if (thisid.station() == 2 && thatid.station() == 3) {
m_dt13_resid[thisid.wheel()+2][thisid.sector()-1][1]->Fill(qoverpt, residdiff);
m_dt13_slope[thisid.wheel()+2][thisid.sector()-1][1]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_posdt13_resid[thisid.wheel()+2][thisid.sector()-1][1]->Fill(residdiff);
m_posdt13_slope[thisid.wheel()+2][thisid.sector()-1][1]->Fill(slopediff);
}
else {
m_negdt13_resid[thisid.wheel()+2][thisid.sector()-1][1]->Fill(residdiff);
m_negdt13_slope[thisid.wheel()+2][thisid.sector()-1][1]->Fill(slopediff);
}
}
else if (thisid.station() == 3 && thatid.station() == 4) {
m_dt13_resid[thisid.wheel()+2][thisid.sector()-1][2]->Fill(qoverpt, residdiff);
m_dt13_slope[thisid.wheel()+2][thisid.sector()-1][2]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_posdt13_resid[thisid.wheel()+2][thisid.sector()-1][2]->Fill(residdiff);
m_posdt13_slope[thisid.wheel()+2][thisid.sector()-1][2]->Fill(slopediff);
}
else {
m_negdt13_resid[thisid.wheel()+2][thisid.sector()-1][2]->Fill(residdiff);
m_negdt13_slope[thisid.wheel()+2][thisid.sector()-1][2]->Fill(slopediff);
}
}
} // end other numhits
} // end this near other
} // end other is DT
} // end loop over other
} // end if DT13
if (dt2 != NULL && dt2->numHits() >= m_minDT2Hits) {
DTChamberId thisid(chamberId->rawId());
for (std::vector<DetId>::const_iterator otherId = chamberIds.begin(); otherId != chamberIds.end(); ++otherId) {
if (otherId->det() == DetId::Muon && otherId->subdetId() == MuonSubdetId::DT) {
DTChamberId thatid(otherId->rawId());
if (thisid.rawId() != thatid.rawId() && thisid.wheel() == thatid.wheel() && thisid.sector() == thatid.sector()) {
MuonChamberResidual *dt2other = muonResidualsFromTrack.chamberResidual(*otherId, MuonChamberResidual::kDT2);
if (dt2other != NULL && dt2other->numHits() >= m_minDT2Hits) {
double slopediff = 1000. * (dt2->global_resslope() - dt2other->global_resslope());
double length = dt2->chamberAlignable()->surface().toGlobal(LocalPoint(0,0,0)).perp() - dt2other->chamberAlignable()->surface().toGlobal(LocalPoint(0,0,0)).perp();
double residdiff = 10. * (dt2->global_residual() + length*dt2->global_resslope() - dt2other->global_residual());
if (thisid.station() == 1 && thatid.station() == 2) {
m_dt2_resid[thisid.wheel()+2][thisid.sector()-1][0]->Fill(qoverpt, residdiff);
m_dt2_slope[thisid.wheel()+2][thisid.sector()-1][0]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_posdt2_resid[thisid.wheel()+2][thisid.sector()-1][0]->Fill(residdiff);
m_posdt2_slope[thisid.wheel()+2][thisid.sector()-1][0]->Fill(slopediff);
}
else {
m_negdt2_resid[thisid.wheel()+2][thisid.sector()-1][0]->Fill(residdiff);
m_negdt2_slope[thisid.wheel()+2][thisid.sector()-1][0]->Fill(slopediff);
}
}
else if (thisid.station() == 2 && thatid.station() == 3) {
m_dt2_resid[thisid.wheel()+2][thisid.sector()-1][1]->Fill(qoverpt, residdiff);
m_dt2_slope[thisid.wheel()+2][thisid.sector()-1][1]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_posdt2_resid[thisid.wheel()+2][thisid.sector()-1][1]->Fill(residdiff);
m_posdt2_slope[thisid.wheel()+2][thisid.sector()-1][1]->Fill(slopediff);
}
else {
m_negdt2_resid[thisid.wheel()+2][thisid.sector()-1][1]->Fill(residdiff);
m_negdt2_slope[thisid.wheel()+2][thisid.sector()-1][1]->Fill(slopediff);
}
}
else if (thisid.station() == 3 && thatid.station() == 4) {
m_dt2_resid[thisid.wheel()+2][thisid.sector()-1][2]->Fill(qoverpt, residdiff);
m_dt2_slope[thisid.wheel()+2][thisid.sector()-1][2]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_posdt2_resid[thisid.wheel()+2][thisid.sector()-1][2]->Fill(residdiff);
m_posdt2_slope[thisid.wheel()+2][thisid.sector()-1][2]->Fill(slopediff);
}
else {
m_negdt2_resid[thisid.wheel()+2][thisid.sector()-1][2]->Fill(residdiff);
m_negdt2_slope[thisid.wheel()+2][thisid.sector()-1][2]->Fill(slopediff);
}
}
} // end other numhits
} // end this near other
} // end other is DT
} // end loop over other
} // end if DT2
} // end if DT
else if (chamberId->det() == DetId::Muon && chamberId->subdetId() == MuonSubdetId::CSC) {
MuonChamberResidual *csc = muonResidualsFromTrack.chamberResidual(*chamberId, MuonChamberResidual::kCSC);
if (csc->numHits() >= m_minCSCHits) {
CSCDetId thisid(chamberId->rawId());
for (std::vector<DetId>::const_iterator otherId = chamberIds.begin(); otherId != chamberIds.end(); ++otherId) {
if (otherId->det() == DetId::Muon && otherId->subdetId() == MuonSubdetId::CSC) {
CSCDetId thatid(otherId->rawId());
if (thisid.rawId() != thatid.rawId() && thisid.endcap() == thatid.endcap()) {
MuonChamberResidual *cscother = muonResidualsFromTrack.chamberResidual(*otherId, MuonChamberResidual::kCSC);
if (cscother != NULL && cscother->numHits() >= m_minCSCHits) {
double slopediff = 1000. * (csc->global_resslope() - cscother->global_resslope());
double length = csc->chamberAlignable()->surface().toGlobal(LocalPoint(0,0,0)).z() - cscother->chamberAlignable()->surface().toGlobal(LocalPoint(0,0,0)).z();
double residdiff = 10. * (csc->global_residual() + length*csc->global_resslope() - cscother->global_residual());
int thischamber = thisid.chamber();
int thisring = thisid.ring();
if (thisid.station() == 1 && (thisring == 1 || thisring == 4)) {
thischamber = (thischamber - 1) / 2 + 1;
thisring = 1;
}
if (thisring == thatid.ring() && thischamber == thatid.chamber()) {
if (thisring == 2 && thisid.station() == 1 && thatid.station() == 2) {
m_cscouter_resid[thisid.endcap()-1][thischamber-1][0]->Fill(qoverpt, residdiff);
m_cscouter_slope[thisid.endcap()-1][thischamber-1][0]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_poscscouter_resid[thisid.endcap()-1][thischamber-1][0]->Fill(residdiff);
m_poscscouter_slope[thisid.endcap()-1][thischamber-1][0]->Fill(slopediff);
}
else {
m_negcscouter_resid[thisid.endcap()-1][thischamber-1][0]->Fill(residdiff);
m_negcscouter_slope[thisid.endcap()-1][thischamber-1][0]->Fill(slopediff);
}
}
else if (thisring == 2 && thisid.station() == 2 && thatid.station() == 3) {
m_cscouter_resid[thisid.endcap()-1][thischamber-1][1]->Fill(qoverpt, residdiff);
m_cscouter_slope[thisid.endcap()-1][thischamber-1][1]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_poscscouter_resid[thisid.endcap()-1][thischamber-1][1]->Fill(residdiff);
m_poscscouter_slope[thisid.endcap()-1][thischamber-1][1]->Fill(slopediff);
}
else {
m_negcscouter_resid[thisid.endcap()-1][thischamber-1][1]->Fill(residdiff);
m_negcscouter_slope[thisid.endcap()-1][thischamber-1][1]->Fill(slopediff);
}
}
else if (thisring == 1 && thisid.station() == 1 && thatid.station() == 2) {
m_cscinner_resid[thisid.endcap()-1][thischamber-1][0]->Fill(qoverpt, residdiff);
m_cscinner_slope[thisid.endcap()-1][thischamber-1][0]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_poscscinner_resid[thisid.endcap()-1][thischamber-1][0]->Fill(residdiff);
m_poscscinner_slope[thisid.endcap()-1][thischamber-1][0]->Fill(slopediff);
}
else {
m_negcscinner_resid[thisid.endcap()-1][thischamber-1][0]->Fill(residdiff);
m_negcscinner_slope[thisid.endcap()-1][thischamber-1][0]->Fill(slopediff);
}
}
else if (thisring == 1 && thisid.station() == 2 && thatid.station() == 3) {
m_cscinner_resid[thisid.endcap()-1][thischamber-1][1]->Fill(qoverpt, residdiff);
m_cscinner_slope[thisid.endcap()-1][thischamber-1][1]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_poscscinner_resid[thisid.endcap()-1][thischamber-1][1]->Fill(residdiff);
m_poscscinner_slope[thisid.endcap()-1][thischamber-1][1]->Fill(slopediff);
}
else {
m_negcscinner_resid[thisid.endcap()-1][thischamber-1][1]->Fill(residdiff);
m_negcscinner_slope[thisid.endcap()-1][thischamber-1][1]->Fill(slopediff);
}
}
else if (thisring == 1 && thisid.station() == 3 && thatid.station() == 4) {
m_cscinner_resid[thisid.endcap()-1][thischamber-1][2]->Fill(qoverpt, residdiff);
m_cscinner_slope[thisid.endcap()-1][thischamber-1][2]->Fill(qoverpt, slopediff);
if (qoverpt > 0) {
m_poscscinner_resid[thisid.endcap()-1][thischamber-1][2]->Fill(residdiff);
m_poscscinner_slope[thisid.endcap()-1][thischamber-1][2]->Fill(slopediff);
}
else {
m_negcscinner_resid[thisid.endcap()-1][thischamber-1][2]->Fill(residdiff);
m_negcscinner_slope[thisid.endcap()-1][thischamber-1][2]->Fill(slopediff);
}
}
}
} // end other numhits
} // end this near other
} // end other is DT
} // end loop over other
} // end if csc
} // end if CSC
} // end loop over chamberIds
} // end if refit is okay
} // end if track pT is within range
} // end loop over tracks
}
bool AlignmentMonitorSegmentDifferences::m_allowTIDTEC [private] |
Definition at line 48 of file AlignmentMonitorSegmentDifferences.cc.
Referenced by event().
TProfile* AlignmentMonitorSegmentDifferences::m_cscinner_resid[2][18][3] [private] |
Definition at line 70 of file AlignmentMonitorSegmentDifferences.cc.
TProfile* AlignmentMonitorSegmentDifferences::m_cscinner_slope[2][18][3] [private] |
Definition at line 71 of file AlignmentMonitorSegmentDifferences.cc.
TProfile* AlignmentMonitorSegmentDifferences::m_cscouter_resid[2][36][2] [private] |
Definition at line 68 of file AlignmentMonitorSegmentDifferences.cc.
TProfile* AlignmentMonitorSegmentDifferences::m_cscouter_slope[2][36][2] [private] |
Definition at line 69 of file AlignmentMonitorSegmentDifferences.cc.
TProfile* AlignmentMonitorSegmentDifferences::m_dt13_resid[5][12][3] [private] |
Definition at line 54 of file AlignmentMonitorSegmentDifferences.cc.
TProfile* AlignmentMonitorSegmentDifferences::m_dt13_slope[5][12][3] [private] |
Definition at line 55 of file AlignmentMonitorSegmentDifferences.cc.
TProfile* AlignmentMonitorSegmentDifferences::m_dt2_resid[5][12][2] [private] |
Definition at line 56 of file AlignmentMonitorSegmentDifferences.cc.
TProfile* AlignmentMonitorSegmentDifferences::m_dt2_slope[5][12][2] [private] |
Definition at line 57 of file AlignmentMonitorSegmentDifferences.cc.
double AlignmentMonitorSegmentDifferences::m_maxTrackerRedChi2 [private] |
Definition at line 47 of file AlignmentMonitorSegmentDifferences.cc.
Referenced by event().
int AlignmentMonitorSegmentDifferences::m_minCSCHits [private] |
Definition at line 51 of file AlignmentMonitorSegmentDifferences.cc.
Referenced by event().
int AlignmentMonitorSegmentDifferences::m_minDT13Hits [private] |
Definition at line 49 of file AlignmentMonitorSegmentDifferences.cc.
Referenced by event().
int AlignmentMonitorSegmentDifferences::m_minDT2Hits [private] |
Definition at line 50 of file AlignmentMonitorSegmentDifferences.cc.
Referenced by event().
int AlignmentMonitorSegmentDifferences::m_minTrackerHits [private] |
Definition at line 46 of file AlignmentMonitorSegmentDifferences.cc.
Referenced by event().
double AlignmentMonitorSegmentDifferences::m_minTrackPt [private] |
Definition at line 45 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_negcscinner_resid[2][18][3] [private] |
Definition at line 78 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_negcscinner_slope[2][18][3] [private] |
Definition at line 79 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_negcscouter_resid[2][36][2] [private] |
Definition at line 76 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_negcscouter_slope[2][36][2] [private] |
Definition at line 77 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_negdt13_resid[5][12][3] [private] |
Definition at line 62 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_negdt13_slope[5][12][3] [private] |
Definition at line 63 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_negdt2_resid[5][12][2] [private] |
Definition at line 64 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_negdt2_slope[5][12][2] [private] |
Definition at line 65 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_poscscinner_resid[2][18][3] [private] |
Definition at line 74 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_poscscinner_slope[2][18][3] [private] |
Definition at line 75 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_poscscouter_resid[2][36][2] [private] |
Definition at line 72 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_poscscouter_slope[2][36][2] [private] |
Definition at line 73 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_posdt13_resid[5][12][3] [private] |
Definition at line 58 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_posdt13_slope[5][12][3] [private] |
Definition at line 59 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_posdt2_resid[5][12][2] [private] |
Definition at line 60 of file AlignmentMonitorSegmentDifferences.cc.
TH1F* AlignmentMonitorSegmentDifferences::m_posdt2_slope[5][12][2] [private] |
Definition at line 61 of file AlignmentMonitorSegmentDifferences.cc.
1.7.3