void lasers() {
  gROOT->SetStyle("Plain");

  TCanvas *c1=new TCanvas("c1","Canvas",600,600);
  gStyle->SetOptStat(0);

  for(unsigned i(0);i<2;i++) {
    if(i==0) TFile::Open("MpsAnalysisLaserThresholdScan470873.root");
    if(i==1) TFile::Open("MpsAnalysisLaserThresholdScan470853.root");

    X021Y126Profile->SetTitle(";Threshold (DTU);Number of hits/4 DTU");

    //Trim0Quad0Means->SetLineColor(2);
    X021Y126Profile->SetLineWidth(2);
    X021Y126Profile->GetYaxis()->SetTitleOffset(1.4);
    //X021Y126Profile->GetXaxis()->SetRangeUser(-100,100);
    X021Y126Profile->Draw();

    if(i==1) {
      TF1 *erfSig=new TF1("erfSig","0.5*[0]*(1.0+TMath::Erf(([1]-x)/(sqrt(2.0)*[2])))",-10000,10000);
      erfSig->SetLineColor(2);
      erfSig->SetParameters(100,160,10);
      X021Y126Profile->Fit("erfSig","","",50,300);
    }

    c1->Update();
    if(i==0) c1->Print("../lasers1.png");
    if(i==1) c1->Print("../lasers2.png");
  }
}