FSQDQM Class Reference

#include <FSQDQM.h>

Inheritance diagram for FSQDQM:

Public Member Functions
	FSQDQM (const edm::ParameterSet &ps)
	~FSQDQM () override
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override

Private Member Functions
void	bookHistograms (DQMStore::IBooker &bei, edm::Run const &, edm::EventSetup const &) override
void	bookHistos (DQMStore *bei)

Private Attributes
unsigned int	bxNumber_
MonitorElement *	CastorJetMulti
MonitorElement *	CastorJetphi
unsigned int	eventNumber_
MonitorElement *	h_leadingtrkpt_ntrk_away
MonitorElement *	h_leadingtrkpt_ntrk_towards
MonitorElement *	h_leadingtrkpt_ntrk_transverse
MonitorElement *	h_leadingtrkpt_ptsum_away
MonitorElement *	h_leadingtrkpt_ptsum_towards
MonitorElement *	h_leadingtrkpt_ptsum_transverse
MonitorElement *	h_ntracks
MonitorElement *	h_ntracks_away
MonitorElement *	h_ntracks_towards
MonitorElement *	h_ntracks_transverse
MonitorElement *	h_ptsum_away
MonitorElement *	h_ptsum_towards
MonitorElement *	h_ptsum_transverse
MonitorElement *	h_trkptsum
std::vector< int >	hltresults
std::string	labelBS_
std::string	labelCastorJet_
std::string	labelPFJet_
std::string	labelTrack_
unsigned int	lumiNumber_
MonitorElement *	NPV
MonitorElement *	PFJeteta
MonitorElement *	PFJetMulti
MonitorElement *	PFJetphi
MonitorElement *	PFJetpt
MonitorElement *	PFJetRapidity
MonitorElement *	PV_chi2
MonitorElement *	PV_d0
MonitorElement *	PV_numTrks
MonitorElement *	PV_sumTrks
edm::EDGetTokenT< edm::View< reco::Vertex > >	pvs_
unsigned int	runNumber_
edm::EDGetTokenT< reco::BasicJetCollection >	tok_castorjet_
edm::EDGetTokenT< reco::PFJetCollection >	tok_pfjet_
edm::EDGetTokenT< reco::TrackCollection >	tok_track_
MonitorElement *	Track_HP_dxyvtx_over_dxyerror
MonitorElement *	Track_HP_dzvtx_over_dzerr
MonitorElement *	Track_HP_Eta
MonitorElement *	Track_HP_Phi
MonitorElement *	Track_HP_Pt
MonitorElement *	Track_HP_ptErr_over_pt
edm::InputTag	vertex_

Detailed Description

Definition at line 80 of file FSQDQM.h.

Constructor & Destructor Documentation

FSQDQM::FSQDQM

(

const edm::ParameterSet &

ps

)

Definition at line 90 of file FSQDQM.cc.

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

{
  edm::LogInfo("FSQDQM") << " Creating FSQDQM "
                 << "\n";
  pvs_                        =consumes<edm::View<reco::Vertex> >(
                                 iConfig.getParameter<edm::InputTag>("pvs"));
  labelPFJet_                 = iConfig.getParameter<std::string>("LabelPFJet");
  labelCastorJet_             = iConfig.getParameter<std::string>("LabelCastorJet");
  labelTrack_                 = iConfig.getParameter<std::string>("LabelTrack");
  tok_pfjet_                  = consumes<reco::PFJetCollection>(labelPFJet_);
  tok_track_                  = consumes<reco::TrackCollection>(labelTrack_);
  tok_castorjet_              = consumes<reco::BasicJetCollection>(labelCastorJet_);

}

FSQDQM::~FSQDQM ( )

override

Definition at line 106 of file FSQDQM.cc.

 {

   edm::LogInfo("FSQDQM") << " Deleting FSQDQM "
                << "\n";
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)

 }

Member Function Documentation

void FSQDQM::analyze ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

overrideprotected

Definition at line 162 of file FSQDQM.cc.

References funct::abs(), edm::EventBase::bunchCrossing(), hiPixelPairStep_cff::deltaPhi, edm::EventID::event(), edm::Event::getByToken(), edm::EventBase::id(), reco::Vertex::isFake(), edm::HandleBase::isValid(), gen::k, edm::EventID::luminosityBlock(), reco::Vertex::normalizedChi2(), reco::Vertex::position(), edm::Handle< T >::product(), reco::TrackBase::qualityByName(), RecoPFJets_cff::recoPFJets, edm::EventID::run(), sistrip::runNumber_, mathSSE::sqrt(), reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), reco::Vertex::tracksSize(), reco::Vertex::x(), reco::Vertex::xError(), reco::Vertex::y(), reco::Vertex::yError(), reco::Vertex::z(), and reco::Vertex::zError().

 {
   using namespace edm;
   using namespace std;
   using namespace reco;

    runNumber_   = iEvent.id().run();
    eventNumber_ = iEvent.id().event();
    lumiNumber_  = iEvent.id().luminosityBlock();
    bxNumber_ = iEvent.bunchCrossing();

    edm::Handle<edm::View<reco::Vertex> > privtxs;
    //    if(!iEvent.getByToken(pvs_, privtxs)) return;
    iEvent.getByToken(pvs_, privtxs);
    double bestvz=-999.9, bestvx=-999.9, bestvy=-999.9;
    double bestvzError=-999.9, bestvxError=-999.9, bestvyError=-999.9;
    if(privtxs.isValid()){
      const reco::Vertex& pvtx=privtxs->front();
      NPV->Fill(privtxs->size());
      if(privtxs->begin() !=privtxs->end() && !(pvtx.isFake()) && pvtx.position().Rho() <= 2. && fabs(pvtx.position().z()) <= 24){
    bestvz = pvtx.z(); 
    bestvx = pvtx.x(); 
    bestvy = pvtx.y();
    bestvzError = pvtx.zError(); 
    bestvxError = pvtx.xError(); 
    bestvyError = pvtx.yError();
    PV_chi2->Fill(pvtx.normalizedChi2());
    PV_d0->Fill(sqrt(pvtx.x() * pvtx.x() + pvtx.y() * pvtx.y()));
    PV_numTrks->Fill(pvtx.tracksSize());
    double vertex_sumTrks = 0.0;
    for(reco::Vertex::trackRef_iterator iTrack= pvtx.tracks_begin(); iTrack != pvtx.tracks_end(); iTrack++)
      {
        vertex_sumTrks += (*iTrack)->pt();
      }
    PV_sumTrks->Fill(vertex_sumTrks);
      }
    }
    
    std::vector<Jet> recoPFJets;
     recoPFJets.clear();
     int  nPFCHSJet=0;
     edm::Handle<PFJetCollection> pfjetchscoll;
     //  if(!iEvent.getByToken(tok_pfjet_, pfjetchscoll)) return;
     iEvent.getByToken(tok_pfjet_, pfjetchscoll);
     if(pfjetchscoll.isValid()){
     const reco::PFJetCollection *pfchsjets = pfjetchscoll.product();
     reco::PFJetCollection::const_iterator pfjetchsclus = pfchsjets->begin();
     for(pfjetchsclus = pfchsjets->begin(); pfjetchsclus!= pfchsjets->end() ; ++pfjetchsclus){
       PFJetpt->Fill( pfjetchsclus->pt());
       PFJeteta->Fill( pfjetchsclus->eta());
       PFJetphi->Fill( pfjetchsclus->phi());
       PFJetRapidity->Fill( pfjetchsclus->rapidity());
       nPFCHSJet++;
     }
     PFJetMulti->Fill( nPFCHSJet);
     }
 
     std::vector<Jet> recoCastorJets;
         recoCastorJets.clear();
     edm::Handle<BasicJetCollection> castorJets;
     //  if(!iEvent.getByToken(tok_castorjet_, castorJets)) return;
     iEvent.getByToken(tok_castorjet_, castorJets);
     if(castorJets.isValid()){
     for (unsigned ijet=0; ijet<castorJets->size();ijet++) {
       recoCastorJets.push_back((*castorJets)[ijet]);
     }
     for (unsigned ijet=0; ijet<recoCastorJets.size(); ijet++) {
       //      cout<<recoCastorJets[ijet].pt()<<endl;
       CastorJetphi->Fill( recoCastorJets[ijet].phi());

       CastorJetMulti->Fill(recoCastorJets.size());
     }
     }
     edm::Handle<reco::TrackCollection> itracks;
     //  if(!iEvent.getByToken(tok_track_, itracks)) return;
     iEvent.getByToken(tok_track_, itracks);
     if(itracks.isValid()){
     reco::TrackBase::TrackQuality hiPurity = reco::TrackBase::qualityByName("highPurity");
     std::vector<TLorentzVector> T_trackRec_P4;

     int ntracks = 0;
     int ntracks_towards = 0;
     int ntracks_transverse = 0;
     int ntracks_away = 0;
     float ptsum = 0;
     float dphi=0;
     float ptsum_towards = 0;
     float ptsum_transverse = 0;
     float ptsum_away = 0;
     
     T_trackRec_P4.clear();
     for(reco::TrackCollection::const_iterator iT = itracks->begin(); iT != itracks->end(); ++iT){
         if(iT->quality(hiPurity)){
           math::XYZPoint bestvtx(bestvx,bestvy,bestvz);
           double dzvtx = iT->dz(bestvtx);
           double dxyvtx = iT->dxy(bestvtx);
           double dzerror = sqrt(iT->dzError()*iT->dzError()+bestvzError*bestvzError);
           double dxyerror = sqrt(iT->d0Error()*iT->d0Error()+bestvxError*bestvyError);
           if((iT->ptError())/iT->pt() < 0.05 && dzvtx < 3.0 && dxyvtx < 3.0){
           TLorentzVector trk; 
           trk.SetPtEtaPhiE(iT->pt(),iT->eta(),iT->phi(),iT->p());
           T_trackRec_P4.push_back(trk);
               Track_HP_Eta->Fill(iT->eta());
           Track_HP_Phi->Fill(iT->phi());
           Track_HP_Pt->Fill(iT->pt());
           Track_HP_ptErr_over_pt->Fill((iT->ptError())/iT->pt());
           Track_HP_dzvtx_over_dzerr->Fill(dzvtx/dzerror);
           Track_HP_dxyvtx_over_dxyerror->Fill(dxyvtx/dxyerror); 
           }
         }
     }

     float highest_pt_track=-999;
     int index=-999;
     for(unsigned int k=0;k<T_trackRec_P4.size();k++){
       if(T_trackRec_P4.at(k).Pt() > highest_pt_track){
         highest_pt_track=T_trackRec_P4.at(k).Pt();
         index = k;
           }

     }
     unsigned int finalid=abs(index);
     //  if(T_trackRec_P4.at(index).Pt()!=0.){
     if(finalid < T_trackRec_P4.size()){
       //    std::sort(T_trackRec_P4.begin(), T_trackRec_P4.end(), SortByPt());
     for(unsigned int itrk=0;itrk<T_trackRec_P4.size();itrk++){
       ++ntracks;
         ptsum= ptsum + T_trackRec_P4.at(itrk).Pt();
         dphi = deltaPhi(T_trackRec_P4.at(itrk).Phi(),T_trackRec_P4.at(index).Phi());
         if(fabs(dphi) < 1.05){
           ++ntracks_towards;
           ptsum_towards = ptsum_towards + T_trackRec_P4.at(itrk).Pt();}
         if(fabs(dphi) > 1.05 && fabs(dphi) < 2.09){
           ++ntracks_transverse;
           ptsum_transverse = ptsum_transverse + T_trackRec_P4.at(itrk).Pt();}
         if(fabs(dphi) > 2.09){
           ++ntracks_away; 
           ptsum_away = ptsum_away + T_trackRec_P4.at(itrk).Pt();}
       }       
     
     
     
     h_ntracks->Fill(ntracks);
         h_trkptsum->Fill(ptsum);
     h_ptsum_towards->Fill(ptsum_towards);
     h_ptsum_transverse->Fill(ptsum_transverse);
     h_ptsum_away->Fill(ptsum_away);
     h_ntracks_towards->Fill(ntracks_towards);
     h_ntracks_transverse->Fill(ntracks_transverse);
     h_ntracks_away->Fill(ntracks_away);

       if(!T_trackRec_P4.empty()){
     h_leadingtrkpt_ntrk_towards->Fill(T_trackRec_P4.at(index).Pt(),ntracks_towards/8.37);
     h_leadingtrkpt_ntrk_transverse->Fill(T_trackRec_P4.at(index).Pt(),ntracks_transverse/8.37);
     h_leadingtrkpt_ntrk_away->Fill(T_trackRec_P4.at(index).Pt(),ntracks_away/8.37);
     h_leadingtrkpt_ptsum_towards->Fill(T_trackRec_P4.at(index).Pt(),ptsum_towards/8.37);
     h_leadingtrkpt_ptsum_transverse->Fill(T_trackRec_P4.at(index).Pt(),ptsum_transverse/8.37);
     h_leadingtrkpt_ptsum_away->Fill(T_trackRec_P4.at(index).Pt(),ptsum_away/8.37);
       }
     }
     }
 }//analyze

void FSQDQM::bookHistograms ( DQMStore::IBooker &  bei, edm::Run const &  , edm::EventSetup const &    )

overrideprivate

Definition at line 115 of file FSQDQM.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::bookProfile(), and DQMStore::IBooker::setCurrentFolder().

                                       {
   bei.setCurrentFolder("Physics/FSQ");
   PFJetpt = bei.book1D("PFJetpt",";p_{T}(PFJet)", 100,0.0 , 100);
   PFJeteta = bei.book1D("PFJeteta", ";#eta(PFJet)", 50, -2.5, 2.5);
   PFJetphi = bei.book1D("PFJetphi", ";#phi(PFJet)", 50, -3.14,3.14);
   PFJetMulti         = bei.book1D("PFJetMulti", ";No. of PFJets", 10, -0.5, 9.5);
   PFJetRapidity      = bei.book1D("PFJetRapidity",";PFJetRapidity", 50, -6.0,6.0);
   CastorJetphi = bei.book1D("CastorJetphi", ";#phi(CastorJet)", 50, -3.14,3.14);
   CastorJetMulti         = bei.book1D("CastorJetMulti", ";No. of CastorJets", 10, -0.5, 9.5);
   Track_HP_Phi       =bei.book1D("Track_HP_Phi",";#phi(HPtrack)", 50, -3.14,3.14);
   Track_HP_Eta       =bei.book1D("Track_HP_Eta", ";#eta(HPtrack)", 50, -2.5, 2.5);
   Track_HP_Pt        =bei.book1D("Track_HP_Pt",  ";p_{T}(HPtrack)",500, 0.0 , 50);
   Track_HP_ptErr_over_pt=bei.book1D("Track_HP_ptErr_over_pt",";{p_{T}Err}/p_{T}",100,0,0.1);
   Track_HP_dzvtx_over_dzerr=bei.book1D("Track_HP_dzvtx_over_dzerr",";dZerr/dZ",100,-10,10);
   Track_HP_dxyvtx_over_dxyerror =bei.book1D("Track_HP_dxyvtx_over_dxyerror",";dxyErr/dxy",100,-10,10);
   NPV      = bei.book1D("NPV",";NPV",10, -0.5, 9.5);
   PV_chi2  = bei.book1D("PV_chi2",";PV_chi2",100, 0.0, 2.0);
   PV_d0    = bei.book1D("PV_d0",";PV_d0",100, -10.0, 10.0);
   PV_numTrks         = bei.book1D("PV_numTrks",";PV_numTrks",100, -0.5, 99.5);
   PV_sumTrks=bei.book1D("PV_sumTrks",";PV_sumTrks",100,0,100);
   h_ptsum_towards      = bei.book1D("h_ptsum_towards",";h_ptsum_towards",100,0,100);
   h_ptsum_transverse      = bei.book1D("h_ptsum_transverse",";h_ptsum_transverse",100,0,100);

   h_ntracks      = bei.book1D("h_ntracks",";h_ntracks",50,-0.5,49.5);
   h_trkptsum      = bei.book1D("h_trkptsum",";h_trkptsum",100,0,100);
   h_ptsum_away      = bei.book1D("h_ptsum_away",";h_ptsum_away",100,0,100);
   h_ntracks_towards      = bei.book1D("h_ntracks_towards",";h_ntracks_towards",50,-0.5,49.5);
   h_ntracks_transverse      = bei.book1D("h_ntracks_transverse",";h_ntracks_transverse",50,-0.5,49.5);
   h_ntracks_away      = bei.book1D("h_ntracks_away",";h_ntracks_away",50,-0.5,49.5);
   
   h_leadingtrkpt_ntrk_away =bei.bookProfile("h_leadingtrkpt_ntrk_away","h_leadingtrkpt_ntrk_away",50,0,50,0,30," ");
   h_leadingtrkpt_ntrk_towards =bei.bookProfile("h_leadingtrkpt_ntrk_towards","h_leadingtrkpt_ntrk_towards",50,0,50,0,30," ");
   h_leadingtrkpt_ntrk_transverse =bei.bookProfile("h_leadingtrkpt_ntrk_transverse","h_leadingtrkpt_ntrk_transverse",50,0,50,0,30," ");
   h_leadingtrkpt_ptsum_away =bei.bookProfile("h_leadingtrkpt_ptsum_away","h_leadingtrkpt_ptsum_away",50,0,50,0,30," ");
   h_leadingtrkpt_ptsum_towards =bei.bookProfile("h_leadingtrkpt_ptsum_towards","h_leadingtrkpt_ptsum_towards",50,0,50,0,30," ");
   h_leadingtrkpt_ptsum_transverse =bei.bookProfile("h_leadingtrkpt_ptsum_transverse","h_leadingtrkpt_ptsum_transverse",50,0,50,0,30," ");
     
 }

void FSQDQM::bookHistos ( DQMStore *  bei )

private

Member Data Documentation

unsigned int FSQDQM::bxNumber_

private

Definition at line 103 of file FSQDQM.h.

MonitorElement* FSQDQM::CastorJetMulti

private

Definition at line 112 of file FSQDQM.h.

MonitorElement* FSQDQM::CastorJetphi

private

Definition at line 111 of file FSQDQM.h.

unsigned int FSQDQM::eventNumber_

private

Definition at line 103 of file FSQDQM.h.

MonitorElement* FSQDQM::h_leadingtrkpt_ntrk_away

private

Definition at line 135 of file FSQDQM.h.

MonitorElement* FSQDQM::h_leadingtrkpt_ntrk_towards

private

Definition at line 136 of file FSQDQM.h.

MonitorElement* FSQDQM::h_leadingtrkpt_ntrk_transverse

private

Definition at line 137 of file FSQDQM.h.

MonitorElement* FSQDQM::h_leadingtrkpt_ptsum_away

private

Definition at line 138 of file FSQDQM.h.

MonitorElement* FSQDQM::h_leadingtrkpt_ptsum_towards

private

Definition at line 139 of file FSQDQM.h.

MonitorElement* FSQDQM::h_leadingtrkpt_ptsum_transverse

private

Definition at line 140 of file FSQDQM.h.

MonitorElement* FSQDQM::h_ntracks

private

Definition at line 133 of file FSQDQM.h.

MonitorElement* FSQDQM::h_ntracks_away

private

Definition at line 131 of file FSQDQM.h.

MonitorElement* FSQDQM::h_ntracks_towards

private

Definition at line 129 of file FSQDQM.h.

MonitorElement* FSQDQM::h_ntracks_transverse

private

Definition at line 130 of file FSQDQM.h.

MonitorElement* FSQDQM::h_ptsum_away

private

Definition at line 128 of file FSQDQM.h.

MonitorElement* FSQDQM::h_ptsum_towards

private

Definition at line 126 of file FSQDQM.h.

MonitorElement* FSQDQM::h_ptsum_transverse

private

Definition at line 127 of file FSQDQM.h.

MonitorElement* FSQDQM::h_trkptsum

private

Definition at line 132 of file FSQDQM.h.

std::vector<int> FSQDQM::hltresults

private

Definition at line 102 of file FSQDQM.h.

std::string FSQDQM::labelBS_

private

Definition at line 96 of file FSQDQM.h.

std::string FSQDQM::labelCastorJet_

private

Definition at line 96 of file FSQDQM.h.

std::string FSQDQM::labelPFJet_

private

Definition at line 96 of file FSQDQM.h.

std::string FSQDQM::labelTrack_

private

Definition at line 96 of file FSQDQM.h.

unsigned int FSQDQM::lumiNumber_

private

Definition at line 103 of file FSQDQM.h.

MonitorElement* FSQDQM::NPV

private

Definition at line 121 of file FSQDQM.h.

MonitorElement* FSQDQM::PFJeteta

private

Definition at line 107 of file FSQDQM.h.

MonitorElement* FSQDQM::PFJetMulti

private

Definition at line 113 of file FSQDQM.h.

MonitorElement* FSQDQM::PFJetphi

private

Definition at line 108 of file FSQDQM.h.

MonitorElement* FSQDQM::PFJetpt

private

Definition at line 106 of file FSQDQM.h.

MonitorElement* FSQDQM::PFJetRapidity

private

Definition at line 114 of file FSQDQM.h.

MonitorElement* FSQDQM::PV_chi2

private

Definition at line 122 of file FSQDQM.h.

MonitorElement* FSQDQM::PV_d0

private

Definition at line 123 of file FSQDQM.h.

MonitorElement* FSQDQM::PV_numTrks

private

Definition at line 124 of file FSQDQM.h.

MonitorElement* FSQDQM::PV_sumTrks

private

Definition at line 125 of file FSQDQM.h.

edm::EDGetTokenT<edm::View<reco::Vertex> > FSQDQM::pvs_

private

Definition at line 97 of file FSQDQM.h.

unsigned int FSQDQM::runNumber_

private

Definition at line 103 of file FSQDQM.h.

edm::EDGetTokenT<reco::BasicJetCollection> FSQDQM::tok_castorjet_

private

Definition at line 100 of file FSQDQM.h.

edm::EDGetTokenT<reco::PFJetCollection> FSQDQM::tok_pfjet_

private

Definition at line 99 of file FSQDQM.h.

edm::EDGetTokenT<reco::TrackCollection> FSQDQM::tok_track_

private

Definition at line 98 of file FSQDQM.h.

MonitorElement* FSQDQM::Track_HP_dxyvtx_over_dxyerror

private

Definition at line 120 of file FSQDQM.h.

MonitorElement* FSQDQM::Track_HP_dzvtx_over_dzerr

private

Definition at line 119 of file FSQDQM.h.

MonitorElement* FSQDQM::Track_HP_Eta

private

Definition at line 116 of file FSQDQM.h.

MonitorElement* FSQDQM::Track_HP_Phi

private

Definition at line 115 of file FSQDQM.h.

MonitorElement* FSQDQM::Track_HP_Pt

private

Definition at line 117 of file FSQDQM.h.

MonitorElement* FSQDQM::Track_HP_ptErr_over_pt

private

Definition at line 118 of file FSQDQM.h.

edm::InputTag FSQDQM::vertex_

private

Definition at line 95 of file FSQDQM.h.