|
|
|
Public Member Functions | |
| void | afterAlignment (const edm::EventSetup &iSetup) |
| AlignmentMonitorMuonVsCurvature (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. | |
| ~AlignmentMonitorMuonVsCurvature () | |
Private Types | |
| enum | { kWheelMinus2 = 0, kWheelMinus1, kWheelZero, kWheelPlus1, kWheelPlus2, kWheelMEm11, kWheelMEm12, kWheelMEm13, kWheelMEm14, kWheelMEp11, kWheelMEp12, kWheelMEp13, kWheelMEp14, kNumWheels } |
| enum | { kEndcapMEm11 = 0, kEndcapMEm12, kEndcapMEm13, kEndcapMEm14, kEndcapMEp11, kEndcapMEp12, kEndcapMEp13, kEndcapMEp14, kNumEndcap } |
| enum | { kDeltaX = 0, kDeltaDxDz, kPtErr, kCurvErr, kNumComponents } |
Private Attributes | |
| bool | m_allowTIDTEC |
| int | m_layer |
| double | m_maxDxy |
| double | m_maxTrackerRedChi2 |
| int | m_minCSCHits |
| int | m_minDT13Hits |
| int | m_minDT2Hits |
| int | m_minTrackerHits |
| double | m_minTrackPt |
| std::string | m_propagator |
| TH1F * | th1f_trackerRedChi2 |
| TH1F * | th1f_trackerRedChi2Diff |
| TH2F * | th2f_endcap [kNumEndcap][36][kNumComponents] |
| TH2F * | th2f_evens [kNumEndcap][kNumComponents] |
| TH2F * | th2f_odds [kNumEndcap][kNumComponents] |
| TH2F * | th2f_wheelsector [kNumWheels][12][kNumComponents] |
| TProfile * | tprofile_endcap [kNumEndcap][36][kNumComponents] |
| TProfile * | tprofile_evens [kNumEndcap][kNumComponents] |
| TProfile * | tprofile_odds [kNumEndcap][kNumComponents] |
| TProfile * | tprofile_wheelsector [kNumWheels][12][kNumComponents] |
Definition at line 43 of file AlignmentMonitorMuonVsCurvature.cc.
anonymous enum [private] |
| kWheelMinus2 | |
| kWheelMinus1 | |
| kWheelZero | |
| kWheelPlus1 | |
| kWheelPlus2 | |
| kWheelMEm11 | |
| kWheelMEm12 | |
| kWheelMEm13 | |
| kWheelMEm14 | |
| kWheelMEp11 | |
| kWheelMEp12 | |
| kWheelMEp13 | |
| kWheelMEp14 | |
| kNumWheels |
Definition at line 64 of file AlignmentMonitorMuonVsCurvature.cc.
anonymous enum [private] |
| kEndcapMEm11 | |
| kEndcapMEm12 | |
| kEndcapMEm13 | |
| kEndcapMEm14 | |
| kEndcapMEp11 | |
| kEndcapMEp12 | |
| kEndcapMEp13 | |
| kEndcapMEp14 | |
| kNumEndcap |
Definition at line 81 of file AlignmentMonitorMuonVsCurvature.cc.
anonymous enum [private] |
Definition at line 93 of file AlignmentMonitorMuonVsCurvature.cc.
{
kDeltaX = 0,
kDeltaDxDz,
kPtErr,
kCurvErr,
kNumComponents
};
| AlignmentMonitorMuonVsCurvature::AlignmentMonitorMuonVsCurvature | ( | const edm::ParameterSet & | cfg | ) |
Definition at line 128 of file AlignmentMonitorMuonVsCurvature.cc.
: AlignmentMonitorBase(cfg, "AlignmentMonitorMuonVsCurvature") , 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_maxDxy(cfg.getParameter<double>("maxDxy")) , m_minDT13Hits(cfg.getParameter<int>("minDT13Hits")) , m_minDT2Hits(cfg.getParameter<int>("minDT2Hits")) , m_minCSCHits(cfg.getParameter<int>("minCSCHits")) , m_layer(cfg.getParameter<int>("layer")) , m_propagator(cfg.getParameter<std::string>("propagator")) {}
| AlignmentMonitorMuonVsCurvature::~AlignmentMonitorMuonVsCurvature | ( | ) | [inline] |
Definition at line 46 of file AlignmentMonitorMuonVsCurvature.cc.
{};
| void AlignmentMonitorMuonVsCurvature::afterAlignment | ( | const edm::EventSetup & | iSetup | ) | [virtual] |
Called after updating AlignableTracker and AlignableMuon (by "endOfLoop()"): may be reimplemented
Reimplemented from AlignmentMonitorBase.
Definition at line 561 of file AlignmentMonitorMuonVsCurvature.cc.
{}
| void AlignmentMonitorMuonVsCurvature::book | ( | ) | [virtual] |
Book or retrieve histograms; MUST be reimplemented.
Implements AlignmentMonitorBase.
Definition at line 142 of file AlignmentMonitorMuonVsCurvature.cc.
References AlignmentMonitorBase::book1D(), AlignmentMonitorBase::book2D(), AlignmentMonitorBase::bookProfile(), Reference_intrackfit_cff::endcap, kCurvErr, kDeltaDxDz, kDeltaX, kEndcapMEm11, kEndcapMEm12, kEndcapMEm13, kEndcapMEm14, kEndcapMEp11, kEndcapMEp12, kEndcapMEp13, kEndcapMEp14, kNumComponents, kNumEndcap, kNumWheels, kPtErr, kWheelMEm11, kWheelMEm12, kWheelMEm13, kWheelMEm14, kWheelMEp11, kWheelMEp12, kWheelMEp13, kWheelMEp14, kWheelMinus1, kWheelMinus2, kWheelPlus1, kWheelPlus2, kWheelZero, th1f_trackerRedChi2, th1f_trackerRedChi2Diff, th2f_endcap, th2f_evens, th2f_odds, th2f_wheelsector, tprofile_endcap, tprofile_evens, tprofile_odds, and tprofile_wheelsector.
{
for (int wheel = 0; wheel < kNumWheels; wheel++) {
std::stringstream wheelname;
if (wheel == kWheelMinus2) wheelname << "wheelm2_";
else if (wheel == kWheelMinus1) wheelname << "wheelm1_";
else if (wheel == kWheelZero) wheelname << "wheelz_";
else if (wheel == kWheelPlus1) wheelname << "wheelp1_";
else if (wheel == kWheelPlus2) wheelname << "wheelp2_";
else if (wheel == kWheelMEm11) wheelname << "wheelmem11_";
else if (wheel == kWheelMEm12) wheelname << "wheelmem12_";
else if (wheel == kWheelMEm13) wheelname << "wheelmem13_";
else if (wheel == kWheelMEm14) wheelname << "wheelmem14_";
else if (wheel == kWheelMEp11) wheelname << "wheelmep11_";
else if (wheel == kWheelMEp12) wheelname << "wheelmep12_";
else if (wheel == kWheelMEp13) wheelname << "wheelmep13_";
else if (wheel == kWheelMEp14) wheelname << "wheelmep14_";
for (int sector = 0; sector < 12; sector++) {
std::stringstream sectorname;
if (sector == 0) sectorname << "sector01_";
else if (sector == 1) sectorname << "sector02_";
else if (sector == 2) sectorname << "sector03_";
else if (sector == 3) sectorname << "sector04_";
else if (sector == 4) sectorname << "sector05_";
else if (sector == 5) sectorname << "sector06_";
else if (sector == 6) sectorname << "sector07_";
else if (sector == 7) sectorname << "sector08_";
else if (sector == 8) sectorname << "sector09_";
else if (sector == 9) sectorname << "sector10_";
else if (sector == 10) sectorname << "sector11_";
else if (sector == 11) sectorname << "sector12_";
for (int component = 0; component < kNumComponents; component++) {
std::stringstream th2f_name, tprofile_name;
th2f_name << "th2f_" << wheelname.str() << sectorname.str();
tprofile_name << "tprofile_" << wheelname.str() << sectorname.str();
double minmax = 15.;
if (component == kDeltaX) {
th2f_name << "deltax";
tprofile_name << "deltax";
minmax = 15.;
}
else if (component == kDeltaDxDz) {
th2f_name << "deltadxdz";
tprofile_name << "deltadxdz";
minmax = 15.;
}
else if (component == kPtErr) {
th2f_name << "pterr";
tprofile_name << "pterr";
minmax = 15.;
}
else if (component == kCurvErr) {
th2f_name << "curverr";
tprofile_name << "curverr";
minmax = 0.0015;
}
if (component == kPtErr) {
th2f_wheelsector[wheel][sector][component] = book2D("/iterN/", th2f_name.str().c_str(), "", 25, -200., 200., 30, -minmax, minmax);
tprofile_wheelsector[wheel][sector][component] = bookProfile("/iterN/", tprofile_name.str().c_str(), "", 25, -200., 200.);
}
else {
th2f_wheelsector[wheel][sector][component] = book2D("/iterN/", th2f_name.str().c_str(), "", 25, -0.05, 0.05, 30, -minmax, minmax);
tprofile_wheelsector[wheel][sector][component] = bookProfile("/iterN/", tprofile_name.str().c_str(), "", 25, -0.05, 0.05);
}
}
}
}
for (int endcap = 0; endcap < kNumEndcap; endcap++) {
std::stringstream endcapname;
if (endcap == kEndcapMEm11) endcapname << "endcapmem11_";
else if (endcap == kEndcapMEm12) endcapname << "endcapmem12_";
else if (endcap == kEndcapMEm13) endcapname << "endcapmem13_";
else if (endcap == kEndcapMEm14) endcapname << "endcapmem14_";
else if (endcap == kEndcapMEp11) endcapname << "endcapmep11_";
else if (endcap == kEndcapMEp12) endcapname << "endcapmep12_";
else if (endcap == kEndcapMEp13) endcapname << "endcapmep13_";
else if (endcap == kEndcapMEp14) endcapname << "endcapmep14_";
for (int component = 0; component < kNumComponents; component++) {
std::stringstream componentname;
double minmax = 15.;
if (component == kDeltaX) {
componentname << "deltax";
minmax = 15.;
}
else if (component == kDeltaDxDz) {
componentname << "deltadxdz";
minmax = 15.;
}
else if (component == kPtErr) {
componentname << "pterr";
minmax = 15.;
}
else if (component == kCurvErr) {
componentname << "curverr";
minmax = 0.0015;
}
std::stringstream th2f_evens_name, th2f_odds_name, tprofile_evens_name, tprofile_odds_name;
th2f_evens_name << "th2f_" << endcapname.str() << "evens_" << componentname.str();
th2f_odds_name << "th2f_" << endcapname.str() << "odds_" << componentname.str();
tprofile_evens_name << "tprofile_" << endcapname.str() << "evens_" << componentname.str();
tprofile_odds_name << "tprofile_" << endcapname.str() << "odds_" << componentname.str();
if (component == kPtErr) {
th2f_evens[endcap][component] = book2D("/iterN/", th2f_evens_name.str().c_str(), "", 25, -200., 200., 30, -minmax, minmax);
th2f_odds[endcap][component] = book2D("/iterN/", th2f_odds_name.str().c_str(), "", 25, -200., 200., 30, -minmax, minmax);
tprofile_evens[endcap][component] = bookProfile("/iterN/", tprofile_evens_name.str().c_str(), "", 25, -200., 200.);
tprofile_odds[endcap][component] = bookProfile("/iterN/", tprofile_odds_name.str().c_str(), "", 25, -200., 200.);
}
else {
th2f_evens[endcap][component] = book2D("/iterN/", th2f_evens_name.str().c_str(), "", 25, -0.05, 0.05, 30, -minmax, minmax);
th2f_odds[endcap][component] = book2D("/iterN/", th2f_odds_name.str().c_str(), "", 25, -0.05, 0.05, 30, -minmax, minmax);
tprofile_evens[endcap][component] = bookProfile("/iterN/", tprofile_evens_name.str().c_str(), "", 25, -0.05, 0.05);
tprofile_odds[endcap][component] = bookProfile("/iterN/", tprofile_odds_name.str().c_str(), "", 25, -0.05, 0.05);
}
for (int chamber = 0; chamber < 36; chamber++) {
std::stringstream chambername;
if (chamber == 0) chambername << "chamber01_";
else if (chamber == 1) chambername << "chamber02_";
else if (chamber == 2) chambername << "chamber03_";
else if (chamber == 3) chambername << "chamber04_";
else if (chamber == 4) chambername << "chamber05_";
else if (chamber == 5) chambername << "chamber06_";
else if (chamber == 6) chambername << "chamber07_";
else if (chamber == 7) chambername << "chamber08_";
else if (chamber == 8) chambername << "chamber09_";
else if (chamber == 9) chambername << "chamber10_";
else if (chamber == 10) chambername << "chamber11_";
else if (chamber == 11) chambername << "chamber12_";
else if (chamber == 12) chambername << "chamber13_";
else if (chamber == 13) chambername << "chamber14_";
else if (chamber == 14) chambername << "chamber15_";
else if (chamber == 15) chambername << "chamber16_";
else if (chamber == 16) chambername << "chamber17_";
else if (chamber == 17) chambername << "chamber18_";
else if (chamber == 18) chambername << "chamber19_";
else if (chamber == 19) chambername << "chamber20_";
else if (chamber == 20) chambername << "chamber21_";
else if (chamber == 21) chambername << "chamber22_";
else if (chamber == 22) chambername << "chamber23_";
else if (chamber == 23) chambername << "chamber24_";
else if (chamber == 24) chambername << "chamber25_";
else if (chamber == 25) chambername << "chamber26_";
else if (chamber == 26) chambername << "chamber27_";
else if (chamber == 27) chambername << "chamber28_";
else if (chamber == 28) chambername << "chamber29_";
else if (chamber == 29) chambername << "chamber30_";
else if (chamber == 30) chambername << "chamber31_";
else if (chamber == 31) chambername << "chamber32_";
else if (chamber == 32) chambername << "chamber33_";
else if (chamber == 33) chambername << "chamber34_";
else if (chamber == 34) chambername << "chamber35_";
else if (chamber == 35) chambername << "chamber36_";
std::stringstream th2f_name, tprofile_name;
th2f_name << "th2f_" << endcapname.str() << chambername.str() << componentname.str();
tprofile_name << "tprofile_" << endcapname.str() << chambername.str() << componentname.str();
if (component == kPtErr) {
th2f_endcap[endcap][chamber][component] = book2D("/iterN/", th2f_name.str().c_str(), "", 25, -200., 200., 30, -minmax, minmax);
tprofile_endcap[endcap][chamber][component] = bookProfile("/iterN/", tprofile_name.str().c_str(), "", 25, -200., 200.);
}
else {
th2f_endcap[endcap][chamber][component] = book2D("/iterN/", th2f_name.str().c_str(), "", 25, -0.05, 0.05, 30, -minmax, minmax);
tprofile_endcap[endcap][chamber][component] = bookProfile("/iterN/", tprofile_name.str().c_str(), "", 25, -0.05, 0.05);
}
}
}
}
th1f_trackerRedChi2 = book1D("/iterN/", "trackerRedChi2", "Refit tracker reduced chi^2", 100, 0., 30.);
th1f_trackerRedChi2Diff = book1D("/iterN/", "trackerRedChi2Diff", "Fit-minus-refit tracker reduced chi^2", 100, -5., 5.);
}
| void AlignmentMonitorMuonVsCurvature::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 323 of file AlignmentMonitorMuonVsCurvature.cc.
References MuonResidualsFromTrack::chamberIds(), MuonResidualsFromTrack::chamberResidual(), TrajectoryStateOnSurface::charge(), reco::TrackBase::charge(), MuonResidualsFromTrack::contains_TIDTEC(), CSC(), GeomDetEnumerators::DT, reco::TrackBase::dxy(), Reference_intrackfit_cff::endcap, edm::EventSetup::get(), MuonChamberResidual::global_hitresid(), MuonChamberResidual::global_resslope(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), MuonChamberResidual::hitresid(), TrajectoryStateOnSurface::isValid(), MuonChamberResidual::kCSC, kCurvErr, kDeltaDxDz, kDeltaX, MuonChamberResidual::kDT13, kEndcapMEm11, kEndcapMEm12, kEndcapMEm13, kEndcapMEm14, kEndcapMEp11, kEndcapMEp12, kEndcapMEp13, kEndcapMEp14, kPtErr, kWheelMEm11, kWheelMEm12, kWheelMEm13, kWheelMEm14, kWheelMEp11, kWheelMEp12, kWheelMEp13, kWheelMEp14, kWheelMinus1, kWheelMinus2, kWheelPlus1, kWheelPlus2, kWheelZero, MuonChamberResidual::localid(), TrajectoryStateOnSurface::localPosition(), m_allowTIDTEC, m_layer, m_maxDxy, m_maxTrackerRedChi2, m_minCSCHits, m_minDT13Hits, m_minTrackerHits, m_minTrackPt, m_propagator, PV3DBase< T, PVType, FrameType >::mag(), Trajectory::measurements(), DetId::Muon, reco::TrackBase::normalizedChi2(), NULL, MuonChamberResidual::numHits(), PV3DBase< T, PVType, FrameType >::perp(), AlignmentMonitorBase::pNavigator(), pos, LargeD0_PixelPairStep_cff::propagator, reco::TrackBase::pt(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), MuonChamberResidual::resslope(), mathSSE::sqrt(), th1f_trackerRedChi2, th1f_trackerRedChi2Diff, th2f_endcap, th2f_evens, th2f_odds, th2f_wheelsector, tprofile_endcap, tprofile_evens, tprofile_odds, tprofile_wheelsector, MuonResidualsFromTrack::trackerNumHits(), MuonResidualsFromTrack::trackerRedChi2(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
edm::ESHandle<GlobalTrackingGeometry> globalGeometry;
iSetup.get<GlobalTrackingGeometryRecord>().get(globalGeometry);
edm::ESHandle<Propagator> propagator;
iSetup.get<TrackingComponentsRecord>().get(m_propagator, propagator);
edm::ESHandle<MagneticField> magneticField;
iSetup.get<IdealMagneticFieldRecord>().get(magneticField);
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 && fabs(track->dxy()) < m_maxDxy) {
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();
double qoverpt = track->charge() / track->pt();
double px = track->px();
double py = track->py();
double pz = track->pz();
th1f_trackerRedChi2->Fill(muonResidualsFromTrack.trackerRedChi2());
th1f_trackerRedChi2Diff->Fill(track->normalizedChi2() - muonResidualsFromTrack.trackerRedChi2());
for (std::vector<DetId>::const_iterator chamberId = chamberIds.begin(); chamberId != chamberIds.end(); ++chamberId) {
if (chamberId->det() == DetId::Muon && chamberId->subdetId() == MuonSubdetId::DT) {
DTChamberId dtid(chamberId->rawId());
if (dtid.station() == 1) {
MuonChamberResidual *dt13 = muonResidualsFromTrack.chamberResidual(*chamberId, MuonChamberResidual::kDT13);
if (dt13 != NULL && dt13->numHits() >= m_minDT13Hits) {
int wheel = -1;
if (dtid.wheel() == -2) wheel = kWheelMinus2;
else if (dtid.wheel() == -1) wheel = kWheelMinus1;
else if (dtid.wheel() == 0) wheel = kWheelZero;
else if (dtid.wheel() == 1) wheel = kWheelPlus1;
else if (dtid.wheel() == 2) wheel = kWheelPlus2;
int sector = dtid.sector() - 1;
double resid_x = dt13->global_hitresid(m_layer);
double resid_dxdz = dt13->global_resslope();
if (fabs(resid_x) < 10.) {
th2f_wheelsector[wheel][sector][kDeltaX]->Fill(qoverpt, resid_x*10.);
tprofile_wheelsector[wheel][sector][kDeltaX]->Fill(qoverpt, resid_x*10.);
th2f_wheelsector[wheel][sector][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
tprofile_wheelsector[wheel][sector][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
}
// derivatives are in local coordinates, so these should be, too
resid_x = dt13->hitresid(m_layer);
resid_dxdz = dt13->resslope();
// calculate derivative
TrajectoryStateOnSurface last_tracker_tsos;
double last_tracker_R = 0.;
std::vector<TrajectoryMeasurement> measurements = traj->measurements();
for (std::vector<TrajectoryMeasurement>::const_iterator im = measurements.begin(); im != measurements.end(); ++im) {
TrajectoryStateOnSurface tsos = im->forwardPredictedState();
if (tsos.isValid()) {
GlobalPoint pos = tsos.globalPosition();
if (pos.perp() < 200. && fabs(pos.z()) < 400.) { // if in tracker (cheap, I know...)
if (pos.perp() > last_tracker_R) {
last_tracker_tsos = tsos;
last_tracker_R = pos.perp();
}
}
}
}
if (last_tracker_R > 0.) {
FreeTrajectoryState ts_rebuilt(last_tracker_tsos.globalPosition(),
last_tracker_tsos.globalMomentum(),
last_tracker_tsos.charge(),
&*magneticField);
double factor = (last_tracker_tsos.globalMomentum().mag() + 1.) / last_tracker_tsos.globalMomentum().mag();
FreeTrajectoryState ts_plus1GeV(last_tracker_tsos.globalPosition(),
GlobalVector(factor*last_tracker_tsos.globalMomentum().x(), factor*last_tracker_tsos.globalMomentum().y(), factor*last_tracker_tsos.globalMomentum().z()),
last_tracker_tsos.charge(),
&*magneticField);
TrajectoryStateOnSurface extrapolation_rebuilt = propagator->propagate(ts_rebuilt, globalGeometry->idToDet(dt13->localid(m_layer))->surface());
TrajectoryStateOnSurface extrapolation_plus1GeV = propagator->propagate(ts_plus1GeV, globalGeometry->idToDet(dt13->localid(m_layer))->surface());
if (extrapolation_rebuilt.isValid() && extrapolation_plus1GeV.isValid()) {
double rebuiltx = extrapolation_rebuilt.localPosition().x();
double plus1x = extrapolation_plus1GeV.localPosition().x();
if (fabs(resid_x) < 10.) {
th2f_wheelsector[wheel][sector][kPtErr]->Fill(1./qoverpt, resid_x/(rebuiltx-plus1x)/sqrt(1. + pz*pz/(px*px + py*py))*fabs(qoverpt)*100.);
tprofile_wheelsector[wheel][sector][kPtErr]->Fill(1./qoverpt, resid_x/(rebuiltx-plus1x)/sqrt(1. + pz*pz/(px*px + py*py))*fabs(qoverpt)*100.);
th2f_wheelsector[wheel][sector][kCurvErr]->Fill(qoverpt, resid_x/(rebuiltx-plus1x)*qoverpt/sqrt(px*px + py*py + pz*pz));
tprofile_wheelsector[wheel][sector][kCurvErr]->Fill(qoverpt, resid_x/(rebuiltx-plus1x)*qoverpt/sqrt(px*px + py*py + pz*pz));
}
}
}
} // if it's a good segment
} // if on our chamber
} // if DT
if (chamberId->det() == DetId::Muon && chamberId->subdetId() == MuonSubdetId::CSC) {
CSCDetId cscid(chamberId->rawId());
if (cscid.station() == 1) {
MuonChamberResidual *csc = muonResidualsFromTrack.chamberResidual(*chamberId, MuonChamberResidual::kCSC);
if (csc != NULL && csc->numHits() >= m_minCSCHits) {
int wheel = -1;
int endcap = -1;
if (cscid.endcap() == 1 && cscid.ring() == 1) { wheel = kWheelMEp11; endcap = kEndcapMEp11; }
else if (cscid.endcap() == 1 && cscid.ring() == 2) { wheel = kWheelMEp12; endcap = kEndcapMEp12; }
else if (cscid.endcap() == 1 && cscid.ring() == 3) { wheel = kWheelMEp13; endcap = kEndcapMEp13; }
else if (cscid.endcap() == 1 && cscid.ring() == 4) { wheel = kWheelMEp14; endcap = kEndcapMEp14; }
else if (cscid.endcap() != 1 && cscid.ring() == 1) { wheel = kWheelMEm11; endcap = kEndcapMEm11; }
else if (cscid.endcap() != 1 && cscid.ring() == 2) { wheel = kWheelMEm12; endcap = kEndcapMEm12; }
else if (cscid.endcap() != 1 && cscid.ring() == 3) { wheel = kWheelMEm13; endcap = kEndcapMEm13; }
else if (cscid.endcap() != 1 && cscid.ring() == 4) { wheel = kWheelMEm14; endcap = kEndcapMEm14; }
int chamber = cscid.chamber() - 1;
int sector = cscid.chamber() / 3;
if (cscid.chamber() == 36) sector = 0;
double resid_x = csc->global_hitresid(m_layer);
double resid_dxdz = csc->global_resslope();
if (fabs(resid_x) < 10.) {
th2f_wheelsector[wheel][sector][kDeltaX]->Fill(qoverpt, resid_x*10.);
tprofile_wheelsector[wheel][sector][kDeltaX]->Fill(qoverpt, resid_x*10.);
th2f_wheelsector[wheel][sector][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
tprofile_wheelsector[wheel][sector][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
th2f_endcap[endcap][chamber][kDeltaX]->Fill(qoverpt, resid_x*10.);
tprofile_endcap[endcap][chamber][kDeltaX]->Fill(qoverpt, resid_x*10.);
th2f_endcap[endcap][chamber][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
tprofile_endcap[endcap][chamber][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
if (cscid.chamber() % 2 == 0) {
th2f_evens[endcap][kDeltaX]->Fill(qoverpt, resid_x*10.);
tprofile_evens[endcap][kDeltaX]->Fill(qoverpt, resid_x*10.);
th2f_evens[endcap][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
tprofile_evens[endcap][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
}
else {
th2f_odds[endcap][kDeltaX]->Fill(qoverpt, resid_x*10.);
tprofile_odds[endcap][kDeltaX]->Fill(qoverpt, resid_x*10.);
th2f_odds[endcap][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
tprofile_odds[endcap][kDeltaDxDz]->Fill(qoverpt, resid_dxdz*1000.);
}
}
// derivatives are in local coordinates, so these should be, too
resid_x = csc->hitresid(m_layer);
resid_dxdz = csc->resslope();
// calculate derivative
TrajectoryStateOnSurface last_tracker_tsos;
double last_tracker_absZ = 0.;
std::vector<TrajectoryMeasurement> measurements = traj->measurements();
for (std::vector<TrajectoryMeasurement>::const_iterator im = measurements.begin(); im != measurements.end(); ++im) {
TrajectoryStateOnSurface tsos = im->forwardPredictedState();
if (tsos.isValid()) {
GlobalPoint pos = tsos.globalPosition();
if (pos.perp() < 200. && fabs(pos.z()) < 400.) {
if (fabs(pos.z()) > last_tracker_absZ) {
last_tracker_tsos = tsos;
last_tracker_absZ = fabs(pos.z());
}
}
}
}
if (last_tracker_absZ > 0.) {
FreeTrajectoryState ts_rebuilt(last_tracker_tsos.globalPosition(),
last_tracker_tsos.globalMomentum(),
last_tracker_tsos.charge(),
&*magneticField);
double factor = (last_tracker_tsos.globalMomentum().mag() + 1.) / last_tracker_tsos.globalMomentum().mag();
FreeTrajectoryState ts_plus1GeV(last_tracker_tsos.globalPosition(),
GlobalVector(factor*last_tracker_tsos.globalMomentum().x(), factor*last_tracker_tsos.globalMomentum().y(), factor*last_tracker_tsos.globalMomentum().z()),
last_tracker_tsos.charge(),
&*magneticField);
TrajectoryStateOnSurface extrapolation_rebuilt = propagator->propagate(ts_rebuilt, globalGeometry->idToDet(csc->localid(m_layer))->surface());
TrajectoryStateOnSurface extrapolation_plus1GeV = propagator->propagate(ts_plus1GeV, globalGeometry->idToDet(csc->localid(m_layer))->surface());
if (extrapolation_rebuilt.isValid() && extrapolation_plus1GeV.isValid()) {
double rebuiltx = extrapolation_rebuilt.localPosition().x();
double plus1x = extrapolation_plus1GeV.localPosition().x();
if (fabs(resid_x) < 10.) {
double pterroverpt = resid_x/(rebuiltx-plus1x)/sqrt(1. + pz*pz/(px*px + py*py))*fabs(qoverpt)*100.;
double curverror = resid_x/(rebuiltx-plus1x)*qoverpt/sqrt(px*px + py*py + pz*pz);
th2f_wheelsector[wheel][sector][kPtErr]->Fill(1./qoverpt, pterroverpt);
tprofile_wheelsector[wheel][sector][kPtErr]->Fill(1./qoverpt, pterroverpt);
th2f_wheelsector[wheel][sector][kCurvErr]->Fill(qoverpt, curverror);
tprofile_wheelsector[wheel][sector][kCurvErr]->Fill(qoverpt, curverror);
th2f_endcap[endcap][chamber][kPtErr]->Fill(1./qoverpt, pterroverpt);
tprofile_endcap[endcap][chamber][kPtErr]->Fill(1./qoverpt, pterroverpt);
th2f_endcap[endcap][chamber][kCurvErr]->Fill(qoverpt, curverror);
tprofile_endcap[endcap][chamber][kCurvErr]->Fill(qoverpt, curverror);
if (cscid.chamber() % 2 == 0) {
th2f_evens[endcap][kPtErr]->Fill(1./qoverpt, pterroverpt);
tprofile_evens[endcap][kPtErr]->Fill(1./qoverpt, pterroverpt);
th2f_evens[endcap][kCurvErr]->Fill(qoverpt, curverror);
tprofile_evens[endcap][kCurvErr]->Fill(qoverpt, curverror);
}
else {
th2f_odds[endcap][kPtErr]->Fill(1./qoverpt, pterroverpt);
tprofile_odds[endcap][kPtErr]->Fill(1./qoverpt, pterroverpt);
th2f_odds[endcap][kCurvErr]->Fill(qoverpt, curverror);
tprofile_odds[endcap][kCurvErr]->Fill(qoverpt, curverror);
}
}
}
}
} // if it's a good segment
} // if on our chamber
} // if CSC
} // end loop over chamberIds
} // end if refit is okay
} // end if track pT is within range
} // end loop over tracks
}
bool AlignmentMonitorMuonVsCurvature::m_allowTIDTEC [private] |
Definition at line 56 of file AlignmentMonitorMuonVsCurvature.cc.
Referenced by event().
int AlignmentMonitorMuonVsCurvature::m_layer [private] |
Definition at line 61 of file AlignmentMonitorMuonVsCurvature.cc.
Referenced by event().
double AlignmentMonitorMuonVsCurvature::m_maxDxy [private] |
Definition at line 57 of file AlignmentMonitorMuonVsCurvature.cc.
Referenced by event().
double AlignmentMonitorMuonVsCurvature::m_maxTrackerRedChi2 [private] |
Definition at line 55 of file AlignmentMonitorMuonVsCurvature.cc.
Referenced by event().
int AlignmentMonitorMuonVsCurvature::m_minCSCHits [private] |
Definition at line 60 of file AlignmentMonitorMuonVsCurvature.cc.
Referenced by event().
int AlignmentMonitorMuonVsCurvature::m_minDT13Hits [private] |
Definition at line 58 of file AlignmentMonitorMuonVsCurvature.cc.
Referenced by event().
int AlignmentMonitorMuonVsCurvature::m_minDT2Hits [private] |
Definition at line 59 of file AlignmentMonitorMuonVsCurvature.cc.
int AlignmentMonitorMuonVsCurvature::m_minTrackerHits [private] |
Definition at line 54 of file AlignmentMonitorMuonVsCurvature.cc.
Referenced by event().
double AlignmentMonitorMuonVsCurvature::m_minTrackPt [private] |
Definition at line 53 of file AlignmentMonitorMuonVsCurvature.cc.
Referenced by event().
std::string AlignmentMonitorMuonVsCurvature::m_propagator [private] |
Definition at line 62 of file AlignmentMonitorMuonVsCurvature.cc.
Referenced by event().
TH1F* AlignmentMonitorMuonVsCurvature::th1f_trackerRedChi2 [private] |
Definition at line 112 of file AlignmentMonitorMuonVsCurvature.cc.
TH1F* AlignmentMonitorMuonVsCurvature::th1f_trackerRedChi2Diff [private] |
Definition at line 113 of file AlignmentMonitorMuonVsCurvature.cc.
TH2F* AlignmentMonitorMuonVsCurvature::th2f_endcap[kNumEndcap][36][kNumComponents] [private] |
Definition at line 109 of file AlignmentMonitorMuonVsCurvature.cc.
TH2F* AlignmentMonitorMuonVsCurvature::th2f_evens[kNumEndcap][kNumComponents] [private] |
Definition at line 104 of file AlignmentMonitorMuonVsCurvature.cc.
TH2F* AlignmentMonitorMuonVsCurvature::th2f_odds[kNumEndcap][kNumComponents] [private] |
Definition at line 105 of file AlignmentMonitorMuonVsCurvature.cc.
TH2F* AlignmentMonitorMuonVsCurvature::th2f_wheelsector[kNumWheels][12][kNumComponents] [private] |
Definition at line 101 of file AlignmentMonitorMuonVsCurvature.cc.
TProfile* AlignmentMonitorMuonVsCurvature::tprofile_endcap[kNumEndcap][36][kNumComponents] [private] |
Definition at line 110 of file AlignmentMonitorMuonVsCurvature.cc.
TProfile* AlignmentMonitorMuonVsCurvature::tprofile_evens[kNumEndcap][kNumComponents] [private] |
Definition at line 106 of file AlignmentMonitorMuonVsCurvature.cc.
TProfile* AlignmentMonitorMuonVsCurvature::tprofile_odds[kNumEndcap][kNumComponents] [private] |
Definition at line 107 of file AlignmentMonitorMuonVsCurvature.cc.
TProfile* AlignmentMonitorMuonVsCurvature::tprofile_wheelsector[kNumWheels][12][kNumComponents] [private] |
Definition at line 102 of file AlignmentMonitorMuonVsCurvature.cc.
1.7.3