TestOutliers Class Reference

#include <RecoTracker/DebugTools/plugins/TestOutliers.cc>

Inheritance diagram for TestOutliers:

Public Member Functions
	TestOutliers (const edm::ParameterSet &)
	~TestOutliers () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () 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
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginRun (edm::Run const &run, const edm::EventSetup &) override
void	endJob () override

Private Attributes
TH1F *	badcluster
TH1F *	badhittype
TH1F *	badlayer
TH1F *	badprocess
TH1F *	countOldT
TH1F *	countOutA
TH1F *	countOutT
TH1F *	deltahits
TH1F *	deltahitsAssocGained
TH1F *	deltahitsAssocLost
TH1F *	deltahitsNO
TH1F *	deltahitsOK
TH1F *	energyLoss
TH1F *	energyLossMax
TH1F *	energyLossRatio
TFile *	file
TH1F *	gainedhits
TH1F *	gainedhits2
TH1F *	goodbadhits
TH1F *	goodbadmerged
TH1F *	goodbadmergedGained
TH1F *	goodbadmergedLost
TH1F *	goodcluster
TH1F *	goodhittype
TH1F *	goodhittype_clgteq4
TH1F *	goodhittype_cllt4
TH1F *	goodhittype_simvecbig
TH1F *	goodhittype_simvecsmall
TH1F *	goodlayer
TH1F *	goodlayer_clgteq4
TH1F *	goodlayer_cllt4
TH1F *	goodlayer_simvecbig
TH1F *	goodlayer_simvecsmall
TH1F *	goodpix_clustersize
TH1F *	goodpix_simvecsize
TH1F *	goodpixgteq4_simvecsize
TH1F *	goodpixlt4_simvecsize
TH1F *	goodprj_clustersize
TH1F *	goodprj_simvecsize
TH1F *	goodprjgteq4_simvecsize
TH1F *	goodprjlt4_simvecsize
TH1F *	goodprocess
TH1F *	goodst1_clustersize
TH1F *	goodst1_simvecsize
TH1F *	goodst1gteq4_simvecsize
TH1F *	goodst1lt4_simvecsize
TH1F *	goodst2_clustersize
TH1F *	goodst2_simvecsize
TH1F *	goodst2gteq4_simvecsize
TH1F *	goodst2lt4_simvecsize
TH1F *	histoD0Old
TH1F *	histoD0Out
TH1F *	histoDzOld
TH1F *	histoDzOut
TH1F *	histoLambdaOld
TH1F *	histoLambdaOut
TH1F *	histoPhiOld
TH1F *	histoPhiOut
TH1F *	histoPtOld
TH1F *	histoPtOut
TH1F *	histoQoverpOld
TH1F *	histoQoverpOut
TH1F *	hitsPerTrackAssocGained
TH1F *	hitsPerTrackAssocLost
TH1F *	hitsPerTrackGained
TH1F *	hitsPerTrackLost
TH1F *	mergedcluster
TH1F *	mergedhittype
TH1F *	mergedlayer
TH1F *	mergedPull
TH1F *	nOfTrackIds
TH1F *	okcutsOld
TH1F *	okcutsOut
std::string	out
TH2F *	posxy
TH2F *	poszr
TH1F *	probXbad
TH1F *	probXdelta
TH1F *	probXgood
TH1F *	probXnoshare
TH1F *	probXshared
TH1F *	probYbad
TH1F *	probYdelta
TH1F *	probYgood
TH1F *	probYnoshare
TH1F *	probYshared
TH1F *	process
TH1F *	sizeOld
TH1F *	sizeOldT
TH1F *	sizeOut
TH1F *	sizeOutT
edm::EDGetTokenT< reco::TrackToTrackingParticleAssociator >	theAssociatorOldToken
edm::EDGetTokenT< reco::TrackToTrackingParticleAssociator >	theAssociatorOutToken
edm::ESHandle< TrackerGeometry >	theG
edm::InputTag	tpTags_
TrackerHitAssociator::Config	trackerHitAssociatorConfig_
TH1F *	tracks
edm::InputTag	trackTagsOld_
edm::InputTag	trackTagsOut_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
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)

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 55 of file TestOutliers.cc.

Constructor & Destructor Documentation

TestOutliers::TestOutliers ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 120 of file TestOutliers.cc.

References LogTrace.

  :
  trackTagsOut_(iConfig.getUntrackedParameter<edm::InputTag>("tracksOut")),
  trackTagsOld_(iConfig.getUntrackedParameter<edm::InputTag>("tracksOld")),
  tpTags_(iConfig.getUntrackedParameter<edm::InputTag>("tp")),
  trackerHitAssociatorConfig_(consumesCollector()),
  theAssociatorOldToken(consumes<reco::TrackToTrackingParticleAssociator>(iConfig.getUntrackedParameter<edm::InputTag>("TrackAssociatorByHitsOld"))),
  theAssociatorOutToken(consumes<reco::TrackToTrackingParticleAssociator>(iConfig.getUntrackedParameter<edm::InputTag>("TrackAssociatorByHitsOut"))),
  out(iConfig.getParameter<std::string>("out"))
{
  LogTrace("TestOutliers") <<"constructor";
//   ParameterSet cuts = iConfig.getParameter<ParameterSet>("RecoTracksCuts");
//   selectRecoTracks = RecoTrackSelector(cuts.getParameter<double>("ptMin"),
//                     cuts.getParameter<double>("minRapidity"),
//                     cuts.getParameter<double>("maxRapidity"),
//                     cuts.getParameter<double>("tip"),
//                     cuts.getParameter<double>("lip"),
//                     cuts.getParameter<int>("minHit"),
//                     cuts.getParameter<double>("maxChi2"));
  
  LogTrace("TestOutliers") <<"end constructor";
}

TestOutliers::~TestOutliers ( )

override

Definition at line 144 of file TestOutliers.cc.

{
 
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void TestOutliers::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 159 of file TestOutliers.cc.

References funct::abs(), TrackerHitAssociator::associateHit(), TrackerHitAssociator::associateHitId(), reco::TrackToTrackingParticleAssociator::associateRecoToSim(), badcluster, badhittype, badlayer, badprocess, ecalDrivenElectronSeedsParameters_cff::beamSpot, TrajectoryStateClosestToPoint::charge(), reco::TrackBase::charge(), funct::cos(), KineDebug3::count(), countOldT, countOutA, countOutT, reco::TrackBase::d0(), reco::TrackBase::d0Error(), delta, deltahits, deltahitsAssocGained, deltahitsAssocLost, deltahitsNO, deltahitsOK, reco::TrackBase::dsz(), reco::TrackBase::dxy(), reco::TrackBase::dz(), reco::TrackBase::dzError(), edm::AssociationMap< Tag >::end(), energyLoss, energyLossMax, energyLossRatio, reco::TrackBase::eta(), spr::find(), edm::AssociationMap< Tag >::find(), plotBeamSpotDB::first, gainedhits, gainedhits2, edm::EventSetup::get(), edm::RefToBase< T >::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), edm::Event::getByToken(), goodbadhits, goodbadmerged, goodbadmergedGained, goodbadmergedLost, goodcluster, goodhittype, goodhittype_clgteq4, goodhittype_cllt4, goodhittype_simvecbig, goodhittype_simvecsmall, goodlayer, goodlayer_clgteq4, goodlayer_cllt4, goodlayer_simvecbig, goodlayer_simvecsmall, goodpix_clustersize, goodpix_simvecsize, goodpixgteq4_simvecsize, goodpixlt4_simvecsize, goodprj_clustersize, goodprj_simvecsize, goodprjgteq4_simvecsize, goodprjlt4_simvecsize, goodprocess, goodst1_clustersize, goodst1_simvecsize, goodst1gteq4_simvecsize, goodst1lt4_simvecsize, goodst2_clustersize, goodst2_simvecsize, goodst2gteq4_simvecsize, goodst2lt4_simvecsize, histoD0Old, histoD0Out, histoDzOld, histoDzOut, histoLambdaOld, histoLambdaOut, histoPhiOld, histoPhiOut, histoPtOld, histoPtOut, histoQoverpOld, histoQoverpOut, TrackQuality_cff::hitAssociator, hitsPerTrackAssocGained, hitsPerTrackAssocLost, hitsPerTrackGained, hitsPerTrackLost, edm::EventBase::id(), TrackerGeometry::idToDet(), findQualityFiles::jj, gen::k, reco::TrackBase::lambda(), TrackerTopology::layer(), LogTrace, funct::m, M_PI, PV3DBase< T, PVType, FrameType >::mag(), makeMuonMisalignmentScenario::matrix, mergedcluster, mergedhittype, mergedlayer, mergedPull, CoreSimTrack::momentum(), FreeTrajectoryState::momentum(), nOfTrackIds, reco::TrackBase::normalizedChi2(), nullptr, reco::TrackBase::numberOfValidHits(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), reco::TrackBase::phi(), reco::TrackBase::phiError(), reco::BeamSpot::position(), FreeTrajectoryState::position(), posxy, poszr, SiPixelRecHit::probabilityX(), SiPixelRecHit::probabilityY(), probXbad, probXdelta, probXgood, probXnoshare, probXshared, probYbad, probYdelta, probYgood, probYnoshare, probYshared, process, PSimHit::processType(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), reco::TrackBase::pt(), reco::TrackBase::ptError(), reco::TrackBase::qoverp(), reco::TrackBase::qoverpError(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), reco::Track::seedRef(), funct::sin(), edm::AssociationMap< Tag >::size(), sizeOld, sizeOldT, sizeOut, sizeOutT, TrackingRecHit::some, jetUpdater_cfi::sort, mathSSE::sqrt(), theAssociatorOldToken, theAssociatorOutToken, theG, TrajectoryStateClosestToPoint::theState(), PV3DBase< T, PVType, FrameType >::theta(), reco::TrackBase::thetaError(), tmp, GeomDet::toGlobal(), tpTags_, trackerHitAssociatorConfig_, tracks, trackTagsOld_, trackTagsOut_, ErrorFrameTransformer::transform(), findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                         {
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopo;
  iSetup.get<TrackerTopologyRcd>().get(tTopo);



  using namespace edm;
  using namespace std;
  using reco::TrackCollection;

  LogTrace("TestOutliers") <<"analyze event #" << iEvent.id();
  
  edm::Handle<edm::View<reco::Track> > tracksOut;
  iEvent.getByLabel(trackTagsOut_,tracksOut);
  edm::Handle<edm::View<reco::Track> > tracksOld;
  iEvent.getByLabel(trackTagsOld_,tracksOld);
  Handle<TrackingParticleCollection> tps;
  iEvent.getByLabel(tpTags_,tps);
  edm::Handle<reco::BeamSpot> beamSpot;
  iEvent.getByLabel("offlineBeamSpot",beamSpot); 

  TrackerHitAssociator hitAssociator(iEvent, trackerHitAssociatorConfig_);

  edm::Handle<reco::TrackToTrackingParticleAssociator> hAssociatorOld;
  iEvent.getByToken(theAssociatorOldToken, hAssociatorOld);
  const reco::TrackToTrackingParticleAssociator *theAssociatorOld = hAssociatorOld.product();

  edm::Handle<reco::TrackToTrackingParticleAssociator> hAssociatorOut;
  iEvent.getByToken(theAssociatorOutToken, hAssociatorOut);
  const reco::TrackToTrackingParticleAssociator *theAssociatorOut = hAssociatorOut.product();

  reco::RecoToSimCollection recSimCollOut=theAssociatorOut->associateRecoToSim(tracksOut, tps);
  reco::RecoToSimCollection recSimCollOld=theAssociatorOld->associateRecoToSim(tracksOld, tps);
  sizeOut->Fill(recSimCollOut.size());
  sizeOld->Fill(recSimCollOld.size());
  sizeOutT->Fill(tracksOut->size());
  sizeOldT->Fill(tracksOld->size());

  LogTrace("TestOutliers") << "tOld size=" << tracksOld->size() << " tOut size=" << tracksOut->size() 
               << " aOld size=" << recSimCollOld.size() << " aOut size=" << recSimCollOut.size();

#if 0
  LogTrace("TestOutliers") << "recSimCollOld.size()=" << recSimCollOld.size() ;
  for(reco::TrackCollection::size_type j=0; j<tracksOld->size(); ++j){
    reco::TrackRef trackOld(tracksOld, j);
    //if ( !selectRecoTracks( *trackOld,beamSpot.product() ) ) continue;
    LogTrace("TestOutliers") << "trackOld->pt()=" << trackOld->pt() << " trackOld->numberOfValidHits()=" << trackOld->numberOfValidHits();
    std::vector<pair<TrackingParticleRef, double> > tpOld;
    if(recSimCollOld.find(trackOld) != recSimCollOld.end()){
      tpOld = recSimCollOld[trackOld];
      if (tpOld.size()!=0) LogTrace("TestOutliers") << " associated" ;
      else LogTrace("TestOutliers") << " NOT associated" ;
    } else LogTrace("TestOutliers") << " NOT associated" ;
  }
  LogTrace("TestOutliers") << "recSimCollOut.size()=" << recSimCollOut.size() ;
  for(reco::TrackCollection::size_type j=0; j<tracksOut->size(); ++j){
    reco::TrackRef trackOut(tracksOut, j);
    //if ( !selectRecoTracks( *trackOut,beamSpot.product() ) ) continue;
    LogTrace("TestOutliers") << "trackOut->pt()=" << trackOut->pt() << " trackOut->numberOfValidHits()=" << trackOut->numberOfValidHits();
    std::vector<pair<TrackingParticleRef, double> > tpOut;
    if(recSimCollOut.find(trackOut) != recSimCollOut.end()){
      tpOut = recSimCollOut[trackOut];
      if (tpOut.size()!=0) LogTrace("TestOutliers") << " associated" ;
      else LogTrace("TestOutliers") << " NOT associated" ;
    } else LogTrace("TestOutliers") << " NOT associated" ;
  }
#endif

  LogTrace("TestOutliers") <<"tracksOut->size()="<<tracksOut->size();
  LogTrace("TestOutliers") <<"tracksOld->size()="<<tracksOld->size();

  std::vector<unsigned int> outused;
  for(unsigned int j=0; j<tracksOld->size(); ++j) {
    countOldT->Fill(1);
    edm::RefToBase<reco::Track> trackOld(tracksOld, j);
    LogTrace("TestOutliers") << "now track old with id=" << j << " seed ref=" << trackOld->seedRef().get() << " pt=" << trackOld->pt()
                  << " eta=" <<  trackOld->eta() << " chi2=" <<  trackOld->normalizedChi2()
                  << " tip= " << fabs(trackOld->dxy(beamSpot->position()))
                  << " lip= " << fabs(trackOld->dsz(beamSpot->position()));

    //     if (i==tracksOut->size()) {
    //       if(recSimCollOld.find(trackOld) != recSimCollOld.end()){      
    //  vector<pair<TrackingParticleRef, double> > tpOld;
    //  tpOld = recSimCollOld[trackOld];
    //  if (tpOld.size()!=0) { 
    //    LogTrace("TestOutliers") <<"no match: old associated and out lost! old has #hits=" << trackOld->numberOfValidHits() 
    //                 << " and fraction=" << tpOld.begin()->second;
    //    if (tpOld.begin()->second>0.5) hitsPerTrackLost->Fill(trackOld->numberOfValidHits());
    //  }
    //       }
    //       continue;
    //     }

    std::vector<unsigned int> outtest;//all the out tracks with the same seed ref
    for(unsigned int k=0; k<tracksOut->size(); ++k) {
      edm::RefToBase<reco::Track> tmpOut = edm::RefToBase<reco::Track>(tracksOut, k);
      if ( tmpOut->seedRef() == trackOld->seedRef() ) {
    outtest.push_back(k);
      }      
    }

    edm::RefToBase<reco::Track> trackOut;
    if (outtest.size()==1) {// if only one that's it
      trackOut = edm::RefToBase<reco::Track>(tracksOut, outtest[0]);
      LogTrace("TestOutliers") << "now track out with id=" << outtest[0] << " seed ref=" << trackOut->seedRef().get() << " pt=" << trackOut->pt();  
      outused.push_back(outtest[0]);
    } else if (outtest.size()>1) {//if > 1 take the one that shares all the hits with the old track
      for(unsigned int p=0; p<outtest.size(); ++p) {
    edm::RefToBase<reco::Track> tmpOut = edm::RefToBase<reco::Track>(tracksOut, outtest[p]);    
    bool allhits = true;
    for (trackingRecHit_iterator itOut = tmpOut->recHitsBegin(); itOut!=tmpOut->recHitsEnd(); itOut++) {
      if ((*itOut)->isValid()) { 
        bool thishit = false;
        for (trackingRecHit_iterator itOld = trackOld->recHitsBegin();  itOld!=trackOld->recHitsEnd(); itOld++) {
          if ((*itOld)->isValid()) {
        const TrackingRecHit* kt = &(**itOld);
        if ( (*itOut)->sharesInput(kt,TrackingRecHit::some) ) { 
          thishit = true;
          break;
        }
          }
        }
        if (!thishit) allhits = false;
      }
    }
    if (allhits) {
      trackOut = edm::RefToBase<reco::Track>(tracksOut, outtest[p]);
      LogTrace("TestOutliers") << "now track out with id=" << outtest[p] << " seed ref=" << trackOut->seedRef().get() << " pt=" << trackOut->pt();  
      outused.push_back(outtest[p]);      
    }
      }
    } 

    if (outtest.empty() || trackOut.get()==nullptr ) {//no out track found for the old track
      if(recSimCollOld.find(trackOld) != recSimCollOld.end()){      
    vector<pair<TrackingParticleRef, double> > tpOld;
    tpOld = recSimCollOld[trackOld];
    if (!tpOld.empty()) { 
      LogTrace("TestOutliers") <<"no match: old associated and out lost! old has #hits=" << trackOld->numberOfValidHits() 
                   << " and fraction=" << tpOld.begin()->second;
      if (tpOld.begin()->second>0.5) hitsPerTrackLost->Fill(trackOld->numberOfValidHits());
    }
      }      
      LogTrace("TestOutliers") <<"...skip to next old track";
      continue;
    }

    //look if old and out are associated    
    LogTrace("TestOutliers") <<"trackOut->seedRef()=" << trackOut->seedRef().get() << " trackOld->seedRef()=" << trackOld->seedRef().get();
    bool oldAssoc = recSimCollOld.find(trackOld) != recSimCollOld.end();
    bool outAssoc = recSimCollOut.find(trackOut) != recSimCollOut.end();
    LogTrace("TestOutliers") <<"outAssoc=" << outAssoc <<" oldAssoc=" << oldAssoc;

    //     if ( trackOut->seedRef()!= trackOld->seedRef() ||  
    //   (trackOut->seedRef() == trackOld->seedRef() && trackOut->numberOfValidHits()>trackOld->numberOfValidHits()) ) {
    //       LogTrace("TestOutliers") <<"out and old tracks does not match...";
    //       LogTrace("TestOutliers") <<"old has #hits=" << trackOld->numberOfValidHits();
    //       std::vector<pair<TrackingParticleRef, double> > tpOld;
    //       if(recSimCollOld.find(trackOld) != recSimCollOld.end()) {
    //  tpOld = recSimCollOld[trackOld];
    //  if (tpOld.size()!=0) { 
    //    LogTrace("TestOutliers") <<"old was associated with fraction=" << tpOld.begin()->second;
    //    if (tpOld.begin()->second>0.5) hitsPerTrackLost->Fill(trackOld->numberOfValidHits());
    //  }
    //       }
    //       LogTrace("TestOutliers") <<"...skip to next old track";
    //       continue;
    //     }
    //     //++i;
    countOutT->Fill(1);

    //if ( !selectRecoTracks( *trackOld,beamSpot.product() ) ) continue;//no more cuts

    tracks->Fill(0);//FIXME

    TrackingParticleRef tpr;
    TrackingParticleRef tprOut;
    TrackingParticleRef tprOld;
    double fracOut;
    std::vector<unsigned int> tpids;
    std::vector<std::pair<TrackingParticleRef, double> > tpOut;
    std::vector<pair<TrackingParticleRef, double> > tpOld;
    //contare outliers delle tracce che erano associate e non lo sono piu!!!!

    if(outAssoc) {//save the ids od the tp associate to the out track
      tpOut = recSimCollOut[trackOut];
      if (!tpOut.empty()) {
    countOutA->Fill(1);
    tprOut = tpOut.begin()->first;
    fracOut = tpOut.begin()->second;
        for (TrackingParticle::g4t_iterator g4T=tprOut->g4Track_begin(); g4T!=tprOut->g4Track_end(); ++g4T) {
          tpids.push_back(g4T->trackId());
        }
      }
    }

    if(oldAssoc){//save the ids od the tp associate to the old track
      tpOld = recSimCollOld[trackOld];
      if (!tpOld.empty()) { 
    tprOld = tpOld.begin()->first;
    //  LogTrace("TestOutliers") <<"old associated and out not! old has #hits=" << trackOld->numberOfValidHits() 
    //               << " and fraction=" << tpOld.begin()->second;
    //  if (tpOld.begin()->second>0.5) hitsPerTrackLost->Fill(trackOld->numberOfValidHits());//deve essere un plot diverso tipo LostAssoc
    if (tpOut.empty()) {
      for (TrackingParticle::g4t_iterator g4T=tprOld->g4Track_begin(); g4T!=tprOld->g4Track_end(); ++g4T) {
        tpids.push_back(g4T->trackId());
      }
    }
      }
    }

    if (tprOut.get()!=nullptr || tprOld.get()!=nullptr) { //at least one of the tracks has to be associated

      tpr = tprOut.get()!=nullptr ? tprOut : tprOld;

      const SimTrack * assocTrack = &(*tpr->g4Track_begin());
      
      //if ( trackOut->numberOfValidHits() < trackOld->numberOfValidHits() ) {
      if ( trackOut->numberOfValidHits() != trackOld->numberOfValidHits() || 
       !(*trackOut->recHitsBegin())->sharesInput((*trackOld->recHitsBegin()),TrackingRecHit::some) ||
       !(*(trackOut->recHitsEnd()-1))->sharesInput((*(trackOld->recHitsEnd()-1)),TrackingRecHit::some)  ) 
    { //there are outliers if the number of valid hits is != or if the first and last hit does not match
    LogTrace("TestOutliers") << "outliers for track with #hits=" << trackOut->numberOfValidHits();
    tracks->Fill(1);
    LogTrace("TestOutliers") << "Out->pt=" << trackOut->pt() << " Old->pt=" << trackOld->pt() 
                 << " tp->pt=" << sqrt(tpr->momentum().perp2()) 
      //<< " trackOut->ptError()=" << trackOut->ptError() << " trackOld->ptError()=" << trackOld->ptError() 
                 << " Old->validHits=" << trackOld->numberOfValidHits() << " Out->validHits=" << trackOut->numberOfValidHits()
      /*<< " fracOld=" << fracOld*/ << " fracOut=" << fracOut
                 << " deltaHits=" << trackOld->numberOfValidHits()-trackOut->numberOfValidHits();

    //compute all the track parameters    
    double PtPullOut = (trackOut->pt()-sqrt(tpr->momentum().perp2()))/trackOut->ptError(); 
    double PtPullOld = (trackOld->pt()-sqrt(tpr->momentum().perp2()))/trackOld->ptError();
    histoPtOut->Fill( PtPullOut );
    histoPtOld->Fill( PtPullOld );

    //LogTrace("TestOutliers") << "MagneticField";        
    edm::ESHandle<MagneticField> theMF;
    iSetup.get<IdealMagneticFieldRecord>().get(theMF);
    FreeTrajectoryState 
      ftsAtProduction(GlobalPoint(tpr->vertex().x(),tpr->vertex().y(),tpr->vertex().z()),
              GlobalVector(assocTrack->momentum().x(),assocTrack->momentum().y(),assocTrack->momentum().z()),
              TrackCharge(trackOld->charge()),
              theMF.product());
    TSCPBuilderNoMaterial tscpBuilder;
    TrajectoryStateClosestToPoint tsAtClosestApproach 
      = tscpBuilder(ftsAtProduction,GlobalPoint(0,0,0));//as in TrackProducerAlgorithm
    GlobalPoint v = tsAtClosestApproach.theState().position();
    GlobalVector p = tsAtClosestApproach.theState().momentum();
          
    //LogTrace("TestOutliers") << "qoverpSim";        
    double qoverpSim = tsAtClosestApproach.charge()/p.mag();
    double lambdaSim = M_PI/2-p.theta();
    double phiSim    = p.phi();
    double dxySim    = (-v.x()*sin(p.phi())+v.y()*cos(p.phi()));
    double dszSim    = v.z()*p.perp()/p.mag() - (v.x()*p.x()+v.y()*p.y())/p.perp() * p.z()/p.mag();
    double d0Sim     = -dxySim;
    double dzSim     = dszSim*p.mag()/p.perp();
          
    //LogTrace("TestOutliers") << "qoverpPullOut";        
    double qoverpPullOut=(trackOut->qoverp()-qoverpSim)/trackOut->qoverpError();
    double qoverpPullOld=(trackOld->qoverp()-qoverpSim)/trackOld->qoverpError();
    double lambdaPullOut=(trackOut->lambda()-lambdaSim)/trackOut->thetaError();
    double lambdaPullOld=(trackOld->lambda()-lambdaSim)/trackOld->thetaError();
    double phi0PullOut=(trackOut->phi()-phiSim)/trackOut->phiError();
    double phi0PullOld=(trackOld->phi()-phiSim)/trackOld->phiError();
    double d0PullOut=(trackOut->d0()-d0Sim)/trackOut->d0Error();
    double d0PullOld=(trackOld->d0()-d0Sim)/trackOld->d0Error();
    double dzPullOut=(trackOut->dz()-dzSim)/trackOut->dzError();
    double dzPullOld=(trackOld->dz()-dzSim)/trackOld->dzError();
          
    //LogTrace("TestOutliers") << "histoQoverpOut";       
    histoQoverpOut->Fill(qoverpPullOut);
    histoQoverpOld->Fill(qoverpPullOld);
    histoPhiOut->Fill(phi0PullOut);  
    histoPhiOld->Fill(phi0PullOld);  
    histoD0Out->Fill(d0PullOut);  
    histoD0Old->Fill(d0PullOld);  
    histoDzOut->Fill(dzPullOut);  
    histoDzOld->Fill(dzPullOld);  
    histoLambdaOut->Fill(lambdaPullOut);
    histoLambdaOld->Fill(lambdaPullOld);

    //delta number of valid hits
    LogTrace("TestOutliers") << "deltahits=" << trackOld->numberOfValidHits()-trackOut->numberOfValidHits();          
    deltahits->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());

    if(tprOut.get()!=nullptr && tprOld.get()==nullptr) { //out associated and old not: gained track
      if (!tpOld.empty() && tpOld.begin()->second<=0.5) {
        deltahitsAssocGained->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
        hitsPerTrackAssocGained->Fill(trackOut->numberOfValidHits());
        LogTrace("TestOutliers") << "a) gained (assoc) track out #hits==" << trackOut->numberOfValidHits() << " old #hits=" << trackOld->numberOfValidHits();    
      } else {
        deltahitsAssocGained->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
        hitsPerTrackAssocGained->Fill(trackOut->numberOfValidHits());
        LogTrace("TestOutliers") << "b) gained (assoc) track out #hits==" << trackOut->numberOfValidHits() << " old #hits=" << trackOld->numberOfValidHits();    
      }
    } else if(tprOut.get()==nullptr && tprOld.get()!=nullptr) { //old associated and out not: lost track
      LogTrace("TestOutliers") <<"old associated and out not! old has #hits=" << trackOld->numberOfValidHits() 
                   << " and fraction=" << tpOld.begin()->second;
      if (tpOld.begin()->second>0.5) {      
        hitsPerTrackAssocLost->Fill(trackOld->numberOfValidHits());
        deltahitsAssocLost->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
      }
    }
    
    if ( fabs(PtPullOut) < fabs(PtPullOld) ) 
      deltahitsOK->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
    else 
      deltahitsNO->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
    
    //LogTrace("TestOutliers") << "RecoTrackSelector";        
    //if (selectRecoTracks(*trackOut,beamSpot.product())) okcutsOut->Fill(1); else okcutsOut->Fill(0);
    //if (selectRecoTracks(*trackOld,beamSpot.product())) okcutsOld->Fill(1); else okcutsOld->Fill(0);

    LogTrace("TestOutliers") << "track old";
    for (trackingRecHit_iterator itOld = trackOld->recHitsBegin(); itOld!=trackOld->recHitsEnd() ; itOld++){
      LogTrace("TestOutliers") << (*itOld)->isValid() << " " << (*itOld)->geographicalId().rawId();
    }
    LogTrace("TestOutliers") << "track out";
    for (trackingRecHit_iterator itOut = trackOut->recHitsBegin(); itOut!=trackOut->recHitsEnd() ; itOut++){
      LogTrace("TestOutliers") << (*itOut)->isValid() << " " << (*itOut)->geographicalId().rawId();
    }
    //LogTrace("TestOutliers") << "itOut";        


    vector<pair<int, trackingRecHit_iterator> > gainedlostoutliers;
    //look for gained hits
    for (trackingRecHit_iterator itOut = trackOut->recHitsBegin(); itOut!=trackOut->recHitsEnd(); itOut++){
      bool gained = true;
      if ((*itOut)->isValid()) {
        for (trackingRecHit_iterator itOld = trackOld->recHitsBegin(); itOld!=trackOld->recHitsEnd() ; itOld++){
          if ( (*itOld)->geographicalId().rawId()==(*itOut)->geographicalId().rawId() ) gained = false;
        }
        if (gained) {
          gainedlostoutliers.push_back(pair<int, trackingRecHit_iterator>(1,itOut));
          LogTrace("TestOutliers") << "broken trajectory during old fit... gained hit " << (*itOut)->geographicalId().rawId();
          gainedhits->Fill(1);
        }
      }
    }

    //look for outliers and lost hits
    for (trackingRecHit_iterator itOld = trackOld->recHitsBegin(); itOld!=trackOld->recHitsEnd() ; itOld++){

      bool outlier = false;
      bool lost = true;

      for (trackingRecHit_iterator itOut = trackOut->recHitsBegin(); itOut!=trackOut->recHitsEnd(); itOut++){
        if ( (*itOld)->geographicalId().rawId()==(*itOut)->geographicalId().rawId() ) {
          lost=false;
          if ( (*itOld)->isValid() && !(*itOut)->isValid() && (*itOld)->geographicalId().rawId()==(*itOut)->geographicalId().rawId() ) {
        LogTrace("TestOutliers") << (*itOld)->isValid() << " " << (*itOut)->isValid() << " " 
                     << (*itOld)->geographicalId().rawId() << " " << (*itOut)->geographicalId().rawId();
        outlier=true;
          }
        }
      }
      if (lost) gainedlostoutliers.push_back(pair<int, trackingRecHit_iterator>(2,itOld));
      if (lost) LogTrace("TestOutliers") << "lost";
      else if (outlier) gainedlostoutliers.push_back(pair<int, trackingRecHit_iterator>(3,itOld));
    }

    for (std::vector<pair<int, trackingRecHit_iterator> >::iterator it = gainedlostoutliers.begin(); it!=gainedlostoutliers.end();++it) {
      LogTrace("TestOutliers") << "type of processed hit:" <<it->first;
      trackingRecHit_iterator itHit = it->second;
      bool gained = false;
      bool lost = false;
      bool outlier = false;
      if (it->first==1) gained = true;
      else if (it->first==2) lost = true;
      else if (it->first==3) outlier = true;

      if (outlier||lost||gained) {
      //if (1) {

        if (lost && (*itHit)->isValid()==false) {
          goodbadmergedLost->Fill(0);
          LogTrace("TestOutliers") << "lost invalid";         
          continue;
        }
        else if (gained && (*itHit)->isValid()==false) {
          goodbadmergedGained->Fill(0);
          LogTrace("TestOutliers") << "gained invalid";
          continue;
        }
          
        //LogTrace("TestOutliers") << "vector<SimHitIdpr>";       
        //look if the hit comes from a correct sim track
        std::vector<SimHitIdpr> simTrackIds = hitAssociator.associateHitId(**itHit);
        bool goodhit = false;
        for(size_t j=0; j<simTrackIds.size(); j++){
          for (size_t jj=0; jj<tpids.size(); jj++){
        if (simTrackIds[j].first == tpids[jj]) goodhit = true;
        break;
          }
        }

        //find what kind of hit is it
        int clustersize = 0;
        int hittypeval = 0;
        int layerval = 0 ;
        if (dynamic_cast<const SiPixelRecHit*>(&**itHit)){  
          LogTrace("TestOutliers") << "SiPixelRecHit";        
          clustersize =  ((const SiPixelRecHit*)(&**itHit))->cluster()->size() ;
          hittypeval  = 1;
        }
        else if (dynamic_cast<const SiStripRecHit2D*>(&**itHit)){
          LogTrace("TestOutliers") << "SiStripRecHit2D";          
          clustersize =  ((const SiStripRecHit2D*)(&**itHit))->cluster()->amplitudes().size() ;
          hittypeval  = 2;
        }
        else if (dynamic_cast<const SiStripMatchedRecHit2D*>(&**itHit)){
          LogTrace("TestOutliers") << "SiStripMatchedRecHit2D";       
          int clsize1 = ((const SiStripMatchedRecHit2D*)(&**itHit))->monoCluster().amplitudes().size();
          int clsize2 =  ((const SiStripMatchedRecHit2D*)(&**itHit))->stereoCluster().amplitudes().size();
          if (clsize1>clsize2) clustersize = clsize1;
          else clustersize = clsize2;
          hittypeval  = 3;
        }
        else if (dynamic_cast<const ProjectedSiStripRecHit2D*>(&**itHit)){
          LogTrace("TestOutliers") << "ProjectedSiStripRecHit2D";         
          clustersize =  ((const ProjectedSiStripRecHit2D*)(&**itHit))->originalHit().cluster()->amplitudes().size();
          hittypeval  = 4;
        }
          
        //find the layer of the hit
        int subdetId = (*itHit)->geographicalId().subdetId();
        DetId id = (*itHit)->geographicalId();
        int layerId  = tTopo->layer(id);
        layerval = subdetId*10+layerId;
      
        //LogTrace("TestOutliers") << "gpos";         
        GlobalPoint gpos = theG->idToDet((*itHit)->geographicalId())->surface().toGlobal((*itHit)->localPosition());


        //get the vector of sim hit associated and choose the one with the largest energy loss
        //double delta = 99999;
        //LocalPoint rhitLPv = (*itHit)->localPosition();
        //vector<PSimHit> assSimHits = hitAssociator.associateHit(**itHit);
        //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();
        // }
        //}
        //LogTrace("TestOutliers") << "energyLoss_";       
        double energyLoss_ = 0.;
        unsigned int monoId = 0;
        std::vector<double> energyLossM;
        std::vector<double> energyLossS;
        std::vector<PSimHit> assSimHits = hitAssociator.associateHit(**itHit);
        if (assSimHits.empty()) continue;
        PSimHit shit;
        std::vector<unsigned int> trackIds;
        energyLossS.clear();
        energyLossM.clear();
        //LogTrace("TestOutliers") << "energyLossM";          
        for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){
          if (outlier) energyLoss->Fill(m->energyLoss());
          unsigned int tId = m->trackId();
          if (find(trackIds.begin(),trackIds.end(),tId)==trackIds.end()) trackIds.push_back(tId);
          LogTrace("TestOutliers") << "id=" << tId;
          if (m->energyLoss()>energyLoss_) {
        shit=*m;
        energyLoss_ = m->energyLoss();
          }
          if (hittypeval==3) {
        if (monoId==0) monoId = m->detUnitId();
        if (monoId == m->detUnitId()){
          energyLossM.push_back(m->energyLoss());
        }
        else {
          energyLossS.push_back(m->energyLoss());
        }
        //std::cout << "detUnitId="  << m->detUnitId() << " trackId=" << m->trackId() << " energyLoss=" << m->energyLoss() << std::endl;
          } else {
        energyLossM.push_back(m->energyLoss());
          }
        }
        unsigned int nIds = trackIds.size();

        if (outlier) {
          goodbadhits->Fill(goodhit);
          posxy->Fill(fabs(gpos.x()),fabs(gpos.y()));
          poszr->Fill(fabs(gpos.z()),sqrt(gpos.x()*gpos.x()+gpos.y()*gpos.y()));
          process->Fill(shit.processType());              
          energyLossMax->Fill(energyLoss_);
        }
          
        //look if the hit is shared and if is produced only by ionization processes
        bool shared = true;
        bool ioniOnly = true;
        unsigned int idc = 0;
        for (size_t jj=0; jj<tpids.size(); jj++){
          idc += std::count(trackIds.begin(),trackIds.end(),tpids[jj]);
        }
        if (idc==trackIds.size()) {
          shared = false;
        }
        for(std::vector<PSimHit>::const_iterator m=assSimHits.begin()+1; m<assSimHits.end(); m++){
          if ((m->processType()!=7&&m->processType()!=8&&m->processType()!=9)&&abs(m->particleType())!=11){
        ioniOnly = false;
        break;
          }
        }
        if (ioniOnly&&!shared){
          LogTrace("TestOutliers") << "delta";
        }
        
        if (goodhit) { 
          if (outlier) {
        goodprocess->Fill(shit.processType());
        if (clustersize>=4) {
          goodhittype_clgteq4->Fill(hittypeval);
          goodlayer_clgteq4->Fill(layerval);
        } else {
          goodhittype_cllt4->Fill(hittypeval);
          goodlayer_cllt4->Fill(layerval);
        }
        //LogTrace("TestOutliers") << "hittypeval && clustersize";        
        if (hittypeval==1 && clustersize>=4) goodpixgteq4_simvecsize->Fill(assSimHits.size());
        if (hittypeval==1 && clustersize<4 ) goodpixlt4_simvecsize->Fill(assSimHits.size());
        if (hittypeval==2 && clustersize>=4) goodst1gteq4_simvecsize->Fill(assSimHits.size());
        if (hittypeval==2 && clustersize<4 ) goodst1lt4_simvecsize->Fill(assSimHits.size());
        if (hittypeval==3 && clustersize>=4) goodst2gteq4_simvecsize->Fill(assSimHits.size());
        if (hittypeval==3 && clustersize<4 ) goodst2lt4_simvecsize->Fill(assSimHits.size());
        if (hittypeval==4 && clustersize>=4) goodprjgteq4_simvecsize->Fill(assSimHits.size());
        if (hittypeval==4 && clustersize<4 ) goodprjlt4_simvecsize->Fill(assSimHits.size());
          
        //LogTrace("TestOutliers") << "hittypeval";       
        if (hittypeval==1) goodpix_clustersize->Fill(clustersize);
        if (hittypeval==2) goodst1_clustersize->Fill(clustersize);
        if (hittypeval==3) goodst2_clustersize->Fill(clustersize);
        if (hittypeval==4) goodprj_clustersize->Fill(clustersize);
        if (hittypeval==1) goodpix_simvecsize->Fill(assSimHits.size());
        if (hittypeval==2) goodst1_simvecsize->Fill(assSimHits.size());
        if (hittypeval==3) goodst2_simvecsize->Fill(assSimHits.size());
        if (hittypeval==4) goodprj_simvecsize->Fill(assSimHits.size());
                  
        //LogTrace("TestOutliers") << "nOfTrackIds";          
        nOfTrackIds->Fill(nIds);
        if (hittypeval!=3) { 
          if (energyLossM.size()>1) {
            sort(energyLossM.begin(),energyLossM.end(),greater<double>());
            energyLossRatio->Fill(energyLossM[1]/energyLossM[0]);
          }
        } else {
          //LogTrace("TestOutliers") << "hittypeval==3";          
          if (energyLossM.size()>1&&energyLossS.size()<=1) {
            sort(energyLossM.begin(),energyLossM.end(),greater<double>());
            energyLossRatio->Fill(energyLossM[1]/energyLossM[0]);
          }
          else if (energyLossS.size()>1&&energyLossM.size()<=1) {
            sort(energyLossS.begin(),energyLossS.end(),greater<double>());
            energyLossRatio->Fill(energyLossS[1]/energyLossS[0]);
          }
          else if (energyLossS.size()>1&&energyLossM.size()>1) {
            sort(energyLossM.begin(),energyLossM.end(),greater<double>());
            sort(energyLossS.begin(),energyLossS.end(),greater<double>());
            energyLossM[1]/energyLossM[0] > energyLossS[1]/energyLossS[0] 
              ? energyLossRatio->Fill(energyLossM[1]/energyLossM[0]) 
              : energyLossRatio->Fill(energyLossS[1]/energyLossS[0]);
          }
        }
          
        LogTrace("TestOutliers") << "before merged";          
        const SiStripMatchedRecHit2D* tmp = dynamic_cast<const SiStripMatchedRecHit2D*>(&**itHit);
        LogTrace("TestOutliers") << "tmp=" << tmp; 
        LogTrace("TestOutliers") << "assSimHits.size()=" << assSimHits.size(); 
        if ( (assSimHits.size()>1 && tmp==nullptr) || 
             (assSimHits.size()>2 && tmp!=nullptr) ) {
          //std::cout << "MERGED HIT" << std::endl;
          //LogTrace("TestOutliers") << "merged";         
          mergedlayer->Fill(layerval);
          mergedcluster->Fill(clustersize);
          mergedhittype->Fill(hittypeval);

          for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){
            unsigned int tId = m->trackId();
            LogTrace("TestOutliers") << "component with id=" << tId <<  " eLoss=" << m->energyLoss() 
                         << " proc=" <<  m->processType() << " part=" <<  m->particleType();
            if (find(tpids.begin(),tpids.end(),tId)==tpids.end()) continue;
            if (m->processType()==2) {
              //GlobalPoint gpos = theG->idToDet((*itHit)->geographicalId())->surface().toGlobal((*itHit)->localPosition());
              //GlobalPoint gpr = rhit->globalPosition();
              AlgebraicSymMatrix33 ger = 
            ErrorFrameTransformer().transform((*itHit)->localPositionError(),theG->idToDet((*itHit)->geographicalId())->surface() ).matrix();
              //AlgebraicSymMatrix ger = rhit->globalPositionError().matrix();
              GlobalPoint gps = theG->idToDet((*itHit)->geographicalId())->surface().toGlobal(m->localPosition());
              //LogTrace("TestOutliers") << gpr << " " << gps << " " << ger;
              ROOT::Math::SVector<double,3> delta;
              delta[0]=gpos.x()-gps.x();
              delta[1]=gpos.y()-gps.y();
              delta[2]=gpos.z()-gps.z();
              LogTrace("TestOutliers") << delta << " " << ger ;
              double mpull = sqrt(delta[0]*delta[0]/ger[0][0]+delta[1]*delta[1]/ger[1][1]+delta[2]*delta[2]/ger[2][2]);
              LogTrace("TestOutliers") << "hit pull=" << mpull;//ger.similarity(delta);
              mergedPull->Fill(mpull);
            }
          }
          LogTrace("TestOutliers") << "end merged";
        } else {//not merged=>good
          //LogTrace("TestOutliers") << "goodlayer";          
          goodlayer->Fill(layerval);          
          goodcluster->Fill(clustersize);
          goodhittype->Fill(hittypeval);
        }   
          }//if (outlier)

          const SiPixelRecHit* pix = dynamic_cast<const SiPixelRecHit*>(&**itHit);
          if ((hittypeval!=3 && assSimHits.size()<2)||(hittypeval==3 && assSimHits.size()<3)){
        if (outlier) {
          goodhittype_simvecsmall->Fill(hittypeval);
          goodlayer_simvecsmall->Fill(layerval);
          goodbadmerged->Fill(1);
          if (pix) {
            probXgood->Fill(pix->probabilityX());
            probYgood->Fill(pix->probabilityY());
          }
          LogTrace("TestOutliers") << "out good";         
        }
        else if (lost){
          goodbadmergedLost->Fill(1);
          LogTrace("TestOutliers") << "lost good";        
        }
        else if (gained){
          goodbadmergedGained->Fill(1);
          LogTrace("TestOutliers") << "gained good";          
        }
        LogTrace("TestOutliers") << "good"; 
          } else {
        if (outlier) {
          goodhittype_simvecbig->Fill(hittypeval);
          goodlayer_simvecbig->Fill(layerval);
          if (ioniOnly&&!shared) { 
            goodbadmerged->Fill(3);
            if (pix) {
              probXdelta->Fill(pix->probabilityX());
              probYdelta->Fill(pix->probabilityY());
            }
          } else if(!ioniOnly&&!shared) {
            goodbadmerged->Fill(4);
            if (pix) {
              probXnoshare->Fill(pix->probabilityX());
              probYnoshare->Fill(pix->probabilityY());
            }
          } else {
            goodbadmerged->Fill(5);
            if (pix) {
              probXshared->Fill(pix->probabilityX());
              probYshared->Fill(pix->probabilityY());
            }
          }
          LogTrace("TestOutliers") << "out merged, ioniOnly=" << ioniOnly << " shared=" << shared;
        }
        else if (lost) {
          if (ioniOnly&&!shared) goodbadmergedLost->Fill(3);
          else if(!ioniOnly&&!shared) goodbadmergedLost->Fill(4);
          else goodbadmergedLost->Fill(5);
          LogTrace("TestOutliers") << "lost merged, ioniOnly=" << ioniOnly << " shared=" << shared;
        }
        else if (gained) {
          if (ioniOnly&&!shared) goodbadmergedGained->Fill(3);
          else if(!ioniOnly&&!shared) goodbadmergedGained->Fill(4);
          else goodbadmergedGained->Fill(5);
          LogTrace("TestOutliers") << "gained merged, ioniOnly=" << ioniOnly << " shared=" << shared;
        }
        LogTrace("TestOutliers") << "merged, ioniOnly=" << ioniOnly << " shared=" << shared;
          }
        } //if (goodhit)
        else {//badhit
          //LogTrace("TestOutliers") << "badhit";   
          if (outlier) {
        badcluster->Fill(clustersize);
        badhittype->Fill(hittypeval);
        badlayer->Fill(layerval);
        badprocess->Fill(shit.processType());
        goodbadmerged->Fill(2);
        const SiPixelRecHit* pix = dynamic_cast<const SiPixelRecHit*>(&**itHit);
        if (pix) {
          probXbad->Fill(pix->probabilityX());
          probYbad->Fill(pix->probabilityY());
        }
        LogTrace("TestOutliers") << "out bad";        
          }
          else if (lost) {
        goodbadmergedLost->Fill(2);
        LogTrace("TestOutliers") << "lost bad";       
          }
          else if (gained) {
        goodbadmergedGained->Fill(2);
        LogTrace("TestOutliers") << "gained bad";         
          }
          LogTrace("TestOutliers") << "bad"; 
        }
      }
    }
      } 
      //else if ( trackOut->numberOfValidHits() > trackOld->numberOfValidHits() ) {
      else if ( false ) {
    LogTrace("TestOutliers") << "outliers for track with #hits=" << trackOut->numberOfValidHits();
    tracks->Fill(1);
    LogTrace("TestOutliers") << "Out->pt=" << trackOut->pt() << " Old->pt=" << trackOld->pt() 
                 << " tp->pt=" << sqrt(tpr->momentum().perp2()) 
      //<< " trackOut->ptError()=" << trackOut->ptError() << " trackOld->ptError()=" << trackOld->ptError() 
                 << " Old->validHits=" << trackOld->numberOfValidHits() << " Out->validHits=" << trackOut->numberOfValidHits()
      /*<< " fracOld=" << fracOld*/ << " fracOut=" << fracOut
                 << " deltaHits=" << trackOld->numberOfValidHits()-trackOut->numberOfValidHits();
    LogTrace("TestOutliers") << "track with gained hits";
    gainedhits2->Fill(trackOut->numberOfValidHits()-trackOld->numberOfValidHits()); 
      } else {
    LogTrace("TestOutliers") << "no outliers for track with #hits=" << trackOut->numberOfValidHits();
      }
    }
    LogTrace("TestOutliers") << "end track old #" << j;
  }  

  for (unsigned int k=0;k<tracksOut->size();++k) {
    if ( find(outused.begin(),outused.end(),k)==outused.end() ) {
      edm::RefToBase<reco::Track> trackOut(tracksOut, k);      
      bool outAssoc = recSimCollOut.find(trackOut) != recSimCollOut.end();
      if (outAssoc) {
    hitsPerTrackGained->Fill(trackOut->numberOfValidHits());
    LogTrace("TestOutliers") << "gained track out id=" << k << " #hits==" << trackOut->numberOfValidHits();    

      }
    }    
  }
}

void TestOutliers::beginRun ( edm::Run const &  run, const edm::EventSetup &  es  )

overrideprivate

Definition at line 897 of file TestOutliers.cc.

References badcluster, badhittype, badlayer, badprocess, countOldT, countOutA, countOutT, deltahits, deltahitsAssocGained, deltahitsAssocLost, deltahitsNO, deltahitsOK, energyLoss, energyLossMax, energyLossRatio, file, gainedhits, gainedhits2, edm::EventSetup::get(), goodbadhits, goodbadmerged, goodbadmergedGained, goodbadmergedLost, goodcluster, goodhittype, goodhittype_clgteq4, goodhittype_cllt4, goodhittype_simvecbig, goodhittype_simvecsmall, goodlayer, goodlayer_clgteq4, goodlayer_cllt4, goodlayer_simvecbig, goodlayer_simvecsmall, goodpix_clustersize, goodpix_simvecsize, goodpixgteq4_simvecsize, goodpixlt4_simvecsize, goodprj_clustersize, goodprj_simvecsize, goodprjgteq4_simvecsize, goodprjlt4_simvecsize, goodprocess, goodst1_clustersize, goodst1_simvecsize, goodst1gteq4_simvecsize, goodst1lt4_simvecsize, goodst2_clustersize, goodst2_simvecsize, goodst2gteq4_simvecsize, goodst2lt4_simvecsize, histoD0Old, histoD0Out, histoDzOld, histoDzOut, histoLambdaOld, histoLambdaOut, histoPhiOld, histoPhiOut, histoPtOld, histoPtOut, histoQoverpOld, histoQoverpOut, hitsPerTrackAssocGained, hitsPerTrackAssocLost, hitsPerTrackGained, hitsPerTrackLost, mergedcluster, mergedhittype, mergedlayer, mergedPull, nOfTrackIds, okcutsOld, okcutsOut, out, posxy, poszr, probXbad, probXdelta, probXgood, probXnoshare, probXshared, probYbad, probYdelta, probYgood, probYnoshare, probYshared, process, sizeOld, sizeOldT, sizeOut, sizeOutT, theG, and tracks.

{
  es.get<TrackerDigiGeometryRecord>().get(theG);
  const bool oldAddDir = TH1::AddDirectoryStatus();
  TH1::AddDirectory(true);
  file = new TFile(out.c_str(),"recreate");
  histoPtOut = new TH1F("histoPtOut","histoPtOut",100,-10,10);
  histoPtOld = new TH1F("histoPtOld","histoPtOld",100,-10,10);
  histoQoverpOut = new TH1F("histoQoverpOut","histoQoverpOut",250,-25,25);
  histoQoverpOld = new TH1F("histoQoverpOld","histoQoverpOld",250,-25,25);
  histoPhiOut = new TH1F("histoPhiOut","histoPhiOut",250,-25,25);
  histoPhiOld = new TH1F("histoPhiOld","histoPhiOld",250,-25,25);
  histoD0Out = new TH1F("histoD0Out","histoD0Out",250,-25,25);
  histoD0Old = new TH1F("histoD0Old","histoD0Old",250,-25,25);
  histoDzOut = new TH1F("histoDzOut","histoDzOut",250,-25,25);
  histoDzOld = new TH1F("histoDzOld","histoDzOld",250,-25,25);
  histoLambdaOut = new TH1F("histoLambdaOut","histoLambdaOut",250,-25,25);
  histoLambdaOld = new TH1F("histoLambdaOld","histoLambdaOld",250,-25,25);
  deltahits = new TH1F("deltahits","deltahits",80,-40,40);
  deltahitsOK = new TH1F("deltahitsOK","deltahitsOK",20,0,20);
  deltahitsNO = new TH1F("deltahitsNO","deltahitsNO",20,0,20);
  okcutsOut = new TH1F("okcutsOut","okcutsOut",2,-0.5,1.5);
  okcutsOld = new TH1F("okcutsOld","okcutsOld",2,-0.5,1.5);
  tracks = new TH1F("tracks_","tracks_",2,-0.5,1.5);
  goodbadhits = new TH1F("goodbadhits","goodbadhits",2,-0.5,1.5);
  process = new TH1F("process","process",20,-0.5,19.5);
  goodcluster = new TH1F("goodcluster","goodcluster",40,-0.5,39.5);
  goodprocess = new TH1F("goodprocess","goodprocess",20,-0.5,19.5);
  badcluster = new TH1F("badcluster","badcluster",40,-0.5,39.5);
  badprocess = new TH1F("badprocess","badprocess",20,-0.5,19.5);
  goodhittype = new TH1F("goodhittype","goodhittype",5,-0.5,4.5);
  goodlayer = new TH1F("goodlayer","goodlayer",70,-0.5,69.5);
  goodhittype_clgteq4 = new TH1F("goodhittype_clgteq4","goodhittype_clgteq4",5,-0.5,4.5);
  goodlayer_clgteq4 = new TH1F("goodlayer_clgteq4","goodlayer_clgteq4",70,-0.5,69.5);
  goodhittype_cllt4 = new TH1F("goodhittype_cllt4","goodhittype_cllt4",5,-0.5,4.5);
  goodlayer_cllt4 = new TH1F("goodlayer_cllt4","goodlayer_cllt4",70,-0.5,69.5);
  badhittype = new TH1F("badhittype","badhittype",5,-0.5,4.5);
  badlayer = new TH1F("badlayer","badlayer",70,-0.5,69.5);
  posxy = new TH2F("posxy","posxy",1200,0,120,1200,0,120);
  poszr = new TH2F("poszr","poszr",3000,0,300,1200,0,120);
  goodpixgteq4_simvecsize = new TH1F("goodpixgteq4_simvecsize","goodpixgteq4_simvecsize",40,-0.5,39.5);
  goodpixlt4_simvecsize  = new TH1F("goodpixlt4_simvecsize","goodpixlt4_simvecsize",40,-0.5,39.5);
  goodst1gteq4_simvecsize = new TH1F("goodst1gteq4_simvecsize","goodst1gteq4_simvecsize",40,-0.5,39.5);
  goodst1lt4_simvecsize  = new TH1F("goodst1lt4_simvecsize","goodst1lt4_simvecsize",40,-0.5,39.5);
  goodst2gteq4_simvecsize = new TH1F("goodst2gteq4_simvecsize","goodst2gteq4_simvecsize",40,-0.5,39.5);
  goodst2lt4_simvecsize  = new TH1F("goodst2lt4_simvecsize","goodst2lt4_simvecsize",40,-0.5,39.5);
  goodprjgteq4_simvecsize = new TH1F("goodprjgteq4_simvecsize","goodprjgteq4_simvecsize",40,-0.5,39.5);
  goodprjlt4_simvecsize  = new TH1F("goodprjlt4_simvecsize","goodprjlt4_simvecsize",40,-0.5,39.5);
  goodpix_clustersize = new TH1F("goodpix_clustersize","goodpix_clustersize",40,-0.5,39.5);
  goodst1_clustersize = new TH1F("goodst1_clustersize","goodst1_clustersize",40,-0.5,39.5);
  goodst2_clustersize = new TH1F("goodst2_clustersize","goodst2_clustersize",40,-0.5,39.5);
  goodprj_clustersize = new TH1F("goodprj_clustersize","goodprj_clustersize",40,-0.5,39.5);
  goodpix_simvecsize = new TH1F("goodpix_simvecsize","goodpix_simvecsize",40,-0.5,39.5);
  goodst1_simvecsize = new TH1F("goodst1_simvecsize","goodst1_simvecsize",40,-0.5,39.5);
  goodst2_simvecsize = new TH1F("goodst2_simvecsize","goodst2_simvecsize",40,-0.5,39.5);
  goodprj_simvecsize = new TH1F("goodprj_simvecsize","goodprj_simvecsize",40,-0.5,39.5);
  goodhittype_simvecsmall = new TH1F("goodhittype_simvecsmall","goodhittype_simvecsmall",5,-0.5,4.5);
  goodlayer_simvecsmall = new TH1F("goodlayer_simvecsmall","goodlayer_simvecsmall",70,-0.5,69.5);
  goodhittype_simvecbig = new TH1F("goodhittype_simvecbig","goodhittype_simvecbig",5,-0.5,4.5);
  goodlayer_simvecbig = new TH1F("goodlayer_simvecbig","goodlayer_simvecbig",70,-0.5,69.5);
  goodbadmerged = new TH1F("goodbadmerged","goodbadmerged",5,0.5,5.5);
  goodbadmergedLost = new TH1F("goodbadmergedLost","goodbadmergedLost",5,0.5,5.5);
  goodbadmergedGained = new TH1F("goodbadmergedGained","goodbadmergedGained",5,0.5,5.5);
  energyLoss = new TH1F("energyLoss","energyLoss",1000,0,0.1);
  energyLossMax = new TH1F("energyLossMax","energyLossMax",1000,0,0.1);
  energyLossRatio = new TH1F("energyLossRatio","energyLossRatio",100,0,1);
  nOfTrackIds = new TH1F("nOfTrackIds","nOfTrackIds",10,0,10);
  mergedPull = new TH1F("mergedPull","mergedPull",100,0,10);
  mergedlayer = new TH1F("mergedlayer","mergedlayer",70,-0.5,69.5);
  mergedhittype = new TH1F("mergedhittype","mergedhittype",5,-0.5,4.5);
  mergedcluster = new TH1F("mergedcluster","mergedcluster",40,-0.5,39.5);
  deltahitsAssocGained = new TH1F("deltahitsAssocGained","deltahitsAssocGained",80, -40, 40);
  deltahitsAssocLost = new TH1F("deltahitsAssocLost","deltahitsAssocLost",80, -40, 40);
  hitsPerTrackLost = new TH1F("hitsPerTrackLost","hitsPerTrackLost",40, -0.5, 39.5);
  hitsPerTrackAssocLost = new TH1F("hitsPerTrackAssocLost","hitsPerTrackAssocLost",40, -0.5, 39.5);
  hitsPerTrackGained = new TH1F("hitsPerTrackGained","hitsPerTrackGained",40, -0.5, 39.5);
  hitsPerTrackAssocGained = new TH1F("hitsPerTrackAssocGained","hitsPerTrackAssocGained",40, -0.5, 39.5);
  sizeOut = new TH1F("sizeOut","sizeOut",900,-0.5,899.5);
  sizeOld = new TH1F("sizeOld","sizeOld",900,-0.5,899.5);
  sizeOutT = new TH1F("sizeOutT","sizeOutT",900,-0.5,899.5);
  sizeOldT = new TH1F("sizeOldT","sizeOldT",900,-0.5,899.5);
  countOutA = new TH1F("countOutA","countOutA",2,0,2);
  countOutT = new TH1F("countOutT","countOutT",2,0,2);
  countOldT = new TH1F("countOldT","countOldT",2,0,2);
  gainedhits = new TH1F("gainedhits","gainedhits",2,0,2);
  gainedhits2 = new TH1F("gainedhits2","gainedhits2",30,-0.5,29.5);
  probXgood = new TH1F("probXgood","probXgood",110,0,1.1);
  probXbad = new TH1F("probXbad","probXbad",110,0,1.1);
  probXdelta = new TH1F("probXdelta","probXdelta",110,0,1.1);
  probXshared = new TH1F("probXshared","probXshared",110,0,1.1);
  probXnoshare = new TH1F("probXnoshare","probXnoshare",110,0,1.1);
  probYgood = new TH1F("probYgood","probYgood",110,0,1.1);
  probYbad = new TH1F("probYbad","probYbad",110,0,1.1);
  probYdelta = new TH1F("probYdelta","probYdelta",110,0,1.1);
  probYshared = new TH1F("probYshared","probYshared",110,0,1.1);
  probYnoshare = new TH1F("probYnoshare","probYnoshare",110,0,1.1);
  TH1::AddDirectory(oldAddDir);
}

void TestOutliers::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 997 of file TestOutliers.cc.

References DEFINE_FWK_MODULE, file, and LogTrace.

                     {
  LogTrace("TestOutliers") << "TestOutliers::endJob";
  file->Write();
  LogTrace("TestOutliers") << "outfile written";
  file->Close();
  LogTrace("TestOutliers") << "oufile closed";
  LogTrace("TestOutliers") << "exiting TestOutliers::endJob";
}

Member Data Documentation

TH1F * TestOutliers::badcluster

private

Definition at line 84 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::badhittype

private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::badlayer

private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::badprocess

private

Definition at line 84 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::countOldT

private

Definition at line 102 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::countOutA

private

Definition at line 102 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::countOutT

private

Definition at line 102 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::deltahits

private

Definition at line 82 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::deltahitsAssocGained

private

Definition at line 82 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::deltahitsAssocLost

private

Definition at line 82 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::deltahitsNO

private

Definition at line 82 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::deltahitsOK

private

Definition at line 82 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::energyLoss

private

Definition at line 100 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::energyLossMax

private

Definition at line 100 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::energyLossRatio

private

Definition at line 100 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TFile* TestOutliers::file

private

Definition at line 75 of file TestOutliers.cc.

Referenced by beginRun(), and endJob().

TH1F* TestOutliers::gainedhits

private

Definition at line 103 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::gainedhits2

private

Definition at line 103 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodbadhits

private

Definition at line 84 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodbadmerged

private

Definition at line 99 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodbadmergedGained

private

Definition at line 99 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodbadmergedLost

private

Definition at line 99 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodcluster

private

Definition at line 84 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodhittype

private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodhittype_clgteq4

private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodhittype_cllt4

private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodhittype_simvecbig

private

Definition at line 99 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodhittype_simvecsmall

private

Definition at line 99 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer

private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer_clgteq4

private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer_cllt4

private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer_simvecbig

private

Definition at line 99 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer_simvecsmall

private

Definition at line 99 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodpix_clustersize

private

Definition at line 91 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodpix_simvecsize

private

Definition at line 95 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodpixgteq4_simvecsize

private

Definition at line 87 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodpixlt4_simvecsize

private

Definition at line 87 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodprj_clustersize

private

Definition at line 94 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodprj_simvecsize

private

Definition at line 98 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodprjgteq4_simvecsize

private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodprjlt4_simvecsize

private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodprocess

private

Definition at line 84 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst1_clustersize

private

Definition at line 92 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst1_simvecsize

private

Definition at line 96 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst1gteq4_simvecsize

private

Definition at line 88 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodst1lt4_simvecsize

private

Definition at line 88 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst2_clustersize

private

Definition at line 93 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst2_simvecsize

private

Definition at line 97 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst2gteq4_simvecsize

private

Definition at line 89 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodst2lt4_simvecsize

private

Definition at line 89 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoD0Old

private

Definition at line 79 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoD0Out

private

Definition at line 79 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoDzOld

private

Definition at line 80 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoDzOut

private

Definition at line 80 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoLambdaOld

private

Definition at line 81 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoLambdaOut

private

Definition at line 81 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoPhiOld

private

Definition at line 78 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoPhiOut

private

Definition at line 78 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoPtOld

private

Definition at line 76 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoPtOut

private

Definition at line 76 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoQoverpOld

private

Definition at line 77 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoQoverpOut

private

Definition at line 77 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::hitsPerTrackAssocGained

private

Definition at line 82 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::hitsPerTrackAssocLost

private

Definition at line 82 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::hitsPerTrackGained

private

Definition at line 82 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::hitsPerTrackLost

private

Definition at line 82 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::mergedcluster

private

Definition at line 100 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::mergedhittype

private

Definition at line 100 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::mergedlayer

private

Definition at line 100 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::mergedPull

private

Definition at line 100 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::nOfTrackIds

private

Definition at line 100 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::okcutsOld

private

Definition at line 83 of file TestOutliers.cc.

Referenced by beginRun().

TH1F* TestOutliers::okcutsOut

private

Definition at line 83 of file TestOutliers.cc.

Referenced by beginRun().

std::string TestOutliers::out

private

Definition at line 74 of file TestOutliers.cc.

Referenced by beginRun().

TH2F* TestOutliers::posxy

private

Definition at line 86 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH2F * TestOutliers::poszr

private

Definition at line 86 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probXbad

private

Definition at line 104 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probXdelta

private

Definition at line 104 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::probXgood

private

Definition at line 104 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probXnoshare

private

Definition at line 104 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probXshared

private

Definition at line 104 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probYbad

private

Definition at line 105 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probYdelta

private

Definition at line 105 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::probYgood

private

Definition at line 105 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probYnoshare

private

Definition at line 105 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probYshared

private

Definition at line 105 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::process

private

Definition at line 84 of file TestOutliers.cc.

Referenced by analyze(), beginRun(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::dumpPython(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::open(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::setProcess(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::setProperty().

TH1F * TestOutliers::sizeOld

private

Definition at line 101 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::sizeOldT

private

Definition at line 101 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::sizeOut

private

Definition at line 101 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::sizeOutT

private

Definition at line 101 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

edm::EDGetTokenT<reco::TrackToTrackingParticleAssociator> TestOutliers::theAssociatorOldToken

private

Definition at line 71 of file TestOutliers.cc.

Referenced by analyze().

edm::EDGetTokenT<reco::TrackToTrackingParticleAssociator> TestOutliers::theAssociatorOutToken

private

Definition at line 72 of file TestOutliers.cc.

Referenced by analyze().

edm::ESHandle<TrackerGeometry> TestOutliers::theG

private

Definition at line 73 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

edm::InputTag TestOutliers::tpTags_

private

Definition at line 69 of file TestOutliers.cc.

Referenced by analyze().

TrackerHitAssociator::Config TestOutliers::trackerHitAssociatorConfig_

private

Definition at line 70 of file TestOutliers.cc.

Referenced by analyze().

TH1F* TestOutliers::tracks

private

Definition at line 84 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

edm::InputTag TestOutliers::trackTagsOld_

private

Definition at line 68 of file TestOutliers.cc.

Referenced by analyze().

edm::InputTag TestOutliers::trackTagsOut_

private

Definition at line 67 of file TestOutliers.cc.

Referenced by analyze().