Subversion Repositories f9daq

Rev

Rev 40 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 40 Rev 50
Line 97... Line 97...
97 
  if (strstr(plopt, "bias") != NULL) {
 97 
  if (strstr(plopt, "bias") != NULL) {
98 
    TCanvas* canvas01 = new TCanvas("canvas01","",800, 800);
 98 
    TCanvas* canvas01 = new TCanvas("canvas01","",800, 800);
99 
    canvas01->cd();
 99 
    canvas01->cd();
100 
    TH2F* h_bias = (TH2F*) rootfile->Get("h_bias");
 100 
    TH2F* h_bias = (TH2F*) rootfile->Get("h_bias");
101 
    h_bias->Draw("colz");
 101 
    h_bias->Draw("colz");
- 
 
 102 
 
- 
 
 103 
    TCanvas* canvas02 = new TCanvas("canvas02","",800, 800);
- 
 
 104 
    //canvas02->Divide(2,2);
- 
 
 105 
    canvas02->cd();
- 
 
 106 
    TH3F* h_biasThreshold = (TH3F*) rootfile->Get("h_biasThreshold");
- 
 
 107 
    //h_biasThreshold->Draw();
- 
 
 108 
    //canvas02->cd(2);
- 
 
 109 
    //h_biasThreshold->GetXaxis()->SetRange(1,5);
- 
 
 110 
    h_biasThreshold->Project3D("zy");
- 
 
 111 
    h_biasThreshold_zy->DrawCopy("colz");
- 
 
 112 
 
- 
 
 113 
    TCanvas* canvas03 = new TCanvas("canvas03","",1200, 1200);
- 
 
 114 
    canvas03->Divide(4,4);
- 
 
 115 
    int x1=0;
- 
 
 116 
    int x2=3;
- 
 
 117 
    int y1=0;
- 
 
 118 
    int y2=3;
- 
 
 119 
    for(int j=x1; j<=x2; j++){
- 
 
 120 
        for(int k=y1; k<=y2; k++){
- 
 
 121 
                int channel = map[j][k];
- 
 
 122 
                int canvasPosition = 4*(4-k-1) + 4 - j;
- 
 
 123 
                canvas03->cd(canvasPosition);
- 
 
 124 
                h_biasThreshold->GetXaxis()->SetRange(channel+1,channel+1);
- 
 
 125 
                char histoName[128];
- 
 
 126 
                sprintf(histoName, "Ch %d", channel);
- 
 
 127 
                h_biasThreshold->Project3D("zy");
- 
 
 128 
                //sprintf(option, "h_biasThresholsd_zy%f", channel);
- 
 
 129 
                h_biasThreshold_zy->SetTitle(histoName);
- 
 
 130 
                h_biasThreshold_zy->DrawCopy("colz");
- 
 
 131 
 
- 
 
 132 
        }
- 
 
 133 
    }
- 
 
 134 
    TCanvas* canvas04 = new TCanvas("canvas04","",1200, 1200);
- 
 
 135 
    canvas04->Divide(4,4);
- 
 
 136 
    int x1=0;
- 
 
 137 
        int x2=3;
- 
 
 138 
        int y1=4;
- 
 
 139 
        int y2=7;
- 
 
 140 
        for(int j=x1; j<=x2; j++){
- 
 
 141 
                for(int k=y1; k<=y2; k++){
- 
 
 142 
                        int channel = map[j][k];
- 
 
 143 
                        int canvasPosition = 4*(8-k-1) + 4 - j;
- 
 
 144 
                        canvas04->cd(canvasPosition);
- 
 
 145 
                        h_biasThreshold->GetXaxis()->SetRange(channel+1,channel+1);
- 
 
 146 
                        char histoName[128];
- 
 
 147 
                        sprintf(histoName, "Ch %d", channel);
- 
 
 148 
                        h_biasThreshold->Project3D("zy");
- 
 
 149 
                        //sprintf(option, "h_biasThresholsd_zy%f", channel);
- 
 
 150 
                        h_biasThreshold_zy->SetTitle(histoName);
- 
 
 151 
                        h_biasThreshold_zy->DrawCopy("colz");
- 
 
 152 
 
- 
 
 153 
                }
- 
 
 154 
        }
- 
 
 155 
        TCanvas* canvas05 = new TCanvas("canvas05","",1200, 1200);
- 
 
 156 
        canvas05->Divide(4,4);
- 
 
 157 
        int x1=4;
- 
 
 158 
        int x2=7;
- 
 
 159 
        int y1=0;
- 
 
 160 
        int y2=3;
- 
 
 161 
        for(int j=x1; j<=x2; j++){
- 
 
 162 
                for(int k=y1; k<=y2; k++){
- 
 
 163 
                        int channel = map[j][k];
- 
 
 164 
                        int canvasPosition = 4*(4-k-1) + 8 - j;
- 
 
 165 
                        canvas05->cd(canvasPosition);
- 
 
 166 
                        h_biasThreshold->GetXaxis()->SetRange(channel+1,channel+1);
- 
 
 167 
                        char histoName[128];
- 
 
 168 
                        sprintf(histoName, "Ch %d", channel);
- 
 
 169 
                        h_biasThreshold->Project3D("zy");
- 
 
 170 
                        //sprintf(option, "h_biasThresholsd_zy%f", channel);
- 
 
 171 
                        h_biasThreshold_zy->SetTitle(histoName);
- 
 
 172 
                        h_biasThreshold_zy->DrawCopy("colz");
- 
 
 173 
 
- 
 
 174 
                }
- 
 
 175 
        }
- 
 
 176 
 
- 
 
 177 
        TCanvas* canvas06 = new TCanvas("canvas06","",1200, 1200);
- 
 
 178 
        canvas06->Divide(4,4);
- 
 
 179 
        int x1=4;
- 
 
 180 
        int x2=7;
- 
 
 181 
        int y1=4;
- 
 
 182 
        int y2=7;
- 
 
 183 
        for(int j=x1; j<=x2; j++){
- 
 
 184 
                for(int k=y1; k<=y2; k++){
- 
 
 185 
                        int channel = map[j][k];
- 
 
 186 
                        int canvasPosition = 4*(8-k-1) + 8 - j;
- 
 
 187 
                        canvas06->cd(canvasPosition);
- 
 
 188 
                        h_biasThreshold->GetXaxis()->SetRange(channel+1,channel+1);
- 
 
 189 
                        char histoName[128];
- 
 
 190 
                        sprintf(histoName, "Ch %d", channel);
- 
 
 191 
                        h_biasThreshold->Project3D("zy");
- 
 
 192 
                        //sprintf(option, "h_biasThresholsd_zy%f", channel);
- 
 
 193 
                        h_biasThreshold_zy->SetTitle(histoName);
- 
 
 194 
                        h_biasThreshold_zy->DrawCopy("colz");
- 
 
 195 
 
- 
 
 196 
                }
- 
 
 197 
        }
- 
 
 198 
 
102 
  }
 199 
  }
103 
 
 200 
 
104 
  if (strstr(plopt, "th") != NULL) {
 201 
  if (strstr(plopt, "th") != NULL) {
105 
   TCanvas* canvas1 = new TCanvas("canvas1","",800,800);
 202 
   TCanvas* canvas1 = new TCanvas("canvas1","",800,800);
106 
   canvas1->cd();
 203 
   canvas1->cd();