profilehist.cpp 3.14 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
void profilehist(){
  double row, col;
  TCanvas* c1 = new TCanvas("c1");
  // std::string hist_name = "h" + std::to_string(row) + "_" + std::to_string(col);
  TH2D* histo_bl = new TH2D("histo_bl", "Baseline distribution",256,0.,256.,256,0.,256.);
  TH2D* histo_rms = new TH2D("histo_rms", "Noise distribution",256,0.,256.,256,0.,256.);
  
  for (row = 0; row < 256; row++){
    for(col = 0; col < 256; col++){
      stringstream foldername, filename;

       foldername << std::setw(3) << std::setfill('0') << row;
       // std::cout << " Foldername : " << foldername.str()<< std::endl;

       filename << std::setw(3) << std::setfill('0') << row << "_"
		<< std::setw(3) << std::setfill('0') << col << ".dat";
       //cout << "Filename : " << filename.str() << std::endl;

     
      
      std::ifstream infile(foldername.str()+"/"+filename.str());
      if(infile.is_open()){
	//	std::cout<<"The data file  is opened..."<<std::endl;
	 std::string line;
	 double avr=0.,rms=0.;
	 
	 while(!infile.eof()){
	   std::getline(infile,line);
	   if(line.rfind("#",0)!=0){
	     continue;
	   }
	   //  cout<<"original line : " <<line<<std::endl;
	   std::string newline = " ";
	   for(int i=0;i<line.length();i++){
	     if(line[i]!=':'){
	       newline=newline+line[i];
	     }
	     else{
	       newline=newline+line[i]+" ";
	     }
	   }

	   // cout<<"The new line with spaces " <<newline<<std::endl;
      
	   std::stringstream s;
	   s<<newline;
	   //cout<<"The line is stored in the stringstream object..."<<std::endl;
	   string temp;
	   double num=0.;
	   double arr[2]={ };
	   int i=0;
	   while(!s.eof()){
	     s>>temp;
	     if(std::stringstream(temp)>>num){
	       //cout<<num<<std::endl;
	       arr[i] =num;
	       i++;
	     }
	     temp = " ";
	   }
	   avr = arr[0];
	   rms = arr[1];
	 }

	 // cout << " avr : " << avr << " " << " rms : " << rms << std::endl;

	 if(avr >= 0 && rms >=0){
	   histo_bl->Fill(row,col,avr);
	   histo_rms->Fill(row,col,rms);
	 }
	 /*	 else{
	   cout<< " Foldername : " << foldername.str() << " " <<  " filename : " << filename.str() << " "
	       <<  " avr : " << avr << " " << " rms : " << rms <<  std::endl;
	       }*/
       
      }
      else{
	cout<<"Error..."<<std::endl;
      }
    }
  }

  c1->Divide(2,3);
  c1->cd(1);
  histo_bl->SetStats(0);
  // gStyle->SetOptStat(111100);
  histo_bl->Draw("COLZ");
  
  
  histo_bl->GetXaxis()->SetTitle("Row");
  histo_bl->GetYaxis()->SetTitle("Column");
  histo_bl->GetZaxis()->SetTitle("Baseline");
  histo_bl->GetZaxis()->SetRangeUser(1000,1200);

  c1->cd(2);
   histo_rms->SetStats(0);
  // gStyle->SetOptStat(111100);
  histo_rms->Draw("COLZ");
  
  
  histo_rms->GetXaxis()->SetTitle("Row");
  histo_rms->GetYaxis()->SetTitle("Column");
  histo_rms->GetZaxis()->SetTitle("Noise");
  histo_rms->GetZaxis()->SetRangeUser(5,25);

  c1->cd(3);
  TProfile *px = histo_rms->ProfileX("px",0.,256.);
  px->Draw("HIST");

  c1->cd(4);
  TProfile *py = histo_rms->ProfileY("py",0.,256.);
  py->Draw("HIST");

  c1->cd(5);
  TH1D *projx = histo_rms->ProjectionX("projx", 0., 256.); 
  projx->Draw("HIST");

  c1->cd(6);
  TH1D *projy = histo_rms->ProjectionY("projy", 0., 256.); 
  projy->Draw("HIST");

}