32 _analysisname(analysisname),
34 _hltProcessName(pset.getParameter<std::
string>(
"hltProcessName")),
35 _genParticleLabel(pset.getParameter<std::
string>(
"genParticleLabel")),
36 _trigResultsLabel(
"TriggerResults",
"", _hltProcessName),
37 _beamSpotLabel(pset.getParameter<std::
string>(
"beamSpotLabel")),
38 _parametersEta(pset.getParameter<std::vector<double> >(
"parametersEta")),
39 _parametersPhi(pset.getParameter<std::vector<double> >(
"parametersPhi")),
40 _parametersTurnOn(pset.getParameter<std::vector<double> >(
"parametersTurnOn")),
41 _parametersDxy(pset.getParameter<std::vector<double> >(
"parametersDxy")),
43 _recMuonTrkSelector(0),
50 _recCaloMETSelector(0),
53 _recPhotonSelector(0),
55 _recCaloJetSelector(0)
58 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::constructor()";
67 if (anpset.
exists(
"parametersTurnOn")) {
71 if (anpset.
exists(
"parametersEta")) {
75 if (anpset.
exists(
"parametersPhi")) {
79 if (anpset.
exists(
"parametersDxy")) {
93 for (std::map<unsigned int, edm::InputTag>::const_iterator it =
_recLabels.begin();
98 auto const genCutParam = objStr +
"_genCut_leading";
99 if (pset.
exists(genCutParam)) {
104 auto const recCutParam = objStr +
"_recCut_leading";
105 if (pset.
exists(recCutParam)) {
113 for (std::map<unsigned int, edm::InputTag>::const_iterator it =
_recLabels.begin();
117 auto const genCutParam = objStr +
"_genCut";
122 auto const recCutParam = objStr +
"_recCut";
190 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::subAnalysisBookHistos()";
197 for (std::map<unsigned int, edm::InputTag>::const_iterator it =
_recLabels.begin();
200 std::vector<std::string> sources(2);
204 for (
size_t i = 0;
i < sources.size();
i++) {
207 if ( source ==
"gen" ) {
208 if ( TString(objStr).Contains(
"MET") ||
209 TString(objStr).Contains(
"MHT") ||
210 TString(objStr).Contains(
"Jet") ) {
213 bookHist(iBooker, source, objStr,
"MaxPt1");
214 bookHist(iBooker, source, objStr,
"MaxPt2");
215 bookHist(iBooker, source, objStr,
"MaxPt3");
216 bookHist(iBooker, source, objStr,
"Eta");
217 bookHist(iBooker, source, objStr,
"Phi");
223 bookHist(iBooker, source, objStr,
"Dxy");
227 if ( TString(objStr).Contains(
"MET") ||
228 TString(objStr).Contains(
"MHT") ) {
229 bookHist(iBooker, source, objStr,
"MaxPt1");
230 bookHist(iBooker, source, objStr,
"SumEt");
232 bookHist(iBooker, source, objStr,
"MaxPt1");
233 bookHist(iBooker, source, objStr,
"MaxPt2");
234 bookHist(iBooker, source, objStr,
"MaxPt3");
235 bookHist(iBooker, source, objStr,
"Eta");
236 bookHist(iBooker, source, objStr,
"Phi");
242 bookHist(iBooker, source, objStr,
"Dxy");
252 for (std::vector<HLTExoticaPlotter>::iterator it =
_plotters.begin();
254 it->plotterBookHistos(iBooker, iRun, iSetup);
261 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::beginRun()";
267 bool changedConfig(
true);
269 edm::LogError(
"ExoticaValidation") <<
"HLTExoticaSubAnalysis::constructor(): "
270 <<
"Initialization of HLTConfigProvider failed!";
283 if (TString(thetriggername).Contains(pattern)) {
288 std::cout <<
"--- TRIGGER PATH : " << thetriggername << std::endl;
293 edm::LogWarning(
"ExoticaValidation") <<
"HLTExoticaSubAnalysis::constructor(): In "
302 <<
"\nHLT Trigger Paths found >>>";
303 for (std::set<std::string>::const_iterator
iter =
_hltPaths.begin();
305 LogTrace(
"ExoticaValidation") << (*iter) <<
"\n";
310 for (std::set<std::string>::iterator iPath =
_hltPaths.begin();
315 if (path.rfind(
"_v") < path.length()) {
316 shortpath = path.substr(0, path.rfind(
"_v"));
324 std::vector<unsigned int> objsNeedHLT;
325 for (std::map<unsigned int, edm::InputTag>::iterator it =
_recLabels.begin() ;
327 objsNeedHLT.push_back(it->first);
352 LogTrace(
"ExoticaValidation") <<
" --- " << shortpath;
365 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::analyze()";
387 std::vector<reco::LeafCandidate> matchesGen; matchesGen.clear();
388 std::vector<reco::LeafCandidate> matchesReco; matchesReco.clear();
389 std::map<int , double> theSumEt;
390 std::map<int, std::vector<const reco::Track*> > trkObjs;
397 for (std::map<unsigned int, edm::InputTag>::iterator it =
_recLabels.begin();
412 if ( TString(objTypeStr).Contains(
"MET") ||
413 TString(objTypeStr).Contains(
"MHT") ||
414 TString(objTypeStr).Contains(
"Jet") )
continue;
428 matchesGen.push_back(
v);
444 for (std::map<unsigned int, edm::InputTag>::iterator it =
_recLabels.begin();
454 <<
" sumEt=" << theSumEt[it->first] << std::endl;
476 if(
verbose>2)
std::cerr <<
"### matchesGen.size() = " << matchesGen.size() << std::endl;
489 std::map<unsigned int, int> countobjects;
490 for (std::map<unsigned int, edm::InputTag>::iterator co =
_recLabels.begin();
493 countobjects.insert(std::pair<unsigned int, int>(co->first, 0));
499 int totalobjectssize3 = 3 * countobjects.size();
502 bool isPassedLeadingCut =
true;
504 for (
size_t j = 0;
j != matchesGen.size(); ++
j) {
505 const unsigned int objType = matchesGen[
j].pdgId();
508 if ( !
select( matchesGen[
j] ) ) {
509 isPassedLeadingCut =
false;
515 std::vector<float> dxys; dxys.clear();
517 for (
size_t j = 0; (
j != matchesGen.size() ) && isPassedLeadingCut; ++
j) {
518 const unsigned int objType = matchesGen[
j].pdgId();
522 float pt = matchesGen[
j].pt();
524 if (countobjects[objType] == 0) {
525 this->
fillHist(
"gen", objTypeStr,
"MaxPt1", pt);
526 ++(countobjects[objType]);
529 else if (countobjects[objType] == 1) {
530 this->
fillHist(
"gen", objTypeStr,
"MaxPt2", pt);
531 ++(countobjects[objType]);
534 else if (countobjects[objType] == 2) {
535 this->
fillHist(
"gen", objTypeStr,
"MaxPt3", pt);
536 ++(countobjects[objType]);
541 if (counttotal == totalobjectssize3) {
542 size_t max_size = matchesGen.size();
543 for (
size_t jj =
j;
jj < max_size;
jj++ ) {
544 matchesGen.erase(matchesGen.end());
550 float eta = matchesGen[
j].eta();
551 float phi = matchesGen[
j].phi();
553 this->
fillHist(
"gen", objTypeStr,
"Eta", eta);
554 this->
fillHist(
"gen", objTypeStr,
"Phi", phi);
561 float momphi = matchesGen[
j].momentum().phi();
562 float dxyGen = (-(vtx.x()-cols->
bs->
x0())*
sin(momphi)+(vtx.y()-cols->
bs->
y0())*
cos(momphi));
563 dxys.push_back(dxyGen);
564 this->
fillHist(
"gen", objTypeStr,
"Dxy", dxyGen);
571 for (std::vector<HLTExoticaPlotter>::iterator an =
_plotters.begin(); an !=
_plotters.end(); ++an) {
574 LogDebug(
"ExoticaValidation") <<
" preparing to call the plotters analysis";
575 an->analyze(ispassTrigger,
"gen", matchesGen, theSumEt, dxys);
576 LogDebug(
"ExoticaValidation") <<
" called the plotter";
583 if(
verbose>2)
std::cerr <<
"### matchesReco.size() = " << matchesReco.size() << std::endl;
599 std::map<unsigned int, int> countobjects;
600 for (std::map<unsigned int, edm::InputTag>::iterator co =
_recLabels.begin();
602 countobjects.insert(std::pair<unsigned int, int>(co->first, 0));
608 int totalobjectssize3 = 3 * countobjects.size();
613 std::vector<float> dxys; dxys.clear();
615 bool isPassedLeadingCut =
true;
617 for (
size_t j = 0;
j != matchesReco.size(); ++
j) {
618 const unsigned int objType = matchesReco[
j].pdgId();
621 if ( !
select( matchesReco[
j] ) ) {
622 isPassedLeadingCut =
false;
628 for (
size_t j = 0; (
j != matchesReco.size() ) && isPassedLeadingCut; ++
j) {
629 const unsigned int objType = matchesReco[
j].pdgId();
633 float pt = matchesReco[
j].pt();
635 if (countobjects[objType] == 0) {
636 this->
fillHist(
"rec", objTypeStr,
"MaxPt1", pt);
637 ++(countobjects[objType]);
640 else if (countobjects[objType] == 1) {
641 if( ! ( TString(objTypeStr).Contains(
"MET") || TString(objTypeStr).Contains(
"MHT") ) ) {
642 this->
fillHist(
"rec", objTypeStr,
"MaxPt2", pt);
644 ++(countobjects[objType]);
647 else if (countobjects[objType] == 2) {
648 if( ! ( TString(objTypeStr).Contains(
"MET") || TString(objTypeStr).Contains(
"MHT") ) ) {
649 this->
fillHist(
"rec", objTypeStr,
"MaxPt3", pt);
651 ++(countobjects[objType]);
656 if (counttotal == totalobjectssize3) {
657 size_t max_size = matchesReco.size();
658 for (
size_t jj =
j;
jj < max_size;
jj++ ) {
659 matchesReco.erase(matchesReco.end());
665 float eta = matchesReco[
j].eta();
666 float phi = matchesReco[
j].phi();
668 if ( !( TString(objTypeStr).Contains(
"MET") || TString(objTypeStr).Contains(
"MHT") ) ) {
669 this->
fillHist(
"rec", objTypeStr,
"Eta", eta);
670 this->
fillHist(
"rec", objTypeStr,
"Phi", phi);
673 this->
fillHist(
"rec", objTypeStr,
"SumEt", theSumEt[objType]);
676 if (trkObjs[objType].
size()>=
j+1) {
677 float dxyRec = trkObjs[objType].at(
j)->dxy(cols->
bs->
position());
678 this->
fillHist(
"rec", objTypeStr,
"Dxy", dxyRec);
679 dxys.push_back(dxyRec);
689 for (std::vector<HLTExoticaPlotter>::iterator an =
_plotters.begin(); an !=
_plotters.end(); ++an) {
692 LogDebug(
"ExoticaValidation") <<
" preparing to call the plotters analysis";
693 an->analyze(ispassTrigger,
"rec", matchesReco, theSumEt, dxys);
694 LogDebug(
"ExoticaValidation") <<
" called the plotter";
704 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::getObjectsType()";
706 static const unsigned int objSize = 14;
707 static const unsigned int objtriggernames[] = {
724 std::set<unsigned int> objsType;
726 for (
unsigned int i = 0;
i < objSize; ++
i) {
730 if (! TString(hltPath).Contains(objTypeStr)) {
734 objsType.insert(objtriggernames[i]);
737 return std::vector<unsigned int>(objsType.begin(), objsType.end());
743 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::getNamesOfObjects()";
745 if (anpset.
exists(
"recMuonLabel")) {
749 if (anpset.
exists(
"recMuonTrkLabel")) {
753 if (anpset.
exists(
"recTrackLabel")) {
757 if (anpset.
exists(
"recElecLabel")) {
761 if (anpset.
exists(
"recPhotonLabel")) {
765 if (anpset.
exists(
"recMETLabel")) {
769 if (anpset.
exists(
"recPFMETLabel")) {
773 if (anpset.
exists(
"recPFMHTLabel")) {
777 if (anpset.
exists(
"genMETLabel")) {
781 if (anpset.
exists(
"recCaloMETLabel")) {
785 if (anpset.
exists(
"hltMETLabel")) {
789 if (anpset.
exists(
"l1METLabel")) {
793 if (anpset.
exists(
"recPFTauLabel")) {
797 if (anpset.
exists(
"recPFJetLabel")) {
801 if (anpset.
exists(
"recCaloJetLabel")) {
807 edm::LogError(
"ExoticaValidation") <<
"HLTExoticaSubAnalysis::getNamesOfObjects, "
808 <<
"Not included any object (recMuonLabel, recElecLabel, ...) "
832 for (std::map<unsigned int, edm::InputTag>::iterator it =
_recLabels.begin();
905 edm::LogError(
"ExoticaValidation") <<
"HLTExoticaSubAnalysis::registerConsumes"
906 <<
" NOT IMPLEMENTED (yet) ERROR: '" << it->second.label() <<
"'";
915 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::getHandlesToObjects()";
923 LogDebug(
"ExoticaValidation") <<
"Added handle to triggerResults";
931 LogDebug(
"ExoticaValidation") <<
"Added handle to genParticles";
943 LogDebug(
"ExoticaValidation") <<
"We have got " <<
_tokens.size() <<
"tokens";
944 for (std::map<unsigned int, edm::EDGetToken>::iterator it =
_tokens.begin();
1017 edm::LogError(
"ExoticaValidation") <<
"HLTExoticaSubAnalysis::getHandlesToObjects "
1018 <<
" NOT IMPLEMENTED (yet) ERROR: '" << it->first <<
"'";
1029 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::bookHist()";
1031 sourceUpper[0] = std::toupper(sourceUpper[0]);
1035 if (variable.find(
"SumEt") != std::string::npos) {
1038 float *
edges =
new float[nBins + 1];
1039 for (
size_t i = 0;
i < nBins + 1;
i++) {
1042 h =
new TH1F(name.c_str(), title.c_str(), nBins,
edges);
1046 else if (variable.find(
"Dxy") != std::string::npos) {
1051 h =
new TH1F(name.c_str(), title.c_str(), nBins,
min,
max);
1054 else if (variable.find(
"MaxPt") != std::string::npos) {
1055 std::string desc = (variable ==
"MaxPt1") ?
"Leading" :
"Next-to-Leading";
1058 float *
edges =
new float[nBins + 1];
1059 for (
size_t i = 0;
i < nBins + 1;
i++) {
1062 h =
new TH1F(name.c_str(), title.c_str(), nBins,
edges);
1067 std::string symbol = (variable ==
"Eta") ?
"#eta" :
"#phi";
1070 int nBins = (int)params[0];
1071 double min = params[1];
1072 double max = params[2];
1073 h =
new TH1F(name.c_str(), title.c_str(), nBins,
min,
max);
1088 const float &
value)
1091 sourceUpper[0] = toupper(sourceUpper[0]);
1094 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::fillHist() " << name <<
" " <<
value;
1096 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::fillHist() " << name <<
" worked";
1103 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::initSelector()";
1142 std::map<int,double> & theSumEt,
std::map<
int,std::vector<const reco::Track*> > & trkObjs )
1145 LogDebug(
"ExoticaValidation") <<
"In HLTExoticaSubAnalysis::insertCandidates()";
1147 theSumEt[objType] = -1;
1150 for (
size_t i = 0;
i < cols->
muons->size();
i++) {
1151 LogDebug(
"ExoticaValidation") <<
"Inserting muon " <<
i ;
1154 matches->push_back(
m);
1157 trkObjs[objType].push_back(cols->
muons->at(
i).bestTrack());
1161 for (
size_t i = 0;
i < cols->
tracks->size();
i++) {
1162 LogDebug(
"ExoticaValidation") <<
"Inserting muonTrack " <<
i ;
1164 ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double>> mom4;
1166 mom4.SetXYZT(mom3.x(),mom3.y(),mom3.z(),mom3.r());
1168 matches->push_back(
m);
1172 for (
size_t i = 0;
i < cols->
tracks->size();
i++) {
1173 LogDebug(
"ExoticaValidation") <<
"Inserting Track " <<
i ;
1175 ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double>> mom4;
1177 mom4.SetXYZT(mom3.x(),mom3.y(),mom3.z(),mom3.r());
1179 matches->push_back(
m);
1184 LogDebug(
"ExoticaValidation") <<
"Inserting electron " <<
i ;
1187 matches->push_back(
m);
1190 trkObjs[objType].push_back(cols->
electrons->at(
i).bestTrack());
1194 for (
size_t i = 0;
i < cols->
photons->size();
i++) {
1195 LogDebug(
"ExoticaValidation") <<
"Inserting photon " <<
i ;
1198 matches->push_back(
m);
1205 for (
size_t i = 0;
i < cols->
pfMETs->size();
i++) {
1206 LogDebug(
"ExoticaValidation") <<
"Inserting PFMET " <<
i ;
1209 matches->push_back(
m);
1210 if(
i==0) theSumEt[objType] = cols->
pfMETs->at(
i).sumEt();
1214 for (
size_t i = 0;
i < cols->
pfMHTs->size();
i++) {
1215 LogDebug(
"ExoticaValidation") <<
"Inserting PFMHT " <<
i ;
1218 matches->push_back(
m);
1219 if(
i==0) theSumEt[objType] = cols->
pfMHTs->at(
i).sumEt();
1223 for (
size_t i = 0;
i < cols->
genMETs->size();
i++) {
1224 LogDebug(
"ExoticaValidation") <<
"Inserting GENMET " <<
i ;
1227 matches->push_back(
m);
1228 if(
i==0) theSumEt[objType] = cols->
genMETs->at(
i).sumEt();
1232 for (
size_t i = 0;
i < cols->
caloMETs->size();
i++) {
1233 LogDebug(
"ExoticaValidation") <<
"Inserting CALOMET " <<
i ;
1236 matches->push_back(
m);
1237 if(
i==0) theSumEt[objType] = cols->
caloMETs->at(
i).sumEt();
1241 for (
size_t i = 0;
i < cols->
l1METs->size();
i++) {
1242 LogDebug(
"ExoticaValidation") <<
"Inserting L1MET " <<
i ;
1245 matches->push_back(
m);
1246 if(
i==0) theSumEt[objType] = cols->
l1METs->at(
i).etTotal();
1250 for (
size_t i = 0;
i < cols->
pfTaus->size();
i++) {
1251 LogDebug(
"ExoticaValidation") <<
"Inserting PFtau " <<
i ;
1254 matches->push_back(
m);
1258 for (
size_t i = 0;
i < cols->
pfJets->size();
i++) {
1259 LogDebug(
"ExoticaValidation") <<
"Inserting jet " <<
i ;
1262 matches->push_back(
m);
1266 for (
size_t i = 0;
i < cols->
caloJets->size();
i++) {
1267 LogDebug(
"ExoticaValidation") <<
"Inserting jet " <<
i ;
1270 matches->push_back(
m);
const std::vector< reco::Muon > * muons
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
StringCutObjectSelector< reco::Track > * _recMuonTrkSelector
std::vector< GenParticle > GenParticleCollection
collection of GenParticles
T getParameter(std::string const &) const
unsigned int _minCandidates
The minimum number of reco/gen candidates needed by the analysis.
Entry const & retrieve(char const *) const
T getUntrackedParameter(std::string const &, T const &) const
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
void getNamesOfObjects(const edm::ParameterSet &anpset)
Creates the maps that map which collection should come from which label.
std::vector< PFTau > PFTauCollection
collection of PFTau objects
bool existsAs(std::string const ¶meterName, bool trackiness=true) const
checks if a parameter exists as a given type
const std::vector< reco::CaloMET > * caloMETs
bool getByToken(EDGetToken token, Handle< PROD > &result) const
std::vector< reco::GenMET > GenMETCollection
collection of GenMET objects
std::map< unsigned int, std::string > _genCut
gen/rec objects cuts
Sin< T >::type sin(const T &t)
bool accept() const
Has at least one path accepted the event?
const std::vector< std::string > & triggerNames() const
names of trigger paths
std::vector< Track > TrackCollection
collection of Tracks
std::vector< HLTExoticaPlotter > _plotters
The plotters: managers of each hlt path where the plots are done.
bool exists(std::string const ¶meterName) const
checks if a parameter exists
StringCutObjectSelector< reco::CaloJet > * _recCaloJetSelector
StringCutObjectSelector< reco::GenMET > * _genMETSelector
std::vector< double > _parametersPhi
const std::vector< reco::CaloJet > * caloJets
void insert(bool ok_to_replace, char const *, Entry const &)
void fillHist(const std::string &source, const std::string &objType, const std::string &variable, const float &value)
const reco::GenParticleCollection * genParticles
container with all the objects needed
void analyze(const edm::Event &iEvent, const edm::EventSetup &iEventSetup, EVTColContainer *cols)
edm::InputTag _trigResultsLabel
StringCutObjectSelector< reco::MET > * _recMETSelector
edm::EDGetTokenT< edm::TriggerResults > _trigResultsToken
std::vector< reco::MET > METCollection
collection of MET objects
std::vector< GsfElectron > GsfElectronCollection
collection of GsfElectron objects
const std::vector< reco::PFTau > * pfTaus
std::vector< Muon > MuonCollection
collection of Muon objects
std::map< unsigned int, std::string > _genCut_leading
gen/rec pt-leading objects cuts
StringCutObjectSelector< reco::PFJet > * _recPFJetSelector
void setPFMHT(const reco::PFMETCollection *v)
std::map< unsigned int, edm::InputTag > _recLabels
StringCutObjectSelector< reco::GsfElectron > * _recElecSelector
void set(const reco::MuonCollection *v)
Setter: multiple overloaded function.
std::map< unsigned int, std::string > _recCut_leading
const std::vector< reco::PFMET > * pfMHTs
StringCutObjectSelector< reco::Muon > * _recMuonSelector
tuple path
else: Piece not in the list, fine.
unsigned int triggerIndex(std::string const &name) const
const std::vector< reco::GenMET > * genMETs
edm::ParameterSet _pset
Internal, working copy of the PSet passed from above.
const reco::BeamSpot * bs
void subAnalysisBookHistos(DQMStore::IBooker &iBooker, const edm::Run &iRun, const edm::EventSetup &iSetup)
std::vector< double > _parametersEta
Some kinematical parameters.
std::map< unsigned int, std::string > _recCut
std::vector< double > _parametersDxy
Cos< T >::type cos(const T &t)
const std::vector< unsigned int > getObjectsType(const std::string &hltpath) const
closes analyze method
StringCutObjectSelector< reco::Track > * _recTrackSelector
MonitorElement * book1D(Args &&...args)
const std::vector< reco::PFJet > * pfJets
HLTConfigProvider _hltConfig
Interface to the HLT information.
const std::vector< reco::Track > * tracks
const std::vector< reco::Photon > * photons
virtual const Point & vertex() const =0
vertex position
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
~HLTExoticaSubAnalysis()
End Constructor.
std::vector< std::string > _hltPathsToCheck
The hlt paths to check for.
std::string _hltProcessName
The labels of the object collections to be used in this analysis.
static const char *const trigNames[]
Transform3DPJ::Vector XYZVector
const std::vector< l1extra::L1EtMissParticle > * l1METs
void setCurrentFolder(const std::string &fullpath)
std::vector< reco::CaloMET > CaloMETCollection
collection of CaloMET objects
void registerConsumes(edm::ConsumesCollector &consCollector)
Registers consumption of objects.
T const * product() const
void getHandlesToObjects(const edm::Event &iEvent, EVTColContainer *col)
Gets the collections themselves.
XYZPointD XYZPoint
point in space with cartesian internal representation
StringCutObjectSelector< reco::PFMET > * _recPFMETSelector
std::string _analysisname
The name of this sub-analysis.
StringCutObjectSelector< reco::Photon > * _recPhotonSelector
std::vector< Photon > PhotonCollection
collectin of Photon objects
std::map< unsigned int, StringCutObjectSelector< reco::GenParticle > * > _genSelectorMap
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d'tor
void beginRun(const edm::Run &iRun, const edm::EventSetup &iEventSetup)
StringCutObjectSelector< l1extra::L1EtMissParticle > * _l1METSelector
const std::vector< reco::GsfElectron > * electrons
std::vector< double > _parametersTurnOn
std::map< unsigned int, edm::EDGetToken > _tokens
std::map< std::string, MonitorElement * > _elements
Structure of the MonitorElements.
StringCutObjectSelector< reco::CaloMET > * _recCaloMETSelector
std::vector< PFJet > PFJetCollection
collection of PFJet objects
HLTExoticaSubAnalysis(const edm::ParameterSet &pset, const std::string &analysisname, edm::ConsumesCollector &&consCollector)
Constructor.
std::vector< reco::PFMET > PFMETCollection
collection of PFMET objects
double y0() const
y coordinate
edm::EDGetTokenT< reco::GenParticleCollection > _genParticleToken
And also the tokens to get the object collections.
const Point & position() const
position
void initSelector(const unsigned int &objtype)
Initializes the selectors of the objects based on which object it is.
edm::InputTag _genParticleLabel
const std::vector< reco::PFMET > * pfMETs
std::set< std::string > _hltPaths
The hlt paths found in the hltConfig.
static const std::string getTypeString(const unsigned int &objtype)
Tranform types into strings.
edm::InputTag _beamSpotLabel
const edm::TriggerResults * triggerResults
StringCutObjectSelector< reco::PFMET > * _recPFMHTSelector
static std::string const source
edm::EDGetTokenT< reco::BeamSpot > _bsToken
StringCutObjectSelector< reco::PFTau > * _recPFTauSelector
tuple size
Write out results.
void bookHist(DQMStore::IBooker &iBooker, const std::string &source, const std::string &objType, const std::string &variable)
The internal functions to book and fill histograms.
std::vector< CaloJet > CaloJetCollection
collection of CaloJet objects
std::map< std::string, std::string > _shortpath2long
Relation between the short and long versions of the path.
void insertCandidates(const unsigned int &objtype, const EVTColContainer *col, std::vector< reco::LeafCandidate > *matches, std::map< int, double > &theSumEt, std::map< int, std::vector< const reco::Track * > > &trkObjs)
double x0() const
x coordinate
virtual const LorentzVector & p4() const =0
four-momentum Lorentz vector