C++ ROOT:当我尝试打开文件并从中读取输入时,ROOT 关闭

C++ ROOT: ROOT closes when I try to open a file and read input from it

我正在使用 CERN ROOT 6.22/00(根据 class 的要求)。我正在尝试读取包含两列数据的“输入文件”,如下所示:

40000 1397251483
40000 1397251484
40010 1397251485
40012 1397251486
40004 1397251487
40003 1397251488
40014 1397251489

这是产生错误的我的代码的最小可重现版本:

# include <iostream> // Header that defines the standard input/output stream objects:
# include <fstream> // Input/output stream class to operate on files.
# include <math.h> // Header declares a set of functions to compute common mathematical operations and transformations.
# include <iomanip> // Header providing parametric manipulators.

using namespace std;
int main()
{

    ifstream inFile;
        cout << "TEST";
    
    int NumOfRows = 1131635;
    char inputFileName[30] = "input.dat"; //File with the model parameters
    char outputFileName[30] = "output.dat";
    
    const int nArray = NumOfRows + 1;
    
    double paramOne[nArray];
    double T[nArray];
    
    //Reading input parameters from file into arrays//
    
    inFile.open(inputFileName,ios::in);
    
    return 0;

}

然而,每当我继续 运行 这段代码时,ROOT 退出,我又回到了终端。我也尝试过 运行 使用 g++ 的代码,但我收到错误:

Segmentation fault (core dumped)

有什么建议吗?

编辑:所以我继续将数组转换为向量,如下所示:

// V2: converted arrays to vectors to avoid memory problems

# include <iostream> // Header that defines the standard input/output stream objects:
# include <fstream> // Input/output stream class to operate on files.
# include <math.h> // Header declares a set of functions to compute common mathematical operations and transformations.
# include <iomanip> // Header providing parametric manipulators.

using namespace std;
int main()
{

    ifstream inFile;
    
    int NumOfRows = 10;
    char inputFileName[30] = "input.dat"; //File with the model parameters
    char outputFileName[30] = "output.dat";
    
    vector<int> TDC;
    vector<int> T;
    
    //Reading input parameters from file into arrays//
    
    inFile.open(inputFileName,ios::in);
    

    // Warning if file cant be opened
    if(!inFile.is_open()){ 
        cout << "Error opening file. \n";
        //cout << "Giving Retry... \n";
    }
    if(inFile.is_open()){
        cout<<"Input File was opened successfully"<<endl;
    }
    if(inFile.good()){
        cout<<"Input File is ready for reading"<<endl;
    }
    
    cout<< fixed;

    int rejects = 0;
    
    
    //reading file
    if(inFile.is_open()){
    
        // Putting cursor at start of file
        inFile.clear();
        
        //Reading first line
        inFile >> TDC[0] >> T[0];
        cout<<"TDC time"<<setw(15)<<"timestamp"<<endl;
        cout<<TDC[0]<<setw(20)<<T[0]<<endl;
}
    return 0;

}

我目前仍然面临内存问题:

Input File was opened successfully
Input File is ready for reading

 *** Break *** segmentation violation



===========================================================
There was a crash.
This is the entire stack trace of all threads:
===========================================================
#0  0x00007fdc11e1d6e7 in waitpid () from /lib/x86_64-linux-gnu/libc.so.6
#1  0x00007fdc11d88107 in ?? () from /lib/x86_64-linux-gnu/libc.so.6
#2  0x00007fdc129bfed3 in TUnixSystem::StackTrace() () from /home/nick/root/lib/libCore.so.6.22
#3  0x00007fdc129c29c5 in TUnixSystem::DispatchSignals(ESignals) () from /home/nick/root/lib/libCore.so.6.22
#4  <signal handler called>
#5  0x00007fdc1243d8c8 in std::istream::operator>>(int&) () from /usr/lib/x86_64-linux-gnu/libstdc++.so.6
#6  0x00007fdc1320eab8 in ?? ()
#7  0x0000558f085c6f00 in ?? ()
#8  0x0ab62774202a8500 in ?? ()
#9  0x0000000000000000 in ?? ()
===========================================================


The lines below might hint at the cause of the crash.
You may get help by asking at the ROOT forum http://root.cern.ch/forum
Only if you are really convinced it is a bug in ROOT then please submit a
report at http://root.cern.ch/bugs Please post the ENTIRE stack trace
from above as an attachment in addition to anything else
that might help us fixing this issue.
===========================================================
#5  0x00007fdc1243d8c8 in std::istream::operator>>(int&) () from /usr/lib/x86_64-linux-gnu/libstdc++.so.6
#6  0x00007fdc1320eab8 in ?? ()
#7  0x0000558f085c6f00 in ?? ()
#8  0x0ab62774202a8500 in ?? ()
#9  0x0000000000000000 in ?? ()
===========================================================

编辑 2:使用建议的解决方案:

// V2: converted arrays to vectors to avoid memory problems

# include <iostream> // Header that defines the standard input/output stream objects:
# include <fstream> // Input/output stream class to operate on files.
# include <math.h> // Header declares a set of functions to compute common mathematical operations and transformations.
# include <iomanip> // Header providing parametric manipulators.

using namespace std;
int main()
{

    ifstream inFile;
    
    int NumOfRows = 1131636;
    char inputFileName[30] = "input.dat"; //File with the model parameters
    char outputFileName[30] = "output.dat";
    
    size_t reasonableTDCSize = 1131635;
    
    vector<int> TDC(NumOfRows);
    vector<int> T(NumOfRows);

    
    //Reading input parameters from file into arrays//
    
    inFile.open(inputFileName,ios::in);
    

    // Warning if file cant be opened
    if(!inFile.is_open()){ 
        cout << "Error opening file. \n";
        //cout << "Giving Retry... \n";
    }
    if(inFile.is_open()){
        cout<<"Input File was opened successfully"<<endl;
    }
    if(inFile.good()){
        cout<<"Input File is ready for reading"<<endl;
    }
    
    cout<< fixed;

    int rejects = 0;
    int tempTDC = 0;
    int tempT = 0;
    
    //reading file
    if(inFile.is_open()){
    
        // Putting cursor at start of file
        inFile.clear();
        
        //Reading first line
        inFile >> tempTDC >> tempT;
        TDC.push_back(tempTDC);
        T.push_back(tempT);
        cout<<"TDC time"<<setw(15)<<"timestamp"<<endl;
        cout<<TDC[0]<<setw(20)<<T[0]<<endl;
        
        
        for (int a = 1; a < NumOfRows; a++){
            inFile >> tempTDC >> tempT;
            if ( tempTDC >= 40000 )
            {
               ++rejects;
               break;
            }

            cout<<tempTDC<<setw(20)<<tempT<<endl;
            // Reading rest of file
            TDC.push_back(tempTDC);
            T.push_back(tempT);
            //cout<<Mod[a]<<setw(15)<<z[a]<<setw(15)<<x[a]<<setw(15)<<M[a]<<setw(15)<<L[a]<<setw(15)<<T[a]<<endl;
                
        }
        
        //To show last and first index only, have this line uncommmented and the cout line in loop commented
        cout<<TDC[NumOfRows-1]<<setw(20)<<T[NumOfRows-1]<<endl;
        
        // Close the file.
        inFile.close(); 
    }
    
    /*
    cout<< "Lines remaining " << NumOfRows - rejects << endl;
            
    if(!inFile.is_open()){
        cout<<"Input File closed successfully"<<endl;
    }
    
    cout<< "Timestamp difference between first and last line is: " << T[NumOfRows-1] - T[0] << endl;
    
    
    cout<<"Creating output file"<<endl;
    ofstream outFile(outputFileName);
    outFile<<"TDC time"<<setw(15)<<"timestamp"<<endl; //Header
        for (int a = 1; a < NumOfRows; a++){
            
            // Reading rest of file
            outFile << TDC[a] << T[a];

                
        }
    outFile<<""<<endl;
    
        // Warning if file cant be opened
    if(!outFile.is_open()){ 
        cout << "Error opening file. \n";
        //cout << "Giving Retry... \n";
    }
    if(outFile.is_open()){
        cout<<"Output File was opened successfully"<<endl;
    }
    if(outFile.good()){
        cout<<"Output File is ready for reading"<<endl;
    }
    outFile.close();
    
    if(!outFile.is_open()){
        cout<<"Output File closed successfully"<<endl;
    }
    
    */
    return 0;

}

在改进版本(“V2”)中,您需要适当调整向量的大小。

您可以这样做:

size_t reasonableTDCSize = /* some value */;
size_t reasonableTSize = /* some value */;

vector<int> TDC(reasonableTDCSize);
vector<int> T(reasonableTSize);

我不太明白你想做什么,但是你可以,例如,做:

vector<int> TDC(NumOfRows);
vector<int> T(NumOfRows);

... if 这在语义上是有道理的。 一旦您正确调整了它们的大小,您的代码就可以工作了。