Inheritance diagram for TrackingRecoMaterialAnalyser:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override

	TrackingRecoMaterialAnalyser (const edm::ParameterSet &)

	~TrackingRecoMaterialAnalyser () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
bool	isDoubleSided (SiStripDetId strip_id)

Private Attributes
const edm::EDGetTokenT< reco::BeamSpot >	beamspotToken_

MonitorElement *	deltaP_in_out_vs_eta_

MonitorElement *	deltaP_in_out_vs_eta_2d_

MonitorElement *	deltaP_in_out_vs_eta_vs_phi_2d_

MonitorElement *	deltaP_in_out_vs_z_

MonitorElement *	deltaP_in_out_vs_z_2d_

MonitorElement *	deltaPl_in_out_vs_eta_

MonitorElement *	deltaPl_in_out_vs_z_

MonitorElement *	deltaPt_in_out_2d_

MonitorElement *	deltaPt_in_out_vs_eta_

MonitorElement *	deltaPt_in_out_vs_z_

std::string	folder_

MonitorElement *	histo_RZ_

MonitorElement *	histo_RZ_Ori_

std::unordered_map< std::string, MonitorElement * >	histosEta_

std::unordered_map< std::string, MonitorElement * >	histosOriEta_

MonitorElement *	P_vs_eta_2d_

TrackTransformer	refitter_

const edm::EDGetTokenT< reco::TrackCollection >	tracksToken_

bool	usePV_

const edm::EDGetTokenT< reco::VertexCollection >	verticesToken_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 35 of file TrackingRecoMaterialAnalyser.cc.

Constructor & Destructor Documentation

TrackingRecoMaterialAnalyser::TrackingRecoMaterialAnalyser ( const edm::ParameterSet & iPSet )

explicit

Definition at line 70 of file TrackingRecoMaterialAnalyser.cc.

     : refitter_(iPSet),
       tracksToken_(consumes<reco::TrackCollection>(iPSet.getParameter<edm::InputTag>("tracks"))),
       beamspotToken_(consumes<reco::BeamSpot>(iPSet.getParameter<edm::InputTag>("beamspot"))),
       verticesToken_(mayConsume<reco::VertexCollection>(iPSet.getParameter<edm::InputTag>("vertices"))),
       usePV_(iPSet.getParameter<bool>("usePV")),
       folder_(iPSet.getParameter<std::string>("folder")),
       histo_RZ_(nullptr),
       histo_RZ_Ori_(nullptr),
       deltaPt_in_out_2d_(nullptr),
       deltaP_in_out_vs_eta_(nullptr),
       deltaP_in_out_vs_z_(nullptr),
       deltaP_in_out_vs_eta_2d_(nullptr),
       deltaP_in_out_vs_eta_vs_phi_2d_(nullptr),
       deltaP_in_out_vs_z_2d_(nullptr),
       deltaPt_in_out_vs_eta_(nullptr),
       deltaPt_in_out_vs_z_(nullptr),
       deltaPl_in_out_vs_eta_(nullptr),
       deltaPl_in_out_vs_z_(nullptr),
       P_vs_eta_2d_(nullptr) {}

TrackingRecoMaterialAnalyser::~TrackingRecoMaterialAnalyser ( void )

override

Definition at line 92 of file TrackingRecoMaterialAnalyser.cc.

92 {}

Member Function Documentation

void TrackingRecoMaterialAnalyser::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

override

Definition at line 173 of file TrackingRecoMaterialAnalyser.cc.

References funct::abs(), pwdgSkimBPark_cfi::beamSpot, beamspotToken_, DEFINE_FWK_MODULE, deltaP_in_out_vs_eta_, deltaP_in_out_vs_eta_2d_, deltaP_in_out_vs_eta_vs_phi_2d_, deltaP_in_out_vs_z_, deltaP_in_out_vs_z_2d_, deltaPl_in_out_vs_eta_, deltaPl_in_out_vs_z_, deltaPt_in_out_2d_, deltaPt_in_out_vs_eta_, deltaPt_in_out_vs_z_, reco::TrackBase::dxy(), PV3DBase< T, PVType, FrameType >::eta(), reco::TrackBase::eta(), dqm::impl::MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::EventSetup::get(), TrajectoryStateOnSurface::globalPosition(), histo_RZ_, histo_RZ_Ori_, histosEta_, histosOriEta_, reco::TrackBase::hitPattern(), SurfaceOrientation::inner, reco::Track::innerMomentum(), isDoubleSided(), reco::Vertex::isFake(), MediumProperties::isValid(), TrajectoryStateOnSurface::isValid(), edm::HandleBase::isValid(), TrackerTopology::layer(), TrajectoryStateOnSurface::localMomentum(), LogTrace, Surface::mediumProperties(), TrackingRecHit::missing, reco::Vertex::ndof(), reco::HitPattern::numberOfLostTrackerHits(), SurfaceOrientation::outer, reco::Track::outerMomentum(), reco::Track::outerPosition(), reco::Track::outerZ(), reco::TrackBase::p(), p2, P_vs_eta_2d_, PV3DBase< T, PVType, FrameType >::perp(), reco::BeamSpot::position(), reco::Vertex::position(), MetAnalyzer::pv(), reco::TrackBase::quality(), reco::TrackBase::qualityByName(), MediumProperties::radLen(), fastSimProducer_cff::radLen, refitter_, sDETS, TrackTransformer::setServices(), mathSSE::sqrt(), DetId::subdetId(), TrajectoryStateOnSurface::surface(), HLT_2018_cff::track, PDWG_EXOHSCP_cff::tracks, tracksToken_, TrackTransformer::transform(), usePV_, findQualityFiles::v, TrackingRecHit::valid, pwdgSkimBPark_cfi::vertices, verticesToken_, hybridSuperClusters_cfi::xi, MediumProperties::xi(), PV3DBase< T, PVType, FrameType >::z(), and reco::Vertex::z().

                                                                                             {
   using namespace edm;
   using namespace reco;
   using namespace std;
 
   refitter_.setServices(setup);
 
   Handle<TrackCollection> tracks;
   Handle<VertexCollection> vertices;
   ESHandle<TrackerTopology> trk_topology;
 
   // Get the TrackerTopology
   setup.get<TrackerTopologyRcd>().get(trk_topology);
 
   // Get Tracks
   event.getByToken(tracksToken_, tracks);
   if (!tracks.isValid() || tracks->empty()) {
     LogInfo("TrackingRecoMaterialAnalyser") << "Invalid or empty track collection" << endl;
     return;
   }
 
   // Track Selector
   auto selector = [](const Track &track, const reco::Vertex::Point &pv) -> bool {
     return (track.quality(track.qualityByName("highPurity")) && track.dxy(pv) < 0.01 &&
             track.hitPattern().numberOfLostTrackerHits(HitPattern::MISSING_OUTER_HITS) == 0);
   };
 
   // Get BeamSpot
   Handle<BeamSpot> beamSpot;
   event.getByToken(beamspotToken_, beamSpot);
   // Bail out if missing
   if (!beamSpot.isValid())
     return;
 
   reco::Vertex::Point pv(beamSpot->position());
   if (usePV_) {
     event.getByToken(verticesToken_, vertices);
     if (vertices.isValid() && !vertices->empty()) {
       // Since we need to use eta and Z information from the tracks, in
       // order not to have the reco material distribution washed out due
       // to geometrical effects, we need to confine the PV to some small
       // region.
       const Vertex &v = (*vertices)[0];
       if (!v.isFake() && v.ndof() > 4 && std::fabs(v.z()) < 24 && std::fabs(v.position().rho()) < 2)
         pv = v.position();
     }
   }
 
   // Main idea:
   // * select first good tracks in input, according to reasonable criteria
   // * refit the tracks so that we have access to the TrajectoryMeasurements
   //   that internally have all the information about the TSOS on all
   //   crossed layers. We need the refit track and not the original one so
   //   that we are able to correctly compute the path travelled by the track
   //   within the detector, using its updated TSOS. The material description
   //   can in principle be derived also directly from the rechits, via the
   //   det()[(->)GeomDet *]->mediumProperties chain, but that would simply give the
   //   face values, not the "real" ones used while propagating the track.
   // * Loop on all measurements, extract the information about the TSOS,
   //   its surface and its mediumProperties
   // * Make plots for the untouched material properties, but also for the
   //   ones corrected by the track direction, since the material properties,
   //   according to the documentation, should refer to normal incidence of
   //   the track, which is seldom the case, according to the current direction
   TrajectoryStateOnSurface current_tsos;
   DetId current_det;
   for (auto const track : *tracks) {
     if (!selector(track, pv))
       continue;
     auto const &inner = track.innerMomentum();
     auto const &outer = track.outerMomentum();
     deltaP_in_out_vs_eta_->Fill(inner.eta(), inner.R() - outer.R());
     deltaP_in_out_vs_z_->Fill(track.outerZ(), inner.R() - outer.R());
     deltaP_in_out_vs_eta_2d_->Fill(inner.eta(), inner.R() - outer.R());
     deltaP_in_out_vs_eta_vs_phi_2d_->Fill(inner.eta(), inner.phi(), inner.R() - outer.R());
     deltaP_in_out_vs_z_2d_->Fill(track.outerZ(), inner.R() - outer.R());
     deltaPt_in_out_vs_eta_->Fill(inner.eta(), inner.rho() - outer.rho());
     deltaPt_in_out_vs_z_->Fill(track.outerZ(), inner.rho() - outer.rho());
     deltaPl_in_out_vs_eta_->Fill(inner.eta(), inner.z() - outer.z());
     deltaPl_in_out_vs_z_->Fill(track.outerZ(), inner.z() - outer.z());
     deltaPt_in_out_2d_->Fill(track.outerZ(), track.outerPosition().rho(), inner.rho() - outer.rho());
     P_vs_eta_2d_->Fill(track.eta(), track.p());
     vector<Trajectory> traj = refitter_.transform(track);
     if (traj.size() > 1 || traj.empty())
       continue;
     for (auto const &tm : traj.front().measurements()) {
       if (tm.recHit().get() &&
           (tm.recHitR().type() == TrackingRecHit::valid || tm.recHitR().type() == TrackingRecHit::missing)) {
         current_tsos = tm.updatedState().isValid() ? tm.updatedState() : tm.forwardPredictedState();
         auto const &localP = current_tsos.localMomentum();
         current_det = tm.recHit()->geographicalId();
         const Surface &surface = current_tsos.surface();
         assert(tm.recHit()->surface() == &surface);
         if (!surface.mediumProperties().isValid()) {
           LogError("TrackingRecoMaterialAnalyser")
               << "Medium properties for material linked to detector"
               << " are invalid at: " << current_tsos.globalPosition() << " " << (SiStripDetId)current_det << endl;
           assert(0);
           continue;
         }
         float p2 = localP.mag2();
         float xf = std::abs(std::sqrt(p2) / localP.z());
         float ori_xi = surface.mediumProperties().xi();
         float ori_radLen = surface.mediumProperties().radLen();
         float xi = ori_xi * xf;
         float radLen = ori_radLen * xf;
 
         // NOTA BENE: THIS ASSUMES THAT THE TRACKS HAVE BEEN PRODUCED
         // WITH SPLITTING, AS IS THE CASE FOR THE generalTracks
         // collection.
 
         // Since there are double-sided (glued) modules all over the
         // tracker, the material budget has been internally
         // partitioned in two equal components, so that each single
         // layer will receive half of the correct radLen. For this
         // reason, only for the double-sided components, we rescale
         // the obtained radLen by 2.
 
         // In particular see code here:
         // http://cmslxr.fnal.gov/dxr/CMSSW_8_0_5/source/Geometry/TrackerGeometryBuilder/src/TrackerGeomBuilderFromGeometricDet.cc#213
         // where, in the SiStrip Tracker, if the module has a partner
         // (i.e. it's a glued detector) the plane is built with a
         // scaling of 0.5. The actual plane is built few lines below:
         // http://cmslxr.fnal.gov/dxr/CMSSW_8_0_5/source/Geometry/TrackerGeometryBuilder/src/TrackerGeomBuilderFromGeometricDet.cc#287
 
         if (isDoubleSided(current_det)) {
           LogTrace("TrackingRecoMaterialAnalyser")
               << "Eta: " << track.eta() << " " << sDETS[current_det.subdetId()] << trk_topology->layer(current_det)
               << " has ori_radLen: " << ori_radLen << " and ori_xi: " << xi << " and has radLen: " << radLen
               << "  and xi: " << xi << endl;
           ori_radLen *= 2.;
           radLen *= 2.;
         }
 
         histosOriEta_[sDETS[current_det.subdetId()] + to_string(trk_topology->layer(current_det))]->Fill(
             current_tsos.globalPosition().eta(), ori_radLen);
         histosEta_[sDETS[current_det.subdetId()] + to_string(trk_topology->layer(current_det))]->Fill(
             current_tsos.globalPosition().eta(), radLen);
         histo_RZ_Ori_->Fill(current_tsos.globalPosition().z(), current_tsos.globalPosition().perp(), ori_radLen);
         histo_RZ_->Fill(current_tsos.globalPosition().z(), current_tsos.globalPosition().perp(), radLen);
         LogInfo("TrackingRecoMaterialAnalyser")
             << "Eta: " << track.eta() << " " << sDETS[current_det.subdetId()] << trk_topology->layer(current_det)
             << " has ori_radLen: " << ori_radLen << " and ori_xi: " << xi << " and has radLen: " << radLen
             << "  and xi: " << xi << endl;
       }
     }
   }
 }

void TrackingRecoMaterialAnalyser::bookHistograms	(	DQMStore::IBooker &	i,
		edm::Run const &	,
		edm::EventSetup const &	setup
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 95 of file TrackingRecoMaterialAnalyser.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile2D(), data, deltaP_in_out_vs_eta_, deltaP_in_out_vs_eta_2d_, deltaP_in_out_vs_eta_vs_phi_2d_, deltaP_in_out_vs_z_, deltaP_in_out_vs_z_2d_, deltaPl_in_out_vs_eta_, deltaPl_in_out_vs_z_, deltaPt_in_out_2d_, deltaPt_in_out_vs_eta_, deltaPt_in_out_vs_z_, folder_, edm::EventSetup::get(), histo_RZ_, histo_RZ_Ori_, histosEta_, histosOriEta_, crabWrapper::key, TrackerGeometry::numberOfLayers(), P_vs_eta_2d_, sDETS, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), and overlapproblemtsosanalyzer_cfi::title.

                                                                               {
   using namespace std;
   edm::ESHandle<TrackerGeometry> trackerGeometry;
   setup.get<TrackerDigiGeometryRecord>().get(trackerGeometry);
 
   ibook.setCurrentFolder(folder_);
 
   // Histogram to store the radiation length map, in the R-Z plane,
   // gathering numbers directly from the trackerRecoMaterial.xml
   // file. The numbers are not corrected for the track angle.
   histo_RZ_Ori_ = ibook.bookProfile2D("OriRadLen", "Original_RadLen", 600, -300., 300, 120, 0., 120., 0., 1.);
 
   // Histogram to store the radiation length map, as before, but
   // correcting the numbers with the track angle. This represents the
   // real material seen by the track.
   histo_RZ_ = ibook.bookProfile2D("RadLen", "RadLen", 600, -300., 300, 120, 0., 120., 0., 1.);
 
   // Histogram to show the deltaP Out-In in the eta-phi plane.
   deltaP_in_out_vs_eta_vs_phi_2d_ = ibook.bookProfile2D(
       "DeltaP_in_out_vs_eta_vs_phi_2d", "DeltaP_in_out_vs_eta_vs_phi_2d", 100, -3.0, 3.0, 100, -3.15, 3.15, 0., 100.);
 
   // Histogram to show the deltaP Out-In vs eta.
   deltaP_in_out_vs_eta_2d_ =
       ibook.book2D("DeltaP_in_out_vs_eta_2d", "DeltaP_in_out_vs_eta_2d", 100, -3.0, 3.0, 100, 0., 1);
 
   // Histogram to show the deltaP Out-In vs Z. The Z coordinate is the
   // one computed at the outermost hit.
   deltaP_in_out_vs_z_2d_ = ibook.book2D("DeltaP_in_out_vs_z_2d", "DeltaP_in_out_vs_z_2d", 600, -300, 300, 200., -1, 1.);
 
   // Histogram to show the deltaP Out-In vs eta. The eta is the one
   // computed at the innermost hit.
   deltaP_in_out_vs_eta_ =
       ibook.bookProfile("DeltaP_in_out_vs_eta", "DeltaP_in_out_vs_eta", 100, -3.0, 3.0, -100., 100.);
   deltaP_in_out_vs_z_ = ibook.bookProfile("DeltaP_in_out_vs_z", "DeltaP_in_out_vs_z", 600, -300, 300, -100., 100.);
 
   // Histogram to show the delta_Pt Out-In vs eta. The eta is the one
   // computed at the innermost hit.
   deltaPt_in_out_vs_eta_ =
       ibook.bookProfile("DeltaPt_in_out_vs_eta", "DeltaPt_in_out_vs_eta", 100, -3.0, 3.0, -100., 100.);
 
   // Histogram to show the delta_Pt Out-In vs Z. The Z is the one
   // computed at the outermost hit.
   deltaPt_in_out_vs_z_ = ibook.bookProfile("DeltaPt_in_out_vs_z", "DeltaPt_in_out_vs_z", 600, -300, 300, -100., 100);
 
   // Histogram to show the delta_Pl Out-In vs eta. The eta is the one
   // computed at the innermost hit.
   deltaPl_in_out_vs_eta_ =
       ibook.bookProfile("DeltaPz_in_out_vs_eta", "DeltaPz_in_out_vs_eta", 100, -3.0, 3.0, -100., 100.);
 
   // Histogram to show the delta_Pl Out-In vs Z. The Z is the one
   // computed at the outermost hit.
   deltaPl_in_out_vs_z_ = ibook.bookProfile("DeltaPz_in_out_vs_z", "DeltaPz_in_out_vs_z", 600, -300, 300, -100., 100.);
 
   // Histogram to show the delta_Pt Out-In in the Z-R plane. Z and R
   // are related to the outermost hit.
   deltaPt_in_out_2d_ = ibook.bookProfile2D("DeltaPt 2D", "DeltaPt 2D", 600, -300., 300, 120, 0., 120., -100., 100.);
 
   // Histogram to show the distribution of p vs eta for all tracks.
   P_vs_eta_2d_ = ibook.book2D("P_vs_eta_2d", "P_vs_eta_2d", 100, -3.0, 3.0, 100., 0., 5.);
   char title[50];
   char key[20];
   for (unsigned int det = 1; det < sDETS.size(); ++det) {
     for (unsigned int sub_det = 1; sub_det <= trackerGeometry->numberOfLayers(det); ++sub_det) {
       memset(title, 0, sizeof(title));
       snprintf(title, sizeof(title), "Original_RadLen_vs_Eta_%s%d", sDETS[det].data(), sub_det);
       snprintf(key, sizeof(key), "%s%d", sDETS[det].data(), sub_det);
       histosOriEta_.insert(
           make_pair<string, MonitorElement *>(key, ibook.bookProfile(title, title, 250, -5.0, 5.0, 0., 1.)));
       snprintf(title, sizeof(title), "RadLen_vs_Eta_%s%d", sDETS[det].data(), sub_det);
       histosEta_.insert(
           make_pair<string, MonitorElement *>(key, ibook.bookProfile(title, title, 250, -5.0, 5.0, 0., 1.)));
     }
   }
 }

bool TrackingRecoMaterialAnalyser::isDoubleSided ( SiStripDetId strip_id )

inlineprivate

Definition at line 43 of file TrackingRecoMaterialAnalyser.cc.

References SiStripDetId::glued(), SiStripDetId::subDetector(), and SiStripDetId::UNKNOWN.

Referenced by analyze().

                                                    {
     return ((strip_id.subDetector() != SiStripDetId::UNKNOWN) && strip_id.glued());
   }

Member Data Documentation

const edm::EDGetTokenT<reco::BeamSpot> TrackingRecoMaterialAnalyser::beamspotToken_

private

Definition at line 48 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze().

MonitorElement* TrackingRecoMaterialAnalyser::deltaP_in_out_vs_eta_

private

Definition at line 57 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::deltaP_in_out_vs_eta_2d_

private

Definition at line 59 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::deltaP_in_out_vs_eta_vs_phi_2d_

private

Definition at line 60 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::deltaP_in_out_vs_z_

private

Definition at line 58 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::deltaP_in_out_vs_z_2d_

private

Definition at line 61 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::deltaPl_in_out_vs_eta_

private

Definition at line 64 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::deltaPl_in_out_vs_z_

private

Definition at line 65 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::deltaPt_in_out_2d_

private

Definition at line 56 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::deltaPt_in_out_vs_eta_

private

Definition at line 62 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::deltaPt_in_out_vs_z_

private

Definition at line 63 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

std::string TrackingRecoMaterialAnalyser::folder_

private

Definition at line 51 of file TrackingRecoMaterialAnalyser.cc.

Referenced by bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::histo_RZ_

private

Definition at line 54 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::histo_RZ_Ori_

private

Definition at line 55 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

std::unordered_map<std::string, MonitorElement *> TrackingRecoMaterialAnalyser::histosEta_

private

Definition at line 53 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

std::unordered_map<std::string, MonitorElement *> TrackingRecoMaterialAnalyser::histosOriEta_

private

Definition at line 52 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackingRecoMaterialAnalyser::P_vs_eta_2d_

private

Definition at line 66 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze(), and bookHistograms().

TrackTransformer TrackingRecoMaterialAnalyser::refitter_

private

Definition at line 46 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze().

const edm::EDGetTokenT<reco::TrackCollection> TrackingRecoMaterialAnalyser::tracksToken_

private

Definition at line 47 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze().

bool TrackingRecoMaterialAnalyser::usePV_

private

Definition at line 50 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze().

const edm::EDGetTokenT<reco::VertexCollection> TrackingRecoMaterialAnalyser::verticesToken_

private

Definition at line 49 of file TrackingRecoMaterialAnalyser.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation