void DumpPlots_Digis(){ gStyle->SetOptStat(0); //gStyle->SetOptTitle(0); TCanvas *c1 = new TCanvas("c1", "c1", 600, 600); c1->Range(0, 0, 1, 1); c1->SetFillColor(0); c1->SetBorderMode(0); c1->SetBorderSize(2); c1->SetFrameBorderMode(0); char filename[200], dirname[200], histname[200]; int run = 1; int adcs= 2; const int layers = 4; const int sensorsize = 128; const int Sensor_Iu = 0; const int Sensor_Iv = 0; HGCalTBTopology IsCellValid; TFile *F; for(int iii=1; iii<=run; iii++){ for(int kkk=1; kkk<= layers; kkk++){ for(int jjj=0;jjj<adcs;jjj++){ sprintf(filename,"test_DigiAndRechitPlotter_TB_Run%i.root",iii); F = TFile::Open(filename); F->cd("hgcaltbdigisplotter"); sprintf(histname,"FullLayer_ADC%i_Layer%i",jjj,kkk); TH2Poly* Overview_Digis_Layer = (TH2Poly*) F->FindObjectAny(histname); Overview_Digis_Layer->Scale(1./512,""); Overview_Digis_Layer->GetYaxis()->SetTitle("Y[cm]"); Overview_Digis_Layer->GetYaxis()->CenterTitle(1); Overview_Digis_Layer->Draw("colztext"); ReverseXAxis(Overview_Digis_Layer); sprintf(dirname,"/home/rchatter/shervinTest/CMSSW_7_6_3_patch2/src/HGCal/Digis/Run%i/FullLayer_Digis_ADC%i_Layer%i.png",iii,jjj,kkk); c1->SaveAs(dirname); sprintf(histname,"FullLayer_ADC%i_Layer%i_summed",jjj,kkk); TH1F* Overview_Digis_Layer_Summed = (TH1F*) F->FindObjectAny(histname); Overview_Digis_Layer_Summed->Draw(); sprintf(dirname,"/home/rchatter/shervinTest/CMSSW_7_6_3_patch2/src/HGCal/Digis/Run%i/FullLayer_Digis_ADC%i_Layer%i_Summed.png",iii,jjj,kkk); c1->SaveAs(dirname); sprintf(histname,"FullLayer_ADC%i_Layer%i_profile",jjj,kkk); TProfile* Overview_Digis_Layer_Profile = (TProfile*) F->FindObjectAny(histname); Overview_Digis_Layer_Profile->Draw(); sprintf(dirname,"/home/rchatter/shervinTest/CMSSW_7_6_3_patch2/src/HGCal/Digis/Run%i/FullLayer_Digis_ADC%i_Layer%i_Profile.png",iii,jjj,kkk); c1->SaveAs(dirname); }// loop over adcs ends here for(int iv = -7; iv < 8; iv++){ for(int iu = -7; iu < 8; iu++){ for(int jjj=0;jjj<adcs;jjj++){ if(!IsCellValid.iu_iv_valid(layers, Sensor_Iu, Sensor_Iv, iu, iv, sensorsize)) continue; sprintf(histname,"Cell_u_%i_v_%i_ADC%i_Layer%i",iu,iv,jjj,kkk); TH1F* Cell_Hist_U_V = (TH1F*) F->FindObjectAny(histname); sprintf(dirname,"/home/rchatter/shervinTest/CMSSW_7_6_3_patch2/src/HGCal/Digis/Run%i/Detailed/Cell_u_%i_v_%i_ADC%i_Layer_%i.png",iii,iu,iv,jjj,kkk); Cell_Hist_U_V->Draw(); c1->SaveAs(dirname); } }//loop over iu ends here }//loop over iv ends here }//loop over layer ends here }//loop over run ends her }
void DrawEff_DS(){ gStyle->SetOptStat(kFALSE); gStyle->SetErrorX(0); //bool isweight=kFALSE; bool isweight=kTRUE; bool isHF=kTRUE; //bool isHF=kFALSE; if(isweight) TString dir="weight"; else TString dir="noweight"; const int Nhist=4; const int NGen=1; const int Ncent=6; //const double HFbin[]={0,0.1,0.2,0.4,0.6,0.8,1,2,4,6,9,12,15,20,30,50,75,100}; const double xerr[Ncent]={0,0,0,0,0,0}; /* const double HFbin_cent[Nhist][Ncent+1]={ {0,2.92,10.1,19.3,23.7,30.2,100}, {0,1.03,5.00,11.6,15.0,20.3,100}, {0,2.7,12.6,28.6,36.9,49.6,100}}; */ const double HFbin_cent[Nhist][Ncent+1]={ {0,2.725,10.065,19.565,23.975,30.455,100}, {0,0.975,5.005,11.775,15.285,20.725,100}, {0,2.515,12.465,28.855,37.235,50.085,100}}; const double cent[Ncent]={95,75,45,25,15,5}; const int nHFbin=sizeof(HFbin)/sizeof(double)-1; //double color[NGen][Nhist]={{2,1,4},{2,1,4}}; //double Style[NGen][Nhist]={{20,24,28},{20,24,28}}; double color[NGen][Nhist]={{2,1,4}}; double Style[NGen][Nhist]={{20,24,28}}; //TString Gen[NGen]={"Hijing","Epos"}; TString Gen[NGen]={"Hijing"}; TString histoname[Nhist]={"HFEnergy4","HFEnergyPlus4","HFEnergyPlus","Ntrk"}; TString etarange[Nhist]={" HF E_{T} |#eta|>4"," HF E_{T} #eta>4"," HF E_{T} #eta>0",""}; TString filename[NGen]; TFile* f[NGen]; TH1F* h_or[NGen][Nhist]; TH1F* reh_or[NGen][Nhist]; TH1F* h_tr[NGen][Nhist]; TH1F* reh_tr[NGen][Nhist]; TH1F* h_ds[NGen][Nhist]; TH1F* reh_ds[NGen][Nhist]; TH1F* h_eff[NGen][Nhist]; TH1F* reh_eff[NGen][Nhist]; TLatex* T1[Nhist]; TLatex* T2; TLatex* T3[NGen]; TF1 *fitEff[NGen][Nhist]; //TGraphAsymmErrors *gEff[NGen][Nhist]; TGraph *gEff[NGen][Nhist]; c1 = new TCanvas("c1"," ",800,500); c2 = new TCanvas("c2"," ",800,700); makeMultiPanelCanvas(c1,NGen,1,0.03,0.03,0.14,0.14,0.03); makeMultiPanelCanvas(c2,1,1,0.03,0.03,0.03,0.03,0.03); //for |eta|>4 HFEnergy: /* for(int k=0;k<NGen;k++) for(int j=0;j<Nhist;j++) fitEff[k][j]= new TF1(Form("fitEff_%s%s",Gen[k].Data(),histoname[j]),"[0]-[1]*exp([2]*x)-[3]*exp(-[4]*x)",0,300); fitEff[0][0]->SetParameters(0.999876,0.44648,-1.78911,0.242036,0.549932); fitEff[0][1]->SetParameters(1.,0.367894,-3.36629,0.104605,0.690272); fitEff[0][2]->SetParameters(0.99981,0.697501,-2.59569,0.15323,0.411757); fitEff[1][0]->SetParameters(0.999918,0.691843,-0.785042,0.0352768,0.449521); fitEff[1][1]->SetParameters(0.999929,0.284,-2.25636,0.264114,0.951221); fitEff[1][2]->SetParameters(0.999887,0.5,-2.39395,0.373215,0.550851); */ for(int k=0;k<NGen;k++){ T3[k]=new TLatex(0.45,0.93,Gen[k]); T3[k]->SetNDC(); T3[k]->SetTextAlign(12); T3[k]->SetTextSize(0.05); T3[k]->SetTextColor(1); T3[k]->SetTextFont(42); } for(int j=0;j<Nhist;j++){ T1[j]=new TLatex(0.4,0.85,etarange[j]); T1[j]->SetNDC(); T1[j]->SetTextAlign(12); T1[j]->SetTextSize(0.05); T1[j]->SetTextColor(1); T1[j]->SetTextFont(42); } T2=new TLatex(0.2,0.25,"Double Side"); T2->SetNDC(); T2->SetTextAlign(12); T2->SetTextSize(0.06); T2->SetTextColor(1); T2->SetTextFont(42); j=3; for(int k=0;k<NGen;k++){ //for(int j=0;j<Nhist;j++){ filename[k]=Form("output/%s/pPbHist_%s_DS",dir.Data(),Gen[k].Data()); f[k]=TFile::Open(Form("%s.root",filename[k].Data())); h_or[k][j]=(TH1F*)f[k]->Get(Form("h%s",histoname[j].Data())); h_tr[k][j]=(TH1F*)f[k]->Get(Form("h%s_tr",histoname[j].Data())); h_ds[k][j]=(TH1F*)f[k]->Get(Form("h%s_DS",histoname[j].Data())); h_eff[k][j]=(TH1F*)f[k]->Get(Form("hEff_%s",histoname[j].Data())); fixedFontHist(h_or[k][j],1.8,2); fixedFontHist(h_tr[k][j],1.8,2); fixedFontHist(h_ds[k][j],1.8,2); //reh_or[k][j]=(TH1F*)h_or[k][j]->Clone(Form("re%s%s",Gen[k].Data(),histoname[j].Data())); //reh_tr[k][j]=(TH1F*)h_tr[k][j]->Clone(Form("re%s%s_tr",Gen[k].Data(),histoname[j].Data())); //reh_ds[k][j]=(TH1F*)h_ds[k][j]->Clone(Form("re%s%s_ds",Gen[k].Data(),histoname[j].Data())); //reh_or[k][j]=(TH1F*)h_or[k][j]->Rebin(Ncent,Form("re%s%s",Gen[k].Data(),histoname[j].Data()),HFbin_cent[j]); //reh_tr[k][j]=(TH1F*)h_tr[k][j]->Rebin(Ncent,Form("re%s%s_tr",Gen[k].Data(),histoname[j].Data()),HFbin_cent[j]); //reh_ds[k][j]=(TH1F*)h_ds[k][j]->Rebin(Ncent,Form("re%s%s_ds",Gen[k].Data(),histoname[j].Data()),HFbin_cent[j]); reh_or[k][j]=(TH1F*)h_or[k][j]->Rebin(nHFbin,Form("re%s%s",Gen[k].Data(),histoname[j].Data()),HFbin); reh_tr[k][j]=(TH1F*)h_tr[k][j]->Rebin(nHFbin,Form("re%s%s_tr",Gen[k].Data(),histoname[j].Data()),HFbin); reh_ds[k][j]=(TH1F*)h_ds[k][j]->Rebin(nHFbin,Form("re%s%s_ds",Gen[k].Data(),histoname[j].Data()),HFbin); //normalizeByBinWidth(reh_or[k][j]); //normalizeByBinWidth(reh_tr[k][j]); //normalizeByBinWidth(reh_ds[k][j]); //gEff[k][j]=new TGraphAsymmErrors(); gEff[k][j]=new TGraph(h_eff[k][j]); gEff[k][j]->SetPoint(0,0,1); gEff[k][j]->SetName(Form("gEff_%s%s",Gen[k].Data(),histoname[j].Data())); //gEff[k][j]->BayesDivide(reh_tr[k][j],reh_ds[k][j]); gEff[k][j]->SetMarkerColor(color[k][j]); gEff[k][j]->SetLineColor(color[k][j]); gEff[k][j]->SetMarkerStyle(Style[k][j]); gEff[k][j]->SetMarkerSize(1.2); gEff[k][j]->GetXaxis()->SetLimits(0.01,99.9); gEff[k][j]->SetTitle(""); //} } TH1D * hFrame = new TH1D("","",100,0,100.); if(!isHF){ hFrame->GetXaxis()->SetLabelOffset(999); hFrame->GetXaxis()->SetTickLength(0);} hFrame->GetYaxis()->SetTitle("Event Selection Efficiency"); hFrame->SetMaximum(1.09); hFrame->SetMinimum(0); if(isHF) hFrame->GetXaxis()->SetTitle("HF E_{T} (GeV/c)"); if(j==3) hFrame->GetXaxis()->SetTitle("track Mult"); hFrame->GetYaxis()->SetTitleSize(0.05); fixedFontHist(hFrame); hFrame->GetYaxis()->SetTitleOffset(1.5); for(k=0;k<NGen;k++){ if(isHF) c1->cd(k+1)->SetLogx(); else c1->cd(k+1); TLegend *leg=new TLegend(0.60,0.2,0.75,0.4); if(k!=0){ hFrame->GetYaxis()->SetTitle(""); hFrame->GetYaxis()->SetLabelSize(0); } hFrame->DrawCopy(); //for(j=0;j<Nhist;j++){ if(isHF) gEff[k][j]->Draw("Psame"); else{ ReverseXAxis(gEff[k][j],0.05); ReverseXGraph(gEff[k][j]); } //fitEff[k][j]->SetLineStyle(color[k][j]); //fitEff[k][j]->SetLineColor(color[k][j]); //fitEff[k][j]->Draw("same"); leg->AddEntry(gEff[k][j],etarange[j],"lp"); //} leg->SetTextSize(0.05); leg->SetBorderSize(0); leg->SetFillColor(0); TLine *l=new TLine(0,1,100,1); l->SetLineStyle(2); l->Draw("same"); leg->Draw("same"); T3[k]->Draw("same"); } if(isHF) c1->Print(Form("pic/OwnBinning/pPbHist_%s_DS_HF.png",dir.Data())); else c1->Print(Form("pic/OwnBinning/pPbHist_%s_DS_cent.png",dir.Data())); //------------------------------------------------------------------------------------------ if(isHF) //c2->cd()->SetLogx(); //else c2->cd(); TLegend *leg=new TLegend(0.7,0.2,0.85,0.4); hFrame->GetYaxis()->SetTitle("Event Selection Efficiency"); hFrame->GetYaxis()->SetLabelSize(18); hFrame->GetYaxis()->SetTitleOffset(1); hFrame->DrawCopy(); for(k=0;k<NGen;k++){ if(k==0){ gEff[k][j]->SetMarkerColor(1); gEff[k][j]->SetLineColor(1); gEff[k][j]->SetMarkerStyle(20); gEff[k][j]->SetMarkerSize(1.6); } else{ gEff[k][j]->SetMarkerColor(2); gEff[k][j]->SetLineColor(2); gEff[k][j]->SetMarkerStyle(24); gEff[k][j]->SetMarkerSize(1.6); } if(isHF) gEff[k][j]->Draw("Psame"); else{ ReverseXAxis(gEff[k][j],0.05); ReverseXGraph(gEff[k][j]); } leg->SetTextSize(0.05); leg->SetBorderSize(0); leg->SetFillColor(0); leg->AddEntry(gEff[k][j],Gen[k],"p"); TLine *l=new TLine(0,1,100,1); l->SetLineStyle(2); } T1[j]->Draw("same"); leg->Draw("same"); l->Draw("same"); if(isHF) c2->Print(Form("pic/OwnBinning/pPbHist_%s_DS_%s_HF.png",dir.Data(),histoname[j].Data())); else c2->Print(Form("pic/OwnBinning/pPbHist_%s_DS_%s_cent.png",dir.Data(),histoname[j].Data())); }
void DumpPlotsReco(TString inputFileName, TString outputFolder, Int_t runNumber, Int_t nSpills) { std::cout << "DumpPlotsReco called for inputFileName: " << inputFileName << ", outputFolder: " << outputFolder << ", runNumber: " << runNumber << ", nSpills: " << nSpills << std::endl; gStyle->SetOptStat(0); //gStyle->SetOptTitle(0); TCanvas *plotCanvas = new TCanvas("plotCanvas", "plotCanvas", 600, 600); plotCanvas->Range(0, 0, 1, 1); plotCanvas->SetFillColor(0); plotCanvas->SetBorderMode(0); plotCanvas->SetBorderSize(2); plotCanvas->SetFrameBorderMode(0); const int nLayers = 8; int nSkirocs = nLayers*2; // const int nChannels = 64; int nEvents = nSpills*EVENTSPERSPILL*nLayers; TFile *inputFile = TFile::Open(inputFileName); inputFile->cd("hgcaltbrechitsplotter_highgain_correlation_cm"); TString outputPlotFullPath, plotName; plotName = Form("Noise_2D_Profile_Layer"); TH2F* h_Noise_2D_Profile_Layer = (TH2F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_Noise_2D_Profile_Layer)) { h_Noise_2D_Profile_Layer->Draw("colz"); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_Noise_2D_Profile_Layer; plotName = Form("AllCells_Ped"); TH1F* h_AllCells_Ped = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_AllCells_Ped)) { h_AllCells_Ped->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_AllCells_Ped; plotName = Form("AllCells_CM"); TH1F* h_AllCells_CM = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_AllCells_CM)) { h_AllCells_CM->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_AllCells_CM; for(int layerCounter=0; layerCounter < nLayers; ++layerCounter){ plotName = Form("Full_Cell_Layer_%d", layerCounter); TH1F* h_Full_Cell = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_Full_Cell)) { h_Full_Cell->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_Full_Cell; plotName = Form("Half_Cell_Layer_%d", layerCounter); TH1F* h_Half_Cell = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_Half_Cell)) { h_Half_Cell->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_Half_Cell; plotName = Form("MB_Cell_Layer_%d", layerCounter); TH1F* h_MB_Cell = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_MB_Cell)) { h_MB_Cell->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_MB_Cell; plotName = Form("Calib_Pads_Layer_%d", layerCounter); TH1F* h_Calib_Pads = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_Calib_Pads)) { h_Calib_Pads->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_Calib_Pads; plotName = Form("Merged_Cell_Layer_%d", layerCounter); TH1F* h_Merged_Cell = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_Merged_Cell)) { h_Merged_Cell->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_Merged_Cell; // for (int skirocCounter = 1; skirocCounter <= nSkirocs; ++skirocCounter) { // for (int channelCounter = 0; channelCounter <= nChannels; ++channelCounter) { // plotName = Form("Ski_%d_Channel_%d_Layer_%d", skirocCounter, channelCounter, layerCounter); // TH1F* h_Ski = (TH1F*) inputFile->FindObjectAny(plotName); // if (not(NULL == h_Ski)) { // h_Ski->Draw(); // outputPlotFullPath = (outputFolder + Form("/Detailed/")) + plotName + Form("_%06d.png", runNumber); // plotCanvas->SaveAs(outputPlotFullPath); // } // delete h_Ski; // } // ends loop over channels // } // ends loop over skirocs } //loop over layers ends here inputFile->cd("hgcaltbrechitsplotter_highgain_new"); plotName = Form("Covar_hist"); TH2F* h_Covar_hist = (TH2F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_Covar_hist)) { h_Covar_hist->Draw("colz"); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_Covar_hist; plotName = Form("Correl_hist"); TH2F* h_Correl_hist = (TH2F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_Correl_hist)) { h_Correl_hist->Draw("colz"); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_Correl_hist; plotName = Form("DiffIJ_hist"); TH2F* h_DiffIJ_hist = (TH2F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_DiffIJ_hist)) { h_DiffIJ_hist->Draw("colz"); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_DiffIJ_hist; plotName = Form("CG_X"); TH1F* h_CG_X = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_CG_X)) { h_CG_X->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_CG_X; plotName = Form("CG_Y"); TH1F* h_CG_Y = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_CG_Y)) { h_CG_Y->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_CG_Y; // plotName = Form("AllCells_Ped"); // TH1F* h_AllCells_Ped = (TH1F*) inputFile->FindObjectAny(plotName); // if (not(NULL == h_AllCells_Ped)) { // h_AllCells_Ped->Draw(); // outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); // plotCanvas->SaveAs(outputPlotFullPath); // } // delete h_AllCells_Ped; // plotName = Form("AllCells_CM"); // TH1F* h_AllCells_CM = (TH1F*) inputFile->FindObjectAny(plotName); // if (not(NULL == h_AllCells_CM)) { // h_AllCells_CM->Draw(); // outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); // plotCanvas->SaveAs(outputPlotFullPath); // } // delete h_AllCells_CM; plotName = Form("Sum_Cluster_ADC"); TH1F* h_Sum_Cluster_ADC = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_Sum_Cluster_ADC)) { h_Sum_Cluster_ADC->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_Sum_Cluster_ADC; plotName = Form("Sum_Cluster_Max"); TH1F* h_Sum_Cluster_Max = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_Sum_Cluster_Max)) { h_Sum_Cluster_Max->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_Sum_Cluster_Max; for(int layerCounter=1; layerCounter <= nLayers; ++layerCounter){ plotName = Form("FullLayer_RecHits_Layer%d_Summed", layerCounter); TH1F* h_FullLayer_RecHits = (TH1F*) inputFile->FindObjectAny(plotName); if (not(NULL == h_FullLayer_RecHits)) { h_FullLayer_RecHits->Draw(); outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_FullLayer_RecHits; } //loop over layers ends here for (int eventCenturyCounter = 0; eventCenturyCounter < int(nEvents/100); ++eventCenturyCounter) { int eventNumber = eventCenturyCounter*100; int layerNumber = 1+((int(eventNumber/EVENTSPERSPILL))%nLayers); plotName = Form("FullLayer_ADC0_Layer%d_Event%d", layerNumber, eventNumber); TH2Poly* h_EventDisplay = (TH2Poly*) inputFile->FindObjectAny(plotName); if (not(NULL == h_EventDisplay)) { h_EventDisplay->GetYaxis()->SetTitle("Y[cm]"); h_EventDisplay->GetYaxis()->CenterTitle(1); h_EventDisplay->Draw("colztext"); ReverseXAxis(h_EventDisplay); eventNumber = EVENTSPERSPILL*int(eventNumber/(EVENTSPERSPILL*nLayers)) + eventNumber%EVENTSPERSPILL; plotName = Form("EventDisplay_Event%d_Layer%d", eventNumber, layerNumber); outputPlotFullPath = (outputFolder + Form("/Detailed/")) + plotName + Form("_%06d.png", runNumber); plotCanvas->SaveAs(outputPlotFullPath); } delete h_EventDisplay; } // loop over events ends here delete inputFile; delete plotCanvas; // }//loop over spill ends here // }//loop over run ends here // OLD: // char filename[50], dirname[100], histname[50]; // int run = 2; // int spill = 3; // for(int iii=0; iii<=run; iii++){ // for(int jjj=1; jjj<=spill;jjj++){ // if(iii != 0) sprintf(filename,"test_DigiAndRechitPlotter_TB_Run%i_Spill%i.root",iii,jjj); // else sprintf(filename,"test_DigiAndRechitPlotter_TB_Cumulative.root"); // const int sensorsize = 128; // const int Sensor_Iu = 0; // const int Sensor_Iv = 0; // HGCalTBTopology IsCellValid; // plotName = Form("FullLayer_RecHits_Layer%i",layerCounter); // // sprintf(histname,"FullLayer_RecHits_Layer%i",layerCounter); // TH2Poly* Overview_Rechits_Layer = (TH2Poly*) inputFile->FindObjectAny(plotName); // Overview_Rechits_Layer->GetYaxis()->SetTitle("Y[cm]"); // Overview_Rechits_Layer->GetYaxis()->CenterTitle(1); // Overview_Rechits_Layer->Draw("colztext"); // ReverseXAxis(Overview_Rechits_Layer); // // if(iii != 0) sprintf(dirname,"/afs/cern.ch/work/r/rchatter/CMSSW_7_6_3_patch2/src/HGCal/DQM_Plots/Run%i/Spill%i/Overview/FullLayer_RecHits_Layer%i.png",iii,jjj,layerCounter); // // else sprintf(dirname,"/afs/cern.ch/work/r/rchatter/CMSSW_7_6_3_patch2/src/HGCal/DQM_Plots/Cumulative/Overview/FullLayer_RecHits_Layer%i.png",layerCounter); // outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); // plotCanvas->SaveAs(outputPlotFullPath); // delete Overview_Rechits_Layer; // plotName = Form("FullLayer_Occupancy_Layer%i", layerCounter); // TH2Poly* Overview_Occupancy_Layer = (TH2Poly*) inputFile->FindObjectAny(plotName); // Overview_Occupancy_Layer->GetYaxis()->SetTitle("Y[cm]"); // Overview_Occupancy_Layer->GetYaxis()->CenterTitle(1); // Overview_Occupancy_Layer->Draw("colztext"); // ReverseXAxis(Overview_Occupancy_Layer); // // if(iii != 0) sprintf(dirname,"/afs/cern.ch/work/r/rchatter/CMSSW_7_6_3_patch2/src/HGCal/DQM_Plots/Run%i/Spill%i/Overview/FullLayer_Occupancy_Layer%i.png",iii,jjj,layerCounter); // // else sprintf(dirname,"/afs/cern.ch/work/r/rchatter/CMSSW_7_6_3_patch2/src/HGCal/DQM_Plots/Cumulative/Overview/FullLayer_Occupancy_Layer%i.png",layerCounter); // outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); // plotCanvas->SaveAs(outputPlotFullPath); // delete Overview_Occupancy_Layer; // // plotName = Form("FullLayer_RecHits_Layer%i_Summed",layerCounter); // // TH1F* Overview_Rechits_Layer_Summed = (TH1F*) inputFile->FindObjectAny(plotName); // // Overview_Rechits_Layer_Summed->Draw(); // // // if(iii != 0) sprintf(outputPlotFullPath,"/afs/cern.ch/work/r/rchatter/CMSSW_7_6_3_patch2/src/HGCal/DQM_Plots/Run%i/Spill%i/Overview/FullLayer_RecHits_Layer%i_Summed.png",iii,jjj,layerCounter); // // // else sprintf(outputPlotFullPath,"/afs/cern.ch/work/r/rchatter/CMSSW_7_6_3_patch2/src/HGCal/DQM_Plots/Cumulative/Overview/FullLayer_RecHits_Layer%i_Summed.png",layerCounter); // // outputPlotFullPath = (outputFolder + Form("/")) + plotName + Form("_%06d.png", runNumber); // // plotCanvas->SaveAs(outputPlotFullPath); // // delete Overview_Rechits_Layer_Summed; // // // for(int iv = -7; iv < 8; iv++){ // // // for(int iu = -7; iu < 8; iu++){ // // // if(!IsCellValid.iu_iv_valid(nLayers, Sensor_Iu, Sensor_Iv, iu, iv, sensorsize)) continue; // // // plotName = Form("Cell_RecHits_u_%i_v_%i_Layer%i",iu,iv,layerCounter); // // // TH1F* Cell_Hist_U_V = (TH1F*) inputFile->FindObjectAny(plotName); // // // // if(iii != 0) sprintf(outputPlotFullPath,"/afs/cern.ch/work/r/rchatter/CMSSW_7_6_3_patch2/src/HGCal/DQM_Plots/Run%i/Spill%i/Detailed/Cell_RecHits_u_%i_v_%i_Layer_%i.png",iii,jjj,iu,iv,layerCounter); // // // // else sprintf(outputPlotFullPath,"/afs/cern.ch/work/r/rchatter/CMSSW_7_6_3_patch2/src/HGCal/DQM_Plots/Cumulative/Detailed/Cell_RecHits_u_%i_v_%i_Layer_%i.png",iu,iv,layerCounter); // // // outputPlotFullPath = (outputFolder + Form("/Detailed/")) + plotName + Form("_%06d.png", runNumber); // // // Cell_Hist_U_V->Draw(); // // // plotCanvas->SaveAs(outputPlotFullPath); // // // delete Cell_Hist_U_V; // // // }//loop over iu ends here // // // }//loop over iv ends here }