CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
CentralityTableProducer Class Reference
Inheritance diagram for CentralityTableProducer:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 CentralityTableProducer (const edm::ParameterSet &)
 
 ~CentralityTableProducer () override
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueueglobalLuminosityBlocksQueue ()
 
SerialTaskQueueglobalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (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 analyze (const edm::Event &, const edm::EventSetup &) override
 
void endJob () override
 
void printBin (const CentralityTable::CBin *)
 

Private Attributes

const CentralityBinsCB
 
CentralityTableCT
 
bool debug_
 
bool firstRunOnly_
 
edm::Service< TFileServicefs
 
edm::ESHandle< CentralityTableinputDB_
 
TFile * inputTFile_
 
string inputTFileName_
 
bool makeDBFromTFile_
 
bool makeTFileFromDB_
 
string rootTag_
 
unsigned int runnum_
 
ofstream text_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- 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)
 

Detailed Description

Definition at line 40 of file CentralityTableProducer.cc.

Constructor & Destructor Documentation

CentralityTableProducer::CentralityTableProducer ( const edm::ParameterSet iConfig)
explicit

Definition at line 81 of file CentralityTableProducer.cc.

References gather_cfg::cout, debug_, firstRunOnly_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), inputTFile_, inputTFileName_, makeDBFromTFile_, makeTFileFromDB_, and rootTag_.

81  : text_("bins.txt"), runnum_(0) {
82  //now do what ever initialization is needed
83  makeDBFromTFile_ = iConfig.getUntrackedParameter<bool>("makeDBFromTFile", true);
84  makeTFileFromDB_ = iConfig.getUntrackedParameter<bool>("makeTFileFromDB", false);
85  firstRunOnly_ = iConfig.getUntrackedParameter<bool>("isMC", false);
86  debug_ = iConfig.getUntrackedParameter<bool>("debug", false);
87  if (makeDBFromTFile_) {
88  inputTFileName_ = iConfig.getParameter<string>("inputTFile");
89  rootTag_ = iConfig.getParameter<string>("rootTag");
90  inputTFile_ = new TFile(inputTFileName_.data(), "read");
91  cout << inputTFileName_.data() << endl;
92  }
93 }
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
CentralityTableProducer::~CentralityTableProducer ( )
override

Definition at line 95 of file CentralityTableProducer.cc.

95  {
96  // do anything here that needs to be done at desctruction time
97  // (e.g. close files, deallocate resources etc.)
98 }

Member Function Documentation

void CentralityTableProducer::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
)
overrideprivate

Definition at line 105 of file CentralityTableProducer.cc.

References CB, gather_cfg::cout, firstRunOnly_, fs, edm::EventBase::id(), TFileDirectory::make(), makeTFileFromDB_, TFileService::mkdir(), edm::EventID::run(), runnum_, and PostProcessorHcalIsoTrack_cfi::subDir.

105  {
106  if (!makeTFileFromDB_)
107  return;
108  if ((!firstRunOnly_ && runnum_ != iEvent.id().run()) || (firstRunOnly_ && runnum_ == 0)) {
109  runnum_ = iEvent.id().run();
110  cout << "Adding table for run : " << runnum_ << endl;
111  TFileDirectory subDir = fs->mkdir(Form("run%d", runnum_));
112  CB = subDir.make<CentralityBins>();
113  }
114 }
RunNumber_t run() const
Definition: EventID.h:38
T * make(const Args &...args) const
make new ROOT object
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
edm::Service< TFileService > fs
edm::EventID id() const
Definition: EventBase.h:59
void CentralityTableProducer::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 117 of file CentralityTableProducer.cc.

References cond::service::PoolDBOutputService::appendSinceTime(), CentralityTable::CBin::b, cond::service::PoolDBOutputService::beginOfTime(), CentralityTable::CBin::bin_edge, CentralityBins::bMeanOfBin(), CentralityBins::bSigmaOfBin(), CB, gather_cfg::cout, cond::service::PoolDBOutputService::createNewIOV(), CT, cond::service::PoolDBOutputService::currentTime(), cond::service::PoolDBOutputService::endOfTime(), CentralityBins::getNbins(), inputTFile_, edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), dqmiolumiharvest::j, CentralityBins::lowEdgeOfBin(), CentralityTable::m_table, makeDBFromTFile_, CentralityTable::BinValues::mean, CentralityTable::CBin::n_coll, CentralityTable::CBin::n_hard, CentralityTable::CBin::n_part, CentralityBins::NcollMeanOfBin(), CentralityBins::NcollSigmaOfBin(), CentralityBins::NhardMeanOfBin(), CentralityBins::NhardSigmaOfBin(), CentralityBins::NpartMeanOfBin(), CentralityBins::NpartSigmaOfBin(), dqmiodatasetharvest::pool, printBin(), rootTag_, runnum_, text_, and CentralityTable::BinValues::var.

117  {
118  if (makeDBFromTFile_) {
119  runnum_ = 1;
120  // Get values from root file
121  CB = (CentralityBins*)inputTFile_->Get(Form("%s/run%d", rootTag_.data(), runnum_));
122  cout << rootTag_.data() << endl;
123  CT = new CentralityTable();
124  CT->m_table.reserve(CB->getNbins());
125 
126  text_ << "# BinEdge NpartMean NpartVar NcollMean NcollVar NhardMean NhardVar bMean bVar" << endl;
127  for (int j = 0; j < CB->getNbins(); j++) {
129  thisBin->bin_edge = CB->lowEdgeOfBin(j);
130  thisBin->n_part.mean = CB->NpartMeanOfBin(j);
131  thisBin->n_part.var = CB->NpartSigmaOfBin(j);
132  thisBin->n_coll.mean = CB->NcollMeanOfBin(j);
133  thisBin->n_coll.var = CB->NcollSigmaOfBin(j);
134  thisBin->n_hard.mean = CB->NhardMeanOfBin(j);
135  thisBin->n_hard.var = CB->NhardSigmaOfBin(j);
136  thisBin->b.mean = CB->bMeanOfBin(j);
137  thisBin->b.var = CB->bSigmaOfBin(j);
138  printBin(thisBin);
139  CT->m_table.push_back(*thisBin);
140  if (thisBin)
141  delete thisBin;
142  }
143 
145  if (pool.isAvailable()) {
146  if (pool->isNewTagRequest("HeavyIonRcd")) {
147  pool->createNewIOV<CentralityTable>(CT, pool->beginOfTime(), pool->endOfTime(), "HeavyIonRcd");
148  } else {
149  pool->appendSinceTime<CentralityTable>(CT, pool->currentTime(), "HeavyIonRcd");
150  }
151  }
152  }
153 }
void printBin(const CentralityTable::CBin *)
std::vector< CBin > m_table
float NpartSigmaOfBin(int bin) const
float lowEdgeOfBin(int bin) const
void appendSinceTime(T *payloadObj, cond::Time_t sinceTime, const std::string &recordName, bool withlogging=false)
float NcollMeanOfBin(int bin) const
bool isNewTagRequest(const std::string &recordName)
bool isAvailable() const
Definition: Service.h:40
float bSigmaOfBin(int bin) const
float NhardMeanOfBin(int bin) const
void createNewIOV(T *firstPayloadObj, cond::Time_t firstSinceTime, cond::Time_t firstTillTime, const std::string &recordName, bool withlogging=false)
float NpartMeanOfBin(int bin) const
int getNbins() const
float bMeanOfBin(int bin) const
float NhardSigmaOfBin(int bin) const
float NcollSigmaOfBin(int bin) const
void CentralityTableProducer::printBin ( const CentralityTable::CBin thisBin)
private

Definition at line 155 of file CentralityTableProducer.cc.

References CentralityTable::CBin::b, CentralityTable::CBin::bin_edge, gather_cfg::cout, DEFINE_FWK_MODULE, CentralityTable::BinValues::mean, CentralityTable::CBin::n_coll, CentralityTable::CBin::n_hard, CentralityTable::CBin::n_part, text_, and CentralityTable::BinValues::var.

Referenced by endJob().

155  {
156  cout << "HF Cut = " << thisBin->bin_edge << endl;
157  cout << "Npart = " << thisBin->n_part.mean << endl;
158  cout << "sigma = " << thisBin->n_part.var << endl;
159  cout << "Ncoll = " << thisBin->n_coll.mean << endl;
160  cout << "sigma = " << thisBin->n_coll.var << endl;
161  cout << "B = " << thisBin->b.mean << endl;
162  cout << "sigma = " << thisBin->b.var << endl;
163  text_ << Form("%0.2f\t%0.2f\t%0.2f\t%0.2f\t%0.2f\t%0.2f\t%0.2f\t%0.2f\t%0.2f",
164  (Float_t)thisBin->bin_edge,
165  (Float_t)thisBin->n_part.mean,
166  (Float_t)thisBin->n_part.var,
167  (Float_t)thisBin->n_coll.mean,
168  (Float_t)thisBin->n_coll.var,
169  (Float_t)thisBin->n_hard.mean,
170  (Float_t)thisBin->n_hard.var,
171  (Float_t)thisBin->b.mean,
172  (Float_t)thisBin->b.var)
173  << endl;
174  cout << "__________________________________________________" << endl;
175 }

Member Data Documentation

const CentralityBins* CentralityTableProducer::CB
private

Definition at line 65 of file CentralityTableProducer.cc.

Referenced by analyze(), and endJob().

CentralityTable* CentralityTableProducer::CT
private

Definition at line 64 of file CentralityTableProducer.cc.

Referenced by endJob().

bool CentralityTableProducer::debug_
private

Definition at line 54 of file CentralityTableProducer.cc.

Referenced by CentralityTableProducer().

bool CentralityTableProducer::firstRunOnly_
private

Definition at line 53 of file CentralityTableProducer.cc.

Referenced by analyze(), and CentralityTableProducer().

edm::Service<TFileService> CentralityTableProducer::fs
private

Definition at line 59 of file CentralityTableProducer.cc.

Referenced by analyze().

edm::ESHandle<CentralityTable> CentralityTableProducer::inputDB_
private

Definition at line 56 of file CentralityTableProducer.cc.

TFile* CentralityTableProducer::inputTFile_
private

Definition at line 57 of file CentralityTableProducer.cc.

Referenced by CentralityTableProducer(), and endJob().

string CentralityTableProducer::inputTFileName_
private

Definition at line 58 of file CentralityTableProducer.cc.

Referenced by CentralityTableProducer().

bool CentralityTableProducer::makeDBFromTFile_
private

Definition at line 51 of file CentralityTableProducer.cc.

Referenced by CentralityTableProducer(), and endJob().

bool CentralityTableProducer::makeTFileFromDB_
private

Definition at line 52 of file CentralityTableProducer.cc.

Referenced by analyze(), and CentralityTableProducer().

string CentralityTableProducer::rootTag_
private

Definition at line 61 of file CentralityTableProducer.cc.

Referenced by CentralityTableProducer(), and endJob().

unsigned int CentralityTableProducer::runnum_
private

Definition at line 67 of file CentralityTableProducer.cc.

Referenced by analyze(), and endJob().

ofstream CentralityTableProducer::text_
private

Definition at line 62 of file CentralityTableProducer.cc.

Referenced by endJob(), and printBin().