#include <CommonTools/UtilAlgos/src/ConfigurableAnalysis.cc>

Inheritance diagram for ConfigurableAnalysis:

Public Member Functions
	ConfigurableAnalysis (const edm::ParameterSet &)

	~ConfigurableAnalysis ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	endJob ()

virtual bool	filter (edm::Event &, const edm::EventSetup &)

Private Attributes
std::vector< std::string >	flows_

NTupler *	ntupler_

Plotter *	plotter_

Selections *	selections_

bool	workAsASelector_

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

typedef WorkerT< EDFilter >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 45 of file ConfigurableAnalysis.cc.

Constructor & Destructor Documentation

ConfigurableAnalysis::ConfigurableAnalysis ( const edm::ParameterSet & iConfig )

explicit

Definition at line 74 of file ConfigurableAnalysis.cc.

References edm::ParameterSet::empty(), edm::ParameterSet::exists(), flows_, reco::get(), edm::ParameterSet::getParameter(), init, moduleLabel(), ntupler_, plotter_, NTupler::registerleaves(), selections_, and workAsASelector_.

                                                                          :
   selections_(0), plotter_(0), ntupler_(0)
 {
 
   std::string moduleLabel = iConfig.getParameter<std::string>("@module_label");
 
   //configure inputag distributor
   if (iConfig.exists("InputTags"))
     edm::Service<InputTagDistributorService>()->init(moduleLabel,iConfig.getParameter<edm::ParameterSet>("InputTags"));
 
   //configure the variable helper
   edm::Service<VariableHelperService>()->init(moduleLabel,iConfig.getParameter<edm::ParameterSet>("Variables"));
 
   //list of selections
   selections_ = new Selections(iConfig.getParameter<edm::ParameterSet>("Selections"));
 
   //plotting device
   edm::ParameterSet plotPset = iConfig.getParameter<edm::ParameterSet>("Plotter");
   if (!plotPset.empty()){
     std::string plotterName = plotPset.getParameter<std::string>("ComponentName");
     plotter_ = PlotterFactory::get()->create(plotterName, plotPset);
   }
   else
     plotter_ = 0;
 
   //ntupling device
   edm::ParameterSet ntPset = iConfig.getParameter<edm::ParameterSet>("Ntupler");
   if (!ntPset.empty()){
     std::string ntuplerName=ntPset.getParameter<std::string>("ComponentName");
     ntupler_ = NTuplerFactory::get()->create(ntuplerName, ntPset);
   }
   else ntupler_=0;
   
   flows_ = iConfig.getParameter<std::vector<std::string> >("flows");
   workAsASelector_ = iConfig.getParameter<bool>("workAsASelector");
 
   //vector of passed selections
   produces<std::vector<bool> >();
 
   //ntupler needs to register its products
   if (ntupler_) ntupler_->registerleaves(this);
 }

ConfigurableAnalysis::~ConfigurableAnalysis ( )

Definition at line 117 of file ConfigurableAnalysis.cc.

References selections_.

 {
   delete selections_;
 }

Member Function Documentation

void ConfigurableAnalysis::beginJob ( void )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 217 of file ConfigurableAnalysis.cc.

218 {

219 }

void ConfigurableAnalysis::endJob ( void )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 223 of file ConfigurableAnalysis.cc.

References Plotter::complete(), plotter_, Selections::print(), and selections_.

                              {
   //print summary tables
   selections_->print();
   if (plotter_) plotter_->complete();
 }

bool ConfigurableAnalysis::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDFilter.

Definition at line 128 of file ConfigurableAnalysis.cc.

References accept(), Selections::begin(), Selections::end(), Plotter::fill(), NTupler::fill(), flows_, iEvent, Filter::name(), ntupler_, plotter_, edm::Event::put(), elec_selection::selection(), selections_, Plotter::setDir(), createPayload::skip, and workAsASelector_.

 {
   using namespace edm;
 
   //will the filter pass or not.
   bool majorGlobalAccept=false;
 
   std::auto_ptr<std::vector<bool> > passedProduct(new std::vector<bool>(flows_.size(),false));
   bool filledOnce=false;  
 
   // loop the requested selections
   for (Selections::iterator selection=selections_->begin(); selection!=selections_->end();++selection){
     //was this flow of filter actually asked for
     bool skip=true;
     unsigned int iFlow=0;
     for (;iFlow!=flows_.size();++iFlow){if (flows_[iFlow]==selection->name()){skip=false; break;}}
     if (skip) continue;
 
     //make a specific direction in the plotter
     if (plotter_) plotter_->setDir(selection->name());
     
     // apply individual filters on the event
     std::map<std::string, bool> accept=selection->accept(iEvent);
     
     bool globalAccept=true;
     std::string separator="";
     std::string cumulative="";
     std::string allButOne="allBut_";
     std::string fullAccept="fullAccept";
     std::string fullContent="fullContent";
 
     if (selection->makeContentPlots() && plotter_)
       plotter_->fill(fullContent,iEvent);
 
     //loop the filters to make cumulative and allButOne job
     for (Selection::iterator filterIt=selection->begin(); filterIt!=selection->end();++filterIt){
       Filter & filter=(**filterIt);
       //      bool lastCut=((filterIt+1)==selection->end());
 
       //increment the directory name
       cumulative+=separator+filter.name(); separator="_";
 
       if (accept[filter.name()]){
     //  if (globalAccept && selection->makeCumulativePlots() && !lastCut)
     if (globalAccept && selection->makeCumulativePlots() && plotter_)
       plotter_->fill(cumulative,iEvent);
       }
       else{
     globalAccept=false;
     // did all the others filter fire
     bool goodForAllButThisOne=true;
     for (std::map<std::string,bool>::iterator decision=accept.begin(); decision!=accept.end();++decision){
       if (decision->first==filter.name()) continue;
       if (!decision->second) {
         goodForAllButThisOne=false;
         break;}
     }
     if (goodForAllButThisOne && selection->makeAllButOnePlots() && plotter_){
       plotter_->fill(allButOne+filter.name(),iEvent);
     }
       }
       
     }// loop over the filters in this selection
 
     if (globalAccept){
       (*passedProduct)[iFlow]=true;
       majorGlobalAccept=true;
       //make final plots only if no cumulative plots
       if (selection->makeFinalPlots() && !selection->makeCumulativePlots() && plotter_)
     plotter_->fill(fullAccept,iEvent);
 
       //make the ntuple and put it in the event
       if (selection->ntuplize() && !filledOnce && ntupler_){
     ntupler_->fill(iEvent);
     filledOnce=true;}
     }
     
   }//loop the different filter order/number: loop the Selections
 
   iEvent.put(passedProduct);
   if (workAsASelector_)
     return majorGlobalAccept;
   else
     return true;
 }