DMRtrends.cc File Reference

#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
#include <map>
#include <iomanip>
#include <experimental/filesystem>
#include "TPad.h"
#include "TCanvas.h"
#include "TGraph.h"
#include "TGraphErrors.h"
#include "TMultiGraph.h"
#include "TH1.h"
#include "THStack.h"
#include "TROOT.h"
#include "TFile.h"
#include "TColor.h"
#include "TLegend.h"
#include "TLegendEntry.h"
#include "TMath.h"
#include "TRegexp.h"
#include "TPaveLabel.h"
#include "TPaveText.h"
#include "TStyle.h"
#include "TLine.h"

Go to the source code of this file.

Classes
class	Geometry
	Class Geometry Contains vector for fit parameters (mean, sigma, etc.) obtained from multiple IOVs See Structure Point for description of the parameters. More...
struct	Point
	Structure Point Contains parameters of Gaussian fits to DMRs. More...

Macros
#define	DMRFactor   10000.
#define	DUMMY   -999.
#define	lumiFactor   1000.

Functions
bool	checkrunlist (vector< int > runs, vector< int > IOVlist={}, TString Year="2018")
void	compileDMRTrends (vector< int > IOVlist, vector< string >labels={"MB"}, TString Year="2018", string myValidation="/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/results/acardini/DMRs/", vector< string > geometries={"GT","SG","MP pix LBL","PIX HLS+ML STR fix"}, bool showlumi=false, bool FORCE=false)
TH1F *	ConvertToHist (TGraphErrors *g)
void	DMRtrends (vector< int > IOVlist, vector< string >labels={"MB"}, TString Year="2018", string myValidation="/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/results/acardini/DMRs/", vector< string > geometries={"GT","SG","MP pix LBL","PIX HLS+ML STR fix"}, vector< Color_t > colours={kBlue, kRed, kGreen, kCyan}, TString outputdir="/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/alignmentObjects/acardini/DMRsTrends/", bool pixelupdate=false, vector< int > pixelupdateruns={314881, 316758, 317527, 318228, 320377}, bool showlumi=false, bool FORCE=false)
double	getintegratedlumiuptorun (int run, TString Year="2018", double min=0.)
TString	getName (TString structure, int layer, TString geometry)
TString	lumifileperyear (TString Year="2018", string RunOrIOV="IOV")
int	main (int argc, char *argv[])
void	PixelUpdateLines (TCanvas *c, TString Year="2018", bool showlumi=false, vector< int >pixelupdateruns={314881, 316758, 317527, 318228, 320377})
void	PlotDMRTrends (vector< string >labels={"MB"}, TString Year="2018", string myValidation="/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/results/acardini/DMRs/", vector< string > geometries={"GT","SG","MP pix LBL","PIX HLS+ML STR fix"}, vector< Color_t > colours={kBlue, kRed, kGreen, kCyan}, TString outputdir="/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/alignmentObjects/acardini/DMRsTrends/", bool pixelupdate=false, vector< int > pixelupdateruns={314881, 316758, 317527, 318228, 320377}, bool showlumi=false)
vector< int >	runlistfromlumifile (TString Year="2018")
void	scalebylumi (TGraphErrors *g, TString Year="2018", double min=0.)

Macro Definition Documentation

#define DMRFactor   10000.

Definition at line 41 of file DMRtrends.cc.

Referenced by Point::GetDeltaMu(), Point::GetMu(), Point::GetMuMinus(), Point::GetMuPlus(), Point::GetSigma(), Point::GetSigmaDeltaMu(), Point::GetSigmaMinus(), and Point::GetSigmaPlus().

#define DUMMY   -999.

Definition at line 33 of file DMRtrends.cc.

Referenced by Point::GetDeltaMu(), Point::GetSigmaDeltaMu(), and scalebylumi().

#define lumiFactor   1000.

Definition at line 37 of file DMRtrends.cc.

Referenced by getintegratedlumiuptorun(), and scalebylumi().

Function Documentation

bool checkrunlist	(	vector< int >	runs,
		vector< int >	IOVlist = {},
		TString	Year = "2018"
	)

Definition at line 211 of file DMRtrends.cc.

References gather_cfg::cout, spr::find(), writedatasetfile::run, lumiPlot::runlist, and runlistfromlumifile().

Referenced by compileDMRTrends(), Point::GetSigmaDeltaMu(), and PlotDMRTrends().

                                                                     {
    vector<int> runlist = runlistfromlumifile(Year);
    vector<int> missingruns; //runs for which the luminosity is not found
    vector<int> lostruns; //IOVs for which the DMR were not found
    bool problemfound=false;
    for(int run : runs){
        if(find(runlist.begin(),runlist.end(),run)==runlist.end()){
        problemfound=true;
        missingruns.push_back(run);
    }
    }
    if(!IOVlist.empty()) for(int IOV : IOVlist){
        if(find(runs.begin(),runs.end(),IOV)==runs.end()){
        problemfound=true;
        lostruns.push_back(IOV);
    }
    }
    std::sort(missingruns.begin(),missingruns.end());
    if(problemfound){
      if(!lostruns.empty()){ 
        cout << "WARNING: some IOVs where not found among the list of available DMRs" << endl << "List of missing IOVs:" << endl;
        for (int lostrun : lostruns) cout << to_string(lostrun) << " ";
        cout << endl;
    }
      if(!missingruns.empty()){
        cout << "WARNING: some runs are missing in the run/luminosity txt file" << endl << "List of missing runs:" << endl;
        for (int missingrun : missingruns) cout << to_string(missingrun) << " ";
        cout << endl;
    }
    }
    return problemfound;

} 

void compileDMRTrends	(	vector< int >	IOVlist,
		vector< string >	labels = {"MB"},
		TString	Year = "2018",
		string	myValidation = "/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/results/acardini/DMRs/",
		vector< string >	geometries = {"GT","SG", "MP pix LBL","PIX HLS+ML STR fix"},
		bool	showlumi = false,
		bool	FORCE = false
	)

Definition at line 260 of file DMRtrends.cc.

References checkrunlist(), gather_cfg::cout, Geometry::DeltaMu(), mps_splice::entry, cmsRelvalreport::exit, f, corrVsCorr::filename, cmsPerfCommons::filenames, groupFilesInBlocks::fout, g, relativeConstraints::geom, getName(), trackerHits::histo, mps_fire::i, utils::is_empty(), label, Geometry::Mu(), Geometry::MuMinus(), Geometry::MuPlus(), gen::n, dataset::name, nlayers, point, fetchall_from_DQM_v2::regexp, Geometry::Run(), writedatasetfile::runs, Geometry::Sigma(), Geometry::SigmaDeltaMu(), Geometry::SigmaMinus(), Geometry::SigmaPlus(), and objects.autophobj::variables.

Referenced by DMRtrends(), and Point::GetSigmaDeltaMu().

                                                                                                                                                          {
    gROOT->SetBatch();
    vector<int>RunNumbers;
    vector<TString> filenames;
    TRegexp regexp("[0-9][0-9][0-9][0-9][0-9][0-9]");
    for (const auto & entry : fs::recursive_directory_iterator(myValidation)){
        bool found_all_labels = true;
    for(string label : labels){
        if(entry.path().string().find(label)==std::string::npos) found_all_labels = false;
    }
    if ((entry.path().string().find("ExtendedOfflineValidation_Images/OfflineValidationSummary.root")!=std::string::npos)&&found_all_labels){
        if(fs::is_empty(entry.path())) cout << "ERROR: Empty file " << entry.path() << endl;
        else{
            TString filename(entry.path().string());
        filenames.push_back(filename);
        TString runstring(filename(regexp));
        if(runstring.IsFloat()){
            int runN=runstring.Atoi();
            RunNumbers.push_back(runN);
        }
        }
    }
    }
    if(checkrunlist(RunNumbers,IOVlist,Year)&&showlumi&&!FORCE){
        cout << "Please fix the run/luminosities file!" << endl;
    exit(EXIT_FAILURE);
    }

    vector<TString> structures { "BPIX", "BPIX_y", "FPIX", "FPIX_y", "TIB", "TID", "TOB", "TEC"};

    const map<TString,int> nlayers{ {"BPIX", 4}, {"FPIX", 3}, {"TIB", 4}, {"TID", 3}, {"TOB", 6}, {"TEC", 9} };


    map<pair<pair<TString,int>,TString>,Geometry> mappoints; // pair = (structure, layer), geometry
    
    std::sort(filenames.begin(),filenames.end());//order the files in alphabetical order
    for (TString filename: filenames){
        int runN;
        TString runstring(filename(regexp));
        if(runstring.IsFloat()){
      runN=runstring.Atoi();
    }
        else{
      cout << "ERROR: run number not retrieved for file " << filename << endl;
      continue;
    }

        TFile * f = new TFile(filename, "READ");

        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 (string geometry: geometries) {
                    TString name = getName(structure, layer, geometry);
                    TH1F *histo      = dynamic_cast<TH1F*>(f->Get( name));
                    //Geometry *geom =nullptr;
                    Point * point = nullptr;
            // 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);
                    }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, histo);
                        }else point= new Point(runN, histo, histoplus, histominus);

                    }else point= new Point(runN, histo);
                    mappoints[make_pair(make_pair(structure,layer),geometry)].points.push_back(*point);
                }
            }
        }
        f->Close();
    }
    TString outname=myValidation+"DMRtrends";
    for(TString label : labels){outname+="_"; outname+=label;}
    outname+=".root";
    TFile * fout = TFile::Open(outname, "RECREATE");
    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 (string geometry : geometries) {
                TString name = 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);
                    TGraphErrors *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;
                    TGraphErrors *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();

}

TH1F * ConvertToHist

(

TGraphErrors *

g

)

Definition at line 526 of file DMRtrends.cc.

References create_public_pileup_plots::bins, trackerHits::histo, mps_fire::i, N, and funct::pow().

Referenced by Point::GetSigmaDeltaMu(), and PlotDMRTrends().

                                    { 
    size_t N=g->GetN();
    double* x=g->GetX();
    double* y=g->GetY();
    double* xerr=g->GetEX();
    vector<float> bins;
    bins.push_back(x[0]-xerr[0]);
    for(size_t i=1;i<N;i++){
      if((x[i-1]+xerr[i-1]) > (bins.back() + pow(10,-6))) bins.push_back(x[i-1]+xerr[i-1]);
      if((x[i]-xerr[i])     > (bins.back() + pow(10,-6))) bins.push_back(x[i]-xerr[i]);
      
    }
    bins.push_back(x[N-1]+xerr[N-1]);
    TString histoname="histo_";
    histoname+=g->GetName();
    TH1F *histo = new TH1F(histoname,g->GetTitle(),bins.size()-1,bins.data());
    for(size_t i=0;i<N;i++){
      histo->Fill(x[i],y[i]);
      histo->SetBinError(histo->FindBin(x[i]),pow(10,-6));
    }
  return histo;
}

void DMRtrends	(	vector< int >	IOVlist,
		vector< string >	labels = {"MB"},
		TString	Year = "2018",
		string	myValidation = "/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/results/acardini/DMRs/",
		vector< string >	geometries = {"GT","SG", "MP pix LBL","PIX HLS+ML STR fix"},
		vector< Color_t >	colours = {kBlue, kRed, kGreen, kCyan},
		TString	outputdir = "/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/alignmentObjects/acardini/DMRsTrends/",
		bool	pixelupdate = false,
		vector< int >	pixelupdateruns = {314881, 316758, 317527, 318228, 320377},
		bool	showlumi = false,
		bool	FORCE = false
	)

Definition at line 249 of file DMRtrends.cc.

References compileDMRTrends(), gather_cfg::cout, and PlotDMRTrends().

Referenced by Point::GetSigmaDeltaMu(), and main().

                                                                                                                                                                                                                                             {
  compileDMRTrends(IOVlist, labels, Year, myValidation, geometries, showlumi, FORCE);
    cout<< "Begin plotting"<<endl;
    PlotDMRTrends(labels, Year, myValidation, geometries, colours, outputdir, pixelupdate, pixelupdateruns, showlumi);

};

double getintegratedlumiuptorun	(	int	run,
		TString	Year = "2018",
		double	min = 0.
	)

Definition at line 451 of file DMRtrends.cc.

References csvLumiCalc::lumi, lumiFactor, lumifileperyear(), min(), Scenarios_cff::scale, and matplotRender::yscale.

Referenced by Point::GetSigmaDeltaMu(), and PixelUpdateLines().

                                                                  {
    TGraph * scale = new TGraph((lumifileperyear(Year,"run")).Data());
    int Nscale=scale->GetN();
    double *xscale=scale->GetX();
    double *yscale=scale->GetY();


    double lumi=min;
    int index=-1;
    for(int j=0;j<Nscale;j++){
        lumi+=yscale[j];
        if(run>=xscale[j]){
            index=j;
            continue;
        }
    }
    lumi=min;  
    for(int j=0;j<index;j++)lumi+=yscale[j]/lumiFactor;

    return lumi;

}

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

Definition at line 167 of file DMRtrends.cc.

References dataset::name.

Referenced by MuonAlignment::analyze(), CaloTowersClient::CaloTowersEndjob(), CaloTowersDQMClient::CaloTowersEndjob(), ESDaqInfoTask::cleanup(), jsoncollector::FastMonitor::commit(), compileDMRTrends(), FWGUIManager::createHelpGLPopup(), FWGUIManager::createShortcutPopup(), HGCalDigiClient::digisEndjob(), JetMETHLTOfflineClient::dqmEndJob(), MonitorElementsDb::endJob(), HGCalGeometryClient::geometryEndjob(), HGCalHitClient::geometryEndjob(), cscdqm::HistoDef::getFullPath(), Point::GetSigmaDeltaMu(), HcalDigisClient::HcalDigisEndjob(), HcalRecHitsDQMClient::HcalRecHitsEndjob(), NoiseRatesClient::NoiseRatesEndjob(), HcalNoiseRatesClient::NoiseRatesEndjob(), GenericMVAComputer::TaggingVariableIterator< Iter_t >::Value::operator PhysicsTools::Variable::Value(), PlotDMRTrends(), HGCalRecHitsClient::recHitsEndjob(), HcalSimHitsClient::SimHitsEndjob(), HGCalSimHitsClient::simHitsEndjob(), and PhysicsTools::TrainerMonitoring::RootObject< T >::write().

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

    return name;
};

TString lumifileperyear	(	TString	Year = "2018",
		string	RunOrIOV = "IOV"
	)

Definition at line 183 of file DMRtrends.cc.

References gather_cfg::cout.

Referenced by getintegratedlumiuptorun(), Point::GetSigmaDeltaMu(), runlistfromlumifile(), and scalebylumi().

                                                      {
  TString LumiFile=getenv("CMSSW_BASE"); LumiFile +="/src/Alignment/OfflineValidation/data/lumiper";
  if(RunOrIOV!="run"&&RunOrIOV!="IOV") cout << "ERROR: Please specify \"run\" or \"IOV\" to retrieve the luminosity run by run or for each IOV"<<endl;
  LumiFile+=RunOrIOV;
  if(Year!="2017"&&Year!="2018") cout << "ERROR: Only 2017 and 2018 lumi-per-run/IOV are available, please check!"<<endl;
  LumiFile+=Year;
  LumiFile+=".txt";
  return LumiFile;

};

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 701 of file DMRtrends.cc.

References gather_cfg::cout, DMRtrends(), mps_fire::i, tablePrinter::labels, and lumiPlot::runlist.

                                   { 
    if (argc == 1) {

            vector<int>IOVlist={290543, 296702, 296966, 297224, 297281, 297429, 297467, 297484, 297494, 297503, 297557, 297599, 297620, 297660, 297670, 298678, 298996, 299062, 299096, 299184, 299327, 299368, 299381, 299443, 299480, 299592, 299594, 299649, 300087, 300155, 300233, 300237, 300280, 300364, 300389, 300399, 300459, 300497, 300515, 300538, 300551, 300574, 300636, 300673, 300780, 300806, 300812, 301046, 301417, 302131, 302573, 302635, 303825, 303998, 304170, 304505, 304672, 305040, 305081};//UL17
            //vector<int> pixelupdateruns {316758, 317527,317661,317664,318227, 320377};//2018
        vector<int> pixelupdateruns {290543, 297281, 298653, 299443, 300389, 302131, 303790, 304911, 305898};//2017

            cout << "WARNING: Running function with arguments specified in DMRtrends.cc" << endl << "If you want to specify the arguments from command line run the macro as follows:" << endl << "DMRtrends labels pathtoDMRs geometriesandcolourspairs outputdirectory showpixelupdate showlumi FORCE" << endl;

        //Example provided for a currently working set of parameters
        DMRtrends(IOVlist,{"vUL17","MB"},"2017", "/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/results/acardini/DMRs/EOY17/", {"EOY17","full ML pixel + strip","SG-mp2607","mp2993"}, {kRed, kBlack, kBlue, kGreen}, "/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/alignmentObjects/acardini/DMRsTrends/", true, pixelupdateruns, true, true); 
        
        
        return 0;
    }
    else if (argc < 11) {
        cout << "DMRtrends IOVlist labels Year pathtoDMRs geometriesandcolourspairs outputdirectory pixelupdatelist showpixelupdate showlumi FORCE" << endl;
        
        return 1;
    }

    TString runlist = argv[1],
            all_labels = argv[2],
            Year = argv[3],
        pathtoDMRs = argv[4],
        geometrieandcolours = argv[5], //name1:title1:color1,name2:title2:color2,name3:title3:color3
            outputdirectory = argv[6],
            pixelupdatelist = argv[7];
    bool    showpixelupdate = argv[8],
        showlumi = argv[9],
        FORCE = argv[10];
    TObjArray *labelarray = all_labels.Tokenize(",");
    vector<string> labels;
    for(int i=0; i < labelarray->GetEntries(); i++)labels.push_back((string)(labelarray->At(i)->GetName()));
    TObjArray *IOVarray = runlist.Tokenize(",");
    vector<int> IOVlist;
    for(int i=0; i < IOVarray->GetEntries(); i++)IOVlist.push_back(stoi(IOVarray->At(i)->GetName()));
    vector<int> pixelupdateruns;
    TObjArray *PIXarray = pixelupdatelist.Tokenize(",");
    for(int i=0; i < PIXarray->GetEntries(); i++)pixelupdateruns.push_back(stoi(PIXarray->At(i)->GetName()));
    vector<string> geometries;
    vector<Color_t> colours;
    TObjArray *geometrieandcolourspairs = geometrieandcolours.Tokenize(",");
    for (int i=0; i < geometrieandcolourspairs->GetEntries(); i++) {
            TObjArray *geomandcolourvec = TString(geometrieandcolourspairs->At(i)->GetName()).Tokenize(":");
        geometries.push_back(geomandcolourvec->At(0)->GetName());
        colours.push_back((Color_t)(atoi(geomandcolourvec->At(1)->GetName())));
    }
    DMRtrends(IOVlist,labels,Year,pathtoDMRs.Data(),geometries,colours,outputdirectory.Data(),showpixelupdate, pixelupdateruns,showlumi,FORCE);

    
    return 0; 
}

void PixelUpdateLines	(	TCanvas *	c,
		TString	Year = "2018",
		bool	showlumi = false,
		vector< int >	pixelupdateruns = {314881, 316758, 317527, 318228, 320377}
	)

Definition at line 393 of file DMRtrends.cc.

References getintegratedlumiuptorun(), label, tablePrinter::labels, mps_splice::line, csvLumiCalc::lumi, Min(), and corrVsCorr::rx.

Referenced by Point::GetSigmaDeltaMu(), and PlotDMRTrends().

                                                                                          {
    vector<TPaveText*> labels;
    double lastlumi=0.;
    c->cd();
    size_t index = 0;
    for(int pixelupdaterun : pixelupdateruns){
           double lumi=0.;
           if(showlumi)lumi=getintegratedlumiuptorun(pixelupdaterun,Year); //The vertical line needs to be drawn at the beginning of the run where the pixel update was implemented, thus only the integrated luminosity up to that run is required.
           else lumi=pixelupdaterun;
           TLine *line = new TLine (lumi,c->GetUymin(),lumi,c->GetUymax());
           line->SetLineColor(kBlue);
           line->SetLineStyle(9);
           line->Draw();
           //Due to the way the coordinates within the Canvas are set, the following steps are required to draw the TPaveText:
           // Compute the gPad coordinates in TRUE normalized space (NDC)
           int ix1;
           int ix2;
           int iw = gPad->GetWw();
           int ih = gPad->GetWh();
           double x1p,y1p,x2p,y2p;
           gPad->GetPadPar(x1p,y1p,x2p,y2p);
           ix1 = (Int_t)(iw*x1p);
           ix2 = (Int_t)(iw*x2p);
           double wndc  = TMath::Min(1.,(double)iw/(double)ih);
           double rw    = wndc/(double)iw;
           double x1ndc = (double)ix1*rw;
           double x2ndc = (double)ix2*rw;
           // Ratios to convert user space in TRUE normalized space (NDC)
           double rx1,ry1,rx2,ry2;
           gPad->GetRange(rx1,ry1,rx2,ry2);
           double rx = (x2ndc-x1ndc)/(rx2-rx1);
           double _sx;
           // Left limit of the TPaveText
           _sx = rx*(lumi-rx1)+x1ndc;
           // To avoid an overlap between the TPaveText a vertical shift is done when the IOVs are too close
           if( _sx < lastlumi ){
         index++;
           }else index=0;
           TPaveText *box= new TPaveText(_sx+0.0015,0.86-index*0.04,_sx+0.035,0.89-index*0.04,"blNDC");
           box->SetFillColor(10);
           box->SetBorderSize(1);
           box->SetLineColor(kBlack);
           TText *textRun = box->AddText(Form("%i",int(pixelupdaterun)));
           textRun->SetTextSize(0.025);
           labels.push_back(box);
           lastlumi=_sx+0.035;
    }
    //Drawing in a separate loop to ensure that the labels are drawn on top of the lines
    for(auto label: labels){
           label->Draw("same");
    }
    c->Update();
}

void PlotDMRTrends	(	vector< string >	labels = {"MB"},
		TString	Year = "2018",
		string	myValidation = "/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/results/acardini/DMRs/",
		vector< string >	geometries = {"GT","SG", "MP pix LBL","PIX HLS+ML STR fix"},
		vector< Color_t >	colours = {kBlue, kRed, kGreen, kCyan},
		TString	outputdir = "/afs/cern.ch/cms/CAF/CMSALCA/ALCA_TRACKERALIGN/data/commonValidation/alignmentObjects/acardini/DMRsTrends/",
		bool	pixelupdate = false,
		vector< int >	pixelupdateruns = {314881, 316758, 317527, 318228, 320377},
		bool	showlumi = false
	)

Definition at line 554 of file DMRtrends.cc.

References EnergyCorrector::c, checkrunlist(), ConvertToHist(), corrVsCorr::filename, g, utils::gErrorIgnoreLevel, getName(), h, mps_fire::i, recoMuon::in, fwlog::kWarning, label, listHistos::legend, SiStripPI::max, min(), dataset::name, nlayers, PixelUpdateLines(), scalebylumi(), taus_updatedMVAIds_cff::variable, and objects.autophobj::variables.

Referenced by DMRtrends(), and Point::GetSigmaDeltaMu().

                                                                                                                                                                                                                 {
    gErrorIgnoreLevel = kWarning;
    checkrunlist(pixelupdateruns,{},Year);
    vector<TString> structures { "BPIX", "BPIX_y", "FPIX", "FPIX_y", "TIB", "TID", "TOB", "TEC"};

    const map<TString,int> nlayers{ {"BPIX", 4}, {"FPIX", 3}, {"TIB", 4}, {"TID", 3}, {"TOB", 6}, {"TEC", 9} };


    TString filename=myValidation+"DMRtrends";
    for(TString label : labels){ filename+="_"; filename+=label;}
    filename+=".root";
    TFile *in= new TFile(filename);
    for (TString& structure: structures) {
        TString structname = structure;
        structname.ReplaceAll("_y", "");
        int layersnumber=nlayers.at(structname);
        for (int layer=0; layer<=layersnumber;layer++){
            vector<TString> variables {"mu", "sigma", "muplus", "sigmaplus", "muminus", "sigmaminus", "deltamu", "sigmadeltamu", "musigma", "muplussigmaplus", "muminussigmaminus", "deltamusigmadeltamu"};
            vector<string> YaxisNames { "#mu [#mum]", "#sigma_{#mu} [#mum]", "#mu outward [#mum]", "#sigma_{#mu outward} [#mum]", "#mu inward [#mum]", "#sigma_{#mu inward} [#mum]", "#Delta#mu [#mum]", "#sigma_{#Delta#mu} [#mum]", "#mu [#mum]", "#mu outward [#mum]", "#mu inward [#mum]", "#Delta#mu [#mum]",}; 
        //For debugging purposes we still might want to have a look at plots for a variable without errors, once ready for the PR those variables will be removed and the iterator will start from 0
            for(size_t i=8; i < variables.size(); i++){
                TString variable= variables.at(i);
                TCanvas * c = new TCanvas("dummy","",2000,800);
                vector<Color_t>::iterator colour = colours.begin();

                TMultiGraph *mg = new TMultiGraph(structure,structure);
                THStack *mh = new THStack(structure,structure);
        size_t igeom=0;
                for (string geometry: geometries) {
                    TString name = getName(structure, layer, geometry);
                    TGraphErrors *g = dynamic_cast<TGraphErrors*> (in->Get(name+"_"+variables.at(i)));
            g->SetName(name+"_"+variables.at(i));
                    if(i>=8){
                        g->SetLineWidth(1);
            g->SetLineColor(*colour);
                        g->SetFillColorAlpha(*colour,0.2);
                    }
            vector<vector<double>> vectors; 
            //if(showlumi&&i<8)scalebylumi(dynamic_cast<TGraph*>(g));
            if(showlumi)scalebylumi(g,Year);
            g->SetLineColor(*colour);
                    g->SetMarkerColor(*colour);
            TH1F *h = ConvertToHist(g); 
            h->SetLineColor(*colour);
            h->SetMarkerColor(*colour);
            h->SetMarkerSize(0);
            h->SetLineWidth(1);

                    if(i<8){
              mg->Add(g,"PL");
              mh->Add(h,"E");
                    }else{
              mg->Add(g,"2");
              mh->Add(h,"E");
                    }
                    ++colour;
            ++igeom;
                }
                if(i<8){
          mg->Draw("a");
                }else{
          mg->Draw("a2");
        }
                double max=6;
                double min=-4;
                double range=max-min;
                if(((variable=="sigma"||variable=="sigmaplus"||variable=="sigmaminus"||variable=="sigmadeltamu")&&range>=2)){
                    mg->SetMaximum(4);
                    mg->SetMinimum(-2);
                }else{
                    mg->SetMaximum(max+range*0.1);
                    mg->SetMinimum(min-range*0.3);
                }

                char* Ytitle= (char *)YaxisNames.at(i).c_str();
                mg->GetYaxis()->SetTitle(Ytitle);
                mg->GetXaxis()->SetTitle(showlumi ? "Integrated lumi [1/fb]" : "IOV number");
                mg->GetXaxis()->CenterTitle(true);
                mg->GetYaxis()->CenterTitle(true);
                mg->GetYaxis()->SetTitleOffset(.5);
                mg->GetYaxis()->SetTitleSize(.05);
                mg->GetXaxis()->SetTitleSize(.04);
        if(showlumi) mg->GetXaxis()->SetLimits(0.,mg->GetXaxis()->GetXmax());
                gStyle->SetOptTitle(0); // TODO
        gPad->SetTickx();
        gPad->SetTicky();
        //c->SetLeftMargin(0.11);
        
                c->Update();

        //gStyle->SetLegendBorderSize(0);
        gStyle->SetLegendTextSize(0.025);

                //TLegend *legend = c->BuildLegend();
                TLegend *legend = c->BuildLegend(0.3,0.15,0.3,0.15);
                    int Ngeom=geometries.size();
                legend->SetNColumns(Ngeom);
                //legend->SetTextSize(0.05);
                TString structtitle = structure;
                if(layer!=0){
                    if(structure=="TID"||structure=="TEC"||structure=="FPIX"||structure=="FPIX_y")structtitle+="_disc";
                    else structtitle+="_layer";
                    structtitle+=layer;
                }
        legend->SetHeader("#scale[1.2]{#bf{CMS} Work in progress}");
                //TLegendEntry *header = (TLegendEntry*)legend->GetListOfPrimitives()->First();
                //header->SetTextSize(.04);
        legend->AddEntry((TObject*)nullptr,structtitle.Data(),"h");
               
        //TLegendEntry *str = (TLegendEntry*)legend->GetListOfPrimitives()->Last();
                //str->SetTextSize(.03);
                PixelUpdateLines(c, Year, showlumi, pixelupdateruns);

        legend->Draw();
        mh->Draw("nostack same");
                c->Update();
                TString structandlayer = getName(structure,layer,"");
                TString printfile=outputdir;
        for(TString label : labels){printfile+=label;printfile+="_";}
        printfile+=variable+structandlayer;
        c->SaveAs(printfile+".pdf");
                c->SaveAs(printfile+".eps");
                c->SaveAs(printfile+".png");
        c->Destructor();
            }

        }

    }
    in->Close();
}

vector< int > runlistfromlumifile

(

TString

Year = "2018"

)

Definition at line 198 of file DMRtrends.cc.

References mps_fire::i, lumifileperyear(), N, writedatasetfile::runs, and Scenarios_cff::scale.

Referenced by checkrunlist(), and Point::GetSigmaDeltaMu().

                                             {
    TGraph * scale = new TGraph((lumifileperyear(Year,"run")).Data());
    double *xscale = scale->GetX();
    size_t N = scale->GetN();
    vector<int> runs;
    for(size_t i=0;i<N;i++)runs.push_back(xscale[i]);
    return runs;
} 

void scalebylumi	(	TGraphErrors *	g,
		TString	Year = "2018",
		double	min = 0.
	)

Definition at line 477 of file DMRtrends.cc.

References DUMMY, mps_fire::i, lumiFactor, lumifileperyear(), min(), N, writedatasetfile::run, Scenarios_cff::scale, and matplotRender::yscale.

Referenced by Point::GetSigmaDeltaMu(), and PlotDMRTrends().

                                                           { 
    size_t N=g->GetN();
    vector<double> x,y,xerr,yerr;

    TGraph * scale = new TGraph((lumifileperyear(Year,"IOV")).Data());
    size_t Nscale=scale->GetN();
    double *xscale=scale->GetX();
    double *yscale=scale->GetY();

    size_t i=0;
    while(i<N){
        double run,yvalue;
    g->GetPoint(i,run,yvalue);
        size_t index=-1;
        for(size_t j=0;j<Nscale;j++){
            if(run==xscale[j]){
                index=j;
                continue;
            }else if(run>xscale[j]) continue;
        }
        if(yscale[index]==0||index<0.){
            N=N-1;
            g->RemovePoint(i);
        }else{
        double xvalue=min;
            for(size_t j=0;j<index;j++)xvalue+=yscale[j]/lumiFactor;
        x   .push_back(xvalue+(yscale[index]/(lumiFactor*2.)));
        if(yvalue<=DUMMY){
          y.push_back(DUMMY);
          yerr.push_back(0.);
        }else{
          y.push_back(yvalue);
          yerr.push_back(g->GetErrorY(i));
        }
        xerr.push_back(yscale[index]/(lumiFactor*2.));
            i=i+1;
        }

    } 
    g->GetHistogram()->Delete(); 
    g->SetHistogram(nullptr);
    for(size_t i=0;i<N;i++){ g->SetPoint(i, x.at(i),y.at(i)); g->SetPointError(i, xerr.at(i),yerr.at(i));}

}