{
//////////////////////////////////////////////////
// tof.C : Root Macro
//
//  Example of the Hall D Time-of-Flight
//
// The reaction gamma p -> p K- K+ pi- pi+ was
// generated using "genr8" and processed through
// HDFast to create the data file test1.rdt 
//  
//*-- Author :    Paul Eugenio 1-Jan-99
//*-- CMZ : PME 14-Jan-99
//////////////////////////////////////////////////

//
//Begin_Html
/*
<a href= "smearedTOF.gif" >
Here is the output using smeared momentum in the calculation</a>
<a href= "tracedTOF.gif" >
Here is the output for using the trace point momentum   </a>

*/
//End_Html

gROOT->Reset();
gSystem->Load("libTMCFast.so"); 

// This file was generated by HDFast
TFile f("test1.rdt");

TTree *tree = (TTree*)f.Get("T");

TBranch *b1 = tree->GetBranch("hepevt");
TBranch *b2 = tree->GetBranch("tof_trace");
TBranch *b3 = tree->GetBranch("offtrk");

TMCFastHepEvt *hepevt;
TMCFastTOF *trace; 
TMCFastOfflineTrack *offtrk;  
hepevt = new TMCFastHepEvt();
trace = new TMCFastTOF();
offtrk = new TMCFastOfflineTrack();

b1->SetAddress(&hepevt);
b2->SetAddress(&trace);
b3->SetAddress(&offtrk);

// create some histograms and some place to put them

TH1F *h1 = new TH1F("h1", "CTOF mass kaon",100,0.1,0.8);
TH1F *h2 = new TH1F("h2", "CTOF mass pion",100,0.0,0.5);
TH1F *h3 = new TH1F("h3", "CTOF mass proton",100,0.2,1.2);
TH1F *h4 = new TH1F("h4", "FTOF mass kaon",100,0.1,0.8);
TH1F *h5 = new TH1F("h5", "FTOF mass pion",100,0.0,0.5);
TH1F *h6 = new TH1F("h6", "FTOF mass proton",100,0.2,1.2);

// now create some pads to put the histo's
c1 = new TCanvas("c1","Halld Time-of-Flight",100,10,700,900);
c1->SetFillColor(18);

pad1 = new TPad("pad1","The pad with the function",0.05,0.66,0.45,0.96,21);
pad2 = new TPad("pad2","The pad with the histogram",0.05,0.34,0.45,0.64,21);
pad3 = new TPad("pad3","The pad with the histogram",0.05,0.02,0.45,0.32,21);
pad4 = new TPad("pad4","The pad with the function",0.55,0.66,0.95,0.96,21);
pad5 = new TPad("pad5","The pad with the histogram",0.55,0.34,0.95,0.64,21);
pad6 = new TPad("pad6","The pad with the histogram",0.55,0.02,0.95,0.32,21);

pad1->Draw();
pad2->Draw();
pad3->Draw();
pad4->Draw();
pad5->Draw();
pad6->Draw();

Int_t nentries = (Int_t)tree->GetEntries();

for (Int_t ev = 0; ev < nentries; ev++) {
  
  tree->GetEvent(ev);        //read  event in memory

    /**** Using the smeared momentum *
    h1->Fill(trace->CTOFmass(1,*offtrk));// kaon
    h1->Fill(trace->CTOFmass(4,*offtrk));// kaon
    h2->Fill(trace->CTOFmass(2,*offtrk));// pion
    h2->Fill(trace->CTOFmass(3,*offtrk));// pion
    h3->Fill(trace->CTOFmass(5,*offtrk));// proton
    h4->Fill(trace->FTOFmass(1,*offtrk));// kaon
    h4->Fill(trace->FTOFmass(4,*offtrk));// kaon
    h5->Fill(trace->FTOFmass(2,*offtrk));// pion
    h5->Fill(trace->FTOFmass(3,*offtrk));// pion
    h6->Fill(trace->FTOFmass(5,*offtrk));// proton
    *******/

  /******* Using the trace point momentum */
    h1->Fill(trace->CTOFmass(1));// kaon
    h1->Fill(trace->CTOFmass(4));// kaon
    h2->Fill(trace->CTOFmass(2));// pion
    h2->Fill(trace->CTOFmass(3));// pion
    h3->Fill(trace->CTOFmass(5));// proton
    h4->Fill(trace->FTOFmass(1));// kaon
    h4->Fill(trace->FTOFmass(4));// kaon
    h5->Fill(trace->FTOFmass(2));// pion
    h5->Fill(trace->FTOFmass(3));// pion
    h6->Fill(trace->FTOFmass(5));// proton
    /************/

  hepevt->Clear();
  trace->Clear();
  offtrk->Clear();
}

// Now draw the histograms
pad1->cd();
h1->Draw();
pad2->cd();
h2->Draw();
pad3->cd();
h3->Draw();
pad4->cd();
h4->Draw();
pad5->cd();
h5->Draw();
pad6->cd();
h6->Draw();
return 0;
}