#include<iostream>
#include<vector>
#include<string>
#include "TH2D.h"


using namespace std;

void ThresholdPlotter()
{
    cout<<"start"<<endl;
    string fileName="noise_Mamba.dat";
    
    ifstream l_file(fileName.c_str());
    int channelsNumber=512;
    std::vector<int> thresholVector(channelsNumber,0);

    for(int channel=0;channel<channelsNumber;channel++)
    {
        double l_noiseFromFile=0;
        l_file >> l_noiseFromFile;
        thresholVector[channel]=l_noiseFromFile;
        cout<<"NoiseRetreiver===> channel: "<< channel <<"noise: "<<l_noiseFromFile<<endl;
    }
   
 
    cout<<"create histograms" <<endl;
    TH2D* noiseHistogram= new TH2D("Noise Histogram","Noise (1#sigma)",100,-0.5,channelsNumber+0.5 , 100, 0, 200);
    TH2D*  lowThresholdHistogram=new  TH2D("low Thresholds Histogram","low Thresholds ",100,-0.5,channelsNumber+0.5 , 100, 0, 200);
    TH2D* highThresholdHistogram=new  TH2D("high Thresholds Histogram","high Thresholds ",100,-0.5,channelsNumber+0.5 , 100, 0, 200);

    int lawThresholdMultiplicity=2.5;
    int highThresholdMultiplicity=3;

    cout<<"Fill histograms" <<endl;

    for(int channel=0;channel<channelsNumber;channel++)
    {
        noiseHistogram->Fill(channel, thresholVector[channel]);
        lowThresholdHistogram->Fill(channel,lawThresholdMultiplicity* thresholVector[channel]);
        highThresholdHistogram->Fill(channel,highThresholdMultiplicity* thresholVector[channel]);
    }

    TCanvas * c1 = new TCanvas("c", "c", 600, 800);
    
    noiseHistogram->GetYaxis()->SetTitle("[ADC]");
    noiseHistogram->GetXaxis()->SetTitle("channel");
    noiseHistogram->SetMarkerStyle(2);
    noiseHistogram->SetLineColor(kWhite);

    lowThresholdHistogram->SetMarkerStyle(2);
    lowThresholdHistogram->SetMarkerColor(kGreen);
    lowThresholdHistogram->SetLineColor(kWhite);

    highThresholdHistogram->SetMarkerStyle(2);
    highThresholdHistogram->SetMarkerColor(kRed);
    highThresholdHistogram->SetLineColor(kWhite);

   
    noiseHistogram->Draw();
    lowThresholdHistogram->Draw("same");
    highThresholdHistogram->Draw("same");
    
    c1->BuildLegend(0.7296238,0.8274793,0.9004702,0.8997934);
    gStyle->SetOptStat(0);

}