plot_raw(){


  TCut XX = "flag<=2";   // flag=2 intime; 0 not intime
  
  TCanvas * c1 =  new TCanvas("c1","raw data: elementID", 10, 20, 800, 800);

  TCanvas * c2 =  new TCanvas("c2","raw data: TDCtime", 400, 200, 800, 800);

  TCanvas * c3 =  new TCanvas("c3","raw data: Global Information", 800, 400, 800, 800);

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

  ///DP stat01
  c1->cd(1);
  save->Draw("elementID","detectorID==55"+XX);
  
  c1->cd(2);
  save->Draw("elementID","detectorID==56"+XX);

  c1->cd(3);
  save->Draw("elementID","detectorID==57"+XX);

  c1->cd(4);
  save->Draw("elementID","detectorID==58"+XX);


  //stat-2
    c1->cd(5);
    save->Draw("elementID","detectorID==59"+XX);

    c1->cd(6);
    save->Draw("elementID","detectorID==60"+XX);

    c1->cd(7);
    save->Draw("elementID","detectorID==61"+XX);

    c1->cd(8);
    save->Draw("elementID","detectorID==62"+XX);


  //TDC time
    c2->Divide(2,4);
    
    c2->cd(1);
    save->Draw("tdcTime","detectorID==55"+XX);

    c2->cd(2);
    save->Draw("tdcTime","detectorID==56"+XX);
    
    c2->cd(3);
    save->Draw("tdcTime","detectorID==57"+XX);

    c2->cd(4);
    save->Draw("tdcTime","detectorID==58"+XX);



    //-- c3 global informaiton

    c3->cd();
    c3->Divide(3,3);


    c3->cd(1);
    save->Draw("fRunID");

    c3->cd(2);
    save->Draw("fEventID");

    c3->cd(3);
    save->Draw("fSpillID");

    c3->cd(4);
    save->Draw("fTriggerBits");
    
    c3->cd(5);
    save->Draw("fTargtePos");

    c3->cd(6);
    save->Draw("fIntensity[0]");  //intensity[33]

    c3->cd(7);
    save->Draw("fNRoads[1]");  // fNRoads[4]

    c3->cd(8);
    save->Draw("fNHits[55]");  // Hits[63]  // DP triggers = DetectorID = 55-62

}