DQMPFCandidateAnalyzer Class Reference

#include <DQMPFCandidateAnalyzer.h>

Inheritance diagram for DQMPFCandidateAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	Get the analysis. More...
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	Inizialize parameters for histo binning. More...
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	Initialize run-based parameters. More...
	DQMPFCandidateAnalyzer (const edm::ParameterSet &)
	Constructor. More...
void	endRun (const edm::Run &, const edm::EventSetup &) override
	Finish up a run. More...
	~DQMPFCandidateAnalyzer () override
	Destructor. More...
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Static Private Member Functions
static bool	jetSortingRule (reco::Jet x, reco::Jet y)

Private Attributes
bool	bypassAllDCSChecks_
bool	bypassAllPVChecks_
std::string	candidateType_
edm::ParameterSet	cleaningParameters_
std::vector< int >	countsPFCand_
std::vector< int >	countsPFCandRECO_
JetMETDQMDCSFilter *	DCSFilter_
std::vector< edm::ParameterSet >	diagnosticsParameters_
std::vector< double >	etaMaxPFCand_
std::vector< double >	etaMaxPFCandRECO_
std::vector< double >	etaMinPFCand_
std::vector< double >	etaMinPFCandRECO_
HLTConfigProvider	FilterhltConfig_
bool	hbhenoifilterdecision
edm::InputTag	hbheNoiseFilterResultTag_
edm::EDGetTokenT< bool >	hbheNoiseFilterResultToken_
std::string	HBHENoiseStringMiniAOD
double	hcalMin_
bool	isMiniAO_
bool	isMiniAOD_
int	LSBegin_
int	LSEnd_
MonitorElement *	m_HOverTrackP_Barrel_hPt_10_20
MonitorElement *	m_HOverTrackP_Barrel_hPt_1_10
MonitorElement *	m_HOverTrackP_Barrel_hPt_20_50
MonitorElement *	m_HOverTrackP_Barrel_hPt_50
MonitorElement *	m_HOverTrackP_EndCap_hPt_10_20
MonitorElement *	m_HOverTrackP_EndCap_hPt_1_10
MonitorElement *	m_HOverTrackP_EndCap_hPt_20_50
MonitorElement *	m_HOverTrackP_EndCap_hPt_50
MonitorElement *	m_HOverTrackP_trackPtVsEta
MonitorElement *	m_HOverTrackPVsEta_hPt_10_20
MonitorElement *	m_HOverTrackPVsEta_hPt_1_10
MonitorElement *	m_HOverTrackPVsEta_hPt_20_50
MonitorElement *	m_HOverTrackPVsEta_hPt_50
MonitorElement *	m_HOverTrackPVsTrackP_Barrel
MonitorElement *	m_HOverTrackPVsTrackP_EndCap
MonitorElement *	m_HOverTrackPVsTrackPt_Barrel
MonitorElement *	m_HOverTrackPVsTrackPt_EndCap
std::map< std::string, MonitorElement * >	map_of_MEs
edm::InputTag	METFilterMiniAODLabel2_
edm::InputTag	METFilterMiniAODLabel_
edm::EDGetTokenT< edm::TriggerResults >	METFilterMiniAODToken2_
edm::EDGetTokenT< edm::TriggerResults >	METFilterMiniAODToken_
int	miniaodfilterdec
int	miniaodfilterindex
edm::InputTag	mInputCollection_
MonitorElement *	mProfileIsoPFChHad_EcalOccupancyCentral
MonitorElement *	mProfileIsoPFChHad_EcalOccupancyEndcap
MonitorElement *	mProfileIsoPFChHad_EMPtCentral
MonitorElement *	mProfileIsoPFChHad_EMPtEndcap
MonitorElement *	mProfileIsoPFChHad_HadPtCentral
MonitorElement *	mProfileIsoPFChHad_HadPtEndcap
MonitorElement *	mProfileIsoPFChHad_HcalOccupancyCentral
MonitorElement *	mProfileIsoPFChHad_HcalOccupancyEndcap
MonitorElement *	mProfileIsoPFChHad_TrackOccupancy
MonitorElement *	mProfileIsoPFChHad_TrackPt
std::vector< MonitorElement * >	multiplicityPFCand_
std::vector< std::string >	multiplicityPFCand_name_
std::vector< std::string >	multiplicityPFCand_nameRECO_
std::vector< MonitorElement * >	multiplicityPFCandRECO_
int	numPV_
std::vector< MonitorElement * >	occupancyPFCand_
std::vector< std::string >	occupancyPFCand_name_
std::vector< std::string >	occupancyPFCand_name_puppiNolepWeight_
std::vector< std::string >	occupancyPFCand_nameRECO_
std::vector< MonitorElement * >	occupancyPFCand_puppiNolepWeight_
std::vector< MonitorElement * >	occupancyPFCandRECO_
edm::EDGetTokenT< std::vector< pat::PackedCandidate > >	pflowPackedToken_
edm::EDGetTokenT< std::vector< reco::PFCandidate > >	pflowToken_
double	ptMinCand_
std::vector< MonitorElement * >	ptPFCand_
std::vector< std::string >	ptPFCand_name_
std::vector< std::string >	ptPFCand_name_puppiNolepWeight_
std::vector< std::string >	ptPFCand_nameRECO_
std::vector< MonitorElement * >	ptPFCand_puppiNolepWeight_
std::vector< MonitorElement * >	ptPFCandRECO_
edm::InputTag	theTriggerResultsLabel_
std::vector< int >	typePFCand_
std::vector< int >	typePFCandRECO_
int	verbose_
edm::InputTag	vertexTag_
edm::EDGetTokenT< std::vector< reco::Vertex > >	vertexToken_

Detailed Description

DQM jetMET analysis monitoring for PFCandidates

    Jan. '16: by

    M. Artur Weber

Definition at line 66 of file DQMPFCandidateAnalyzer.h.

Constructor & Destructor Documentation

DQMPFCandidateAnalyzer::DQMPFCandidateAnalyzer

(

const edm::ParameterSet &

pSet

)

Constructor.

Definition at line 60 of file DQMPFCandidateAnalyzer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

{

  miniaodfilterdec=-1;

  candidateType_= pSet.getUntrackedParameter<std::string>("CandType");
  //here only choice between miniaod or reco

  LSBegin_     = pSet.getParameter<int>("LSBegin");
  LSEnd_       = pSet.getParameter<int>("LSEnd");

  isMiniAOD_= (std::string("Packed") == candidateType_);

  mInputCollection_       = pSet.getParameter<edm::InputTag>("PFCandidateLabel");

  if(isMiniAOD_){
    pflowPackedToken_ = consumes<std::vector<pat::PackedCandidate> >(mInputCollection_);
  }else{
    pflowToken_ = consumes<std::vector<reco::PFCandidate> >(mInputCollection_);
  }
  
  miniaodfilterdec=-1;

 // Smallest track pt
  ptMinCand_ = pSet.getParameter<double>("ptMinCand");

  // Smallest raw HCAL energy linked to the track
  hcalMin_ = pSet.getParameter<double>("hcalMin");

  diagnosticsParameters_ = pSet.getParameter<std::vector<edm::ParameterSet> >("METDiagonisticsParameters");

  edm::ConsumesCollector iC  = consumesCollector();
  //DCS
  DCSFilter_ = new JetMETDQMDCSFilter(pSet.getParameter<ParameterSet>("DCSFilter"), iC );
  if(isMiniAOD_){
    METFilterMiniAODLabel_=pSet.getParameter<edm::InputTag>("FilterResultsLabelMiniAOD");
    METFilterMiniAODToken_=consumes<edm::TriggerResults>(METFilterMiniAODLabel_);
    
    METFilterMiniAODLabel2_=pSet.getParameter<edm::InputTag>("FilterResultsLabelMiniAOD2");
    METFilterMiniAODToken2_=consumes<edm::TriggerResults>(METFilterMiniAODLabel2_);


    HBHENoiseStringMiniAOD = pSet.getParameter<std::string>("HBHENoiseLabelMiniAOD");
  }

  if(!isMiniAOD_){
    hbheNoiseFilterResultTag_    = pSet.getParameter<edm::InputTag>("HBHENoiseFilterResultLabel");
    hbheNoiseFilterResultToken_=consumes<bool>(hbheNoiseFilterResultTag_);
  }
  //jet cleanup parameters
  cleaningParameters_ = pSet.getParameter<ParameterSet>("CleaningParameters");

  //Vertex requirements
  bypassAllPVChecks_    = cleaningParameters_.getParameter<bool>("bypassAllPVChecks");
  bypassAllDCSChecks_    = cleaningParameters_.getParameter<bool>("bypassAllDCSChecks");
  vertexTag_    = cleaningParameters_.getParameter<edm::InputTag>("vertexCollection");
  vertexToken_  = consumes<std::vector<reco::Vertex> >(edm::InputTag(vertexTag_));
  
  verbose_= pSet.getParameter<int>("verbose");
}  

DQMPFCandidateAnalyzer::~DQMPFCandidateAnalyzer ( )

override

Destructor.

Definition at line 125 of file DQMPFCandidateAnalyzer.cc.

References LogTrace.

                                                {

  delete DCSFilter_;
  LogTrace("DQMPFCandidateAnalyzer")<<"[DQMPFCandidateAnalyzer] Saving the histos";
}

Member Function Documentation

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

override

Get the analysis.

Definition at line 411 of file DQMPFCandidateAnalyzer.cc.

References funct::abs(), edm::HLTGlobalStatus::accept(), EnergyCorrector::c, reco::PFBlockElementCluster::clusterRef(), gather_cfg::cout, cmsPerfPublish::DirName, ECAL, allElectronIsolations_cfi::elements, reco::PFCandidate::elementsInBlocks(), stringResolutionProvider_cfi::et, reco::LeafCandidate::eta(), pat::PackedCandidate::eta(), HcalObjRepresent::Fill(), edm::Event::getByToken(), HCAL, reco::PFCandidate::hcalEnergy(), mps_fire::i, cuy::ii, edm::HandleBase::isValid(), LogDebug, edm::EventBase::luminosityBlock(), nTracks(), pfLinker_cff::particleFlow, reco::PFCandidate::particleId(), pat::PackedCandidate::pdgId(), reco::LeafCandidate::phi(), pat::PackedCandidate::phi(), edm::Handle< T >::product(), reco::LeafCandidate::pt(), pat::PackedCandidate::pt(), pat::PackedCandidate::puppiWeight(), pat::PackedCandidate::puppiWeightNoLep(), reco::PFCandidate::rawEcalEnergy(), reco::PFCandidate::rawHcalEnergy(), edm::OwnVector< T, P >::size(), AlCaHLTBitMon_QueryRunRegistry::string, reco::PFBlockElementTrack::trackRef(), and particleFlowSuperClusterECAL_cfi::vertexCollection.

                                                                                        {


  //Vertex information
  Handle<VertexCollection> vertexHandle;
  iEvent.getByToken(vertexToken_, vertexHandle);

  if (!vertexHandle.isValid()) {
    LogDebug("") << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
    if (verbose_) std::cout << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
  }
  numPV_ = 0;
  if ( vertexHandle.isValid() ){
    VertexCollection vertexCollection = *(vertexHandle.product());
    numPV_  = vertexCollection.size();
  }
  bool bPrimaryVertex = (bypassAllPVChecks_ || (numPV_>0));


  int myLuminosityBlock;
  myLuminosityBlock = iEvent.luminosityBlock();

  if (myLuminosityBlock<LSBegin_) return;
  if (myLuminosityBlock>LSEnd_ && LSEnd_>0) return;


  if (verbose_) std::cout << "METAnalyzer analyze" << std::endl;

  std::string DirName = "JetMET/PFCandidates/"+mInputCollection_.label();

  bool hbhenoifilterdecision=true;
  if(!isMiniAOD_){//not checked for MiniAOD -> for miniaod decision filled as "triggerResults" bool
    edm::Handle<bool> HBHENoiseFilterResultHandle;
    iEvent.getByToken(hbheNoiseFilterResultToken_, HBHENoiseFilterResultHandle);
    if (!HBHENoiseFilterResultHandle.isValid()) {
      LogDebug("") << "METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
      if (verbose_) std::cout << "METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
    }
    hbhenoifilterdecision= *HBHENoiseFilterResultHandle;
  }else{//need to check if we go for version 1 or version 2
    edm::Handle<edm::TriggerResults> metFilterResults;
    iEvent.getByToken(METFilterMiniAODToken_, metFilterResults);
    if(metFilterResults.isValid()){
      if(miniaodfilterindex!=-1){
    hbhenoifilterdecision = metFilterResults->accept(miniaodfilterindex);
      }
    }else{
      iEvent.getByToken(METFilterMiniAODToken2_, metFilterResults);
      if(metFilterResults.isValid()){
    if(miniaodfilterindex!=-1){
      hbhenoifilterdecision = metFilterResults->accept(miniaodfilterindex);
    }
      }
    }
  }

  //DCS Filter
  bool bDCSFilter = (bypassAllDCSChecks_ || DCSFilter_->filter(iEvent, iSetup));

  for (unsigned int i=0;i<countsPFCand_.size();i++) {
    countsPFCand_[i]=0;
  }
  if(bDCSFilter && hbhenoifilterdecision && bPrimaryVertex){
    if(isMiniAOD_){
      edm::Handle<std::vector<pat::PackedCandidate> > packedParticleFlow;
      iEvent.getByToken(pflowPackedToken_, packedParticleFlow);
      //11, 13, 22 for el/mu/gamma, 211 chargedHadron, 130 neutralHadrons, 1 and 2 hadHF and EGammaHF
      for (unsigned int i = 0; i < packedParticleFlow->size(); i++) {
    const pat::PackedCandidate& c = packedParticleFlow->at(i);
    for (unsigned int j=0; j<typePFCand_.size(); j++) {
      if (abs(c.pdgId())==typePFCand_[j]) {
        //second check for endcap, if inside barrel Max and Min symmetric around 0
        if ( ((c.eta()>etaMinPFCand_[j]) && (c.eta()<etaMaxPFCand_[j])) || ((c.eta()> (-etaMaxPFCand_[j])) && (c.eta()< (-etaMinPFCand_[j]))) ){
          countsPFCand_[j]+=1;
          ptPFCand_[j]   = map_of_MEs[DirName + "/"+ptPFCand_name_[j]];
          if ( ptPFCand_[j]       && ptPFCand_[j]->getRootObject()) ptPFCand_[j]->Fill(c.eta(), c.phi(), c.pt()*c.puppiWeight());
          occupancyPFCand_[j]   = map_of_MEs[DirName + "/"+occupancyPFCand_name_[j]];
          if ( occupancyPFCand_[j]       && occupancyPFCand_[j]->getRootObject()) occupancyPFCand_[j]->Fill(c.eta(), c.phi(),c.puppiWeight());
          ptPFCand_puppiNolepWeight_[j]   = map_of_MEs[DirName + "/"+ptPFCand_name_puppiNolepWeight_[j]];
          if ( ptPFCand_puppiNolepWeight_[j]       && ptPFCand_puppiNolepWeight_[j]->getRootObject()) ptPFCand_puppiNolepWeight_[j]->Fill(c.eta(), c.phi(), c.pt()*c.puppiWeightNoLep());
          occupancyPFCand_puppiNolepWeight_[j]   = map_of_MEs[DirName + "/"+occupancyPFCand_name_puppiNolepWeight_[j]];
          if ( occupancyPFCand_puppiNolepWeight_[j]       && occupancyPFCand_puppiNolepWeight_[j]->getRootObject()){ occupancyPFCand_puppiNolepWeight_[j]->Fill(c.eta(), c.phi(),c.puppiWeightNoLep());
          }
        }
      }
    }
      }
      for (unsigned int j=0; j<countsPFCand_.size(); j++) {
    multiplicityPFCand_[j]   = map_of_MEs[DirName + "/"+multiplicityPFCand_name_[j]];
    if(multiplicityPFCand_[j] && multiplicityPFCand_[j]->getRootObject()){  multiplicityPFCand_[j]->Fill(countsPFCand_[j]);
    }
      }
    }else{
      edm::Handle<std::vector<reco::PFCandidate> > particleFlow;
      iEvent.getByToken(pflowToken_, particleFlow);
       for (unsigned int i = 0; i < particleFlow->size(); i++) {
    const reco::PFCandidate& c = particleFlow->at(i);
    for (unsigned int j=0; j<typePFCandRECO_.size(); j++) {
      if (c.particleId()==typePFCandRECO_[j]) {
        //second check for endcap, if inside barrel Max and Min symmetric around 0
        if ( ((c.eta()>etaMinPFCandRECO_[j]) && (c.eta()<etaMaxPFCandRECO_[j])) || ((c.eta()> (-etaMaxPFCandRECO_[j])) && (c.eta()< (-etaMinPFCandRECO_[j]))) ){
          countsPFCandRECO_[j]+=1;
          ptPFCandRECO_[j]   = map_of_MEs[DirName + "/"+ptPFCand_nameRECO_[j]];
          if ( ptPFCandRECO_[j]       && ptPFCandRECO_[j]->getRootObject()) ptPFCandRECO_[j]->Fill(c.eta(), c.phi(), c.pt());
          occupancyPFCandRECO_[j]   = map_of_MEs[DirName + "/"+occupancyPFCand_nameRECO_[j]];
          if ( occupancyPFCandRECO_[j]       && occupancyPFCandRECO_[j]->getRootObject()) occupancyPFCandRECO_[j]->Fill(c.eta(), c.phi());
        }
        //fill quantities for isolated charged hadron quantities
        //only for charged hadrons
        if ( c.particleId() == 1 &&  c.pt() > ptMinCand_ ){
          // At least 1 GeV in HCAL
          double ecalRaw = c.rawEcalEnergy();
          double hcalRaw = c.rawHcalEnergy();
          if ( (ecalRaw + hcalRaw) > hcalMin_ ){
        const PFCandidate::ElementsInBlocks& theElements = c.elementsInBlocks();
        if( theElements.empty() ) continue;
        unsigned int iTrack=-999;
        std::vector<unsigned int> iECAL;// =999;
        std::vector<unsigned int> iHCAL;// =999;
        const reco::PFBlockRef blockRef = theElements[0].first;
        const edm::OwnVector<reco::PFBlockElement>& elements = blockRef->elements();
        // Check that there is only one track in the block.
        unsigned int nTracks = 0;
        for(unsigned int iEle=0; iEle<elements.size(); iEle++) {             
          // Find the tracks in the block
          PFBlockElement::Type type = elements[iEle].type();
          switch( type ) {
          case PFBlockElement::TRACK:
            iTrack = iEle;
            nTracks++;
            break;
          case PFBlockElement::ECAL:
            iECAL.push_back( iEle );
            break;
          case PFBlockElement::HCAL:
            iHCAL.push_back( iEle );
            break;
          default:
            continue;
          } 
        }
        if ( nTracks == 1 ){
          // Characteristics of the track
          const reco::PFBlockElementTrack& et = dynamic_cast<const reco::PFBlockElementTrack &>( elements[iTrack] );
          mProfileIsoPFChHad_TrackOccupancy=map_of_MEs[DirName+"/"+"IsoPfChHad_Track_profile"];
          if (mProfileIsoPFChHad_TrackOccupancy  && mProfileIsoPFChHad_TrackOccupancy->getRootObject()) mProfileIsoPFChHad_TrackOccupancy->Fill(et.trackRef()->eta(),et.trackRef()->phi());
          mProfileIsoPFChHad_TrackPt=map_of_MEs[DirName+"/"+"IsoPfChHad_TrackPt"];
          if (mProfileIsoPFChHad_TrackPt  && mProfileIsoPFChHad_TrackPt->getRootObject()) mProfileIsoPFChHad_TrackPt->Fill(et.trackRef()->eta(),et.trackRef()->phi(),et.trackRef()->pt());
          if(c.rawEcalEnergy()==0){//isolated hadron, nothing in ECAL
            //right now take corrected hcalEnergy, do we want the rawHcalEnergy instead
            m_HOverTrackP_trackPtVsEta=map_of_MEs[DirName+"/"+"HOverTrackP_trackPtVsEta"];
            if(m_HOverTrackP_trackPtVsEta  && m_HOverTrackP_trackPtVsEta->getRootObject()) m_HOverTrackP_trackPtVsEta->Fill(c.pt(), c.eta(), c.hcalEnergy()/et.trackRef()->p());
            if(c.pt()>1 && c.pt()<10){
              m_HOverTrackPVsEta_hPt_1_10=map_of_MEs[DirName+"/"+"HOverTrackPVsEta_hPt_1_10"];
              if(m_HOverTrackPVsEta_hPt_1_10  && m_HOverTrackPVsEta_hPt_1_10->getRootObject()) m_HOverTrackPVsEta_hPt_1_10->Fill(c.eta(), c.hcalEnergy()/et.trackRef()->p());
            }else if(c.pt()>10 && c.pt()<20){
              m_HOverTrackPVsEta_hPt_10_20=map_of_MEs[DirName+"/"+"HOverTrackPVsEta_hPt_10_20"];
              if(m_HOverTrackPVsEta_hPt_10_20  && m_HOverTrackPVsEta_hPt_10_20->getRootObject()) m_HOverTrackPVsEta_hPt_10_20->Fill(c.eta(), c.hcalEnergy()/et.trackRef()->p());
            }else if(c.pt()>20 && c.pt()<50){
              m_HOverTrackPVsEta_hPt_20_50=map_of_MEs[DirName+"/"+"HOverTrackPVsEta_hPt_20_50"];
              if(m_HOverTrackPVsEta_hPt_20_50  && m_HOverTrackPVsEta_hPt_20_50->getRootObject()) m_HOverTrackPVsEta_hPt_20_50->Fill(c.eta(), c.hcalEnergy()/et.trackRef()->p());
            }else if(c.pt()>50){
              m_HOverTrackPVsEta_hPt_50=map_of_MEs[DirName+"/"+"HOverTrackPVsEta_hPt_50"];
              if(m_HOverTrackPVsEta_hPt_50  && m_HOverTrackPVsEta_hPt_50->getRootObject()) m_HOverTrackPVsEta_hPt_50->Fill(c.eta(), c.hcalEnergy()/et.trackRef()->p());
            }
            if(fabs(c.eta()<1.392)){
              if(c.pt()>1 && c.pt()<10){
            m_HOverTrackP_Barrel_hPt_1_10=map_of_MEs[DirName+"/"+"HOverTrackP_Barrel_hPt_1_10"];
            if(m_HOverTrackP_Barrel_hPt_1_10 && m_HOverTrackP_Barrel_hPt_1_10->getRootObject()) m_HOverTrackP_Barrel_hPt_1_10->Fill(c.hcalEnergy()/et.trackRef()->p());
              }else if(c.pt()>10 && c.pt()<20){
            m_HOverTrackP_Barrel_hPt_10_20=map_of_MEs[DirName+"/"+"HOverTrackP_Barrel_hPt_10_20"];
            if(m_HOverTrackP_Barrel_hPt_10_20 && m_HOverTrackP_Barrel_hPt_10_20->getRootObject()) m_HOverTrackP_Barrel_hPt_10_20->Fill(c.hcalEnergy()/et.trackRef()->p());
              }else if(c.pt()>20 && c.pt()<50){
            m_HOverTrackP_Barrel_hPt_20_50=map_of_MEs[DirName+"/"+"HOverTrackP_Barrel_hPt_20_50"];
            if(m_HOverTrackP_Barrel_hPt_20_50 && m_HOverTrackP_Barrel_hPt_20_50->getRootObject()) m_HOverTrackP_Barrel_hPt_20_50->Fill(c.hcalEnergy()/et.trackRef()->p());
              }else if(c.pt()>50){
            m_HOverTrackP_Barrel_hPt_50=map_of_MEs[DirName+"/"+"HOverTrackP_Barrel_hPt_50"];
            if(m_HOverTrackP_Barrel_hPt_50 && m_HOverTrackP_Barrel_hPt_50->getRootObject()) m_HOverTrackP_Barrel_hPt_50->Fill(c.hcalEnergy()/et.trackRef()->p());
              }
              m_HOverTrackPVsTrackP_Barrel=map_of_MEs[DirName+"/"+"HOverTrackPVsTrackP_Barrel"];
              if(m_HOverTrackPVsTrackP_Barrel  && m_HOverTrackPVsTrackP_Barrel->getRootObject()) m_HOverTrackPVsTrackP_Barrel->Fill(et.trackRef()->p(), c.hcalEnergy()/et.trackRef()->p());
              m_HOverTrackPVsTrackPt_Barrel=map_of_MEs[DirName+"/"+"HOverTrackPVsTrackPt_Barrel"];
              if(m_HOverTrackPVsTrackPt_Barrel  && m_HOverTrackPVsTrackPt_Barrel->getRootObject()) m_HOverTrackPVsTrackPt_Barrel->Fill(et.trackRef()->pt(), c.hcalEnergy()/et.trackRef()->p());
            }else{
              m_HOverTrackPVsTrackP_EndCap=map_of_MEs[DirName+"/"+"HOverTrackPVsTrackP_EndCap"];
              if(m_HOverTrackPVsTrackP_EndCap  && m_HOverTrackPVsTrackP_EndCap->getRootObject()) m_HOverTrackPVsTrackP_EndCap->Fill(et.trackRef()->p(), c.hcalEnergy()/et.trackRef()->p());
              m_HOverTrackPVsTrackPt_EndCap=map_of_MEs[DirName+"/"+"HOverTrackPVsTrackPt_EndCap"];
              if(m_HOverTrackPVsTrackPt_EndCap  && m_HOverTrackPVsTrackPt_EndCap->getRootObject())m_HOverTrackPVsTrackPt_EndCap->Fill(et.trackRef()->pt(), c.hcalEnergy()/et.trackRef()->p());
              if(c.pt()>1 && c.pt()<10){
            m_HOverTrackP_EndCap_hPt_1_10=map_of_MEs[DirName+"/"+"HOverTrackP_EndCap_hPt_1_10"];
            if(m_HOverTrackP_EndCap_hPt_1_10 && m_HOverTrackP_EndCap_hPt_1_10->getRootObject()) m_HOverTrackP_EndCap_hPt_1_10->Fill(c.hcalEnergy()/et.trackRef()->p());
              }else if(c.pt()>10 && c.pt()<20){
            m_HOverTrackP_EndCap_hPt_10_20=map_of_MEs[DirName+"/"+"HOverTrackP_EndCap_hPt_10_20"];
            if(m_HOverTrackP_EndCap_hPt_10_20 && m_HOverTrackP_EndCap_hPt_10_20->getRootObject()) m_HOverTrackP_EndCap_hPt_10_20->Fill(c.hcalEnergy()/et.trackRef()->p());
              }else if(c.pt()>20 && c.pt()<50){
            m_HOverTrackP_EndCap_hPt_20_50=map_of_MEs[DirName+"/"+"HOverTrackP_EndCap_hPt_20_50"];
            if(m_HOverTrackP_EndCap_hPt_20_50 && m_HOverTrackP_EndCap_hPt_20_50->getRootObject()) m_HOverTrackP_EndCap_hPt_20_50->Fill(c.hcalEnergy()/et.trackRef()->p());
              }else if(c.pt()>50){
            m_HOverTrackP_EndCap_hPt_50=map_of_MEs[DirName+"/"+"HOverTrackP_EndCap_hPt_50"];
            if(m_HOverTrackP_EndCap_hPt_50 && m_HOverTrackP_EndCap_hPt_50->getRootObject()) m_HOverTrackP_EndCap_hPt_50->Fill(c.hcalEnergy()/et.trackRef()->p());
              }
            }     
          }
          //ECAL element
          for(unsigned int ii=0;ii<iECAL.size();ii++) {
            const reco::PFBlockElementCluster& eecal = dynamic_cast<const reco::PFBlockElementCluster &>( elements[ iECAL[ii] ] );
            if(fabs(eecal.clusterRef()->eta())<1.479){
              mProfileIsoPFChHad_EcalOccupancyCentral=map_of_MEs[DirName+"/"+"IsoPfChHad_ECAL_profile_central"];
              if (mProfileIsoPFChHad_EcalOccupancyCentral  && mProfileIsoPFChHad_EcalOccupancyCentral->getRootObject()) mProfileIsoPFChHad_EcalOccupancyCentral->Fill(eecal.clusterRef()->eta(),eecal.clusterRef()->phi());
              mProfileIsoPFChHad_EMPtCentral=map_of_MEs[DirName+"/"+"IsoPfChHad_EMPt_central"];
              if (mProfileIsoPFChHad_EMPtCentral  && mProfileIsoPFChHad_EMPtCentral->getRootObject()) mProfileIsoPFChHad_EMPtCentral->Fill(eecal.clusterRef()->eta(),eecal.clusterRef()->phi(),eecal.clusterRef()->pt());
            }else{
              mProfileIsoPFChHad_EcalOccupancyEndcap=map_of_MEs[DirName+"/"+"IsoPfChHad_ECAL_profile_endcap"];
              if (mProfileIsoPFChHad_EcalOccupancyEndcap  && mProfileIsoPFChHad_EcalOccupancyEndcap->getRootObject()) mProfileIsoPFChHad_EcalOccupancyEndcap->Fill(eecal.clusterRef()->eta(),eecal.clusterRef()->phi());
              mProfileIsoPFChHad_EMPtEndcap=map_of_MEs[DirName+"/"+"IsoPfChHad_EMPt_endcap"];
              if (mProfileIsoPFChHad_EMPtEndcap  && mProfileIsoPFChHad_EMPtEndcap->getRootObject()) mProfileIsoPFChHad_EMPtEndcap->Fill(eecal.clusterRef()->eta(),eecal.clusterRef()->phi(),eecal.clusterRef()->pt());
            }


          } 
          //HCAL element
          for(unsigned int ii=0;ii<iHCAL.size();ii++) {
                const reco::PFBlockElementCluster& ehcal = dynamic_cast<const reco::PFBlockElementCluster &>( elements[ iHCAL[ii] ] );
            if(fabs(ehcal.clusterRef()->eta())<1.740){
              mProfileIsoPFChHad_HcalOccupancyCentral=map_of_MEs[DirName+"/"+"IsoPfChHad_HCAL_profile_central"];
              if (mProfileIsoPFChHad_HcalOccupancyCentral  && mProfileIsoPFChHad_HcalOccupancyCentral->getRootObject()) mProfileIsoPFChHad_HcalOccupancyCentral->Fill(ehcal.clusterRef()->eta(),ehcal.clusterRef()->phi());
              mProfileIsoPFChHad_HadPtCentral=map_of_MEs[DirName+"/"+"IsoPfChHad_HadPt_central"];
              if (mProfileIsoPFChHad_HadPtCentral  && mProfileIsoPFChHad_HadPtCentral->getRootObject()) mProfileIsoPFChHad_HadPtCentral->Fill(ehcal.clusterRef()->eta(),ehcal.clusterRef()->phi(),ehcal.clusterRef()->pt());
            }else{
              mProfileIsoPFChHad_HcalOccupancyEndcap=map_of_MEs[DirName+"/"+"IsoPfChHad_HCAL_profile_endcap"];
              if (mProfileIsoPFChHad_HcalOccupancyEndcap  && mProfileIsoPFChHad_HcalOccupancyEndcap->getRootObject()) mProfileIsoPFChHad_HcalOccupancyEndcap->Fill(ehcal.clusterRef()->eta(),ehcal.clusterRef()->phi());
              mProfileIsoPFChHad_HadPtEndcap=map_of_MEs[DirName+"/"+"IsoPfChHad_HadPt_endcap"];
              if (mProfileIsoPFChHad_HadPtEndcap  && mProfileIsoPFChHad_HadPtEndcap->getRootObject()) mProfileIsoPFChHad_HadPtEndcap->Fill(ehcal.clusterRef()->eta(),ehcal.clusterRef()->phi(),ehcal.clusterRef()->pt());
            }   
          }
        }
          }
        }
      } 
    }
      }
      for (unsigned int j=0; j<countsPFCandRECO_.size(); j++) {
    multiplicityPFCandRECO_[j]   = map_of_MEs[DirName + "/"+multiplicityPFCand_nameRECO_[j]];
    if(multiplicityPFCandRECO_[j] && multiplicityPFCandRECO_[j]->getRootObject()){  multiplicityPFCandRECO_[j]->Fill(countsPFCandRECO_[j]);
    }
      }
    }//candidate loop for both miniaod and reco
  }
}

void DQMPFCandidateAnalyzer::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  iRun, edm::EventSetup const &    )

override

Inizialize parameters for histo binning.

Definition at line 132 of file DQMPFCandidateAnalyzer.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), cmsPerfPublish::DirName, JetChargeProducer_cfi::exp, objects.autophobj::float, mps_fire::i, cmsBatch::log, M_PI, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and findQualityFiles::v.

                                             {

  ibooker.setCurrentFolder("JetMET/PFCandidates/"+mInputCollection_.label());
  std::string DirName = "JetMET/PFCandidates/"+mInputCollection_.label();

  if(!isMiniAOD_){
    if(!occupancyPFCandRECO_.empty())occupancyPFCandRECO_.clear();
    if(!occupancyPFCand_nameRECO_.empty())occupancyPFCand_nameRECO_.clear();
    if(!etaMinPFCandRECO_.empty())etaMinPFCandRECO_.clear();
    if(!etaMaxPFCandRECO_.empty())etaMaxPFCandRECO_.clear();
    if(!typePFCandRECO_.empty())typePFCandRECO_.clear();
    if(!countsPFCandRECO_.empty())countsPFCandRECO_.clear();  
    if(!ptPFCandRECO_.empty())ptPFCandRECO_.clear();
    if(!ptPFCand_nameRECO_.empty())ptPFCand_nameRECO_.clear();
    if(!multiplicityPFCandRECO_.empty())multiplicityPFCandRECO_.clear();
    if(!multiplicityPFCand_nameRECO_.empty())multiplicityPFCand_nameRECO_.clear();
    for (std::vector<edm::ParameterSet>::const_iterator v = diagnosticsParameters_.begin(); v!=diagnosticsParameters_.end(); v++) {
      int etaNBinsPFCand = v->getParameter<int>("etaNBins");
      double etaMinPFCand = v->getParameter<double>("etaMin");
      double etaMaxPFCand = v->getParameter<double>("etaMax");
      int phiNBinsPFCand = v->getParameter<int>("phiNBins");
      double phiMinPFCand = v->getParameter<double>("phiMin");
      double phiMaxPFCand = v->getParameter<double>("phiMax");
      int nMinPFCand = v->getParameter<int>("nMin");
      int nMaxPFCand = v->getParameter<int>("nMax");
      int nbinsPFCand = v->getParameter<double>("nbins");      
      etaMinPFCandRECO_.push_back(etaMinPFCand);
      etaMaxPFCandRECO_.push_back(etaMaxPFCand);
      typePFCandRECO_.push_back(v->getParameter<int>("type"));
      countsPFCandRECO_.push_back(0);
      multiplicityPFCandRECO_.push_back(ibooker.book1D(std::string(v->getParameter<std::string>("name")).append("_multiplicity_").c_str(),std::string(v->getParameter<std::string>("name"))+"multiplicity", nbinsPFCand, nMinPFCand, nMaxPFCand));
      multiplicityPFCand_nameRECO_.push_back(std::string(v->getParameter<std::string>("name")).append("_multiplicity_"));
      map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+ multiplicityPFCand_nameRECO_[multiplicityPFCand_nameRECO_.size()-1], multiplicityPFCandRECO_[multiplicityPFCandRECO_.size()-1]));
      
      //push back names first, we need to create histograms with the name and fill it for endcap plots later
      occupancyPFCand_nameRECO_.push_back(std::string(v->getParameter<std::string>("name")).append("_occupancy_"));
      
      ptPFCand_nameRECO_.push_back(std::string(v->getParameter<std::string>("name")).append("_pt_"));
      //special booking for endcap plots, merge plots for eminus and eplus into one plot, using variable binning
      //barrel plots have eta-boundaries on minus and plus side
      //parameters start on minus side
      if(etaMinPFCand*etaMaxPFCand<0){//barrel plots, plot only in barrel region
    occupancyPFCandRECO_.push_back(ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_occupancy_").c_str(),std::string(v->getParameter<std::string>("name"))+"occupancy", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
    ptPFCandRECO_.push_back(ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_pt_").c_str(),std::string(v->getParameter<std::string>("name"))+"pt", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
      }else{//endcap or forward plots, 
    const int nbins_eta_endcap=2*(etaNBinsPFCand+1);
    double eta_limits_endcap[nbins_eta_endcap];
    for(int i=0;i<nbins_eta_endcap;i++){
      if(i<(etaNBinsPFCand+1)){
        eta_limits_endcap[i]=etaMinPFCand+i*(etaMaxPFCand-etaMinPFCand)/(double)etaNBinsPFCand;
      }else{
        eta_limits_endcap[i]= -etaMaxPFCand +(i- (etaNBinsPFCand+1) )*(etaMaxPFCand-etaMinPFCand)/(double)etaNBinsPFCand;
      }
    }
    TH2F* hist_temp_occup = new TH2F((occupancyPFCand_nameRECO_[occupancyPFCand_nameRECO_.size()-1]).c_str(),"occupancy",nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
    occupancyPFCandRECO_.push_back(ibooker.book2D(occupancyPFCand_nameRECO_[occupancyPFCand_nameRECO_.size()-1],hist_temp_occup));
    TH2F* hist_temp_pt = new TH2F((ptPFCand_nameRECO_[ptPFCand_nameRECO_.size()-1]).c_str(),"pt",nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
    ptPFCandRECO_.push_back(ibooker.book2D(ptPFCand_nameRECO_[ptPFCand_nameRECO_.size()-1], hist_temp_pt));
      }
      
      map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+ occupancyPFCand_nameRECO_[occupancyPFCand_nameRECO_.size()-1], occupancyPFCandRECO_[occupancyPFCandRECO_.size()-1]));
      map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+ ptPFCand_nameRECO_[ptPFCand_nameRECO_.size()-1], ptPFCandRECO_[ptPFCandRECO_.size()-1]));
    }
    
    mProfileIsoPFChHad_TrackOccupancy=ibooker.book2D("IsoPfChHad_Track_profile","Isolated PFChHadron Tracker_occupancy", 108, -2.7, 2.7, 160, -M_PI,M_PI);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_Track_profile"        ,mProfileIsoPFChHad_TrackOccupancy));
    mProfileIsoPFChHad_TrackPt=ibooker.book2D("IsoPfChHad_TrackPt","Isolated PFChHadron TrackPt", 108, -2.7, 2.7, 160, -M_PI,M_PI);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_TrackPt"        ,mProfileIsoPFChHad_TrackPt));
    
    
    mProfileIsoPFChHad_EcalOccupancyCentral  = ibooker.book2D("IsoPfChHad_ECAL_profile_central","IsolatedPFChHa ECAL occupancy (Barrel)", 180, -1.479, 1.479, 125, -M_PI,M_PI);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_ECAL_profile_central"        ,mProfileIsoPFChHad_EcalOccupancyCentral));
    mProfileIsoPFChHad_EMPtCentral=ibooker.book2D("IsoPfChHad_EMPt_central","Isolated PFChHadron HadPt_central", 180, -1.479, 1.479, 360, -M_PI,M_PI);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_EMPt_central"        ,mProfileIsoPFChHad_EMPtCentral));
    
    mProfileIsoPFChHad_EcalOccupancyEndcap  = ibooker.book2D("IsoPfChHad_ECAL_profile_endcap","IsolatedPFChHa ECAL occupancy (Endcap)", 110, -2.75, 2.75, 125, -M_PI,M_PI);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_ECAL_profile_endcap"        ,mProfileIsoPFChHad_EcalOccupancyEndcap));
    mProfileIsoPFChHad_EMPtEndcap=ibooker.book2D("IsoPfChHad_EMPt_endcap","Isolated PFChHadron EMPt_endcap", 110, -2.75, 2.75, 125, -M_PI,M_PI);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_EMPt_endcap"        ,mProfileIsoPFChHad_EMPtEndcap));
    
    const int nbins_eta=16;
    
    double eta_limits[nbins_eta]={-2.650,-2.500,-2.322,-2.172,-2.043,-1.930,-1.830,-1.740,1.740,1.830,1.930,2.043,2.172,2.3122,2.500,2.650};
    
    TH2F* hist_temp_HCAL =new TH2F("IsoPfChHad_HCAL_profile_endcap","IsolatedPFChHa HCAL occupancy (outer endcap)",nbins_eta-1,eta_limits, 36, -M_PI,M_PI);
    TH2F* hist_tempPt_HCAL=(TH2F*)hist_temp_HCAL->Clone("Isolated PFCHHadron HadPt (outer endcap)");
    
    mProfileIsoPFChHad_HcalOccupancyCentral  = ibooker.book2D("IsoPfChHad_HCAL_profile_central","IsolatedPFChHa HCAL occupancy (Central Part)", 40, -1.740, 1.740, 72, -M_PI,M_PI);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_HCAL_profile_central"        ,mProfileIsoPFChHad_HcalOccupancyCentral));
    mProfileIsoPFChHad_HadPtCentral=ibooker.book2D("IsoPfChHad_HadPt_central","Isolated PFChHadron HadPt_central", 40, -1.740, 1.740, 72, -M_PI,M_PI);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_HadPt_central"        ,mProfileIsoPFChHad_HadPtCentral));
    
    mProfileIsoPFChHad_HcalOccupancyEndcap  = ibooker.book2D("IsoPfChHad_HCAL_profile_endcap",hist_temp_HCAL);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_HCAL_profile_endcap"        ,mProfileIsoPFChHad_HcalOccupancyEndcap));
    mProfileIsoPFChHad_HadPtEndcap=ibooker.book2D("IsoPfChHad_HadPt_endcap",hist_tempPt_HCAL);
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"IsoPfChHad_HadPt_endcap"        ,mProfileIsoPFChHad_HadPtEndcap));

    //actual HCAL segmentation in pseudorapidity -> reduce by a factor of two
    //const int nbins_eta_hcal_total=54;   
    //double eta_limits_hcal_total[nbins_eta_hcal_total]=
    //   {-2.650,-2.500,-2.322,-2.172,-2.043,-1.930,-1.830,-1.740,-1.653,-1.566,-1.479,-1.392,-1.305, 
    //   -1.218,-1.131,-1.044,-0.957,-0.870,-0.783,-0.696,-0.609,-0.522,-0.435,-0.348,-0.261,-0.174,-0.087,0.0,
    //   0.087, 0.174, 0.261, 0.348, 0.435, 0.522, 0.609, 0.783, 0.870, 0.957, 1.044, 1.131, 1.218
    //   1.305, 1.392, 1.479, 1.566, 1.653, 1.740, 1.830, 1.930, 2.043, 2.172, 2.322, 2.500, 2.650}
    //
    
    const int nbins_eta_hcal_total=28;   
    double eta_limits_hcal_total[nbins_eta_hcal_total]={-2.650,-2.322,-2.043,-1.830,-1.653,-1.479,-1.305,-1.131,-0.957,-0.783,-0.609,-0.435,-0.261,-0.087,
                                0.087, 0.261, 0.435, 0.609, 0.783, 0.957, 1.131, 1.305, 1.479, 1.653, 1.830, 2.043, 2.322, 2.650};
    float eta_limits_hcal_total_f[nbins_eta_hcal_total];
    float log_bin_spacing = log(200.)/40.;
    const int nbins_pt_total_hcal= 41;
    double pt_limits_hcal[nbins_pt_total_hcal];
    float pt_limits_hcal_f[nbins_pt_total_hcal];
    for(int i=0;i<nbins_pt_total_hcal;i++){
      pt_limits_hcal[i]=exp(i*log_bin_spacing);
      pt_limits_hcal_f[i]=exp(i*log_bin_spacing);
    }
    for(int i=0;i<nbins_eta_hcal_total;i++){
      eta_limits_hcal_total_f[i]=(float)eta_limits_hcal_total[i];
    }
    m_HOverTrackP_trackPtVsEta= ibooker.book2D("HOverTrackP_trackPtVsEta","HOverTrackP_trackPtVsEta",nbins_pt_total_hcal-1,pt_limits_hcal_f,nbins_eta_hcal_total-1,eta_limits_hcal_total_f);
    m_HOverTrackPVsTrackP_Barrel = ibooker.bookProfile("HOverTrackPVsTrackP_Barrel","HOverTrackPVsTrackP_Barrel",nbins_pt_total_hcal-1,pt_limits_hcal, 0, 4, " ");
    m_HOverTrackPVsTrackP_EndCap = ibooker.bookProfile("HOverTrackPVsTrackP_EndCap","HOverTrackPVsTrackP_EndCap",nbins_pt_total_hcal-1,pt_limits_hcal, 0, 4, " ");
    m_HOverTrackPVsTrackPt_Barrel = ibooker.bookProfile("HOverTrackPVsTrackPt_Barrel","HOverTrackPVsTrackPt_Barrel",nbins_pt_total_hcal-1,pt_limits_hcal, 0, 4, " ");
    m_HOverTrackPVsTrackPt_EndCap = ibooker.bookProfile("HOverTrackPVsTrackPt_EndCap","HOverTrackPVsTrackPt_EndCap",nbins_pt_total_hcal-1,pt_limits_hcal, 0, 4, " ");
 
    m_HOverTrackPVsEta_hPt_1_10 = ibooker.bookProfile("HOverTrackPVsEta_hPt_1_10","HOverTrackPVsEta, 1<hPt<10 GeV",nbins_eta_hcal_total-1, eta_limits_hcal_total, 0, 4, " ");
    m_HOverTrackPVsEta_hPt_10_20 = ibooker.bookProfile("HOverTrackPVsEta_hPt_10_20","HOverTrackPVsEta, 10<hPt<20 GeV",nbins_eta_hcal_total-1, eta_limits_hcal_total, 0, 4, " ");
    m_HOverTrackPVsEta_hPt_20_50 = ibooker.bookProfile("HOverTrackPVsEta_hPt_20_50","HOverTrackPVsEta, 20<hPt<50 GeV",nbins_eta_hcal_total-1, eta_limits_hcal_total, 0, 4, " ");
    m_HOverTrackPVsEta_hPt_50 = ibooker.bookProfile("HOverTrackPVsEta_hPt_50","HOverTrackPVsEta, hPt>50 GeV",nbins_eta_hcal_total-1, eta_limits_hcal_total, 0, 4, " ");

    m_HOverTrackP_Barrel_hPt_1_10= ibooker.book1D("HOverTrackP_Barrel_hPt_1_10","HOverTrackP_B, 1<hPt<10 GeV",50,0,4);
    m_HOverTrackP_Barrel_hPt_10_20= ibooker.book1D("HOverTrackP_Barrel_hPt_10_20","HOverTrackP_B, 10<hPt<20 GeV",50,0,4);
    m_HOverTrackP_Barrel_hPt_20_50= ibooker.book1D("HOverTrackP_Barrel_hPt_20_50","HOverTrackP_B, 20<hPt<50 GeV",50,0,4);
    m_HOverTrackP_Barrel_hPt_50= ibooker.book1D("HOverTrackP_Barrel_hPt_50","HOverTrackP_B, hPt>50 GeV",50,0,4);

    m_HOverTrackP_EndCap_hPt_1_10= ibooker.book1D("HOverTrackP_EndCap_hPt_1_10","HOverTrackP_E, 1<hPt<10 GeV",50,0,4);
    m_HOverTrackP_EndCap_hPt_10_20= ibooker.book1D("HOverTrackP_EndCap_hPt_10_20","HOverTrackP_E, 10<hPt<20 GeV",50,0,4);
    m_HOverTrackP_EndCap_hPt_20_50= ibooker.book1D("HOverTrackP_EndCap_hPt_20_50","HOverTrackP_E, 20<hPt<50 GeV",50,0,4);
    m_HOverTrackP_EndCap_hPt_50= ibooker.book1D("HOverTrackP_EndCap_hPt_50","HOverTrackP_E, hPt>50 GeV",50,0,4);

    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackP_trackPtVsEta"        ,m_HOverTrackP_trackPtVsEta));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackPVsTrackP_Barrel"        ,m_HOverTrackPVsTrackP_Barrel));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackPVsTrackP_EndCap"        ,m_HOverTrackPVsTrackP_EndCap));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackPVsTrackPt_Barrel"        ,m_HOverTrackPVsTrackPt_Barrel));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackPVsTrackPt_EndCap"        ,m_HOverTrackPVsTrackPt_EndCap));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackPVsEta_hPt_1_10"        ,m_HOverTrackPVsEta_hPt_1_10));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackPVsEta_hPt_10_20"        ,m_HOverTrackPVsEta_hPt_10_20));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackPVsEta_hPt_20_50"        ,m_HOverTrackPVsEta_hPt_20_50));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackPVsEta_hPt_50"        ,m_HOverTrackPVsEta_hPt_50));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackP_Barrel_hPt_1_10"        ,m_HOverTrackP_Barrel_hPt_1_10));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackP_Barrel_hPt_10_20"        ,m_HOverTrackP_Barrel_hPt_10_20));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackP_Barrel_hPt_20_50"        ,m_HOverTrackP_Barrel_hPt_20_50));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackP_Barrel_hPt_50"        ,m_HOverTrackP_Barrel_hPt_50));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackP_EndCap_hPt_1_10"        ,m_HOverTrackP_EndCap_hPt_1_10));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackP_EndCap_hPt_10_20"        ,m_HOverTrackP_EndCap_hPt_10_20));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackP_EndCap_hPt_20_50"        ,m_HOverTrackP_EndCap_hPt_20_50));
    map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+"HOverTrackP_EndCap_hPt_50"        ,m_HOverTrackP_EndCap_hPt_50));
  }else{//MiniAOD workflow
    if(!occupancyPFCand_.empty())occupancyPFCand_.clear();
    if(!occupancyPFCand_name_.empty())occupancyPFCand_name_.clear();
    if(!occupancyPFCand_puppiNolepWeight_.empty())occupancyPFCand_puppiNolepWeight_.clear();
    if(!occupancyPFCand_name_puppiNolepWeight_.empty())occupancyPFCand_name_puppiNolepWeight_.clear();
    if(!etaMinPFCand_.empty())etaMinPFCand_.clear();
    if(!etaMaxPFCand_.empty())etaMaxPFCand_.clear();
    if(!typePFCand_.empty())typePFCand_.clear();
    if(!countsPFCand_.empty())countsPFCand_.clear();  
    if(!ptPFCand_.empty())ptPFCand_.clear();
    if(!ptPFCand_name_.empty())ptPFCand_name_.clear();
    if(!ptPFCand_puppiNolepWeight_.empty())ptPFCand_puppiNolepWeight_.clear();
    if(!ptPFCand_name_puppiNolepWeight_.empty())ptPFCand_name_puppiNolepWeight_.clear();
    if(!multiplicityPFCand_.empty())multiplicityPFCand_.clear();
    if(!multiplicityPFCand_name_.empty())multiplicityPFCand_name_.clear();
    for (std::vector<edm::ParameterSet>::const_iterator v = diagnosticsParameters_.begin(); v!=diagnosticsParameters_.end(); v++) {
      int etaNBinsPFCand = v->getParameter<int>("etaNBins");
      double etaMinPFCand = v->getParameter<double>("etaMin");
      double etaMaxPFCand = v->getParameter<double>("etaMax");
      int phiNBinsPFCand = v->getParameter<int>("phiNBins");
      double phiMinPFCand = v->getParameter<double>("phiMin");
      double phiMaxPFCand = v->getParameter<double>("phiMax");
      int nMinPFCand = v->getParameter<int>("nMin");
      int nMaxPFCand = v->getParameter<int>("nMax");
      int nbinsPFCand = v->getParameter<double>("nbins");
      
      // etaNBins_.push_back(etaNBins);
      etaMinPFCand_.push_back(etaMinPFCand);
      etaMaxPFCand_.push_back(etaMaxPFCand);
      typePFCand_.push_back(v->getParameter<int>("type"));
      countsPFCand_.push_back(0);
      
      multiplicityPFCand_.push_back(ibooker.book1D(std::string(v->getParameter<std::string>("name")).append("_multiplicity_").c_str(),std::string(v->getParameter<std::string>("name"))+"multiplicity", nbinsPFCand, nMinPFCand, nMaxPFCand));
      multiplicityPFCand_name_.push_back(std::string(v->getParameter<std::string>("name")).append("_multiplicity_"));
      map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+ multiplicityPFCand_name_[multiplicityPFCand_name_.size()-1], multiplicityPFCand_[multiplicityPFCand_.size()-1]));
      //push back names first, we need to create histograms with the name and fill it for endcap plots later
      occupancyPFCand_name_.push_back(std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiWeight_"));
      ptPFCand_name_.push_back(std::string(v->getParameter<std::string>("name")).append("_pt_puppiWeight_"));
      //push back names first, we need to create histograms with the name and fill it for endcap plots later
      occupancyPFCand_name_puppiNolepWeight_.push_back(std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiNolepWeight_"));     
      ptPFCand_name_puppiNolepWeight_.push_back(std::string(v->getParameter<std::string>("name")).append("_pt_puppiNolepWeight_"));
      //special booking for endcap plots, merge plots for eminus and eplus into one plot, using variable binning
      //barrel plots have eta-boundaries on minus and plus side
      //parameters start on minus side
      if(etaMinPFCand*etaMaxPFCand<0){//barrel plots, plot only in barrel region
    occupancyPFCand_.push_back(ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiWeight_").c_str(),std::string(v->getParameter<std::string>("name"))+"occupancy_puppiWeight_", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
    ptPFCand_.push_back(ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_pt_puppiWeight_").c_str(),std::string(v->getParameter<std::string>("name"))+"pt_puppiWeight_", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
    occupancyPFCand_puppiNolepWeight_.push_back(ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiNolepWeight_").c_str(),std::string(v->getParameter<std::string>("name"))+"occupancy_puppiNolepWeight", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
    ptPFCand_puppiNolepWeight_.push_back(ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_pt_puppiNolepWeight_").c_str(),std::string(v->getParameter<std::string>("name"))+"pt_puppiNolepWeight", etaNBinsPFCand, etaMinPFCand, etaMaxPFCand, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand));
      }else{//endcap or forward plots, 
    const int nbins_eta_endcap=2*(etaNBinsPFCand+1);
    double eta_limits_endcap[nbins_eta_endcap];
    for(int i=0;i<nbins_eta_endcap;i++){
      if(i<(etaNBinsPFCand+1)){
        eta_limits_endcap[i]=etaMinPFCand+i*(etaMaxPFCand-etaMinPFCand)/(double)etaNBinsPFCand;
      }else{
        eta_limits_endcap[i]= -etaMaxPFCand +(i- (etaNBinsPFCand+1) )*(etaMaxPFCand-etaMinPFCand)/(double)etaNBinsPFCand;
      }
    }
    TH2F* hist_temp_occup = new TH2F((occupancyPFCand_name_[occupancyPFCand_name_.size()-1]).c_str(),(occupancyPFCand_name_[occupancyPFCand_name_.size()-1]).c_str(),nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
    occupancyPFCand_.push_back(ibooker.book2D(occupancyPFCand_name_[occupancyPFCand_name_.size()-1],hist_temp_occup));
    TH2F* hist_temp_pt = new TH2F((ptPFCand_name_[ptPFCand_name_.size()-1]).c_str(),(ptPFCand_name_[ptPFCand_name_.size()-1]).c_str(),nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
    ptPFCand_.push_back(ibooker.book2D(ptPFCand_name_[ptPFCand_name_.size()-1], hist_temp_pt));
    TH2F* hist_temp_occup_puppiNolepWeight = new TH2F((occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size()-1]).c_str(),(occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size()-1]).c_str(),nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
    occupancyPFCand_puppiNolepWeight_.push_back(ibooker.book2D(occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size()-1],hist_temp_occup_puppiNolepWeight));
    TH2F* hist_temp_pt_puppiNolepWeight = new TH2F((ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size()-1]).c_str(),(ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size()-1]).c_str(),nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
    ptPFCand_puppiNolepWeight_.push_back(ibooker.book2D(ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size()-1], hist_temp_pt_puppiNolepWeight));
      }
      map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+ occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size()-1], occupancyPFCand_puppiNolepWeight_[occupancyPFCand_puppiNolepWeight_.size()-1]));
      map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+ ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size()-1], ptPFCand_puppiNolepWeight_[ptPFCand_puppiNolepWeight_.size()-1]));
      map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+ occupancyPFCand_name_[occupancyPFCand_name_.size()-1], occupancyPFCand_[occupancyPFCand_.size()-1]));
      map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+"/"+ ptPFCand_name_[ptPFCand_name_.size()-1], ptPFCand_[ptPFCand_.size()-1]));
    }
  }   


}   

void DQMPFCandidateAnalyzer::dqmBeginRun ( const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

override

Initialize run-based parameters.

Definition at line 372 of file DQMPFCandidateAnalyzer.cc.

References runEdmFileComparison::found, mps_fire::i, cond::persistency::search(), electronDataDiscovery::search2(), and AlCaHLTBitMon_QueryRunRegistry::string.

{

  miniaodfilterindex=-1;

  if(isMiniAOD_){
    bool changed_filter=true;
    if (FilterhltConfig_.init(iRun,iSetup,METFilterMiniAODLabel_.process(),changed_filter)){
      miniaodfilterdec=0;
      for(unsigned int i=0;i<FilterhltConfig_.size();i++){
    std::string search=FilterhltConfig_.triggerName(i).substr(5);//actual label of filter, the first 5 items are Flag_, so stripped off
    std::string search2=HBHENoiseStringMiniAOD;//all filters end with DQM, which is not in the flag --> ONLY not for HBHEFilters
    std::size_t found=search2.find(search);
    if(found!=std::string::npos){
      miniaodfilterindex=i;
    }
      }
    }else if(FilterhltConfig_.init(iRun,iSetup,METFilterMiniAODLabel2_.process(),changed_filter)){
      miniaodfilterdec=1;
      for(unsigned int i=0;i<FilterhltConfig_.size();i++){
    std::string search=FilterhltConfig_.triggerName(i).substr(5);//actual label of filter, the first 5 items are Flag_, so stripped off
    std::string search2=HBHENoiseStringMiniAOD;//all filters end with DQM, which is not in the flag --> ONLY not for HBHEFilters
    std::size_t found=search2.find(search);
    if(found!=std::string::npos){
      miniaodfilterindex=i;
    }
      }
    }else{
      edm::LogWarning("MiniAOD MET Filter HLT OBject version")<<"nothing found with both RECO and reRECO label"<<std::endl;
    }
  }
}

void DQMPFCandidateAnalyzer::endRun ( const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

override

Finish up a run.

Definition at line 406 of file DQMPFCandidateAnalyzer.cc.

{
}

static bool DQMPFCandidateAnalyzer::jetSortingRule ( reco::Jet  x, reco::Jet  y  )

inlinestaticprivate

Definition at line 91 of file DQMPFCandidateAnalyzer.h.

References reco::LeafCandidate::pt().

{return x.pt() > y.pt();}

Member Data Documentation

bool DQMPFCandidateAnalyzer::bypassAllDCSChecks_

private

Definition at line 187 of file DQMPFCandidateAnalyzer.h.

bool DQMPFCandidateAnalyzer::bypassAllPVChecks_

private

Definition at line 186 of file DQMPFCandidateAnalyzer.h.

std::string DQMPFCandidateAnalyzer::candidateType_

private

Definition at line 118 of file DQMPFCandidateAnalyzer.h.

edm::ParameterSet DQMPFCandidateAnalyzer::cleaningParameters_

private

Definition at line 125 of file DQMPFCandidateAnalyzer.h.

std::vector<int> DQMPFCandidateAnalyzer::countsPFCand_

private

Definition at line 172 of file DQMPFCandidateAnalyzer.h.

std::vector<int> DQMPFCandidateAnalyzer::countsPFCandRECO_

private

Definition at line 178 of file DQMPFCandidateAnalyzer.h.

JetMETDQMDCSFilter* DQMPFCandidateAnalyzer::DCSFilter_

private

Definition at line 123 of file DQMPFCandidateAnalyzer.h.

std::vector<edm::ParameterSet> DQMPFCandidateAnalyzer::diagnosticsParameters_

private

Definition at line 126 of file DQMPFCandidateAnalyzer.h.

std::vector<double> DQMPFCandidateAnalyzer::etaMaxPFCand_

private

Definition at line 171 of file DQMPFCandidateAnalyzer.h.

std::vector<double> DQMPFCandidateAnalyzer::etaMaxPFCandRECO_

private

Definition at line 177 of file DQMPFCandidateAnalyzer.h.

std::vector<double> DQMPFCandidateAnalyzer::etaMinPFCand_

private

Definition at line 171 of file DQMPFCandidateAnalyzer.h.

std::vector<double> DQMPFCandidateAnalyzer::etaMinPFCandRECO_

private

Definition at line 177 of file DQMPFCandidateAnalyzer.h.

HLTConfigProvider DQMPFCandidateAnalyzer::FilterhltConfig_

private

Definition at line 104 of file DQMPFCandidateAnalyzer.h.

bool DQMPFCandidateAnalyzer::hbhenoifilterdecision

private

Definition at line 110 of file DQMPFCandidateAnalyzer.h.

edm::InputTag DQMPFCandidateAnalyzer::hbheNoiseFilterResultTag_

private

Definition at line 100 of file DQMPFCandidateAnalyzer.h.

edm::EDGetTokenT<bool> DQMPFCandidateAnalyzer::hbheNoiseFilterResultToken_

private

Definition at line 101 of file DQMPFCandidateAnalyzer.h.

std::string DQMPFCandidateAnalyzer::HBHENoiseStringMiniAOD

private

Definition at line 102 of file DQMPFCandidateAnalyzer.h.

double DQMPFCandidateAnalyzer::hcalMin_

private

Definition at line 130 of file DQMPFCandidateAnalyzer.h.

bool DQMPFCandidateAnalyzer::isMiniAO_

private

Definition at line 120 of file DQMPFCandidateAnalyzer.h.

bool DQMPFCandidateAnalyzer::isMiniAOD_

private

Definition at line 191 of file DQMPFCandidateAnalyzer.h.

int DQMPFCandidateAnalyzer::LSBegin_

private

Definition at line 183 of file DQMPFCandidateAnalyzer.h.

int DQMPFCandidateAnalyzer::LSEnd_

private

Definition at line 184 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_10_20

private

Definition at line 145 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_1_10

private

Definition at line 144 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_20_50

private

Definition at line 146 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_50

private

Definition at line 147 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_10_20

private

Definition at line 150 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_1_10

private

Definition at line 149 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_20_50

private

Definition at line 151 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_50

private

Definition at line 152 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_trackPtVsEta

private

Definition at line 132 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_10_20

private

Definition at line 140 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_1_10

private

Definition at line 139 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_20_50

private

Definition at line 141 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_50

private

Definition at line 142 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackP_Barrel

private

Definition at line 133 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackP_EndCap

private

Definition at line 134 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackPt_Barrel

private

Definition at line 136 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackPt_EndCap

private

Definition at line 137 of file DQMPFCandidateAnalyzer.h.

std::map< std::string,MonitorElement* > DQMPFCandidateAnalyzer::map_of_MEs

private

Definition at line 189 of file DQMPFCandidateAnalyzer.h.

edm::InputTag DQMPFCandidateAnalyzer::METFilterMiniAODLabel2_

private

Definition at line 107 of file DQMPFCandidateAnalyzer.h.

edm::InputTag DQMPFCandidateAnalyzer::METFilterMiniAODLabel_

private

Definition at line 105 of file DQMPFCandidateAnalyzer.h.

edm::EDGetTokenT<edm::TriggerResults> DQMPFCandidateAnalyzer::METFilterMiniAODToken2_

private

Definition at line 108 of file DQMPFCandidateAnalyzer.h.

edm::EDGetTokenT<edm::TriggerResults> DQMPFCandidateAnalyzer::METFilterMiniAODToken_

private

Definition at line 106 of file DQMPFCandidateAnalyzer.h.

int DQMPFCandidateAnalyzer::miniaodfilterdec

private

Definition at line 113 of file DQMPFCandidateAnalyzer.h.

int DQMPFCandidateAnalyzer::miniaodfilterindex

private

Definition at line 111 of file DQMPFCandidateAnalyzer.h.

edm::InputTag DQMPFCandidateAnalyzer::mInputCollection_

private

Definition at line 98 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EcalOccupancyCentral

private

Definition at line 162 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EcalOccupancyEndcap

private

Definition at line 163 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EMPtCentral

private

Definition at line 156 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EMPtEndcap

private

Definition at line 157 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HadPtCentral

private

Definition at line 154 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HadPtEndcap

private

Definition at line 155 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HcalOccupancyCentral

private

Definition at line 160 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HcalOccupancyEndcap

private

Definition at line 161 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_TrackOccupancy

private

Definition at line 164 of file DQMPFCandidateAnalyzer.h.

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_TrackPt

private

Definition at line 158 of file DQMPFCandidateAnalyzer.h.

std::vector<MonitorElement* > DQMPFCandidateAnalyzer::multiplicityPFCand_

private

Definition at line 167 of file DQMPFCandidateAnalyzer.h.

std::vector<std::string> DQMPFCandidateAnalyzer::multiplicityPFCand_name_

private

Definition at line 168 of file DQMPFCandidateAnalyzer.h.

std::vector<std::string> DQMPFCandidateAnalyzer::multiplicityPFCand_nameRECO_

private

Definition at line 176 of file DQMPFCandidateAnalyzer.h.

std::vector<MonitorElement* > DQMPFCandidateAnalyzer::multiplicityPFCandRECO_

private

Definition at line 175 of file DQMPFCandidateAnalyzer.h.

int DQMPFCandidateAnalyzer::numPV_

private

Definition at line 180 of file DQMPFCandidateAnalyzer.h.

std::vector<MonitorElement* > DQMPFCandidateAnalyzer::occupancyPFCand_

private

Definition at line 167 of file DQMPFCandidateAnalyzer.h.

std::vector<std::string> DQMPFCandidateAnalyzer::occupancyPFCand_name_

private

Definition at line 168 of file DQMPFCandidateAnalyzer.h.

std::vector<std::string> DQMPFCandidateAnalyzer::occupancyPFCand_name_puppiNolepWeight_

private

Definition at line 170 of file DQMPFCandidateAnalyzer.h.

std::vector<std::string> DQMPFCandidateAnalyzer::occupancyPFCand_nameRECO_

private

Definition at line 176 of file DQMPFCandidateAnalyzer.h.

std::vector<MonitorElement* > DQMPFCandidateAnalyzer::occupancyPFCand_puppiNolepWeight_

private

Definition at line 169 of file DQMPFCandidateAnalyzer.h.

std::vector<MonitorElement* > DQMPFCandidateAnalyzer::occupancyPFCandRECO_

private

Definition at line 175 of file DQMPFCandidateAnalyzer.h.

edm::EDGetTokenT<std::vector<pat::PackedCandidate> > DQMPFCandidateAnalyzer::pflowPackedToken_

private

Definition at line 116 of file DQMPFCandidateAnalyzer.h.

edm::EDGetTokenT<std::vector<reco::PFCandidate> > DQMPFCandidateAnalyzer::pflowToken_

private

Definition at line 115 of file DQMPFCandidateAnalyzer.h.

double DQMPFCandidateAnalyzer::ptMinCand_

private

Definition at line 128 of file DQMPFCandidateAnalyzer.h.

std::vector<MonitorElement* > DQMPFCandidateAnalyzer::ptPFCand_

private

Definition at line 167 of file DQMPFCandidateAnalyzer.h.

std::vector<std::string> DQMPFCandidateAnalyzer::ptPFCand_name_

private

Definition at line 168 of file DQMPFCandidateAnalyzer.h.

std::vector<std::string> DQMPFCandidateAnalyzer::ptPFCand_name_puppiNolepWeight_

private

Definition at line 170 of file DQMPFCandidateAnalyzer.h.

std::vector<std::string> DQMPFCandidateAnalyzer::ptPFCand_nameRECO_

private

Definition at line 176 of file DQMPFCandidateAnalyzer.h.

std::vector<MonitorElement* > DQMPFCandidateAnalyzer::ptPFCand_puppiNolepWeight_

private

Definition at line 169 of file DQMPFCandidateAnalyzer.h.

std::vector<MonitorElement* > DQMPFCandidateAnalyzer::ptPFCandRECO_

private

Definition at line 175 of file DQMPFCandidateAnalyzer.h.

edm::InputTag DQMPFCandidateAnalyzer::theTriggerResultsLabel_

private

Definition at line 99 of file DQMPFCandidateAnalyzer.h.

std::vector<int> DQMPFCandidateAnalyzer::typePFCand_

private

Definition at line 172 of file DQMPFCandidateAnalyzer.h.

std::vector<int> DQMPFCandidateAnalyzer::typePFCandRECO_

private

Definition at line 178 of file DQMPFCandidateAnalyzer.h.

int DQMPFCandidateAnalyzer::verbose_

private

Definition at line 181 of file DQMPFCandidateAnalyzer.h.

edm::InputTag DQMPFCandidateAnalyzer::vertexTag_

private

Definition at line 95 of file DQMPFCandidateAnalyzer.h.

edm::EDGetTokenT<std::vector<reco::Vertex> > DQMPFCandidateAnalyzer::vertexToken_

private

Definition at line 96 of file DQMPFCandidateAnalyzer.h.