/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/SimMuon/Neutron/src/AsciiNeutronReader.cc

Go to the documentation of this file.

#include "SimMuon/Neutron/src/AsciiNeutronReader.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/EDMException.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include <sstream>
#include<iostream>

using namespace std;

AsciiNeutronReader::AsciiNeutronReader(string fileNameBase) :
  theFileNameBase(fileNameBase),
  //theStreamPos(nChamberTypes+1, 0) //TODO WON"T WORK!  replace with a map 
  theStreamPos(11,0)
{
}


void AsciiNeutronReader::readNextEvent(int chamberType, edm::PSimHitContainer & result) {
  stringstream fileName;
  fileName << theFileNameBase << chamberType;
  ifstream fin(fileName.str().c_str(), ios::in);
  if(!fin.is_open()) {
    throw cms::Exception("NeutronReader") 
      << "Muon neutron noise file missing " << fileName.str();
  }

  int nhits = read_nhits(fin, chamberType);
  for(int ihit = 0; ihit < nhits; ++ihit) {
    float entryX, entryY, entryZ, exitX, exitY, exitZ;
    float p, tof, eloss, theta, phi;
    int type, layer, track;

    fin >> entryX >> entryY >> entryZ
        >> exitX  >> exitY  >> exitZ
        >> p >> tof >> eloss >> type >> layer >> track >> theta >> phi;
    LocalPoint entry(entryX, entryY, entryZ);
    LocalPoint exit (exitX,  exitY,  exitZ);
    PSimHit phit(entry, exit, p, tof, eloss, type, layer, track, theta, phi);
    result.push_back(phit);
  }
  theStreamPos[chamberType] = fin.tellg();
}


int AsciiNeutronReader::read_nhits(ifstream & fin, int chamberType) {
  int nhits = 0;
    // go to the last saved place
  fin.seekg(theStreamPos[chamberType]);
  if(fin.eof()) {
    resetStreampos(fin, chamberType);
  }
  LogDebug("NeutronReader") << "starting from pos " << theStreamPos[chamberType] 
                  << " EOF " << fin.eof();
  fin >> nhits;
  if(fin.eof()) {
    resetStreampos(fin, chamberType);
    fin >> nhits;
  }
  LogDebug("NeutronReader") << "returning nhits " << nhits;
  return nhits;
}


void AsciiNeutronReader::resetStreampos(ifstream & fin, int chamberType) {
  LogDebug("NeutronReader") << "reached EOF, resetting streampos ";
  theStreamPos[chamberType] = 0;
  fin.clear();
  fin.seekg(0, ios::beg);
}