CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
TestOutliers Class Reference

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

Inheritance diagram for TestOutliers:
edm::EDAnalyzer

Public Member Functions

 TestOutliers (const edm::ParameterSet &)
 
 ~TestOutliers ()
 
- Public Member Functions inherited from edm::EDAnalyzer
 EDAnalyzer ()
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 

Private Member Functions

virtual void analyze (const edm::Event &, const edm::EventSetup &)
 
virtual void beginRun (edm::Run &run, const edm::EventSetup &)
 
virtual void endJob ()
 

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
 
TrackerHitAssociator * hitAssociator
 
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
 
edm::ParameterSet psetold
 
edm::ParameterSet psetout
 
TH1F * sizeOld
 
TH1F * sizeOldT
 
TH1F * sizeOut
 
TH1F * sizeOutT
 
TrackAssociatorBase * theAssociatorOld
 
TrackAssociatorBase * theAssociatorOut
 
edm::ESHandle< TrackerGeometry > theG
 
edm::InputTag tpTags_
 
TH1F * tracks
 
edm::InputTag trackTagsOld_
 
edm::InputTag trackTagsOut_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
typedef WorkerT< EDAnalyzer > WorkerType
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const * currentContext () const
 

Detailed Description

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

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

Definition at line 62 of file TestOutliers.cc.

Constructor & Destructor Documentation

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

Definition at line 129 of file TestOutliers.cc.

References edm::ParameterSet::getParameter(), LogTrace, psetold, and psetout.

130  :
131  trackTagsOut_(iConfig.getUntrackedParameter<edm::InputTag>("tracksOut")),
132  trackTagsOld_(iConfig.getUntrackedParameter<edm::InputTag>("tracksOld")),
133  tpTags_(iConfig.getUntrackedParameter<edm::InputTag>("tp")),
134  out(iConfig.getParameter<std::string>("out"))
135 {
136  LogTrace("TestOutliers") <<"constructor";
137 // ParameterSet cuts = iConfig.getParameter<ParameterSet>("RecoTracksCuts");
138 // selectRecoTracks = RecoTrackSelector(cuts.getParameter<double>("ptMin"),
139 // cuts.getParameter<double>("minRapidity"),
140 // cuts.getParameter<double>("maxRapidity"),
141 // cuts.getParameter<double>("tip"),
142 // cuts.getParameter<double>("lip"),
143 // cuts.getParameter<int>("minHit"),
144 // cuts.getParameter<double>("maxChi2"));
145 
146  psetold = iConfig.getParameter<ParameterSet>("TrackAssociatorByHitsPSetOld");
147  psetout = iConfig.getParameter<ParameterSet>("TrackAssociatorByHitsPSetOut");
148  LogTrace("TestOutliers") <<"end constructor";
149 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
TestOutliers::psetold
edm::ParameterSet psetold
Definition: TestOutliers.cc:114
TestOutliers::out
std::string out
Definition: TestOutliers.cc:82
TestOutliers::psetout
edm::ParameterSet psetout
Definition: TestOutliers.cc:114
TestOutliers::trackTagsOut_
edm::InputTag trackTagsOut_
Definition: TestOutliers.cc:74
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:517
TestOutliers::tpTags_
edm::InputTag tpTags_
Definition: TestOutliers.cc:76
TestOutliers::trackTagsOld_
edm::InputTag trackTagsOld_
Definition: TestOutliers.cc:75
TestOutliers::~TestOutliers ( )

Definition at line 152 of file TestOutliers.cc.

153 {
154 
155  // do anything here that needs to be done at desctruction time
156  // (e.g. close files, deallocate resources etc.)
157 
158 }

Member Function Documentation

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

Implements edm::EDAnalyzer.

Definition at line 167 of file TestOutliers.cc.

References abs, TrackerHitAssociator::associateHit(), TrackerHitAssociator::associateHitId(), TrackAssociatorBase::associateRecoToSim(), badcluster, badhittype, badlayer, badprocess, reco::TrackBase::charge(), TrajectoryStateClosestToPoint::charge(), funct::cos(), prof2calltree::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(), energyLoss, energyLossMax, energyLossRatio, reco::TrackBase::eta(), spr::find(), first, gainedhits, gainedhits2, edm::EventSetup::get(), edm::RefToBase< T >::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), 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, hitAssociator, hitsPerTrackAssocGained, hitsPerTrackAssocLost, hitsPerTrackGained, hitsPerTrackLost, edm::EventBase::id(), j, gen::k, reco::TrackBase::lambda(), LogTrace, m, M_PI, PV3DBase< T, PVType, FrameType >::mag(), mergedcluster, mergedhittype, mergedlayer, mergedPull, CoreSimTrack::momentum(), FreeTrajectoryState::momentum(), nOfTrackIds, reco::TrackBase::normalizedChi2(), reco::TrackBase::numberOfValidHits(), L1TEmulatorMonitor_cff::p, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), reco::TrackBase::phi(), reco::TrackBase::phiError(), FreeTrajectoryState::position(), posxy, poszr, SiPixelRecHit::probabilityX(), SiPixelRecHit::probabilityY(), probXbad, probXdelta, probXgood, probXnoshare, probXshared, probYbad, probYdelta, probYgood, probYnoshare, probYshared, process, PSimHit::processType(), edm::ESHandle< class >::product(), psetold, psetout, reco::TrackBase::pt(), reco::TrackBase::ptError(), reco::TrackBase::qoverp(), reco::TrackBase::qoverpError(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), reco::Track::seedRef(), funct::sin(), sizeOld, sizeOldT, sizeOut, sizeOutT, TrackingRecHit::some, python.multivaluedict::sort(), mathSSE::sqrt(), theAssociatorOld, theAssociatorOut, theG, TrajectoryStateClosestToPoint::theState(), PV3DBase< T, PVType, FrameType >::theta(), reco::TrackBase::thetaError(), tmp, tpTags_, TrackAssociatorByHits_cfi::TrackAssociatorByHits, tracks, trackTagsOld_, trackTagsOut_, ErrorFrameTransformer::transform(), v, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

167  {
168 
169  using namespace edm;
170  using namespace std;
171  using reco::TrackCollection;
172 
173  LogTrace("TestOutliers") <<"analyze event #" << iEvent.id();
174 
175  edm::Handle<edm::View<reco::Track> > tracksOut;
176  iEvent.getByLabel(trackTagsOut_,tracksOut);
177  edm::Handle<edm::View<reco::Track> > tracksOld;
178  iEvent.getByLabel(trackTagsOld_,tracksOld);
179  Handle<TrackingParticleCollection> tps;
180  iEvent.getByLabel(tpTags_,tps);
181  edm::Handle<reco::BeamSpot> beamSpot;
182  iEvent.getByLabel("offlineBeamSpot",beamSpot);
183 
184  hitAssociator = new TrackerHitAssociator(iEvent);
185 
186  theAssociatorOld = new TrackAssociatorByHits(psetold);
187  theAssociatorOut = new TrackAssociatorByHits(psetout);
188  reco::RecoToSimCollection recSimCollOut=theAssociatorOut->associateRecoToSim(tracksOut, tps, &iEvent);
189  reco::RecoToSimCollection recSimCollOld=theAssociatorOld->associateRecoToSim(tracksOld, tps, &iEvent);
190  sizeOut->Fill(recSimCollOut.size());
191  sizeOld->Fill(recSimCollOld.size());
192  sizeOutT->Fill(tracksOut->size());
193  sizeOldT->Fill(tracksOld->size());
194 
195  LogTrace("TestOutliers") << "tOld size=" << tracksOld->size() << " tOut size=" << tracksOut->size()
196  << " aOld size=" << recSimCollOld.size() << " aOut size=" << recSimCollOut.size();
197 
198 #if 0
199  LogTrace("TestOutliers") << "recSimCollOld.size()=" << recSimCollOld.size() ;
200  for(reco::TrackCollection::size_type j=0; j<tracksOld->size(); ++j){
201  reco::TrackRef trackOld(tracksOld, j);
202  //if ( !selectRecoTracks( *trackOld,beamSpot.product() ) ) continue;
203  LogTrace("TestOutliers") << "trackOld->pt()=" << trackOld->pt() << " trackOld->numberOfValidHits()=" << trackOld->numberOfValidHits();
204  std::vector<pair<TrackingParticleRef, double> > tpOld;
205  if(recSimCollOld.find(trackOld) != recSimCollOld.end()){
206  tpOld = recSimCollOld[trackOld];
207  if (tpOld.size()!=0) LogTrace("TestOutliers") << " associated" ;
208  else LogTrace("TestOutliers") << " NOT associated" ;
209  } else LogTrace("TestOutliers") << " NOT associated" ;
210  }
211  LogTrace("TestOutliers") << "recSimCollOut.size()=" << recSimCollOut.size() ;
212  for(reco::TrackCollection::size_type j=0; j<tracksOut->size(); ++j){
213  reco::TrackRef trackOut(tracksOut, j);
214  //if ( !selectRecoTracks( *trackOut,beamSpot.product() ) ) continue;
215  LogTrace("TestOutliers") << "trackOut->pt()=" << trackOut->pt() << " trackOut->numberOfValidHits()=" << trackOut->numberOfValidHits();
216  std::vector<pair<TrackingParticleRef, double> > tpOut;
217  if(recSimCollOut.find(trackOut) != recSimCollOut.end()){
218  tpOut = recSimCollOut[trackOut];
219  if (tpOut.size()!=0) LogTrace("TestOutliers") << " associated" ;
220  else LogTrace("TestOutliers") << " NOT associated" ;
221  } else LogTrace("TestOutliers") << " NOT associated" ;
222  }
223 #endif
224 
225  LogTrace("TestOutliers") <<"tracksOut->size()="<<tracksOut->size();
226  LogTrace("TestOutliers") <<"tracksOld->size()="<<tracksOld->size();
227 
228  std::vector<unsigned int> outused;
229  for(unsigned int j=0; j<tracksOld->size(); ++j) {
230  countOldT->Fill(1);
231  edm::RefToBase<reco::Track> trackOld(tracksOld, j);
232  LogTrace("TestOutliers") << "now track old with id=" << j << " seed ref=" << trackOld->seedRef().get() << " pt=" << trackOld->pt()
233  << " eta=" << trackOld->eta() << " chi2=" << trackOld->normalizedChi2()
234  << " tip= " << fabs(trackOld->dxy(beamSpot->position()))
235  << " lip= " << fabs(trackOld->dsz(beamSpot->position()));
236 
237  // if (i==tracksOut->size()) {
238  // if(recSimCollOld.find(trackOld) != recSimCollOld.end()){
239  // vector<pair<TrackingParticleRef, double> > tpOld;
240  // tpOld = recSimCollOld[trackOld];
241  // if (tpOld.size()!=0) {
242  // LogTrace("TestOutliers") <<"no match: old associated and out lost! old has #hits=" << trackOld->numberOfValidHits()
243  // << " and fraction=" << tpOld.begin()->second;
244  // if (tpOld.begin()->second>0.5) hitsPerTrackLost->Fill(trackOld->numberOfValidHits());
245  // }
246  // }
247  // continue;
248  // }
249 
250  std::vector<unsigned int> outtest;//all the out tracks with the same seed ref
251  for(unsigned int k=0; k<tracksOut->size(); ++k) {
252  edm::RefToBase<reco::Track> tmpOut = edm::RefToBase<reco::Track>(tracksOut, k);
253  if ( tmpOut->seedRef() == trackOld->seedRef() ) {
254  outtest.push_back(k);
255  }
256  }
257 
258  edm::RefToBase<reco::Track> trackOut;
259  if (outtest.size()==1) {// if only one that's it
260  trackOut = edm::RefToBase<reco::Track>(tracksOut, outtest[0]);
261  LogTrace("TestOutliers") << "now track out with id=" << outtest[0] << " seed ref=" << trackOut->seedRef().get() << " pt=" << trackOut->pt();
262  outused.push_back(outtest[0]);
263  } else if (outtest.size()>1) {//if > 1 take the one that shares all the hits with the old track
264  for(unsigned int p=0; p<outtest.size(); ++p) {
265  edm::RefToBase<reco::Track> tmpOut = edm::RefToBase<reco::Track>(tracksOut, outtest[p]);
266  bool allhits = true;
267  for (trackingRecHit_iterator itOut = tmpOut->recHitsBegin(); itOut!=tmpOut->recHitsEnd(); itOut++) {
268  if ((*itOut)->isValid()) {
269  bool thishit = false;
270  for (trackingRecHit_iterator itOld = trackOld->recHitsBegin(); itOld!=trackOld->recHitsEnd(); itOld++) {
271  if ((*itOld)->isValid()) {
272  const TrackingRecHit* kt = &(**itOld);
273  if ( (*itOut)->sharesInput(kt,TrackingRecHit::some) ) {
274  thishit = true;
275  break;
276  }
277  }
278  }
279  if (!thishit) allhits = false;
280  }
281  }
282  if (allhits) {
283  trackOut = edm::RefToBase<reco::Track>(tracksOut, outtest[p]);
284  LogTrace("TestOutliers") << "now track out with id=" << outtest[p] << " seed ref=" << trackOut->seedRef().get() << " pt=" << trackOut->pt();
285  outused.push_back(outtest[p]);
286  }
287  }
288  }
289 
290  if (outtest.size()==0 || trackOut.get()==0 ) {//no out track found for the old track
291  if(recSimCollOld.find(trackOld) != recSimCollOld.end()){
292  vector<pair<TrackingParticleRef, double> > tpOld;
293  tpOld = recSimCollOld[trackOld];
294  if (tpOld.size()!=0) {
295  LogTrace("TestOutliers") <<"no match: old associated and out lost! old has #hits=" << trackOld->numberOfValidHits()
296  << " and fraction=" << tpOld.begin()->second;
297  if (tpOld.begin()->second>0.5) hitsPerTrackLost->Fill(trackOld->numberOfValidHits());
298  }
299  }
300  LogTrace("TestOutliers") <<"...skip to next old track";
301  continue;
302  }
303 
304  //look if old and out are associated
305  LogTrace("TestOutliers") <<"trackOut->seedRef()=" << trackOut->seedRef().get() << " trackOld->seedRef()=" << trackOld->seedRef().get();
306  bool oldAssoc = recSimCollOld.find(trackOld) != recSimCollOld.end();
307  bool outAssoc = recSimCollOut.find(trackOut) != recSimCollOut.end();
308  LogTrace("TestOutliers") <<"outAssoc=" << outAssoc <<" oldAssoc=" << oldAssoc;
309 
310  // if ( trackOut->seedRef()!= trackOld->seedRef() ||
311  // (trackOut->seedRef() == trackOld->seedRef() && trackOut->numberOfValidHits()>trackOld->numberOfValidHits()) ) {
312  // LogTrace("TestOutliers") <<"out and old tracks does not match...";
313  // LogTrace("TestOutliers") <<"old has #hits=" << trackOld->numberOfValidHits();
314  // std::vector<pair<TrackingParticleRef, double> > tpOld;
315  // if(recSimCollOld.find(trackOld) != recSimCollOld.end()) {
316  // tpOld = recSimCollOld[trackOld];
317  // if (tpOld.size()!=0) {
318  // LogTrace("TestOutliers") <<"old was associated with fraction=" << tpOld.begin()->second;
319  // if (tpOld.begin()->second>0.5) hitsPerTrackLost->Fill(trackOld->numberOfValidHits());
320  // }
321  // }
322  // LogTrace("TestOutliers") <<"...skip to next old track";
323  // continue;
324  // }
325  // //++i;
326  countOutT->Fill(1);
327 
328  //if ( !selectRecoTracks( *trackOld,beamSpot.product() ) ) continue;//no more cuts
329 
330  tracks->Fill(0);//FIXME
331 
332  bool hasOut = false;
333 
334  TrackingParticleRef tpr;
335  TrackingParticleRef tprOut;
336  TrackingParticleRef tprOld;
337  double fracOut;
338  std::vector<unsigned int> tpids;
339  std::vector<std::pair<TrackingParticleRef, double> > tpOut;
340  std::vector<pair<TrackingParticleRef, double> > tpOld;
341  //contare outliers delle tracce che erano associate e non lo sono piu!!!!
342 
343  if(outAssoc) {//save the ids od the tp associate to the out track
344  tpOut = recSimCollOut[trackOut];
345  if (tpOut.size()!=0) {
346  countOutA->Fill(1);
347  tprOut = tpOut.begin()->first;
348  fracOut = tpOut.begin()->second;
349  for (TrackingParticle::g4t_iterator g4T=tprOut->g4Track_begin(); g4T!=tprOut->g4Track_end(); ++g4T) {
350  tpids.push_back(g4T->trackId());
351  }
352  }
353  }
354 
355  if(oldAssoc){//save the ids od the tp associate to the old track
356  tpOld = recSimCollOld[trackOld];
357  if (tpOld.size()!=0) {
358  tprOld = tpOld.begin()->first;
359  // LogTrace("TestOutliers") <<"old associated and out not! old has #hits=" << trackOld->numberOfValidHits()
360  // << " and fraction=" << tpOld.begin()->second;
361  // if (tpOld.begin()->second>0.5) hitsPerTrackLost->Fill(trackOld->numberOfValidHits());//deve essere un plot diverso tipo LostAssoc
362  if (tpOut.size()==0) {
363  for (TrackingParticle::g4t_iterator g4T=tprOld->g4Track_begin(); g4T!=tprOld->g4Track_end(); ++g4T) {
364  tpids.push_back(g4T->trackId());
365  }
366  }
367  }
368  }
369 
370  if (tprOut.get()!=0 || tprOld.get()!=0) { //at least one of the tracks has to be associated
371 
372  tpr = tprOut.get()!=0 ? tprOut : tprOld;
373 
374  const SimTrack * assocTrack = &(*tpr->g4Track_begin());
375 
376  //if ( trackOut->numberOfValidHits() < trackOld->numberOfValidHits() ) {
377  if ( trackOut->numberOfValidHits() != trackOld->numberOfValidHits() ||
378  !trackOut->recHitsBegin()->get()->sharesInput(trackOld->recHitsBegin()->get(),TrackingRecHit::some) ||
379  !(trackOut->recHitsEnd()-1)->get()->sharesInput((trackOld->recHitsEnd()-1)->get(),TrackingRecHit::some) )
380  { //there are outliers if the number of valid hits is != or if the first and last hit does not match
381  hasOut=true;
382  LogTrace("TestOutliers") << "outliers for track with #hits=" << trackOut->numberOfValidHits();
383  tracks->Fill(1);
384  LogTrace("TestOutliers") << "Out->pt=" << trackOut->pt() << " Old->pt=" << trackOld->pt()
385  << " tp->pt=" << sqrt(tpr->momentum().perp2())
386  //<< " trackOut->ptError()=" << trackOut->ptError() << " trackOld->ptError()=" << trackOld->ptError()
387  << " Old->validHits=" << trackOld->numberOfValidHits() << " Out->validHits=" << trackOut->numberOfValidHits()
388  /*<< " fracOld=" << fracOld*/ << " fracOut=" << fracOut
389  << " deltaHits=" << trackOld->numberOfValidHits()-trackOut->numberOfValidHits();
390 
391  //compute all the track parameters
392  double PtPullOut = (trackOut->pt()-sqrt(tpr->momentum().perp2()))/trackOut->ptError();
393  double PtPullOld = (trackOld->pt()-sqrt(tpr->momentum().perp2()))/trackOld->ptError();
394  histoPtOut->Fill( PtPullOut );
395  histoPtOld->Fill( PtPullOld );
396 
397  //LogTrace("TestOutliers") << "MagneticField";
398  edm::ESHandle<MagneticField> theMF;
399  iSetup.get<IdealMagneticFieldRecord>().get(theMF);
400  FreeTrajectoryState
401  ftsAtProduction(GlobalPoint(tpr->vertex().x(),tpr->vertex().y(),tpr->vertex().z()),
402  GlobalVector(assocTrack->momentum().x(),assocTrack->momentum().y(),assocTrack->momentum().z()),
403  TrackCharge(trackOld->charge()),
404  theMF.product());
405  TSCPBuilderNoMaterial tscpBuilder;
406  TrajectoryStateClosestToPoint tsAtClosestApproach
407  = tscpBuilder(ftsAtProduction,GlobalPoint(0,0,0));//as in TrackProducerAlgorithm
408  GlobalPoint v = tsAtClosestApproach.theState().position();
409  GlobalVector p = tsAtClosestApproach.theState().momentum();
410 
411  //LogTrace("TestOutliers") << "qoverpSim";
412  double qoverpSim = tsAtClosestApproach.charge()/p.mag();
413  double lambdaSim = M_PI/2-p.theta();
414  double phiSim = p.phi();
415  double dxySim = (-v.x()*sin(p.phi())+v.y()*cos(p.phi()));
416  double dszSim = v.z()*p.perp()/p.mag() - (v.x()*p.x()+v.y()*p.y())/p.perp() * p.z()/p.mag();
417  double d0Sim = -dxySim;
418  double dzSim = dszSim*p.mag()/p.perp();
419 
420  //LogTrace("TestOutliers") << "qoverpPullOut";
421  double qoverpPullOut=(trackOut->qoverp()-qoverpSim)/trackOut->qoverpError();
422  double qoverpPullOld=(trackOld->qoverp()-qoverpSim)/trackOld->qoverpError();
423  double lambdaPullOut=(trackOut->lambda()-lambdaSim)/trackOut->thetaError();
424  double lambdaPullOld=(trackOld->lambda()-lambdaSim)/trackOld->thetaError();
425  double phi0PullOut=(trackOut->phi()-phiSim)/trackOut->phiError();
426  double phi0PullOld=(trackOld->phi()-phiSim)/trackOld->phiError();
427  double d0PullOut=(trackOut->d0()-d0Sim)/trackOut->d0Error();
428  double d0PullOld=(trackOld->d0()-d0Sim)/trackOld->d0Error();
429  double dzPullOut=(trackOut->dz()-dzSim)/trackOut->dzError();
430  double dzPullOld=(trackOld->dz()-dzSim)/trackOld->dzError();
431 
432  //LogTrace("TestOutliers") << "histoQoverpOut";
433  histoQoverpOut->Fill(qoverpPullOut);
434  histoQoverpOld->Fill(qoverpPullOld);
435  histoPhiOut->Fill(phi0PullOut);
436  histoPhiOld->Fill(phi0PullOld);
437  histoD0Out->Fill(d0PullOut);
438  histoD0Old->Fill(d0PullOld);
439  histoDzOut->Fill(dzPullOut);
440  histoDzOld->Fill(dzPullOld);
441  histoLambdaOut->Fill(lambdaPullOut);
442  histoLambdaOld->Fill(lambdaPullOld);
443 
444  //delta number of valid hits
445  LogTrace("TestOutliers") << "deltahits=" << trackOld->numberOfValidHits()-trackOut->numberOfValidHits();
446  deltahits->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
447 
448  if(tprOut.get()!=0 && tprOld.get()==0) { //out associated and old not: gained track
449  if (tpOld.size()!=0 && tpOld.begin()->second<=0.5) {
450  deltahitsAssocGained->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
451  hitsPerTrackAssocGained->Fill(trackOut->numberOfValidHits());
452  LogTrace("TestOutliers") << "a) gained (assoc) track out #hits==" << trackOut->numberOfValidHits() << " old #hits=" << trackOld->numberOfValidHits();
453  } else {
454  deltahitsAssocGained->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
455  hitsPerTrackAssocGained->Fill(trackOut->numberOfValidHits());
456  LogTrace("TestOutliers") << "b) gained (assoc) track out #hits==" << trackOut->numberOfValidHits() << " old #hits=" << trackOld->numberOfValidHits();
457  }
458  } else if(tprOut.get()==0 && tprOld.get()!=0) { //old associated and out not: lost track
459  LogTrace("TestOutliers") <<"old associated and out not! old has #hits=" << trackOld->numberOfValidHits()
460  << " and fraction=" << tpOld.begin()->second;
461  if (tpOld.begin()->second>0.5) {
462  hitsPerTrackAssocLost->Fill(trackOld->numberOfValidHits());
463  deltahitsAssocLost->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
464  }
465  }
466 
467  if ( fabs(PtPullOut) < fabs(PtPullOld) )
468  deltahitsOK->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
469  else
470  deltahitsNO->Fill(trackOld->numberOfValidHits()-trackOut->numberOfValidHits());
471 
472  //LogTrace("TestOutliers") << "RecoTrackSelector";
473  //if (selectRecoTracks(*trackOut,beamSpot.product())) okcutsOut->Fill(1); else okcutsOut->Fill(0);
474  //if (selectRecoTracks(*trackOld,beamSpot.product())) okcutsOld->Fill(1); else okcutsOld->Fill(0);
475 
476  LogTrace("TestOutliers") << "track old";
477  for (trackingRecHit_iterator itOld = trackOld->recHitsBegin(); itOld!=trackOld->recHitsEnd() ; itOld++){
478  LogTrace("TestOutliers") << (*itOld)->isValid() << " " << (*itOld)->geographicalId().rawId();
479  }
480  LogTrace("TestOutliers") << "track out";
481  for (trackingRecHit_iterator itOut = trackOut->recHitsBegin(); itOut!=trackOut->recHitsEnd() ; itOut++){
482  LogTrace("TestOutliers") << (*itOut)->isValid() << " " << (*itOut)->geographicalId().rawId();
483  }
484  //LogTrace("TestOutliers") << "itOut";
485 
486 
487  vector<pair<int, trackingRecHit_iterator> > gainedlostoutliers;
488  //look for gained hits
489  for (trackingRecHit_iterator itOut = trackOut->recHitsBegin(); itOut!=trackOut->recHitsEnd(); itOut++){
490  bool gained = true;
491  if ((*itOut)->isValid()) {
492  for (trackingRecHit_iterator itOld = trackOld->recHitsBegin(); itOld!=trackOld->recHitsEnd() ; itOld++){
493  if ( (*itOld)->geographicalId().rawId()==(*itOut)->geographicalId().rawId() ) gained = false;
494  }
495  if (gained) {
496  gainedlostoutliers.push_back(pair<int, trackingRecHit_iterator>(1,itOut));
497  LogTrace("TestOutliers") << "broken trajectory during old fit... gained hit " << (*itOut)->geographicalId().rawId();
498  gainedhits->Fill(1);
499  }
500  }
501  }
502 
503  //look for outliers and lost hits
504  for (trackingRecHit_iterator itOld = trackOld->recHitsBegin(); itOld!=trackOld->recHitsEnd() ; itOld++){
505 
506  bool outlier = false;
507  bool lost = true;
508 
509  for (trackingRecHit_iterator itOut = trackOut->recHitsBegin(); itOut!=trackOut->recHitsEnd(); itOut++){
510  if ( (*itOld)->geographicalId().rawId()==(*itOut)->geographicalId().rawId() ) {
511  lost=false;
512  if ( (*itOld)->isValid() && !(*itOut)->isValid() && (*itOld)->geographicalId().rawId()==(*itOut)->geographicalId().rawId() ) {
513  LogTrace("TestOutliers") << (*itOld)->isValid() << " " << (*itOut)->isValid() << " "
514  << (*itOld)->geographicalId().rawId() << " " << (*itOut)->geographicalId().rawId();
515  outlier=true;
516  }
517  }
518  }
519  if (lost) gainedlostoutliers.push_back(pair<int, trackingRecHit_iterator>(2,itOld));
520  if (lost) LogTrace("TestOutliers") << "lost";
521  else if (outlier) gainedlostoutliers.push_back(pair<int, trackingRecHit_iterator>(3,itOld));
522  }
523 
524  for (std::vector<pair<int, trackingRecHit_iterator> >::iterator it = gainedlostoutliers.begin(); it!=gainedlostoutliers.end();++it) {
525  LogTrace("TestOutliers") << "type of processed hit:" <<it->first;
526  trackingRecHit_iterator itHit = it->second;
527  bool gained = false;
528  bool lost = false;
529  bool outlier = false;
530  if (it->first==1) gained = true;
531  else if (it->first==2) lost = true;
532  else if (it->first==3) outlier = true;
533 
534  if (outlier||lost||gained) {
535  //if (1) {
536 
537  if (lost && (*itHit)->isValid()==false) {
538  goodbadmergedLost->Fill(0);
539  LogTrace("TestOutliers") << "lost invalid";
540  continue;
541  }
542  else if (gained && (*itHit)->isValid()==false) {
543  goodbadmergedGained->Fill(0);
544  LogTrace("TestOutliers") << "gained invalid";
545  continue;
546  }
547 
548  //LogTrace("TestOutliers") << "vector<SimHitIdpr>";
549  //look if the hit comes from a correct sim track
550  std::vector<SimHitIdpr> simTrackIds = hitAssociator->associateHitId(**itHit);
551  bool goodhit = false;
552  for(size_t j=0; j<simTrackIds.size(); j++){
553  for (size_t jj=0; jj<tpids.size(); jj++){
554  if (simTrackIds[j].first == tpids[jj]) goodhit = true;
555  break;
556  }
557  }
558 
559  //find what kind of hit is it
560  int clustersize = 0;
561  int hittypeval = 0;
562  int layerval = 0 ;
563  if (dynamic_cast<const SiPixelRecHit*>(&**itHit)){
564  LogTrace("TestOutliers") << "SiPixelRecHit";
565  clustersize = ((const SiPixelRecHit*)(&**itHit))->cluster()->size() ;
566  hittypeval = 1;
567  }
568  else if (dynamic_cast<const SiStripRecHit2D*>(&**itHit)){
569  LogTrace("TestOutliers") << "SiStripRecHit2D";
570  clustersize = ((const SiStripRecHit2D*)(&**itHit))->cluster()->amplitudes().size() ;
571  hittypeval = 2;
572  }
573  else if (dynamic_cast<const SiStripMatchedRecHit2D*>(&**itHit)){
574  LogTrace("TestOutliers") << "SiStripMatchedRecHit2D";
575  int clsize1 = ((const SiStripMatchedRecHit2D*)(&**itHit))->monoHit()->cluster()->amplitudes().size();
576  int clsize2 = ((const SiStripMatchedRecHit2D*)(&**itHit))->stereoHit()->cluster()->amplitudes().size();
577  if (clsize1>clsize2) clustersize = clsize1;
578  else clustersize = clsize2;
579  hittypeval = 3;
580  }
581  else if (dynamic_cast<const ProjectedSiStripRecHit2D*>(&**itHit)){
582  LogTrace("TestOutliers") << "ProjectedSiStripRecHit2D";
583  clustersize = ((const ProjectedSiStripRecHit2D*)(&**itHit))->originalHit().cluster()->amplitudes().size();
584  hittypeval = 4;
585  }
586 
587  //find the layer of the hit
588  int subdetId = (*itHit)->geographicalId().subdetId();
589  int layerId = 0;
590  DetId id = (*itHit)->geographicalId();
591  if (id.subdetId()==3) layerId = ((TIBDetId)(id)).layer();
592  if (id.subdetId()==5) layerId = ((TOBDetId)(id)).layer();
593  if (id.subdetId()==1) layerId = ((PXBDetId)(id)).layer();
594  if (id.subdetId()==4) layerId = ((TIDDetId)(id)).wheel();
595  if (id.subdetId()==6) layerId = ((TECDetId)(id)).wheel();
596  if (id.subdetId()==2) layerId = ((PXFDetId)(id)).disk();
597  layerval = subdetId*10+layerId;
598 
599  //LogTrace("TestOutliers") << "gpos";
600  GlobalPoint gpos = theG->idToDet((*itHit)->geographicalId())->surface().toGlobal((*itHit)->localPosition());
601 
602 
603  //get the vector of sim hit associated and choose the one with the largest energy loss
604  //double delta = 99999;
605  //LocalPoint rhitLPv = (*itHit)->localPosition();
606  //vector<PSimHit> assSimHits = hitAssociator->associateHit(**itHit);
607  //if (assSimHits.size()==0) continue;
608  //PSimHit shit;
609  //for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){
610  //if ((m->localPosition()-rhitLPv).mag()<delta) {
611  // shit=*m;
612  // delta = (m->localPosition()-rhitLPv).mag();
613  // }
614  //}
615  //LogTrace("TestOutliers") << "energyLoss_";
616  double energyLoss_ = 0.;
617  unsigned int monoId = 0;
618  std::vector<double> energyLossM;
619  std::vector<double> energyLossS;
620  std::vector<PSimHit> assSimHits = hitAssociator->associateHit(**itHit);
621  if (assSimHits.size()==0) continue;
622  PSimHit shit;
623  std::vector<unsigned int> trackIds;
624  energyLossS.clear();
625  energyLossM.clear();
626  //LogTrace("TestOutliers") << "energyLossM";
627  for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){
628  if (outlier) energyLoss->Fill(m->energyLoss());
629  unsigned int tId = m->trackId();
630  if (find(trackIds.begin(),trackIds.end(),tId)==trackIds.end()) trackIds.push_back(tId);
631  LogTrace("TestOutliers") << "id=" << tId;
632  if (m->energyLoss()>energyLoss_) {
633  shit=*m;
634  energyLoss_ = m->energyLoss();
635  }
636  if (hittypeval==3) {
637  if (monoId==0) monoId = m->detUnitId();
638  if (monoId == m->detUnitId()){
639  energyLossM.push_back(m->energyLoss());
640  }
641  else {
642  energyLossS.push_back(m->energyLoss());
643  }
644  //std::cout << "detUnitId=" << m->detUnitId() << " trackId=" << m->trackId() << " energyLoss=" << m->energyLoss() << std::endl;
645  } else {
646  energyLossM.push_back(m->energyLoss());
647  }
648  }
649  unsigned int nIds = trackIds.size();
650 
651  if (outlier) {
652  goodbadhits->Fill(goodhit);
653  posxy->Fill(fabs(gpos.x()),fabs(gpos.y()));
654  poszr->Fill(fabs(gpos.z()),sqrt(gpos.x()*gpos.x()+gpos.y()*gpos.y()));
655  process->Fill(shit.processType());
656  energyLossMax->Fill(energyLoss_);
657  }
658 
659  //look if the hit is shared and if is produced only by ionization processes
660  bool shared = true;
661  bool ioniOnly = true;
662  unsigned int idc = 0;
663  for (size_t jj=0; jj<tpids.size(); jj++){
664  idc += std::count(trackIds.begin(),trackIds.end(),tpids[jj]);
665  }
666  if (idc==trackIds.size()) {
667  shared = false;
668  }
669  for(std::vector<PSimHit>::const_iterator m=assSimHits.begin()+1; m<assSimHits.end(); m++){
670  if ((m->processType()!=7&&m->processType()!=8&&m->processType()!=9)&&abs(m->particleType())!=11){
671  ioniOnly = false;
672  break;
673  }
674  }
675  if (ioniOnly&&!shared){
676  LogTrace("TestOutliers") << "delta";
677  }
678 
679  if (goodhit) {
680  if (outlier) {
681  goodprocess->Fill(shit.processType());
682  if (clustersize>=4) {
683  goodhittype_clgteq4->Fill(hittypeval);
684  goodlayer_clgteq4->Fill(layerval);
685  } else {
686  goodhittype_cllt4->Fill(hittypeval);
687  goodlayer_cllt4->Fill(layerval);
688  }
689  //LogTrace("TestOutliers") << "hittypeval && clustersize";
690  if (hittypeval==1 && clustersize>=4) goodpixgteq4_simvecsize->Fill(assSimHits.size());
691  if (hittypeval==1 && clustersize<4 ) goodpixlt4_simvecsize->Fill(assSimHits.size());
692  if (hittypeval==2 && clustersize>=4) goodst1gteq4_simvecsize->Fill(assSimHits.size());
693  if (hittypeval==2 && clustersize<4 ) goodst1lt4_simvecsize->Fill(assSimHits.size());
694  if (hittypeval==3 && clustersize>=4) goodst2gteq4_simvecsize->Fill(assSimHits.size());
695  if (hittypeval==3 && clustersize<4 ) goodst2lt4_simvecsize->Fill(assSimHits.size());
696  if (hittypeval==4 && clustersize>=4) goodprjgteq4_simvecsize->Fill(assSimHits.size());
697  if (hittypeval==4 && clustersize<4 ) goodprjlt4_simvecsize->Fill(assSimHits.size());
698 
699  //LogTrace("TestOutliers") << "hittypeval";
700  if (hittypeval==1) goodpix_clustersize->Fill(clustersize);
701  if (hittypeval==2) goodst1_clustersize->Fill(clustersize);
702  if (hittypeval==3) goodst2_clustersize->Fill(clustersize);
703  if (hittypeval==4) goodprj_clustersize->Fill(clustersize);
704  if (hittypeval==1) goodpix_simvecsize->Fill(assSimHits.size());
705  if (hittypeval==2) goodst1_simvecsize->Fill(assSimHits.size());
706  if (hittypeval==3) goodst2_simvecsize->Fill(assSimHits.size());
707  if (hittypeval==4) goodprj_simvecsize->Fill(assSimHits.size());
708 
709  //LogTrace("TestOutliers") << "nOfTrackIds";
710  nOfTrackIds->Fill(nIds);
711  if (hittypeval!=3) {
712  if (energyLossM.size()>1) {
713  sort(energyLossM.begin(),energyLossM.end(),greater<double>());
714  energyLossRatio->Fill(energyLossM[1]/energyLossM[0]);
715  }
716  } else {
717  //LogTrace("TestOutliers") << "hittypeval==3";
718  if (energyLossM.size()>1&&energyLossS.size()<=1) {
719  sort(energyLossM.begin(),energyLossM.end(),greater<double>());
720  energyLossRatio->Fill(energyLossM[1]/energyLossM[0]);
721  }
722  else if (energyLossS.size()>1&&energyLossM.size()<=1) {
723  sort(energyLossS.begin(),energyLossS.end(),greater<double>());
724  energyLossRatio->Fill(energyLossS[1]/energyLossS[0]);
725  }
726  else if (energyLossS.size()>1&&energyLossM.size()>1) {
727  sort(energyLossM.begin(),energyLossM.end(),greater<double>());
728  sort(energyLossS.begin(),energyLossS.end(),greater<double>());
729  energyLossM[1]/energyLossM[0] > energyLossS[1]/energyLossS[0]
730  ? energyLossRatio->Fill(energyLossM[1]/energyLossM[0])
731  : energyLossRatio->Fill(energyLossS[1]/energyLossS[0]);
732  }
733  }
734 
735  LogTrace("TestOutliers") << "before merged";
736  const SiStripMatchedRecHit2D* tmp = dynamic_cast<const SiStripMatchedRecHit2D*>(&**itHit);
737  LogTrace("TestOutliers") << "tmp=" << tmp;
738  LogTrace("TestOutliers") << "assSimHits.size()=" << assSimHits.size();
739  if ( (assSimHits.size()>1 && tmp==0) ||
740  (assSimHits.size()>2 && tmp!=0) ) {
741  //std::cout << "MERGED HIT" << std::endl;
742  //LogTrace("TestOutliers") << "merged";
743  mergedlayer->Fill(layerval);
744  mergedcluster->Fill(clustersize);
745  mergedhittype->Fill(hittypeval);
746 
747  for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){
748  unsigned int tId = m->trackId();
749  LogTrace("TestOutliers") << "component with id=" << tId << " eLoss=" << m->energyLoss()
750  << " proc=" << m->processType() << " part=" << m->particleType();
751  if (find(tpids.begin(),tpids.end(),tId)==tpids.end()) continue;
752  if (m->processType()==2) {
753  //GlobalPoint gpos = theG->idToDet((*itHit)->geographicalId())->surface().toGlobal((*itHit)->localPosition());
754  //GlobalPoint gpr = rhit->globalPosition();
755  AlgebraicSymMatrix ger =
756  ErrorFrameTransformer().transform((*itHit)->localPositionError(),theG->idToDet((*itHit)->geographicalId())->surface() ).matrix();
757  //AlgebraicSymMatrix ger = rhit->globalPositionError().matrix();
758  GlobalPoint gps = theG->idToDet((*itHit)->geographicalId())->surface().toGlobal(m->localPosition());
759  //LogTrace("TestOutliers") << gpr << " " << gps << " " << ger;
760  CLHEP::HepVector delta(3);
761  delta[0]=gpos.x()-gps.x();
762  delta[1]=gpos.y()-gps.y();
763  delta[2]=gpos.z()-gps.z();
764  LogTrace("TestOutliers") << delta << " " << ger ;
765  double mpull = sqrt(delta[0]*delta[0]/ger[0][0]+delta[1]*delta[1]/ger[1][1]+delta[2]*delta[2]/ger[2][2]);
766  LogTrace("TestOutliers") << "hit pull=" << mpull;//ger.similarity(delta);
767  mergedPull->Fill(mpull);
768  }
769  }
770  LogTrace("TestOutliers") << "end merged";
771  } else {//not merged=>good
772  //LogTrace("TestOutliers") << "goodlayer";
773  goodlayer->Fill(layerval);
774  goodcluster->Fill(clustersize);
775  goodhittype->Fill(hittypeval);
776  }
777  }//if (outlier)
778 
779  const SiPixelRecHit* pix = dynamic_cast<const SiPixelRecHit*>(&**itHit);
780  if ((hittypeval!=3 && assSimHits.size()<2)||(hittypeval==3 && assSimHits.size()<3)){
781  if (outlier) {
782  goodhittype_simvecsmall->Fill(hittypeval);
783  goodlayer_simvecsmall->Fill(layerval);
784  goodbadmerged->Fill(1);
785  if (pix) {
786  probXgood->Fill(pix->probabilityX());
787  probYgood->Fill(pix->probabilityY());
788  }
789  LogTrace("TestOutliers") << "out good";
790  }
791  else if (lost){
792  goodbadmergedLost->Fill(1);
793  LogTrace("TestOutliers") << "lost good";
794  }
795  else if (gained){
796  goodbadmergedGained->Fill(1);
797  LogTrace("TestOutliers") << "gained good";
798  }
799  LogTrace("TestOutliers") << "good";
800  } else {
801  if (outlier) {
802  goodhittype_simvecbig->Fill(hittypeval);
803  goodlayer_simvecbig->Fill(layerval);
804  if (ioniOnly&&!shared) {
805  goodbadmerged->Fill(3);
806  if (pix) {
807  probXdelta->Fill(pix->probabilityX());
808  probYdelta->Fill(pix->probabilityY());
809  }
810  } else if(!ioniOnly&&!shared) {
811  goodbadmerged->Fill(4);
812  if (pix) {
813  probXnoshare->Fill(pix->probabilityX());
814  probYnoshare->Fill(pix->probabilityY());
815  }
816  } else {
817  goodbadmerged->Fill(5);
818  if (pix) {
819  probXshared->Fill(pix->probabilityX());
820  probYshared->Fill(pix->probabilityY());
821  }
822  }
823  LogTrace("TestOutliers") << "out merged, ioniOnly=" << ioniOnly << " shared=" << shared;
824  }
825  else if (lost) {
826  if (ioniOnly&&!shared) goodbadmergedLost->Fill(3);
827  else if(!ioniOnly&&!shared) goodbadmergedLost->Fill(4);
828  else goodbadmergedLost->Fill(5);
829  LogTrace("TestOutliers") << "lost merged, ioniOnly=" << ioniOnly << " shared=" << shared;
830  }
831  else if (gained) {
832  if (ioniOnly&&!shared) goodbadmergedGained->Fill(3);
833  else if(!ioniOnly&&!shared) goodbadmergedGained->Fill(4);
834  else goodbadmergedGained->Fill(5);
835  LogTrace("TestOutliers") << "gained merged, ioniOnly=" << ioniOnly << " shared=" << shared;
836  }
837  LogTrace("TestOutliers") << "merged, ioniOnly=" << ioniOnly << " shared=" << shared;
838  }
839  } //if (goodhit)
840  else {//badhit
841  //LogTrace("TestOutliers") << "badhit";
842  if (outlier) {
843  badcluster->Fill(clustersize);
844  badhittype->Fill(hittypeval);
845  badlayer->Fill(layerval);
846  badprocess->Fill(shit.processType());
847  goodbadmerged->Fill(2);
848  const SiPixelRecHit* pix = dynamic_cast<const SiPixelRecHit*>(&**itHit);
849  if (pix) {
850  probXbad->Fill(pix->probabilityX());
851  probYbad->Fill(pix->probabilityY());
852  }
853  LogTrace("TestOutliers") << "out bad";
854  }
855  else if (lost) {
856  goodbadmergedLost->Fill(2);
857  LogTrace("TestOutliers") << "lost bad";
858  }
859  else if (gained) {
860  goodbadmergedGained->Fill(2);
861  LogTrace("TestOutliers") << "gained bad";
862  }
863  LogTrace("TestOutliers") << "bad";
864  }
865  }
866  }
867  }
868  //else if ( trackOut->numberOfValidHits() > trackOld->numberOfValidHits() ) {
869  else if ( 0 ) {
870  LogTrace("TestOutliers") << "outliers for track with #hits=" << trackOut->numberOfValidHits();
871  tracks->Fill(1);
872  LogTrace("TestOutliers") << "Out->pt=" << trackOut->pt() << " Old->pt=" << trackOld->pt()
873  << " tp->pt=" << sqrt(tpr->momentum().perp2())
874  //<< " trackOut->ptError()=" << trackOut->ptError() << " trackOld->ptError()=" << trackOld->ptError()
875  << " Old->validHits=" << trackOld->numberOfValidHits() << " Out->validHits=" << trackOut->numberOfValidHits()
876  /*<< " fracOld=" << fracOld*/ << " fracOut=" << fracOut
877  << " deltaHits=" << trackOld->numberOfValidHits()-trackOut->numberOfValidHits();
878  LogTrace("TestOutliers") << "track with gained hits";
879  gainedhits2->Fill(trackOut->numberOfValidHits()-trackOld->numberOfValidHits());
880  } else {
881  LogTrace("TestOutliers") << "no outliers for track with #hits=" << trackOut->numberOfValidHits();
882  }
883  }
884  LogTrace("TestOutliers") << "end track old #" << j;
885  }
886 
887  for (unsigned int k=0;k<tracksOut->size();++k) {
888  if ( find(outused.begin(),outused.end(),k)==outused.end() ) {
889  edm::RefToBase<reco::Track> trackOut(tracksOut, k);
890  bool outAssoc = recSimCollOut.find(trackOut) != recSimCollOut.end();
891  if (outAssoc) {
892  hitsPerTrackGained->Fill(trackOut->numberOfValidHits());
893  LogTrace("TestOutliers") << "gained track out id=" << k << " #hits==" << trackOut->numberOfValidHits();
894 
895  }
896  }
897  }
898  delete hitAssociator;
899  delete theAssociatorOld;
900  delete theAssociatorOut;
901 }
TestOutliers::histoLambdaOut
TH1F * histoLambdaOut
Definition: TestOutliers.cc:89
reco::TrackBase::qoverp
double qoverp() const
q/p
Definition: TrackBase.h:115
delta
dbl * delta
Definition: mlp_gen.cc:36
TestOutliers::hitsPerTrackGained
TH1F * hitsPerTrackGained
Definition: TestOutliers.cc:90
reco::TrackBase::d0Error
double d0Error() const
error on d0
Definition: TrackBase.h:211
TestOutliers::goodhittype_cllt4
TH1F * goodhittype_cllt4
Definition: TestOutliers.cc:93
TestOutliers::mergedcluster
TH1F * mergedcluster
Definition: TestOutliers.cc:108
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:66
TestOutliers::histoDzOld
TH1F * histoDzOld
Definition: TestOutliers.cc:88
TestOutliers::goodprj_clustersize
TH1F * goodprj_clustersize
Definition: TestOutliers.cc:102
TestOutliers::probXshared
TH1F * probXshared
Definition: TestOutliers.cc:112
reco::TrackBase::d0
double d0() const
dxy parameter in perigee convention (d0 = - dxy)
Definition: TrackBase.h:123
TestOutliers::goodpixlt4_simvecsize
TH1F * goodpixlt4_simvecsize
Definition: TestOutliers.cc:95
ErrorFrameTransformer::transform
GlobalError transform(const LocalError &le, const Surface &surf) const
Definition: ErrorFrameTransformer.h:13
TestOutliers::goodhittype
TH1F * goodhittype
Definition: TestOutliers.cc:93
TestOutliers::histoQoverpOut
TH1F * histoQoverpOut
Definition: TestOutliers.cc:85
TestOutliers::hitsPerTrackAssocLost
TH1F * hitsPerTrackAssocLost
Definition: TestOutliers.cc:90
TrajectoryStateClosestToPoint::theState
const FreeTrajectoryState & theState() const
Definition: TrajectoryStateClosestToPoint.h:124
TestOutliers::goodst1_simvecsize
TH1F * goodst1_simvecsize
Definition: TestOutliers.cc:104
TestOutliers::goodst2lt4_simvecsize
TH1F * goodst2lt4_simvecsize
Definition: TestOutliers.cc:97
TestOutliers::goodbadmergedLost
TH1F * goodbadmergedLost
Definition: TestOutliers.cc:107
TestOutliers::goodst1gteq4_simvecsize
TH1F * goodst1gteq4_simvecsize
Definition: TestOutliers.cc:96
TestOutliers::goodprj_simvecsize
TH1F * goodprj_simvecsize
Definition: TestOutliers.cc:106
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:63
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:10
TrackAssociatorByHits_cfi.TrackAssociatorByHits
tuple TrackAssociatorByHits
Definition: TrackAssociatorByHits_cfi.py:3
PV3DBase::y
T y() const
Definition: PV3DBase.h:57
abs
#define abs(x)
Definition: mlp_lapack.h:159
TestOutliers::deltahitsNO
TH1F * deltahitsNO
Definition: TestOutliers.cc:90
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:139
TestOutliers::psetold
edm::ParameterSet psetold
Definition: TestOutliers.cc:114
TestOutliers::energyLossRatio
TH1F * energyLossRatio
Definition: TestOutliers.cc:108
TestOutliers::histoPtOut
TH1F * histoPtOut
Definition: TestOutliers.cc:84
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:21
TrackerHitAssociator::associateHitId
std::vector< SimHitIdpr > associateHitId(const TrackingRecHit &thit)
Definition: TrackerHitAssociator.cc:280
TestOutliers::hitsPerTrackAssocGained
TH1F * hitsPerTrackAssocGained
Definition: TestOutliers.cc:90
TestOutliers::process
TH1F * process
Definition: TestOutliers.cc:92
TestOutliers::goodst2_simvecsize
TH1F * goodst2_simvecsize
Definition: TestOutliers.cc:105
TestOutliers::probYdelta
TH1F * probYdelta
Definition: TestOutliers.cc:113
TestOutliers::energyLossMax
TH1F * energyLossMax
Definition: TestOutliers.cc:108
TestOutliers::probYgood
TH1F * probYgood
Definition: TestOutliers.cc:113
TestOutliers::badhittype
TH1F * badhittype
Definition: TestOutliers.cc:93
PV3DBase::theta
Geom::Theta< T > theta() const
Definition: PV3DBase.h:69
TestOutliers::goodlayer_clgteq4
TH1F * goodlayer_clgteq4
Definition: TestOutliers.cc:93
TestOutliers::mergedhittype
TH1F * mergedhittype
Definition: TestOutliers.cc:108
TrackCharge
int TrackCharge
Definition: TrackCharge.h:4
TestOutliers::energyLoss
TH1F * energyLoss
Definition: TestOutliers.cc:108
TestOutliers::deltahits
TH1F * deltahits
Definition: TestOutliers.cc:90
TestOutliers::histoD0Old
TH1F * histoD0Old
Definition: TestOutliers.cc:87
TestOutliers::probYshared
TH1F * probYshared
Definition: TestOutliers.cc:113
TestOutliers::goodprocess
TH1F * goodprocess
Definition: TestOutliers.cc:92
TestOutliers::goodhittype_simvecsmall
TH1F * goodhittype_simvecsmall
Definition: TestOutliers.cc:107
TestOutliers::goodpix_clustersize
TH1F * goodpix_clustersize
Definition: TestOutliers.cc:99
TestOutliers::goodprjgteq4_simvecsize
TH1F * goodprjgteq4_simvecsize
Definition: TestOutliers.cc:98
PV3DBase::mag
T mag() const
Definition: PV3DBase.h:61
TestOutliers::badlayer
TH1F * badlayer
Definition: TestOutliers.cc:93
TestOutliers::goodst1_clustersize
TH1F * goodst1_clustersize
Definition: TestOutliers.cc:100
TestOutliers::goodpixgteq4_simvecsize
TH1F * goodpixgteq4_simvecsize
Definition: TestOutliers.cc:95
TestOutliers::histoPhiOld
TH1F * histoPhiOld
Definition: TestOutliers.cc:86
TestOutliers::badcluster
TH1F * badcluster
Definition: TestOutliers.cc:92
TestOutliers::probXgood
TH1F * probXgood
Definition: TestOutliers.cc:112
TestOutliers::countOutA
TH1F * countOutA
Definition: TestOutliers.cc:110
TestOutliers::goodprjlt4_simvecsize
TH1F * goodprjlt4_simvecsize
Definition: TestOutliers.cc:98
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:131
PV3DBase::z
T z() const
Definition: PV3DBase.h:58
TestOutliers::goodpix_simvecsize
TH1F * goodpix_simvecsize
Definition: TestOutliers.cc:103
TestOutliers::goodlayer
TH1F * goodlayer
Definition: TestOutliers.cc:93
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
TestOutliers::badprocess
TH1F * badprocess
Definition: TestOutliers.cc:92
reco::TrackBase::ptError
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:194
TestOutliers::histoPhiOut
TH1F * histoPhiOut
Definition: TestOutliers.cc:86
reco::TrackBase::phiError
double phiError() const
error on phi
Definition: TrackBase.h:207
j
int j
Definition: DBlmapReader.cc:9
reco::TrackBase::lambda
double lambda() const
Lambda angle.
Definition: TrackBase.h:119
TestOutliers::goodlayer_simvecsmall
TH1F * goodlayer_simvecsmall
Definition: TestOutliers.cc:107
TestOutliers::psetout
edm::ParameterSet psetout
Definition: TestOutliers.cc:114
TestOutliers::goodhittype_simvecbig
TH1F * goodhittype_simvecbig
Definition: TestOutliers.cc:107
SiPixelRecHit::probabilityY
float probabilityY() const
Definition: SiPixelRecHit.h:82
reco::TrackBase::numberOfValidHits
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:232
reco::Track::recHitsBegin
trackingRecHit_iterator recHitsBegin() const
Iterator to first hit on the track.
Definition: Track.h:63
first
bool first
Definition: L1TdeRCT.cc:79
TestOutliers::trackTagsOut_
edm::InputTag trackTagsOut_
Definition: TestOutliers.cc:74
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:355
TestOutliers::probXnoshare
TH1F * probXnoshare
Definition: TestOutliers.cc:112
FreeTrajectoryState::momentum
GlobalVector momentum() const
Definition: FreeTrajectoryState.h:81
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:517
TrackAssociatorBase::associateRecoToSim
virtual reco::RecoToSimCollection associateRecoToSim(edm::Handle< edm::View< reco::Track > > &tCH, edm::Handle< TrackingParticleCollection > &tPCH, const edm::Event *event=0, const edm::EventSetup *setup=0) const
compare reco to sim the handle of reco::Track and TrackingParticle collections
Definition: TrackAssociatorBase.h:44
TestOutliers::probYnoshare
TH1F * probYnoshare
Definition: TestOutliers.cc:113
gen::k
int k[5][pyjets_maxn]
Definition: Pythia6Hadronizer.cc:86
reco::TrackBase::qoverpError
double qoverpError() const
error on signed transverse curvature
Definition: TrackBase.h:192
TestOutliers::mergedlayer
TH1F * mergedlayer
Definition: TestOutliers.cc:108
reco::TrackBase::dz
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:127
reco::TrackBase::dzError
double dzError() const
error on dz
Definition: TrackBase.h:215
TrackingParticle::g4t_iterator
std::vector< SimTrack >::const_iterator g4t_iterator
Definition: TrackingParticle.h:36
TestOutliers::deltahitsOK
TH1F * deltahitsOK
Definition: TestOutliers.cc:90
DetId
Definition: DetId.h:20
FreeTrajectoryState::position
GlobalPoint position() const
Definition: FreeTrajectoryState.h:78
TestOutliers::goodst1lt4_simvecsize
TH1F * goodst1lt4_simvecsize
Definition: TestOutliers.cc:96
TestOutliers::histoPtOld
TH1F * histoPtOld
Definition: TestOutliers.cc:84
TestOutliers::goodhittype_clgteq4
TH1F * goodhittype_clgteq4
Definition: TestOutliers.cc:93
reco::Track::seedRef
edm::RefToBase< TrajectorySeed > seedRef() const
Definition: Track.h:112
M_PI
#define M_PI
Definition: BFit3D.cc:3
TestOutliers::tpTags_
edm::InputTag tpTags_
Definition: TestOutliers.cc:76
TestOutliers::deltahitsAssocGained
TH1F * deltahitsAssocGained
Definition: TestOutliers.cc:90
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
TestOutliers::theAssociatorOld
TrackAssociatorBase * theAssociatorOld
Definition: TestOutliers.cc:77
TestOutliers::goodlayer_cllt4
TH1F * goodlayer_cllt4
Definition: TestOutliers.cc:93
L1TEmulatorMonitor_cff.p
tuple p
Definition: L1TEmulatorMonitor_cff.py:48
TestOutliers::histoDzOut
TH1F * histoDzOut
Definition: TestOutliers.cc:88
TestOutliers::theAssociatorOut
TrackAssociatorBase * theAssociatorOut
Definition: TestOutliers.cc:78
TestOutliers::probXdelta
TH1F * probXdelta
Definition: TestOutliers.cc:112
TestOutliers::goodst2_clustersize
TH1F * goodst2_clustersize
Definition: TestOutliers.cc:101
PSimHit::processType
unsigned short processType() const
Definition: PSimHit.h:118
tmp
std::vector< std::vector< double > > tmp
Definition: MVATrainer.cc:100
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
TestOutliers::histoLambdaOld
TH1F * histoLambdaOld
Definition: TestOutliers.cc:89
TestOutliers::countOldT
TH1F * countOldT
Definition: TestOutliers.cc:110
TestOutliers::countOutT
TH1F * countOutT
Definition: TestOutliers.cc:110
CoreSimTrack::momentum
const math::XYZTLorentzVectorD & momentum() const
particle info...
Definition: CoreSimTrack.h:36
TestOutliers::goodst2gteq4_simvecsize
TH1F * goodst2gteq4_simvecsize
Definition: TestOutliers.cc:97
TestOutliers::goodcluster
TH1F * goodcluster
Definition: TestOutliers.cc:92
TestOutliers::histoD0Out
TH1F * histoD0Out
Definition: TestOutliers.cc:87
TestOutliers::theG
edm::ESHandle< TrackerGeometry > theG
Definition: TestOutliers.cc:81
AlgebraicSymMatrix
CLHEP::HepSymMatrix AlgebraicSymMatrix
Definition: AlgebraicObjects.h:14
TrackerHitAssociator::associateHit
std::vector< PSimHit > associateHit(const TrackingRecHit &thit)
Definition: TrackerHitAssociator.cc:120
TestOutliers::mergedPull
TH1F * mergedPull
Definition: TestOutliers.cc:108
TestOutliers::goodbadmerged
TH1F * goodbadmerged
Definition: TestOutliers.cc:107
TestOutliers::deltahitsAssocLost
TH1F * deltahitsAssocLost
Definition: TestOutliers.cc:90
TestOutliers::goodlayer_simvecbig
TH1F * goodlayer_simvecbig
Definition: TestOutliers.cc:107
TestOutliers::histoQoverpOld
TH1F * histoQoverpOld
Definition: TestOutliers.cc:85
SiPixelRecHit::probabilityX
float probabilityX() const
Definition: SiPixelRecHit.h:79
TestOutliers::nOfTrackIds
TH1F * nOfTrackIds
Definition: TestOutliers.cc:108
TestOutliers::trackTagsOld_
edm::InputTag trackTagsOld_
Definition: TestOutliers.cc:75
TestOutliers::hitsPerTrackLost
TH1F * hitsPerTrackLost
Definition: TestOutliers.cc:90
TestOutliers::gainedhits
TH1F * gainedhits
Definition: TestOutliers.cc:111
edm::Ref::get
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:241
TestOutliers::hitAssociator
TrackerHitAssociator * hitAssociator
Definition: TestOutliers.cc:80
PV3DBase::x
T x() const
Definition: PV3DBase.h:56
edm::RefToBase::get
value_type const * get() const
Definition: RefToBase.h:212
v
mathSSE::Vec4< T > v
Definition: newBasic3DVector.h:343
TestOutliers::goodbadhits
TH1F * goodbadhits
Definition: TestOutliers.cc:92
GlobalVector
Global3DVector GlobalVector
Definition: GlobalVector.h:10
reco::TrackBase::thetaError
double thetaError() const
error on theta
Definition: TrackBase.h:201
TestOutliers::goodbadmergedGained
TH1F * goodbadmergedGained
Definition: TestOutliers.cc:107
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:27
TestOutliers::gainedhits2
TH1F * gainedhits2
Definition: TestOutliers.cc:111
reco::Track::recHitsEnd
trackingRecHit_iterator recHitsEnd() const
Iterator to last hit on the track.
Definition: Track.h:65
void TestOutliers::beginRun ( edm::Run &  run,
const edm::EventSetup &  es 
)
privatevirtual

Definition at line 906 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.

907 {
908  es.get<TrackerDigiGeometryRecord>().get(theG);
909  const bool oldAddDir = TH1::AddDirectoryStatus();
910  TH1::AddDirectory(true);
911  file = new TFile(out.c_str(),"recreate");
912  histoPtOut = new TH1F("histoPtOut","histoPtOut",100,-10,10);
913  histoPtOld = new TH1F("histoPtOld","histoPtOld",100,-10,10);
914  histoQoverpOut = new TH1F("histoQoverpOut","histoQoverpOut",250,-25,25);
915  histoQoverpOld = new TH1F("histoQoverpOld","histoQoverpOld",250,-25,25);
916  histoPhiOut = new TH1F("histoPhiOut","histoPhiOut",250,-25,25);
917  histoPhiOld = new TH1F("histoPhiOld","histoPhiOld",250,-25,25);
918  histoD0Out = new TH1F("histoD0Out","histoD0Out",250,-25,25);
919  histoD0Old = new TH1F("histoD0Old","histoD0Old",250,-25,25);
920  histoDzOut = new TH1F("histoDzOut","histoDzOut",250,-25,25);
921  histoDzOld = new TH1F("histoDzOld","histoDzOld",250,-25,25);
922  histoLambdaOut = new TH1F("histoLambdaOut","histoLambdaOut",250,-25,25);
923  histoLambdaOld = new TH1F("histoLambdaOld","histoLambdaOld",250,-25,25);
924  deltahits = new TH1F("deltahits","deltahits",80,-40,40);
925  deltahitsOK = new TH1F("deltahitsOK","deltahitsOK",20,0,20);
926  deltahitsNO = new TH1F("deltahitsNO","deltahitsNO",20,0,20);
927  okcutsOut = new TH1F("okcutsOut","okcutsOut",2,-0.5,1.5);
928  okcutsOld = new TH1F("okcutsOld","okcutsOld",2,-0.5,1.5);
929  tracks = new TH1F("tracks_","tracks_",2,-0.5,1.5);
930  goodbadhits = new TH1F("goodbadhits","goodbadhits",2,-0.5,1.5);
931  process = new TH1F("process","process",20,-0.5,19.5);
932  goodcluster = new TH1F("goodcluster","goodcluster",40,-0.5,39.5);
933  goodprocess = new TH1F("goodprocess","goodprocess",20,-0.5,19.5);
934  badcluster = new TH1F("badcluster","badcluster",40,-0.5,39.5);
935  badprocess = new TH1F("badprocess","badprocess",20,-0.5,19.5);
936  goodhittype = new TH1F("goodhittype","goodhittype",5,-0.5,4.5);
937  goodlayer = new TH1F("goodlayer","goodlayer",70,-0.5,69.5);
938  goodhittype_clgteq4 = new TH1F("goodhittype_clgteq4","goodhittype_clgteq4",5,-0.5,4.5);
939  goodlayer_clgteq4 = new TH1F("goodlayer_clgteq4","goodlayer_clgteq4",70,-0.5,69.5);
940  goodhittype_cllt4 = new TH1F("goodhittype_cllt4","goodhittype_cllt4",5,-0.5,4.5);
941  goodlayer_cllt4 = new TH1F("goodlayer_cllt4","goodlayer_cllt4",70,-0.5,69.5);
942  badhittype = new TH1F("badhittype","badhittype",5,-0.5,4.5);
943  badlayer = new TH1F("badlayer","badlayer",70,-0.5,69.5);
944  posxy = new TH2F("posxy","posxy",1200,0,120,1200,0,120);
945  poszr = new TH2F("poszr","poszr",3000,0,300,1200,0,120);
946  goodpixgteq4_simvecsize = new TH1F("goodpixgteq4_simvecsize","goodpixgteq4_simvecsize",40,-0.5,39.5);
947  goodpixlt4_simvecsize = new TH1F("goodpixlt4_simvecsize","goodpixlt4_simvecsize",40,-0.5,39.5);
948  goodst1gteq4_simvecsize = new TH1F("goodst1gteq4_simvecsize","goodst1gteq4_simvecsize",40,-0.5,39.5);
949  goodst1lt4_simvecsize = new TH1F("goodst1lt4_simvecsize","goodst1lt4_simvecsize",40,-0.5,39.5);
950  goodst2gteq4_simvecsize = new TH1F("goodst2gteq4_simvecsize","goodst2gteq4_simvecsize",40,-0.5,39.5);
951  goodst2lt4_simvecsize = new TH1F("goodst2lt4_simvecsize","goodst2lt4_simvecsize",40,-0.5,39.5);
952  goodprjgteq4_simvecsize = new TH1F("goodprjgteq4_simvecsize","goodprjgteq4_simvecsize",40,-0.5,39.5);
953  goodprjlt4_simvecsize = new TH1F("goodprjlt4_simvecsize","goodprjlt4_simvecsize",40,-0.5,39.5);
954  goodpix_clustersize = new TH1F("goodpix_clustersize","goodpix_clustersize",40,-0.5,39.5);
955  goodst1_clustersize = new TH1F("goodst1_clustersize","goodst1_clustersize",40,-0.5,39.5);
956  goodst2_clustersize = new TH1F("goodst2_clustersize","goodst2_clustersize",40,-0.5,39.5);
957  goodprj_clustersize = new TH1F("goodprj_clustersize","goodprj_clustersize",40,-0.5,39.5);
958  goodpix_simvecsize = new TH1F("goodpix_simvecsize","goodpix_simvecsize",40,-0.5,39.5);
959  goodst1_simvecsize = new TH1F("goodst1_simvecsize","goodst1_simvecsize",40,-0.5,39.5);
960  goodst2_simvecsize = new TH1F("goodst2_simvecsize","goodst2_simvecsize",40,-0.5,39.5);
961  goodprj_simvecsize = new TH1F("goodprj_simvecsize","goodprj_simvecsize",40,-0.5,39.5);
962  goodhittype_simvecsmall = new TH1F("goodhittype_simvecsmall","goodhittype_simvecsmall",5,-0.5,4.5);
963  goodlayer_simvecsmall = new TH1F("goodlayer_simvecsmall","goodlayer_simvecsmall",70,-0.5,69.5);
964  goodhittype_simvecbig = new TH1F("goodhittype_simvecbig","goodhittype_simvecbig",5,-0.5,4.5);
965  goodlayer_simvecbig = new TH1F("goodlayer_simvecbig","goodlayer_simvecbig",70,-0.5,69.5);
966  goodbadmerged = new TH1F("goodbadmerged","goodbadmerged",5,0.5,5.5);
967  goodbadmergedLost = new TH1F("goodbadmergedLost","goodbadmergedLost",5,0.5,5.5);
968  goodbadmergedGained = new TH1F("goodbadmergedGained","goodbadmergedGained",5,0.5,5.5);
969  energyLoss = new TH1F("energyLoss","energyLoss",1000,0,0.1);
970  energyLossMax = new TH1F("energyLossMax","energyLossMax",1000,0,0.1);
971  energyLossRatio = new TH1F("energyLossRatio","energyLossRatio",100,0,1);
972  nOfTrackIds = new TH1F("nOfTrackIds","nOfTrackIds",10,0,10);
973  mergedPull = new TH1F("mergedPull","mergedPull",100,0,10);
974  mergedlayer = new TH1F("mergedlayer","mergedlayer",70,-0.5,69.5);
975  mergedhittype = new TH1F("mergedhittype","mergedhittype",5,-0.5,4.5);
976  mergedcluster = new TH1F("mergedcluster","mergedcluster",40,-0.5,39.5);
977  deltahitsAssocGained = new TH1F("deltahitsAssocGained","deltahitsAssocGained",80, -40, 40);
978  deltahitsAssocLost = new TH1F("deltahitsAssocLost","deltahitsAssocLost",80, -40, 40);
979  hitsPerTrackLost = new TH1F("hitsPerTrackLost","hitsPerTrackLost",40, -0.5, 39.5);
980  hitsPerTrackAssocLost = new TH1F("hitsPerTrackAssocLost","hitsPerTrackAssocLost",40, -0.5, 39.5);
981  hitsPerTrackGained = new TH1F("hitsPerTrackGained","hitsPerTrackGained",40, -0.5, 39.5);
982  hitsPerTrackAssocGained = new TH1F("hitsPerTrackAssocGained","hitsPerTrackAssocGained",40, -0.5, 39.5);
983  sizeOut = new TH1F("sizeOut","sizeOut",900,-0.5,899.5);
984  sizeOld = new TH1F("sizeOld","sizeOld",900,-0.5,899.5);
985  sizeOutT = new TH1F("sizeOutT","sizeOutT",900,-0.5,899.5);
986  sizeOldT = new TH1F("sizeOldT","sizeOldT",900,-0.5,899.5);
987  countOutA = new TH1F("countOutA","countOutA",2,0,2);
988  countOutT = new TH1F("countOutT","countOutT",2,0,2);
989  countOldT = new TH1F("countOldT","countOldT",2,0,2);
990  gainedhits = new TH1F("gainedhits","gainedhits",2,0,2);
991  gainedhits2 = new TH1F("gainedhits2","gainedhits2",30,-0.5,29.5);
992  probXgood = new TH1F("probXgood","probXgood",110,0,1.1);
993  probXbad = new TH1F("probXbad","probXbad",110,0,1.1);
994  probXdelta = new TH1F("probXdelta","probXdelta",110,0,1.1);
995  probXshared = new TH1F("probXshared","probXshared",110,0,1.1);
996  probXnoshare = new TH1F("probXnoshare","probXnoshare",110,0,1.1);
997  probYgood = new TH1F("probYgood","probYgood",110,0,1.1);
998  probYbad = new TH1F("probYbad","probYbad",110,0,1.1);
999  probYdelta = new TH1F("probYdelta","probYdelta",110,0,1.1);
1000  probYshared = new TH1F("probYshared","probYshared",110,0,1.1);
1001  probYnoshare = new TH1F("probYnoshare","probYnoshare",110,0,1.1);
1002  TH1::AddDirectory(oldAddDir);
1003 }
TestOutliers::histoLambdaOut
TH1F * histoLambdaOut
Definition: TestOutliers.cc:89
TestOutliers::hitsPerTrackGained
TH1F * hitsPerTrackGained
Definition: TestOutliers.cc:90
TestOutliers::goodhittype_cllt4
TH1F * goodhittype_cllt4
Definition: TestOutliers.cc:93
TestOutliers::mergedcluster
TH1F * mergedcluster
Definition: TestOutliers.cc:108
TestOutliers::histoDzOld
TH1F * histoDzOld
Definition: TestOutliers.cc:88
TestOutliers::goodprj_clustersize
TH1F * goodprj_clustersize
Definition: TestOutliers.cc:102
TestOutliers::probXshared
TH1F * probXshared
Definition: TestOutliers.cc:112
TestOutliers::goodpixlt4_simvecsize
TH1F * goodpixlt4_simvecsize
Definition: TestOutliers.cc:95
TestOutliers::goodhittype
TH1F * goodhittype
Definition: TestOutliers.cc:93
TestOutliers::histoQoverpOut
TH1F * histoQoverpOut
Definition: TestOutliers.cc:85
TestOutliers::hitsPerTrackAssocLost
TH1F * hitsPerTrackAssocLost
Definition: TestOutliers.cc:90
TestOutliers::goodst1_simvecsize
TH1F * goodst1_simvecsize
Definition: TestOutliers.cc:104
TestOutliers::goodst2lt4_simvecsize
TH1F * goodst2lt4_simvecsize
Definition: TestOutliers.cc:97
TestOutliers::goodbadmergedLost
TH1F * goodbadmergedLost
Definition: TestOutliers.cc:107
TestOutliers::goodst1gteq4_simvecsize
TH1F * goodst1gteq4_simvecsize
Definition: TestOutliers.cc:96
TestOutliers::goodprj_simvecsize
TH1F * goodprj_simvecsize
Definition: TestOutliers.cc:106
TestOutliers::okcutsOut
TH1F * okcutsOut
Definition: TestOutliers.cc:91
TestOutliers::file
TFile * file
Definition: TestOutliers.cc:83
TestOutliers::deltahitsNO
TH1F * deltahitsNO
Definition: TestOutliers.cc:90
TestOutliers::energyLossRatio
TH1F * energyLossRatio
Definition: TestOutliers.cc:108
TestOutliers::histoPtOut
TH1F * histoPtOut
Definition: TestOutliers.cc:84
TestOutliers::hitsPerTrackAssocGained
TH1F * hitsPerTrackAssocGained
Definition: TestOutliers.cc:90
TestOutliers::process
TH1F * process
Definition: TestOutliers.cc:92
TestOutliers::goodst2_simvecsize
TH1F * goodst2_simvecsize
Definition: TestOutliers.cc:105
TestOutliers::probYdelta
TH1F * probYdelta
Definition: TestOutliers.cc:113
TestOutliers::energyLossMax
TH1F * energyLossMax
Definition: TestOutliers.cc:108
TestOutliers::probYgood
TH1F * probYgood
Definition: TestOutliers.cc:113
TestOutliers::badhittype
TH1F * badhittype
Definition: TestOutliers.cc:93
TestOutliers::goodlayer_clgteq4
TH1F * goodlayer_clgteq4
Definition: TestOutliers.cc:93
TestOutliers::mergedhittype
TH1F * mergedhittype
Definition: TestOutliers.cc:108
TestOutliers::energyLoss
TH1F * energyLoss
Definition: TestOutliers.cc:108
TestOutliers::deltahits
TH1F * deltahits
Definition: TestOutliers.cc:90
TestOutliers::histoD0Old
TH1F * histoD0Old
Definition: TestOutliers.cc:87
TestOutliers::probYshared
TH1F * probYshared
Definition: TestOutliers.cc:113
TestOutliers::goodprocess
TH1F * goodprocess
Definition: TestOutliers.cc:92
TestOutliers::goodhittype_simvecsmall
TH1F * goodhittype_simvecsmall
Definition: TestOutliers.cc:107
TestOutliers::goodpix_clustersize
TH1F * goodpix_clustersize
Definition: TestOutliers.cc:99
TestOutliers::goodprjgteq4_simvecsize
TH1F * goodprjgteq4_simvecsize
Definition: TestOutliers.cc:98
TestOutliers::badlayer
TH1F * badlayer
Definition: TestOutliers.cc:93
TestOutliers::goodst1_clustersize
TH1F * goodst1_clustersize
Definition: TestOutliers.cc:100
TestOutliers::goodpixgteq4_simvecsize
TH1F * goodpixgteq4_simvecsize
Definition: TestOutliers.cc:95
TestOutliers::histoPhiOld
TH1F * histoPhiOld
Definition: TestOutliers.cc:86
TestOutliers::badcluster
TH1F * badcluster
Definition: TestOutliers.cc:92
TestOutliers::probXgood
TH1F * probXgood
Definition: TestOutliers.cc:112
TestOutliers::countOutA
TH1F * countOutA
Definition: TestOutliers.cc:110
TestOutliers::goodprjlt4_simvecsize
TH1F * goodprjlt4_simvecsize
Definition: TestOutliers.cc:98
TestOutliers::goodpix_simvecsize
TH1F * goodpix_simvecsize
Definition: TestOutliers.cc:103
TestOutliers::goodlayer
TH1F * goodlayer
Definition: TestOutliers.cc:93
TestOutliers::badprocess
TH1F * badprocess
Definition: TestOutliers.cc:92
TestOutliers::histoPhiOut
TH1F * histoPhiOut
Definition: TestOutliers.cc:86
TestOutliers::out
std::string out
Definition: TestOutliers.cc:82
TestOutliers::goodlayer_simvecsmall
TH1F * goodlayer_simvecsmall
Definition: TestOutliers.cc:107
TestOutliers::goodhittype_simvecbig
TH1F * goodhittype_simvecbig
Definition: TestOutliers.cc:107
TestOutliers::probXnoshare
TH1F * probXnoshare
Definition: TestOutliers.cc:112
TestOutliers::probYnoshare
TH1F * probYnoshare
Definition: TestOutliers.cc:113
TestOutliers::mergedlayer
TH1F * mergedlayer
Definition: TestOutliers.cc:108
TestOutliers::deltahitsOK
TH1F * deltahitsOK
Definition: TestOutliers.cc:90
TestOutliers::goodst1lt4_simvecsize
TH1F * goodst1lt4_simvecsize
Definition: TestOutliers.cc:96
TestOutliers::histoPtOld
TH1F * histoPtOld
Definition: TestOutliers.cc:84
TestOutliers::goodhittype_clgteq4
TH1F * goodhittype_clgteq4
Definition: TestOutliers.cc:93
TestOutliers::deltahitsAssocGained
TH1F * deltahitsAssocGained
Definition: TestOutliers.cc:90
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
TestOutliers::goodlayer_cllt4
TH1F * goodlayer_cllt4
Definition: TestOutliers.cc:93
TestOutliers::histoDzOut
TH1F * histoDzOut
Definition: TestOutliers.cc:88
TestOutliers::probXdelta
TH1F * probXdelta
Definition: TestOutliers.cc:112
TestOutliers::goodst2_clustersize
TH1F * goodst2_clustersize
Definition: TestOutliers.cc:101
TestOutliers::histoLambdaOld
TH1F * histoLambdaOld
Definition: TestOutliers.cc:89
TestOutliers::countOldT
TH1F * countOldT
Definition: TestOutliers.cc:110
TestOutliers::countOutT
TH1F * countOutT
Definition: TestOutliers.cc:110
TestOutliers::goodst2gteq4_simvecsize
TH1F * goodst2gteq4_simvecsize
Definition: TestOutliers.cc:97
TestOutliers::goodcluster
TH1F * goodcluster
Definition: TestOutliers.cc:92
TestOutliers::histoD0Out
TH1F * histoD0Out
Definition: TestOutliers.cc:87
TestOutliers::theG
edm::ESHandle< TrackerGeometry > theG
Definition: TestOutliers.cc:81
TestOutliers::okcutsOld
TH1F * okcutsOld
Definition: TestOutliers.cc:91
TestOutliers::mergedPull
TH1F * mergedPull
Definition: TestOutliers.cc:108
TestOutliers::goodbadmerged
TH1F * goodbadmerged
Definition: TestOutliers.cc:107
TestOutliers::deltahitsAssocLost
TH1F * deltahitsAssocLost
Definition: TestOutliers.cc:90
TestOutliers::goodlayer_simvecbig
TH1F * goodlayer_simvecbig
Definition: TestOutliers.cc:107
TestOutliers::histoQoverpOld
TH1F * histoQoverpOld
Definition: TestOutliers.cc:85
TestOutliers::nOfTrackIds
TH1F * nOfTrackIds
Definition: TestOutliers.cc:108
TestOutliers::hitsPerTrackLost
TH1F * hitsPerTrackLost
Definition: TestOutliers.cc:90
TestOutliers::gainedhits
TH1F * gainedhits
Definition: TestOutliers.cc:111
TestOutliers::goodbadhits
TH1F * goodbadhits
Definition: TestOutliers.cc:92
TestOutliers::goodbadmergedGained
TH1F * goodbadmergedGained
Definition: TestOutliers.cc:107
TestOutliers::gainedhits2
TH1F * gainedhits2
Definition: TestOutliers.cc:111
void TestOutliers::endJob ( void  )
privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1006 of file TestOutliers.cc.

References file, and LogTrace.

1006  {
1007  LogTrace("TestOutliers") << "TestOutliers::endJob";
1008  file->Write();
1009  LogTrace("TestOutliers") << "outfile written";
1010  file->Close();
1011  LogTrace("TestOutliers") << "oufile closed";
1012  LogTrace("TestOutliers") << "exiting TestOutliers::endJob";
1013 }
TestOutliers::file
TFile * file
Definition: TestOutliers.cc:83
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:517

Member Data Documentation

TH1F * TestOutliers::badcluster
private

Definition at line 92 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::badhittype
private

Definition at line 93 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::badlayer
private

Definition at line 93 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::badprocess
private

Definition at line 92 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::countOldT
private

Definition at line 110 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::countOutA
private

Definition at line 110 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::countOutT
private

Definition at line 110 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::deltahits
private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::deltahitsAssocGained
private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::deltahitsAssocLost
private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::deltahitsNO
private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::deltahitsOK
private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::energyLoss
private

Definition at line 108 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::energyLossMax
private

Definition at line 108 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::energyLossRatio
private

Definition at line 108 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TFile* TestOutliers::file
private

Definition at line 83 of file TestOutliers.cc.

Referenced by beginRun(), and endJob().

TH1F* TestOutliers::gainedhits
private

Definition at line 111 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::gainedhits2
private

Definition at line 111 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodbadhits
private

Definition at line 92 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodbadmerged
private

Definition at line 107 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodbadmergedGained
private

Definition at line 107 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodbadmergedLost
private

Definition at line 107 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodcluster
private

Definition at line 92 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodhittype
private

Definition at line 93 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodhittype_clgteq4
private

Definition at line 93 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodhittype_cllt4
private

Definition at line 93 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodhittype_simvecbig
private

Definition at line 107 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodhittype_simvecsmall
private

Definition at line 107 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer
private

Definition at line 93 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer_clgteq4
private

Definition at line 93 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer_cllt4
private

Definition at line 93 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer_simvecbig
private

Definition at line 107 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodlayer_simvecsmall
private

Definition at line 107 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodpix_clustersize
private

Definition at line 99 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodpix_simvecsize
private

Definition at line 103 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodpixgteq4_simvecsize
private

Definition at line 95 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodpixlt4_simvecsize
private

Definition at line 95 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodprj_clustersize
private

Definition at line 102 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodprj_simvecsize
private

Definition at line 106 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodprjgteq4_simvecsize
private

Definition at line 98 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodprjlt4_simvecsize
private

Definition at line 98 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodprocess
private

Definition at line 92 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst1_clustersize
private

Definition at line 100 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst1_simvecsize
private

Definition at line 104 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst1gteq4_simvecsize
private

Definition at line 96 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodst1lt4_simvecsize
private

Definition at line 96 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst2_clustersize
private

Definition at line 101 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst2_simvecsize
private

Definition at line 105 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::goodst2gteq4_simvecsize
private

Definition at line 97 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::goodst2lt4_simvecsize
private

Definition at line 97 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoD0Old
private

Definition at line 87 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoD0Out
private

Definition at line 87 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoDzOld
private

Definition at line 88 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoDzOut
private

Definition at line 88 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoLambdaOld
private

Definition at line 89 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoLambdaOut
private

Definition at line 89 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoPhiOld
private

Definition at line 86 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoPhiOut
private

Definition at line 86 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoPtOld
private

Definition at line 84 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoPtOut
private

Definition at line 84 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::histoQoverpOld
private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::histoQoverpOut
private

Definition at line 85 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TrackerHitAssociator* TestOutliers::hitAssociator
private

Definition at line 80 of file TestOutliers.cc.

Referenced by analyze().

TH1F * TestOutliers::hitsPerTrackAssocGained
private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::hitsPerTrackAssocLost
private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::hitsPerTrackGained
private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::hitsPerTrackLost
private

Definition at line 90 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::mergedcluster
private

Definition at line 108 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::mergedhittype
private

Definition at line 108 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::mergedlayer
private

Definition at line 108 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::mergedPull
private

Definition at line 108 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::nOfTrackIds
private

Definition at line 108 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::okcutsOld
private

Definition at line 91 of file TestOutliers.cc.

Referenced by beginRun().

TH1F* TestOutliers::okcutsOut
private

Definition at line 91 of file TestOutliers.cc.

Referenced by beginRun().

std::string TestOutliers::out
private

Definition at line 82 of file TestOutliers.cc.

Referenced by beginRun().

TH2F* TestOutliers::posxy
private

Definition at line 94 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH2F * TestOutliers::poszr
private

Definition at line 94 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probXbad
private

Definition at line 112 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probXdelta
private

Definition at line 112 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::probXgood
private

Definition at line 112 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probXnoshare
private

Definition at line 112 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probXshared
private

Definition at line 112 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probYbad
private

Definition at line 113 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probYdelta
private

Definition at line 113 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::probYgood
private

Definition at line 113 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probYnoshare
private

Definition at line 113 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::probYshared
private

Definition at line 113 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::process
private

Definition at line 92 of file TestOutliers.cc.

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

edm::ParameterSet TestOutliers::psetold
private

Definition at line 114 of file TestOutliers.cc.

Referenced by analyze(), and TestOutliers().

edm::ParameterSet TestOutliers::psetout
private

Definition at line 114 of file TestOutliers.cc.

Referenced by analyze(), and TestOutliers().

TH1F * TestOutliers::sizeOld
private

Definition at line 109 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::sizeOldT
private

Definition at line 109 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F* TestOutliers::sizeOut
private

Definition at line 109 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TH1F * TestOutliers::sizeOutT
private

Definition at line 109 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

TrackAssociatorBase* TestOutliers::theAssociatorOld
private

Definition at line 77 of file TestOutliers.cc.

Referenced by analyze().

TrackAssociatorBase* TestOutliers::theAssociatorOut
private

Definition at line 78 of file TestOutliers.cc.

Referenced by analyze().

edm::ESHandle<TrackerGeometry> TestOutliers::theG
private

Definition at line 81 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

edm::InputTag TestOutliers::tpTags_
private

Definition at line 76 of file TestOutliers.cc.

Referenced by analyze().

TH1F* TestOutliers::tracks
private

Definition at line 92 of file TestOutliers.cc.

Referenced by analyze(), and beginRun().

edm::InputTag TestOutliers::trackTagsOld_
private

Definition at line 75 of file TestOutliers.cc.

Referenced by analyze().

edm::InputTag TestOutliers::trackTagsOut_
private

Definition at line 74 of file TestOutliers.cc.

Referenced by analyze().

Generated for CMSSW Reference Manual by doxygen 1.8.5