#include <PrepareDMRTrends.h>

Public Member Functions
void	compileDMRTrends (std::vector< int > IOVlist, TString Variable, std::vector< std::string > inputFiles, std::vector< TString > structures, const std::map< TString, int > nlayers, bool FORCE=false)

TString	getName (TString structure, int layer, TString geometry)

	PrepareDMRTrends (const char *outputFileName, boost::property_tree::ptree &json)

	~PrepareDMRTrends ()

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

const char *	outputFileName_

Detailed Description

Definition at line 148 of file PrepareDMRTrends.h.

Constructor & Destructor Documentation

◆ PrepareDMRTrends()

PrepareDMRTrends::PrepareDMRTrends	(	const char *	outputFileName,
		boost::property_tree::ptree &	json
	)

Definition at line 7 of file PrepareDMRTrends.cc.

References geometries, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                             : outputFileName_(outputFileName) {
   geometries.clear();
   for (const auto &childTree : json) {
     geometries.push_back(childTree.second.get<std::string>("title"));
   }
 }

◆ ~PrepareDMRTrends()

PrepareDMRTrends::~PrepareDMRTrends ( )

inline

Definition at line 151 of file PrepareDMRTrends.h.

151 {}

Member Function Documentation

◆ compileDMRTrends()

void PrepareDMRTrends::compileDMRTrends	(	std::vector< int >	IOVlist,
		TString	Variable,
		std::vector< std::string >	inputFiles,
		std::vector< TString >	structures,
		const std::map< TString, int >	nlayers,
		bool	FORCE = `false`
	)

Definition at line 36 of file PrepareDMRTrends.cc.

References gather_cfg::cout, Geometry::DeltaMu(), DMRFactor, f, groupFilesInBlocks::fout, g, relativeConstraints::geom, geometries, getName(), timingPdfMaker::histo, mps_fire::i, crabTemplate::inputFiles, utils::is_empty(), nano_mu_digi_cff::layer, Geometry::Mu(), Geometry::MuMinus(), Geometry::MuPlus(), dqmiodumpmetadata::n, Skims_PA_cff::name, nlayers, outputFileName_, point, writedatasetfile::runs, Geometry::Sigma(), Geometry::SigmaDeltaMu(), Geometry::SigmaMinus(), Geometry::SigmaPlus(), trends(), and hcalRecHitTable_cff::variables.

Referenced by trends().

                                                     {
   gROOT->SetBatch();
 
   float ScaleFactor = DMRFactor;
   if (Variable == "DrmsNR")
     ScaleFactor = 1;
 
   map<pair<pair<TString, int>, TString>, Geometry> mappoints;  // pair = (structure, layer), geometry
   Point *point = nullptr;
   TFile *f = nullptr;
 
   for (unsigned int i = 0; i < inputFiles.size(); ++i) {
     if (fs::is_empty(inputFiles.at(i).c_str())) {
       cout << "ERROR: Empty file " << inputFiles.at(i).c_str() << endl;
       continue;
     }
 
     int runN = IOVlist.at(i);
 
     f = new TFile(inputFiles.at(i).c_str(), "READ");
     std::cout << inputFiles.at(i) << std::endl;
 
     for (TString &structure : structures) {
       TString structname = structure;
       structname.ReplaceAll("_y", "");
       size_t layersnumber = nlayers.at(structname);
       for (size_t layer = 0; layer <= layersnumber; layer++) {
         for (const string &geometry : geometries) {
           TString name = Variable + "_" + getName(structure, layer, geometry);
           TH1F *histo = dynamic_cast<TH1F *>(f->Get(name));
 
           // Three possibilities:
           //  - All histograms are produced correctly
           //  - Only the non-split histograms are produced
           //  - No histogram is produced correctly
           //  FORCE means that the Point is not added to the points collection in the chosen geometry for that structure
           //  If FORCE is not enabled a default value for the Point is used (-9999) which will appear in the plots
           if (!histo) {
             //cout << "Run" << runN << " Histogram: " << name << " not found" << endl;
             if (FORCE)
               continue;
             point = new Point(runN, ScaleFactor);
           } else if (structure != "TID" && structure != "TEC") {
             TH1F *histoplus = dynamic_cast<TH1F *>(f->Get((name + "_plus")));
             TH1F *histominus = dynamic_cast<TH1F *>(f->Get((name + "_minus")));
             if (!histoplus || !histominus) {
               //cout << "Run" << runN << " Histogram: " << name << " plus or minus not found" << endl;
               if (FORCE)
                 continue;
               point = new Point(runN, ScaleFactor, histo);
             } else
               point = new Point(runN, ScaleFactor, histo, histoplus, histominus);
 
           } else
             point = new Point(runN, ScaleFactor, histo);
           mappoints[make_pair(make_pair(structure, layer), geometry)].points.push_back(*point);
         }
       }
     }
     f->Close();
   }
 
   TFile *fout = TFile::Open(outputFileName_, "RECREATE");
   TGraphErrors *g = nullptr;
   for (TString &structure : structures) {
     TString structname = structure;
     structname.ReplaceAll("_y", "");
     size_t layersnumber = nlayers.at(structname);
     for (size_t layer = 0; layer <= layersnumber; layer++) {
       for (const string &geometry : geometries) {
         TString name = Variable + "_" + getName(structure, layer, geometry);
         Geometry geom = mappoints[make_pair(make_pair(structure, layer), geometry)];
         using Trend = vector<float> (Geometry::*)() const;
         vector<Trend> trends{&Geometry::Mu,
                              &Geometry::Sigma,
                              &Geometry::MuPlus,
                              &Geometry::SigmaPlus,
                              &Geometry::MuMinus,
                              &Geometry::SigmaMinus,
                              &Geometry::DeltaMu,
                              &Geometry::SigmaDeltaMu};
         vector<TString> variables{
             "mu", "sigma", "muplus", "sigmaplus", "muminus", "sigmaminus", "deltamu", "sigmadeltamu"};
         vector<float> runs = geom.Run();
         size_t n = runs.size();
         vector<float> emptyvec;
         for (size_t i = 0; i < runs.size(); i++)
           emptyvec.push_back(0.);
         for (size_t iVar = 0; iVar < variables.size(); iVar++) {
           Trend trend = trends.at(iVar);
           g = new TGraphErrors(n, runs.data(), (geom.*trend)().data(), emptyvec.data(), emptyvec.data());
           g->SetTitle(geometry.c_str());
           g->Write(name + "_" + variables.at(iVar));
         }
         vector<pair<Trend, Trend>> trendspair{make_pair(&Geometry::Mu, &Geometry::Sigma),
                                               make_pair(&Geometry::MuPlus, &Geometry::SigmaPlus),
                                               make_pair(&Geometry::MuMinus, &Geometry::SigmaMinus),
                                               make_pair(&Geometry::DeltaMu, &Geometry::SigmaDeltaMu)};
         vector<pair<TString, TString>> variablepairs{make_pair("mu", "sigma"),
                                                      make_pair("muplus", "sigmaplus"),
                                                      make_pair("muminus", "sigmaminus"),
                                                      make_pair("deltamu", "sigmadeltamu")};
         for (size_t iVar = 0; iVar < variablepairs.size(); iVar++) {
           Trend meantrend = trendspair.at(iVar).first;
           Trend sigmatrend = trendspair.at(iVar).second;
           g = new TGraphErrors(
               n, runs.data(), (geom.*meantrend)().data(), emptyvec.data(), (geom.*sigmatrend)().data());
           g->SetTitle(geometry.c_str());
           TString graphname = name + "_" + variablepairs.at(iVar).first;
           graphname += variablepairs.at(iVar).second;
           g->Write(graphname);
         }
       }
     }
   }
   fout->Close();
 
   delete point;
   delete f;
   delete g;
 }

◆ getName()

TString PrepareDMRTrends::getName	(	TString	structure,
		int	layer,
		TString	geometry
	)

Definition at line 18 of file PrepareDMRTrends.cc.

References nano_mu_digi_cff::layer, and Skims_PA_cff::name.

Referenced by compileDMRTrends(), and plotting.Plot::draw().

                                                                                 {
   geometry.ReplaceAll(" ", "_");
   TString name = geometry + "_" + structure;
   if (layer != 0) {
     if (structure == "TID" || structure == "TEC")
       name += "_disc";
     else
       name += "_layer";
     name += layer;
   }
 
   return name;
 };

Member Data Documentation

◆ geometries

std::vector<std::string> PrepareDMRTrends::geometries

private

Definition at line 163 of file PrepareDMRTrends.h.

Referenced by compileDMRTrends(), and PrepareDMRTrends().

◆ outputFileName_

const char* PrepareDMRTrends::outputFileName_