void DrawDmDiff2(TString path="/panda/pandaroot/macro/lmd/testPixel/mom_1_5/"){
  //  double dP[11]={-1e-1,-1e-2,-1e-3,-1e-4,-1e-5,0,1e-5,1e-4,1e-3,1e-2,1e-1};
  TString dPstr[11]={"0","-1e-1", "-1e-2", "-1e-3", "-1e-4", "-1e-5", "1e-5", "1e-4", "1e-3", "1e-2", "1e-1"};
  TH2D *hthphLMD[11]; TH2D *hthphIP[11]; 
  TH2D *hthIPthLMD[11]; TH2D *hphIPphLMD[11]; 
  TCanvas *c1 = new TCanvas();
  TString name_out = path+"/PlotsAccSummary_noCuts_PrecIsPbeam.pdf";
  name_out_o = name_out+"(";
  name_out_c = name_out+")";
  TLegend *leg = new TLegend(0.75,0.6,0.98,0.99);
  leg->SetHeader("P_{sim}=P_{beam}+dP, P_{rec}=P_{beam}");
  c1->Divide(2,2);
  for(int i=0;i<11;i++){
    TString path_loc = path+"ShiftMom_";
    path_loc+=dPstr[i];
    TString file_in = path_loc+"/Lumi_EffTrk_0.root";
    cout<<file_in<<endl;
    TFile *f_in = new TFile(file_in.Data(),"READ");
    if(f_in==0) continue;
    TH2 *hresPhi_phi = (TH2D*)f_in->Get("hresPhiPhi");
    if(hresPhi_phi==0) continue;
    // hthphLMD[i] = new TH2D();     hthphIP[i] = new TH2D();
    hthphLMD[i] = (TH2D*)f_in->Get("hphthLMD");
    TString nameAcc = "dP =";
    nameAcc+=dPstr[i];
    hthphLMD[i]->SetTitle(nameAcc.Data());
    c1->cd(1);
    if(i>0) hthphLMD[i]->Divide(hthphLMD[0]);
    hthphLMD[i]->SetStats(kFALSE);
    hthphLMD[i]->Draw("colz");
    //    c1->Print(name_out_o);
    hthphIP[i] =  (TH2D*)f_in->Get("hphthIP");
    hthphIP[i]->SetTitle(nameAcc.Data());
    c1->cd(2);
    if(i>0) hthphIP[i]->Divide(hthphIP[0]);
    hthphIP[i]->SetStats(kFALSE);
    hthphIP[i]->Draw("colz");

    hthIPthLMD[i] = (TH2D*)f_in->Get("hthIPthLMD");
    hthIPthLMD[i]->SetTitle(nameAcc.Data());
    c1->cd(3);
    if(i>0) hthIPthLMD[i]->Divide(hthIPthLMD[0]);
    hthIPthLMD[i]->SetStats(kFALSE);
    hthIPthLMD[i]->Draw("colz");

    hphIPphLMD[i] = (TH2D*)f_in->Get("hphIPphLMD");
    hphIPphLMD[i]->SetTitle(nameAcc.Data());
    c1->cd(4);
    if(i>0) hphIPphLMD[i]->Divide(hphIPphLMD[0]);
    hphIPphLMD[i]->SetStats(kFALSE);
    hphIPphLMD[i]->Draw("colz");
    c1->Print(name_out_o);
    // if(i==10) c1->Print(name_out_c);
    // else c1->Print(name_out_o);
    //   leg->AddEntry(hresPhi_phi_mean[i],Pleg.Data(),"ep");
    // f_in->Close();
  }
  // leg->Draw();
  c1->Clear();
  c1->Print(name_out_c);




}