---

PhysicsHistograms.cc

Go to the documentation of this file.

//------------------------------------------------------------------------
// -*- C++ -*-
//
//
// Original Author:  Petar Maksimovic
//         Created:  Christmas 2007
// $Id: PhysicsHistograms.cc,v 1.5 2008/07/07 20:06:56 srappocc Exp $
//
// Revision History:
//------------------------------------------------------------------------


#include "PhysicsTools/StarterKit/interface/PhysicsHistograms.h"

#include <string>
#include <sstream>

//------------------------------------------------------------------------
//------------------------------------------------------------------------
PhysicsHistograms::PhysicsHistograms( KinAxisLimits const & muonAxis, 
                                      KinAxisLimits const & electronAxis, 
                                      KinAxisLimits const & tauAxis, 
                                      KinAxisLimits const & jetAxis, 
                                      KinAxisLimits const & METAxis, 
                                      KinAxisLimits const & photonAxis, 
                                      KinAxisLimits const & trackAxis
                                      )
{
  //--- Initialize histogram objects
  std::cout << "PhysicsHistograms: Creating muon histograms" << std::endl;
  muonHistograms_     = 
    new pat::HistoMuon    ( "muon", "Muon", "mu", muonAxis.pt1, muonAxis.pt2, muonAxis.m1, muonAxis.m2 );
  std::cout << "PhysicsHistograms: Creating electron histograms" << std::endl;
  electronHistograms_ = 
    new pat::HistoElectron( "electron", "Electron", "e", electronAxis.pt1, electronAxis.pt2, electronAxis.m1, electronAxis.m2 );
  std::cout << "PhysicsHistograms: Creating tau histograms" << std::endl;
  tauHistograms_      = 
    new pat::HistoTau     ( "tau", "Tau", "tau", tauAxis.pt1, tauAxis.pt2, tauAxis.m1, tauAxis.m2 );
  std::cout << "PhysicsHistograms: Creating jet histograms" << std::endl;
  jetHistograms_      =
    new pat::HistoJet     ( "jet", "Jet", "jet", jetAxis.pt1, jetAxis.pt2, jetAxis.m1, jetAxis.m2 );
  std::cout << "PhysicsHistograms: Creating met histograms" << std::endl;
  metHistograms_      = 
    new pat::HistoMET     ( "met", "MET", "met", METAxis.pt1, METAxis.pt2, METAxis.m1, METAxis.m2 );
  std::cout << "PhysicsHistograms: Creating photon histograms" << std::endl;
  photonHistograms_   = 
    new pat::HistoPhoton  ( "photon", "Photon", "photon", photonAxis.pt1, photonAxis.pt2, photonAxis.m1, photonAxis.m2 );
  std::cout << "PhysicsHistograms: Creating track histograms" << std::endl;
  trackHistograms_    = 
    new pat::HistoTrack   ( "track", "Track", "track", trackAxis.pt1, trackAxis.pt2, trackAxis.m1, trackAxis.m2 );
}


//------------------------------------------------------------------------
//------------------------------------------------------------------------
PhysicsHistograms::~PhysicsHistograms()
{
  delete muonHistograms_     ;
  delete electronHistograms_ ;
  delete tauHistograms_      ;
  delete jetHistograms_      ;
  delete metHistograms_      ;
  delete photonHistograms_   ;
  delete trackHistograms_    ;

  for ( unsigned int i = 0; i < allVarHistos_.size(); ++i ) {
    if ( allVarHistos_[i] ) delete allVarHistos_[i];
  }

  outputFile_.close();
}



//------------------------------------------------------------------------
//------------------------------------------------------------------------
void
PhysicsHistograms::configure( std::string & histos_to_disable,   // comma separated list of names
                              std::string & histos_to_enable )   // comma separated list of names
{
  std::cout << "PhysicsHistograms:: configuring..."
            << "\n   First disabling: " << histos_to_disable
            << "\n   Then  enabling : " << histos_to_enable
            << std::endl;


  //--- Pass this information to histogramGroups
  muonHistograms_    ->configure( histos_to_disable, histos_to_enable ) ;
  electronHistograms_->configure( histos_to_disable, histos_to_enable ) ;
  tauHistograms_     ->configure( histos_to_disable, histos_to_enable ) ;
  metHistograms_     ->configure( histos_to_disable, histos_to_enable ) ;
  jetHistograms_     ->configure( histos_to_disable, histos_to_enable ) ;
  photonHistograms_  ->configure( histos_to_disable, histos_to_enable ) ;
  trackHistograms_   ->configure( histos_to_disable, histos_to_enable ) ;

}


//------------------------------------------------------------------------
//------------------------------------------------------------------------
void
PhysicsHistograms::select( std::string  vars_to_select,   // comma separated list of names
                           std::vector< pat::PhysVarHisto * > & selectedVars )
{
  std::cout << "PhysicsHistograms:: selecting the following variables:\n\t"
            << vars_to_select
            << std::endl;


  //--- Pass this information to histogramGroups
  muonHistograms_    ->select( vars_to_select, selectedVars ) ;
  electronHistograms_->select( vars_to_select, selectedVars ) ;
  tauHistograms_     ->select( vars_to_select, selectedVars ) ;
  metHistograms_     ->select( vars_to_select, selectedVars ) ;
  jetHistograms_     ->select( vars_to_select, selectedVars ) ;
  photonHistograms_  ->select( vars_to_select, selectedVars ) ;
  trackHistograms_   ->select( vars_to_select, selectedVars ) ;
  
  std::vector<pat::PhysVarHisto*>::iterator i = allVarHistos_.begin();
  std::vector<pat::PhysVarHisto*>::iterator end = allVarHistos_.end();
  std::string temp = "," + vars_to_select + ",";
  for ( ; i != end; ++i  ) {
    std::string test = "," + (*i)->name() + ",";
    std::cout << "testing " << test << std::endl;
    if ( temp.find( test ) != std::string::npos || temp == ",all," ) {
      std::cout << "FOUND!" << std::endl;
      selectedVars.push_back ( *i );
    }
  }

  std::cout << "PhysicsHistograms:: selected " << selectedVars.size()
            << " variables." << std::endl;
}

void PhysicsHistograms::clearVec()
{
  muonHistograms_    ->clearVec() ;
  electronHistograms_->clearVec() ;
  tauHistograms_     ->clearVec() ;
  metHistograms_     ->clearVec() ;
  jetHistograms_     ->clearVec() ;
  photonHistograms_  ->clearVec() ;
  trackHistograms_   ->clearVec() ;
  for ( uint i = 0; i < allVarHistos_.size(); i++ ) {
    allVarHistos_[i]->clearVec();
  }
}

//------------------------------------------------------------------------
//------------------------------------------------------------------------
void
PhysicsHistograms::beginJob()
{
  // Dummy for now
}


//------------------------------------------------------------------------
//------------------------------------------------------------------------
void
PhysicsHistograms::endJob()
{
  // Dummy for now
}



//------------------------------------------------------------------------
// &&& Design suggestion: this should go into HistoGroup<> instead.
//------------------------------------------------------------------------
std::ostream & operator<<( std::ostream & out, const reco::Candidate & cand )
{
  char buff[1000];
  sprintf( buff, "Pt, Eta, Phi, M = (%6.2f, %6.2f, %6.2f, %6.2f)",
           cand.pt(), cand.eta(), cand.phi(), cand.mass() );
  out << buff;
  return out;
}



//--- Code graveyard:
//
//   // Muons
//   typedef std::vector<Muon>::const_iterator muonIter ;
//   outputFile_ << "Muons: " << muons.size() << endl;
//   int imuon = 0;
//   for (muonIter muon = muons.begin(); muon != muons.end();++muon, ++imuon) {
//     outputFile_ << setw(6) << imuon << " : " << *muon << endl;
//     histoMuon_->fill( *muon, imuon+1 );
//   }

---