|
|
| Line 1: |
Line 1: |
| − | #include "TROOT.h"
| + | == t == |
| − | #include "TCanvas.h"
| + | |
| − | #include "TEfficiency.h"
| + | |
| − | #include "TF1.h"
| + | |
| − | #include "TFile.h"
| + | |
| − | #include "TH1D.h"
| + | |
| − | #include "TH2D.h"
| + | |
| − | #include "TLatex.h"
| + | |
| − | #include "TList.h"
| + | |
| − | #include "TPad.h"
| + | |
| − | #include "TPaveStats.h"
| + | |
| − | #include "TPaveText.h"
| + | |
| − | #include "TProfile.h"
| + | |
| − | #include "TString.h"
| + | |
| − | #include "TStyle.h"
| + | |
| − | #include "TText.h"
| + | |
| − | #include "TTree.h"
| + | |
| − | #include "TGraph.h"
| + | |
| − | | + | |
| − | #include <cstring>
| + | |
| − | #include <iostream>
| + | |
| − | #include <fstream>
| + | |
| − | #include <stdio.h>
| + | |
| − | #include <stdlib.h>
| + | |
| − | | + | |
| − | using namespace std;
| + | |
| − | | + | |
| − | static const Int_t NUMPMT = 14;
| + | |
| − | static const Int_t NUMPIXEL = 16;
| + | |
| − | static const Int_t NUMPADDLE = 14;
| + | |
| − | static const Int_t NUMPADDLES = NUMPMT*NUMPIXEL;
| + | |
| − | | + | |
| − | static const float adc_charge = 50*1e-15; // 50 fC, corresponding to an ADC chan
| + | |
| − | static const float e = 1.6e-19; // C, electron charge
| + | |
| − | static const Int_t xcanvas = 800; // width of canvases
| + | |
| − | static const Int_t ycanvas = 800; // height of canvases
| + | |
| − | | + | |
| − | TTree *T; // tree from root file, raw adc/tdc/hit data
| + | |
| − | // local arrays to hold data per event from tree
| + | |
| − | Double_t adc[NUMPADDLES];
| + | |
| − | Double_t adc_c[NUMPADDLES];
| + | |
| − | Double_t tdcl[NUMPADDLES];
| + | |
| − | Double_t tdcl_c[NUMPADDLES];
| + | |
| − | Double_t tdct[NUMPADDLES];
| + | |
| − | Double_t tdct_c[NUMPADDLES];
| + | |
| − | Double_t nahit;
| + | |
| − | Double_t nthit;
| + | |
| − | Double_t nhit;
| + | |
| − | | + | |
| − | //Latest map
| + | |
| − | Int_t pixel1[NUMPMT]={4, 4, 2, 1, 4, 3, 4,13, 4,13,13,13, 5,13};
| + | |
| − | Int_t pixel2[NUMPMT]={5, 8,11,13,14,13,16,16,16,16,16,16,16,16};
| + | |
| − | //Int_t pixel1[NUMPMT]={4, 4, 2, 1, 4, 3, 6,13, 4,13,13,13, 5,13};
| + | |
| − | //Int_t pixel2[NUMPMT]={5, 8,11,13,14,13,16,16,16,16,16,16,16,16};
| + | |
| − | //Original map of Ralph/Tommy
| + | |
| − | //Int_t pixel1[NUMPMT]={4, 4, 2, 1, 4, 3, 4,13, 4,13,13,13, 5,13};
| + | |
| − | //Int_t pixel2[NUMPMT]={5, 8,11,13,14,13,16,16,16,16,16,16,16,16};
| + | |
| − | //Map to test for missing pixels
| + | |
| − | //Int_t pixel1[NUMPMT]={0,0,0,0,0,0,0,0,0,0,0,0,0,0};
| + | |
| − | //Int_t pixel2[NUMPMT]={0,0,0,0,0,0,0,0,0,0,0,0,0,0};
| + | |
| − | | + | |
| − | Int_t paddleindex[NUMPADDLES];
| + | |
| − | | + | |
| − | Int_t run; // run number used in titles of plots
| + | |
| − | Int_t n_events_to_analyze; // number of events to analyze ... -1 = all.
| + | |
| − | TStyle *MyStyle = new TStyle("MyStyle","MyStyle");
| + | |
| − | | + | |
| − | void plot(Int_t runno, Int_t n_events=-1){
| + | |
| − |
| + | |
| − | TString filename;
| + | |
| − | filename.Form("scint_%d.root",runno);
| + | |
| − | TFile *_file0 = TFile::Open(filename);
| + | |
| − | | + | |
| − | run=runno;
| + | |
| − | // Set Canvas preferences
| + | |
| − | MyStyle->SetTitleFontSize(0.08);
| + | |
| − | MyStyle->SetTitleX(0.15);
| + | |
| − | MyStyle->SetTitleY(0.99);
| + | |
| − | MyStyle->SetStatW(0.9);
| + | |
| − | MyStyle->SetMarkerStyle(6);
| + | |
| − | gStyle->SetCanvasDefH(xcanvas);
| + | |
| − | gStyle->SetCanvasDefW(ycanvas);
| + | |
| − | gStyle->SetPalette(1);
| + | |
| − | gROOT->SetStyle("MyStyle");
| + | |
| − | | + | |
| − | T = (TTree *)_file0->Get("T");
| + | |
| − | // Setting addresses so we can call them later
| + | |
| − | T->SetBranchAddress("C.cdetm1r.adc",&adc);
| + | |
| − | T->SetBranchAddress("C.cdetm1r.adc_c",&adc_c);
| + | |
| − | T->SetBranchAddress("C.cdetm1r.tdcl",&tdcl);
| + | |
| − | T->SetBranchAddress("C.cdetm1r.tdct",&tdct);
| + | |
| − | T->SetBranchAddress("C.cdetm1r.tdcl_c",&tdcl_c);
| + | |
| − | T->SetBranchAddress("C.cdetm1r.tdct_c",&tdct_c);
| + | |
| − | T->SetBranchAddress("C.cdetm1r.nhit",&nhit);
| + | |
| − | T->SetBranchAddress("C.cdetm1r.nahit",&nahit);
| + | |
| − | T->SetBranchAddress("C.cdetm1r.nthit",&nthit);
| + | |
| − | | + | |
| − | Int_t n_entries = T->GetEntries(); // Checks how many total entries in T
| + | |
| − | cout << "Found " << n_entries << " events"<<endl;
| + | |
| − | if (n_events==-1){
| + | |
| − | n_events_to_analyze = n_entries;
| + | |
| − | cout << "Analyzing all events." << endl;
| + | |
| − | } else {
| + | |
| − | if (n_events < n_entries) {
| + | |
| − | n_events_to_analyze = n_events;
| + | |
| − | cout << "Analyzing " << n_events << " events." << endl;
| + | |
| − | }else{
| + | |
| − | n_events_to_analyze = n_entries;
| + | |
| − | cout << "Analyzing " << n_entries << " events." << endl;
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | }
| + | |
| − | | + | |
| − | | + | |
| − | void setPaddleIndices(){
| + | |
| − | // Setting the location of paddles, taking into account the missing pixels in each PMT
| + | |
| − | for(Int_t pmt=0; pmt<NUMPMT;pmt++){
| + | |
| − | Int_t ipaddle = (pmt+1)*NUMPADDLE+1;
| + | |
| − | for(Int_t pixel=0; pixel<NUMPIXEL;pixel++){
| + | |
| − | Int_t index = pmt*NUMPIXEL+pixel;
| + | |
| − | if (pixel!=pixel1[pmt]-1&&pixel!=pixel2[pmt]-1){
| + | |
| − | ipaddle--;
| + | |
| − | paddleindex[ipaddle]=index;
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | return;
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t getPaddleIndex(Int_t pmt, Int_t ipaddle){
| + | |
| − | //pmt numbers are from 0 to 13 as are ipaddle numbers and the code accounts for removed pixels
| + | |
| − | setPaddleIndices();
| + | |
| − | | + | |
| − | Int_t localpaddle = (pmt+1)*NUMPADDLE-ipaddle;
| + | |
| − | | + | |
| − | return paddleindex[localpaddle];
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t getPaddlePixel(Int_t pmt, Int_t ipaddle){
| + | |
| − | //gives the pixel number for the given paddle number of the pmt, again paddle and pmt numbers are from 0 to 13
| + | |
| − | setPaddleIndices();
| + | |
| − | | + | |
| − | Int_t localpaddle = (pmt+1)*NUMPADDLE-ipaddle;
| + | |
| − | | + | |
| − | return paddleindex[localpaddle]-(pmt)*NUMPIXEL;
| + | |
| − | }
| + | |
| − | | + | |
| − | void print_event(Int_t adc_cut=50, Int_t start_event=1, Int_t num_events=1, Int_t tdc_min=750, Int_t tdc_width=300){
| + | |
| − | | + | |
| − | for (Int_t id=start_event;id<start_event+num_events;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − | cout << "Event " << id << endl;
| + | |
| − | for (Int_t pmt=1; pmt<=NUMPMT;pmt++){
| + | |
| − | for (Int_t index=1; index<=NUMPIXEL; index++){
| + | |
| − | Int_t ipixel = (pmt-1)*NUMPIXEL+index-1;
| + | |
| − | if(tdcl[ipixel]>tdc_min&&tdcl[ipixel]<tdc_min+tdc_width){
| + | |
| − | cout << "Leading Edge TDC Hit on global pixel " << ipixel << " = " << tdcl[ipixel] << endl;
| + | |
| − | }
| + | |
| − | if(tdct[ipixel]>tdc_min&&tdct[ipixel]<tdc_min+tdc_width){
| + | |
| − | cout << "Trailing Edge TDC Hit on global pixel " << ipixel << " = " << tdct[ipixel] << endl;
| + | |
| − | }
| + | |
| − | if(adc_c[ipixel]>adc_cut){
| + | |
| − | cout << "ADC Hit on global pixel " << ipixel << " = " << adc_c[ipixel] << endl;
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − |
| + | |
| − | return;
| + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *plot_adc_fit(Int_t pmt=1, Int_t tdc_min=750, Int_t tdc_width=300, Int_t adc_neighbor_cut=60, Int_t adc_cut=50){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title, grtitle;
| + | |
| − | title.Form("run_%d_ADC_Fit",run);
| + | |
| − | TCanvas *cADCFit = new TCanvas("cADCFit",title,xcanvas,ycanvas);
| + | |
| − | title.Form("run_%d_ADC_Mean_Fit",run);
| + | |
| − | TCanvas *cADCMeanFit= new TCanvas("cADCMeanFit",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | // Setting up different arrays and objects used to create the two canvases.
| + | |
| − | | + | |
| − | TH1D *htmp[NUMPIXEL];
| + | |
| − | TF1 *function;
| + | |
| − | TF1 *function1;
| + | |
| − | Int_t histmax_cut[16], gaus_cut_plus[16];
| + | |
| − | Float_t gaus_cut_minus[16];
| + | |
| − | | + | |
| − | Double_t upixel[NUMPIXEL]={1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16};
| + | |
| − | Double_t epixel[NUMPIXEL]={0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
| + | |
| − | Double_t errors[NUMPIXEL]={0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
| + | |
| − | | + | |
| − | Double_t constants[16], means[16], sigmas[16];
| + | |
| − | | + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.6);
| + | |
| − | MyStyle->SetStatY(0.4);
| + | |
| − | | + | |
| − | setPaddleIndices(); //Setting the geometric paddle locations
| + | |
| − | | + | |
| − | Int_t nbin=90;
| + | |
| − | Int_t min=-100, max=800;
| + | |
| − | | + | |
| − | for(Int_t i = 1; i <= NUMPIXEL; i++)
| + | |
| − | {
| + | |
| − | tmpentry.Form("htmp%d",i);
| + | |
| − | htmp[i-1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | htmp[i-1]->SetLineColor(kRed);
| + | |
| − | title.Form("Run %d ADC Fit pmt %d, paddle %d: %d < tdc < %d",run,pmt,i,tdc_min,tdc_min+tdc_width);
| + | |
| − | htmp[i-1]->SetTitle(title);
| + | |
| − | }
| + | |
| − | | + | |
| − | // Filling the histograms using only adc data with good tdc and good adc neighbor cut.
| + | |
| − | | + | |
| − | Int_t nentries=n_events_to_analyze;
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++)
| + | |
| − | {
| + | |
| − | T->GetEntry(id);
| + | |
| − | Int_t ipaddle = (pmt)*NUMPADDLE+1;
| + | |
| − | for (Int_t pixel=0; pixel < NUMPIXEL; pixel++)
| + | |
| − | {
| + | |
| − | Int_t index = (pmt-1)*NUMPIXEL+pixel;
| + | |
| − | if (pixel!=pixel1[pmt-1]-1 && pixel!=pixel2[pmt-1]-1)
| + | |
| − | {
| + | |
| − | ipaddle--;
| + | |
| − | // if(tdcl[index] > tdc_min && tdcl[index] < tdc_min+tdc_width)
| + | |
| − | // {
| + | |
| − | if (ipaddle < 2)
| + | |
| − | {
| + | |
| − | if (adc_c[paddleindex[ipaddle+1]] < 10)
| + | |
| − | {
| + | |
| − | htmp[pixel]->Fill(adc_c[index]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | else if (ipaddle > NUMPMT*NUMPADDLE) // should this be > or = ?
| + | |
| − | {
| + | |
| − | if (adc_c[paddleindex[ipaddle-1]] < 10)
| + | |
| − | {
| + | |
| − | htmp[pixel]->Fill(adc_c[index]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | else
| + | |
| − | {
| + | |
| − | if (adc_c[paddleindex[ipaddle-1]] < 10 && adc_c[paddleindex[ipaddle+1]] < 10)
| + | |
| − | {
| + | |
| − | htmp[pixel]->Fill(adc_c[index]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | // }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | | + | |
| − | // Creating the canvas of adc data with good tdc and fitting the adc with a gaussian or "landau" function.
| + | |
| − | | + | |
| − | cADCFit->Clear();
| + | |
| − | cADCMeanFit->Clear();
| + | |
| − | cADCFit->Divide(4,4);
| + | |
| − | | + | |
| − | Int_t count = 0;
| + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++)
| + | |
| − | {
| + | |
| − | if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1)
| + | |
| − | {
| + | |
| − | // means[i]=0.0;
| + | |
| − | // int check = 0;
| + | |
| − | cADCFit->cd(count+1);
| + | |
| − | gPad->SetLogy();
| + | |
| − | cADCFit->Update();
| + | |
| − |
| + | |
| − | Int_t entries = htmp[i]->GetEntries();
| + | |
| − | //float mean = htmp[i]->GetMean(1);
| + | |
| − | //float RMS = htmp[i]->GetRMS(1);
| + | |
| − |
| + | |
| − | htmp[i]->SetStats(0);
| + | |
| − | htmp[i]->Draw();
| + | |
| − | | + | |
| − | Double_t bc =0.0;
| + | |
| − | Double_t bn = 0.0;
| + | |
| − | | + | |
| − | for (Int_t j=15; j<100 ;j++){
| + | |
| − | bc = htmp[i]->GetBinContent(j);
| + | |
| − | if( bc == 0.0 || (bc <htmp[i]->GetBinContent(j+1) && bc < htmp[i]->GetBinContent(j+2) && bc <htmp[i]->GetBinContent(j-1) && bc <htmp[i]->GetBinContent(j-2) && bc <htmp[i]->GetBinContent(j-3) && bc <htmp[i]->GetBinContent(j-4)) ){
| + | |
| − | bn = htmp[i]->GetBinCenter(j);
| + | |
| − | htmp[i]->Fit("gaus","","",bn,bn+600);
| + | |
| − | break;
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | //htmp[i]->Fit("landau","","", 0, 250);
| + | |
| − | function = htmp[i]->GetFunction("gaus");
| + | |
| − | function->SetLineColor(1);
| + | |
| − | | + | |
| − | constants[i] = function->GetParameter(0);
| + | |
| − | means[i] = function->GetParameter(1);
| + | |
| − | sigmas[i] = function->GetParameter(2);
| + | |
| − |
| + | |
| − | count++;
| + | |
| − |
| + | |
| − | }
| + | |
| − | }
| + | |
| − |
| + | |
| − | | + | |
| − | // Printing out means and sigmas to screen with corresponding pixel number.
| + | |
| − | | + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++)
| + | |
| − | {
| + | |
| − | cout << "Pixel number: " << i+1 << " \t mean = " << means[i] << " \t sigma = " << sigmas[i] << endl;
| + | |
| − | }
| + | |
| − | | + | |
| − | title.Form("run_%d_ADC_pmt_%d_tdc_min_%d_max_%d.png",run,pmt,tdc_min,tdc_min+tdc_width);
| + | |
| − | cADCFit->Print(title);
| + | |
| − | cADCFit->cd(0);
| + | |
| − | | + | |
| − | // End of work on canvas with adc data fitted with a function
| + | |
| − | | + | |
| − | // Creating the canvas of mean adc data with good tdc data.
| + | |
| − | | + | |
| − | cADCMeanFit->cd();
| + | |
| − | | + | |
| − | Double_t mmean[2]={0,0};
| + | |
| − | Double_t merror[2]={1,1};
| + | |
| − | Double_t mpixel[2]={0,0};
| + | |
| − | Double_t mepixel[2]={0,0};
| + | |
| − | mpixel[0]=pixel1[pmt-1];
| + | |
| − | mpixel[1]=pixel2[pmt-1];
| + | |
| − | | + | |
| − | Double_t yline[2]={100,100};
| + | |
| − | Double_t xline[2]={0,17};
| + | |
| − | | + | |
| − | // Calculating errors on the mean adc data and printing out number of entries, means, and errors
| + | |
| − | // to screen with corresponding pixel number.
| + | |
| − | | + | |
| − | count = 0;
| + | |
| − | for (Int_t i = 0; i < NUMPIXEL; i++)
| + | |
| − | {
| + | |
| − | if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1)
| + | |
| − | {
| + | |
| − | Int_t entries = htmp[i]->GetEntries();
| + | |
| − | errors[i] = sigmas[i]/sqrt(entries);
| + | |
| − | if(means[i] > 1)
| + | |
| − | {
| + | |
| − | cout << "Pixel number: " << i+1 << " \t entries: " << entries << "\t mean: "
| + | |
| − | << means[i] << "\t error: " << errors[i] << endl;
| + | |
| − | }
| + | |
| − | count++;
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | // Creating three different graphs for used pixels, missing pixels, and general formatting.
| + | |
| − | | + | |
| − | gr = new TGraphErrors(NUMPIXEL,upixel,means,epixel,errors);
| + | |
| − | gr->SetMarkerStyle(21);
| + | |
| − | gr->GetXaxis()->SetTitle("Pixel Number");
| + | |
| − | gr->GetYaxis()->SetTitle("Mean ADC Fit (Good TDC)");
| + | |
| − | gr->GetYaxis()->SetTitleOffset(1.4);
| + | |
| − | grtitle.Form("run_%d_pmt_%d_adc_mean",run,pmt);
| + | |
| − | gr->SetTitle(grtitle);
| + | |
| − | gr2 = new TGraphErrors(2,mpixel,mmean,mepixel,merror);
| + | |
| − | gr2->SetMarkerStyle(21);
| + | |
| − | gr2->SetMarkerColor(2);
| + | |
| − | gr2->SetTitle("");
| + | |
| − | gr3 = new TGraph(2,xline,yline);
| + | |
| − | gr3->SetLineColor(2);
| + | |
| − | gr3->SetLineWidth(2);
| + | |
| − | gr3->SetLineStyle(2);
| + | |
| − | gr3->SetTitle("");
| + | |
| − | | + | |
| − | // Drawing the three different graphs to the adc mean values canvas.
| + | |
| − | | + | |
| − | gr->Draw("AP");
| + | |
| − | gr2->Draw("P");
| + | |
| − | gr3->Draw("L");
| + | |
| − |
| + | |
| − | cADCMeanFit->Update();
| + | |
| − | cADCMeanFit->cd(0);
| + | |
| − | | + | |
| − | // End of work on canvas with adc mean data.
| + | |
| − | | + | |
| − | // Creating a file to output parameters of function fitted to the ADC data with good TDC.
| + | |
| − | TString filename;
| + | |
| − | filename.Form("run_%d_pmt_%d_adc_fitresults.txt", run, pmt);
| + | |
| − | ofstream outfile(filename);
| + | |
| − | | + | |
| − | outfile << "Run #: " << run << "\t PMT #: " << pmt << " \t Paramters of fitted ADC data with good TDC. \n";
| + | |
| − | outfile << "Pixel # \t Mean Value \t Error \t Sigma \n" << endl;
| + | |
| − | for(Int_t i = 0; i < NUMPIXEL; i++)
| + | |
| − | {
| + | |
| − | outfile << i+1 << "\t" << means[i] << "\t" << errors[i] << "\t" << sigmas[i] << endl;
| + | |
| − | }
| + | |
| − | cout << "Created output file for parameters." << endl;
| + | |
| − | outfile.close();
| + | |
| − |
| + | |
| − | return cADCMeanFit;
| + | |
| − | return cADCFit;
| + | |
| − | }
| + | |
| − | | + | |
| − |
| + | |
| − | | + | |
| − | TCanvas *plot_adc(Int_t pmt=1, Int_t tdc_min=750, Int_t tdc_width=300){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_ADC",run);
| + | |
| − | TCanvas *cADC= new TCanvas("cADC",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | TH1D *htmpa[NUMPIXEL];//=new TH1D("htmp","htmp",nbin,min,max);
| + | |
| − | TH1D *htmpb[NUMPIXEL];//=new TH1D("htmp1","htmp1",nbin,min,max);
| + | |
| − | | + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.6);
| + | |
| − | MyStyle->SetStatW(0.4);
| + | |
| − | | + | |
| − | Int_t nbin=600;
| + | |
| − | Int_t min=-100, max=500;
| + | |
| − | for(Int_t icounter=1;icounter<=NUMPIXEL;icounter++){
| + | |
| − | tmpentry.Form("htmpa%d", icounter);
| + | |
| − | htmpa[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | tmpentry.Form("htmpb%d", icounter);
| + | |
| − | htmpb[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | htmpa[icounter - 1]->SetLineColor(kBlue);
| + | |
| − | htmpb[icounter - 1]->SetLineColor(kRed);
| + | |
| − | title.Form("Run %d ADC pmt %d, paddle %d: %d < tdc < %d",run,pmt,icounter,tdc_min,tdc_min+tdc_width);
| + | |
| − | htmpa[icounter - 1]->SetTitle(title);
| + | |
| − | htmpb[icounter - 1]->SetTitle(title);
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t nentries=n_events_to_analyze;
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − | for (Int_t index=1; index<=NUMPIXEL; index++){
| + | |
| − | Int_t ipaddle = (pmt-1)*NUMPIXEL+index-1;
| + | |
| − | htmpa[index-1]->Fill(adc_c[ipaddle]);
| + | |
| − | if(tdcl[ipaddle]>tdc_min&&tdcl[ipaddle]<tdc_min+tdc_width){
| + | |
| − | htmpb[index-1]->Fill(adc_c[ipaddle]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | cADC->Clear();
| + | |
| − | cADC->Divide(4,4);
| + | |
| − | | + | |
| − | //plot histos
| + | |
| − | Int_t icount=0;
| + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++){
| + | |
| − | | + | |
| − | if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1) {
| + | |
| − | | + | |
| − | //cout<<"into loop 2, i = " << i << endl;
| + | |
| − | | + | |
| − | cADC->cd( icount + 1 );
| + | |
| − | gPad->SetLogy();
| + | |
| − | | + | |
| − | cADC->Update();
| + | |
| − |
| + | |
| − | Int_t entries = htmpa[i]->GetEntries();
| + | |
| − | float mean = htmpa[i]->GetMean(1);
| + | |
| − | float RMS = htmpa[i]->GetRMS(1);
| + | |
| − | //cout << entries <<" "<< mean <<" "<< RMS <<endl;
| + | |
| − |
| + | |
| − | htmpa[i]->SetStats(0);
| + | |
| − | //current->Modified();
| + | |
| − | htmpa[i]->Draw();
| + | |
| − | htmpb[i]->Draw("same");
| + | |
| − | | + | |
| − | icount++;
| + | |
| − |
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | title.Form("run_%d_ADC_pmt_%d_tdc_min_%d_max_%d.png",
| + | |
| − | run,pmt,tdc_min,tdc_min+tdc_width);
| + | |
| − | cADC->Print(title);
| + | |
| − | cADC->cd(0);
| + | |
| − | return cADC;
| + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *plot_tdc(Int_t pmt=1, Int_t adc_cut=40, Int_t tdc_min=750, Int_t tdc_width=300){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_TDC",run);
| + | |
| − | TCanvas *cTDC= new TCanvas("cTDC",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | TH1D *htmpa[NUMPIXEL];//=new TH1D("htmp","htmp",nbin,min,max);
| + | |
| − | TH1D *htmpb[NUMPIXEL];//=new TH1D("htmp1","htmp1",nbin,min,max);
| + | |
| − | | + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.6);
| + | |
| − | MyStyle->SetStatW(0.4);
| + | |
| − | | + | |
| − | Int_t nbin=tdc_width;
| + | |
| − | Int_t min=tdc_min, max=tdc_min+tdc_width;
| + | |
| − | for(Int_t icounter=1;icounter<=NUMPIXEL;icounter++){
| + | |
| − | tmpentry.Form("htmpa%d", icounter);
| + | |
| − | htmpa[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | tmpentry.Form("htmpb%d", icounter);
| + | |
| − | htmpb[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | htmpa[icounter - 1]->SetLineColor(kBlue);
| + | |
| − | htmpb[icounter - 1]->SetLineColor(kRed);
| + | |
| − | title.Form("Run %d ADC pmt %d, paddle %d: %d < tdc < %d",run,pmt,icounter,tdc_min,tdc_min+tdc_width);
| + | |
| − | htmpa[icounter - 1]->SetTitle(title);
| + | |
| − | htmpb[icounter - 1]->SetTitle(title);
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t nentries = n_events_to_analyze;
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − | for (Int_t index=1; index<=NUMPIXEL; index++){
| + | |
| − | Int_t ipaddle = (pmt-1)*NUMPIXEL+index-1;
| + | |
| − | htmpa[index-1]->Fill(tdcl[ipaddle]);
| + | |
| − | if(adc_c[ipaddle]>adc_cut){
| + | |
| − | htmpb[index-1]->Fill(tdcl[ipaddle]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | cTDC->Clear();
| + | |
| − | cTDC->Divide(4,4);
| + | |
| − | | + | |
| − | //plot histos
| + | |
| − | Int_t icount=0;
| + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++){
| + | |
| − | | + | |
| − | if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1) {
| + | |
| − | | + | |
| − | //cout<<"into loop 2, i = " << i << endl;
| + | |
| − | | + | |
| − | cTDC->cd( icount + 1 );
| + | |
| − | //gPad->SetLogy();
| + | |
| − | | + | |
| − | cTDC->Update();
| + | |
| − |
| + | |
| − | Int_t entries = htmpa[i]->GetEntries();
| + | |
| − | float mean = htmpa[i]->GetMean(1);
| + | |
| − | float RMS = htmpa[i]->GetRMS(1);
| + | |
| − | //cout << entries <<" "<< mean <<" "<< RMS <<endl;
| + | |
| − |
| + | |
| − | htmpa[i]->SetStats(0);
| + | |
| − | //current->Modified();
| + | |
| − | htmpa[i]->Draw();
| + | |
| − | htmpb[i]->Draw("same");
| + | |
| − | | + | |
| − | icount++;
| + | |
| − |
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | title.Form("run_%d_TDC_pmt_%d_adc_cut_%d.png",
| + | |
| − | run,pmt,adc_cut);
| + | |
| − | cTDC->Print(title);
| + | |
| − | cTDC->cd(0);
| + | |
| − | return cTDC;
| + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *plot_tdc_adc(Int_t pmt=1, Int_t tdc_min=1, Int_t tdc_width=100){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_TDC",run);
| + | |
| − | TCanvas *cTDCA= new TCanvas("cTDCA",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | TH2D *htmpa[NUMPIXEL];//=new TH1D("htmp","htmp",nbin,min,max);
| + | |
| − | | + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.6);
| + | |
| − | MyStyle->SetStatW(0.4);
| + | |
| − | | + | |
| − | Int_t nbinx=tdc_width;
| + | |
| − | Int_t minx=tdc_min, maxx=tdc_min+tdc_width;
| + | |
| − | Int_t nbiny=350;
| + | |
| − | Int_t miny=0, maxy=350;
| + | |
| − | for(Int_t icounter=1;icounter<=NUMPIXEL;icounter++){
| + | |
| − | tmpentry.Form("htmpa%d", icounter);
| + | |
| − | htmpa[icounter - 1] = new TH2D(tmpentry,tmpentry,nbinx,minx,maxx,nbiny,miny,maxy);
| + | |
| − | htmpa[icounter - 1]->SetLineColor(kBlue);
| + | |
| − | title.Form("Run %d ADC pmt %d, paddle %d: %d < tdc < %d",run,pmt,icounter,tdc_min,tdc_min+tdc_width);
| + | |
| − | htmpa[icounter - 1]->SetTitle(title);
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t nentries = n_events_to_analyze;
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − | for (Int_t index=1; index<=NUMPIXEL; index++){
| + | |
| − | Int_t ipaddle = (pmt-1)*NUMPIXEL+index-1;
| + | |
| − | htmpa[index-1]->Fill(tdcl[ipaddle]-tdct[ipaddle],adc_c[ipaddle]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | cTDCA->Clear();
| + | |
| − | cTDCA->Divide(4,4);
| + | |
| − | | + | |
| − | //plot histos
| + | |
| − | Int_t icount=0;
| + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++){
| + | |
| − | | + | |
| − | if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1) {
| + | |
| − | | + | |
| − | //cout<<"into loop 2, i = " << i << endl;
| + | |
| − | | + | |
| − | cTDCA->cd( icount + 1 );
| + | |
| − | //gPad->SetLogy();
| + | |
| − | | + | |
| − | cTDCA->Update();
| + | |
| − |
| + | |
| − | Int_t entries = htmpa[i]->GetEntries();
| + | |
| − | float mean = htmpa[i]->GetMean(1);
| + | |
| − | float RMS = htmpa[i]->GetRMS(1);
| + | |
| − | //cout << entries <<" "<< mean <<" "<< RMS <<endl;
| + | |
| − |
| + | |
| − | htmpa[i]->SetStats(0);
| + | |
| − | //current->Modified();
| + | |
| − | htmpa[i]->Draw("COLZ");
| + | |
| − | | + | |
| − | icount++;
| + | |
| − |
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | title.Form("run_%d_TDC_pmt_%d.png",
| + | |
| − | run,pmt);
| + | |
| − | cTDCA->Print(title);
| + | |
| − | cTDCA->cd(0);
| + | |
| − | return cTDCA;
| + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *plot_tdc_diff(Int_t pmt=1, Int_t adc_cut=40, Int_t tdc_min=1, Int_t tdc_width=100){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_TDC",run);
| + | |
| − | TCanvas *cTDCT= new TCanvas("cTDCT",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | TH1D *htmpa[NUMPIXEL];//=new TH1D("htmp","htmp",nbin,min,max);
| + | |
| − | TH1D *htmpb[NUMPIXEL];//=new TH1D("htmp1","htmp1",nbin,min,max);
| + | |
| − | | + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.6);
| + | |
| − | MyStyle->SetStatW(0.4);
| + | |
| − | | + | |
| − | Int_t nbin=tdc_width;
| + | |
| − | Int_t min=tdc_min, max=tdc_min+tdc_width;
| + | |
| − | for(Int_t icounter=1;icounter<=NUMPIXEL;icounter++){
| + | |
| − | tmpentry.Form("htmpa%d", icounter);
| + | |
| − | htmpa[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | tmpentry.Form("htmpb%d", icounter);
| + | |
| − | htmpb[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | htmpa[icounter - 1]->SetLineColor(kBlue);
| + | |
| − | htmpb[icounter - 1]->SetLineColor(kRed);
| + | |
| − | title.Form("Run %d ADC pmt %d, paddle %d: %d < tdc < %d",run,pmt,icounter,tdc_min,tdc_min+tdc_width);
| + | |
| − | htmpa[icounter - 1]->SetTitle(title);
| + | |
| − | htmpb[icounter - 1]->SetTitle(title);
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t nentries = n_events_to_analyze;
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − | for (Int_t index=1; index<=NUMPIXEL; index++){
| + | |
| − | Int_t ipaddle = (pmt-1)*NUMPIXEL+index-1;
| + | |
| − | htmpa[index-1]->Fill(tdcl[ipaddle]-tdct[ipaddle]);
| + | |
| − | if(adc_c[ipaddle]>adc_cut){
| + | |
| − | htmpb[index-1]->Fill(tdcl[ipaddle]-tdct[ipaddle]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | cTDCT->Clear();
| + | |
| − | cTDCT->Divide(4,4);
| + | |
| − | | + | |
| − | //plot histos
| + | |
| − | Int_t icount=0;
| + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++){
| + | |
| − | | + | |
| − | if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1) {
| + | |
| − | | + | |
| − | //cout<<"into loop 2, i = " << i << endl;
| + | |
| − | | + | |
| − | cTDCT->cd( icount + 1 );
| + | |
| − | //gPad->SetLogy();
| + | |
| − | | + | |
| − | cTDCT->Update();
| + | |
| − |
| + | |
| − | Int_t entries = htmpa[i]->GetEntries();
| + | |
| − | float mean = htmpa[i]->GetMean(1);
| + | |
| − | float RMS = htmpa[i]->GetRMS(1);
| + | |
| − | //cout << entries <<" "<< mean <<" "<< RMS <<endl;
| + | |
| − |
| + | |
| − | htmpa[i]->SetStats(0);
| + | |
| − | //current->Modified();
| + | |
| − | htmpa[i]->Draw();
| + | |
| − | htmpb[i]->Draw("same");
| + | |
| − | | + | |
| − | icount++;
| + | |
| − |
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | title.Form("run_%d_TDC_pmt_%d_adc_cut_%d.png",
| + | |
| − | run,pmt,adc_cut);
| + | |
| − | cTDCT->Print(title);
| + | |
| − | cTDCT->cd(0);
| + | |
| − | return cTDCT;
| + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *plot_occupancy(Int_t adc_cut=40, Int_t adc_neighbor_cut=10000, Int_t multiplicity_cut=12, Int_t tdc_min=750, Int_t tdc_width=300){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_Occupancy",run);
| + | |
| − | | + | |
| − | Int_t nbin=196;
| + | |
| − | Int_t min=1, max=197;
| + | |
| − | Int_t nbinm=11;
| + | |
| − | Int_t minm=-1, maxm=10;
| + | |
| − | TH1D *hoccupancy = new TH1D("hoccupancy","hoccupancy",nbin,min,max);
| + | |
| − | TH1D *hmultiplicity = new TH1D("hmultiplicity","hmultiplicity",nbinm,minm,maxm);
| + | |
| − | TH2D *heatmap = new TH2D("heatmap","heatmap",nbin,min,max,nbinm,minm,maxm);
| + | |
| − | | + | |
| − | TCanvas *cOCCUPANCY= new TCanvas("cOCCUPANCY",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | Int_t nentries=n_events_to_analyze;
| + | |
| − |
| + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.9);
| + | |
| − | MyStyle->SetStatW(0.2);
| + | |
| − | | + | |
| − | setPaddleIndices();
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − |
| + | |
| − | Int_t nmultiplicity=0;
| + | |
| − | Int_t good_paddle[500];
| + | |
| − | for(Int_t icount=0;icount<500;icount++){good_paddle[icount]=-1;}
| + | |
| − | | + | |
| − | for(Int_t pmt=0; pmt<NUMPMT;pmt++){
| + | |
| − | Int_t ipaddle = (pmt+1)*NUMPADDLE+1;
| + | |
| − | for(Int_t pixel=0; pixel<NUMPIXEL;pixel++){
| + | |
| − | Int_t index = pmt*NUMPIXEL+pixel;
| + | |
| − | if (pixel!=pixel1[pmt]-1&&pixel!=pixel2[pmt]-1){
| + | |
| − | ipaddle--;
| + | |
| − | if (tdcl[index]>tdc_min&&tdcl[index]<tdc_min+tdc_width){
| + | |
| − | nmultiplicity++;
| + | |
| − | if (ipaddle == 1){
| + | |
| − | if (adc_c[index]>adc_cut && adc_c[paddleindex[ipaddle+1]] < adc_neighbor_cut ){
| + | |
| − | good_paddle[nmultiplicity-1]=ipaddle;
| + | |
| − | }
| + | |
| − | }else if (ipaddle == NUMPMT*NUMPADDLE){
| + | |
| − | if (adc_c[index]>adc_cut && adc_c[paddleindex[ipaddle-1]] < adc_neighbor_cut ){
| + | |
| − | good_paddle[nmultiplicity-1]=ipaddle;
| + | |
| − | }
| + | |
| − | }else{
| + | |
| − | if (adc_c[index]>adc_cut && adc_c[paddleindex[ipaddle-1]] < adc_neighbor_cut && adc_c[paddleindex[ipaddle+1]] < adc_neighbor_cut ){
| + | |
| − | good_paddle[nmultiplicity-1]=ipaddle;
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | if(nmultiplicity>0&&nmultiplicity<=multiplicity_cut){
| + | |
| − | for(Int_t icount=0;icount<nmultiplicity;icount++){
| + | |
| − | hoccupancy->Fill(good_paddle[icount]);
| + | |
| − | heatmap->Fill(good_paddle[icount],nmultiplicity);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | hmultiplicity->Fill(nmultiplicity);
| + | |
| − | }
| + | |
| − | | + | |
| − | cOCCUPANCY->Clear();
| + | |
| − | cOCCUPANCY->Divide(1,2) ;
| + | |
| − | | + | |
| − | title.Form("run_%d_OCCUPANCY_tdc_min_%d_max_%d.png",
| + | |
| − | run,tdc_min,tdc_min+tdc_width);
| + | |
| − | cOCCUPANCY->Print(title);
| + | |
| − | cOCCUPANCY->cd(1);
| + | |
| − | //gPad->SetLogy();
| + | |
| − | hoccupancy->Draw();
| + | |
| − | hoccupancy->GetXaxis()->SetNdivisions(NUMPMT,NUMPMT,0,0);
| + | |
| − | hoccupancy->SetLineColor(kBlue);
| + | |
| − | gPad->SetGridx();
| + | |
| − | TVirtualPad *c1_1 = cOCCUPANCY->cd(2);
| + | |
| − | c1_1->Divide(2,1);
| + | |
| − | c1_1->cd(1);
| + | |
| − | hmultiplicity->Draw();
| + | |
| − | hmultiplicity->SetLineColor(kBlue);
| + | |
| − | c1_1->cd(2);
| + | |
| − | heatmap->Draw("COLZ");
| + | |
| − | heatmap->GetXaxis()->SetNdivisions(NUMPMT,NUMPMT,0,0);
| + | |
| − | gPad->SetGridx();
| + | |
| − | return cOCCUPANCY;
| + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *plot_occupancy_single(Int_t event_id = 1, Int_t adc_cut=40, Int_t multiplicity_cut=12, Int_t tdc_min=750, Int_t tdc_width=200){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_Occupancy",run);
| + | |
| − | | + | |
| − | Int_t nbin=196;
| + | |
| − | Int_t min=1, max=197;
| + | |
| − | Int_t nbinm=11;
| + | |
| − | Int_t minm=-1, maxm=10;
| + | |
| − | TH1D *hoccupancy = new TH1D("hoccupancy","hoccupancy",nbin,min,max);
| + | |
| − | TH1D *hmultiplicity = new TH1D("hmultiplicity","hmultiplicity",nbinm,minm,maxm);
| + | |
| − | | + | |
| − | TCanvas *cOCCUPANCYE= new TCanvas("cOCCUPANCYE",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | Int_t nentries=n_events_to_analyze;
| + | |
| − |
| + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.9);
| + | |
| − | MyStyle->SetStatW(0.2);
| + | |
| − | | + | |
| − | //for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(event_id);
| + | |
| − |
| + | |
| − | Int_t nmultiplicity=0;
| + | |
| − | Int_t good_paddle[500];
| + | |
| − | for(Int_t icount=0;icount<500;icount++){good_paddle[icount]=-1;}
| + | |
| − | | + | |
| − | for(Int_t pmt=0; pmt<NUMPMT;pmt++){
| + | |
| − | Int_t ipaddle = (pmt+1)*NUMPADDLE+1;
| + | |
| − | for(Int_t pixel=0; pixel<NUMPIXEL;pixel++){
| + | |
| − | Int_t index = pmt*NUMPIXEL+pixel;
| + | |
| − | if (pixel!=pixel1[pmt]-1&&pixel!=pixel2[pmt]-1){
| + | |
| − | ipaddle--;
| + | |
| − | if (tdcl[index]>tdc_min&&tdcl[index]<tdc_min+tdc_width){
| + | |
| − | nmultiplicity++;
| + | |
| − | if (adc_c[index]>adc_cut){
| + | |
| − | good_paddle[nmultiplicity-1]=ipaddle;
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | if(nmultiplicity>0&&nmultiplicity<=multiplicity_cut){
| + | |
| − | for(Int_t icount=0;icount<nmultiplicity;icount++){hoccupancy->Fill(good_paddle[icount]);}
| + | |
| − | }
| + | |
| − | hmultiplicity->Fill(nmultiplicity);
| + | |
| − | //}
| + | |
| − | | + | |
| − | cOCCUPANCYE->Clear();
| + | |
| − | cOCCUPANCYE->Divide(1,2) ;
| + | |
| − | | + | |
| − | title.Form("run_%d_OCCUPANCYE_tdc_min_%d_max_%d.png",
| + | |
| − | run,tdc_min,tdc_min+tdc_width);
| + | |
| − | cOCCUPANCYE->Print(title);
| + | |
| − | cOCCUPANCYE->cd(1);
| + | |
| − | hoccupancy->Draw();
| + | |
| − | hoccupancy->GetXaxis()->SetNdivisions(NUMPMT,NUMPMT,0,0);
| + | |
| − | hoccupancy->SetLineColor(kBlue);
| + | |
| − | gPad->SetGridx();
| + | |
| − | cOCCUPANCYE->cd(2);
| + | |
| − | hmultiplicity->Draw();
| + | |
| − | hmultiplicity->SetLineColor(kBlue);
| + | |
| − | return cOCCUPANCYE;
| + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *plot_ratio(Int_t pmt=1, Int_t tdc_min=750, Int_t tdc_width=200){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_ADCRATIO",run);
| + | |
| − | TCanvas *cADCRATIO= new TCanvas("cADCRATIO",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | TH1D *htmpa[NUMPIXEL];//=new TH1D("htmp","htmp",nbin,min,max);
| + | |
| − | TH1D *htmpb[NUMPIXEL];//=new TH1D("htmp1","htmp1",nbin,min,max);
| + | |
| − | TH1D *htmpc[NUMPIXEL];//=new TH1D("htmp1","htmp1",nbin,min,max);
| + | |
| − | | + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.6);
| + | |
| − | MyStyle->SetStatW(0.4);
| + | |
| − | | + | |
| − | Int_t nbin=25;
| + | |
| − | Int_t min=0, max=100;
| + | |
| − | for(Int_t icounter=1;icounter<=NUMPIXEL;icounter++){
| + | |
| − | tmpentry.Form("htmpa%d", icounter);
| + | |
| − | htmpa[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | tmpentry.Form("htmpb%d", icounter);
| + | |
| − | htmpb[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | tmpentry.Form("htmpc%d", icounter);
| + | |
| − | htmpc[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | htmpa[icounter - 1]->SetLineColor(kBlue);
| + | |
| − | htmpb[icounter - 1]->SetLineColor(kRed);
| + | |
| − | htmpc[icounter - 1]->SetLineColor(kRed);
| + | |
| − | title.Form("Run %d ADC pmt %d, paddle %d: %d < tdc < %d",run,pmt,icounter,tdc_min,tdc_min+tdc_width);
| + | |
| − | htmpa[icounter - 1]->SetTitle(title);
| + | |
| − | htmpb[icounter - 1]->SetTitle(title);
| + | |
| − | htmpc[icounter - 1]->SetTitle(title);
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t nentries=n_events_to_analyze;
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − | for (Int_t index=1; index<=NUMPIXEL; index++){
| + | |
| − | Int_t ipaddle = (pmt-1)*NUMPIXEL+index-1;
| + | |
| − | htmpa[index-1]->Fill(adc_c[ipaddle]);
| + | |
| − | if(tdcl[ipaddle]>tdc_min&&tdcl[ipaddle]<tdc_min+tdc_width){
| + | |
| − | htmpb[index-1]->Fill(adc_c[ipaddle]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | cADCRATIO->Clear();
| + | |
| − | cADCRATIO->Divide(4,4);
| + | |
| − | | + | |
| − | //plot histos
| + | |
| − | Int_t icount=0;
| + | |
| − | TF1 *myfit = new TF1("myfit","1.0-0.5*ROOT::Math::erfc((x-[0])/[1])",0,1);
| + | |
| − | myfit->SetParName(0,"Mean");
| + | |
| − | myfit->SetParName(1,"Width");
| + | |
| − | | + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++){
| + | |
| − | | + | |
| − | if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1) {
| + | |
| − | | + | |
| − | //cout<<"into loop 2, i = " << i << endl;
| + | |
| − | | + | |
| − | cADCRATIO->cd( icount + 1 );
| + | |
| − | //gPad->SetLogy();
| + | |
| − | | + | |
| − | //cADC->Update();
| + | |
| − |
| + | |
| − | Int_t entries = htmpa[i]->GetEntries();
| + | |
| − | float mean = htmpa[i]->GetMean(1);
| + | |
| − | float RMS = htmpa[i]->GetRMS(1);
| + | |
| − | //cout << entries <<" "<< mean <<" "<< RMS <<endl;
| + | |
| − |
| + | |
| − | htmpb[i]->SetStats(0);
| + | |
| − | //current->Modified();
| + | |
| − | | + | |
| − | myfit->SetParameter(0,40.0);
| + | |
| − | myfit->SetParameter(1,10.0);
| + | |
| − | htmpc[i] = (TH1D*)htmpb[i]->Clone();
| + | |
| − | htmpc[i]->Divide(htmpa[i]);
| + | |
| − | htmpc[i]->Fit("myfit");
| + | |
| − | | + | |
| − | // write the mean and sigma of the fit
| + | |
| − | // on the canvas itself
| + | |
| − | | + | |
| − | float fitmean = myfit->GetParameter(0);
| + | |
| − | float fitsig = myfit->GetParameter(1);
| + | |
| − | | + | |
| − | // create strings to write on plot
| + | |
| − | TString lmean = Form("mean = %1.1f",fitmean);
| + | |
| − | TString lsig = Form("sigma = %1.1f",fitsig);
| + | |
| − |
| + | |
| − | // create the latex object to write
| + | |
| − | TLatex tl;
| + | |
| − | tl->SetTextSize(0.07);
| + | |
| − | tl->SetTextAngle(0);
| + | |
| − | tl->SetTextColor(1);
| + | |
| − | | + | |
| − | // NOTE: first two numbers are the x and y VALUES
| + | |
| − | // on the individual graph/histogram
| + | |
| − | tl->DrawLatex(30,0.85,lmean);
| + | |
| − | tl->DrawLatex(30,0.75,lsig);
| + | |
| − |
| + | |
| − | icount++;
| + | |
| − |
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | title.Form("run_%d_ADCRATIO_pmt_%d_tdc_min_%d_max_%d.png",
| + | |
| − | run,pmt,tdc_min,tdc_min+tdc_width);
| + | |
| − | cADCRATIO->Print(title);
| + | |
| − | cADCRATIO->cd(0);
| + | |
| − | return cADCRATIO;
| + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *plot_mean_ratio(Int_t pmt=1, Int_t tdc_min=750, Int_t tdc_width=200){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_ADCMEANRATIO",run);
| + | |
| − | TCanvas *cADCMEANRATIO= new TCanvas("cADCMEANRATIO",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | TH1D *htmpa[NUMPIXEL];//=new TH1D("htmp","htmp",nbin,min,max);
| + | |
| − | TH1D *htmpb[NUMPIXEL];//=new TH1D("htmp1","htmp1",nbin,min,max);
| + | |
| − | TH1D *htmpc[NUMPIXEL];//=new TH1D("htmp1","htmp1",nbin,min,max);
| + | |
| − | | + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.6);
| + | |
| − | MyStyle->SetStatW(0.4);
| + | |
| − | | + | |
| − | Int_t nbin=25;
| + | |
| − | Int_t min=0, max=100;
| + | |
| − | for(Int_t icounter=1;icounter<=NUMPIXEL;icounter++){
| + | |
| − | tmpentry.Form("htmpa%d", icounter);
| + | |
| − | htmpa[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | tmpentry.Form("htmpb%d", icounter);
| + | |
| − | htmpb[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | tmpentry.Form("htmpc%d", icounter);
| + | |
| − | htmpc[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | htmpa[icounter - 1]->SetLineColor(kBlue);
| + | |
| − | htmpb[icounter - 1]->SetLineColor(kRed);
| + | |
| − | htmpc[icounter - 1]->SetLineColor(kRed);
| + | |
| − | title.Form("Run %d ADC pmt %d, paddle %d: %d < tdc < %d",run,pmt,icounter,tdc_min,tdc_min+tdc_width);
| + | |
| − | htmpa[icounter - 1]->SetTitle(title);
| + | |
| − | htmpb[icounter - 1]->SetTitle(title);
| + | |
| − | htmpc[icounter - 1]->SetTitle(title);
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t nentries=n_events_to_analyze;
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − |
| + | |
| − | for (Int_t index=1; index<=NUMPIXEL; index++){
| + | |
| − | Int_t ipaddle = (pmt-1)*NUMPIXEL+index-1;
| + | |
| − | htmpa[index-1]->Fill(adc_c[ipaddle]);
| + | |
| − | if(tdcl[ipaddle]>tdc_min&&tdcl[ipaddle]<tdc_min+tdc_width){
| + | |
| − | htmpb[index-1]->Fill(adc_c[ipaddle]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | cADCMEANRATIO->Clear();
| + | |
| − | | + | |
| − | //plot histos
| + | |
| − | Int_t icount=0;
| + | |
| − | cADCMEANRATIO->cd();
| + | |
| − | | + | |
| − | Double_t mean[NUMPADDLE]={0,0,0,0,0,0,0,0,0,0,0,0,0,0};
| + | |
| − | Double_t sigma[NUMPADDLE]={1,1,1,1,1,1,1,1,1,1,1,1,1,1};
| + | |
| − | Double_t paddle[NUMPADDLE]={1,2,3,4,5,6,7,8,9,10,11,12,13,14};
| + | |
| − | Double_t epaddle[NUMPADDLE]={0,0,0,0,0,0,0,0,0,0,0,0,0,0};
| + | |
| − | | + | |
| − | TF1 *myfit = new TF1("myfit","1.0-0.5*ROOT::Math::erfc((x-[0])/[1])",0,1);
| + | |
| − | myfit->SetParName(0,"Mean");
| + | |
| − | myfit->SetParName(1,"Width");
| + | |
| − | | + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++){
| + | |
| − | | + | |
| − | if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1) {
| + | |
| − | | + | |
| − | htmpb[i]->SetStats(0);
| + | |
| − | | + | |
| − | myfit->SetParameter(0,40.0);
| + | |
| − | myfit->SetParameter(1,10.0);
| + | |
| − | htmpc[i] = (TH1D*)htmpb[i]->Clone();
| + | |
| − | htmpc[i]->Divide(htmpa[i]);
| + | |
| − | htmpc[i]->Fit("myfit");
| + | |
| − | | + | |
| − | Int_t numentries = htmpc[i]->GetEntries();
| + | |
| − | mean[icount]=myfit->GetParameter(0);
| + | |
| − | //sigma[icount]=myfit->GetParameter(1)/sqrt(numentries);
| + | |
| − | sigma[icount]=myfit->GetParameter(1);
| + | |
| − | | + | |
| − | icount++;
| + | |
| − |
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | gr = new TGraphErrors(NUMPADDLE,paddle,mean,epaddle,sigma);
| + | |
| − | gr->SetMarkerStyle(21);
| + | |
| − | gr->GetXaxis()->SetTitle("Paddle Number");
| + | |
| − | gr->GetYaxis()->SetTitle("50% Threshold (Good TDC)");
| + | |
| − | gr->GetHistogram()->SetMaximum(100);
| + | |
| − | gr->GetHistogram()->SetMinimum(0);
| + | |
| − | | + | |
| − | gr->Draw("AP");
| + | |
| − | | + | |
| − | cADCMEANRATIO->Update();
| + | |
| − | | + | |
| − | title.Form("run_%d_ADCMEANRATIO_pmt_%d_tdc_min_%d_max_%d.png",
| + | |
| − | run,pmt,tdc_min,tdc_min+tdc_width);
| + | |
| − | cADCMEANRATIO->Print(title);
| + | |
| − | cADCMEANRATIO->cd();
| + | |
| − | return cADCMEANRATIO;
| + | |
| − | | + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *plot_mean_adc(Int_t pmt=1, Int_t tdc_min=750, Int_t tdc_width=200){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_ADCMean",run);
| + | |
| − | TCanvas *cADCMean= new TCanvas("cADCMean",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | TH1D *htmpa[NUMPIXEL];//=new TH1D("htmp","htmp",nbin,min,max);
| + | |
| − | | + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.6);
| + | |
| − | MyStyle->SetStatW(0.4);
| + | |
| − | | + | |
| − | Int_t nbin=600;
| + | |
| − | Int_t min=-100, max=1000;
| + | |
| − | for(Int_t icounter=1;icounter<=NUMPIXEL;icounter++){
| + | |
| − | tmpentry.Form("htmpa%d", icounter);
| + | |
| − | htmpa[icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | title.Form("Run %d ADC pmt %d, paddle %d: %d < tdc < %d",run,pmt,icounter,tdc_min,tdc_min+tdc_width);
| + | |
| − | htmpa[icounter - 1]->SetTitle(title);
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t nentries=n_events_to_analyze;
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − | for (Int_t index=1; index<=NUMPIXEL; index++){
| + | |
| − | Int_t ipaddle = (pmt-1)*NUMPIXEL+index-1;
| + | |
| − | if(tdcl[ipaddle]>tdc_min&&tdcl[ipaddle]<tdc_min+tdc_width){
| + | |
| − | htmpa[index-1]->Fill(adc_c[ipaddle]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | cADCMean->Clear();
| + | |
| − | | + | |
| − | //plot histos
| + | |
| − | Int_t icount=0;
| + | |
| − | cADCMean->cd();
| + | |
| − | | + | |
| − | Double_t mean[NUMPIXEL]={0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
| + | |
| − | Double_t sigma[NUMPIXEL]={1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1};
| + | |
| − | Double_t pixel[NUMPIXEL]={1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16};
| + | |
| − | Double_t epixel[NUMPIXEL]={0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
| + | |
| − | | + | |
| − | Double_t mmean[2]={0,0};
| + | |
| − | Double_t msigma[2]={1,1};
| + | |
| − | Double_t mpixel[2]={0,0};
| + | |
| − | Double_t mepixel[2]={0,0};
| + | |
| − | mpixel[0]=pixel1[pmt-1];
| + | |
| − | mpixel[1]=pixel2[pmt-1];
| + | |
| − | | + | |
| − | Double_t yline[2]={100,100};
| + | |
| − | Double_t xline[2]={0,17};
| + | |
| − |
| + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++){
| + | |
| − | | + | |
| − | if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1) {
| + | |
| − | | + | |
| − | Int_t entries = htmpa[i]->GetEntries();
| + | |
| − | mean[i] = htmpa[i]->GetMean(1);
| + | |
| − | sigma[i] = htmpa[i]->GetRMS(1)/sqrt(entries);
| + | |
| − | cout << i+1 << " "<< entries <<" "<< mean[i] <<" "<< sigma[i] <<endl;
| + | |
| − |
| + | |
| − | icount++;
| + | |
| − |
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | gr = new TGraphErrors(NUMPIXEL,pixel,mean,epixel,sigma);
| + | |
| − | gr->SetMarkerStyle(21);
| + | |
| − | gr->GetXaxis()->SetTitle("Pixel Number");
| + | |
| − | gr->GetYaxis()->SetTitle("Mean ADC (Good TDC)");
| + | |
| − | gr2 = new TGraphErrors(2,mpixel,mmean,mepixel,msigma);
| + | |
| − | gr2->SetMarkerStyle(21);
| + | |
| − | gr2->SetMarkerColor(3);
| + | |
| − | gr3 = new TGraph(2,xline,yline);
| + | |
| − | gr3->SetLineColor(3);
| + | |
| − | gr3->SetLineWidth(4);
| + | |
| − | | + | |
| − | gr->Draw("AP");
| + | |
| − | gr2->Draw("P");
| + | |
| − | gr3->Draw("L");
| + | |
| − |
| + | |
| − | cADCMean->Update();
| + | |
| − | | + | |
| − | title.Form("run_%d_ADCMean_pmt_%d_tdc_min_%d_max_%d.png",
| + | |
| − | run,pmt,tdc_min,tdc_min+tdc_width);
| + | |
| − | cADCMean->Print(title);
| + | |
| − | cADCMean->cd(0);
| + | |
| − | return cADCMean;
| + | |
| − | }
| + | |
| − | | + | |
| − | TCanvas *calibrate_adc(Int_t ped = 1, Int_t tdc_min=750, Int_t tdc_width=200){
| + | |
| − | | + | |
| − | TString cut, draw, draw1, title;
| + | |
| − | title.Form("run_%d_ADCCalibrate",run);
| + | |
| − | TCanvas *cADCCalib= new TCanvas("cADCCalib",title,xcanvas,ycanvas);
| + | |
| − | | + | |
| − | TH1D *htmpa[NUMPMT][NUMPIXEL];//=new TH1D("htmp","htmp",nbin,min,max);
| + | |
| − | TH1D *htmpb[NUMPMT][NUMPIXEL];//=new TH1D("htmp1","htmp1",nbin,min,max);
| + | |
| − | TH1D *htmpc[NUMPMT][NUMPIXEL];//=new TH1D("htmp1","htmp1",nbin,min,max);
| + | |
| − | | + | |
| − | TString tmpentry;
| + | |
| − | MyStyle->SetStatX(0.9);
| + | |
| − | MyStyle->SetStatY(0.6);
| + | |
| − | MyStyle->SetStatW(0.4);
| + | |
| − | | + | |
| − | Int_t nbin=600;
| + | |
| − | Int_t min=-100, max=1000;
| + | |
| − | for (Int_t pmt=1; pmt<=NUMPMT; pmt++){
| + | |
| − | for(Int_t icounter=1;icounter<=NUMPIXEL;icounter++){
| + | |
| − | tmpentry.Form("htmpa_%d_%d", pmt,icounter);
| + | |
| − | htmpa[pmt-1][icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | tmpentry.Form("htmpb_%d_%d", pmt,icounter);
| + | |
| − | htmpb[pmt-1][icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | tmpentry.Form("htmpc_%d_%d", pmt,icounter);
| + | |
| − | htmpc[pmt-1][icounter - 1] = new TH1D(tmpentry,tmpentry,nbin,min,max);
| + | |
| − | htmpa[pmt-1][icounter - 1]->SetLineColor(kBlue);
| + | |
| − | htmpb[pmt-1][icounter - 1]->SetLineColor(kRed);
| + | |
| − | htmpc[pmt-1][icounter - 1]->SetLineColor(kGreen);
| + | |
| − | title.Form("Run %d ADC pmt %d, paddle %d: %d < tdc < %d",run,pmt,icounter,tdc_min,tdc_min+tdc_width);
| + | |
| − | htmpa[pmt-1][icounter - 1]->SetTitle(title);
| + | |
| − | htmpb[pmt-1][icounter - 1]->SetTitle(title);
| + | |
| − | htmpc[pmt-1][icounter - 1]->SetTitle(title);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | Int_t nentries=n_events_to_analyze;
| + | |
| − | | + | |
| − | for (Int_t id=1;id<=nentries;id++){
| + | |
| − | T->GetEntry(id);
| + | |
| − | for (Int_t pmt=1; pmt<=NUMPMT; pmt++){
| + | |
| − | for (Int_t index=1; index<=NUMPIXEL; index++){
| + | |
| − | Int_t ipaddle = (pmt-1)*NUMPIXEL+index-1;
| + | |
| − | htmpa[pmt-1][index-1]->Fill(adc[ipaddle]);
| + | |
| − | if(tdcl[ipaddle]<tdc_min||tdcl[ipaddle]>tdc_min+tdc_width){
| + | |
| − | htmpb[pmt-1][index-1]->Fill(adc[ipaddle]);
| + | |
| − | }else{
| + | |
| − | htmpc[pmt-1][index-1]->Fill(adc[ipaddle]);
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | }
| + | |
| − | | + | |
| − | cADCCalib->Clear();
| + | |
| − | cADCCalib->Divide(4,4);
| + | |
| − | | + | |
| − | //plot histos
| + | |
| − | for (Int_t pmt=1; pmt<=NUMPMT; pmt++){
| + | |
| − | Int_t icount=0;
| + | |
| − | for (Int_t i=0; i<NUMPIXEL; i++){
| + | |
| − | | + | |
| − | //if(i != pixel1[pmt-1]-1 && i != pixel2[pmt-1]-1) {
| + | |
| − | | + | |
| − | //cout<<"into loop 2, i = " << i << endl;
| + | |
| − | | + | |
| − | cADCCalib->cd( icount + 1 );
| + | |
| − | gPad->SetLogy();
| + | |
| − | | + | |
| − | cADCCalib->Update();
| + | |
| − |
| + | |
| − | Int_t entries = htmpb[pmt-1][i]->GetEntries();
| + | |
| − | float mean = htmpb[pmt-1][i]->GetMean(1);
| + | |
| − | float RMS = htmpb[pmt-1][i]->GetRMS(1);
| + | |
| − | Int_t imean = mean;
| + | |
| − | float mean2 = htmpc[pmt-1][i]->GetMean(1);
| + | |
| − | float gain;
| + | |
| − | if (mean2-mean > 2.0){
| + | |
| − | gain = 100.0/(mean2-mean);
| + | |
| − | }else{
| + | |
| − | gain = 1.0;
| + | |
| − | }
| + | |
| − | //cout << entries <<" "<< mean <<" "<< RMS <<endl;
| + | |
| − | if (ped == 1){
| + | |
| − | if(pmt==1&&i==0) cout << "C.cdetm1r.ped = ";
| + | |
| − | cout << imean << " ";
| + | |
| − | }else{
| + | |
| − | if(pmt==1&&i==0) cout << "C.cdetm1r.gain = ";
| + | |
| − | cout << gain << " ";
| + | |
| − | }
| + | |
| − |
| + | |
| − |
| + | |
| − | htmpa[pmt-1][i]->SetStats(0);
| + | |
| − | //current->Modified();
| + | |
| − | htmpa[pmt-1][i]->Draw();
| + | |
| − | htmpb[pmt-1][i]->Draw("same");
| + | |
| − | htmpc[pmt-1][i]->Draw("same");
| + | |
| − | | + | |
| − | icount++;
| + | |
| − |
| + | |
| − | //}
| + | |
| − | }
| + | |
| − | cout << endl;
| + | |
| − | }
| + | |
| − | | + | |
| − | title.Form("run_%d_ADCCalib_tdc_min_%d_max_%d.png",
| + | |
| − | run,tdc_min,tdc_min+tdc_width);
| + | |
| − | cADCCalib->Print(title);
| + | |
| − | cADCCalib->cd(0);
| + | |
| − | return cADCCalib;
| + | |
| − | }
| + | |
