#include "TString.h"
#include "TFile.h"
#include "TProfile.h"
#include "TH1.h"
#include "THStack.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TLegend.h"
#include "TLatex.h"
#include <iostream>
#include <map>

Functions
bool	getHists (std::map< TString, TH1D * > &hists, TFile *&file, TString var, int var_val, float r)

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

void	setStack (THStack &stack, std::map< TString, TH1D > &hists)

void	setStyle ()

void	split (const TString &str, TString &sub1, TString &sub2)

Function Documentation

bool getHists	(	std::map< TString, TH1D * > &	hists,
		TFile *&	file,
		TString	var,
		int	var_val,
		float	r
	)

Definition at line 145 of file offsetStack.cc.

References gather_cfg::cout, mps_fire::i, dqmiodumpmetadata::n, mergeVDriftHistosByStation::name, AlCaHLTBitMon_ParallelJobs::p, and fw3dlego::xbins.

Referenced by main().

                                                                                               {
   for (auto& pair : hists) {
     TString name = Form("p_offset_eta_%s%i_%s", var.Data(), var_val, pair.first.Data());
     TProfile* p = (TProfile*)file->FindObjectAny(name);
     if (!p) {
       std::cout << "Could not find " + name << std::endl;
       return false;
     }
     pair.second = p->ProjectionX(pair.first);
     pair.second->Scale(r * r / 2 / var_val);
 
     const double* xbins = p->GetXaxis()->GetXbins()->GetArray();
     for (int i = 1, n = p->GetNbinsX(); i <= n; i++) {
       pair.second->SetBinContent(i, pair.second->GetBinContent(i) / (xbins[i] - xbins[i - 1]));
       pair.second->SetBinError(i, pair.second->GetBinError(i) / (xbins[i] - xbins[i - 1]));
     }
   }
   return true;
 }

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

===============================================================================================================================================================================================

variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

CALO JETS

PF JETS

Definition at line 19 of file offsetStack.cc.

References c, gather_cfg::cout, timingPdfMaker::file1, getHists(), h, label, alignmentValidation::outdir, alignCSCRings::r, setStack(), setStyle(), submitPVValidationJobs::split(), runonSM::text, and isotrackApplyRegressor::var.

                                  {
   if (argc < 5) {
     std::cout << "Please provide a label:filename, variable, deltaR value, and an out directory." << std::endl;
     std::cout << "You may also include a second file at the end." << std::endl;
     std::cout << "Example: offsetStack label1:file1.root npv 0.4 plots" << std::endl;
     std::cout << "Example: offsetStack label1:file1.root npv 0.4 plots label2:file2.root" << std::endl;
     return -1;
   }
 
   TString label1, fname1;
   split(argv[1], label1, fname1);
   TString label = label1;
 
   TString var = argv[2];
   float r = std::stof(argv[3]);
   TString outdir = argv[4];
 
   TFile* file1 = TFile::Open(fname1);
   if (!file1) {
     std::cout << "Invalid file1: " << fname1 << std::endl;
     return -1;
   }
 
   TH1F* h = (TH1F*)file1->FindObjectAny(var);
   if (!h) {
     std::cout << "Could not find " + var + " hist" << std::endl;
     return -1;
   }
   int avg = h->GetMean() + 0.5;
 
   std::map<TString, TH1D*> hists1 = {{"chm", 0}, {"chu", 0}, {"nh", 0}, {"ne", 0}, {"hfh", 0}, {"hfe", 0}, {"lep", 0}};
   if (!getHists(hists1, file1, var, avg, r))
     return -1;
 
   setStyle();
   TCanvas c("c", "c", 600, 600);
 
   THStack* stack1 = new THStack(
       "stack1", Form(";#eta;<Offset Energy_{T}(#pi%.1f^{2})>/<%s> [GeV]", r, var == "npv" ? "N_{PV}" : "#mu"));
   setStack(stack1, hists1);
   stack1->Draw("hist");
   TString legPF = "F";
 
   TLegend* leg = new TLegend(.4, .67, .65, .92);
 
   //file2 included//
   if (argc > 5) {
     TString label2, fname2;
     split(argv[5], label2, fname2);
 
     TFile* file2 = TFile::Open(fname2);
     if (!file2) {
       std::cout << "Invalid file2: " << fname2 << std::endl;
       return -1;
     }
 
     TH1F* h2 = (TH1F*)file2->FindObjectAny(var);
     if (!h2) {
       std::cout << "Could not find " + var + " hist in file2" << std::endl;
       return -1;
     }
     int avg2 = h2->GetMean() + 0.5;
 
     std::map<TString, TH1D*> hists2 = hists1;
     if (!getHists(hists2, file2, var, avg2, r))
       return -1;
 
     THStack* stack2 = new THStack("stack2", "stack2");
     setStack(stack2, hists2);
     stack2->Draw("samepe");
 
     legPF = "PF";
     leg->SetHeader(Form("#bf{Markers: %s, Histograms: %s}", label2.Data(), label1.Data()));
     label = label1 + "_vs_" + label2;
 
     //Draw Markers for EM Deposits and Hadronic Deposits in two separate regions//
     TH1D* hfe_clone = (TH1D*)hists2["hfe"]->Clone("hfe_clone");
     TH1D* hfh_clone = (TH1D*)hists2["hfh"]->Clone("hfh_clone");
 
     THStack* cloneStack = new THStack("cloneStack", "cloneStack");
     cloneStack->Add(hists2["ne"]);
     cloneStack->Add(hfe_clone);
     cloneStack->Add(hists2["nh"]);
     cloneStack->Add(hfh_clone);
     cloneStack->Draw("samepe");
 
     hists2["ne"]->SetAxisRange(-2.9, 2.9);
     hists2["hfe"]->SetAxisRange(-5, -2.6);
     hists2["nh"]->SetAxisRange(-2.9, 2.9);
     hists2["hfh"]->SetAxisRange(-5, -2.6);
     hists2["chu"]->SetAxisRange(-2.9, 2.9);
     hists2["chm"]->SetAxisRange(-2.9, 2.9);
 
     hfe_clone->SetAxisRange(2.6, 5);
     hfh_clone->SetAxisRange(2.6, 5);
   }
 
   leg->SetBorderSize(0);
   leg->SetFillColor(0);
   leg->SetFillStyle(0);
   leg->SetTextSize(0.04);
   leg->SetTextFont(42);
 
   leg->AddEntry(hists1["ne"], "Photons", legPF);
   leg->AddEntry(hists1["hfe"], "EM Deposits", legPF);
   leg->AddEntry(hists1["nh"], "Neutral Hadrons", legPF);
   leg->AddEntry(hists1["hfh"], "Hadronic Deposits", legPF);
   leg->AddEntry(hists1["chu"], "Unassoc. Charged Hadrons", legPF);
   leg->AddEntry(hists1["chm"], "Assoc. Charged Hadrons", legPF);
 
   leg->Draw();
 
   TLatex text;
   text.SetNDC();
 
   text.SetTextSize(0.065);
   text.SetTextFont(61);
   text.DrawLatex(0.2, 0.87, "CMS");
 
   text.SetTextSize(0.045);
   text.SetTextFont(42);
   text.DrawLatex(1 - label.Length() / 41., 0.96, label);
 
   c.Print(outdir + "/stack_" + label + ".pdf");
 }

void setStack	(	THStack *&	stack,
		std::map< TString, TH1D * > &	hists
	)

Definition at line 187 of file offsetStack.cc.

Referenced by main().

                                                               {
   stack->Add(hists["ne"]);
   stack->Add(hists["hfe"]);
   stack->Add(hists["nh"]);
   stack->Add(hists["hfh"]);
   stack->Add(hists["chu"]);
   stack->Add(hists["chm"]);
 
   hists["ne"]->SetMarkerStyle(kMultiply);
   hists["hfe"]->SetMarkerStyle(kOpenStar);
   hists["nh"]->SetMarkerStyle(kOpenDiamond);
   hists["hfh"]->SetMarkerStyle(kOpenTriangleUp);
   hists["chu"]->SetMarkerStyle(kOpenCircle);
   hists["chm"]->SetMarkerStyle(kOpenCircle);
 
   hists["ne"]->SetFillColor(kBlue);
   hists["hfe"]->SetFillColor(kViolet + 2);
   hists["nh"]->SetFillColor(kGreen);
   hists["hfh"]->SetFillColor(kPink + 6);
   hists["chu"]->SetFillColor(kRed - 9);
   hists["chm"]->SetFillColor(kRed);
 
   hists["ne"]->SetLineColor(kBlack);
   hists["hfe"]->SetLineColor(kBlack);
   hists["nh"]->SetLineColor(kBlack);
   hists["hfh"]->SetLineColor(kBlack);
   hists["chu"]->SetLineColor(kBlack);
   hists["chm"]->SetLineColor(kBlack);
 }

void setStyle ( )

set the style for the printout

Definition at line 165 of file offsetStack.cc.

Referenced by getGlobalCanvas(), and main().

                 {
   gStyle->SetPadTopMargin(0.05);
   gStyle->SetPadBottomMargin(0.1);
   gStyle->SetPadLeftMargin(0.16);
   gStyle->SetPadRightMargin(0.02);
 
   gStyle->SetOptStat(0);
   gStyle->SetOptTitle(0);
 
   gStyle->SetTitleFont(42, "XYZ");
   gStyle->SetTitleSize(0.05, "XYZ");
   gStyle->SetTitleXOffset(0.9);
   gStyle->SetTitleYOffset(1.4);
 
   gStyle->SetLabelFont(42, "XYZ");
   gStyle->SetLabelOffset(0.007, "XYZ");
   gStyle->SetLabelSize(0.04, "XYZ");
 
   gStyle->SetPadTickX(1);
   gStyle->SetPadTickY(1);
 }

void split	(	const TString &	str,
		TString &	sub1,
		TString &	sub2
	)

Definition at line 217 of file offsetStack.cc.

References str.

                                                              {
   int pos = str.First(':');
   sub1 = str(0, pos);
   sub2 = str(pos + 1, str.Length());
 }

Functions

Function Documentation

CALO JETS

PF JETS