- Timestamp:
- 05/04/12 16:28:06 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
fact/tools/rootmacros/PulseTemplates/pixel.C
r13544 r13559 97 97 if (mVerbosityLevel > 2) cout << endl << "...book pixel histograms" << endl; 98 98 TString histo_name; 99 TString histo_title; 99 100 100 101 for (int order = 0; order < mMaxPulseOrder; order++) … … 104 105 histo_name += "_"; 105 106 histo_name += order; 107 108 histo_title = "Overlay of detected pulses of one pulse order for one Pixel "; 109 histo_title += "[Pixel_Order] "; 110 histo_title += mChid; 111 histo_title += "_"; 112 histo_title += order; 113 106 114 if (mVerbosityLevel > 3) cout << "\t...booking " << histo_name << endl; 107 115 hMaxOverlay[order]=new TH2F( 108 116 histo_name, 109 "Overlay of detected pulses of one pulse order for one Pixel",117 histo_title, 110 118 mPixelOverlayXaxisLeft + mPixelOverlayXaxisRight , 111 119 (-1*mPixelOverlayXaxisLeft)-0.5, … … 145 153 histo_name += "_"; 146 154 histo_name += order; 155 156 histo_title = "Overlay at rising edge of detected pulses of one pulse order for one Pixel "; 157 histo_title += "[Pixel_Order] "; 158 histo_title += mChid; 159 histo_title += "_"; 160 histo_title += order; 161 147 162 if (mVerbosityLevel > 3) cout << "\t...booking " << histo_name << endl; 148 163 hEdgeOverlay[order] = new TH2F( 149 164 histo_name, 150 "Overlay at rising edge of detected pulses of one pulse order for one Pixel ",165 histo_title, 151 166 mPixelOverlayXaxisLeft + mPixelOverlayXaxisRight , 152 167 (-1*mPixelOverlayXaxisLeft)-0.5, … … 178 193 histo_name += "_"; 179 194 histo_name += order; 195 196 histo_title = "Mean value of each slice in overlay plot (Tprofile) "; 197 histo_title += "[Pixel_Order] "; 198 histo_title += mChid; 199 histo_title += "_"; 200 histo_title += order; 201 180 202 if (mVerbosityLevel > 3) cout << "\t...booking " << histo_name << endl; 181 203 hMaxProfile[order] = new TProfile( 182 204 histo_name, 183 "Mean value of each slice in overlay plot (Tprofile)",205 histo_title, 184 206 mPixelOverlayXaxisLeft + mPixelOverlayXaxisRight ,//nbinsx 185 207 (-1*mPixelOverlayXaxisLeft)-0.5, //xlow … … 212 234 histo_name += "_"; 213 235 histo_name += order; 236 237 histo_title = "Mean value of each slice in overlay plot (Tprofile) "; 238 histo_title += "[Pixel_Order] "; 239 histo_title += mChid; 240 histo_title += "_"; 241 histo_title += order; 242 214 243 if (mVerbosityLevel > 3) cout << "\t...booking " << histo_name << endl; 215 244 hEdgeProfile[order] = new TProfile( 216 245 histo_name, 217 "Mean value of each slice in overlay plot (Tprofile)",246 histo_title, 218 247 mPixelOverlayXaxisLeft + mPixelOverlayXaxisRight ,//nbinsx 219 248 (-1*mPixelOverlayXaxisLeft)-0.5, //xlow … … 242 271 void 243 272 Pixel::DrawHistograms( 244 TCanvas* pixelCanvas,273 TCanvas** pixelCanvas, 245 274 int* histoFrameNr 246 275 ) … … 249 278 for (int pulse_order = 0; pulse_order < mMaxPulseOrder; pulse_order++) 250 279 { 251 pixelCanvas[pulse_order] .cd( histoFrameNr[0] );280 pixelCanvas[pulse_order]->cd( histoFrameNr[0] ); 252 281 hMaxOverlay[pulse_order]->Draw("COLZ"); 253 pixelCanvas[pulse_order] .cd( histoFrameNr[2] );254 hMaxProfile[pulse_order]->Draw( "COLZ");282 pixelCanvas[pulse_order]->cd( histoFrameNr[2] ); 283 hMaxProfile[pulse_order]->Draw(); 255 284 hEdgeProfile[pulse_order]->Draw("SAME"); 256 285 257 pixelCanvas[pulse_order] .cd( histoFrameNr[1] );286 pixelCanvas[pulse_order]->cd( histoFrameNr[1] ); 258 287 hEdgeOverlay[pulse_order]->Draw("COLZ"); 259 288 } … … 269 298 { 270 299 TString histo_name; 300 TString histo_title; 271 301 for (int order =0; order < mMaxPulseOrder; order++) 272 302 { … … 281 311 histo_name += "_"; 282 312 histo_name += order; 313 histo_title = "Distribution of rising edge [Pixel_Order]"; 314 histo_title += mChid; 315 histo_title += "_"; 316 histo_title += order; 283 317 if (mVerbosityLevel > 3) cout << "\t...booking " << histo_name << endl; 284 318 hSlopeRisingEdge[order] = new TH1F ( 285 319 histo_name, 286 "Distribution of rising edge",320 histo_title, 287 321 600, 288 322 -10.1, … … 300 334 histo_name += "_"; 301 335 histo_name += order; 336 histo_title = "Deviation of rising edge and maximum [Pixel_Order]"; 337 histo_title += mChid; 338 histo_title += "_"; 339 histo_title += order; 302 340 if (mVerbosityLevel > 3) cout << "\t...booking " << histo_name << endl; 303 341 hRisingEdgeToMax[order] = new TH1F ( 304 342 histo_name, 305 "Deviation of rising edge and maximum",343 histo_title, 306 344 600, 307 345 -10.1, … … 319 357 histo_name += "_"; 320 358 histo_name += order; 359 histo_title = "Deviation of rising edge and maximum [Pixel_Order]"; 360 histo_title += mChid; 361 histo_title += "_"; 362 histo_title += order; 321 363 if (mVerbosityLevel > 3) cout << "\t...booking " << histo_name << endl; 322 364 hPosOfMax[order] = new TH1F ( 323 365 histo_name, 324 "Deviation of rising edge and maximum",366 histo_title, 325 367 600, 326 368 -10.1, … … 339 381 void 340 382 Pixel::DrawDistributionHistograms( 341 TCanvas* pixelCanvas,383 TCanvas** pixelCanvas, 342 384 int* histoFrameNr 343 385 ) … … 346 388 for (int pulse_order = 0; pulse_order < mMaxPulseOrder; pulse_order++) 347 389 { 348 pixelCanvas[pulse_order] .cd( histoFrameNr[0] );390 pixelCanvas[pulse_order]->cd( histoFrameNr[0] ); 349 391 hSlopeRisingEdge[pulse_order]->Draw(); 350 pixelCanvas[pulse_order] .cd( histoFrameNr[1] );392 pixelCanvas[pulse_order]->cd( histoFrameNr[1] ); 351 393 hRisingEdgeToMax[pulse_order]->Draw(); 352 394 353 pixelCanvas[pulse_order] .cd( histoFrameNr[2] );395 pixelCanvas[pulse_order]->cd( histoFrameNr[2] ); 354 396 hPosOfMax[pulse_order]->Draw(); 355 397 }
Note:
See TracChangeset
for help on using the changeset viewer.