Back to home page

Project CMSSW displayed by LXR
 
 

    

File indexing completed on 2023-03-17 11:27:44 

0001 #include <fstream>
0002 #include <TFile.h>
0003 #include <TString.h>
0004 #include <TCanvas.h>
0005 
0006 
0007 
0008 using namespace std;
0009 
0010 void prn(TString s0, TString s1) {
0011   std::cout << "\t>> " << s0 << ": " << s1 << std::endl;
0012 }
0013 
0014 void prn(TString s0, double d) {
0015   std::cout << "\t>> " << s0 << ": " << d << std::endl;
0016 }
0017 
0018 
0019 void ProcessRelVal(TFile &ref_file, TFile &val_file, ifstream &recstr, const int nHist1, const int nHist2, const int nHits2D, const int nProfInd, const int nHistTot, TString ref_vers, TString val_vers, int harvest = 0, bool bRBX = false);
0020 
0021 void RelValMacro(TString ref_vers = "218", TString val_vers = "218", TString rfname, TString vfname, TString InputStream = "InputRelVal.txt", int harvest = 0) {
0022 
0023   ifstream RelValStream;
0024 
0025   RelValStream.open(InputStream);
0026 
0027   TFile Ref_File(rfname);
0028   TFile Val_File(vfname);
0029 
0030   /*
0031     A note about MC histograms: 22 of them are not included in this current implementation.
0032     Two things must be done to include them again:
0033     1. Uncomment the appropriate lines below in the service variables detailing how many
0034        histograms ProcessRelVal will expect from the Input txt files.
0035     2. Change the aforementioned Input txt files DrawSwitch flags from 0 to 1 (the second 
0036        number after the histogram name). 
0037   */
0038 
0039 
0040   //Service variables
0041 
0042   //SimHits
0043   const int SH_nHistTot = 154;
0044   const int SH_nHist1 = 54;
0045   const int SH_nHist2 = 0;
0046   const int SH_nHist2D = 100;
0047   const int SH_nProfInd = 0;
0048   
0049   ProcessRelVal(Ref_File, Val_File, RelValStream, SH_nHist1, SH_nHist2, SH_nHist2D, SH_nProfInd, SH_nHistTot, ref_vers, val_vers, harvest, false, true);
0050   
0051   Ref_File.Close();
0052   Val_File.Close();
0053   
0054   return;
0055 }
0056 
0057 void ProcessRelVal(TFile &ref_file, TFile &val_file, ifstream &recstr, const int nHist1, const int nHist2, const int nHist2D, const int nProfInd, const int nHistTot, TString ref_vers, TString val_vers, int harvest, bool bRBX, bool bHD = false) {
0058   
0059   TString RefHistDir, ValHistDir;
0060   
0061   if (bRBX) {
0062     if (harvest == 11) {
0063       RefHistDir = "DQMData/Run 1/NoiseRatesV/Run summary/NoiseRatesTask";
0064       ValHistDir = "DQMData/Run 1/NoiseRatesV/Run summary/NoiseRatesTask";
0065     } else if (harvest == 10) {
0066       RefHistDir = "DQMData/NoiseRatesV/NoiseRatesTask";
0067       ValHistDir = "DQMData/Run 1/NoiseRatesV/Run summary/NoiseRatesTask";
0068     } else if (harvest == 1) {
0069       RefHistDir = "DQMData/Run 1/NoiseRatesV/Run summary/NoiseRatesTask";
0070       ValHistDir = "DQMData/NoiseRatesV/NoiseRatesTask";
0071     } else {
0072       RefHistDir = "DQMData/NoiseRatesV/NoiseRatesTask";
0073       ValHistDir = "DQMData/NoiseRatesV/NoiseRatesTask";
0074     }
0075   } else if (bHD) {
0076     if (harvest == 11) {
0077       RefHistDir = "DQMData/Run 2/HcalHitsV/Run summary/SimHitsValidationHcal";
0078       ValHistDir = "DQMData/Run 2/HcalHitsV/Run summary/SimHitsValidationHcal";
0079     } else if (harvest == 10) {
0080       RefHistDir = "DQMData/HcalDigisV/HcalDigiTask";
0081       ValHistDir = "DQMData/Run 1/HcalDigisV/Run summary/HcalDigiTask";
0082     } else if (harvest == 1) {
0083       RefHistDir = "DQMData/Run 1/HcalDigisV/Run summary/HcalDigiTask";
0084       ValHistDir = "DQMData/HcalDigisV/HcalDigiTask";
0085     } else {
0086       RefHistDir = "DQMData/HcalDigisV/HcalDigiTask";
0087       ValHistDir = "DQMData/HcalDigisV/HcalDigiTask";
0088     }
0089   } else {
0090     if (harvest == 11) {
0091       RefHistDir = "DQMData/Run 1/HcalRecHitsV/Run summary/HcalRecHitTask";
0092       ValHistDir = "DQMData/Run 1/HcalRecHitsV/Run summary/HcalRecHitTask";
0093     } else if (harvest == 10) {
0094       RefHistDir = "DQMData/HcalRecHitsV/HcalRecHitTask";
0095       ValHistDir = "DQMData/Run 1/HcalRecHitsV/Run summary/HcalRecHitTask";
0096     } else if (harvest == 1) {
0097       RefHistDir = "DQMData/Run 1/HcalRecHitsV/Run summary/HcalRecHitTask";
0098       ValHistDir = "DQMData/HcalRecHitsV/HcalRecHitTask";
0099     } else {
0100       RefHistDir = "DQMData/HcalRecHitsV/HcalRecHitTask";
0101       ValHistDir = "DQMData/HcalRecHitsV/HcalRecHitTask";
0102     }
0103   }
0104   
0105   TCanvas* myc = 0;
0106   TLegend* leg = 0;
0107   TPaveText* ptchi2 = 0;
0108   TPaveStats *ptstats_r = 0;
0109   TPaveStats *ptstats_v = 0;
0110   
0111   TH1F * ref_hist1[nHist1];
0112   TH2F * ref_hist2[nHist2];
0113   TProfile * ref_prof[nProfInd];
0114   TH1D * ref_fp[nProfInd];
0115   TH2F * ref_hist2D[nHist2D];
0116   
0117   TH1F * val_hist1[nHist1];
0118   TH2F * val_hist2[nHist2];
0119   TProfile * val_prof[nProfInd];
0120   TH1D * val_fp[nProfInd];
0121   TH2F * val_hist2D[nHist2D];
0122   
0123   int i;
0124   int DrawSwitch;
0125   TString StatSwitch, Chi2Switch, LogSwitch, DimSwitch;
0126   int RefCol, ValCol;
0127   TString HistName, HistName2;
0128   char xAxisTitle[200];
0129   int nRebin;
0130   float xAxisMin, xAxisMax, yAxisMin, yAxisMax;
0131   TString OutLabel, ProfileLabel;
0132   string xTitleCheck;
0133   
0134   float hmax = 0;
0135 
0136   int nh1 = 0;
0137   int nh2 = 0;
0138   int npr = 0;
0139   int npi = 0;
0140   int n2D = 0;
0141   
0142   for (i = 0; i < nHistTot; i++) {
0143     
0144     //Read in 1/0 switch saying whether this histogram is used
0145     //Skip it if not used, otherwise get output file label, histogram
0146     //axis ranges and title
0147     //ALTERED: Reads in all inputs and then uses 1/0 switch to skip
0148     //See below
0149     recstr >> HistName >> DrawSwitch;
0150     prn("HistName", HistName);
0151     //        if (DrawSwitch == 0) continue;
0152       
0153     recstr >> OutLabel >> nRebin;
0154     recstr >> xAxisMin >> xAxisMax >> yAxisMin >> yAxisMax;
0155     recstr >> DimSwitch >> StatSwitch >> Chi2Switch >> LogSwitch;
0156     recstr >> RefCol >> ValCol;
0157     recstr.getline(xAxisTitle, 200);
0158     
0159     //Make sure extra Profile info is also taken care of
0160     if (DrawSwitch == 0) {
0161       if (DimSwitch == "TM") 
0162     recstr >> ProfileLabel;
0163       continue;
0164     }
0165     
0166     
0167     // Nasty trick:
0168     // recovering colors redefined in rootlogon.C (for "rainbow" Palette)
0169     Float_t r, g, b;
0170     Float_t saturation = 1;
0171     Float_t lightness = 0.5;
0172     Float_t maxHue = 280;
0173     Float_t minHue = 0;
0174     Int_t maxPretty = 50;
0175     Float_t hue;
0176 
0177     for (int j = 0; j < maxPretty; j++) {
0178       hue = maxHue - (j + 1)*((maxHue - minHue) / maxPretty);
0179       TColor::HLStoRGB(hue, lightness, saturation, r, g, b);
0180       TColor *color = (TColor*) (gROOT->GetListOfColors()->At(j + 51));
0181       color->SetRGB(r, g, b);
0182     }
0183     gStyle->SetPalette(1);
0184 
0185 
0186     //Format canvas
0187     if (DimSwitch == "PRwide") {
0188       gStyle->SetPadLeftMargin(0.06);
0189       gStyle->SetPadRightMargin(0.03);
0190       myc = new TCanvas("myc", "", 1200, 600);
0191     } else myc = new TCanvas("myc", "", 800, 600);
0192     myc->SetGrid();
0193 
0194     xTitleCheck = xAxisTitle;
0195     xTitleCheck = xTitleCheck.substr(1, 7);
0196 
0197     //Format pad
0198     if (LogSwitch == "Log") myc->SetLogy(1);
0199     else myc->SetLogy(0);
0200 
0201     //AF
0202     if (LogSwitch == "Log" && DimSwitch == "2D"){
0203       myc->SetLogy(0);
0204       myc->SetLogz(1);
0205     }
0206 
0207     if (DimSwitch == "1D") {
0208       //Get histograms from files
0209       ref_file.cd(RefHistDir);
0210       ref_hist1[nh1] = (TH1F*) gDirectory->Get(HistName);
0211       
0212       val_file.cd(ValHistDir);
0213       val_hist1[nh1] = (TH1F*) gDirectory->Get(HistName);
0214 
0215       // HACK to change what is embedded in DQM histos
0216       ref_hist1[nh1]->GetXaxis()->SetLabelSize(0.04);
0217       val_hist1[nh1]->GetXaxis()->SetLabelSize(0.04);
0218       ref_hist1[nh1]->GetYaxis()->SetLabelSize(0.04);
0219       val_hist1[nh1]->GetYaxis()->SetLabelSize(0.04);
0220       ref_hist1[nh1]->GetXaxis()->SetTitleSize(0.045);
0221       val_hist1[nh1]->GetXaxis()->SetTitleSize(0.045);
0222 
0223       ref_hist1[nh1]->GetXaxis()->SetTickLength(-0.015);
0224       val_hist1[nh1]->GetXaxis()->SetTickLength(-0.015);
0225       ref_hist1[nh1]->GetYaxis()->SetTickLength(-0.015);
0226       val_hist1[nh1]->GetYaxis()->SetTickLength(-0.015);
0227 
0228       ref_hist1[nh1]->GetXaxis()->SetLabelOffset(0.012);
0229       val_hist1[nh1]->GetXaxis()->SetLabelOffset(0.012);
0230       ref_hist1[nh1]->GetYaxis()->SetLabelOffset(0.012);
0231       val_hist1[nh1]->GetYaxis()->SetLabelOffset(0.012);
0232 
0233       ref_hist1[nh1]->GetXaxis()->SetTitleOffset(0.8);
0234       val_hist1[nh1]->GetXaxis()->SetTitleOffset(0.8);
0235 
0236 
0237       //Rebin histograms -- has to be done first
0238       if (nRebin != 1) {
0239     ref_hist1[nh1]->Rebin(nRebin);
0240     val_hist1[nh1]->Rebin(nRebin);
0241       }
0242 
0243       //Set the colors, styles, titles, stat boxes and format axes for the histograms
0244       ref_hist1[nh1]->SetStats(kTRUE);
0245       val_hist1[nh1]->SetStats(kTRUE);
0246 
0247       if (StatSwitch != "Stat" && StatSwitch != "Statrv") {
0248     ref_hist1[nh1]->SetStats(kFALSE);
0249     val_hist1[nh1]->SetStats(kFALSE);
0250       }
0251 
0252       //Min/Max Convetion: Default AxisMin = 0. Default AxisMax = -1.
0253       //xAxis
0254       if (xAxisMin == 0) xAxisMin = ref_hist1[nh1]->GetXaxis()->GetXmin();
0255       if (xAxisMax < 0) xAxisMax = ref_hist1[nh1]->GetXaxis()->GetXmax();
0256 
0257       if (xAxisMax > 0 || xAxisMin != 0) {
0258     ref_hist1[nh1]->GetXaxis()->SetRangeUser(xAxisMin, xAxisMax);
0259     val_hist1[nh1]->GetXaxis()->SetRangeUser(xAxisMin, xAxisMax);
0260       }
0261       //yAxis
0262       if (yAxisMin != 0) ref_hist1[nh1]->SetMinimum(yAxisMin);
0263       if (yAxisMax > 0) ref_hist1[nh1]->SetMaximum(yAxisMax);
0264       else if (ref_hist1[nh1]->GetMaximum() < val_hist1[nh1]->GetMaximum() &&
0265            val_hist1[nh1]->GetMaximum() > 0) {
0266     if (LogSwitch == "Log") ref_hist1[nh1]->SetMaximum(2 * val_hist1[nh1]->GetMaximum());
0267     else ref_hist1[nh1]->SetMaximum(1.05 * val_hist1[nh1]->GetMaximum());
0268       }
0269 
0270       //Title
0271       if (xTitleCheck != "NoTitle") ref_hist1[nh1]->GetXaxis()->SetTitle(xAxisTitle);
0272 
0273       //Different histo colors and styles
0274       ref_hist1[nh1]->SetTitle("");
0275       ref_hist1[nh1]->SetLineColor(RefCol);
0276       ref_hist1[nh1]->SetLineStyle(1);
0277       ref_hist1[nh1]->SetMarkerSize(0.02);
0278       if (StatSwitch != "Stat" && StatSwitch != "Statrv") ref_hist1[nh1]->SetLineWidth(2);
0279 
0280       val_hist1[nh1]->SetTitle("");
0281       val_hist1[nh1]->SetLineColor(ValCol);
0282       val_hist1[nh1]->SetLineStyle(2);
0283       val_hist1[nh1]->SetMarkerSize(0.02);
0284       if (StatSwitch != "Stat" && StatSwitch != "Statrv") val_hist1[nh1]->SetLineWidth(2);
0285 
0286       //Legend
0287       TLegend *leg = new TLegend(0.50, 0.91, 0.84, 0.99, "", "brNDC");
0288       leg->SetBorderSize(2);
0289       leg->SetFillStyle(1001);
0290       leg->AddEntry(ref_hist1[nh1], "CMSSW_" + ref_vers, "l");
0291       leg->AddEntry(val_hist1[nh1], "CMSSW_" + val_vers, "l");
0292 
0293       if (Chi2Switch == "Chi2") {
0294     //Draw and save histograms
0295     ref_hist1[nh1]->SetFillColor(40);//42 Originally, now 40 which is lgiht brown
0296     ref_hist1[nh1]->Draw("hist");
0297     val_hist1[nh1]->SetLineStyle(1);
0298     if (StatSwitch == "Statrv") val_hist1[nh1]->Draw("sames e0");
0299     else val_hist1[nh1]->Draw("same e0");
0300 
0301     //Get p-value from chi2 test
0302     const float NCHI2MIN = 0.01;
0303     
0304     float pval;
0305     stringstream mystream;
0306     char tempbuff[30];
0307 
0308     pval = ref_hist1[nh1]->Chi2Test(val_hist1[nh1]);
0309     
0310     sprintf(tempbuff, "Chi2 p-value: %6.3E%c", pval, '\0');
0311     mystream << tempbuff;
0312 
0313     ptchi2 = new TPaveText(0.05, 0.92, 0.35, 0.99, "NDC");
0314 
0315     if (pval > NCHI2MIN) ptchi2->SetFillColor(kGreen);
0316     else ptchi2->SetFillColor(kRed);
0317 
0318     ptchi2->SetTextSize(0.03);
0319     ptchi2->AddText(mystream.str().c_str());
0320     ptchi2->Draw();
0321       } else {
0322     //Draw and save histograms
0323     ref_hist1[nh1]->Draw("hist");
0324     if (StatSwitch == "Statrv") val_hist1[nh1]->Draw("hist sames");
0325     else val_hist1[nh1]->Draw("hist same");
0326       }
0327 
0328       //Stat Box where required
0329       if (StatSwitch == "Stat" || StatSwitch == "Statrv") {
0330     ptstats_r = new TPaveStats(0.85, 0.86, 0.98, 0.98, "brNDC");
0331     ptstats_r->SetTextColor(RefCol);
0332     ref_hist1[nh1]->GetListOfFunctions()->Add(ptstats_r);
0333     ptstats_r->SetParent(ref_hist1[nh1]->GetListOfFunctions());
0334     ptstats_v = new TPaveStats(0.85, 0.74, 0.98, 0.86, "brNDC");
0335     ptstats_v->SetTextColor(ValCol);
0336     val_hist1[nh1]->GetListOfFunctions()->Add(ptstats_v);
0337     ptstats_v->SetParent(val_hist1[nh1]->GetListOfFunctions());
0338 
0339     ptstats_r->Draw();
0340     ptstats_v->Draw();
0341       }
0342 
0343       leg->Draw();
0344 
0345       myc->SaveAs(OutLabel);
0346       nh1++;
0347     } else if (DimSwitch == "PR" || DimSwitch == "PRwide") {
0348       //Profiles not associated with histograms
0349       //Get profiles from files
0350       ref_file.cd(RefHistDir);
0351       ref_prof[npi] = (TProfile*) gDirectory->Get(HistName);
0352       
0353       val_file.cd(ValHistDir);
0354       val_prof[npi] = (TProfile*) gDirectory->Get(HistName);
0355 
0356       // HACK to change what is embedded in DQM histos
0357       ref_prof[npi]->GetXaxis()->SetLabelSize(0.04);
0358       val_prof[npi]->GetXaxis()->SetLabelSize(0.04);
0359       ref_prof[npi]->GetYaxis()->SetLabelSize(0.04);
0360       val_prof[npi]->GetYaxis()->SetLabelSize(0.04);
0361       ref_prof[npi]->GetXaxis()->SetTitleSize(0.045);
0362       val_prof[npi]->GetXaxis()->SetTitleSize(0.045);
0363 
0364       ref_prof[npi]->GetXaxis()->SetTickLength(-0.015);
0365       val_prof[npi]->GetXaxis()->SetTickLength(-0.015);
0366       ref_prof[npi]->GetYaxis()->SetTickLength(-0.015);
0367       val_prof[npi]->GetYaxis()->SetTickLength(-0.015);
0368         
0369       ref_prof[npi]->GetXaxis()->SetLabelOffset(0.02);
0370       val_prof[npi]->GetXaxis()->SetLabelOffset(0.02);
0371       ref_prof[npi]->GetYaxis()->SetLabelOffset(0.02);
0372       val_prof[npi]->GetYaxis()->SetLabelOffset(0.02);
0373       
0374       ref_prof[npi]->GetXaxis()->SetTitleOffset(1.3);
0375       val_prof[npi]->GetXaxis()->SetTitleOffset(1.3);
0376 
0377 
0378       //Legend
0379       leg = new TLegend(0.50, 0.91, 0.84, 0.99, "", "brNDC");
0380       leg->SetBorderSize(2);
0381       leg->SetFillStyle(1001);
0382 
0383       //Ordinary profiles
0384       if (DimSwitch == "PR") {
0385     ref_prof[npi]->SetTitle("");
0386     ref_prof[npi]->SetErrorOption("");
0387 
0388     val_prof[npi]->SetTitle("");
0389     val_prof[npi]->SetErrorOption("");
0390 
0391     ref_prof[npi]->GetXaxis()->SetTitle(xAxisTitle);
0392 
0393     if (StatSwitch != "Stat" && StatSwitch != "Statrv") {
0394       ref_prof[npi]->SetStats(kFALSE);
0395       val_prof[npi]->SetStats(kFALSE);
0396     }
0397 
0398     ref_prof[npi]->SetLineColor(41);
0399     ref_prof[npi]->SetLineStyle(1);
0400     ref_prof[npi]->SetLineWidth(1);
0401     ref_prof[npi]->SetMarkerColor(41);
0402     ref_prof[npi]->SetMarkerStyle(21);
0403     ref_prof[npi]->SetMarkerSize(0.8);
0404 
0405     val_prof[npi]->SetLineColor(43);
0406     val_prof[npi]->SetLineStyle(1);
0407     val_prof[npi]->SetLineWidth(1);
0408     val_prof[npi]->SetMarkerColor(43);
0409     val_prof[npi]->SetMarkerStyle(22);
0410     val_prof[npi]->SetMarkerSize(1.0);
0411 
0412     if (ref_prof[npi]->GetMaximum() < val_prof[npi]->GetMaximum() &&
0413         val_prof[npi]->GetMaximum() > 0) {
0414       if (LogSwitch == "Log") ref_prof[npi]->SetMaximum(2 * val_prof[npi]->GetMaximum());
0415       else ref_prof[npi]->SetMaximum(1.05 * val_prof[npi]->GetMaximum());
0416     }
0417 
0418     ref_prof[npi]->Draw("hist pl");
0419     val_prof[npi]->Draw("hist pl same");
0420 
0421     leg->AddEntry(ref_prof[npi], "CMSSW_" + ref_vers, "pl");
0422     leg->AddEntry(val_prof[npi], "CMSSW_" + val_vers, "pl");
0423       } else if (DimSwitch == "PRwide") {
0424     //Wide profiles
0425     TString temp = HistName + "_px_v";
0426     ref_fp[npi] = ref_prof[npi]->ProjectionX();
0427     val_fp[npi] = val_prof[npi]->ProjectionX(temp.Data());
0428     
0429     ref_fp[npi]->SetTitle("");
0430     val_fp[npi]->SetTitle("");
0431 
0432     ref_fp[npi]->GetXaxis()->SetTitle(xAxisTitle);
0433     
0434     if (StatSwitch != "Stat" && StatSwitch != "Statrv") {
0435       ref_fp[npi]->SetStats(kFALSE);
0436       val_fp[npi]->SetStats(kFALSE);
0437     }
0438 
0439     int nbins = ref_fp[npi]->GetNbinsX();
0440     for (int j = 1; j < nbins; j++) {
0441       ref_fp[npi]->SetBinError(j, 0.);
0442       val_fp[npi]->SetBinError(j, 0.);
0443     }
0444     ref_fp[npi]->SetLineWidth(0);
0445     ref_fp[npi]->SetLineColor(0); // 5 yellow
0446     ref_fp[npi]->SetLineStyle(1);
0447     ref_fp[npi]->SetMarkerColor(2);
0448     ref_fp[npi]->SetMarkerStyle(20);
0449     ref_fp[npi]->SetMarkerSize(0.5);
0450 
0451     val_fp[npi]->SetLineWidth(0);
0452     val_fp[npi]->SetLineColor(0); // 45 blue
0453     val_fp[npi]->SetLineStyle(2);
0454     val_fp[npi]->SetMarkerColor(4);
0455     val_fp[npi]->SetMarkerStyle(22);
0456     val_fp[npi]->SetMarkerSize(0.5);
0457 
0458     if (ref_fp[npi]->GetMaximum() < val_fp[npi]->GetMaximum() &&
0459         val_fp[npi]->GetMaximum() > 0) {
0460       if (LogSwitch == "Log") ref_fp[npi]->SetMaximum(2 * val_fp[npi]->GetMaximum());
0461       else ref_fp[npi]->SetMaximum(1.05 * val_fp[npi]->GetMaximum());
0462     }
0463 
0464     ref_fp[npi]->Draw("p9");
0465     val_fp[npi]->Draw("p9same");
0466 
0467     leg->AddEntry(ref_fp[npi], "CMSSW_" + ref_vers, "lp");
0468     leg->AddEntry(val_fp[npi], "CMSSW_" + val_vers, "lp");
0469 
0470       }
0471 
0472       leg->Draw("");
0473 
0474       myc->SaveAs(OutLabel);
0475       
0476       npi++;
0477     } else if (DimSwitch == "TM") {
0478       //Timing Histograms (special: read two lines at once)
0479       recstr >> HistName2;
0480 
0481       ref_file.cd(RefHistDir);
0482 
0483       ref_hist2[nh2] = (TH2F*) gDirectory->Get(HistName);
0484       ref_prof[npi] = (TProfile*) gDirectory->Get(HistName2);
0485 
0486       ref_hist2[nh2]->SetMarkerStyle(21);
0487       ref_prof[npi] ->SetMarkerStyle(21);
0488       ref_hist2[nh2]->SetMarkerSize(0.02);
0489       ref_prof[npi] ->SetMarkerSize(0.02);
0490 
0491       val_file.cd(ValHistDir);
0492       
0493       val_hist2[nh2] = (TH2F*) gDirectory->Get(HistName);
0494       val_prof[npi] = (TProfile*) gDirectory->Get(HistName2);
0495 
0496       val_hist2[nh2]->SetMarkerStyle(21);
0497       val_prof[npi] ->SetMarkerStyle(21);
0498       val_hist2[nh2]->SetMarkerSize(0.02);
0499       val_prof[npi] ->SetMarkerSize(0.02);
0500 
0501       // HACK to change what is embedded in DQM histos
0502       ref_hist2[nh2]->GetXaxis()->SetLabelSize(0.04);
0503       val_hist2[nh2]->GetXaxis()->SetLabelSize(0.04);
0504       ref_hist2[nh2]->GetYaxis()->SetLabelSize(0.04);
0505       val_hist2[nh2]->GetYaxis()->SetLabelSize(0.04);
0506       ref_hist2[nh2]->GetXaxis()->SetTitleSize(0.045);
0507       val_hist2[nh2]->GetXaxis()->SetTitleSize(0.045);
0508       
0509       ref_hist2[nh2]->GetXaxis()->SetTickLength(-0.015);
0510       val_hist2[nh2]->GetXaxis()->SetTickLength(-0.015);
0511       ref_hist2[nh2]->GetYaxis()->SetTickLength(-0.015);
0512       val_hist2[nh2]->GetYaxis()->SetTickLength(-0.015);
0513 
0514       ref_hist2[nh2]->GetXaxis()->SetLabelOffset(0.02);
0515       val_hist2[nh2]->GetXaxis()->SetLabelOffset(0.02);
0516       ref_hist2[nh2]->GetYaxis()->SetLabelOffset(0.02);
0517       val_hist2[nh2]->GetYaxis()->SetLabelOffset(0.02);
0518 
0519       ref_hist2[nh2]->GetXaxis()->SetTitleOffset(1.3);
0520       val_hist2[nh2]->GetXaxis()->SetTitleOffset(1.3);
0521 
0522 
0523       //Min/Max Convetion: Default AxisMin = 0. Default AxisMax = -1.
0524       //xAxis
0525       if (xAxisMin == 0) xAxisMin = ref_hist2[nh2]->GetXaxis()->GetXmin();
0526       if (xAxisMax < 0) xAxisMax = ref_hist2[nh2]->GetXaxis()->GetXmax();
0527 
0528       if (xAxisMax > 0 || xAxisMin != 0) {
0529     ref_hist2[nh2]->GetXaxis()->SetRangeUser(xAxisMin, xAxisMax);
0530     val_hist2[nh2]->GetXaxis()->SetRangeUser(xAxisMin, xAxisMax);
0531       }
0532       //yAxis
0533       if (yAxisMin != 0) ref_hist2[nh2]->SetMinimum(yAxisMin);
0534       if (yAxisMax > 0) ref_hist2[nh2]->SetMaximum(yAxisMax);
0535       else if (ref_hist2[nh2]->GetMaximum() < val_hist2[nh2]->GetMaximum() &&
0536            val_hist2[nh2]->GetMaximum() > 0) {
0537     if (LogSwitch == "Log") ref_hist2[nh2]->SetMaximum(2 * val_hist2[nh2]->GetMaximum());
0538     else ref_hist2[nh2]->SetMaximum(1.05 * val_hist2[nh2]->GetMaximum());
0539       }
0540 
0541       //AF
0542       if (yAxisMax > 0 || yAxisMin != 0) {
0543     ref_hist2[nh2]->GetYaxis()->SetRangeUser(yAxisMin, yAxisMax);
0544     val_hist2[nh2]->GetYaxis()->SetRangeUser(yAxisMin, yAxisMax);
0545       }
0546 
0547       //Legend
0548       leg = new TLegend(0.50, 0.91, 0.84, 0.99, "", "brNDC");
0549       leg->SetBorderSize(2);
0550       leg->SetFillStyle(1001);
0551 
0552       ref_hist2[nh2]->GetXaxis()->SetTitle(xAxisTitle);
0553       ref_hist2[nh2]->SetStats(kFALSE);
0554 
0555       ref_hist2[nh2]->SetTitle("");
0556       val_hist2[nh2]->SetTitle("");
0557 
0558       ref_hist2[nh2]->SetMarkerColor(RefCol); // rose
0559       ref_hist2[nh2]->Draw();
0560       ref_prof[npi]->SetLineColor(41);
0561       ref_prof[npi]->Draw("same");
0562         
0563       val_hist2[nh2]->SetMarkerColor(ValCol);
0564       val_hist2[nh2]->Draw("same");
0565       val_prof[npi]->SetLineColor(45);
0566       val_prof[npi]->Draw("same");
0567 
0568       leg->AddEntry(ref_prof[npi], "CMSSW_" + ref_vers, "pl");
0569       leg->AddEntry(val_prof[npi], "CMSSW_" + val_vers, "pl");
0570 
0571       leg->Draw("");
0572       
0573       myc->SaveAs(OutLabel);
0574 
0575       npi++;
0576       nh2++;
0577       i++;
0578     } else if (DimSwitch == "2D") {
0579 
0580       myc->SetGrid(0, 0);
0581 
0582       //Get histograms from files
0583       ref_file.cd(RefHistDir);
0584       ref_hist2D[n2D] = (TH2F*) gDirectory->Get(HistName);
0585         
0586       val_file.cd(ValHistDir);
0587       val_hist2D[n2D] = (TH2F*) gDirectory->Get(HistName);
0588 
0589       ref_hist2D[n2D]->SetStats(kFALSE);
0590       val_hist2D[n2D]->SetStats(kFALSE);
0591 
0592       // HACK to change what is embedded in DQM histos
0593       ref_hist2D[n2D]->GetXaxis()->SetLabelSize(0.04);
0594       val_hist2D[n2D]->GetXaxis()->SetLabelSize(0.04);
0595       ref_hist2D[n2D]->GetYaxis()->SetLabelSize(0.04);
0596       val_hist2D[n2D]->GetYaxis()->SetLabelSize(0.04);
0597       ref_hist2D[n2D]->GetXaxis()->SetTitleSize(0.045);
0598       val_hist2D[n2D]->GetXaxis()->SetTitleSize(0.045);
0599 
0600       ref_hist2D[n2D]->GetXaxis()->SetTickLength(-0.015);
0601       val_hist2D[n2D]->GetXaxis()->SetTickLength(-0.015);
0602       ref_hist2D[n2D]->GetYaxis()->SetTickLength(-0.015);
0603       val_hist2D[n2D]->GetYaxis()->SetTickLength(-0.015);
0604 
0605       ref_hist2D[n2D]->GetXaxis()->SetLabelOffset(0.02);
0606       val_hist2D[n2D]->GetXaxis()->SetLabelOffset(0.02);
0607       ref_hist2D[n2D]->GetYaxis()->SetLabelOffset(0.02);
0608       val_hist2D[n2D]->GetYaxis()->SetLabelOffset(0.02);
0609 
0610       ref_hist2D[n2D]->GetXaxis()->SetTitleOffset(1.0);
0611       val_hist2D[n2D]->GetXaxis()->SetTitleOffset(1.0);
0612 
0613       ref_hist2D[n2D]->SetTitle("");
0614       val_hist2D[n2D]->SetTitle("");
0615 
0616       // special zoom on HB/HE depth1
0617       if (n2D == 1) {
0618     ref_hist2D[n2D]->GetXaxis()->SetRangeUser(-29., 28.);
0619     val_hist2D[n2D]->GetXaxis()->SetRangeUser(-29., 28.);
0620       }
0621       
0622       //AF
0623       //Min/Max Convetion: Default AxisMin = 0. Default AxisMax = -1.
0624       //xAxis
0625       //if (xAxisMax > 0 || xAxisMin != 0) {
0626       ref_hist2D[n2D]->GetXaxis()->SetRangeUser(xAxisMin, xAxisMax);
0627       val_hist2D[n2D]->GetXaxis()->SetRangeUser(xAxisMin, xAxisMax);
0628       //}
0629       //yAxis
0630       //if (yAxisMax > 0 || yAxisMin != 0) {
0631       ref_hist2D[n2D]->GetYaxis()->SetRangeUser(yAxisMin, yAxisMax);
0632       val_hist2D[n2D]->GetYaxis()->SetRangeUser(yAxisMin, yAxisMax);
0633       //}
0634 
0635       TLegend *leg1 = new TLegend(0.50, 0.91, 0.84, 0.99, "", "brNDC");
0636       leg1->SetBorderSize(2);
0637       leg1->SetFillStyle(1001);
0638       leg1->AddEntry(ref_hist2D[n2D], "CMSSW_" + ref_vers, "l");
0639 
0640       if (xTitleCheck != "NoTitle") ref_hist2D[n2D]->GetXaxis()->SetTitle(xAxisTitle);
0641       ref_hist2D[n2D]->Draw("colz");
0642       leg1->Draw();
0643       myc->SaveAs("ref_" + OutLabel);
0644 
0645 
0646       TLegend *leg2 = new TLegend(0.50, 0.91, 0.84, 0.99, "", "brNDC");
0647       leg2->SetBorderSize(2);
0648       leg2->SetFillStyle(1001);
0649       leg2->AddEntry(val_hist2D[n2D], "CMSSW_" + val_vers, "l");
0650 
0651       if (xTitleCheck != "NoTitle") val_hist2D[n2D]->GetXaxis()->SetTitle(xAxisTitle);
0652       val_hist2D[n2D]->Draw("colz");
0653       leg2->Draw();
0654       myc->SaveAs("val_" + OutLabel);
0655       
0656       n2D++;
0657     }
0658 
0659 
0660     if (myc) delete myc;
0661     if (leg) delete leg;
0662     if (ptchi2) delete ptchi2;
0663     if (ptstats_r) delete ptstats_r;
0664     if (ptstats_v) delete ptstats_v;
0665   }
0666   return;
0667 }
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.2.1 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!