#include <RecoTracker/DebugTools/src/TestSmoothHits.cc>
Public Member Functions | |
TestSmoothHits (const edm::ParameterSet &) | |
~TestSmoothHits () | |
Private Member Functions | |
virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
virtual void | beginRun (edm::Run &run, const edm::EventSetup &) |
virtual void | endJob () |
std::pair< LocalPoint, LocalVector > | projectHit (const PSimHit &, const StripGeomDetUnit *, const BoundPlane &) |
Private Attributes | |
std::string | builderName |
const edm::ParameterSet | conf_ |
TFile * | file |
edm::ESHandle< TrajectoryFitter > | fit |
std::string | fname |
TH2F * | hChi2_vs_clsize |
TH2F * | hChi2_vs_Process |
std::map< std::string, TH1F * > | hChi2Increment |
TrackerHitAssociator * | hitAssociator |
std::map< std::string, TH1F * > | hPullGM_X_ts |
std::map< std::string, TH1F * > | hPullGM_X_ts_mono |
std::map< std::string, TH1F * > | hPullGM_X_ts_stereo |
std::map< std::string, TH1F * > | hPullGM_Y_ts |
std::map< std::string, TH1F * > | hPullGM_Y_ts_mono |
std::map< std::string, TH1F * > | hPullGM_Y_ts_stereo |
std::map< std::string, TH1F * > | hPullGM_Z_ts |
std::map< std::string, TH1F * > | hPullGM_Z_ts_mono |
std::map< std::string, TH1F * > | hPullGM_Z_ts_stereo |
std::map< std::string, TH1F * > | hPullGP_X_rs |
std::map< std::string, TH1F * > | hPullGP_X_rs_mono |
std::map< std::string, TH1F * > | hPullGP_X_rs_stereo |
std::map< std::string, TH1F * > | hPullGP_X_tr |
std::map< std::string, TH1F * > | hPullGP_X_tr_mono |
std::map< std::string, TH1F * > | hPullGP_X_tr_stereo |
std::map< std::string, TH1F * > | hPullGP_X_ts |
std::map< std::string, TH1F * > | hPullGP_X_ts_mono |
std::map< std::string, TH1F * > | hPullGP_X_ts_stereo |
std::map< std::string, TH1F * > | hPullGP_Y_rs |
std::map< std::string, TH1F * > | hPullGP_Y_rs_mono |
std::map< std::string, TH1F * > | hPullGP_Y_rs_stereo |
std::map< std::string, TH1F * > | hPullGP_Y_tr |
std::map< std::string, TH1F * > | hPullGP_Y_tr_mono |
std::map< std::string, TH1F * > | hPullGP_Y_tr_stereo |
std::map< std::string, TH1F * > | hPullGP_Y_ts |
std::map< std::string, TH1F * > | hPullGP_Y_ts_mono |
std::map< std::string, TH1F * > | hPullGP_Y_ts_stereo |
std::map< std::string, TH1F * > | hPullGP_Z_rs |
std::map< std::string, TH1F * > | hPullGP_Z_rs_mono |
std::map< std::string, TH1F * > | hPullGP_Z_rs_stereo |
std::map< std::string, TH1F * > | hPullGP_Z_tr |
std::map< std::string, TH1F * > | hPullGP_Z_tr_mono |
std::map< std::string, TH1F * > | hPullGP_Z_tr_stereo |
std::map< std::string, TH1F * > | hPullGP_Z_ts |
std::map< std::string, TH1F * > | hPullGP_Z_ts_mono |
std::map< std::string, TH1F * > | hPullGP_Z_ts_stereo |
TH1F * | hTotChi2Increment |
double | maxeta |
double | mineta |
std::string | propagatorName |
edm::ESHandle< TrajectorySmoother > | smooth |
std::string | sname |
std::string | srcName |
edm::ESHandle < TransientTrackingRecHitBuilder > | theBuilder |
edm::ESHandle< TrackerGeometry > | theG |
edm::ESHandle< MagneticField > | theMF |
edm::ESHandle< Propagator > | thePropagator |
edm::Handle < TrackCandidateCollection > | theTCCollection |
std::stringstream | title |
Description: <one line="" class="" summary>="">
Implementation: <Notes on="" implementation>="">
Definition at line 51 of file TestSmoothHits.h.
TestSmoothHits::TestSmoothHits | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 27 of file TestSmoothHits.cc.
References builderName, conf_, fname, edm::ParameterSet::getParameter(), LogTrace, maxeta, mineta, propagatorName, sname, and srcName.
: conf_(iConfig){ LogTrace("TestSmoothHits") << conf_<< std::endl; propagatorName = conf_.getParameter<std::string>("Propagator"); builderName = conf_.getParameter<std::string>("TTRHBuilder"); srcName = conf_.getParameter<std::string>("src"); fname = conf_.getParameter<std::string>("Fitter"); sname = conf_.getParameter<std::string>("Smoother"); mineta = conf_.getParameter<double>("mineta"); maxeta = conf_.getParameter<double>("maxeta"); }
TestSmoothHits::~TestSmoothHits | ( | ) |
Definition at line 39 of file TestSmoothHits.cc.
{}
void TestSmoothHits::analyze | ( | const edm::Event & | iEvent, |
const edm::EventSetup & | iSetup | ||
) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 193 of file TestSmoothHits.cc.
References TrackerHitAssociator::associateHit(), TrajectoryStateOnSurface::cartesianError(), heavyFlavorValidationHarvestingSequence_cff::combiner, GlobalErrorBase< T, ErrorWeightType >::cxx(), GlobalErrorBase< T, ErrorWeightType >::cyy(), GlobalErrorBase< T, ErrorWeightType >::czz(), delta, PTrajectoryStateOnDet::detId(), fit, edm::Event::getByLabel(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), hChi2_vs_clsize, hChi2_vs_Process, hChi2Increment, hitAssociator, hPullGM_X_ts, hPullGM_X_ts_mono, hPullGM_X_ts_stereo, hPullGM_Y_ts, hPullGM_Y_ts_mono, hPullGM_Y_ts_stereo, hPullGM_Z_ts, hPullGM_Z_ts_mono, hPullGM_Z_ts_stereo, hPullGP_X_rs, hPullGP_X_rs_mono, hPullGP_X_rs_stereo, hPullGP_X_tr, hPullGP_X_tr_mono, hPullGP_X_tr_stereo, hPullGP_X_ts, hPullGP_X_ts_mono, hPullGP_X_ts_stereo, hPullGP_Y_rs, hPullGP_Y_rs_mono, hPullGP_Y_rs_stereo, hPullGP_Y_tr, hPullGP_Y_tr_mono, hPullGP_Y_tr_stereo, hPullGP_Y_ts, hPullGP_Y_ts_mono, hPullGP_Y_ts_stereo, hPullGP_Z_rs, hPullGP_Z_rs_mono, hPullGP_Z_rs_stereo, hPullGP_Z_tr, hPullGP_Z_tr_mono, hPullGP_Z_tr_stereo, hPullGP_Z_ts, hPullGP_Z_ts_mono, hPullGP_Z_ts_stereo, hTotChi2Increment, i, TrajectoryStateOnSurface::isValid(), PSimHit::localPosition(), LogTrace, m, mag(), CartesianTrajectoryError::matrix(), maxeta, mineta, PSimHit::momentumAtEntry(), CartesianTrajectoryError::position(), PSimHit::processType(), edm::ESHandle< T >::product(), projectHit(), TrackCandidate::recHits(), query::result, alignCSCRings::s, TrackCandidate::seed(), smooth, mathSSE::sqrt(), srcName, theBuilder, theG, theMF, thePropagator, theTCCollection, title, Surface::toGlobal(), TrackCandidate::trajectoryStateOnDet(), trajectoryStateTransform::transientState(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{ LogTrace("TestSmoothHits") << "new event" << std::endl; iEvent.getByLabel(srcName,theTCCollection ); hitAssociator = new TrackerHitAssociator(iEvent); TrajectoryStateCombiner combiner; for (TrackCandidateCollection::const_iterator i=theTCCollection->begin(); i!=theTCCollection->end();i++){ LogTrace("TestSmoothHits") << "new candidate" << std::endl; const TrackCandidate * theTC = &(*i); PTrajectoryStateOnDet state = theTC->trajectoryStateOnDet(); const TrackCandidate::range& recHitVec=theTC->recHits(); //convert PTrajectoryStateOnDet to TrajectoryStateOnSurface DetId detId(state.detId()); TrajectoryStateOnSurface theTSOS= trajectoryStateTransform::transientState(state, &(theG->idToDet(detId)->surface()),theMF.product()); if (theTSOS.globalMomentum().eta()>maxeta || theTSOS.globalMomentum().eta()<mineta) continue; //convert the TrackingRecHit vector to a TransientTrackingRecHit vector TransientTrackingRecHit::RecHitContainer hits; for (edm::OwnVector<TrackingRecHit>::const_iterator i=recHitVec.first; i!=recHitVec.second; i++){ hits.push_back(theBuilder->build(&(*i) )); } //call the fitter std::vector<Trajectory> fitted = fit->fit(theTC->seed(), hits, theTSOS); //call the smoother std::vector<Trajectory> result; for(std::vector<Trajectory>::iterator it = fitted.begin(); it != fitted.end(); it++) { std::vector<Trajectory> smoothed = smooth->trajectories(*it); result.insert(result.end(), smoothed.begin(), smoothed.end()); } if (result.size()==0) continue; std::vector<TrajectoryMeasurement> vtm = result[0].measurements(); TSOS lastState = theTSOS; for (std::vector<TrajectoryMeasurement>::iterator tm=vtm.begin(); tm!=vtm.end();tm++){ TransientTrackingRecHit::ConstRecHitPointer rhit = tm->recHit(); if ((rhit)->isValid()==0&&rhit->det()!=0) continue; LogTrace("TestSmoothHits") << "new hit" ; int subdetId = rhit->det()->geographicalId().subdetId(); int layerId = 0; DetId id = rhit->det()->geographicalId(); if (id.subdetId()==3) layerId = ((TIBDetId)(id)).layer(); if (id.subdetId()==5) layerId = ((TOBDetId)(id)).layer(); if (id.subdetId()==1) layerId = ((PXBDetId)(id)).layer(); if (id.subdetId()==4) layerId = ((TIDDetId)(id)).wheel(); if (id.subdetId()==6) layerId = ((TECDetId)(id)).wheel(); if (id.subdetId()==2) layerId = ((PXFDetId)(id)).disk(); LogTrace("TestSmoothHits") << "subdetId=" << subdetId << " layerId=" << layerId ; double delta = 99999; LocalPoint rhitLPv = rhit->localPosition(); std::vector<PSimHit> assSimHits = hitAssociator->associateHit(*(rhit->hit())); if (assSimHits.size()==0) continue; PSimHit shit; for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){ if ((m->localPosition()-rhitLPv).mag()<delta) { shit=*m; delta = (m->localPosition()-rhitLPv).mag(); } } TSOS currentState = combiner(tm->backwardPredictedState(), tm->forwardPredictedState()); TSOS updatedState = tm->updatedState(); //plot chi2 increment double chi2increment = tm->estimate(); LogTrace("TestSmoothHits") << "tm->estimate()=" << tm->estimate(); title.str(""); title << "Chi2Increment_" << subdetId << "-" << layerId; hChi2Increment[title.str()]->Fill( chi2increment ); hTotChi2Increment->Fill( chi2increment ); hChi2_vs_Process->Fill( chi2increment, shit.processType() ); if (dynamic_cast<const SiPixelRecHit*>(rhit->hit())) hChi2_vs_clsize->Fill( chi2increment, ((const SiPixelRecHit*)(rhit->hit()))->cluster()->size() ); if (dynamic_cast<const SiStripRecHit2D*>(rhit->hit())) hChi2_vs_clsize->Fill( chi2increment, ((const SiStripRecHit2D*)(rhit->hit()))->cluster()->amplitudes().size() ); //test hits const Surface * surf = &( (rhit)->det()->surface() ); LocalVector shitLMom; LocalPoint shitLPos; if (dynamic_cast<const SiStripMatchedRecHit2D*>(rhit->hit())) { double rechitmatchedx = rhit->localPosition().x(); double rechitmatchedy = rhit->localPosition().y(); double mindist = 999999; float distx, disty; std::pair<LocalPoint,LocalVector> closestPair; const StripGeomDetUnit* stripDet =(StripGeomDetUnit*) ((const GluedGeomDet *)(rhit)->det())->stereoDet(); const BoundPlane& plane = (rhit)->det()->surface(); for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++) { //project simhit; std::pair<LocalPoint,LocalVector> hitPair = projectHit((*m),stripDet,plane); distx = fabs(rechitmatchedx - hitPair.first.x()); disty = fabs(rechitmatchedy - hitPair.first.y()); double dist = distx*distx+disty*disty; if(sqrt(dist)<mindist){ mindist = dist; closestPair = hitPair; } } shitLPos = closestPair.first; shitLMom = closestPair.second; } else { shitLPos = shit.localPosition(); shitLMom = shit.momentumAtEntry(); } GlobalVector shitGMom = surf->toGlobal(shitLMom); GlobalPoint shitGPos = surf->toGlobal(shitLPos); // if (dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())) { // double rechitmatchedx = rhit->localPosition().x(); // double rechitmatchedy = rhit->localPosition().y(); // double mindist = 999999; // double distx, disty; // std::pair<LocalPoint,LocalVector> closestPair; // const StripGeomDetUnit* stripDet =(StripGeomDetUnit*) ((const GluedGeomDet *)(rhit)->det())->stereoDet(); // const BoundPlane& plane = (rhit)->det()->surface(); // for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){ // const PSimHit& hit = *m; // const StripTopology& topol = stripDet->specificTopology(); // GlobalPoint globalpos = stripDet->surface().toGlobal(hit.localPosition()); // LocalPoint localHit = plane.toLocal(globalpos); // //track direction // LocalVector locdir=hit.localDirection(); // //rotate track in new frame // GlobalVector globaldir= stripDet->surface().toGlobal(locdir); // LocalVector dir=plane.toLocal(globaldir); // float scale = -localHit.z() / dir.z(); // LocalPoint projectedPos = localHit + scale*dir; // float selfAngle = topol.stripAngle( topol.strip( hit.localPosition())); // LocalVector stripDir( sin(selfAngle), cos(selfAngle), 0); // vector along strip in hit frame // LocalVector localStripDir = LocalVector( plane.toLocal(stripDet->surface().toGlobal( stripDir))); // std::pair<LocalPoint,LocalVector> hitPair( projectedPos, localStripDir); // distx = fabs(rechitmatchedx - hitPair.first.x()); // disty = fabs(rechitmatchedy - hitPair.first.y()); // double dist = distx*distx+disty*disty; // if(sqrt(dist)<mindist){ // mindist = dist; // closestPair = hitPair; // } // } // shitLPos = closestPair.first; // shitLMom = closestPair.second; // } else { // shitLPos = shit.localPosition(); // shitLMom = shit.momentumAtEntry(); // } // GlobalVector shitGMom = surf->toGlobal(shitLMom); // GlobalPoint shitGPos = surf->toGlobal(shitLPos); GlobalVector tsosGMom = currentState.globalMomentum(); GlobalError tsosGMEr(currentState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3,3)); GlobalPoint tsosGPos = currentState.globalPosition(); GlobalError tsosGPEr = currentState.cartesianError().position(); GlobalPoint rhitGPos = (rhit)->globalPosition(); GlobalError rhitGPEr = (rhit)->globalPositionError(); double pullGPX_rs = (rhitGPos.x()-shitGPos.x())/sqrt(rhitGPEr.cxx()); double pullGPY_rs = (rhitGPos.y()-shitGPos.y())/sqrt(rhitGPEr.cyy()); double pullGPZ_rs = (rhitGPos.z()-shitGPos.z())/sqrt(rhitGPEr.czz()); //double pullGPX_rs = (rhitGPos.x()-shitGPos.x()); //double pullGPY_rs = (rhitGPos.y()-shitGPos.y()); //double pullGPZ_rs = (rhitGPos.z()-shitGPos.z()); LogTrace("TestSmoothHits") << "rs" << std::endl; title.str(""); title << "PullGP_X_" << subdetId << "-" << layerId << "_rs"; hPullGP_X_rs[title.str()]->Fill( pullGPX_rs ); title.str(""); title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs"; hPullGP_Y_rs[title.str()]->Fill( pullGPY_rs ); title.str(""); title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs"; hPullGP_Z_rs[title.str()]->Fill( pullGPZ_rs ); double pullGPX_tr = (tsosGPos.x()-rhitGPos.x())/sqrt(tsosGPEr.cxx()+rhitGPEr.cxx()); double pullGPY_tr = (tsosGPos.y()-rhitGPos.y())/sqrt(tsosGPEr.cyy()+rhitGPEr.cyy()); double pullGPZ_tr = (tsosGPos.z()-rhitGPos.z())/sqrt(tsosGPEr.czz()+rhitGPEr.czz()); //double pullGPX_tr = (tsosGPos.x()-rhitGPos.x()); //double pullGPY_tr = (tsosGPos.y()-rhitGPos.y()); //double pullGPZ_tr = (tsosGPos.z()-rhitGPos.z()); LogTrace("TestSmoothHits") << "tr" << std::endl; title.str(""); title << "PullGP_X_" << subdetId << "-" << layerId << "_tr"; hPullGP_X_tr[title.str()]->Fill( pullGPX_tr ); title.str(""); title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr"; hPullGP_Y_tr[title.str()]->Fill( pullGPY_tr ); title.str(""); title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr"; hPullGP_Z_tr[title.str()]->Fill( pullGPZ_tr ); double pullGPX_ts = (tsosGPos.x()-shitGPos.x())/sqrt(tsosGPEr.cxx()); double pullGPY_ts = (tsosGPos.y()-shitGPos.y())/sqrt(tsosGPEr.cyy()); double pullGPZ_ts = (tsosGPos.z()-shitGPos.z())/sqrt(tsosGPEr.czz()); //double pullGPX_ts = (tsosGPos.x()-shitGPos.x()); //double pullGPY_ts = (tsosGPos.y()-shitGPos.y()); //double pullGPZ_ts = (tsosGPos.z()-shitGPos.z()); LogTrace("TestSmoothHits") << "ts1" << std::endl; title.str(""); title << "PullGP_X_" << subdetId << "-" << layerId << "_ts"; hPullGP_X_ts[title.str()]->Fill( pullGPX_ts ); title.str(""); title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts"; hPullGP_Y_ts[title.str()]->Fill( pullGPY_ts ); title.str(""); title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts"; hPullGP_Z_ts[title.str()]->Fill( pullGPZ_ts ); double pullGMX_ts = (tsosGMom.x()-shitGMom.x())/sqrt(tsosGMEr.cxx()); double pullGMY_ts = (tsosGMom.y()-shitGMom.y())/sqrt(tsosGMEr.cyy()); double pullGMZ_ts = (tsosGMom.z()-shitGMom.z())/sqrt(tsosGMEr.czz()); //double pullGMX_ts = (tsosGMom.x()-shitGMom.x()); //double pullGMY_ts = (tsosGMom.y()-shitGMom.y()); //double pullGMZ_ts = (tsosGMom.z()-shitGMom.z()); LogTrace("TestSmoothHits") << "ts2" << std::endl; title.str(""); title << "PullGM_X_" << subdetId << "-" << layerId << "_ts"; hPullGM_X_ts[title.str()]->Fill( pullGMX_ts ); title.str(""); title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts"; hPullGM_Y_ts[title.str()]->Fill( pullGMY_ts ); title.str(""); title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts"; hPullGM_Z_ts[title.str()]->Fill( pullGMZ_ts ); if (dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())) { //mono LogTrace("TestSmoothHits") << "MONO HIT" << std::endl; auto m = dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())->monoHit(); CTTRHp tMonoHit = theBuilder->build(&m); if (tMonoHit==0) continue; vector<PSimHit> assMonoSimHits = hitAssociator->associateHit(*tMonoHit->hit()); if (assMonoSimHits.size()==0) continue; const PSimHit sMonoHit = *(assSimHits.begin()); const Surface * monoSurf = &( tMonoHit->det()->surface() ); if (monoSurf==0) continue; TSOS monoState = thePropagator->propagate(lastState,*monoSurf); if (monoState.isValid()==0) continue; LocalVector monoShitLMom = sMonoHit.momentumAtEntry(); GlobalVector monoShitGMom = monoSurf->toGlobal(monoShitLMom); LocalPoint monoShitLPos = sMonoHit.localPosition(); GlobalPoint monoShitGPos = monoSurf->toGlobal(monoShitLPos); GlobalVector monoTsosGMom = monoState.globalMomentum(); GlobalError monoTsosGMEr(monoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3,3)); GlobalPoint monoTsosGPos = monoState.globalPosition(); GlobalError monoTsosGPEr = monoState.cartesianError().position(); GlobalPoint monoRhitGPos = tMonoHit->globalPosition(); GlobalError monoRhitGPEr = tMonoHit->globalPositionError(); double pullGPX_rs_mono = (monoRhitGPos.x()-monoShitGPos.x())/sqrt(monoRhitGPEr.cxx()); double pullGPY_rs_mono = (monoRhitGPos.y()-monoShitGPos.y())/sqrt(monoRhitGPEr.cyy()); double pullGPZ_rs_mono = (monoRhitGPos.z()-monoShitGPos.z())/sqrt(monoRhitGPEr.czz()); //double pullGPX_rs_mono = (monoRhitGPos.x()-monoShitGPos.x()); //double pullGPY_rs_mono = (monoRhitGPos.y()-monoShitGPos.y()); //double pullGPZ_rs_mono = (monoRhitGPos.z()-monoShitGPos.z()); title.str(""); title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_mono"; hPullGP_X_rs_mono[title.str()]->Fill( pullGPX_rs_mono ); title.str(""); title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_mono"; hPullGP_Y_rs_mono[title.str()]->Fill( pullGPY_rs_mono ); title.str(""); title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_mono"; hPullGP_Z_rs_mono[title.str()]->Fill( pullGPZ_rs_mono ); double pullGPX_tr_mono = (monoTsosGPos.x()-monoRhitGPos.x())/sqrt(monoTsosGPEr.cxx()+monoRhitGPEr.cxx()); double pullGPY_tr_mono = (monoTsosGPos.y()-monoRhitGPos.y())/sqrt(monoTsosGPEr.cyy()+monoRhitGPEr.cyy()); double pullGPZ_tr_mono = (monoTsosGPos.z()-monoRhitGPos.z())/sqrt(monoTsosGPEr.czz()+monoRhitGPEr.czz()); //double pullGPX_tr_mono = (monoTsosGPos.x()-monoRhitGPos.x()); //double pullGPY_tr_mono = (monoTsosGPos.y()-monoRhitGPos.y()); //double pullGPZ_tr_mono = (monoTsosGPos.z()-monoRhitGPos.z()); title.str(""); title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_mono"; hPullGP_X_tr_mono[title.str()]->Fill( pullGPX_tr_mono ); title.str(""); title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_mono"; hPullGP_Y_tr_mono[title.str()]->Fill( pullGPY_tr_mono ); title.str(""); title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_mono"; hPullGP_Z_tr_mono[title.str()]->Fill( pullGPZ_tr_mono ); double pullGPX_ts_mono = (monoTsosGPos.x()-monoShitGPos.x())/sqrt(monoTsosGPEr.cxx()); double pullGPY_ts_mono = (monoTsosGPos.y()-monoShitGPos.y())/sqrt(monoTsosGPEr.cyy()); double pullGPZ_ts_mono = (monoTsosGPos.z()-monoShitGPos.z())/sqrt(monoTsosGPEr.czz()); //double pullGPX_ts_mono = (monoTsosGPos.x()-monoShitGPos.x()); //double pullGPY_ts_mono = (monoTsosGPos.y()-monoShitGPos.y()); //double pullGPZ_ts_mono = (monoTsosGPos.z()-monoShitGPos.z()); title.str(""); title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_mono"; hPullGP_X_ts_mono[title.str()]->Fill( pullGPX_ts_mono ); title.str(""); title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_mono"; hPullGP_Y_ts_mono[title.str()]->Fill( pullGPY_ts_mono ); title.str(""); title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_mono"; hPullGP_Z_ts_mono[title.str()]->Fill( pullGPZ_ts_mono ); double pullGMX_ts_mono = (monoTsosGMom.x()-monoShitGMom.x())/sqrt(monoTsosGMEr.cxx()); double pullGMY_ts_mono = (monoTsosGMom.y()-monoShitGMom.y())/sqrt(monoTsosGMEr.cyy()); double pullGMZ_ts_mono = (monoTsosGMom.z()-monoShitGMom.z())/sqrt(monoTsosGMEr.czz()); //double pullGMX_ts_mono = (monoTsosGMom.x()-monoShitGMom.x()); //double pullGMY_ts_mono = (monoTsosGMom.y()-monoShitGMom.y()); //double pullGMZ_ts_mono = (monoTsosGMom.z()-monoShitGMom.z()); title.str(""); title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_mono"; hPullGM_X_ts_mono[title.str()]->Fill( pullGMX_ts_mono ); title.str(""); title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_mono"; hPullGM_Y_ts_mono[title.str()]->Fill( pullGMY_ts_mono ); title.str(""); title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_mono"; hPullGM_Z_ts_mono[title.str()]->Fill( pullGMZ_ts_mono ); //stereo LogTrace("TestSmoothHits") << "STEREO HIT" << std::endl; auto s = dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())->stereoHit(); CTTRHp tStereoHit = theBuilder->build(&s); if (tStereoHit==0) continue; vector<PSimHit> assStereoSimHits = hitAssociator->associateHit(*tStereoHit->hit()); if (assStereoSimHits.size()==0) continue; const PSimHit sStereoHit = *(assSimHits.begin()); const Surface * stereoSurf = &( tStereoHit->det()->surface() ); if (stereoSurf==0) continue; TSOS stereoState = thePropagator->propagate(lastState,*stereoSurf); if (stereoState.isValid()==0) continue; LocalVector stereoShitLMom = sStereoHit.momentumAtEntry(); GlobalVector stereoShitGMom = stereoSurf->toGlobal(stereoShitLMom); LocalPoint stereoShitLPos = sStereoHit.localPosition(); GlobalPoint stereoShitGPos = stereoSurf->toGlobal(stereoShitLPos); GlobalVector stereoTsosGMom = stereoState.globalMomentum(); GlobalError stereoTsosGMEr(stereoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3,3)); GlobalPoint stereoTsosGPos = stereoState.globalPosition(); GlobalError stereoTsosGPEr = stereoState.cartesianError().position(); GlobalPoint stereoRhitGPos = tStereoHit->globalPosition(); GlobalError stereoRhitGPEr = tStereoHit->globalPositionError(); double pullGPX_rs_stereo = (stereoRhitGPos.x()-stereoShitGPos.x())/sqrt(stereoRhitGPEr.cxx()); double pullGPY_rs_stereo = (stereoRhitGPos.y()-stereoShitGPos.y())/sqrt(stereoRhitGPEr.cyy()); double pullGPZ_rs_stereo = (stereoRhitGPos.z()-stereoShitGPos.z())/sqrt(stereoRhitGPEr.czz()); //double pullGPX_rs_stereo = (stereoRhitGPos.x()-stereoShitGPos.x()); //double pullGPY_rs_stereo = (stereoRhitGPos.y()-stereoShitGPos.y()); //double pullGPZ_rs_stereo = (stereoRhitGPos.z()-stereoShitGPos.z()); title.str(""); title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_stereo"; hPullGP_X_rs_stereo[title.str()]->Fill( pullGPX_rs_stereo ); title.str(""); title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_stereo"; hPullGP_Y_rs_stereo[title.str()]->Fill( pullGPY_rs_stereo ); title.str(""); title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_stereo"; hPullGP_Z_rs_stereo[title.str()]->Fill( pullGPZ_rs_stereo ); double pullGPX_tr_stereo = (stereoTsosGPos.x()-stereoRhitGPos.x())/sqrt(stereoTsosGPEr.cxx()+stereoRhitGPEr.cxx()); double pullGPY_tr_stereo = (stereoTsosGPos.y()-stereoRhitGPos.y())/sqrt(stereoTsosGPEr.cyy()+stereoRhitGPEr.cyy()); double pullGPZ_tr_stereo = (stereoTsosGPos.z()-stereoRhitGPos.z())/sqrt(stereoTsosGPEr.czz()+stereoRhitGPEr.czz()); //double pullGPX_tr_stereo = (stereoTsosGPos.x()-stereoRhitGPos.x()); //double pullGPY_tr_stereo = (stereoTsosGPos.y()-stereoRhitGPos.y()); //double pullGPZ_tr_stereo = (stereoTsosGPos.z()-stereoRhitGPos.z()); title.str(""); title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_stereo"; hPullGP_X_tr_stereo[title.str()]->Fill( pullGPX_tr_stereo ); title.str(""); title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_stereo"; hPullGP_Y_tr_stereo[title.str()]->Fill( pullGPY_tr_stereo ); title.str(""); title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_stereo"; hPullGP_Z_tr_stereo[title.str()]->Fill( pullGPZ_tr_stereo ); double pullGPX_ts_stereo = (stereoTsosGPos.x()-stereoShitGPos.x())/sqrt(stereoTsosGPEr.cxx()); double pullGPY_ts_stereo = (stereoTsosGPos.y()-stereoShitGPos.y())/sqrt(stereoTsosGPEr.cyy()); double pullGPZ_ts_stereo = (stereoTsosGPos.z()-stereoShitGPos.z())/sqrt(stereoTsosGPEr.czz()); //double pullGPX_ts_stereo = (stereoTsosGPos.x()-stereoShitGPos.x()); //double pullGPY_ts_stereo = (stereoTsosGPos.y()-stereoShitGPos.y()); //double pullGPZ_ts_stereo = (stereoTsosGPos.z()-stereoShitGPos.z()); title.str(""); title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_stereo"; hPullGP_X_ts_stereo[title.str()]->Fill( pullGPX_ts_stereo ); title.str(""); title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_stereo"; hPullGP_Y_ts_stereo[title.str()]->Fill( pullGPY_ts_stereo ); title.str(""); title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_stereo"; hPullGP_Z_ts_stereo[title.str()]->Fill( pullGPZ_ts_stereo ); double pullGMX_ts_stereo = (stereoTsosGMom.x()-stereoShitGMom.x())/sqrt(stereoTsosGMEr.cxx()); double pullGMY_ts_stereo = (stereoTsosGMom.y()-stereoShitGMom.y())/sqrt(stereoTsosGMEr.cyy()); double pullGMZ_ts_stereo = (stereoTsosGMom.z()-stereoShitGMom.z())/sqrt(stereoTsosGMEr.czz()); //double pullGMX_ts_stereo = (stereoTsosGMom.x()-stereoShitGMom.x()); //double pullGMY_ts_stereo = (stereoTsosGMom.y()-stereoShitGMom.y()); //double pullGMZ_ts_stereo = (stereoTsosGMom.z()-stereoShitGMom.z()); title.str(""); title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_stereo"; hPullGM_X_ts_stereo[title.str()]->Fill( pullGMX_ts_stereo ); title.str(""); title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_stereo"; hPullGM_Y_ts_stereo[title.str()]->Fill( pullGMY_ts_stereo ); title.str(""); title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_stereo"; hPullGM_Z_ts_stereo[title.str()]->Fill( pullGMZ_ts_stereo ); } lastState = updatedState; //#endif } } delete hitAssociator; LogTrace("TestSmoothHits") << "end of event" << std::endl; }
void TestSmoothHits::beginRun | ( | edm::Run & | run, |
const edm::EventSetup & | iSetup | ||
) | [private, virtual] |
Definition at line 41 of file TestSmoothHits.cc.
References builderName, file, fit, fname, edm::EventSetup::get(), hChi2_vs_clsize, hChi2_vs_Process, hChi2Increment, hPullGM_X_ts, hPullGM_X_ts_mono, hPullGM_X_ts_stereo, hPullGM_Y_ts, hPullGM_Y_ts_mono, hPullGM_Y_ts_stereo, hPullGM_Z_ts, hPullGM_Z_ts_mono, hPullGM_Z_ts_stereo, hPullGP_X_rs, hPullGP_X_rs_mono, hPullGP_X_rs_stereo, hPullGP_X_tr, hPullGP_X_tr_mono, hPullGP_X_tr_stereo, hPullGP_X_ts, hPullGP_X_ts_mono, hPullGP_X_ts_stereo, hPullGP_Y_rs, hPullGP_Y_rs_mono, hPullGP_Y_rs_stereo, hPullGP_Y_tr, hPullGP_Y_tr_mono, hPullGP_Y_tr_stereo, hPullGP_Y_ts, hPullGP_Y_ts_mono, hPullGP_Y_ts_stereo, hPullGP_Z_rs, hPullGP_Z_rs_mono, hPullGP_Z_rs_stereo, hPullGP_Z_tr, hPullGP_Z_tr_mono, hPullGP_Z_tr_stereo, hPullGP_Z_ts, hPullGP_Z_ts_mono, hPullGP_Z_ts_stereo, hTotChi2Increment, i, j, propagatorName, smooth, sname, theBuilder, theG, theMF, thePropagator, and title.
{ iSetup.get<TrackerDigiGeometryRecord>().get(theG); iSetup.get<IdealMagneticFieldRecord>().get(theMF); iSetup.get<TrackingComponentsRecord>().get(propagatorName,thePropagator); iSetup.get<TransientRecHitRecord>().get(builderName,theBuilder); iSetup.get<TrajectoryFitter::Record>().get(fname, fit); iSetup.get<TrajectoryFitter::Record>().get(sname, smooth); file = new TFile("testSmoothHits.root","recreate"); for (int i=0; i!=6; i++) for (int j=0; j!=9; j++){ if (i==0 && j>2) break; if (i==1 && j>1) break; if (i==2 && j>3) break; if (i==3 && j>2) break; if (i==4 && j>5) break; if (i==5 && j>8) break; title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts"; hPullGP_X_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts"; hPullGP_Y_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts"; hPullGP_Z_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "Chi2Increment_" << i+1 << "-" << j+1; hChi2Increment[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,0,100); title.str(""); title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts"; hPullGM_X_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts"; hPullGM_Y_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts"; hPullGM_Z_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr"; hPullGP_X_tr[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr"; hPullGP_Y_tr[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr"; hPullGP_Z_tr[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs"; hPullGP_X_rs[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs"; hPullGP_Y_rs[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs"; hPullGP_Z_rs[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); if ( ((i==2||i==4)&&(j==0||j==1)) || (i==3||i==5) ){ //mono title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGP_X_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGP_Y_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGP_Z_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGM_X_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGM_Y_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGM_Z_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr_mono"; hPullGP_X_tr_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr_mono"; hPullGP_Y_tr_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr_mono"; hPullGP_Z_tr_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs_mono"; hPullGP_X_rs_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs_mono"; hPullGP_Y_rs_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs_mono"; hPullGP_Z_rs_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); //stereo title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGP_X_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGP_Y_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGP_Z_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGM_X_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGM_Y_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGM_Z_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr_stereo"; hPullGP_X_tr_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr_stereo"; hPullGP_Y_tr_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr_stereo"; hPullGP_Z_tr_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs_stereo"; hPullGP_X_rs_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs_stereo"; hPullGP_Y_rs_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs_stereo"; hPullGP_Z_rs_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50); } } hTotChi2Increment = new TH1F("TotChi2Increment","TotChi2Increment",1000,0,100); hChi2_vs_Process = new TH2F("Chi2_vs_Process","Chi2_vs_Process",1000,0,100,17,-0.5,16.5); hChi2_vs_clsize = new TH2F("Chi2_vs_clsize","Chi2_vs_clsize",1000,0,100,17,-0.5,16.5); }
void TestSmoothHits::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 640 of file TestSmoothHits.cc.
References file, hChi2_vs_clsize, hChi2_vs_Process, hChi2Increment, hPullGM_X_ts, hPullGM_X_ts_mono, hPullGM_X_ts_stereo, hPullGM_Y_ts, hPullGM_Y_ts_mono, hPullGM_Y_ts_stereo, hPullGM_Z_ts, hPullGM_Z_ts_mono, hPullGM_Z_ts_stereo, hPullGP_X_rs, hPullGP_X_rs_mono, hPullGP_X_rs_stereo, hPullGP_X_tr, hPullGP_X_tr_mono, hPullGP_X_tr_stereo, hPullGP_X_ts, hPullGP_X_ts_mono, hPullGP_X_ts_stereo, hPullGP_Y_rs, hPullGP_Y_rs_mono, hPullGP_Y_rs_stereo, hPullGP_Y_tr, hPullGP_Y_tr_mono, hPullGP_Y_tr_stereo, hPullGP_Y_ts, hPullGP_Y_ts_mono, hPullGP_Y_ts_stereo, hPullGP_Z_rs, hPullGP_Z_rs_mono, hPullGP_Z_rs_stereo, hPullGP_Z_tr, hPullGP_Z_tr_mono, hPullGP_Z_tr_stereo, hPullGP_Z_ts, hPullGP_Z_ts_mono, hPullGP_Z_ts_stereo, hTotChi2Increment, i, j, and title.
{ //file->Write(); TDirectory * chi2i = file->mkdir("Chi2_Increment"); TDirectory * gp_ts = file->mkdir("GP_TSOS-SimHit"); TDirectory * gm_ts = file->mkdir("GM_TSOS-SimHit"); TDirectory * gp_tr = file->mkdir("GP_TSOS-RecHit"); TDirectory * gp_rs = file->mkdir("GP_RecHit-SimHit"); TDirectory * gp_tsx = gp_ts->mkdir("X"); TDirectory * gp_tsy = gp_ts->mkdir("Y"); TDirectory * gp_tsz = gp_ts->mkdir("Z"); TDirectory * gm_tsx = gm_ts->mkdir("X"); TDirectory * gm_tsy = gm_ts->mkdir("Y"); TDirectory * gm_tsz = gm_ts->mkdir("Z"); TDirectory * gp_trx = gp_tr->mkdir("X"); TDirectory * gp_try = gp_tr->mkdir("Y"); TDirectory * gp_trz = gp_tr->mkdir("Z"); TDirectory * gp_rsx = gp_rs->mkdir("X"); TDirectory * gp_rsy = gp_rs->mkdir("Y"); TDirectory * gp_rsz = gp_rs->mkdir("Z"); TDirectory * gp_tsx_mono = gp_ts->mkdir("MONOX"); TDirectory * gp_tsy_mono = gp_ts->mkdir("MONOY"); TDirectory * gp_tsz_mono = gp_ts->mkdir("MONOZ"); TDirectory * gm_tsx_mono = gm_ts->mkdir("MONOX"); TDirectory * gm_tsy_mono = gm_ts->mkdir("MONOY"); TDirectory * gm_tsz_mono = gm_ts->mkdir("MONOZ"); TDirectory * gp_trx_mono = gp_tr->mkdir("MONOX"); TDirectory * gp_try_mono = gp_tr->mkdir("MONOY"); TDirectory * gp_trz_mono = gp_tr->mkdir("MONOZ"); TDirectory * gp_rsx_mono = gp_rs->mkdir("MONOX"); TDirectory * gp_rsy_mono = gp_rs->mkdir("MONOY"); TDirectory * gp_rsz_mono = gp_rs->mkdir("MONOZ"); TDirectory * gp_tsx_stereo = gp_ts->mkdir("STEREOX"); TDirectory * gp_tsy_stereo = gp_ts->mkdir("STEREOY"); TDirectory * gp_tsz_stereo = gp_ts->mkdir("STEREOZ"); TDirectory * gm_tsx_stereo = gm_ts->mkdir("STEREOX"); TDirectory * gm_tsy_stereo = gm_ts->mkdir("STEREOY"); TDirectory * gm_tsz_stereo = gm_ts->mkdir("STEREOZ"); TDirectory * gp_trx_stereo = gp_tr->mkdir("STEREOX"); TDirectory * gp_try_stereo = gp_tr->mkdir("STEREOY"); TDirectory * gp_trz_stereo = gp_tr->mkdir("STEREOZ"); TDirectory * gp_rsx_stereo = gp_rs->mkdir("STEREOX"); TDirectory * gp_rsy_stereo = gp_rs->mkdir("STEREOY"); TDirectory * gp_rsz_stereo = gp_rs->mkdir("STEREOZ"); chi2i->cd(); hTotChi2Increment->Write(); hChi2_vs_Process->Write(); hChi2_vs_clsize->Write(); for (int i=0; i!=6; i++) for (int j=0; j!=9; j++){ if (i==0 && j>2) break; if (i==1 && j>1) break; if (i==2 && j>3) break; if (i==3 && j>2) break; if (i==4 && j>5) break; if (i==5 && j>8) break; chi2i->cd(); title.str(""); title << "Chi2Increment_" << i+1 << "-" << j+1; hChi2Increment[title.str()]->Write(); gp_ts->cd(); gp_tsx->cd(); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts"; hPullGP_X_ts[title.str()]->Write(); gp_tsy->cd(); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts"; hPullGP_Y_ts[title.str()]->Write(); gp_tsz->cd(); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts"; hPullGP_Z_ts[title.str()]->Write(); gm_ts->cd(); gm_tsx->cd(); title.str(""); title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts"; hPullGM_X_ts[title.str()]->Write(); gm_tsy->cd(); title.str(""); title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts"; hPullGM_Y_ts[title.str()]->Write(); gm_tsz->cd(); title.str(""); title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts"; hPullGM_Z_ts[title.str()]->Write(); gp_tr->cd(); gp_trx->cd(); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr"; hPullGP_X_tr[title.str()]->Write(); gp_try->cd(); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr"; hPullGP_Y_tr[title.str()]->Write(); gp_trz->cd(); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr"; hPullGP_Z_tr[title.str()]->Write(); gp_rs->cd(); gp_rsx->cd(); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs"; hPullGP_X_rs[title.str()]->Write(); gp_rsy->cd(); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs"; hPullGP_Y_rs[title.str()]->Write(); gp_rsz->cd(); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs"; hPullGP_Z_rs[title.str()]->Write(); if ( ((i==2||i==4)&&(j==0||j==1)) || (i==3||i==5) ){ //mono gp_ts->cd(); gp_tsx_mono->cd(); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGP_X_ts_mono[title.str()]->Write(); gp_tsy_mono->cd(); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGP_Y_ts_mono[title.str()]->Write(); gp_tsz_mono->cd(); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGP_Z_ts_mono[title.str()]->Write(); gm_ts->cd(); gm_tsx_mono->cd(); title.str(""); title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGM_X_ts_mono[title.str()]->Write(); gm_tsy_mono->cd(); title.str(""); title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGM_Y_ts_mono[title.str()]->Write(); gm_tsz_mono->cd(); title.str(""); title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts_mono"; hPullGM_Z_ts_mono[title.str()]->Write(); gp_tr->cd(); gp_trx_mono->cd(); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr_mono"; hPullGP_X_tr_mono[title.str()]->Write(); gp_try_mono->cd(); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr_mono"; hPullGP_Y_tr_mono[title.str()]->Write(); gp_trz_mono->cd(); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr_mono"; hPullGP_Z_tr_mono[title.str()]->Write(); gp_rs->cd(); gp_rsx_mono->cd(); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs_mono"; hPullGP_X_rs_mono[title.str()]->Write(); gp_rsy_mono->cd(); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs_mono"; hPullGP_Y_rs_mono[title.str()]->Write(); gp_rsz_mono->cd(); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs_mono"; hPullGP_Z_rs_mono[title.str()]->Write(); //stereo gp_ts->cd(); gp_tsx_stereo->cd(); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGP_X_ts_stereo[title.str()]->Write(); gp_tsy_stereo->cd(); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGP_Y_ts_stereo[title.str()]->Write(); gp_tsz_stereo->cd(); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGP_Z_ts_stereo[title.str()]->Write(); gm_ts->cd(); gm_tsx_stereo->cd(); title.str(""); title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGM_X_ts_stereo[title.str()]->Write(); gm_tsy_stereo->cd(); title.str(""); title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGM_Y_ts_stereo[title.str()]->Write(); gm_tsz_stereo->cd(); title.str(""); title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts_stereo"; hPullGM_Z_ts_stereo[title.str()]->Write(); gp_tr->cd(); gp_trx_stereo->cd(); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr_stereo"; hPullGP_X_tr_stereo[title.str()]->Write(); gp_try_stereo->cd(); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr_stereo"; hPullGP_Y_tr_stereo[title.str()]->Write(); gp_trz_stereo->cd(); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr_stereo"; hPullGP_Z_tr_stereo[title.str()]->Write(); gp_rs->cd(); gp_rsx_stereo->cd(); title.str(""); title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs_stereo"; hPullGP_X_rs_stereo[title.str()]->Write(); gp_rsy_stereo->cd(); title.str(""); title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs_stereo"; hPullGP_Y_rs_stereo[title.str()]->Write(); gp_rsz_stereo->cd(); title.str(""); title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs_stereo"; hPullGP_Z_rs_stereo[title.str()]->Write(); } } file->Close(); }
std::pair< LocalPoint, LocalVector > TestSmoothHits::projectHit | ( | const PSimHit & | hit, |
const StripGeomDetUnit * | stripDet, | ||
const BoundPlane & | plane | ||
) | [private] |
Definition at line 884 of file TestSmoothHits.cc.
References funct::cos(), dir, PSimHit::localDirection(), PSimHit::localPosition(), pileupReCalc_HLTpaths::scale, funct::sin(), StripGeomDetUnit::specificTopology(), StripTopology::strip(), StripTopology::stripAngle(), GeomDet::surface(), Surface::toGlobal(), GloballyPositioned< T >::toLocal(), and PV3DBase< T, PVType, FrameType >::z().
Referenced by analyze().
{ const StripTopology& topol = stripDet->specificTopology(); GlobalPoint globalpos= stripDet->surface().toGlobal(hit.localPosition()); LocalPoint localHit = plane.toLocal(globalpos); //track direction LocalVector locdir=hit.localDirection(); //rotate track in new frame GlobalVector globaldir= stripDet->surface().toGlobal(locdir); LocalVector dir=plane.toLocal(globaldir); float scale = -localHit.z() / dir.z(); LocalPoint projectedPos = localHit + scale*dir; float selfAngle = topol.stripAngle( topol.strip( hit.localPosition())); LocalVector stripDir( sin(selfAngle), cos(selfAngle), 0); // vector along strip in hit frame LocalVector localStripDir( plane.toLocal(stripDet->surface().toGlobal( stripDir))); return std::pair<LocalPoint,LocalVector>( projectedPos, localStripDir); }
std::string TestSmoothHits::builderName [private] |
Definition at line 69 of file TestSmoothHits.h.
Referenced by beginRun(), and TestSmoothHits().
const edm::ParameterSet TestSmoothHits::conf_ [private] |
Definition at line 63 of file TestSmoothHits.h.
Referenced by TestSmoothHits().
TFile* TestSmoothHits::file [private] |
Definition at line 82 of file TestSmoothHits.h.
Referenced by beginRun(), and endJob().
edm::ESHandle<TrajectoryFitter> TestSmoothHits::fit [private] |
Definition at line 79 of file TestSmoothHits.h.
Referenced by analyze(), and beginRun().
std::string TestSmoothHits::fname [private] |
Definition at line 71 of file TestSmoothHits.h.
Referenced by beginRun(), and TestSmoothHits().
TH2F * TestSmoothHits::hChi2_vs_clsize [private] |
Definition at line 98 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
TH2F* TestSmoothHits::hChi2_vs_Process [private] |
Definition at line 98 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hChi2Increment [private] |
Definition at line 96 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
Definition at line 64 of file TestSmoothHits.h.
Referenced by analyze().
std::map<std::string,TH1F*> TestSmoothHits::hPullGM_X_ts [private] |
Definition at line 87 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGM_X_ts_mono [private] |
Definition at line 103 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGM_X_ts_stereo [private] |
Definition at line 116 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGM_Y_ts [private] |
Definition at line 88 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGM_Y_ts_mono [private] |
Definition at line 104 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGM_Y_ts_stereo [private] |
Definition at line 117 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGM_Z_ts [private] |
Definition at line 89 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGM_Z_ts_mono [private] |
Definition at line 105 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGM_Z_ts_stereo [private] |
Definition at line 118 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_X_rs [private] |
Definition at line 90 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_X_rs_mono [private] |
Definition at line 106 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_X_rs_stereo [private] |
Definition at line 119 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_X_tr [private] |
Definition at line 93 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_X_tr_mono [private] |
Definition at line 109 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_X_tr_stereo [private] |
Definition at line 122 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_X_ts [private] |
Definition at line 84 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_X_ts_mono [private] |
Definition at line 100 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_X_ts_stereo [private] |
Definition at line 113 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Y_rs [private] |
Definition at line 91 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Y_rs_mono [private] |
Definition at line 107 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Y_rs_stereo [private] |
Definition at line 120 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Y_tr [private] |
Definition at line 94 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Y_tr_mono [private] |
Definition at line 110 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Y_tr_stereo [private] |
Definition at line 123 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Y_ts [private] |
Definition at line 85 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Y_ts_mono [private] |
Definition at line 101 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Y_ts_stereo [private] |
Definition at line 114 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Z_rs [private] |
Definition at line 92 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Z_rs_mono [private] |
Definition at line 108 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Z_rs_stereo [private] |
Definition at line 121 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Z_tr [private] |
Definition at line 95 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Z_tr_mono [private] |
Definition at line 111 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Z_tr_stereo [private] |
Definition at line 124 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Z_ts [private] |
Definition at line 86 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Z_ts_mono [private] |
Definition at line 102 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
std::map<std::string,TH1F*> TestSmoothHits::hPullGP_Z_ts_stereo [private] |
Definition at line 115 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
TH1F* TestSmoothHits::hTotChi2Increment [private] |
Definition at line 97 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().
double TestSmoothHits::maxeta [private] |
Definition at line 66 of file TestSmoothHits.h.
Referenced by analyze(), and TestSmoothHits().
double TestSmoothHits::mineta [private] |
Definition at line 66 of file TestSmoothHits.h.
Referenced by analyze(), and TestSmoothHits().
std::string TestSmoothHits::propagatorName [private] |
Definition at line 68 of file TestSmoothHits.h.
Referenced by beginRun(), and TestSmoothHits().
Definition at line 80 of file TestSmoothHits.h.
Referenced by analyze(), and beginRun().
std::string TestSmoothHits::sname [private] |
Definition at line 72 of file TestSmoothHits.h.
Referenced by beginRun(), and TestSmoothHits().
std::string TestSmoothHits::srcName [private] |
Definition at line 70 of file TestSmoothHits.h.
Referenced by analyze(), and TestSmoothHits().
Definition at line 77 of file TestSmoothHits.h.
Referenced by analyze(), and beginRun().
edm::ESHandle<TrackerGeometry> TestSmoothHits::theG [private] |
Definition at line 74 of file TestSmoothHits.h.
Referenced by analyze(), and beginRun().
edm::ESHandle<MagneticField> TestSmoothHits::theMF [private] |
Definition at line 75 of file TestSmoothHits.h.
Referenced by analyze(), and beginRun().
edm::ESHandle<Propagator> TestSmoothHits::thePropagator [private] |
Definition at line 76 of file TestSmoothHits.h.
Referenced by analyze(), and beginRun().
Definition at line 78 of file TestSmoothHits.h.
Referenced by analyze().
std::stringstream TestSmoothHits::title [private] |
Definition at line 83 of file TestSmoothHits.h.
Referenced by analyze(), beginRun(), and endJob().