CalibMuon/CSCCalibration/test/readGains.cpp - master - githubmirror/cmssw

CalibMuon/CSCCalibration/test/readGains.cpp

45 строк1 KB

Giulio Eulisse

Snapshot of CMSSW_6_2_0_pre8.

26 июн 2013, 18:12

26 июн 2013, 18:12214ab73

#include <iostream>
#include <fstream>
#include <vector>
#include <iostream>
#include <string>

using namespace std;

int main(){
  float gainSlope, gainIntercept,chi2;
  int index,flag,flag1;
  int nrlines=0;

  std::vector<int>   index_id;
  std::vector<float> Slope;
  std::vector<float> Intercept;
  std::vector<float> gainChi2;
  
  std::ifstream dbdata; 
  dbdata.open("FileName",std::ios::in); 
  if(!dbdata) {
    std::cerr <<"Error: FileName -> no such file!"<< std::endl;
    exit(1);
  }

  while (!dbdata.eof() ) { 
    dbdata >> index >> gainSlope >> gainIntercept >>chi2 >>flag >>flag1; 
    index_id.push_back(index);
    Slope.push_back(gainSlope);
    Intercept.push_back(gainIntercept);
    gainChi2.push_back(chi2);
    nrlines++;
  }
  dbdata.close();
  std::ofstream myGainsFile("GoodVals_FileName",std::ios::out);
 
  for(int i=0; i<nrlines-1;++i){
    if (Slope[i]>6.0 && Slope[i]<11.0){
      myGainsFile<<index_id[i]<<"  "<<Slope[i]<<"  "<<Intercept[i]<<"  "<<gainChi2[i]<<std::endl;
      if (flag==1 || flag1==1){
  std::cout<<"Flag not 0: "<<index_id[i]<<" " <<flag<<"  "<<flag1<<std::endl;
      }
    }
  }
}