---

pat::PatKitHelper Class Reference

#include <PhysicsTools/StarterKit/interface/PatKitHelper.h>

List of all members.

Public Member Functions
void	addNtupleVar (edm::EDProducer *prod, std::string name, std::string type)
void	bookHistos (edm::EDProducer *producer)
void	fillHistograms (edm::Event &event, edm::Handle< std::vector< pat::Muon > > &muonHandle, edm::Handle< std::vector< pat::Electron > > &electronHandle, edm::Handle< std::vector< pat::Tau > > &tauHandle, edm::Handle< std::vector< pat::Jet > > &jetHandle, edm::Handle< std::vector< pat::MET > > &METHandle, edm::Handle< std::vector< pat::Photon > > &photonHandle)
PhysicsHistograms::KinAxisLimits	getAxisLimits (std::string name)
void	getHandles (edm::Event &event, edm::Handle< std::vector< pat::Muon > > &muonHandle, edm::Handle< std::vector< pat::Electron > > &electronHandle, edm::Handle< std::vector< pat::Tau > > &tauHandle, edm::Handle< std::vector< pat::Jet > > &jetHandle, edm::Handle< std::vector< pat::MET > > &METHandle, edm::Handle< std::vector< pat::Photon > > &photonHandle)
	PatKitHelper (edm::ParameterSet const &parameters)
void	saveNtuple (edm::Event &event, const std::vector< pat::PhysVarHisto * > &ntvars)
template<class T>
void	saveNtupleVar (edm::Event &event, std::string name, T value)
template<class T>
void	saveNtupleVec (edm::Event &event, std::string name, const vector< T > &value)
template<class T>
void	saveNtupleVec (edm::Event &event, std::string name, const std::vector< T > &invec)
	~PatKitHelper ()
Public Attributes
edm::Service< TFileService >	fs_
pat::PhysVarHisto *	h_eventNumber_
pat::PhysVarHisto *	h_runNumber_
std::vector< pat::PhysVarHisto * >	ntVars_
edm::ParameterSet	parameters_
PhysicsHistograms *	physHistos_
int	verboseLevel_

---

Detailed Description

Definition at line 33 of file PatKitHelper.h.

---

Constructor & Destructor Documentation

PatKitHelper::PatKitHelper

(

edm::ParameterSet const &

parameters

)

Definition at line 7 of file PatKitHelper.cc.

                                                             :
  parameters_( iConfig ),
  physHistos_(0),
  verboseLevel_(0)
{
}

PatKitHelper::~PatKitHelper

(

)

Definition at line 15 of file PatKitHelper.cc.

References physHistos_.

{
  if ( physHistos_ ) delete physHistos_;
}

---

Member Function Documentation

void PatKitHelper::addNtupleVar	(	edm::EDProducer *	prod,
		std::string	name,
		std::string	type
	)

Definition at line 170 of file PatKitHelper.cc.

References GenMuonPlsPt100GeV_cfg::cout, and lat::endl().

Referenced by bookHistos().

{
  if      (type == "D") {
    producer->produces<double>( name ).setBranchAlias( name );
  }
  else if (type == "F") {
    producer->produces<float>( name ).setBranchAlias( name );
  }
  else if (type == "I") {
    producer->produces<int>( name ).setBranchAlias( name );
  }
  else if (type == "i") {
    producer->produces<unsigned int>( name ).setBranchAlias( name );
  }
  else if (type == "S") {
    producer->produces<short>( name ).setBranchAlias( name );
  }
  else if (type == "s") {
    producer->produces<unsigned short>( name ).setBranchAlias( name );
  }
  else if (type == "L") {
    producer->produces<long>( name ).setBranchAlias( name );
  }
  else if (type == "l") {
    producer->produces<unsigned long>( name ).setBranchAlias( name );
  }
  else if (type == "vD") {
    producer->produces<vector<double> >( name ).setBranchAlias( name );
  }
  else if (type == "vF") {
    producer->produces<vector<float> >( name ).setBranchAlias( name );
  }
  else if (type == "vI") {
    producer->produces<vector<int> >( name ).setBranchAlias( name );
  }
  else if (type == "vi") {
    producer->produces<vector<unsigned int> >( name ).setBranchAlias( name );
  }
  else if (type == "vS") {
    producer->produces<vector<short> >( name ).setBranchAlias( name );
  }
  else if (type == "vs") {
    producer->produces<vector<unsigned short> >( name ).setBranchAlias( name );
  }
  else if (type == "vL") {
    producer->produces<vector<long> >( name ).setBranchAlias( name );
  }
  else if (type == "vl") {
    producer->produces<vector<unsigned long> >( name ).setBranchAlias( name );
  }
  // &&& else if (type == "p4") {
  // &&&   producer->produces<math::XYZTLorentzVector> ( name ).setBranchAlias( name );
  // &&& }
  else {
    std::cout << "PatAnalyzerKit::addNtupleVar (ERROR): "
              << "unknown type " << type << std::endl;

    // &&& Throw an exception in order to abort the job!
  }
}

void PatKitHelper::bookHistos

(

edm::EDProducer *

producer

)

Definition at line 35 of file PatKitHelper.cc.

References addNtupleVar(), GenMuonPlsPt100GeV_cfg::cout, lat::endl(), getAxisLimits(), edm::ParameterSet::getParameter(), h_eventNumber_, h_runNumber_, pat::PhysVarHisto::makeTH1(), ntVars_, p, parameters_, physHistos_, and summary.

Referenced by CompositeKit::CompositeKit(), and PatAnalyzerKit::PatAnalyzerKit().

{
    // Initialize TFileService
  edm::Service<TFileService> fs;
  TFileDirectory summary = TFileDirectory( fs->mkdir("summary") );


  physHistos_ = new PhysicsHistograms ( getAxisLimits("muonAxis"),
                                        getAxisLimits("electronAxis"),
                                        getAxisLimits("tauAxis"),
                                        getAxisLimits("jetAxis"),
                                        getAxisLimits("METAxis"),
                                        getAxisLimits("photonAxis"),
                                        getAxisLimits("trackAxis")
                                        );

  // Get list of histograms to enable and disable
  string histos_to_disable =
    parameters_.getParameter<string>    ("disable");
  string histos_to_enable  =
    parameters_.getParameter<string>    ("enable");
  physHistos_->configure( histos_to_disable, histos_to_enable );

  // PhysicsHistograms takes ownership of the memory for these histograms here
  physHistos_->addHisto( h_runNumber_ = 
                        new PhysVarHisto( "runNumber", "Run Number",
                                          10000, 0, 10000, 
                                          &summary ,"", "I") 
                        );
  physHistos_->addHisto( h_eventNumber_ = 
                        new PhysVarHisto( "eventNumber", "Event Number",
                                          10000, 0, 10000, 
                                          &summary ,"", "I") 
                        );

  // Be sure to make the histograms!
  h_runNumber_->makeTH1();
  h_eventNumber_->makeTH1();

  // &&& Ntuple booking begin

  // Now that we know which variables are in the game, we could also
  // decide which ones to ntuplize
  string list_of_ntuple_vars =
    parameters_.getParameter<std::string>    ("ntuplize");

  if (list_of_ntuple_vars != "") {
    //
    //--- Collect all PhysVarHistos which need to store ntuple
    //--- variables and put them in here.
    physHistos_->select( list_of_ntuple_vars, ntVars_ );

    //--- Iterate over the list and "book" them via EDM
    std::vector< PhysVarHisto* >::iterator
      p    = ntVars_.begin(),
      pEnd = ntVars_.end();

    for ( ; p != pEnd; ++p ) {
      cout << "Adding ntuple variable " << (*p)->name() << endl;
      addNtupleVar( producer, (*p)->name(), (*p)->type() );
    }
    //
  } // end if

}

void PatKitHelper::fillHistograms	(	edm::Event &	event,
		edm::Handle< std::vector< pat::Muon > > &	muonHandle,
		edm::Handle< std::vector< pat::Electron > > &	electronHandle,
		edm::Handle< std::vector< pat::Tau > > &	tauHandle,
		edm::Handle< std::vector< pat::Jet > > &	jetHandle,
		edm::Handle< std::vector< pat::MET > > &	METHandle,
		edm::Handle< std::vector< pat::Photon > > &	photonHandle
	)

Definition at line 130 of file PatKitHelper.cc.

References PhysicsHistograms::clearVec(), PhysicsHistograms::fillCollection(), ntVars_, physHistos_, and saveNtuple().

Referenced by CompositeKit::produce(), and PatAnalyzerKit::produce().

{
  physHistos_->clearVec();
  physHistos_->fillCollection(*muonHandle);
  physHistos_->fillCollection(*electronHandle);
  physHistos_->fillCollection(*tauHandle);
  physHistos_->fillCollection(*jetHandle);
  physHistos_->fillCollection(*METHandle);
  physHistos_->fillCollection(*photonHandle);



  // save the list of ntuple varibles to the event record
  saveNtuple( event, ntVars_ );
}

PhysicsHistograms::KinAxisLimits PatKitHelper::getAxisLimits

(

std::string

name

)

Definition at line 20 of file PatKitHelper.cc.

References edm::ParameterSet::getParameter(), and parameters_.

Referenced by bookHistos(), and CompositeKit::CompositeKit().

{
  edm::ParameterSet axisLimitsSet = parameters_.getParameter<edm::ParameterSet>(name);
  
  PhysicsHistograms::KinAxisLimits axisLimits;

  axisLimits.pt1 = axisLimitsSet.getParameter<double>("pt1");
  axisLimits.pt2 = axisLimitsSet.getParameter<double>("pt2");
  axisLimits.m1  = axisLimitsSet.getParameter<double>("m1");
  axisLimits.m2  = axisLimitsSet.getParameter<double>("m2");

  return axisLimits;
}

void PatKitHelper::getHandles	(	edm::Event &	event,
		edm::Handle< std::vector< pat::Muon > > &	muonHandle,
		edm::Handle< std::vector< pat::Electron > > &	electronHandle,
		edm::Handle< std::vector< pat::Tau > > &	tauHandle,
		edm::Handle< std::vector< pat::Jet > > &	jetHandle,
		edm::Handle< std::vector< pat::MET > > &	METHandle,
		edm::Handle< std::vector< pat::Photon > > &	photonHandle
	)

Definition at line 102 of file PatKitHelper.cc.

References edm::ParameterSet::getParameter(), and parameters_.

Referenced by CompositeKit::produce(), and PatAnalyzerKit::produce().

{
  
  edm::InputTag muonName     = parameters_.getParameter<edm::InputTag>("muonSrc"    );
  edm::InputTag electronName = parameters_.getParameter<edm::InputTag>("electronSrc");
  edm::InputTag tauName      = parameters_.getParameter<edm::InputTag>("tauSrc"     );
  edm::InputTag jetName      = parameters_.getParameter<edm::InputTag>("jetSrc"     );
  edm::InputTag METName      = parameters_.getParameter<edm::InputTag>("METSrc"     );
  edm::InputTag photonName   = parameters_.getParameter<edm::InputTag>("photonSrc"  );

  

  event.getByLabel(muonName     , muonHandle);
  event.getByLabel(electronName , electronHandle);
  event.getByLabel(tauName      , tauHandle);
  event.getByLabel(jetName      , jetHandle);
  event.getByLabel(METName      , METHandle);
  event.getByLabel(photonName   , photonHandle);
}

void PatKitHelper::saveNtuple	(	edm::Event &	event,
		const std::vector< pat::PhysVarHisto * > &	ntvars
	)

Definition at line 238 of file PatKitHelper.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), p, and verboseLevel_.

Referenced by fillHistograms().

{
  //  Ntuplization

  if ( verboseLevel_ > 0 )
    cout << "About to save ntuple" << endl;
  if ( ntvars.size() ) {

    //--- Iterate over the list and "fill" them via EDM
    std::vector< PhysVarHisto* >::const_iterator
      p    = ntvars.begin(),
      pEnd = ntvars.end();

    for ( ; p != pEnd; ++p ) {

      if      ((*p)->type() == "D") {
        saveNtupleVar<double>( event, (*p)->name(), (*p)->value() );
      }
      else if ((*p)->type() == "F") {
        saveNtupleVar<float>( event, (*p)->name(), (*p)->value() );
      }
      else if ((*p)->type() == "I") {
        saveNtupleVar<int>( event, (*p)->name(), static_cast<int>((*p)->value()) );
      }
      else if ((*p)->type() == "i") {
        saveNtupleVar<unsigned int>( event, (*p)->name(), static_cast<unsigned int>((*p)->value()) );
      }
      else if ((*p)->type() == "S") {
        saveNtupleVar<short>( event, (*p)->name(), static_cast<short>((*p)->value()) );
      }
      else if ((*p)->type() == "s") {
        saveNtupleVar<unsigned short>( event, (*p)->name(), static_cast<unsigned short>((*p)->value()) );
      }
      else if ((*p)->type() == "L") {
        saveNtupleVar<long>( event, (*p)->name(), static_cast<long>((*p)->value()) );
      }
      else if ((*p)->type() == "l") {
        saveNtupleVar<unsigned long>( event, (*p)->name(), static_cast<unsigned long>((*p)->value()) );
      }
      else if ((*p)->type() == "vD") {
        vector<double> retvec;
        (*p)->vec( retvec );
        saveNtupleVec<double>( event, (*p)->name(), retvec );
      }
      else if ((*p)->type() == "vF") {
        vector<float> retvec;
        (*p)->vec( retvec );
        saveNtupleVec<float>( event, (*p)->name(), retvec );
      }
      else if ((*p)->type() == "vI") {
        vector<int> retvec;
        (*p)->vec( retvec );
        saveNtupleVec<int>( event, (*p)->name(), retvec );
      }
      else if ((*p)->type() == "vi") {
        vector<unsigned int> retvec;
        (*p)->vec( retvec );
        saveNtupleVec<unsigned int>( event, (*p)->name(), retvec );
      }
      else if ((*p)->type() == "vS") {
        vector<short> retvec;
        (*p)->vec( retvec );
        saveNtupleVec<short>( event, (*p)->name(), retvec );
      }
      else if ((*p)->type() == "vs") {
        vector<unsigned short> retvec;
        (*p)->vec( retvec );
        saveNtupleVec<unsigned short>( event, (*p)->name(), retvec );
      }
      else if ((*p)->type() == "vL") {
        vector<long> retvec;
        (*p)->vec( retvec );
        saveNtupleVec<long>( event, (*p)->name(), retvec );
      }
      else if ((*p)->type() == "vl") {
        vector<unsigned long> retvec;
        (*p)->vec( retvec );
        saveNtupleVec<unsigned long>( event, (*p)->name(), retvec );
      }
      // &&& else if (type == "p4") {
      // &&&   produces<math::XYZTLorentzVector> ( name ).setBranchAlias( name );
      // &&& }
      else {
        std::cout << "PatAnalyzerKit::addNtupleVar (ERROR): "
                  << "unknown type " << std::endl;

        // &&& Throw an exception in order to abort the job!
      }


      (*p)->clearVec(); // Clear ntuple cache here as well

    }
    //
  } // end if
}

template<class T>

void PatKitHelper::saveNtupleVar	(	edm::Event &	event,
		std::string	name,
		T	value
	)			[inline]

Definition at line 339 of file PatKitHelper.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), and verboseLevel_.

{
  std::auto_ptr<T> aptr( new T (value ) );
  if ( verboseLevel_ > 0 )
    cout << "Putting variable " << name << " with value " << value << " in event" << endl;
  event.put( aptr, name );
}

template<class T>

void pat::PatKitHelper::saveNtupleVec	(	edm::Event &	event,
		std::string	name,
		const vector< T > &	value
	)			[inline]

Definition at line 351 of file PatKitHelper.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), i, and verboseLevel_.

{
  std::auto_ptr<vector<T> > aptr( new vector<T> ( value ) );
  if ( verboseLevel_ > 0 ) {
    cout << "Putting variable " << name << " with values : " << endl;
    typename vector<T>::const_iterator i = aptr->begin(), iend = aptr->end();
    for ( ; i != iend; ++i ) {
      cout << *i << " ";
    }
    cout << endl;
  }
  event.put( aptr, name );
}

template<class T>

void pat::PatKitHelper::saveNtupleVec	(	edm::Event &	event,
		std::string	name,
		const std::vector< T > &	invec
	)			[inline]

---

Member Data Documentation

edm::Service<TFileService> pat::PatKitHelper::fs_

Definition at line 96 of file PatKitHelper.h.

pat::PhysVarHisto* pat::PatKitHelper::h_eventNumber_

Definition at line 105 of file PatKitHelper.h.

Referenced by bookHistos().

pat::PhysVarHisto* pat::PatKitHelper::h_runNumber_

Definition at line 104 of file PatKitHelper.h.

Referenced by bookHistos().

std::vector< pat::PhysVarHisto* > pat::PatKitHelper::ntVars_

Definition at line 99 of file PatKitHelper.h.

Referenced by bookHistos(), and fillHistograms().

edm::ParameterSet pat::PatKitHelper::parameters_

Definition at line 90 of file PatKitHelper.h.

Referenced by bookHistos(), getAxisLimits(), and getHandles().

PhysicsHistograms* pat::PatKitHelper::physHistos_

Definition at line 93 of file PatKitHelper.h.

Referenced by bookHistos(), fillHistograms(), and ~PatKitHelper().

int pat::PatKitHelper::verboseLevel_

Definition at line 87 of file PatKitHelper.h.

Referenced by saveNtuple(), saveNtupleVar(), and saveNtupleVec().

---

The documentation for this class was generated from the following files:

• PhysicsTools/StarterKit/interface/PatKitHelper.h
• PhysicsTools/StarterKit/src/PatKitHelper.cc

---