TrackerSeedValidator Class Reference

#include <TrackerSeedValidator.h>

Inheritance diagram for TrackerSeedValidator:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	Method called once per event. More...
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	Method called to book the DQM histograms. More...
	TrackerSeedValidator (const edm::ParameterSet &pset)
	Constructor. More...
virtual	~TrackerSeedValidator ()
	Destructor. More...
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
std::vector< edm::EDGetTokenT < reco::TrackToTrackingParticleAssociator > >	associatorTokens
std::string	builderName
CosmicTrackingParticleSelector	cosmictpSelector
std::string	dirName_
std::unique_ptr < MTVHistoProducerAlgoForTracker >	histoProducerAlgo_
edm::ESHandle < TransientTrackingRecHitBuilder >	theTTRHBuilder
TrackingParticleSelector	tpSelector

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
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 Member Functions inherited from MultiTrackValidatorBase
	MultiTrackValidatorBase (const edm::ParameterSet &pset, edm::ConsumesCollector &&iC, bool isSeed=false)
	Constructor. More...
virtual	~MultiTrackValidatorBase ()
	Destructor. More...
Protected Attributes inherited from MultiTrackValidatorBase
std::vector< edm::InputTag >	associators
edm::EDGetTokenT< reco::BeamSpot >	bsSrc
bool	ignoremissingtkcollection_
std::vector< edm::InputTag >	label
edm::EDGetTokenT< std::vector < PileupSummaryInfo > >	label_pileupinfo
edm::EDGetTokenT < TrackingParticleCollection >	label_tp_effic
edm::EDGetTokenT < TrackingParticleCollection >	label_tp_fake
edm::EDGetTokenT < TrackingVertexCollection >	label_tv
std::vector< edm::EDGetTokenT < edm::View< reco::Track > > >	labelToken
std::vector< edm::EDGetTokenT < edm::View< TrajectorySeed > > >	labelTokenSeed
edm::EDGetTokenT < edm::ValueMap < reco::DeDxData > >	m_dEdx1Tag
edm::EDGetTokenT < edm::ValueMap < reco::DeDxData > >	m_dEdx2Tag
std::string	parametersDefiner
std::vector< edm::EDGetTokenT < std::vector< PSimHit > > >	simHitTokens_

Detailed Description

Class that prodecs histrograms to validate Track Reconstruction performances

Author

cerati

Definition at line 18 of file TrackerSeedValidator.h.

Constructor & Destructor Documentation

TrackerSeedValidator::TrackerSeedValidator

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 38 of file TrackerSeedValidator.cc.

References MultiTrackValidatorBase::associators, associatorTokens, builderName, dirName_, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and tpSelector.

                                                                     :
  MultiTrackValidatorBase(pset, consumesCollector(),true),
  histoProducerAlgo_(std::make_unique<MTVHistoProducerAlgoForTracker>(pset.getParameter<ParameterSet>("histoProducerAlgoBlock"), consumesCollector())) {
  dirName_ = pset.getParameter<std::string>("dirName");

  tpSelector = TrackingParticleSelector(pset.getParameter<double>("ptMinTP"),
                    pset.getParameter<double>("minRapidityTP"),
                    pset.getParameter<double>("maxRapidityTP"),
                    pset.getParameter<double>("tipTP"),
                    pset.getParameter<double>("lipTP"),
                    pset.getParameter<int>("minHitTP"),
                    pset.getParameter<bool>("signalOnlyTP"),
                    pset.getParameter<bool>("intimeOnlyTP"),
                    pset.getParameter<bool>("chargedOnlyTP"),
                    pset.getParameter<bool>("stableOnlyTP"),
                    pset.getParameter<std::vector<int> >("pdgIdTP"));

  builderName = pset.getParameter<std::string>("TTRHBuilder");

  for (auto const& associator: associators) {
    associatorTokens.push_back(consumes<reco::TrackToTrackingParticleAssociator>(associator));
  }
}

TrackerSeedValidator::~TrackerSeedValidator ( )

virtual

Destructor.

Definition at line 62 of file TrackerSeedValidator.cc.

{}

Member Function Documentation

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

overridevirtual

Method called once per event.

Implements edm::stream::EDAnalyzerBase.

Definition at line 105 of file TrackerSeedValidator.cc.

References reco::TrackBase::appendHits(), reco::TrackToTrackingParticleAssociator::associateRecoToSim(), reco::TrackToTrackingParticleAssociator::associateSimToReco(), MultiTrackValidatorBase::associators, associatorTokens, MultiTrackValidatorBase::bsSrc, builderName, EncodedEventId::bunchCrossing(), TrackingParticleNumberOfLayers::calculate(), LocalTrajectoryParameters::charge(), funct::cos(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, alignCSCRings::e, edm::AssociationMap< Tag >::end(), EncodedEventId::event(), TrackingParticle::eventId(), edm::AssociationMap< Tag >::find(), plotBeamSpotDB::first, TrajectoryStateOnSurface::freeState(), PerigeeConversions::ftsToPerigeeError(), edm::EventSetup::get(), PileupSummaryInfo::getPU_NumInteractions(), histoProducerAlgo_, i, MultiTrackValidatorBase::label, MultiTrackValidatorBase::label_pileupinfo, MultiTrackValidatorBase::label_tp_effic, MultiTrackValidatorBase::label_tp_fake, MultiTrackValidatorBase::labelTokenSeed, LogTrace, PTrajectoryStateOnDet::parameters(), MultiTrackValidatorBase::parametersDefiner, reco::BeamSpot::position(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), TrajectorySeed::recHits(), fileCollector::seed, MultiTrackValidatorBase::simHitTokens_, funct::sin(), findQualityFiles::size, mathSSE::sqrt(), TrajectorySeed::startingState(), theTTRHBuilder, tpSelector, trajectoryStateTransform::transientState(), w, and cms::Exception::what().

                                                                                   {

  edm::LogInfo("TrackValidator") << "\n====================================================" << "\n"
                 << "Analyzing new event" << "\n"
                 << "====================================================\n" << "\n";

  edm::ESHandle<ParametersDefinerForTP> parametersDefinerTP;
  setup.get<TrackAssociatorRecord>().get(parametersDefiner,parametersDefinerTP);

  edm::ESHandle<TrackerTopology> httopo;
  setup.get<TrackerTopologyRcd>().get(httopo);
  const TrackerTopology& ttopo = *httopo;

  edm::ESHandle<TransientTrackingRecHitBuilder> theTTRHBuilder;
  setup.get<TransientRecHitRecord>().get(builderName,theTTRHBuilder);

  edm::ESHandle<MagneticField> theMF;
  setup.get<IdealMagneticFieldRecord>().get(theMF);

  edm::Handle<TrackingParticleCollection>  TPCollectionHeff ;
  event.getByToken(label_tp_effic,TPCollectionHeff);
  const TrackingParticleCollection tPCeff = *(TPCollectionHeff.product());

  edm::Handle<TrackingParticleCollection>  TPCollectionHfake ;
  event.getByToken(label_tp_fake,TPCollectionHfake);
  const TrackingParticleCollection tPCfake = *(TPCollectionHfake.product());

  if (tPCeff.size()==0) {edm::LogInfo("TrackValidator") << "TP Collection for efficiency studies has size = 0! Skipping Event." ; return;}
  if (tPCfake.size()==0) {edm::LogInfo("TrackValidator") << "TP Collection for fake rate studies has size = 0! Skipping Event." ; return;}

  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  event.getByToken(bsSrc,recoBeamSpotHandle);
  reco::BeamSpot bs = *recoBeamSpotHandle;

  edm::Handle< vector<PileupSummaryInfo> > puinfoH;
  event.getByToken(label_pileupinfo,puinfoH);
  PileupSummaryInfo puinfo;

  for (unsigned int puinfo_ite=0;puinfo_ite<(*puinfoH).size();++puinfo_ite){
    if ((*puinfoH)[puinfo_ite].getBunchCrossing()==0){
      puinfo=(*puinfoH)[puinfo_ite];
      break;
    }
  }

  // Calculate the number of 3D layers for TPs
  //
  // I would have preferred to produce the ValueMap to Event and read
  // it from there, but there would have been quite some number of
  // knock-on effects, and again the fact that we take two TP
  // collections do not support Ref<TP>'s would have complicated that.
  //
  // In principle we could use the SimHitTPAssociationList read above
  // for parametersDefinerIsCosmic_, but since we don't currently
  // support Ref<TP>s, we can't in general use it since eff/fake TP
  // collections can, in general, be different.
  TrackingParticleNumberOfLayers tpNumberOfLayersAlgo(event, simHitTokens_);
  auto nlayers_tPCeff_ptrs = tpNumberOfLayersAlgo.calculate(TPCollectionHeff, setup);
  const auto& nLayers_tPCeff = *(std::get<TrackingParticleNumberOfLayers::nTrackerLayers>(nlayers_tPCeff_ptrs));
  const auto& nPixelLayers_tPCeff = *(std::get<TrackingParticleNumberOfLayers::nPixelLayers>(nlayers_tPCeff_ptrs));
  const auto& nStripMonoAndStereoLayers_tPCeff = *(std::get<TrackingParticleNumberOfLayers::nStripMonoAndStereoLayers>(nlayers_tPCeff_ptrs));

  int w=0;
  for (unsigned int ww=0;ww<associators.size();ww++){
    edm::Handle<reco::TrackToTrackingParticleAssociator> theAssociator;
    event.getByToken(associatorTokens[ww], theAssociator);
    const reco::TrackToTrackingParticleAssociator *associator = theAssociator.product();

    for (unsigned int www=0;www<label.size();www++){
      edm::LogVerbatim("TrackValidator") << "Analyzing "
                     << label[www].process()<<":"
                     << label[www].label()<<":"
                     << label[www].instance()<<" with "
                     << associators[ww] <<"\n";
      //
      //get collections from the event
      //
      edm::Handle<edm::View<TrajectorySeed> > seedCollection;
      event.getByToken(labelTokenSeed[www], seedCollection);
      if (seedCollection->size()==0) {
    edm::LogInfo("TrackValidator") << "SeedCollection size = 0!" ;
    continue;
      }

      //associate seeds
      LogTrace("TrackValidator") << "Calling associateRecoToSim method" << "\n";
      reco::RecoToSimCollectionSeed recSimColl=associator->associateRecoToSim(seedCollection,
                                                                              TPCollectionHfake);
      LogTrace("TrackValidator") << "Calling associateSimToReco method" << "\n";
      reco::SimToRecoCollectionSeed simRecColl=associator->associateSimToReco(seedCollection,
                                                                              TPCollectionHeff);

      //
      //fill simulation histograms
      //compute number of seeds per eta interval
      //
      edm::LogVerbatim("TrackValidator") << "\n# of TrackingParticles: " << tPCeff.size() << "\n";
      int ats(0);     //This counter counts the number of simTracks that are "associated" to recoTracks
      int st(0);          //This counter counts the number of simulated tracks passing the MTV selection (i.e. tpSelector(tp) )
      unsigned sts(0);   //This counter counts the number of simTracks surviving the bunchcrossing cut
      unsigned asts(0);  //This counter counts the number of simTracks that are "associated" to recoTracks surviving the bunchcrossing cut
      for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
    TrackingParticleRef tp(TPCollectionHeff, i);

    if (tp->charge()==0) continue;

    if(! tpSelector(*tp)) continue;

    TrackingParticle::Vector momentumTP = tp->momentum();
    TrackingParticle::Point vertexTP = tp->vertex();
    //Calcualte the impact parameters w.r.t. PCA
    TrackingParticle::Vector momentum = parametersDefinerTP->momentum(event,setup,tp);
    TrackingParticle::Point vertex = parametersDefinerTP->vertex(event,setup,tp);
    double dxySim = (-vertex.x()*sin(momentum.phi())+vertex.y()*cos(momentum.phi()));
    double dzSim = vertex.z() - (vertex.x()*momentum.x()+vertex.y()*momentum.y())/sqrt(momentum.perp2())
      * momentum.z()/sqrt(momentum.perp2());

        if(tp->eventId().bunchCrossing() == 0) {
          st++;
        }

        if(w == 0)
          histoProducerAlgo_->fill_generic_simTrack_histos(momentumTP,vertexTP, tp->eventId().bunchCrossing());

    const TrajectorySeed* matchedSeedPointer=0;
    std::vector<std::pair<edm::RefToBase<TrajectorySeed>, double> > rt;
    if(simRecColl.find(tp) != simRecColl.end()){
      rt = simRecColl[tp];
      if (rt.size()!=0) {
        ats++;
        matchedSeedPointer = rt.begin()->first.get();
        edm::LogVerbatim("TrackValidator") << "TrackingParticle #" << st
                           << " with pt=" << sqrt(tp->momentum().perp2())
                           << " associated with quality:" << rt.begin()->second <<"\n";
      }
    }else{
      edm::LogVerbatim("TrackValidator") << "TrackingParticle #" << st
                         << " with pt=" << sqrt(tp->momentum().perp2())
                         << " NOT associated to any TrajectorySeed" << "\n";
    }

        int nSimHits = tp->numberOfTrackerHits();
        int nSimLayers = nLayers_tPCeff[tp];
        int nSimPixelLayers = nPixelLayers_tPCeff[tp];
        int nSimStripMonoAndStereoLayers = nStripMonoAndStereoLayers_tPCeff[tp];

    //fixme convert seed into track
    reco::Track* matchedTrackPointer = 0;
    if (matchedSeedPointer) {
      TSCBLBuilderNoMaterial tscblBuilder;
      TransientTrackingRecHit::RecHitPointer recHit = theTTRHBuilder->build(&*(matchedSeedPointer->recHits().second-1));
      TrajectoryStateOnSurface state = trajectoryStateTransform::transientState( matchedSeedPointer->startingState(), recHit->surface(), theMF.product());
      TrajectoryStateClosestToBeamLine tsAtClosestApproachSeed = tscblBuilder(*state.freeState(),bs);//as in TrackProducerAlgorithm
      if(!(tsAtClosestApproachSeed.isValid())){
        edm::LogVerbatim("SeedValidator")<<"TrajectoryStateClosestToBeamLine not valid";
        continue;
      }
      const reco::TrackBase::Point vSeed1(tsAtClosestApproachSeed.trackStateAtPCA().position().x(),
                          tsAtClosestApproachSeed.trackStateAtPCA().position().y(),
                          tsAtClosestApproachSeed.trackStateAtPCA().position().z());
      const reco::TrackBase::Vector pSeed(tsAtClosestApproachSeed.trackStateAtPCA().momentum().x(),
                          tsAtClosestApproachSeed.trackStateAtPCA().momentum().y(),
                          tsAtClosestApproachSeed.trackStateAtPCA().momentum().z());
      //GlobalPoint vSeed(vSeed1.x()-bs.x0(),vSeed1.y()-bs.y0(),vSeed1.z()-bs.z0());
      PerigeeTrajectoryError seedPerigeeErrors = PerigeeConversions::ftsToPerigeeError(tsAtClosestApproachSeed.trackStateAtPCA());
      matchedTrackPointer = new reco::Track(0.,0., vSeed1, pSeed, 1, seedPerigeeErrors.covarianceMatrix());
      matchedTrackPointer->appendHits(matchedSeedPointer->recHits().first,matchedSeedPointer->recHits().second, ttopo);
    }

    double dR=0;//fixme: plots vs dR not implemented for now
    histoProducerAlgo_->fill_recoAssociated_simTrack_histos(w,*tp,tp->momentum(),tp->vertex(),dxySim,dzSim,0,0,nSimHits,nSimLayers,nSimPixelLayers,nSimStripMonoAndStereoLayers,
                                matchedTrackPointer,puinfo.getPU_NumInteractions(),dR, nullptr);

    sts++;
    if (matchedTrackPointer) asts++;

      } // End  for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){

      if(w == 0)
        histoProducerAlgo_->fill_simTrackBased_histos(st);

      //
      //fill reconstructed seed histograms
      //
      edm::LogVerbatim("TrackValidator") << "\n# of TrajectorySeeds with "
                     << label[www].process()<<":"
                     << label[www].label()<<":"
                     << label[www].instance()
                     << ": " << seedCollection->size() << "\n";
      int sat(0); //This counter counts the number of recoTracks that are associated to SimTracks from Signal only
      int at(0); //This counter counts the number of recoTracks that are associated to SimTracks
      int rT(0); //This counter counts the number of recoTracks in general

      TSCBLBuilderNoMaterial tscblBuilder;
      for(TrajectorySeedCollection::size_type i=0; i<seedCollection->size(); ++i){
    edm::RefToBase<TrajectorySeed> seed(seedCollection, i);
    rT++;

    //get parameters and errors from the seed state
    TransientTrackingRecHit::RecHitPointer recHit = theTTRHBuilder->build(&*(seed->recHits().second-1));
    TrajectoryStateOnSurface state = trajectoryStateTransform::transientState( seed->startingState(), recHit->surface(), theMF.product());
    TrajectoryStateClosestToBeamLine tsAtClosestApproachSeed = tscblBuilder(*state.freeState(),bs);//as in TrackProducerAlgorithm
    if(!(tsAtClosestApproachSeed.isValid())){
      edm::LogVerbatim("SeedValidator")<<"TrajectoryStateClosestToBeamLine not valid";
      continue;
        }
    const reco::TrackBase::Point vSeed1(tsAtClosestApproachSeed.trackStateAtPCA().position().x(),
                        tsAtClosestApproachSeed.trackStateAtPCA().position().y(),
                        tsAtClosestApproachSeed.trackStateAtPCA().position().z());
    const reco::TrackBase::Vector pSeed(tsAtClosestApproachSeed.trackStateAtPCA().momentum().x(),
                        tsAtClosestApproachSeed.trackStateAtPCA().momentum().y(),
                        tsAtClosestApproachSeed.trackStateAtPCA().momentum().z());
    //GlobalPoint vSeed(vSeed1.x()-bs.x0(),vSeed1.y()-bs.y0(),vSeed1.z()-bs.z0());
    PerigeeTrajectoryError seedPerigeeErrors = PerigeeConversions::ftsToPerigeeError(tsAtClosestApproachSeed.trackStateAtPCA());

    //fixme
    reco::Track* trackFromSeed = new reco::Track(0.,0., vSeed1, pSeed, 1, seedPerigeeErrors.covarianceMatrix());
    trackFromSeed->appendHits(seed->recHits().first,seed->recHits().second, ttopo);

    bool isSigSimMatched(false);
    bool isSimMatched(false);
    bool isChargeMatched(true);
    int numAssocSeeds = 0;
    int nSimHits = 0;
    double sharedFraction = 0.;
    std::vector<std::pair<TrackingParticleRef, double> > tp;
    if(recSimColl.find(seed) != recSimColl.end()) {
      tp = recSimColl[seed];
      if (tp.size()!=0) {

            nSimHits = tp[0].first->numberOfTrackerHits();
            sharedFraction = tp[0].second;
        isSimMatched = true;
        if (tp[0].first->charge() != seed->startingState().parameters().charge()) isChargeMatched = false;
        if(simRecColl.find(tp[0].first) != simRecColl.end()) numAssocSeeds = simRecColl[tp[0].first].size();
        //std::cout << numAssocRecoTracks << std::endl;

        at++;

        for (unsigned int tp_ite=0;tp_ite<tp.size();++tp_ite){
              TrackingParticle trackpart = *(tp[tp_ite].first);
          if ((trackpart.eventId().event() == 0) && (trackpart.eventId().bunchCrossing() == 0)){
            isSigSimMatched = true;
        sat++;
        break;
          }
            }


        edm::LogVerbatim("SeedValidator") << "TrajectorySeed #" << rT << " associated with quality:" << tp.begin()->second <<"\n";
      }
    } else {
      edm::LogVerbatim("SeedValidator") << "TrajectorySeed #" << rT << " NOT associated to any TrackingParticle" << "\n";
    }

    double dR = 0.;//fixme: plots vs dR not implemented for now
    histoProducerAlgo_->fill_generic_recoTrack_histos(w,*trackFromSeed,bs.position(), nullptr, isSimMatched,isSigSimMatched,
                              isChargeMatched, numAssocSeeds, 
                              puinfo.getPU_NumInteractions(),
                              nSimHits, sharedFraction, dR);

    //Fill other histos
    try{
      if (tp.size()==0) continue;

      histoProducerAlgo_->fill_simAssociated_recoTrack_histos(w,*trackFromSeed);

      TrackingParticleRef tpr = tp.begin()->first;

      //compute tracking particle parameters at point of closest approach to the beamline
      TrackingParticle::Vector momentumTP = parametersDefinerTP->momentum(event,setup,tpr);
      TrackingParticle::Point vertexTP = parametersDefinerTP->vertex(event,setup,tpr);

      //      LogTrace("SeedValidatorTEST") << "assocChi2=" << tp.begin()->second << "\n"
      //                     << "" <<  "\n"
      //                     << "ptREC=" << ptSeed << "\n"
      //                     << "etaREC=" << etaSeed << "\n"
      //                     << "qoverpREC=" << qoverpSeed << "\n"
      //                     << "dxyREC=" << dxySeed << "\n"
      //                     << "dzREC=" << dzSeed << "\n"
      //                     << "thetaREC=" << thetaSeed << "\n"
      //                     << "phiREC=" << phiSeed << "\n"
      //                     << "" <<  "\n"
      //                     << "qoverpError()=" << qoverpErrorSeed << "\n"
      //                     << "dxyError()=" << dxyErrorSeed << "\n"
      //                     << "dzError()=" << dzErrorSeed << "\n"
      //                     << "thetaError()=" << lambdaErrorSeed << "\n"
      //                     << "phiError()=" << phiErrorSeed << "\n"
      //                     << "" <<  "\n"
      //                     << "ptSIM=" << sqrt(assocTrack->momentum().perp2()) << "\n"
      //                     << "etaSIM=" << assocTrack->momentum().Eta() << "\n"
      //                     << "qoverpSIM=" << qoverpSim << "\n"
      //                     << "dxySIM=" << dxySim << "\n"
      //                     << "dzSIM=" << dzSim << "\n"
      //                     << "thetaSIM=" << M_PI/2-lambdaSim << "\n"
      //                     << "phiSIM=" << phiSim << "\n"
      //                     << "" << "\n"
      //                     << "contrib_Qoverp=" << contrib_Qoverp << "\n"
      //                     << "contrib_dxy=" << contrib_dxy << "\n"
      //                     << "contrib_dz=" << contrib_dz << "\n"
      //                     << "contrib_theta=" << contrib_theta << "\n"
      //                     << "contrib_phi=" << contrib_phi << "\n"
      //                     << "" << "\n"
      //                     <<"chi2PULL="<<contrib_Qoverp+contrib_dxy+contrib_dz+contrib_theta+contrib_phi<<"\n";

      histoProducerAlgo_->fill_ResoAndPull_recoTrack_histos(w,momentumTP,vertexTP,tpr->charge(),
                                *trackFromSeed,bs.position());


    } catch (cms::Exception e){
      LogTrace("SeedValidator") << "exception found: " << e.what() << "\n";
    }
      }// End of for(TrajectorySeedCollection::size_type i=0; i<seedCollection->size(); ++i)

      histoProducerAlgo_->fill_trackBased_histos(w,at,rT,st);

      edm::LogVerbatim("SeedValidator") << "Total Simulated: " << st << "\n"
                     << "Total Associated (simToReco): " << ats << "\n"
                     << "Total Reconstructed: " << rT << "\n"
                     << "Total Associated (recoToSim): " << at << "\n"
                     << "Total Fakes: " << rT-at << "\n";
      w++;
    }
  }
}

void TrackerSeedValidator::bookHistograms ( DQMStore::IBooker &  ibook, edm::Run const &  , edm::EventSetup const &  setup  )

overridevirtual

Method called to book the DQM histograms.

Implements DQMEDAnalyzer.

Definition at line 64 of file TrackerSeedValidator.cc.

References ecalcalib_dqm_sourceclient-live_cfg::algo, MultiTrackValidatorBase::associators, DQMStore::IBooker::cd(), TrackerOfflineValidation_Dqm_cff::dirName, dirName_, histoProducerAlgo_, edm::InputTag::instance(), edm::InputTag::label(), MultiTrackValidatorBase::label, edm::InputTag::process(), python.rootplot.root2matplotlib::replace(), and DQMStore::IBooker::setCurrentFolder().

                                                                                                           {
  {
    ibook.cd();
    ibook.setCurrentFolder(dirName_ + "simulation");

    //Booking histograms concerning with simulated tracks
    histoProducerAlgo_->bookSimHistos(ibook);
  }

  for (unsigned int ww=0;ww<associators.size();ww++){
    for (unsigned int www=0;www<label.size();www++){

      ibook.cd();
      InputTag algo = label[www];
      string dirName=dirName_;
      if (algo.process()!="")
    dirName+=algo.process()+"_";
      if(algo.label()!="")
    dirName+=algo.label()+"_";
      if(algo.instance()!="")
    dirName+=algo.instance()+"_";
      //      if (dirName.find("Seeds")<dirName.length()){
      //    dirName.replace(dirName.find("Seeds"),6,"");
      //      }
      string assoc= associators[ww].label();;
      if (assoc.find("Track")<assoc.length()){
    assoc.replace(assoc.find("Track"),5,"");
      }
      dirName+=assoc;
      std::replace(dirName.begin(), dirName.end(), ':', '_');

      ibook.setCurrentFolder(dirName.c_str());

      //Booking histograms concerning with reconstructed tracks
      histoProducerAlgo_->bookSimTrackHistos(ibook);
      histoProducerAlgo_->bookRecoHistos(ibook);
    }//end loop www
  }// end loop ww
}

Member Data Documentation

std::vector<edm::EDGetTokenT<reco::TrackToTrackingParticleAssociator> > TrackerSeedValidator::associatorTokens

private

Definition at line 33 of file TrackerSeedValidator.h.

Referenced by analyze(), and TrackerSeedValidator().

std::string TrackerSeedValidator::builderName

private

Definition at line 35 of file TrackerSeedValidator.h.

Referenced by analyze(), and TrackerSeedValidator().

CosmicTrackingParticleSelector TrackerSeedValidator::cosmictpSelector

private

Definition at line 42 of file TrackerSeedValidator.h.

std::string TrackerSeedValidator::dirName_

private

Definition at line 37 of file TrackerSeedValidator.h.

Referenced by bookHistograms(), and TrackerSeedValidator().

std::unique_ptr<MTVHistoProducerAlgoForTracker> TrackerSeedValidator::histoProducerAlgo_

private

Definition at line 43 of file TrackerSeedValidator.h.

Referenced by analyze(), and bookHistograms().

edm::ESHandle<TransientTrackingRecHitBuilder> TrackerSeedValidator::theTTRHBuilder

private

Definition at line 36 of file TrackerSeedValidator.h.

Referenced by analyze().

TrackingParticleSelector TrackerSeedValidator::tpSelector

private

Definition at line 41 of file TrackerSeedValidator.h.

Referenced by analyze(), and TrackerSeedValidator().