readTestVector.cc File Reference

#include <iostream>
#include <iomanip>
#include <fstream>
#include <memory>
#include <map>
#include "DataFormats/L1Trigger/interface/EGamma.h"
#include "DataFormats/L1Trigger/interface/Muon.h"
#include "DataFormats/L1Trigger/interface/Tau.h"
#include "DataFormats/L1Trigger/interface/Jet.h"
#include "DataFormats/L1Trigger/interface/EtSum.h"
#include "L1Trigger/L1TGlobal/src/L1TMenuEditor/L1TriggerMenu.hxx"
#include "FWCore/Utilities/interface/typedefs.h"
#include <boost/program_options.hpp>
#include "TH1.h"
#include "TFile.h"

Go to the source code of this file.

Macros
#define	M_PI   3.14159265358979323846

Functions
double	convertEtaFromHW (int hwEta, double max, double step, int hwMax)
double	convertPhiFromHW (int hwPhi, double step)
double	convertPtFromHW (int hwPt, double max, double step)
int	main (int argc, char **argv)
void	parseAlgo (std::string algo)
void	parseEGs (std::vector< std::string > egs, bool verbose)
void	parseEtSums (std::vector< std::string > etsums, bool verbose)
void	parseJets (std::vector< std::string > jets, bool verbose)
void	parseMuons (std::vector< std::string > muons, bool verbose)
void	parseTaus (std::vector< std::string > taus, bool verbose)
l1t::EGamma	unpackEGs (std::string ieg)
l1t::EtSum	unpackEtSums (std::string ietsum, l1t::EtSum::EtSumType type)
l1t::Jet	unpackJets (std::string ijet)
l1t::Muon	unpackMuons (std::string imu)
l1t::Tau	unpackTaus (std::string itau)

Variables
double	EtaStepCalo_ = 230
double	EtaStepMuon_ = 450
TH1D *	h_l1eg_eta_
TH1D *	h_l1eg_isolation_
TH1D *	h_l1eg_num_
TH1D *	h_l1eg_phi_
TH1D *	h_l1eg_pt_
TH1D *	h_l1et_
TH1D *	h_l1etm_et_
TH1D *	h_l1etm_phi_
TH1D *	h_l1ht_
TH1D *	h_l1htm_et_
TH1D *	h_l1htm_phi_
TH1D *	h_l1jet_eta_
TH1D *	h_l1jet_num_
TH1D *	h_l1jet_phi_
TH1D *	h_l1jet_pt_
TH1D *	h_l1mu_charge_
TH1D *	h_l1mu_eta_
TH1D *	h_l1mu_isolation_
TH1D *	h_l1mu_num_
TH1D *	h_l1mu_phi_
TH1D *	h_l1mu_pt_
TH1D *	h_l1mu_quality_
TH1D *	h_l1tau_eta_
TH1D *	h_l1tau_isolation_
TH1D *	h_l1tau_num_
TH1D *	h_l1tau_phi_
TH1D *	h_l1tau_pt_
std::vector< int >	l1taccepts
std::vector< int >	l1taccepts_evt
std::vector< std::string >	l1tnames
const int	MAX_ALGO_BITS = 512
double	MaxCaloEta_ = 5.0
double	MaxEt_ = 2047
double	MaxJetPt_ = 1023
double	MaxLepPt_ = 255
double	MaxMuonEta_ = 2.45
double	PhiStepCalo_ = 144
double	PhiStepMuon_ = 576
double	PtStep_ = 0.5
std::string	zeroes128 = std::string(128,'0')
std::string	zeroes16 = std::string(16,'0')
std::string	zeroes8 = std::string(8,'0')

Macro Definition Documentation

#define M_PI   3.14159265358979323846

Definition at line 25 of file readTestVector.cc.

Referenced by GflashTrajectory::_refreshCache(), reco::WMuNuCandidate::acop(), reco::WMuNuCandidatePtr::acop(), gen::Pythia6Gun::addAntiParticle(), fireworks::addCircle(), fireworks::addDashedLine(), fireworks::addDoubleLines(), ALIUtils::addPii(), OpticalObject::addPii(), PixelUnpackingRegions::addRegionLocal(), MuonCSCChamberResidual::addResidual(), TrackDetectorAssociator::addTAMuonSegmentMatch(), heppy::FSRWeightAlgo::alphaRatio(), FSRWeightProducer::alphaRatio(), L1TStage2EMTF::analyze(), QcdHighPtDQM::analyze(), EwkMuDQM::analyze(), EwkMuLumiMonitorDQM::analyze(), L1GctValidation::analyze(), CheckHitPattern::analyze(), TrackParameterAnalyzer::analyze(), DTTrigTest::analyze(), MuonGeometrySVGTemplate::analyze(), EcalSimHitsValidation::analyze(), EcalRecHitsValidation::analyze(), MuonDTDigis::analyze(), TestOutliers::analyze(), EcalDigisValidation::analyze(), FFTJetPileupAnalyzer::analyze(), L1TCSCTF::analyze(), MuonTrackValidator::analyze(), B2GDQM::analyzeAllHad(), B2GDQM::analyzeSemiE(), B2GDQM::analyzeSemiMu(), CSCSegAlgoSK::areHitsCloseInGlobalPhi(), CSCSegAlgoTC::areHitsCloseInGlobalPhi(), FWCaloDataHistProxyBuilder::assertCaloDataSlice(), myFastSimVal::beginJob(), JetValidation::beginJob(), JetPlotsExample< Jet >::beginJob(), JetAnaPythia< Jet >::beginJob(), PFTester::beginJob(), DijetMass< Jet >::beginJob(), PatTrackAnalyzer::beginJob(), WMuNuSelector::beginJob(), L1GctValidation::beginJob(), myJetAna::beginJob(), TrackParameterAnalyzer::beginJob(), MuonCSCSeedFromRecHits::bestEndcapHit(), AlignmentMonitorGeneric::book(), AlignmentMonitorMuonSystemMap1D::book(), QcdHighPtDQM::bookHistograms(), EwkMuDQM::bookHistograms(), L1TCSCTF::bookHistograms(), METAnalyzer::bookMonitorElement(), HiEvtPlaneFlatten::bounds(), HiEvtPlaneFlatten::bounds2(), BremsstrahlungSimulator::brem(), MuonBremsstrahlungSimulator::brem(), FWTauProxyBuilderBase::buildBaseTau(), FWECALDetailViewBuilder::buildCaloData(), FFTJetEFlowSmoother::buildKernelConvolver(), FFTJetPileupProcessor::buildKernelConvolver(), FFTJetPatRecoProducer::buildKernelConvolver(), KDTreeLinkerTrackEcal::buildTree(), KDTreeLinkerTrackHcal::buildTree(), FWMETProxyBuilder::buildViewType(), FWJetProxyBuilder::buildViewType(), L1TMuon::GeometryTranslator::calcDTSpecificPoint(), CSCSectorReceiverLUT::calcGlobalPhiME(), GsfVertexWeightCalculator::calculate(), TwoTrackMinimumDistanceHelixLine::calculate(), pat::CaloIsolationEnergy::calculate(), ClusterShapeAlgo::Calculate_EnergyDepTopology(), ALIUtils::CalculateAngleDimensionFactorFromString(), DTTracoChip::calculateAngles(), JetPlusTrackProducerAA::calculateBGtracksJet(), TBPositionCalc::CalculateCMSPos(), CSCPairResidualsConstraint::calculatePhi(), tkDetUtil::calculatePhiWindow(), CSCEfficiency::chamberCandidates(), AlignmentTwoBodyDecayTrackSelector::checkAcoplanarity(), HLTEgammaL1MatchFilterPairs::CheckL1Matching(), AlignmentTwoBodyDecayTrackSelector::checkMETAcoplanarity(), ParticleKinematicLinearizedTrackState::checkParameters(), PerigeeMultiLTS::checkParameters(), PerigeeLinearizedTrackState::checkParameters(), checkPhiInRange(), PairProductionSimulator::compute(), HFShower::compute(), HDShower::compute(), EMShower::compute(), ParticleKinematicLinearizedTrackState::computeChargedJacobians(), PerigeeLinearizedTrackState::computeChargedJacobians(), GaussianSumUtilities1D::computeMode(), HFCherenkov::computeNbOfPhotons(), HFCherenkov::computeNPE(), HFCherenkov::computeNPEinPMT(), HFCherenkov::computeNPhTrapped(), DTDigitizer::computeTime(), EcalHitMaker::configureGeometry(), LaserAlignment::ConvertAngle(), convertPhiFromHW(), CordicXilinx::CordicXilinx(), JetPlusTrackCorrector::correctAA(), CaloDetIdAssociator::crossedElement(), CSCOverlapsAlignmentAlgorithm::CSCOverlapsAlignmentAlgorithm(), MuonShowerInformationFiller::cscPositionToDets(), AntiElectronIDMVA5::dCrackPhi(), AntiElectronIDMVA6::dCrackPhi(), PhiScaleHelper::deg2rad(), degrees(), delPhi(), TtDilepLRSignalSelObservables::delta(), CachedTrajectory::delta(), PFBenchmarkAlgo::deltaPhi(), HLTHemiDPhiFilter::deltaPhi(), CaloSubdetectorGeometry::deltaPhi(), L1MuonPixelTrackFitter::deltaPhi(), JetVertexMain::DeltaR(), PythiaFilterEMJetHeep::deltaR(), IsoTrackCalib::deltaR(), CalibratableTest::deltaR(), reco::deltaR2(), JetAnalyzer_HeavyIons_matching::deltaRR(), FFTJetProducer::determinePileup(), OpticalObject::diff2pi(), pat::helper::ParametrizationHelper::diffToParameters(), reco::isodeposit::Direction::Direction(), MuonSeedOrcaPatternRecognition::discriminator(), L1MuGMTPhiLUT::dphi(), IsoTrackCalib::dPhi(), IsoTrig::dPhi(), CSCPairResidualsConstraint::dphidzFromTrack(), TrackerMap::drawModule(), MuonShowerInformationFiller::dtPositionToDets(), ALIUtils::dumpDimensions(), EcalTrapezoidParameters::EcalTrapezoidParameters(), L1CaloGeometry::emJetPhiIndex(), CordicXilinx::encodeAngle(), OpticalObject::eq2ang(), LandauFluctuationGenerator::ersatzt(), L1CaloGeometry::etSumPhiIndex(), l1t::CorrCondition::evaluateCondition(), AlignmentMonitorTracksFromTrajectories::event(), KinematicPerigeeConversions::extendedPerigeeFromKinematicParameters(), HcalDDDGeometryLoader::fill(), MTVHistoProducerAlgoForTracker::fill_ResoAndPull_recoTrack_histos(), FWECALDetailViewBuilder::fillData(), MaterialBudgetEcalHistos::fillEndTrack(), HCALRecHitAnalyzer::FillGeometry(), ECALRecHitAnalyzer::FillGeometry(), ApeEstimator::fillHistsForAnalyzerMode(), ApeEstimator::fillHitHistsXForAnalyzerMode(), ApeEstimator::fillHitHistsYForAnalyzerMode(), cms::PileupJPTJetIdAlgo::fillJPTBlock(), ColinearityKinematicConstraintT< Dim >::fillValue(), PFFilter::filter(), PFMETFilter::filter(), WMuNuValidator::filter(), ThirdHitPrediction::findArcIntersection(), EgammaHLTTrackIsolation::findIsoTracks(), EgammaHLTTrackIsolation::findIsoTracksWithoutEle(), L1GctMet::floatingPointAlgo(), PerigeeConversions::ftsToPerigeeParameters(), FWLegoEvePFCandidate::FWLegoEvePFCandidate(), PairProductionSimulator::gbteth(), BremsstrahlungSimulator::gbteth(), MuonBremsstrahlungSimulator::gbteth(), gen::Pythia6JetGun::generateEvent(), gen::Pythia6PtYDistGun::generateEvent(), gen::Py8JetGun::generatePartonsAndHadronize(), LaserBeamsBarrel::GeneratePrimaries(), LaserBeamsTEC1::GeneratePrimaries(), LaserBeamsTEC2::GeneratePrimaries(), HcalGeometry::getCells(), EcalEndcapGeometry::getCells(), CaloSubdetectorGeometry::getCells(), EcalBarrelGeometry::getCells(), TrackFitter::getCharge(), HcalDDDGeometry::getClosestCell(), EcalBarrelGeometry::getClosestCell(), HGCalGeometry::getClosestCellIndex(), RPCSeedrecHitFinder::getdeltaPhifromrecHits(), DetIdAssociator::getDetIdsCloseToAPoint(), FittedEntriesManager::GetDifferentBetweenLasers(), ALIUtils::getDimensionValue(), IsoTrig::getDistInCM(), CastorSD::getEnergyDeposit(), EcalClusterToolsT< noZS >::getEnergyDepTopology(), HcalDDDRecConstants::getEtaPhi(), IsolatedPixelTrackCandidateProducer::GetEtaPhiAtEcal(), IsolatedPixelTrackCandidateL1TProducer::GetEtaPhiAtEcal(), IsoTrig::GetEtaPhiAtEcal(), NtupleManager::GetGlobalAngles(), CSCSectorReceiverLUT::getGlobalPhiValue(), GflashTrajectory::getL2DAtR(), CocoaDBMgr::GetOptAlignInfoFromOptO(), fireworks::getPhiRange(), getPhysicalPhi(), AngleConverter::getProcessorPhi(), ThirdHitPrediction::getRanges(), CachedTrajectory::getWideTrajectory(), l1t::CaloTools::gtPhi(), HcalFlexiHardcodeGeometryLoader::HcalFlexiHardcodeGeometryLoader(), HcalHardcodeGeometryLoader::HcalHardcodeGeometryLoader(), HcalTopology::HcalTopology(), HDRShower::HDRShower(), helicityAngle(), HerwigMaxPtPartonFilter::HerwigMaxPtPartonFilter(), MultiHitGeneratorFromChi2::hitSets(), PixelTripletNoTipGenerator::hitTriplets(), PixelTripletLargeTipGenerator::hitTriplets(), PixelTripletHLTGenerator::hitTriplets(), HLTPhi2METFilter::hltFilter(), HLTEgammaDoubleEtDeltaPhiFilter::hltFilter(), HLTMuonDimuonL3Filter::hltFilter(), HLTMuonTrimuonL3Filter::hltFilter(), HLTMuonDimuonL2FromL1TFilter::hltFilter(), HLTMuonDimuonL2Filter::hltFilter(), HLTExclDiJetFilter< T >::hltFilter(), HLTEgammaDoubleEtPhiFilter::hltFilter(), HLTAcoFilter::hltFilter(), HLTDoublet< T1, T2 >::hltFilter(), L1CaloGeometry::htSumPhiIndex(), WMuNuValidator::init_histograms(), GflashTrajectory::initializeTrajectory(), int2dphi(), int2phi(), reco::helper::VirtualJetProducerHelper::intersection(), DetIdAssociator::iPhi(), MuonMesh::isClusteredWith(), AngleConverter::isCSCCounterClockwise(), L1TMuon::GeometryTranslator::isCSCCounterClockwise(), MuonMesh::isDuplicateOf(), CSCSegAlgoShowering::isHitNearSegment(), CSCSegAlgoSK::isHitNearSegment(), CSCSegAlgoTC::isHitNearSegment(), CSCSegAlgoDF::isHitNearSegment(), pat::helper::ParametrizationHelper::isPhysical(), PerigeeConversions::jacobianCurvilinear2Perigee(), jacobianCurvilinearToCartesian(), PerigeeConversions::jacobianPerigee2Curvilinear(), GlobalTrackQualityProducer::kink(), L1CaloGeometry::L1CaloGeometry(), reco::GsfTrack::lambdaMode(), VoronoiAlgorithm::lp_populate(), main(), CaloTowerHardcodeGeometryLoader::makeCell(), CastorHardcodeGeometryLoader::makeCell(), HFClusterAlgo::makeCluster(), HcalFlexiHardcodeGeometryLoader::makeHECells(), Multi5x5BremRecoveryClusterAlgo::makeIslandSuperClusters(), fireworks::makeRhoPhiSuperCluster(), TrackParameterAnalyzer::match(), ThirdHitPredictionFromCircle::HelixRZ::maxCurvature(), OuterDetCompatibility::maximalLocalDisplacement(), MuonDT13ChamberResidual::MuonDT13ChamberResidual(), MuonResidualsFitter_compute_log_convolution(), MuonResidualsFitter_logGaussPowerTails(), MuonResidualsFitter_logPowerLawTails(), MuonResidualsFitter_logPureGaussian(), MuonResidualsFitter_logPureGaussian2D(), MuonTrackDT13ChamberResidual::MuonTrackDT13ChamberResidual(), AntiElectronIDMVA5::MVAValue(), AntiElectronIDMVA6::MVAValue(), SiPixelFedCablingMapBuilder::myprint(), DetIdAssociator::nearElement(), BeamProfileVtxGenerator::newVertex(), SiStripElectronSeedGenerator::normalPhi(), HcalTopology::nPhiBins(), fftjetcms::JetToPeakDistance::operator()(), fftjetcms::JetConvergenceDistance::operator()(), TtDilepLRSignalSelObservables::operator()(), OuterHitPhiPrediction::operator()(), operator<<(), phi2int(), phi2short(), Phi_0_2pi(), JetUtil::Phi_0_2pi(), Phi_mpi_pi(), JetUtil::Phi_mpi_pi(), reco::helper::CastorJetIDHelper::phiangle(), CastorClusterProducer::phiangle(), pat::PackedCandidate::phiAtVtx(), CSCSegAlgoSK::phiAtZ(), CSCSegAlgoTC::phiAtZ(), HcalTopology::phiBin(), SiStripElectronSeedGenerator::phiDiff(), hitfit::phidiff(), InnerDeltaPhi::phiRange(), TrackerMap::phival(), L1MuBMTrackSegPhi::phiValue(), L1MuDTTrackSegPhi::phiValue(), SiStripRegionCabling::position(), SiStripRegionCabling::positionIndex(), hcalCalib::Process(), l1t::Stage2Layer2DemuxEGAlgoFirmwareImp1::processEvent(), l1t::Stage2Layer2DemuxTauAlgoFirmwareImp1::processEvent(), l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent(), l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::processEvent(), AlignmentMonitorMuonSystemMap1D::processMuonResidualsFromTrack(), CastorFastTowerProducer::produce(), FFTJetEFlowSmoother::produce(), JetPlusTrackProducerAA::produce(), JetCrystalsAssociator::produce(), SelectedElectronFEDListProducer< TEle, TCand >::produce(), BaseParticlePropagator::propagate(), spr::propagateCalo(), cms::MuonTCMETValueMapProducer::propagateTrack(), spr::propagateTrackerEnd(), TCMETAlgo::propagateTrackToCalorimeterFace(), proxim(), PhiScaleHelper::rad2deg(), radiants(), radius(), JetUtil::radius(), CSCChamberFitter::radiusCorrection(), ThirdHitPredictionFromCircle::HelixRZ::rAtZ(), SimG4HcalHitJetFinder::rDist(), CSCEfficiency::recHitSegment_Efficiencies(), reco::reduceRange(), SoftLepton::refineJetAxis(), ParticleKinematicLinearizedTrackState::refittedParamFromEquation(), PerigeeMultiLTS::refittedParamFromEquation(), PerigeeLinearizedTrackState::refittedParamFromEquation(), SiStripRegionCabling::regionDimensions(), CastorSD::rotateUnitID(), RecoTracktoTP::s_lambda(), TPtoRecoTrack::s_lambda(), KDTreeLinkerTrackHcal::searchLinks(), KDTreeLinkerTrackEcal::searchLinks(), MuonIdProducer::sectorPhi(), SeedMultiplicityAnalyzer::SeedMultiplicityAnalyzer(), CastorShowerLibrary::select(), DetIdAssociator::selectAllInACone(), PFElectronAlgo::SetCandidates(), HDRShower::setHit(), PFElectronAlgo::SetIDOutputs(), PFElectronAlgo::SetLinks(), SimpleL1MuGMTCand::setPhi(), HFClusterAlgo::setup(), GenericBenchmark::setup(), short2phi(), WeakEffectsWeightProducer::sigma0_qqbarll(), HelixExtrapolatorToLine2Order::solve3rdOrder(), spin(), ThirdHitPrediction::spinCloser(), l1t::Stage1Layer2EtSumAlgorithmImpHI::Stage1Layer2EtSumAlgorithmImpHI(), l1t::Stage1Layer2EtSumAlgorithmImpHW::Stage1Layer2EtSumAlgorithmImpHW(), VoronoiAlgorithm::subtract_momentum(), LinkByRecHit::testTrackAndClusterByRecHit(), L1MuGMTLFPhiProLUT::TheLookupFunction(), L1MuGMTMIAUPhiPro1LUT::TheLookupFunction(), reco::IsolatedPixelTrackCandidate::towerIndex(), l1t::CaloTools::towerPhi(), l1t::CaloTools::towerPhiSize(), TrackCount::TrackCount(), TrackerStablePhiSort(), TtFullHadSignalSel::TtFullHadSignalSel(), HcalDDDSimConstants::unitPhi(), unsafe_acos07(), unsafe_acos71(), unsafe_asin71(), SiStripElectronAnalyzer::unwrapPhi(), SiStripElectronSeedGenerator::unwrapPhi(), SiStripElectronAlgo::unwrapPhi(), LaserSteppingAction::UserSteppingAction(), ColinearityKinematicConstraint::value(), L1GctHtMissLut::value(), reco::isodeposit::NumCrystalEtaPhiVeto::veto(), reco::isodeposit::RectangularEtaPhiVeto::veto(), VoronoiAlgorithm::voronoi_area_incident(), and VirtualJetProducer::writeJets().

Function Documentation

double convertEtaFromHW	(	int	hwEta,
		double	max,
		double	step,
		int	hwMax
	)

Definition at line 641 of file readTestVector.cc.

References eta, and relval_parameters_module::step.

Referenced by parseEGs(), parseJets(), parseMuons(), and parseTaus().

                                                                        {
  hwMax++;
  double binWidth = 2*max/step;
  double eta = ( hwEta<int(hwMax/2+1.001) ) ? double(hwEta)*binWidth+0.5*binWidth : -(double(hwMax-hwEta)*binWidth -0.5*binWidth);
  return eta;
}

double convertPhiFromHW	(	int	hwPhi,
		double	step
	)

Definition at line 647 of file readTestVector.cc.

References M_PI, and phi().

Referenced by parseEGs(), parseEtSums(), parseJets(), parseMuons(), and parseTaus().

                                                 {
  double phi = double(hwPhi)/step * 2 * M_PI;
  return phi;
}

double convertPtFromHW	(	int	hwPt,
		double	max,
		double	step
	)

Definition at line 635 of file readTestVector.cc.

References bookConverter::max, EnergyCorrector::pt, and relval_parameters_module::step.

Referenced by parseEGs(), parseEtSums(), parseJets(), parseMuons(), and parseTaus().

                                                           {
  double pt = double(hwPt) * step;
  if( pt>max ) pt = max;
  return pt;
}

int main	(	int	argc,
		char **	argv
	)

Setup histograms

Definition at line 105 of file readTestVector.cc.

References HLT_FULL_cff::algorithm, gather_cfg::cout, EtaStepCalo_, EtaStepMuon_, cppFunctionSkipper::exception, mergeVDriftHistosByStation::file, groupFilesInBlocks::fin, h_l1eg_eta_, h_l1eg_isolation_, h_l1eg_num_, h_l1eg_phi_, h_l1eg_pt_, h_l1et_, h_l1etm_et_, h_l1etm_phi_, h_l1ht_, h_l1htm_et_, h_l1htm_phi_, h_l1jet_eta_, h_l1jet_num_, h_l1jet_phi_, h_l1jet_pt_, h_l1mu_charge_, h_l1mu_eta_, h_l1mu_isolation_, h_l1mu_num_, h_l1mu_phi_, h_l1mu_pt_, h_l1mu_quality_, h_l1tau_eta_, h_l1tau_isolation_, h_l1tau_num_, h_l1tau_phi_, h_l1tau_pt_, jetmet_cfg::histofile, i, cmsHarvester::index, metsig::jet, l1taccepts, l1taccepts_evt, l1tnames, geometryCSVtoXML::line, M_PI, MAX_ALGO_BITS, MaxCaloEta_, MaxEt_, Data_TkAlMuonIsolated_Run2015B_PromptReco_v1_cff::maxEvents, MaxJetPt_, MaxLepPt_, MaxMuonEta_, RPCpg::mu, mergeVDriftHistosByStation::name, NULL, AlCaHLTBitMon_ParallelJobs::options, convertSQLitetoXML_cfg::output, parseAlgo(), parseEGs(), parseEtSums(), parseJets(), parseMuons(), parseTaus(), PhiStepCalo_, PhiStepMuon_, PtStep_, insertMaterial::readXML(), DTTTrigCorrFirst::run, AlCaHLTBitMon_QueryRunRegistry::string, and metsig::tau.

                                 {
  using namespace boost;
  namespace po = boost::program_options;

  std::string vector_file;
  std::string xml_file;
  std::string histo_file;
  bool dumpEvents;
  int maxEvents;

  po::options_description desc("Main options");
  desc.add_options()
    ("vector_file,i", po::value<std::string>(&vector_file)->default_value(""), "Input file")
    ("menu_file,m",   po::value<std::string>(&xml_file)->default_value(""), "Menu file")
    ("hist_file,o", po::value<std::string>(&histo_file)->default_value(""), "Output histogram file")
    ("dumpEvents,d",  po::value<bool>(&dumpEvents)->default_value(false), "Dump event-by-event information")
    ("maxEvents,n",  po::value<int>(&maxEvents)->default_value(-1), "Number of events (default is all)")
    ("help,h", "Produce help message")
    ;

  po::variables_map vm, vm0;

  // parse the first time, using only common options and allow unregistered options 
  try{
    po::store(po::command_line_parser(argc, argv).options(desc).allow_unregistered().run(), vm0);
    po::notify(vm0);
  } catch(std::exception &ex) {
    std::cout << "Invalid options: " << ex.what() << std::endl;
    std::cout << "Use readTestVector --help to get a list of all the allowed options"  << std::endl;
    return 999;
  } catch(...) {
    std::cout << "Unidentified error parsing options." << std::endl;
    return 1000;
  }

  // if help, print help
  if(vm0.count("help")) {
    std::cout << "Usage: readTestVector [options]\n";
    std::cout << desc;
    return 0;
  }

  if( vector_file=="" ){
    std::cout << "No input file specified" << std::endl;
    return 99;
  }

  bool readXML = true;
  if( xml_file=="" ){
    readXML = false;
    std::cout << "No menu file specified" << std::endl;
  }

  bool output = true;
  if( histo_file=="" ){
    output = false;
  }

  TFile* histofile = NULL;
  if( output ){
    histofile = new TFile(histo_file.c_str(),"RECREATE");
    histofile->cd();
  }

  l1taccepts.resize(MAX_ALGO_BITS);
  l1taccepts_evt.resize(MAX_ALGO_BITS);
  for( unsigned int i=0; i<l1taccepts.size(); i++ ) l1taccepts[i] = 0;

  if( readXML ){
    // Load XML.
    std::auto_ptr<l1t::L1TriggerMenu> tm(l1t::l1TriggerMenu(xml_file));

    l1tnames.resize(MAX_ALGO_BITS);
    l1t::AlgorithmList algorithms = tm->algorithms();
    for( l1t::AlgorithmList::algorithm_const_iterator i = algorithms.algorithm().begin();
         i != algorithms.algorithm().end(); ++i ){

      l1t::Algorithm algorithm = (*i);

      int index = algorithm.index();
      std::string name = algorithm.name();

      l1tnames[index] = name;
    }
  }

  h_l1mu_pt_  = new TH1D("h_l1mu_pt", ";L1 #mu p_{T}", int((MaxLepPt_+PtStep_)/(PtStep_) + 1.001), 0, MaxLepPt_+PtStep_ );
  h_l1mu_eta_ = new TH1D("h_l1mu_eta",";L1 #mu #eta",  int(EtaStepMuon_/2+0.0001), -MaxMuonEta_, MaxMuonEta_ );
  h_l1mu_phi_ = new TH1D("h_l1mu_phi",";L1 #mu #phi",  PhiStepMuon_+1, 0, 2*M_PI );
  h_l1mu_charge_ = new TH1D("h_l1mu_charge_",";L1 #mu charge",  2, 0, 2 );
  h_l1mu_quality_ = new TH1D("h_l1mu_quality_",";L1 #mu quality",  16, 0, 16 );
  h_l1mu_isolation_ = new TH1D("h_l1mu_isolation_",";L1 #mu isolation",  4, 0, 4 );
  h_l1mu_num_ = new TH1D("h_l1mu_num",";L1 Number of #mu",  10, 0, 10 );

  h_l1jet_pt_  = new TH1D("h_l1jet_pt", ";L1 jet p_{T}", int((MaxJetPt_+PtStep_)/(4*PtStep_) + 1.001), 0, MaxJetPt_+PtStep_ );
  h_l1jet_eta_ = new TH1D("h_l1jet_eta",";L1 jet #eta",  int(EtaStepCalo_/2+0.0001), -MaxCaloEta_, MaxCaloEta_ );
  h_l1jet_phi_ = new TH1D("h_l1jet_phi",";L1 jet #phi",  PhiStepCalo_+1, 0, 2*M_PI );
  h_l1jet_num_ = new TH1D("h_l1jet_num",";L1 Number of jets",  13, 0, 13 );

  h_l1eg_pt_  = new TH1D("h_l1eg_pt", ";L1 EG p_{T}", int((MaxLepPt_+PtStep_)/(PtStep_) + 1.001), 0, MaxLepPt_+PtStep_ );
  h_l1eg_eta_ = new TH1D("h_l1eg_eta",";L1 EG #eta",  int(EtaStepCalo_/2+0.0001), -MaxCaloEta_, MaxCaloEta_ );
  h_l1eg_phi_ = new TH1D("h_l1eg_phi",";L1 EG #phi",  PhiStepCalo_+1, 0, 2*M_PI );
  h_l1eg_isolation_ = new TH1D("h_l1eg_isolation_",";L1 EG isolation",  4, 0, 4 );
  h_l1eg_num_ = new TH1D("h_l1eg_num",";L1 Number of EGs",  13, 0, 13 );

  h_l1tau_pt_  = new TH1D("h_l1tau_pt", ";L1 #tau p_{T}", int((MaxLepPt_+PtStep_)/(PtStep_) + 1.001), 0, MaxLepPt_+PtStep_ );
  h_l1tau_eta_ = new TH1D("h_l1tau_eta",";L1 #tau #eta",  int(EtaStepCalo_/2+0.0001), -MaxCaloEta_, MaxCaloEta_ );
  h_l1tau_phi_ = new TH1D("h_l1tau_phi",";L1 #tau #phi",  PhiStepCalo_+1, 0, 2*M_PI );
  h_l1tau_isolation_ = new TH1D("h_l1tau_isolation_",";L1 #tau isolation",  4, 0, 4 );
  h_l1tau_num_ = new TH1D("h_l1tau_num",";L1 Number of #tau",  13, 0, 13 );

  h_l1ht_ = new TH1D("h_l1ht_", ";L1 #SigmaH_{T}", int((MaxEt_+PtStep_)/(16*PtStep_) + 1.001), 0, MaxEt_+PtStep_ );
  h_l1et_ = new TH1D("h_l1et_", ";L1 #SigmaE_{T}", int((MaxEt_+PtStep_)/(16*PtStep_) + 1.001), 0, MaxEt_+PtStep_ );
  h_l1htm_et_ = new TH1D("h_l1htm_et_", ";L1 Missing H_{T}", int((MaxEt_+PtStep_)/(16*PtStep_) + 1.001), 0, MaxEt_+PtStep_ );
  h_l1etm_et_ = new TH1D("h_l1etm_et_", ";L1 Missing E_{T}", int((MaxEt_+PtStep_)/(16*PtStep_) + 1.001), 0, MaxEt_+PtStep_ );
  h_l1htm_phi_ = new TH1D("h_l1htm_phi_", ";L1 Missing H_{T} #phi", PhiStepCalo_+1, 0, 2*M_PI );
  h_l1etm_phi_ = new TH1D("h_l1etm_phi_", ";L1 Missing E_{T} #phi", PhiStepCalo_+1, 0, 2*M_PI );


  std::ifstream file(vector_file);
  std::string line;

  std::vector<std::string> mu;
  std::vector<std::string> eg;
  std::vector<std::string> tau;
  std::vector<std::string> jet;
  std::vector<std::string> etsum;
  mu.resize(8);
  eg.resize(12);
  tau.resize(8);
  jet.resize(12);
  etsum.resize(4);

  std::string bx, ex, alg, fin;

  int evt=0;
  int l1a=0;

  if( dumpEvents ) printf(" ==== Objects ===\n");

  while(std::getline(file, line)){
    evt++;
    std::stringstream linestream(line);
    linestream >> bx 
           >> mu[0] >> mu[1] >> mu[2] >> mu[3] >> mu[4] >> mu[5] >> mu[6] >> mu[7] 
           >> eg[0] >> eg[1] >> eg[2] >> eg[3] >> eg[4] >> eg[5] >> eg[6] >> eg[7] >> eg[8] >> eg[9] >> eg[10] >> eg[11] 
           >> tau[0] >> tau[1] >> tau[2] >> tau[3] >> tau[4] >> tau[5] >> tau[6] >> tau[7] 
           >> jet[0] >> jet[1] >> jet[2] >> jet[3] >> jet[4] >> jet[5] >> jet[6] >> jet[7] >> jet[8] >> jet[9] >> jet[10] >> jet[11] 
           >> etsum[0] >> etsum[1] >> etsum[2] >> etsum[3] 
           >> ex >> alg >> fin;

    if( dumpEvents ) printf("  <==== BX = %s ====>\n",bx.c_str());
    parseMuons(mu, dumpEvents);
    parseEGs(eg, dumpEvents);
    parseTaus(tau, dumpEvents);
    parseJets(jet, dumpEvents);
    parseEtSums(etsum, dumpEvents);

    parseAlgo(alg);

    int finOR = atoi( fin.c_str() );
    l1a += finOR;

    if( dumpEvents ){
      printf("    == Algos ==\n");
      if( finOR ){
    printf(" Triggers with non-zero accepts\n");
    if( readXML && l1tnames.size()>0 ) printf("\t bit\t L1A\t Name\n");
    else                               printf("\t bit\t L1A\n");

    for( int i=0; i<MAX_ALGO_BITS; i++ ){
      if( l1taccepts_evt[i]>0 ){
        if( readXML && l1tnames.size()>0 ) printf("\t %d \t %d \t %s\n", i, l1taccepts_evt[i], l1tnames[i].c_str());
        else                               printf("\t %d \t %d\n", i, l1taccepts_evt[i] );
      }
    }
      }
      else{
    printf("\n No triggers passed\n");
      }
      // extra spacing between bx
      printf("\n\n");
    }

    if( evt==maxEvents ) break;
  }

  printf(" =========== Summary of results ==========\n");
  printf(" There were %d L1A out of %d events (%.1f%%)\n", l1a, evt, float(l1a)/float(evt)*100);
  printf("\n Triggers with non-zero accepts\n");
  if( readXML && l1tnames.size()>0 ) printf("\t bit\t L1A\t Name\n");
  else                               printf("\t bit\t L1A\n");

  for( int i=0; i<MAX_ALGO_BITS; i++ ){
    if( l1taccepts[i]>0 ){
      if( readXML && l1tnames.size()>0 ) printf("\t %d \t %d \t %s\n", i, l1taccepts[i], l1tnames[i].c_str());
      else                               printf("\t %d \t %d\n", i, l1taccepts[i] );
    }
  }

  /*
  printf(" =========== Dump of results ==========\n");
  for( int i=0; i<MAX_ALGO_BITS; i++ ){
     printf("%d\n", l1taccepts[i]);
  }
  */


  /*
  int NSet=13;
  printf("# <algo_id>");
  for( int iSet=0; iSet<NSet; iSet++ ){
    printf(", <set_%d>",iSet+1);
  }
  printf("\n");

  printf("-1");
  for( int iSet=0; iSet<NSet; iSet++ ){
    int useSet = 0;
    if( iSet==0 )      useSet = -3;
    else if( iSet==1 ) useSet = -2;
    else               useSet = iSet-1;
    printf(", %d",useSet);
  }
  printf("\n");

  for( int i=0; i<MAX_ALGO_BITS; i++ ){
    printf("%d",i);
    for( int iSet=0; iSet<NSet; iSet++ ){
      int prescale = 1;
      if( iSet==0 )       prescale = 0;
      else if( iSet==1 )  prescale = 1;
      else if( iSet==2 )  prescale = 1;
      else if( iSet==3 )  prescale = 2;
      else if( iSet==4 )  prescale = 10;
      else if( iSet==5 && i<MAX_ALGO_BITS/2 ) prescale = 2;
      else if( iSet==6 && i>MAX_ALGO_BITS/2 ) prescale = 5;
      else if( iSet==7 && i%2==0 ) prescale = 3;
      else if( iSet==8 && i%2==1 ) prescale = 6;
      else if( iSet==9 && i%3==0 ) prescale = 4;
      else if( iSet==10 ) prescale = 10;
      else if( iSet==11 ) prescale = 2;
      else if( iSet==12 ) prescale = 1;

      printf(", %d",prescale);
    }
    printf("\n");
  }
  */

  if( output ){
    histofile->Write();
    histofile->Close();
  }


  return 0;
}

void parseAlgo

(

std::string

algo

)

Definition at line 520 of file readTestVector.cc.

References accept(), i, recoMuon::in, l1taccepts, l1taccepts_evt, MAX_ALGO_BITS, alignCSCRings::s, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, and zeroes128.

Referenced by main().

                                {

  // clear the algo per event
  l1taccepts_evt.clear();
  l1taccepts_evt.resize(MAX_ALGO_BITS);

  // Return if no triggers fired
  if( algo==zeroes128 ) return;

  std::stringstream ss;
  std::string s;
  int pos, algRes, accept;
  
  for( int i=0; i<MAX_ALGO_BITS; i++ ){
    pos = 127 - int( i/4 );
    std::string in(1,algo[pos]);
    char* endPtr = (char*) in.c_str();
    algRes = strtol( in.c_str(), &endPtr, 16 );
    accept = ( algRes >> (i%4) ) & 1;
    l1taccepts_evt[i] += accept;
    l1taccepts[i] += accept;
  }

  return;

}

void parseEGs	(	std::vector< std::string >	egs,
		bool	verbose
	)

Definition at line 400 of file readTestVector.cc.

References convertEtaFromHW(), convertPhiFromHW(), convertPtFromHW(), eta, EtaStepCalo_, h_l1eg_eta_, h_l1eg_isolation_, h_l1eg_num_, h_l1eg_phi_, h_l1eg_pt_, l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwIso(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), i, MaxCaloEta_, MaxLepPt_, phi(), PhiStepCalo_, EnergyCorrector::pt, PtStep_, AlCaHLTBitMon_QueryRunRegistry::string, unpackEGs(), and zeroes8.

Referenced by main().

                                                       {

  if( verbose) printf("    == EGammas ==\n");
  int neg=0;
  for( unsigned int i=0; i<egs.size(); i++ ){
    std::string ieg = egs[i];
    if( ieg==zeroes8 ) continue;
    neg++;
    l1t::EGamma eg = unpackEGs( ieg );

    double pt = convertPtFromHW( eg.hwPt(), MaxLepPt_, PtStep_ );
    double eta = convertEtaFromHW( eg.hwEta(), MaxCaloEta_, EtaStepCalo_, 0xff );
    double phi = convertPhiFromHW( eg.hwPhi(), PhiStepCalo_ );

    h_l1eg_pt_->Fill( pt );
    h_l1eg_eta_->Fill( eta );
    h_l1eg_phi_->Fill( phi );
    h_l1eg_isolation_->Fill( eg.hwIso() );

    if( verbose) printf(" l1t::EGamma %d:\t pt = %d (%.1f),\t eta = %d (%+.2f),\t phi = %d (%.2f),\t iso = %d\n", i, eg.hwPt(), pt, eg.hwEta(), eta, eg.hwPhi(), phi, eg.hwIso());
  }
  h_l1eg_num_->Fill(neg);

  return;
}

void parseEtSums	(	std::vector< std::string >	etsums,
		bool	verbose
	)

Definition at line 474 of file readTestVector.cc.

References convertPhiFromHW(), convertPtFromHW(), h_l1et_, h_l1etm_et_, h_l1etm_phi_, h_l1ht_, h_l1htm_et_, h_l1htm_phi_, l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), L1Analysis::kMissingEt, L1Analysis::kMissingHt, L1Analysis::kTotalEt, L1Analysis::kTotalHt, MaxEt_, phi(), PhiStepCalo_, EnergyCorrector::pt, PtStep_, AlCaHLTBitMon_QueryRunRegistry::string, unpackEtSums(), and zeroes8.

Referenced by main().

                                                            {

  if( verbose) printf("    == EtSums ==\n");

  //et sum
  std::string iet = etsum[0];
  if( iet!=zeroes8 ){
    l1t::EtSum et = unpackEtSums( iet, l1t::EtSum::EtSumType::kTotalEt );
    double pt = convertPtFromHW( et.hwPt(), MaxEt_, PtStep_ );
    h_l1et_->Fill( pt );
    if( verbose) printf(" l1t::EtSum TotalEt:\t Et = %d\n", et.hwPt());
  }
  //ht sum
  std::string iht = etsum[1];
  if( iht!=zeroes8 ){
    l1t::EtSum ht = unpackEtSums( iht, l1t::EtSum::EtSumType::kTotalHt );
    double pt = convertPtFromHW( ht.hwPt(), MaxEt_, PtStep_ );
    h_l1ht_->Fill( pt );
    if( verbose) printf(" l1t::EtSum TotalHt:\t Ht = %d\n", ht.hwPt());
  }
  //etm
  std::string ietm = etsum[2];
  if( ietm!=zeroes8 ){
    l1t::EtSum etm = unpackEtSums( ietm, l1t::EtSum::EtSumType::kMissingEt );
    double pt = convertPtFromHW( etm.hwPt(), MaxEt_, PtStep_ );
    double phi = convertPhiFromHW( etm.hwPhi(), PhiStepCalo_ );
    h_l1etm_et_->Fill( pt );
    h_l1etm_phi_->Fill( phi );
    if( verbose) printf(" l1t::EtSum MissingEt:\t Et = %d,\t phi = %d\n", etm.hwPt(), etm.hwPhi());
  }
  //htm
  std::string ihtm = etsum[3];
  if( ihtm!=zeroes8 ){
    l1t::EtSum htm = unpackEtSums( ihtm, l1t::EtSum::EtSumType::kMissingHt );
    double pt = convertPtFromHW( htm.hwPt(), MaxEt_, PtStep_ );
    double phi = convertPhiFromHW( htm.hwPhi(), PhiStepCalo_ );
    h_l1htm_et_->Fill( pt );
    h_l1htm_phi_->Fill( phi );
    if( verbose) printf(" l1t::EtSum MissingHt:\t Et = %d,\t phi = %d\n", htm.hwPt(), htm.hwPhi());
  }

  return;
}

void parseJets	(	std::vector< std::string >	jets,
		bool	verbose
	)

Definition at line 450 of file readTestVector.cc.

References convertEtaFromHW(), convertPhiFromHW(), convertPtFromHW(), eta, EtaStepCalo_, h_l1jet_eta_, h_l1jet_num_, h_l1jet_phi_, h_l1jet_pt_, l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), i, metsig::jet, MaxCaloEta_, MaxJetPt_, njet, phi(), PhiStepCalo_, EnergyCorrector::pt, PtStep_, AlCaHLTBitMon_QueryRunRegistry::string, unpackJets(), and zeroes8.

Referenced by main().

                                                         {

  if( verbose) printf("    == Jets ==\n");
  int njet=0;
  for( unsigned int i=0; i<jets.size(); i++ ){
    std::string ijet = jets[i];
    if( ijet==zeroes8 ) continue;
    njet++;
    l1t::Jet jet = unpackJets( ijet );

    double pt = convertPtFromHW( jet.hwPt(), MaxJetPt_, PtStep_ );
    double eta = convertEtaFromHW( jet.hwEta(), MaxCaloEta_, EtaStepCalo_, 0xff );
    double phi = convertPhiFromHW( jet.hwPhi(), PhiStepCalo_ );

    h_l1jet_pt_->Fill( pt );
    h_l1jet_eta_->Fill( eta );
    h_l1jet_phi_->Fill( phi );

    if( verbose) printf(" l1t::Jet %d:\t pt = %d (%.1f),\t eta = %d (%+.2f),\t phi = %d (%.2f)\n", i, jet.hwPt(), pt, jet.hwEta(), eta, jet.hwPhi(), phi);
  }
  h_l1jet_num_->Fill(njet);

  return;
}

void parseMuons	(	std::vector< std::string >	muons,
		bool	verbose
	)

Definition at line 367 of file readTestVector.cc.

References RecoTauCleanerPlugins::charge, convertEtaFromHW(), convertPhiFromHW(), convertPtFromHW(), eta, EtaStepMuon_, h_l1mu_charge_, h_l1mu_eta_, h_l1mu_isolation_, h_l1mu_num_, h_l1mu_phi_, h_l1mu_pt_, h_l1mu_quality_, l1t::Muon::hwCharge(), l1t::Muon::hwChargeValid(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwIso(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), i, MaxLepPt_, MaxMuonEta_, RPCpg::mu, phi(), PhiStepMuon_, EnergyCorrector::pt, PtStep_, AlCaHLTBitMon_QueryRunRegistry::string, unpackMuons(), and zeroes16.

Referenced by main().

                                                           {

  if( verbose) printf("    == Muons ==\n");
  int nmu=0;
  for( unsigned int i=0; i<muons.size(); i++ ){
    std::string imu = muons[i];
    if( imu==zeroes16 ) continue;
    nmu++;
    l1t::Muon mu = unpackMuons( imu );

    double pt = convertPtFromHW( mu.hwPt(), MaxLepPt_, PtStep_ );
    double eta = convertEtaFromHW( mu.hwEta(), MaxMuonEta_, EtaStepMuon_, 0x1ff );
    double phi = convertPhiFromHW( mu.hwPhi(), PhiStepMuon_ );

    int iso = mu.hwIso();
    int qual = mu.hwQual();
    int charge = mu.hwCharge();
    int chargeValid = mu.hwChargeValid();

    h_l1mu_pt_->Fill( pt );
    h_l1mu_eta_->Fill( eta );
    h_l1mu_phi_->Fill( phi );
    h_l1mu_charge_->Fill( charge );

    h_l1mu_quality_->Fill( qual );
    h_l1mu_isolation_->Fill( iso );

    if( verbose) printf(" l1t::Muon %d:\t pt = %d (%.1f),\t eta = %d (%+.2f),\t phi = %d (%.2f),\t iso = %d,\t qual = %d,\t charge = %d,\t chargeValid = %d\n", i, mu.hwPt(), pt, mu.hwEta(), eta, mu.hwPhi(), phi, iso, qual, charge, chargeValid);
  }
  h_l1mu_num_->Fill(nmu);

  return;
}

void parseTaus	(	std::vector< std::string >	taus,
		bool	verbose
	)

Definition at line 425 of file readTestVector.cc.

References convertEtaFromHW(), convertPhiFromHW(), convertPtFromHW(), eta, EtaStepCalo_, h_l1tau_eta_, h_l1tau_isolation_, h_l1tau_num_, h_l1tau_phi_, h_l1tau_pt_, l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwIso(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), i, MaxCaloEta_, MaxLepPt_, phi(), PhiStepCalo_, EnergyCorrector::pt, PtStep_, AlCaHLTBitMon_QueryRunRegistry::string, metsig::tau, unpackTaus(), and zeroes8.

Referenced by main().

                                                         {

  if( verbose) printf("    == Taus ==\n");
  int ntau=0;
  for( unsigned int i=0; i<taus.size(); i++ ){
    std::string itau = taus[i];
    if( itau==zeroes8 ) continue;
    ntau++;
    l1t::Tau tau = unpackTaus( itau );

    double pt = convertPtFromHW( tau.hwPt(), MaxLepPt_, PtStep_ );
    double eta = convertEtaFromHW( tau.hwEta(), MaxCaloEta_, EtaStepCalo_, 0xff );
    double phi = convertPhiFromHW( tau.hwPhi(), PhiStepCalo_ );

    h_l1tau_pt_->Fill( pt );
    h_l1tau_eta_->Fill( eta );
    h_l1tau_phi_->Fill( phi );
    h_l1tau_isolation_->Fill( tau.hwIso() );

    if( verbose) printf(" l1t::Tau %d:\t pt = %d (%.1f),\t eta = %d (%+.2f),\t phi = %d (%.2f),\t iso = %d\n", i, tau.hwPt(), pt, tau.hwEta(), eta, tau.hwPhi(), phi, tau.hwIso());
  }
  h_l1tau_num_->Fill(ntau);

  return;
}

l1t::EGamma unpackEGs

(

std::string

ieg

)

Definition at line 569 of file readTestVector.cc.

References eta, p4, phi(), and EnergyCorrector::pt.

Referenced by parseEGs().

                                    {

  char* endPtr = (char*) ieg.c_str();
  unsigned int packedVal = strtol( ieg.c_str(), &endPtr, 16 );

  int pt  = (packedVal>>0)  & 0x1ff;
  int eta = (packedVal>>9)  & 0xff;
  int phi = (packedVal>>17) & 0xff;
  int iso = (packedVal>>25) & 0x3;
  int qual= (packedVal>>27) & 0x31;

  ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > *p4 = new ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> >();
  l1t::EGamma eg(*p4, pt, eta, phi, qual, iso);

  return eg;
}

l1t::EtSum unpackEtSums	(	std::string	ietsum,
		l1t::EtSum::EtSumType	type
	)

Definition at line 618 of file readTestVector.cc.

References L1Analysis::kMissingEt, L1Analysis::kMissingHt, p4, phi(), and EnergyCorrector::pt.

Referenced by parseEtSums().

                                                                 {

  char* endPtr = (char*) ietsum.c_str();
  unsigned int packedVal = strtol( ietsum.c_str(), &endPtr, 16 );

  int pt  = (packedVal>>0)  & 0xfff;
  int phi = 0;
  if( type==l1t::EtSum::EtSumType::kMissingEt ||
      type==l1t::EtSum::EtSumType::kMissingHt ) phi = (packedVal>>12) & 0xff;
  
  ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > *p4 = new ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> >();
  l1t::EtSum sum(*p4, type, pt, 0,phi, 0); 

  return sum;
}

l1t::Jet unpackJets

(

std::string

ijet

)

Definition at line 602 of file readTestVector.cc.

References eta, metsig::jet, p4, phi(), and EnergyCorrector::pt.

Referenced by parseJets().

                                   {

  char* endPtr = (char*) ijet.c_str();
  unsigned int packedVal = strtol( ijet.c_str(), &endPtr, 16 );

  int pt  = (packedVal>>0)  & 0x7ff;
  int eta = (packedVal>>11) & 0xff;
  int phi = (packedVal>>19) & 0xff;
  int qual= (packedVal>>27) & 0x1;

  ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > *p4 = new ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> >();
  l1t::Jet jet(*p4, pt, eta, phi, qual);

  return jet;
}

l1t::Muon unpackMuons

(

std::string

imu

)

Definition at line 548 of file readTestVector.cc.

References RecoTauCleanerPlugins::charge, eta, RPCpg::mu, p4, phi(), EnergyCorrector::pt, and GlobalPosition_Frontier_DevDB_cff::tag.

Referenced by parseMuons().

                                    {

  char* endPtr = (char*) imu.c_str();
  cms_uint64_t packedVal = strtoull( imu.c_str(), &endPtr, 16 );

  int pt  = (packedVal>>10) & 0x1ff;
  int eta = (packedVal>>23) & 0x1ff;
  int phi = (packedVal>>0)  & 0x3ff;
  int iso = (packedVal>>32) & 0x3;
  int qual= (packedVal>>19) & 0xf;
  int charge = (packedVal>>34) & 0x1;
  int chargeValid = (packedVal>>35) & 0x1;
  int mip = 1;
  int tag = 1;

  ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > *p4 = new ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> >();
  l1t::Muon mu(*p4, pt, eta, phi, qual, charge, chargeValid, iso, mip, tag);

  return mu;

}

l1t::Tau unpackTaus

(

std::string

itau

)

Definition at line 585 of file readTestVector.cc.

References eta, p4, phi(), EnergyCorrector::pt, and metsig::tau.

Referenced by parseTaus().

                                   {

  char* endPtr = (char*) itau.c_str();
  unsigned int packedVal = strtol( itau.c_str(), &endPtr, 16 );

  int pt  = (packedVal>>0)  & 0x1ff;
  int eta = (packedVal>>9)  & 0xff;
  int phi = (packedVal>>17) & 0xff;
  int iso = (packedVal>>25) & 0x3;
  int qual= (packedVal>>27) & 0x31;

  ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > *p4 = new ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> >();
  l1t::Tau tau(*p4, pt, eta, phi, qual, iso);

  return tau;
}

Variable Documentation

double EtaStepCalo_ = 230

Definition at line 98 of file readTestVector.cc.

Referenced by main(), parseEGs(), parseJets(), parseTaus(), and l1t::GenToInputProducer::produce().

double EtaStepMuon_ = 450

Definition at line 99 of file readTestVector.cc.

Referenced by main(), parseMuons(), and l1t::GenToInputProducer::produce().

TH1D* h_l1eg_eta_

Definition at line 69 of file readTestVector.cc.

Referenced by main(), and parseEGs().

TH1D* h_l1eg_isolation_

Definition at line 71 of file readTestVector.cc.

Referenced by main(), and parseEGs().

TH1D* h_l1eg_num_

Definition at line 72 of file readTestVector.cc.

Referenced by main(), and parseEGs().

TH1D* h_l1eg_phi_

Definition at line 70 of file readTestVector.cc.

Referenced by main(), and parseEGs().

TH1D* h_l1eg_pt_

Definition at line 68 of file readTestVector.cc.

Referenced by main(), and parseEGs().

TH1D* h_l1et_

Definition at line 81 of file readTestVector.cc.

Referenced by main(), and parseEtSums().

TH1D* h_l1etm_et_

Definition at line 83 of file readTestVector.cc.

Referenced by main(), and parseEtSums().

TH1D* h_l1etm_phi_

Definition at line 85 of file readTestVector.cc.

Referenced by main(), and parseEtSums().

TH1D* h_l1ht_

Definition at line 80 of file readTestVector.cc.

Referenced by main(), and parseEtSums().

TH1D* h_l1htm_et_

Definition at line 82 of file readTestVector.cc.

Referenced by main(), and parseEtSums().

TH1D* h_l1htm_phi_

Definition at line 84 of file readTestVector.cc.

Referenced by main(), and parseEtSums().

TH1D* h_l1jet_eta_

Definition at line 64 of file readTestVector.cc.

Referenced by main(), and parseJets().

TH1D* h_l1jet_num_

Definition at line 66 of file readTestVector.cc.

Referenced by main(), and parseJets().

TH1D* h_l1jet_phi_

Definition at line 65 of file readTestVector.cc.

Referenced by main(), and parseJets().

TH1D* h_l1jet_pt_

Definition at line 63 of file readTestVector.cc.

Referenced by main(), and parseJets().

TH1D* h_l1mu_charge_

Definition at line 58 of file readTestVector.cc.

Referenced by main(), and parseMuons().

TH1D* h_l1mu_eta_

Definition at line 56 of file readTestVector.cc.

Referenced by main(), and parseMuons().

TH1D* h_l1mu_isolation_

Definition at line 60 of file readTestVector.cc.

Referenced by main(), and parseMuons().

TH1D* h_l1mu_num_

Definition at line 61 of file readTestVector.cc.

Referenced by main(), and parseMuons().

TH1D* h_l1mu_phi_

Definition at line 57 of file readTestVector.cc.

Referenced by main(), and parseMuons().

TH1D* h_l1mu_pt_

Definition at line 55 of file readTestVector.cc.

Referenced by main(), and parseMuons().

TH1D* h_l1mu_quality_

Definition at line 59 of file readTestVector.cc.

Referenced by main(), and parseMuons().

TH1D* h_l1tau_eta_

Definition at line 75 of file readTestVector.cc.

Referenced by main(), and parseTaus().

TH1D* h_l1tau_isolation_

Definition at line 77 of file readTestVector.cc.

Referenced by main(), and parseTaus().

TH1D* h_l1tau_num_

Definition at line 78 of file readTestVector.cc.

Referenced by main(), and parseTaus().

TH1D* h_l1tau_phi_

Definition at line 76 of file readTestVector.cc.

Referenced by main(), and parseTaus().

TH1D* h_l1tau_pt_

Definition at line 74 of file readTestVector.cc.

Referenced by main(), and parseTaus().

std::vector<int> l1taccepts

Definition at line 32 of file readTestVector.cc.

Referenced by main(), and parseAlgo().

std::vector<int> l1taccepts_evt

Definition at line 34 of file readTestVector.cc.

Referenced by main(), and parseAlgo().

std::vector<std::string> l1tnames

Definition at line 33 of file readTestVector.cc.

Referenced by main().

const int MAX_ALGO_BITS = 512

Definition at line 22 of file readTestVector.cc.

Referenced by main(), and parseAlgo().

double MaxCaloEta_ = 5.0

Definition at line 92 of file readTestVector.cc.

Referenced by main(), parseEGs(), parseJets(), parseTaus(), and l1t::GenToInputProducer::produce().

double MaxEt_ = 2047

Definition at line 90 of file readTestVector.cc.

Referenced by main(), parseEtSums(), and l1t::GenToInputProducer::produce().

double MaxJetPt_ = 1023

Definition at line 89 of file readTestVector.cc.

Referenced by main(), parseJets(), and l1t::GenToInputProducer::produce().

double MaxLepPt_ = 255

Definition at line 88 of file readTestVector.cc.

Referenced by main(), parseEGs(), parseMuons(), parseTaus(), and l1t::GenToInputProducer::produce().

double MaxMuonEta_ = 2.45

Definition at line 93 of file readTestVector.cc.

Referenced by main(), parseMuons(), and l1t::GenToInputProducer::produce().

double PhiStepCalo_ = 144

Definition at line 95 of file readTestVector.cc.

Referenced by main(), parseEGs(), parseEtSums(), parseJets(), parseTaus(), and l1t::GenToInputProducer::produce().

double PhiStepMuon_ = 576

Definition at line 96 of file readTestVector.cc.

Referenced by main(), parseMuons(), and l1t::GenToInputProducer::produce().

double PtStep_ = 0.5

Definition at line 101 of file readTestVector.cc.

Referenced by main(), parseEGs(), parseEtSums(), parseJets(), parseMuons(), parseTaus(), and l1t::GenToInputProducer::produce().

std::string zeroes128 = std::string(128,'0')

Definition at line 28 of file readTestVector.cc.

Referenced by parseAlgo().

std::string zeroes16 = std::string(16,'0')

Definition at line 29 of file readTestVector.cc.

Referenced by parseMuons().

std::string zeroes8 = std::string(8,'0')

Definition at line 30 of file readTestVector.cc.

Referenced by parseEGs(), parseEtSums(), parseJets(), and parseTaus().