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 &ev, edm::EventSetup const &es) final

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

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

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

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

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

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

	~DQMEDAnalyzer () override=default

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 () 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

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)

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 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::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 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 34 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 94 of file TrackingRecoMaterialAnalyser.cc.

95 {

96 }

Member Function Documentation

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

override

Definition at line 189 of file TrackingRecoMaterialAnalyser.cc.

References funct::abs(), ecalDrivenElectronSeedsParameters_cff::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(), 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(), refitter_, sDETS, TrackTransformer::setServices(), mathSSE::sqrt(), DetId::subdetId(), TrajectoryStateOnSurface::surface(), HiIsolationCommonParameters_cff::track, l1t::tracks, tracksToken_, TrackTransformer::transform(), usePV_, findQualityFiles::v, TrackingRecHit::valid, electrons_cff::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 99 of file TrackingRecoMaterialAnalyser.cc.

References DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), 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, DQMStore::IBooker::setCurrentFolder(), and fftjetcommon_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