Fill data in outputTree.Root

ROOT
, , ,
1

Hi everybody,
I want to get information from a branch to draw a graph, specifically SignalPlane, I want to get data from the 3 column (Top plane) which is the y axis and x axis which is column 2 (Bot Palne). What should I do to get the information below which is my code

image
image1920×1080 290 KB

#include "TH2F.h"
#include "TFile.h"
#include "TCanvas.h"
#include <vector>

void readplane()
{   
    
    TChain *tIn = new TChain("TofHitsTree");
        {
            tIn->Add(TString::Format("/home/thaitho/tof/output/run16573_outputTree.root"));
        }
    // Khai bao
    Double_t EventTimeOfFlight;
    Int_t EventNSignals, EventNHits;
    vector<int> *SignalType, *HitFeb, *HitSampicChannel, *HitSampic, *HitDaqChannel, *HitPlane, *HitBar, *HitEdge;
    vector<double> *SignalPlane, * SignalBar, *SignalPosition, *SignalTime, *HitCell0Time, *HitTotValue, *HitBaseline, *HitRawPeak, *HitPeak, *HitPeakTime, *HitCfTime_0p1, *HitCfTime_0p9, *HitVoltageIntegral, *HitRawTotValue, *HitUnixTime, *HitSampleLength, *HitPeakSample, *HitChannelOnPlane, *HitSignalIndex, *HitIsSaturated, *HitFitSuccess, *HitCfSuccess;
    vector<vector<double>> *SignalHitsListIndex, *HitWaveform, *HitFitParameter;

    tIn->SetBranchAddress("TofHitsTree.EventTimeOfFlight", &EventTimeOfFlight);
    tIn->SetBranchAddress("TofHitsTree.EventNSignals", &EventNSignals);
    tIn->SetBranchAddress("TofHitsTree.EventNHits", &EventNHits);
    tIn->SetBranchAddress("TofHitsTree.SignalType", &SignalType);
    tIn->SetBranchAddress("TofHitsTree.SignalPlane", &SignalPlane);
    tIn->SetBranchAddress("TofHitsTree.SignalBar", &SignalBar);
    tIn->SetBranchAddress("TofHitsTree.SignalPosition", &SignalPosition);
    tIn->SetBranchAddress("TofHitsTree.SignalTime", &SignalTime);
    tIn->SetBranchAddress("TofHitsTree.SignalHitsListIndex", &SignalHitsListIndex);
    tIn->SetBranchAddress("TofHitsTree.HitFeb", &HitFeb);
    tIn->SetBranchAddress("TofHitsTree.HitSampicChannel", &HitSampicChannel);
    tIn->SetBranchAddress("TofHitsTree.HitCell0Time", &HitCell0Time);
    tIn->SetBranchAddress("TofHitsTree.HitTotValue", &HitTotValue);
    tIn->SetBranchAddress("TofHitsTree.HitWaveform", &HitWaveform);
    tIn->SetBranchAddress("TofHitsTree.HitSampic", &HitSampic);
    tIn->SetBranchAddress("TofHitsTree.HitDaqChannel", &HitDaqChannel);
    tIn->SetBranchAddress("TofHitsTree.HitPlane", &HitPlane);
    tIn->SetBranchAddress("TofHitsTree.HitBar", &HitBar);
    tIn->SetBranchAddress("TofHitsTree.HitEdge", &HitEdge);
    tIn->SetBranchAddress("TofHitsTree.HitBaseline", &HitBaseline);
    tIn->SetBranchAddress("TofHitsTree.HitRawPeak", &HitRawPeak);
    tIn->SetBranchAddress("TofHitsTree.HitPeak", &HitPeak);
    tIn->SetBranchAddress("TofHitsTree.HitPeakSample", &HitPeakSample);
    tIn->SetBranchAddress("TofHitsTree.HitPeakTime", &HitPeakTime);
    tIn->SetBranchAddress("TofHitsTree.HitCfTime_0p1", &HitCfTime_0p1);
    tIn->SetBranchAddress("TofHitsTree.HitCfTime_0p9", &HitCfTime_0p9);
    tIn->SetBranchAddress("TofHitsTree.HitVoltageIntegral", &HitVoltageIntegral);
    tIn->SetBranchAddress("TofHitsTree.HitIsSaturated", &HitIsSaturated);
    tIn->SetBranchAddress("TofHitsTree.HitChannelOnPlane", &HitChannelOnPlane);
    tIn->SetBranchAddress("TofHitsTree.HitRawTotValue", &HitRawTotValue);
    tIn->SetBranchAddress("TofHitsTree.HitUnixTime", &HitUnixTime);
    tIn->SetBranchAddress("TofHitsTree.HitSampleLength", &HitSampleLength);
    tIn->SetBranchAddress("TofHitsTree.HitFitSuccess", &HitFitSuccess);
    tIn->SetBranchAddress("TofHitsTree.HitCfSuccess", &HitCfSuccess);
    tIn->SetBranchAddress("TofHitsTree.HitFitParameter", &HitFitParameter);
    tIn->SetBranchAddress("TofHitsTree.HitSignalIndex", &HitSignalIndex);
  
    //  tao Histograms 2 chieu
    TH2F *hist_top_bot = new TH2F("hist_top_bot", "Top vs Bottom Plane", 100, -50, 50, 100, -50, 50);
    TH2F *hist_top_upstream = new TH2F("hist_top_upstream", "Top vs Upstream", 100, -50, 50, 100, -50, 50);
    TH2F *hist_top_downstream = new TH2F("hist_top_downstream", "Top vs Downstream", 100, -50, 50, 100, -50, 50);


    Long64_t nEntries = tIn->GetEntries();
    for (Long64_t iEntry = 2; iEntry < nEntries; ++iEntry) {
        tIn->GetEntry(iEntry);

        hist_top_bot->Fill(1.0,1.0);
        hist_top_upstream->Fill(1.0, 1.0);
        hist_top_downstream->Fill(1.0,1.0);

    }

    // Draw histograms
    TCanvas *top_vs_planes = new TCanvas("top_vs_plane", "Correlations", 800, 600);

    top_vs_planes->Divide(2, 2);
    top_vs_planes->cd(1);
    hist_top_bot->Draw("colz");
    top_vs_planes->cd(2);
    hist_top_upstream->Draw("colz");
    top_vs_planes->cd(3);
    hist_top_downstream->Draw("colz");
    top_vs_planes->Draw();

    // save to file
    TFile *tfile = new TFile("mychart.root","RECREATE");
    hist_top_upstream->Write();
    hist_top_downstream->Write();
    hist_top_bot->Write();
    tfile->Close();
    delete hist_top_upstream;
    delete tIn;
    delete top_vs_planes;
}
2

Hi,

Let me propose a solution which uses RDataFrame:

ROOT::RDataFrame d("TofHitsTree", "run16573_outputTree.root");
auto myGraph1 = d.Graph("xColumn", "yColumn");
myGraph1->Draw("AP")

Does that work for you?

Best,
D

3

Thank Mr. Danilo, i did it

4

This topic was automatically closed 14 days after the last reply. New replies are no longer allowed.