//
//--- plot DP alignment distributions 
//
//
/*
root [1] save->Show(0)
======> EVENT:0
 quad            = 0                   //0, 1, 2, 3 stands for top-left, top-right, bottom-left, bottom-right 
 fElementID      = 9                   //nearest station-1 hit elementID
 bElementID      = 9                   //nearest station-2 hit elementID
 fElementID_exp  = 10                  //expected station-1 hit elementID
 bElementID_exp  = 10                  //expected station-2 hit elementID
 
 // nearest station-1 hit position  -- relative to the center of the detector panel (yPos - YPos_center), y-offset = 3" ~10cm
 fPos            = -31.3569            
 fPos_exp        = -31.6229            // expected station-1 hit position
 bPos            = -29.654              // nearest station-2 hit position
 bPos_exp        = -30.25               // expected station-2 hit position

 tx              = 0.0419373            //track slope in X-Z
 ty              = 0.00957994            //track slope in Y-Z
 x0              = -20.3776              //track intersection in X-Z
 y0              = 9.55376               //track intersection in Y-Z
*/
//---------
{

  TCut XX1 = "bElementID<100 && fElementID<100 && abs((1.6-y0)/ty)<1000";
  TCut XX2 = "abs(fPos)<100 && abs(bPos)<100";

  TCanvas * c1 = new TCanvas("c1", "DP Alignment Data", 10, 20, 600, 800);
  TCanvas * c2 = new TCanvas("c2", "DP Alignment Data", 600, 100, 600, 800);

  TCanvas * c3 = new TCanvas("c3", "DP Alignment Data", 800, 200, 600, 800);

  TCanvas * c4 = new TCanvas("c4", "DP Alignment Data, Zvtx", 1000, 400, 600, 800);

 
  //--
  c1->cd();

  c1->Divide(4,4);

  c1->cd(1);
  save->Draw("quad",XX1+XX2);

  c1->cd(2);
  save->Draw("fElementID",XX1+XX2);

  c1->cd(3);
  save->Draw("bElementID",XX1+XX2);

  c1->cd(4);
  save->Draw("fPos",XX1+XX2);

  c1->cd(5);
  save->Draw("bPos",XX1+XX2);

 
  c1->cd(6);
  save->Draw("fPos:fElementID",XX1+XX2,"colz");

  c1->cd(7);
  save->Draw("bPos:bElementID",XX1+XX2,"colz");

  c1->cd(8);
  save->Draw("(1.6-y0)/ty",XX1+XX2+"abs((1.6-y0)/ty)<500");

  
  c1->cd(9);
  save->Draw("fElementID:bElementID",XX1+XX2,"colz");

  c1->cd(10);
  save->Draw("fElementID_exp",XX1+XX2);

  c1->cd(11);
  save->Draw("bElementID_exp",XX1+XX2);

  c1->cd(12);
  save->Draw("x0",XX1+XX2);
  
  c1->cd(13);
  save->Draw("y0",XX1+XX2);
   
  c1->cd(14);
  save->Draw("tx",XX1+XX2);
  
  c1->cd(15);
  save->Draw("ty",XX1+XX2);

  c1->cd(16)->SetLogy();
  save->Draw("fPos-fPos_exp",XX1+XX2);
  
   //--
   c2->cd();
   c2->Divide(3,3);

   c2->cd(1);
   save->Draw("y0:x0",XX1+XX2,"colz");

   c2->cd(2);
   save->Draw("y0:x0",XX1+XX2+"quad==0","colz");

   c2->cd(3);
   save->Draw("y0:x0",XX1+XX2+"quad==1","colz");

   c2->cd(4);
   save->Draw("y0:x0",XX1+XX2+"quad==2","colz");

   c2->cd(5);
   save->Draw("y0:x0",XX1+XX2+"quad==3","colz");

   c2->cd(6)->SetLogy();
   save->Draw("fElementID-fElementID_exp",XX1+XX2+"abs(fElementID-fElementID_exp)<10");

   c2->cd(7)->SetLogy();
   save->Draw("bElementID-bElementID_exp",XX1+XX2+"abs(bElementID-bElementID_exp)<10");

   c2->cd(8);
   save->Draw("(fElementID-fElementID_exp):quad",XX1+XX2+"abs(fElementID-fElementID_exp)<10","colz");

   c2->cd(9);
   save->Draw("(bElementID-bElementID_exp):quad",XX1+XX2+"abs(bElementID-bElementID_exp)<5","colz");


   //--
   c3->cd();
   c3->Divide(2,2);

   c3->cd(1);
   save->Draw("(fPos+40+10):fElementID","fElementID<100 && fPos <100 && quad==0","colz");

   c3->cd(2);
   save->Draw("(bPos+50+10):bElementID","bElementID<100 && bPos <100 && quad==0","colz");

   c3->cd(3);
   save->Draw("(fPos-40-10):fElementID","fElementID<100 && fPos <100 && quad==3","colz");

   c3->cd(4);
   save->Draw("(bPos-50-10):bElementID","bElementID<100 && bPos <100 && quad==3","colz");


   //-- z-vtx distributions

   c4->cd();
   c4->Divide(4,2);

   c4->cd(1);
   save->Draw("(1.6-y0)/ty",XX1+XX2+"quad==0");
   
   c4->cd(2);
   save->Draw("(1.6-y0)/ty",XX1+XX2+"quad==1");

   c4->cd(3);
   save->Draw("(1.6-y0)/ty",XX1+XX2+"quad==2");

   c4->cd(4);
   save->Draw("(1.6-y0)/ty",XX1+XX2+"quad==3");
   
   c4->cd(5);
   save->Draw("ty",XX1+XX2+"quad==0");
   
   c4->cd(6);
   save->Draw("ty",XX1+XX2+"quad==1");

   c4->cd(7);
   save->Draw("ty",XX1+XX2+"quad==2");

   c4->cd(8);
   save->Draw("ty",XX1+XX2+"quad==3");

}