Textfile lesen

Zafar

Dieser Beitrag wurde gelöscht!

Zafar

@wob Ich habe die ersten 50 Zeilen genommen und es hat problemlos funktioniert. Nur bei mehr als 10000 Zeilen ignoriert er die ersten Zeilen, voll komisch.

Das Output startet einfach ab Zeile 1279

wob

@Zafar sagte in Textfile lesen:

@wob Ich habe die ersten 50 Zeilen genommen und es hat problemlos funktioniert. Nur bei mehr als 10000 Zeilen ignoriert er die ersten Zeilen, voll komisch.

Das Output startet einfach ab Zeile 1279

Du machst etwas anderes falsch oder die Textdateien enthalten nicht das, was du denkst. Die Funktion weiß doch am Anfang noch gar nicht, wie viele Zeilen es geben wird.

Wie sieht dein vollständiges Programm aus, woher weißt du, dass es die ersten Zeilen nicht korrekt einliest? Hinterfrage jedes Detail deiner Annahmen!

(und zuletzt: bist du sicher, dass du C++ brauchst, geht nicht Python? Geht es um parton density functions von Blei-Blei und Muonen-Transversalimpuls? Wäre meine Vermutung anhand des Dateinamens...)

Zafar

@wob Ich habe was ausprobiert, habe in der while Schreife einen counter eingebaut und dann diesen counter ausgeben lassen. Bei beispielsweise 500 Zeilen gibt der counter mir auch 500 an. Vielleicht liegt es daran, dass mein bash diese große Anzahl an Zeilen net komplett ausgeben kann. Ich benutze zur Info Linux.

Ich benutze ein framework, das sich root nennt. das ist auf c++ basiert

wob

Was hat bash jetzt damit zu tun?! (es gibt auch ein paar Commits von mir im ROOT-Source...)

Rat: Nimm Python (ok, zumindest für soche Sache wie Einlesen)! ROOT hat doch wundervolle Python-Bindings.

Zafar

@wob einfach den Code fürs Einlesen in python umwandeln?

wob

Das ist nur ein genereller Rat. Weil solche Dinge in Python sehr viel einfacher sind. Aber du hast doch schon praktisch den ganzen Einlesecode. Ich kann nicht erraten, was bei dir falsch läuft.

Liest du in einen TTree ein? Nutzt du RDataFrames?

Zafar

@wob Weder das eine noch das andere. Ich lese nur die Datei ein mit ifstream und will daraus ein Histogram erstellen

wob

kannst du die Datei irgendwo hochladen?

Zafar

@wob Hast du Dropbox?

wob

nö

Zafar

https://github.com/zmomtaz/Textfile-.git

Hier ist die Textdatei

wob

Also bei mir kommt da raus:

$ g++ orig.cpp&&./a.out|grep Found|wc -l
11276

Aber warum nicht einfach alles wegschmeißen und es so machen: root starten und eingeben:

auto t = new TTree("t", "t");
t->ReadFile("data.txt", "q2lo/D:q2hi:ylo:yhi:qTlo:qThi:PDF0:uncertainty");
t->Draw("PDF0>>h(100,-20,20)");

oder halt nach dem ReadFile einfach nen TBrowser starten.

Zafar

Aus dem TTree kann ich "yield" einlesen oder?

wob

@Zafar sagte in Textfile lesen:

Aus dem TTree kann ich "yield" einlesen oder?

Das, was du offenbar yield nennst, ist in der Datei doch die vorletzte Spalte. Die scheint doch der Überschrift nach PDF0 zu heißen. Aber den Namen kannst du doch frei vergeben in dem ReadFile.

Generell: ich würde dir gleiche Namen für gleiche Dinge raten. Wenn es in der Überschrift q2lo und q2hi heißt, warum nicht dann auch die Variablen exakt so benennen?

DirkB

@Zafar sagte in Textfile lesen:

@wob Ich habe die ersten 50 Zeilen genommen und es hat problemlos funktioniert. Nur bei mehr als 10000 Zeilen ignoriert er die ersten Zeilen, voll komisch.

Das Output startet einfach ab Zeile 1279

Wie groß ist denn der Buffer deiner Konsole?

Leite die Ausgabe doch mal in eine Datei um.

wob

@DirkB sagte in Textfile lesen:

Wie groß ist denn der Buffer deiner Konsole?

Da die Datei 11276 Datenzeilen hat, rate ich einfach mal: 10.000 Zeilen

Zafar

Der Zeilenpuffer Limit war tatsächlich bei 10000. Es lag tatsächlich an diesem Limit. Es war nur eine kleine Feinjustierung nötig

Trotzdem vielen Dank für die Hilfen!

Zafar

@wob , ich habe versucht meinen Histogram mit dem yield und den Bins zu erstelle. ich bekomme irgendwie eine Endlosschleife , wobei der yield immer 0 ist bei printf. Die Variablen nbinsx und nbinsy habe ich vorher deklariert und sind die Anzahl der Bins in x und y Richtung.

while (in) {
        in_pos = in.tellg();
        if (in >> m_low >> m_high >> y_low >> y_high >> pt_low >> pt_high >> yield >>
            uncertainty) {
			if(y_low > 6) continue;
			yield = yield  * 1e-15 * 2.632e-11;
			for (int i=1;i<=nbinsx;i++)
			{
				for(int j=1;j<=nbinsy;j++)
				
				{	printf("Found entry, yield is %f\n", yield);
					drell->SetBinContent(i,j,yield);
				}
       
			}

wob

Generell: ich nehme an, dass drell vom Typ TH2* ist. Was ist nicht verstehe ist, warum du überhaupt eine Loop über die Bins machst und dort SetBinContent aufrufst. Es ergibt doch gar keinen Sinn, innderhalb der Loop jedes Bin auf denselben Wert zu setzen. Warum nutzt du nicht Fill(...)?

Was soll denn auf die X-Achse und was auf die Y-Achse? Dieses schöne Programm hin oder her, aber wäre es nicht viel viel einfacher, aus dem Tree (siehe mein Dreizeiler oben) ein Histogramm zu erzeugen?

Wert = 0: dir ist schon klar, dass 1e-15 * 2.632e-11 = 2.63e-26 ist? Und dass das verdammt klein ist? Was sind das überhaupt für Werte? Nutze Konstanten und benenne diese Werte!