void mySelector::SlaveTerminate()
{
// The SlaveTerminate() function is called after all entries or objects
// have been processed. When running with PROOF SlaveTerminate() is called
// on each slave server.
cout<<"-------------End of Slave Job----------------"<<endl;
tNow.Set();tNow.Print();
cout<<"Number of Event: "<<fNumberOfEvents<<", elapsed time: "<<tNow.Convert()-tBegin.Convert()<<"seconds,rate: "<<float(fNumberOfEvents)/(tNow.Convert()-tBegin.Convert())<<endl;
}
void mySelector::Begin(TTree * /*tree*/)
{
// The Begin() function is called at the start of the query.
// When running with PROOF Begin() is only called on the client.
// The tree argument is deprecated (on PROOF 0 is passed).
TString option = GetOption();
tBegin.Set();
fNumberOfEvents = 0;
cout<<"-------Begin of job-------"<<endl;
tBegin.Print();
}
void MySelector::Begin(TTree * /*tree*/)
{
// The Begin() function is called at the start of the query.
// When running with PROOF Begin() is only called on the client.
// The tree argument is deprecated (on PROOF 0 is passed).
TString option = GetOption();
// some time measurement
tBegin.Set(); printf("*==* ---------- Begin of Job ----------");
tBegin.Print();
}
void mySelector::Terminate()
{
// The Terminate() function is the last function to be called during
// a query. It always runs on the client, it can be used to present
// the results graphically or save the results to file.
TString option = GetOption();
TString output = "gr_";
output += option;
output += ".root";
TFile *hfile = new TFile(output,"RECREATE","FONLL CCbar cross section");
hfile->cd();
const Int_t npoint = eventnumber;
Float_t x[npoint];
Float_t y[npoint];
Float_t ylow[npoint];
Float_t yup[npoint];
for(int i=0;i<npoint;i++)
{
x[i] = Pt[i];
y[i] = Central[i];
ylow[i] = Down[i];
yup[i] = Up[i];
}
//TGraph *grFONLLD0 = new TGraph(npoint,Pt,Central);
TDirectoryFile *ratioErr = new TDirectoryFile(option,"ratio error of scale pp500 to pp200 ");
ratioErr->SetName(option);
ratioErr->Add(hRatio);
TGraph *grFONLLRatio = new TGraph(npoint,x,y);
grFONLLRatio->SetName("grFONLLRatio");
TGraph *grFONLLRatio_u = new TGraph(npoint,x,yup);
grFONLLRatio_u->SetName("grFONLLRatio_u");
TGraph *grFONLLRatio_d = new TGraph(npoint,x,ylow);
grFONLLRatio_d->SetName("grFONLLRatio_d");
grFONLLRatio->Print();
///grFONLLRatio->Write();
ratioErr->Add(grFONLLRatio);
grFONLLRatio_u->Print();
//grFONLLRatio_u->Write();
ratioErr->Add(grFONLLRatio_u);
grFONLLRatio_d->Print();
//grFONLLRatio_d->Write();
ratioErr->Add(grFONLLRatio_d);
ratioErr->Write();
hfile->Print();
hfile->Close();
tNow.Set();
cout<<"----------End of job----------"<<endl;
tNow.Print();
}
void MySelector::SlaveTerminate()
{
// The SlaveTerminate() function is called after all entries or
// objects have been processed. When running with PROOF
// SlaveTerminate() is called on each slave server.
// some statistics at end of job
printf("\n *==* ---------- End of Slave Job ---------- ");
tNow.Set(); tNow.Print();
printf(
"Number of Events: %i, elapsed time: %i sec, rate: %g evts/sec\n"
,fNumberOfEvents,
tNow.Convert()-tBegin.Convert(),
float(fNumberOfEvents)/(tNow.Convert()-tBegin.Convert()) );
}
void MySelector::Terminate()
{
// The Terminate() function is the last function to be called
// during a query. It always runs on the client, it can be used
// to present the results graphically or save the results to
// file.
// finally, store all output
TFile hfile("MySelector_Result.root","RECREATE","MuonResults");
fOutput->Write();
//Example to retrieve output from output list
h_resistance=
dynamic_cast<TH1F *>(fOutput->FindObject("resistance"));
TCanvas c_result("cresult","Resistance",100,100,300,300);
h_resistance->Draw();
c_result.SaveAs("ResistanceDistribution.png");
tNow.Set(); printf("*==* ---------- End of Job ---------- ");
tNow.Print();
}
