groot/macros/groot100.C

0001 void groot100(int JobNumber, int RunNumber, int SkipEvents, int PRUN_EVTS)
0002 {
0003   cout << "Prepare to skip " << SkipEvents << " Events" << endl;
0004   gSystem->Load("libgroot");
0005
0006   char filename[500];
0007   sprintf(filename, "data/beam_000000%d-0000.evt", RunNumber);
0008   cout << "Analyzing file named: " << filename << endl;
0009   poncsopen(filename);
0010   prun(2);
0011
0012   Quiver::BlobTimeCut      = false;    // default is false
0013   Quiver::PadBlobs         = true;     // default is true
0014   ABlob::RecalibrateCharge = true;     // default is false
0015
0016   Quiver::PaddingLimit = 4;     //default is 1
0017   ABlob::GaussPosition = true;  //default is false
0018
0019   cout << "Prepare to skip " << SkipEvents << " Events" << endl;
0020   setFirstEventNr(SkipEvents);
0021   prun(PRUN_EVTS);
0022   //prun(SkipEvents+10);
0023
0024   groot* Tree = groot::instance();
0025
0026   char outfile[500];
0027   sprintf(outfile,"scratch/Analysis_%010d_%04d.root", RunNumber, JobNumber);
0028   cout << "Outputting Histograms to: " << outfile << endl;
0029   Tree->SaveTheHistograms(outfile);
0030 }