CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
PhotonConversionTrajectorySeedProducerFromSingleLegAlgo.h
Go to the documentation of this file.
1 #ifndef PhotonConversionFinderFromTracks_H
2 #define PhotonConversionFinderFromTracks_H
3 
4 #include "FWCore/Framework/interface/Frameworkfwd.h"
5 #include "FWCore/Framework/interface/EDAnalyzer.h"
6 #include "FWCore/Framework/interface/Event.h"
7 #include "FWCore/ParameterSet/interface/ParameterSet.h"
8 #include "FWCore/MessageLogger/interface/MessageLogger.h"
9 
10 #include "SeedForPhotonConversion1Leg.h"
11 #include "RecoTracker/TkSeedGenerator/interface/FastHelix.h"
12 
13 #include "IdealHelixParameters.h"
14 
15 #include "PrintRecoObjects.h"
16 
17 #include "DataFormats/VertexReco/interface/Vertex.h"
18 #include "DataFormats/VertexReco/interface/VertexFwd.h"
19 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
20 #include "DataFormats/TrackReco/interface/TrackFwd.h"
21 
22 #include "CombinedHitPairGeneratorForPhotonConversion.h"
23 
24 #include "RecoTracker/TkSeedGenerator/interface/ClusterChecker.h"
25 #include "RecoTracker/TkTrackingRegions/plugins/GlobalTrackingRegionProducerFromBeamSpot.h"
26 #include "FWCore/Framework/interface/ConsumesCollector.h"
27 
28 #include <sstream>
29 
30 inline bool lt_(std::pair<double,short> a, std::pair<double,short> b) { return a.first<b.first; }
31 
32 class dso_hidden PhotonConversionTrajectorySeedProducerFromSingleLegAlgo{
33 
34 
35  public:
36 
37  PhotonConversionTrajectorySeedProducerFromSingleLegAlgo(const edm::ParameterSet &,
38  edm::ConsumesCollector && iC);
39  ~PhotonConversionTrajectorySeedProducerFromSingleLegAlgo();
40 
41  void find(const edm::Event & event, const edm::EventSetup & setup, TrajectorySeedCollection & output);
42 
43  IdealHelixParameters* getIdealHelixParameters(){return &_IdealHelixParameters;}
44 
45  private:
46 
47  void loopOnTracks();
48  bool inspectTrack(const reco::Track* track, const TrackingRegion & region, math::XYZPoint& primaryVertexPoint);
49 
50  bool rejectTrack(const reco::Track& track);
51 
52  bool selectPriVtxCompatibleWithTrack(const reco::Track& tk, std::vector<reco::Vertex>& selectedPriVtxCompatibleWithTrack);
53  void loopOnPriVtx(const reco::Track& tk, const std::vector<reco::Vertex>& selectedPriVtxCompatibleWithTrack);
54 
55  //Data Members
56 
57  TrajectorySeedCollection * seedCollection=nullptr;
58 
59  std::unique_ptr<CombinedHitPairGeneratorForPhotonConversion> theHitsGenerator;
60  std::unique_ptr<SeedForPhotonConversion1Leg> theSeedCreator;
61  std::unique_ptr<GlobalTrackingRegionProducerFromBeamSpot> theRegionProducer;
62 
63 
64  ClusterChecker theClusterCheck;
65  bool theSilentOnClusterCheck;
66 
67  double _vtxMinDoF, _maxDZSigmas;
68  size_t _maxNumSelVtx;
69  bool _applyTkVtxConstraint;
70  size_t _countSeedTracks;
71  edm::InputTag _primaryVtxInputTag, _beamSpotInputTag;
72  edm::EDGetTokenT<reco::VertexCollection> token_vertex;
73  edm::EDGetTokenT<reco::BeamSpot> token_bs;
74  edm::EDGetTokenT<reco::TrackCollection> token_refitter;
75 
76  typedef std::vector<TrackingRegion* > Regions;
77  typedef Regions::const_iterator IR;
78  Regions regions;
79 
80  edm::Handle<reco::TrackCollection> trackCollectionH;
81 
82  const edm::EventSetup* myEsetup;
83  const edm::Event* myEvent;
84 
85  const MagneticField* magField;
86  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
87  const reco::BeamSpot * theBeamSpot;
88 
89  IdealHelixParameters _IdealHelixParameters;
90 
91  edm::Handle<reco::VertexCollection> vertexHandle;
92  reco::Vertex primaryVertex;
93 
94  PrintRecoObjects po;
95 
96  std::stringstream ss;
97 
98 };
99 #endif
dso_hidden
#define dso_hidden
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
TrajectorySeedCollection
std::vector< TrajectorySeed > TrajectorySeedCollection
Definition: TrajectorySeedCollection.h:6
event
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
convertSQLitetoXML_cfg.output
tuple output
Definition: convertSQLitetoXML_cfg.py:32
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
b
double b
Definition: hdecay.h:120
PhotonConversionTrajectorySeedProducerFromSingleLegAlgo::theHitsGenerator
std::unique_ptr< CombinedHitPairGeneratorForPhotonConversion > theHitsGenerator
Definition: PhotonConversionTrajectorySeedProducerFromSingleLegAlgo.h:59
a
double a
Definition: hdecay.h:121
lt_
bool lt_(std::pair< double, short > a, std::pair< double, short > b)
Definition: PhotonConversionTrajectorySeedProducerFromSingleLegAlgo.h:30
HcalObjRepresent::setup
void setup(std::vector< TH2F > &depth, std::string name, std::string units="")
Definition: HcalObjRepresent.h:628
Generated for CMSSW Reference Manual by doxygen 1.8.5