SimpleIterate.C

Go to the documentation of this file.

{
//////////////////////////////////////////////////////////
// Simple code for iterating over the FD  MtericsTree
// and printing some DQ metics in the order they 
// appear in the TTree
// Usage: root -b -l -q SimpleIterate.C
//  And in the code below edit the std::vector printConditions
//  to add what you want to print.
//////////////////////////////////////////////////////////


//Reset ROOT and connect tree file
   gROOT->Reset();
   TFile *f = (TFile*)gROOT->GetListOfFiles()->FindObject("MetricsTree_All.root");
   if (!f) {
      f = new TFile("MetricsTree_All.root");
   }
   TTree * mytree;
    f->GetObject("mytree",mytree);

//Declaration of leaves types
   Int_t           run;
   Int_t           subrun;
   Int_t           par;
   Int_t           nevents;
   Int_t           firstsec;
   Int_t           lastsec;
   Int_t           ngoodpix;
   Int_t           ngoodfeb;
   Int_t           ngooddcm;
   Int_t           ngooddb;
   Int_t           ngoodmip;
   Int_t           ngoodcdb;
   Int_t           nactivefeb;
   Int_t           nactivedcm;
   Int_t           nactivedb;
   Int_t           pass_runlen;
   Int_t           pass_empty;
   Int_t           pass_hits;
   Int_t           pass_db;
   Int_t           pass_time;
   Int_t           pass_all;
   Int_t           setsize;
   Double_t        rectimesec;
   Double_t        setlivetime;
   Double_t        hitrate;
   Double_t        midhitrate;
   Double_t        miprate;
   Double_t        midmiprate;
   Double_t        mipratio;
   Double_t        mipasym;
   Double_t        nactivechannels;
   Double_t        emptypercentage;
   Int_t           dbencoded;
   Int_t           dbaencoded;
   Int_t           pass_reco;
   Int_t           pass_slc;
   Int_t           pass_trk;
   Int_t           corrupted;
   Int_t           procsec;
   Double_t        numslc;
   Double_t        trkfrac2D;
   Int_t           preliminary;

   // Set branch addresses.
   mytree->SetBranchAddress("run",&run);
   mytree->SetBranchAddress("subrun",&subrun);
   mytree->SetBranchAddress("par",&par);
   mytree->SetBranchAddress("nevents",&nevents);
   mytree->SetBranchAddress("firstsec",&firstsec);
   mytree->SetBranchAddress("lastsec",&lastsec);
   mytree->SetBranchAddress("ngoodpix",&ngoodpix);
   mytree->SetBranchAddress("ngoodfeb",&ngoodfeb);
   mytree->SetBranchAddress("ngooddcm",&ngooddcm);
   mytree->SetBranchAddress("ngooddb",&ngooddb);
   mytree->SetBranchAddress("ngoodmip",&ngoodmip);
   mytree->SetBranchAddress("ngoodcdb",&ngoodcdb);
   mytree->SetBranchAddress("nactivefeb",&nactivefeb);
   mytree->SetBranchAddress("nactivedcm",&nactivedcm);
   mytree->SetBranchAddress("nactivedb",&nactivedb);
   mytree->SetBranchAddress("pass_runlen",&pass_runlen);
   mytree->SetBranchAddress("pass_empty",&pass_empty);
   mytree->SetBranchAddress("pass_hits",&pass_hits);
   mytree->SetBranchAddress("pass_db",&pass_db);
   mytree->SetBranchAddress("pass_time",&pass_time);
   mytree->SetBranchAddress("pass_all",&pass_all);
   mytree->SetBranchAddress("setsize",&setsize);
   mytree->SetBranchAddress("rectimesec",&rectimesec);
   mytree->SetBranchAddress("setlivetime",&setlivetime);
   mytree->SetBranchAddress("hitrate",&hitrate);
   mytree->SetBranchAddress("midhitrate",&midhitrate);
   mytree->SetBranchAddress("miprate",&miprate);
   mytree->SetBranchAddress("midmiprate",&midmiprate);
   mytree->SetBranchAddress("mipratio",&mipratio);
   mytree->SetBranchAddress("mipasym",&mipasym);
   mytree->SetBranchAddress("nactivechannels",&nactivechannels);
   mytree->SetBranchAddress("emptypercentage",&emptypercentage);
   mytree->SetBranchAddress("dbencoded",&dbencoded);
   mytree->SetBranchAddress("dbaencoded",&dbaencoded);
   mytree->SetBranchAddress("pass_reco",&pass_reco);
   mytree->SetBranchAddress("pass_slc",&pass_slc);
   mytree->SetBranchAddress("pass_trk",&pass_trk);
   mytree->SetBranchAddress("corrupted",&corrupted);
   mytree->SetBranchAddress("procsec",&procsec);
   mytree->SetBranchAddress("numslc",&numslc);
   mytree->SetBranchAddress("trkfrac2D",&trkfrac2D);
   mytree->SetBranchAddress("preliminary",&preliminary);

//     This is the loop skeleton
//       To read only selected branches, Insert statements like:
// mytree->SetBranchStatus("*",0);  // disable all branches
// TTreePlayer->SetBranchStatus("branchname",1);  // activate branchname

   Long64_t nentries = mytree->GetEntries();

   Long64_t nbytes = 0;
   for (Long64_t i=0; i<nentries;i++) {
      nbytes += mytree->GetEntry(i);

      //Define here conditions for printing output
      std::vector<bool> printConditions;
      printConditions.push_back( bool(i%10000==0)               ); //Every 10K rows
      //printConditions.push_back( bool(i>=1405900 && i<=1407000) ); //Range 1
      //printConditions.push_back( bool(i>=1407450 && i<=1409000) ); //Range 2
      printConditions.push_back( bool(i>=1350000 && i<=1360000) ); //Range 2

      //Compose (logical OR) all the conditions 
      bool jointPrintCondition = 0;
      for(int k=0; k<printConditions.size();k++) 
        jointPrintCondition = jointPrintCondition || printConditions[k];
      
      //Actually print
      if (jointPrintCondition){
        std::cout 
          << "Row "
          << std::setw(7) << i
          << " :  "
          << std::setw(5) << run
          << "/" 
          << std::setw(2) << subrun
          << " GoodDB " 
          << std::setw(2) << ngooddb
          << " PassAll " 
          << std::setw(1) << pass_all
          << " Length " 
          << std::setw(4) << lastsec-firstsec
          << " DBdecoded "
          << std::setw(14) << std::bitset<14>(dbencoded)
          << std::endl;
      }

   }

   //Happy ending!
   std::cout<<"Successfully completed the loop over ttree" << std::endl;
}