C   104

realNodpwMip

Guest on 26th June 2022 03:56:18 AM

  1. void realNodpwMip() {
  2.   bool leg(false);
  3.  
  4.   gROOT->SetStyle("Plain");
  5.  
  6.   TCanvas *c1=new TCanvas("c1","Canvas",800,600);
  7.   gStyle->SetOptStat(0);
  8.   gStyle->SetPalette(1);
  9.   //gStyle->SetPadLeftMargin(0.12);
  10.   //gStyle->SetPadRightMargin(0.12);
  11.   //c1->UseCurrentStyle();
  12.  
  13.   TFile::Open("summary_REAL_NODPW.root");
  14.  
  15.   hAmm->SetTitle(";MIP x position (#mum);MIP y position (#mum);Laser signal (DTU)");
  16.   for(int x(1);x<31;x++) {
  17.     for(int y(1);y<31;y++) {
  18.       double b=500.0*hAmm->GetBinContent(x,y);
  19.       if(b<0.01) b=0.01;
  20.       hAmm->SetBinContent(x,y,b);
  21.     }
  22.   }
  23.  
  24.   if(leg) {
  25.     hAmm->GetXaxis()->SetTitleOffset(1.4);
  26.     hAmm->GetYaxis()->SetTitleOffset(1.5);
  27.     hAmm->GetZaxis()->SetTitleOffset(1.3);
  28.     hAmm->GetZaxis()->SetRangeUser(0,287.0);
  29.   } else {
  30.     hAmm->GetXaxis()->SetTitleOffset(1.1);
  31.     hAmm->GetYaxis()->SetTitleOffset(1.1);
  32.     hAmm->GetZaxis()->SetRangeUser(0,299.0);
  33.   }
  34.  
  35.   //Trim0Sigmas2D->GetXaxis()->SetRangeUser(0,83);
  36.   //Trim0Sigmas2D->GetZaxis()->SetRangeUser(0,15);
  37.  
  38.   if(leg) hAmm->Draw("lego2");
  39.   else    hAmm->Draw("colz");
  40.  
  41.   c1->Update();
  42.   c1->Print("../realNodpwMip.png");
  43. }

Raw Paste


Login or Register to edit or fork this paste. It's free.