#include <SusyPostProcessor.h>

Inheritance diagram for SusyPostProcessor:

Public Member Functions
	SusyPostProcessor (const edm::ParameterSet &pSet)

	~SusyPostProcessor () override

Public Member Functions inherited from DQMEDHarvester
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	beginRun (edm::Run const &, edm::EventSetup const &) override

	DQMEDHarvester ()

virtual void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)

void	endJob () final

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &, edm::EventSetup const &) final

void	endRun (edm::Run const &, edm::EventSetup const &) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDHarvester () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

void	QuantilePlots (MonitorElement *&, double, DQMStore::IBooker &)

Private Attributes
double	_quantile

std::vector< std::string >	Dirs

std::vector< MonitorElement * >	histoVector

edm::ParameterSet	iConfig

MonitorElement *	MEx

MonitorElement *	MEy

std::string	SUSYFolder

Static Private Attributes
static const char *	messageLoggerCatregory = "SusyDQMPostProcessor"

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from DQMEDHarvester
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 22 of file SusyPostProcessor.h.

Constructor & Destructor Documentation

SusyPostProcessor::SusyPostProcessor ( const edm::ParameterSet & pSet )

explicit

Definition at line 12 of file SusyPostProcessor.cc.

References edm::ParameterSet::getParameter().

 {
   iConfig = pSet;
 
   SUSYFolder = iConfig.getParameter<string>("folderName");
   _quantile = iConfig.getParameter<double>("quantile");
 
 }

SusyPostProcessor::~SusyPostProcessor ( )

override

Definition at line 21 of file SusyPostProcessor.cc.

21 {

22 }

Member Function Documentation

void SusyPostProcessor::dqmEndJob	(	DQMStore::IBooker &	ibook_,
		DQMStore::IGetter &	iget_
	)

overrideprivatevirtual

Implements DQMEDHarvester.

Definition at line 57 of file SusyPostProcessor.cc.

References DQMStore::IGetter::dirExists(), TrackerOfflineValidation_Dqm_cff::dirName, MonitorElement::DQM_KIND_TH1F, spr::find(), runEdmFileComparison::found, DQMStore::IGetter::get(), DQMStore::IGetter::getContents(), DQMStore::IGetter::getSubdirs(), MonitorElement::getTH1F(), mps_fire::i, createfilelist::int, tablePrinter::prefix, DQMStore::IGetter::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   // MET
   //----------------------------------------------------------------------------
   iget_.setCurrentFolder("JetMET/MET");
 
   Dirs = iget_.getSubdirs();
 
   std::vector<std::string> metFolders;
 
   metFolders.push_back("Uncleaned/");
   metFolders.push_back("Cleaned/");
 
   //Need our own copy for thread safety
   TF1 mygaus("mygaus","gaus");
 
   for (int i=0; i<int(Dirs.size()); i++) {
 
     std::string prefix = "dummy";
 
     if (size_t(Dirs[i].find("met")) != string::npos) prefix = "met";
     if (size_t(Dirs[i].find("pfMet"))   != string::npos) prefix = "pfMET";
 
     for (std::vector<std::string>::const_iterator ic=metFolders.begin();
      ic!=metFolders.end(); ic++) {
 
       std::string dirName = Dirs[i] +"/" + *ic;
 
       MEx = iget_.get(dirName + "/"+"MEx");
       MEy = iget_.get(dirName + "/"+"MEy");
 
       if (MEx && MEx->kind() == MonitorElement::DQM_KIND_TH1F) {
     if (MEx->getTH1F()->GetEntries() > 50) MEx->getTH1F()->Fit(&mygaus, "q");
       }
 
       if (MEy && MEy->kind() == MonitorElement::DQM_KIND_TH1F) {
     if (MEy->getTH1F()->GetEntries() > 50) MEy->getTH1F()->Fit(&mygaus, "q");
       }
     }
   }
 
 
   // SUSY
   //----------------------------------------------------------------------------
   iget_.setCurrentFolder(SUSYFolder);
   Dirs = iget_.getSubdirs();
   for (int i=0; i<int(Dirs.size()); i++)
     {
       size_t found = Dirs[i].find("Alpha");
       if (found!=string::npos) continue;
       if(!iget_.dirExists(Dirs[i])){
     edm::LogError(messageLoggerCatregory)<< "Directory "<<Dirs[i]<<" doesn't exist!!";
     continue;
       }      
       vector<MonitorElement*> histoVector = iget_.getContents(Dirs[i]);
       for (int i=0; i<int(histoVector.size()); i++) {
     QuantilePlots(histoVector[i],_quantile,ibook_);
       } 
     }
 }

void SusyPostProcessor::QuantilePlots	(	MonitorElement *&	ME,
		double	q_value,
		DQMStore::IBooker &	ibooker_
	)

private

Definition at line 26 of file SusyPostProcessor.cc.

References DQMStore::IBooker::book1D(), MonitorElement::Fill(), MonitorElement::getPathname(), MonitorElement::getTH1(), createfilelist::int, SiStripPI::mean, dataset::name, lumiQueryAPI::q, RMS, and DQMStore::IBooker::setCurrentFolder().

 {
   if(ME->getTH1()->GetEntries()>0.)
     {
       Quantile q(static_cast<const TH1*>(ME->getTH1()));
       Float_t mean=q[q_value].first;
       Float_t RMS=q[q_value].second;
       
       Float_t xLow=-5.5;
       Float_t xUp=9.5;
       Int_t NBin=15;
       
       if(mean>0.)
     {
       Float_t DBin=RMS*TMath::Sqrt(12.)/2.;
       xLow=mean-int(mean/DBin+2)*DBin;
       xUp=int(0.2*mean/DBin)*DBin+mean+5*DBin;
       NBin=(xUp-xLow)/DBin;
     }
 
       ibooker_.setCurrentFolder(ME->getPathname());
       TString name=ME->getTH1()->GetName();
       name+="_quant";
       ME=ibooker_.book1D(name,"",NBin, xLow, xUp);
       ME->Fill(mean-RMS);
       ME->Fill(mean+RMS);
     }
 }

Member Data Documentation

double SusyPostProcessor::_quantile

private

Definition at line 38 of file SusyPostProcessor.h.

std::vector<std::string> SusyPostProcessor::Dirs

private

Definition at line 41 of file SusyPostProcessor.h.

std::vector<MonitorElement*> SusyPostProcessor::histoVector

private

Definition at line 40 of file SusyPostProcessor.h.

edm::ParameterSet SusyPostProcessor::iConfig

private

Definition at line 31 of file SusyPostProcessor.h.

const char * SusyPostProcessor::messageLoggerCatregory = "SusyDQMPostProcessor"

staticprivate

Definition at line 35 of file SusyPostProcessor.h.

MonitorElement* SusyPostProcessor::MEx

private

Definition at line 43 of file SusyPostProcessor.h.

MonitorElement* SusyPostProcessor::MEy

private

Definition at line 44 of file SusyPostProcessor.h.

std::string SusyPostProcessor::SUSYFolder

private

Definition at line 37 of file SusyPostProcessor.h.

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation