// // Check for possible errors and shows raw distribution of variables coming from the DSP (no calibrations needed) - Emiliano Mocchiutti // // FCaloQLOOK.c version 1.11 (2006-07-17) // // The only input needed is the path to the directory created by YODA for the data file you want to analyze. // // Changelog: // // 1.09 - 1.11 (2006-07-17): Adapted to flight conditions. // // 1.07 - 1.09 (2006-05-29): Fixed bug in output filename when input is not in the form DW_YYMMDD_NNN. Bug in latchup alarm not coming from the status word, fixed. // // 1.06 - 1.07 (2006-05-23): Don't print warning in case of latchup not recognized from the status word, fixed. Added "expert" figures from 21 to 28 which appears only in the // case of latchup alarm. // // 1.05 - 1.06 (2006-03-22): Add optimize flag in compiling the script! // // 1.04 - 1.05 (2006-03-22): Corrected wrong .C files. // // 1.03 - 1.04 (2006-03-20): Documentation updated. // // 1.02 - 1.03 (2006-03-20): Changed name of shared libraries (from FCaloQLOOK_cxx.so to libFCaloQLOOK.so). // // 1.01 - 1.02 (2006-03-13): Include files from YODA without "event" directory. // // 1.00 - 1.01 (2006-02-28): Works on YODA v6 output (single file), does not require anymore calocommon package. // // 0.00 - 1.00 (2006-02-28): Clone of CaloQLOOK.c (ground software). // #include <iostream> #include <fstream> // #include <TObject.h> #include <TString.h> #include <TFile.h> #include <TCanvas.h> #include <TH1.h> #include <TH1F.h> #include <TH2D.h> #include <TLatex.h> #include <TPad.h> #include <TPolyLine.h> #include <TStyle.h> #include <TSystem.h> // #include <PamelaRun.h> #include <physics/calorimeter/CalorimeterEvent.h> // #include <FCaloQLOOKfun.h> // using namespace std; // void stringcopy(TString& s1, const TString& s2, Int_t from=0, Int_t to=0){ if ( to == 0 ){ Int_t t2length = s2.Length(); s1 = ""; to = t2length; }; for (Int_t i = from; i<to; i++){ s1.Append(s2[i],1); }; } void stringappend(TString& s1, const TString& s2){ Int_t t2length = s2.Length(); for (Int_t i = 0; i<t2length; i++){ s1.Append(s2[i],1); }; } TString getFilename(const TString filename){ // const string fil = gSystem->BaseName(filename.Data()); Int_t posiz = fil.find(".root"); // TString file2; if ( posiz == -1 ){ file2 = gSystem->BaseName(filename.Data()); } else { Int_t posiz2 = 0; stringcopy(file2,gSystem->BaseName(filename.Data()),posiz2,posiz); TString pdat(".dat"); stringappend(file2,pdat); }; return file2; } void FCaloQLOOK(TString filename, Int_t fromevent=0, Int_t toevent=0, TString outDir="", TString saveas = "ps"){ gStyle->SetPaperSize(19.,25.); const char* startingdir = gSystem->WorkingDirectory(); // printf(" basename is %s \n",gSystem->BaseName(filename.Data())); // printf(" getfilename is %s \n",getFilename(filename).Data()); if ( !strcmp(outDir.Data(),"") ) outDir = startingdir; // // ifstream myfile; myfile.open(filename.Data()); if ( !myfile ){ printf(" %s :no such file, exiting...\n\n",filename.Data()); return; }; myfile.close(); // TFile *File = new TFile(filename.Data()); TTree *tr = (TTree*)File->Get("Physics"); if ( !tr ) { printf(" Physics : no such tree in %s \n",filename.Data()); printf(" Exiting, are you sure this is a LEVEL0 not corrupted file? \n\n"); return; }; // // Define variables // Float_t ctshitthr = 0.65; Float_t cbasethr = 0.95; Float_t cdexythr = 0.995; //Float_t cdexythr = 0.90; Float_t cdexycthr = 0.95; Float_t h1rth = 0.90; Int_t calevnth = 13; // pamela::calorimeter::CalorimeterEvent *ce = 0; pamela::EventHeader *eh = 0; pamela::PscuHeader *ph = 0; // UInt_t found; tr->SetBranchStatus("*",0,&found); //disable all branches // found = 0; tr->SetBranchStatus("iev*",1,&found); // printf("enabled %i branches \n",found); found = 0; tr->SetBranchStatus("dexy*",1,&found); // printf("enabled %i branches \n",found); found = 0; tr->SetBranchStatus("stwerr*",1,&found); // printf("enabled %i branches \n",found); found = 0; tr->SetBranchStatus("perror*",1,&found); // printf("enabled %i branches \n",found); found = 0; tr->SetBranchStatus("cal*",1,&found); // printf("enabled %i branches \n",found); found = 0; tr->SetBranchStatus("base*",1,&found); //printf("enabled %i branches \n",found); found = 0; tr->SetBranchStatus("Pscu*",1,&found); //printf("enabled %i branches \n",found); found = 0; tr->SetBranchStatus("Calorimeter*",1,&found); //printf("calo enabled %i branches \n",found); found = 0; tr->SetBranchStatus("Header*",1,&found); //printf("head enabled %i branches \n",found); // tr->SetBranchAddress("Calorimeter", &ce); tr->SetBranchAddress("Header", &eh); // Long64_t nevents = tr->GetEntries(); if ( nevents < 1 ) { printf("The file is empty!\n"); return; }; Int_t const size = nevents; Double_t iev2[size+1]; Int_t ver[4][23], lver[4][23],shit[4][11], rshit[4][11], se,bl, pl, fulldiff[4],cshit, cpre[4],nullsh[4][11], lalarm[4]; se = 0; pl = 0; Float_t allbase = 0.; Int_t tshit, trshit, compdata,rawdata, errorfull, calevn1, calevn2, calevn3, calevn4, pshit[4][11]; tshit = 0; trshit = 0; char *sezione; Bool_t check = false; sezione = ""; stringstream errore; errore.str(""); // TString fififile = getFilename(filename); const char *file = fififile; // // book histograms // TH1F *h1; TH1F *h1r; TH1F *hdiff; TH1F *Baseline; TH1F *Dexy; TH1F *Dexyc; TH1F *Calstriphit; TH1F *calev01; TH1F *calev23; TH1F *calev03; TH1F *lupstw[4]; TH1F *lup[4]; TH2D *calev2; h1 = new TH1F("C7","Strips hit, compress mode",100,0.,200.); h1r = new TH1F("C8","Strips hit, raw mode",100,-1,1); h1r->SetBit(TH1F::kCanRebin); hdiff = new TH1F("C6","Differences in full mode",100,-3,3); hdiff->SetBit(TH1F::kCanRebin); Baseline = new TH1F("C9","baselines",100,0,3); Baseline->SetBit(TH1F::kCanRebin); Dexy = new TH1F("C12","dexy",100,0.,10000.); Dexy->SetBit(TH1F::kCanRebin); Dexyc = new TH1F("C10","dexyc",100,-100.,37000.); // Dexyc->SetBit(TH1F::kCanRebin); Calstriphit = new TH1F("C11","calstriphit[1:4]",100,0.,200.); calev01 = new TH1F("C3","|calevnum(1)-calevnum(2)|",100,0,1); calev01->SetBit(TH1F::kCanRebin); calev23 = new TH1F("C4","|calevnum(3)-calevnum(4)|",100,0,1); calev23->SetBit(TH1F::kCanRebin); calev03 = new TH1F("C5","|calevnum(1)-calevnum(4)|",100,0,1); calev03->SetBit(TH1F::kCanRebin); calev2 = new TH2D("C2","calevnum(2)%iev",3000,0.,2.,100,0.,2.); calev2->SetBit(TH2D::kCanRebin); stringstream oss; stringstream noss; for ( Int_t i=0; i<4; i++){ oss.str(""); oss << "OBT stw latch up section " << i; noss.str(""); noss << "C" << 21+i; lupstw[i] = new TH1F(noss.str().c_str(),oss.str().c_str(),100,-1,1); lupstw[i]->SetBit(TH1F::kCanRebin); oss.str(""); oss << "OBT no stw latch up section " << i; noss.str(""); noss << "C" << 25+i; lup[i] = new TH1F(noss.str().c_str(),oss.str().c_str(),100,-1,1); lup[i]->SetBit(TH1F::kCanRebin); }; // // run over each event and take out some variables // if ( fromevent > toevent && toevent ){ printf("It must be fromevent < toevent \n"); return; }; if ( fromevent > nevents+1 || fromevent < 0 ) { printf("You can choose fromevent between 0 (all) and %i \n",(int)nevents+1); return; }; if ( toevent > nevents+1 || toevent < 0 ) { printf("You can choose toevent between 0 (all) and %i \n",(int)nevents+1); return; }; Int_t minevent = 0; Int_t maxevent = 0; if ( fromevent == 0 ) { minevent = 0; maxevent = nevents -1; } else { minevent = fromevent - 1; if ( toevent > 0 ){ maxevent = toevent - 1; } else { maxevent = fromevent - 1; }; }; for (Int_t k = 0; k < 4; k++ ){ for (Int_t m = 0; m < 23 ; m++ ){ ver[k][m] = 0 ; lver[k][m] = 0 ; cpre[k]=0; if ( m < 11 ) { pshit[k][m] = 0; nullsh[k][m] = 0; }; lalarm[k] = 0; }; }; Int_t pdone, bdone; pdone = 0; bdone = 0; bool isCOMP = 0; bool isFULL = 0; bool isRAW = 0; Int_t alldexy=0; Int_t alldexy2=0; Int_t planebases=0; // Int_t stri=0; Int_t fcheck = 0; Int_t cestw=0; Int_t cmask = 127 ; compdata = 0; rawdata = 0; errorfull = 0; calevn1 = 0; calevn2 = 0; calevn3 = 0; calevn4 = 0; fulldiff[0] = 0; fulldiff[1] = 0; fulldiff[2] = 0; fulldiff[3] = 0; printf("\n Processed events: \n\n"); unsigned short int calev0=0; unsigned short int oldcalev0; unsigned short int calev1=0; unsigned short int oldcalev1; unsigned short int calevv2=0; unsigned short int oldcalev2; unsigned short int calev3=0; unsigned short int oldcalev3; Int_t i = minevent; Float_t headc = 0.; Float_t headco = 0.; Bool_t h1rcheck = false; Int_t h1rin = 0; Int_t h1rout = 0; Int_t intshit = 0; Int_t outtshit = 0; Int_t incshit = 0; Int_t outcshit = 0; Int_t inbase = 0; Int_t outbase = 0; Int_t indexy = 0; Int_t outdexy = 0; Int_t indexyc = 0; Int_t outdexyc = 0; Int_t obt = 0; Int_t minobt[4]; Int_t maxobt[4]; Int_t swminobt[4]; Int_t swmaxobt[4]; Bool_t firstobt[4]; Bool_t swfirstobt[4]; for ( Int_t i = 0; i<4; i++){ minobt[i] = 0; maxobt[i] = 0; swminobt[i] = 0; swmaxobt[i] = 0; firstobt[i] = true; swfirstobt[i] = true; }; while ( i < maxevent+1){ tshit = 0; trshit = 0; tr->GetEntry(i); iev2[i] = ce->iev; // ph = eh->GetPscuHeader(); headco = headc; headc = ph->GetCounter(); obt = ph->GetOrbitalTime(); if ( (int)abs((int)(ce->calevnum[0]-ce->calevnum[1])) ) calevn1++; if ( (int)abs((int)(ce->calevnum[2]-ce->calevnum[3])) ) calevn2++; if ( (int)abs((int)(ce->calevnum[0]-ce->calevnum[3])) ) calevn3++; // calevn1 += (int)abs((int)(ce->calevnum[0]-ce->calevnum[1])); // calevn2 += (int)abs((int)(ce->calevnum[2]-ce->calevnum[3])); // calevn3 += (int)abs((int)(ce->calevnum[0]-ce->calevnum[3])); calev01->Fill(abs((int)(ce->calevnum[0]-ce->calevnum[1]))); calev23->Fill(abs((int)(ce->calevnum[2]-ce->calevnum[3]))); calev03->Fill(abs((int)(ce->calevnum[0]-ce->calevnum[3]))); unsigned short calvnm2 = (unsigned short)ce->calevnum[2]; calev2->Fill(ce->iev,calvnm2); calevn4 += (int)ce->calevnum[2]; oldcalev0 = calev0; calev0 = (int)ce->calevnum[0]; oldcalev1 = calev1; calev1 = (int)ce->calevnum[1]; oldcalev2 = calevv2; calevv2 = (int)ce->calevnum[2]; oldcalev3 = calev3; calev3 = (int)ce->calevnum[3]; goto jumpprintout; if ( (headc - headco -1.) == 0. && ((calev0 - oldcalev0 - 1) > 0 || (calev1 - oldcalev1 - 1) > 0 || (calevv2 - oldcalev2 - 1) > 0 || (calev3 - oldcalev3 - 1) > 0) ){ if ( i != minevent ) { printf(" %f 0 Event %i: %i\n",headco,i,oldcalev0); printf(" %f 0 Event %i: %i\n",headc,i+1,calev0); printf(" %f 1 Event %i: %i\n",headco,i,oldcalev1); printf(" %f 1 Event %i: %i\n",headc,i+1,calev1); printf(" %f 2 Event %i: %i\n",headco,i,oldcalev2); printf(" %f 2 Event %i: %i\n",headc,i+1,calevv2); printf(" %f 3 Event %i: %i\n",headco,i,oldcalev3); printf(" %f 3 Event %i: %i\n",headc,i+1,calev3); isRAW = 0; if ( ce->stwerr[0] & (1 << 3)) isRAW = 1; if ( (calev0 - oldcalev0 - 1) > 0 && !isRAW && ce->perror[0] != 129 && oldcalev0 != 0) { ver[0][10]++; }; isRAW = 0; if ( ce->stwerr[0] & (1 << 3)) isRAW = 1; if ( (calev1 - oldcalev1 - 1) > 0 && !isRAW && ce->perror[1] != 129 && oldcalev1 != 0 ){ ver[1][10]++; }; isRAW = 0; if ( ce->stwerr[0] & (1 << 3)) isRAW = 1; if ( (calevv2 - oldcalev2 - 1) > 0 && !isRAW && ce->perror[2] != 129 && oldcalev2 != 0 ){ ver[2][10]++; }; isRAW = 0; if ( ce->stwerr[0] & (1 << 3)) isRAW = 1; if ( (calev3 - oldcalev3 - 1) > 0 && !isRAW && ce->perror[3] != 129 && oldcalev3 != 0 ){ ver[3][10]++; }; } }; jumpprintout: memcpy(shit, nullsh, sizeof(nullsh)); memcpy(rshit, nullsh, sizeof(nullsh)); for (Int_t l = 0; l < 2; l++ ){ for (Int_t ii = 0; ii < 22; ii++){ if (l == 0 && ii%2 == 0){ se = 3; pl = ii/2; pdone = 0; }; if (l == 0 && ii%2 != 0){ se = 2; pl = (ii-1)/2; pdone = 0; }; if (l == 1 && ii%2 == 0){ se = 0; pl = ii/2; pdone = 0; }; if (l == 1 && ii%2 != 0){ se = 1; pl = (ii-1)/2; pdone = 0; }; isCOMP = 0; isFULL = 0; isRAW = 0; if ( ce->stwerr[se] & (1 << 16) ) isCOMP = 1; if ( ce->stwerr[se] & (1 << 17) ) isFULL = 1; if ( ce->stwerr[se] & (1 << 3) ) isRAW = 1; bl = -1; for (Int_t kk = 0; kk < 96 ; kk++ ){ if ( kk%16 == 0 ){ bdone = 0; bl++; allbase = ce->base[l][ii][bl]; // alldexy=0; // alldexy2=0; // stri=0; // for (Int_t e = 0; e < 16 ; e++ ){ // stri = e + 16 * bl; // alldexy += (int)ce->dexyc[l][ii][stri]; // alldexy2 += (int)ce->dexy[l][ii][stri]; // }; }; // if ( kk == 0 ){ planebases = 0; alldexy=0; alldexy2=0; for (Int_t e = 0; e < 96 ; e++ ){ if ( e < 6 ) planebases += (int)ce->base[l][ii][e]; alldexy += (int)ce->dexyc[l][ii][e]; alldexy2 += (int)ce->dexy[l][ii][e]; }; }; // if ( !isRAW ) { // if ( !pdone ){ // if ( (ce->base[l][ii][bl]>32000 || ce->base[l][ii][bl] == 0 ) && ( alldexy > 512000 || alldexy == 0) && ce->perror[se] == 0 ) { if ( (planebases>192000 || planebases == 0 ) && ( alldexy > 3072000 || alldexy == 0) && ce->perror[se] == 0 ) { pdone = 1; pshit[se][pl]++ ; if ( (ce->stwerr[se] & (1 << 4)) == 0 ) { lalarm[se]++; lup[se]->Fill(obt); if ( firstobt[se] ) minobt[se] = obt; if ( obt > maxobt[se] ) maxobt[se] = obt; }; lver[se][2]++ ; }; } // if ( !bdone ){ Baseline->Fill(ce->base[l][ii][bl]); if ( ce->base[l][ii][bl] > 2000. && ce->base[l][ii][bl] < 4500. ){ inbase++; } else { outbase++; }; bdone = 1; }; // if (ce->dexyc[l][ii][kk] > 0 ) { shit[se][pl]++ ; compdata = 1; }; // if ( isFULL ) { if ( ce->dexy[l][ii][kk] > 0 && ce->dexy[l][ii][kk] < 32000 ) { rshit[se][pl]++ ; rawdata = 1; }; Dexy->Fill(ce->dexy[l][ii][kk]); if ( ce->dexy[l][ii][kk] < 2000. || ce->dexy[l][ii][kk] > 5000. ) { outdexy++; } else { indexy++; }; if ( ce->dexyc[l][ii][kk]>0 && (ce->dexy[l][ii][kk]-ce->dexyc[l][ii][kk]) != 0 && ce->perror[se] != 132 && (ce->base[l][ii][bl] == 0 || ce->base[l][ii][bl] >32000) ){ hdiff->Fill(ce->dexyc[l][ii][kk]-ce->dexy[l][ii][kk]); fulldiff[se]++; errorfull = 1; }; }; // Dexyc->Fill(ce->dexyc[l][ii][kk]); if ( ce->dexyc[l][ii][kk] != 0. && ( ce->dexyc[l][ii][kk] < 2000. || ce->dexyc[l][ii][kk] > 5000. ) ) { outdexyc++; } else { indexyc++; }; // if (ce->dexyc[l][ii][kk] < 0 ) ver[se][21]++ ; // if ( allbase == 0. || allbase > 32000.) { fcheck = 0; for (Int_t nn = bl*16; nn < (bl+1)*16 ; nn++){ if ( ce->dexyc[l][ii][nn] > 0. && ce->dexyc[l][ii][nn] < 32000. ) fcheck++; }; if ( fcheck ) { cpre[se]++; allbase = 1.; }; }; // } else { if ( !pdone ){ if ( ( alldexy2 > 3072000 || alldexy2 == 0) && ce->perror[se] == 0 ) { // if ( (alldexy2>512000 || alldexy2 == 0) && ce->perror[se] == 0 ) { pdone = 1; pshit[se][pl]++ ; if ( (ce->stwerr[se] & (1 << 4)) == 0 ){ lalarm[se]++; lup[se]->Fill(obt); if ( firstobt[se] ) minobt[se] = obt; if ( obt > maxobt[se] ) maxobt[se] = obt; }; lver[se][2]++ ; }; }; if ( ce->dexy[l][ii][kk] > 0 && ce->dexy[l][ii][kk] < 32000 ) { rshit[se][pl]++ ; rawdata = 1; }; Dexy->Fill(ce->dexy[l][ii][kk]); if ( ce->dexy[l][ii][kk] < 2000. || ce->dexy[l][ii][kk] > 5000. ) { outdexy++; } else { indexy++; }; }; }; }; }; cshit = 0; for (Int_t k = 0; k < 4 ; k++ ){ isCOMP = 0; isFULL = 0; if ( ce->stwerr[se] & (1 << 16) ) isCOMP = 1; if ( ce->stwerr[se] & (1 << 17) ) isFULL = 1; if ( isCOMP ) ver[k][16]++; if ( isFULL ) ver[k][17]++; cshit += (int)ce->calstriphit[k]; cestw=0; if ( ce->stwerr[k] ) cestw = ce->stwerr[k] & cmask ; if ( cestw ){ if ( cestw & (1 << 0) ) ver[k][6]++ ; if ( cestw & (1 << 1) ) ver[k][5]++ ; if ( cestw & (1 << 2) ) ver[k][4]++ ; if ( cestw & (1 << 3) ) ver[k][3]++ ; if ( cestw & (1 << 4) ){ ver[k][2]++ ; lupstw[k]->Fill(obt); if ( swfirstobt[k] ) swminobt[k] = obt; if ( obt > swmaxobt[k] ) swmaxobt[k] = obt; }; if ( cestw & (1 << 5) ) ver[k][1]++ ; if ( cestw & (1 << 6) ) ver[k][0]++ ; }; if ( ce->perror[k] != 0. ){ if (ce->perror[k] == 128) ver[k][7]++ ; if (ce->perror[k] == 129) ver[k][8]++ ; if (ce->perror[k] == 132) ver[k][11]++ ; if (ce->perror[k] == 133) ver[k][12]++ ; if (ce->perror[k] == 134) ver[k][13]++ ; if (ce->perror[k] == 135) ver[k][14]++ ; if (ce->perror[k] == 136) ver[k][15]++ ; if (ce->perror[k] == 139) ver[k][18]++ ; if (ce->perror[k] == 140) ver[k][19]++ ; if (ce->perror[k] == 141) ver[k][20]++ ; if (ce->perror[k] == 142) ver[k][22]++ ; }; for (Int_t kk = 0; kk < 11 ; kk++ ){ tshit += shit[k][kk]; trshit += rshit[k][kk]; }; }; if (isCOMP || isFULL){ Calstriphit->Fill(cshit); // if ( (cshit > 0 && cshit < 25) || (cshit > 40 && cshit < 80) ){ if ( (cshit > 10 && cshit < 100) ){ incshit++; } else { outcshit++; }; // if ( tshit>0 ) h1->Fill(tshit); // if ( (tshit > 0 && tshit < 25) || (tshit > 40 && tshit < 80) ){ if ( (tshit > 10 && tshit < 100) ){ intshit++; } else { outtshit++; }; }; if ( trshit>0 ) { h1r->Fill(trshit); if ( trshit < 4210 ){ h1rout++; // printf("ma come... trshit %i \n",trshit); } else { h1rin++; }; }; if ( i%1000 == 0 && i > 0 ) printf("%iK\n",i/1000); i++; }; printf("\n"); if ( (float)h1rout/((float)h1rin+(float)h1rout) > h1rth ){ h1rcheck = true; }; // // output figures, first sheet: // gStyle->SetOptDate(0); gStyle->SetOptStat(1111); TCanvas *figura = new TCanvas("Calorimeter_Detector_Report_1/3", "Calorimeter_Detector_Report_1/3", 1100, 900); figura->SetFillColor(10); figura->Range(0,0,100,100); errore.str(""); errore << "EXPERT -- File: " << file; errore << " "; TLatex *t=new TLatex(); t->SetTextFont(32); t->SetTextColor(1); t->SetTextAlign(12); t->SetTextSize(0.015); t->DrawLatex(2.,99.,errore.str().c_str()); TPad *pd5; TPad *pd6; TPad *pd7; TPad *pd8; TPad *pad1; TPad *pad2; TPad *pad3; TPad *pad4; TPad *pd1; TPad *pd2; TPad *pd3; TPad *pd4; // Double_t h1max; TPolyLine *banda1; TPolyLine *banda2; Float_t ctshit = 0; Float_t ccshit = 0; Float_t cbase = 0; Float_t cdexy = 0; Float_t cdexyc = 0; // if ( compdata && rawdata ){ ctshit = (float)intshit/((float)outtshit + (float)intshit); Int_t pd5col = 10; Int_t pd6col = 10; if ( ctshit < ctshitthr ) { check = true; pd5col = 45; }; if ( h1rcheck ) { check = true; pd6col = 45; }; pd5 = new TPad("pd5","This is pad5",0.51,0.51,0.98,0.98,pd5col); pd6 = new TPad("pd6","This is pad6",0.51,0.02,0.98,0.49,pd6col); pd5->SetTicks(); pd6->SetTicks(); pd5->Range(0,0,100,100); pd6->Range(0,0,100,100); figura->cd(); pd5->Draw(); pd5->cd(); h1->SetXTitle("Number of hit (dexyc>0)"); h1->SetYTitle("Number of events"); h1->Draw(); // h1max = h1->GetMaximum()*1.05; Double_t xc[4] = {10.,100.,100.,10.}; Double_t yc[4] = {0.,0.,h1max,h1max}; banda1 = new TPolyLine(4,xc,yc); banda1->SetLineColor(5); banda1->SetFillColor(5); banda1->SetLineWidth(1); banda1->Draw("fSAME"); // Double_t xd[4] = {40.,80.,80.,40.}; // Double_t yd[4] = {0.,0.,h1max,h1max}; // banda2 = new TPolyLine(4,xd,yd); // banda2->SetLineColor(5); // banda2->SetFillColor(5); // banda2->SetLineWidth(1); //banda2->Draw("fSAME"); h1->Draw("SAME"); // figura->cd(); pd6->Draw(); pd6->cd(); gPad->SetLogy(); h1r->SetXTitle("Number of hit (dexy>0)"); h1r->SetYTitle("Number of events"); h1r->Draw(); h1max = h1r->GetMaximum()*2.05; Double_t xg[4] = {4220.,4230.,4230.,4220.}; Double_t yg[4] = {0.,0.,h1max,h1max}; banda1 = new TPolyLine(4,xg,yg); banda1->SetLineColor(5); banda1->SetFillColor(5); banda1->SetLineWidth(1); banda1->Draw("fSAME"); h1r->Draw("SAME"); }; if ( (compdata && !rawdata) || (!compdata && rawdata) ){ if ( tshit !=0 ) { ctshit = (float)intshit/((float)outtshit + (float)intshit); Int_t pd5col = 10; if ( ctshit < ctshitthr ) { check = true; pd5col = 45; }; pd5 = new TPad("pd5","This is pad5",0.51,0.02,0.98,0.98,pd5col); pd5->Range(0,0,100,100); pd5->SetTicks(); pd5->Draw(); pd5->cd(); h1->SetXTitle("Number of hit (dexyc>0)"); h1->SetYTitle("Number of events"); h1->Draw(); h1max = h1->GetMaximum()*1.05; Double_t xe[4] = {10.,100.,100.,10.}; Double_t ye[4] = {0.,0.,h1max,h1max}; banda1 = new TPolyLine(4,xe,ye); banda1->SetLineColor(5); banda1->SetFillColor(5); banda1->SetLineWidth(1); banda1->Draw("fSAME"); // Double_t xf[4] = {40.,80.,80.,40.}; // Double_t yf[4] = {0.,0.,h1max,h1max}; // banda2 = new TPolyLine(4,xf,yf); // banda2->SetLineColor(5); // banda2->SetFillColor(5); // banda2->SetLineWidth(1); // banda2->Draw("fSAME"); h1->Draw("SAME"); }; if ( trshit !=0 ) { Int_t pd5col = 10; if ( h1rcheck ) { check = true; pd5col = 45; }; pd5 = new TPad("pd5","This is pad5",0.51,0.02,0.98,0.98,pd5col); pd5->Range(0,0,100,100); pd5->SetTicks(); pd5->Draw(); pd5->cd(); gPad->SetLogy(); h1r->SetXTitle("Number of hit (dexyc>0 or dexy>0)"); h1r->SetYTitle("Number of events"); h1r->Draw(); h1max = h1r->GetMaximum()*2.05; Double_t xh[4] = {4220.,4230.,4230.,4220.}; Double_t yh[4] = {0.,0.,h1max,h1max}; banda1 = new TPolyLine(4,xh,yh); banda1->SetLineColor(5); banda1->SetFillColor(5); banda1->SetLineWidth(1); banda1->Draw("fSAME"); h1r->Draw("SAME"); }; }; if ( !errorfull ) { if ( calevn1 > calevnth || calevn2 > calevnth || calevn3 > calevnth ) { pd1 = new TPad("pd1","This is pad1",0.02,0.51,0.24,0.98,45); pd2 = new TPad("pd2","This is pad2",0.26,0.51,0.49,0.98,45); pd3 = new TPad("pd3","This is pad3",0.02,0.02,0.24,0.49,45); pd4 = new TPad("pd4","This is pad4",0.26,0.02,0.49,0.49,45); figura->cd(); // printf("ou2! \n"); check = true; pd1->Range(0,0,100,100); pd2->Range(0,0,100,100); pd3->Range(0,0,100,100); pd4->Range(0,0,100,100); pd1->Draw(); pd2->Draw(); pd3->Draw(); pd4->Draw(); pd1->cd(); gPad->SetLogy(); calev01->SetXTitle("delta calevnum 1-2"); calev01->SetYTitle("Number of events"); calev01->Draw(); pd2->cd(); gPad->SetLogy(); calev23->SetXTitle("delta calevnum 3-4"); calev23->SetYTitle("Number of events"); calev23->Draw(); pd3->cd(); gPad->SetLogy(); calev03->SetXTitle("delta calevnum 1-4"); calev03->SetYTitle("Number of events"); calev03->Draw(); pd4->cd(); calev2->SetYTitle("calevnum 2"); calev2->SetXTitle("iev"); calev2->Draw(); } else { if ( calevn4 ) { pd3 = new TPad("pd3","This is pad3",0.02,0.02,0.49,0.98,10); figura->cd(); pd3->Range(0,0,100,100); pd3->Draw(); pd3->cd(); calev2->SetYTitle("calevnum 2"); calev2->SetXTitle("iev"); calev2->Draw(); } else { figura->cd(); t=new TLatex(); t->SetTextFont(32); t->SetTextColor(1); t->SetTextSize(0.04); t->SetTextAlign(12); t->DrawLatex(10.,49.,"No calevnum infos from DSP"); }; }; } else { figura->cd(); check = true; // printf("nou! \n"); pd1 = new TPad("pd1","This is pad1",0.02,0.02,0.49,0.98,45); pd1->Range(0,0,100,100); pd1->Draw(); pd1->cd(); gPad->SetLogy(); hdiff->SetXTitle("Differences in FULL mode"); hdiff->SetYTitle("Number of events"); hdiff->Draw(); }; // // output figures, second sheet: // TCanvas *figura2 = new TCanvas("Calorimeter_Detector_Report_2/3","Calorimeter_Detector_Report_2/3", 1100, 900); if ( compdata ){ Int_t pd1col = 10; Int_t pd2col = 10; Int_t pd3col = 10; Int_t pd4col = 10; if ( (outdexy+indexy) ) { cdexy = (float)indexy/((float)outdexy + (float)indexy); if ( cdexy < cdexythr ) { check = true; pd4col = 45; }; }; if ( (outcshit+incshit) ) { ccshit = (float)incshit/((float)outcshit + (float)incshit); if ( ccshit < ctshitthr ) { check = true; pd3col = 45; }; }; if ( (outdexyc+indexyc) ) { cdexyc = (float)indexyc/((float)outdexyc + (float)indexyc); if ( cdexyc < cdexycthr ) { check = true; pd2col = 45; }; }; if ( (outbase+inbase) ) { cbase = (float)inbase/((float)outbase + (float)inbase); if ( cbase < cbasethr ) { check = true; pd1col = 45; }; }; figura2->SetFillColor(10); figura2->Range(0,0,100,100); errore.str(""); errore << "EXPERT -- File: " << file; errore << " "; t=new TLatex(); t->SetTextFont(32); t->SetTextColor(1); t->SetTextAlign(12); t->SetTextSize(0.015); t->DrawLatex(2.,99.,errore.str().c_str()); pd1 = new TPad("pd1","This is pad1",0.02,0.51,0.49,0.98,pd1col); pd2 = new TPad("pd2","This is pad2",0.51,0.51,0.98,0.98,pd2col); pd3 = new TPad("pd3","This is pad3",0.02,0.02,0.49,0.49,pd3col); pd4 = new TPad("pd4","This is pad4",0.51,0.02,0.98,0.49,pd4col); figura2->cd(); pd1->Range(0,0,100,100); pd2->Range(0,0,100,100); pd3->Range(0,0,100,100); pd4->Range(0,0,100,100); pd1->SetTicks(); pd2->SetTicks(); pd3->SetTicks(); pd4->SetTicks(); pd1->Draw(); pd2->Draw(); pd3->Draw(); pd4->Draw(); pd1->cd(); Baseline->SetXTitle("ADC channels"); Baseline->SetYTitle("Number of events"); Baseline->Draw(); h1max = Baseline->GetMaximum()*1.05; Double_t xc[4] = {2000.,4500.,4500.,2000.}; Double_t yc[4] = {0.,0.,h1max,h1max}; banda1 = new TPolyLine(4,xc,yc); banda1->SetLineColor(5); banda1->SetFillColor(5); banda1->SetLineWidth(1); banda1->Draw("fSAME"); Baseline->Draw("SAME"); // pd2->cd(); gPad->SetLogy(); Dexyc->SetXTitle("ADC channels"); Dexyc->SetYTitle("Number of events"); Dexyc->Draw(); h1max = Dexyc->GetMaximum()*2.05; Double_t xe[4] = {2000.,5000.,5000.,2000.}; Double_t ye[4] = {0.,0.,h1max,h1max}; banda1 = new TPolyLine(4,xe,ye); banda1->SetLineColor(5); banda1->SetFillColor(5); banda1->SetLineWidth(1); banda1->Draw("fSAME"); Double_t xj[4] = {0.,40.,40.,0.}; Double_t yj[4] = {0.,0.,h1max,h1max}; banda2 = new TPolyLine(4,xj,yj); banda2->SetLineColor(5); banda2->SetFillColor(5); banda2->SetLineWidth(1); banda2->Draw("fSAME"); Dexyc->Draw("SAME"); // pd3->cd(); // gPad->SetLogy(); Calstriphit->SetXTitle("Number of hit"); Calstriphit->SetYTitle("Number of events"); Calstriphit->Draw(); h1max = Calstriphit->GetMaximum()*1.05; Double_t xg[4] = {10.,100.,100.,10.}; Double_t yg[4] = {0.,0.,h1max,h1max}; banda1 = new TPolyLine(4,xg,yg); banda1->SetLineColor(5); banda1->SetFillColor(5); banda1->SetLineWidth(1); banda1->Draw("fSAME"); // Double_t xh[4] = {40.,80.,80.,40.}; // Double_t yh[4] = {0.,0.,h1max,h1max}; // banda2 = new TPolyLine(4,xh,yh); // banda2->SetLineColor(5); // banda2->SetFillColor(5); // banda2->SetLineWidth(1); // banda2->Draw("fSAME"); Calstriphit->Draw("SAME"); // pd4->cd(); Dexy->SetXTitle("ADC channels"); Dexy->SetYTitle("Number of events"); Dexy->Draw(); h1max = Dexy->GetMaximum()*1.05; Double_t xd[4] = {2000.,5000.,5000.,2000.}; Double_t yd[4] = {0.,0.,h1max,h1max}; banda1 = new TPolyLine(4,xd,yd); banda1->SetLineColor(5); banda1->SetFillColor(5); banda1->SetLineWidth(1); banda1->Draw("fSAME"); Dexy->Draw("SAME"); } else { Int_t pd4col = 10; if ( (outdexy+indexy) ) { cdexy = (float)indexy/((float)outdexy + (float)indexy); if ( cdexy < cdexythr ) { check = true; pd4col = 45; }; }; figura2->SetFillColor(10); figura2->Range(0,0,100,100); errore.str(""); errore << "EXPERT -- File: " << file; errore << " "; t=new TLatex(); t->SetTextFont(32); t->SetTextColor(1); t->SetTextAlign(12); t->SetTextSize(0.015); t->DrawLatex(2.,99.,errore.str().c_str()); pd4 = new TPad("pd4","This is pad4",0.02,0.02,0.98,0.98,pd4col); figura2->cd(); pd4->Range(0,0,100,100); pd4->SetTicks(); pd4->Draw(); // pd4->cd(); Dexy->SetXTitle("ADC channels"); Dexy->SetYTitle("Number of events"); Dexy->Draw(); h1max = Dexy->GetMaximum()*1.05; Double_t xd[4] = {2000.,5000.,5000.,2000.}; Double_t yd[4] = {0.,0.,h1max,h1max}; banda1 = new TPolyLine(4,xd,yd); banda1->SetLineColor(5); banda1->SetFillColor(5); banda1->SetLineWidth(1); banda1->Draw("fSAME"); Dexy->Draw("SAME"); }; // TCanvas *figura3 = 0; Bool_t printfigure3 = false; for (Int_t i = 0; i<4; i++){ if ( ver[i][2] || lver[i][2] ){ printfigure3 = true; break; }; }; if ( printfigure3 ){ figura3 = new TCanvas("Calorimeter_Detector_Report_2bis/3","Calorimeter_Detector_Report_2bis/3", 1100, 900); figura3->SetFillColor(10); figura3->Range(0,0,100,100); errore.str(""); errore << "EXPERT -- File: " << file; errore << " "; t=new TLatex(); t->SetTextFont(32); t->SetTextColor(1); t->SetTextAlign(12); t->SetTextSize(0.015); t->DrawLatex(2.,99.,errore.str().c_str()); pd1 = new TPad("pd1","This is pad1",0.02,0.51,0.49,0.73,45); pd5 = new TPad("pd5","This is pad5",0.02,0.76,0.49,0.98,45); pd2 = new TPad("pd2","This is pad2",0.51,0.51,0.98,0.73,45); pd6 = new TPad("pd6","This is pad6",0.51,0.76,0.98,0.98,45); pd3 = new TPad("pd3","This is pad3",0.02,0.02,0.49,0.23,45); pd7 = new TPad("pd7","This is pad7",0.02,0.26,0.49,0.49,45); pd4 = new TPad("pd4","This is pad4",0.51,0.02,0.98,0.23,45); pd8 = new TPad("pd8","This is pad8",0.51,0.26,0.98,0.49,45); figura3->cd(); pd1->Range(0,0,100,100); pd2->Range(0,0,100,100); pd3->Range(0,0,100,100); pd4->Range(0,0,100,100); pd1->SetTicks(); pd2->SetTicks(); pd3->SetTicks(); pd4->SetTicks(); pd1->Draw(); pd2->Draw(); pd3->Draw(); pd4->Draw(); pd5->Range(0,0,100,100); pd6->Range(0,0,100,100); pd7->Range(0,0,100,100); pd8->Range(0,0,100,100); pd5->SetTicks(); pd6->SetTicks(); pd7->SetTicks(); pd8->SetTicks(); pd5->Draw(); pd6->Draw(); pd7->Draw(); pd8->Draw(); pd1->cd(); lup[0]->SetAxisRange((Double_t)min(minobt[0],swminobt[0])*0.9,(Double_t)max(maxobt[0],swmaxobt[0])*1.1,"X"); lup[0]->SetXTitle("OBT"); lup[0]->SetYTitle("Number of events"); lup[0]->Draw(); pd5->cd(); lupstw[0]->SetAxisRange((Double_t)min(minobt[0],swminobt[0])*0.9,(Double_t)max(maxobt[0],swmaxobt[0])*1.1,"X"); // lupstw[0]->SetAxisRange((Double_t)swminobt[0]*0.9,(Double_t)swmaxobt[0]*1.1,"X"); lupstw[0]->SetXTitle("OBT"); lupstw[0]->SetYTitle("Number of events"); lupstw[0]->Draw(); pd2->cd(); lup[1]->SetAxisRange((Double_t)min(minobt[1],swminobt[1])*0.9,(Double_t)max(maxobt[1],swmaxobt[1])*1.1,"X"); // lup[1]->SetAxisRange((Double_t)minobt[1]*0.9,(Double_t)maxobt[1]*1.1,"X"); lup[1]->SetXTitle("OBT"); lup[1]->SetYTitle("Number of events"); lup[1]->Draw(); pd6->cd(); lupstw[1]->SetAxisRange((Double_t)min(minobt[1],swminobt[1])*0.9,(Double_t)max(maxobt[1],swmaxobt[1])*1.1,"X"); // lupstw[1]->SetAxisRange((Double_t)swminobt[1]*0.9,(Double_t)swmaxobt[1]*1.1,"X"); lupstw[1]->SetXTitle("OBT"); lupstw[1]->SetYTitle("Number of events"); lupstw[1]->Draw(); pd3->cd(); lup[2]->SetAxisRange((Double_t)min(minobt[2],swminobt[2])*0.9,(Double_t)max(maxobt[2],swmaxobt[2])*1.1,"X"); // lup[2]->SetAxisRange((Double_t)minobt[2]*0.9,(Double_t)maxobt[2]*1.1,"X"); lup[2]->SetXTitle("OBT"); lup[2]->SetYTitle("Number of events"); lup[2]->Draw(); pd7->cd(); lupstw[2]->SetAxisRange((Double_t)min(minobt[2],swminobt[2])*0.9,(Double_t)max(maxobt[2],swmaxobt[2])*1.1,"X"); // lupstw[2]->SetAxisRange((Double_t)swminobt[2]*0.9,(Double_t)swmaxobt[2]*1.1,"X"); lupstw[2]->SetXTitle("OBT"); lupstw[2]->SetYTitle("Number of events"); lupstw[2]->Draw(); pd4->cd(); lup[3]->SetAxisRange((Double_t)min(minobt[3],swminobt[3])*0.9,(Double_t)max(maxobt[3],swmaxobt[3])*1.1,"X"); // lup[3]->SetAxisRange((Double_t)minobt[3]*0.9,(Double_t)maxobt[3]*1.1,"X"); lup[3]->SetXTitle("OBT"); lup[3]->SetYTitle("Number of events"); lup[3]->Draw(); pd8->cd(); lupstw[3]->SetAxisRange((Double_t)min(minobt[3],swminobt[3])*0.9,(Double_t)max(maxobt[3],swmaxobt[3])*1.1,"X"); // lupstw[3]->SetAxisRange((Double_t)swminobt[3]*0.9,(Double_t)swmaxobt[3]*1.1,"X"); lupstw[3]->SetXTitle("OBT"); lupstw[3]->SetYTitle("Number of events"); lupstw[3]->Draw(); }; // // output figures, report sheet: // TCanvas *rapporto= new TCanvas("Calorimeter_Detector_Report_3/3", "Calorimeter_Detector_Report_3/3", 1100, 900); rapporto->SetFillColor(10); rapporto->Range(0,0,100,100); t=new TLatex(); t->SetTextFont(32); t->SetTextColor(1); t->SetTextSize(0.05); t->SetTextAlign(12); // t->DrawLatex(33.,95.,"Calorimeter quick look: "); errore.str(""); errore << "BASIC: C1 - File: " << file; errore << " "; t->SetTextSize(0.02); t->DrawLatex(2.,99.,errore.str().c_str()); pad1 = new TPad("pad1","This is pad1",0.02,0.47,0.49,0.90,19); pad2 = new TPad("pad2","This is pad2",0.51,0.47,0.98,0.90,19); pad3 = new TPad("pad3","This is pad3",0.02,0.02,0.49,0.45,19); pad4 = new TPad("pad4","This is pad4",0.51,0.02,0.98,0.45,19); pad1->Range(0,0,100,100); pad2->Range(0,0,100,100); pad3->Range(0,0,100,100); pad4->Range(0,0,100,100); pad1->Draw(); pad2->Draw(); pad3->Draw(); pad4->Draw(); for (Int_t i = 0; i < 4; i++){ if (i == 2) { pad1->cd() ; sezione = "** Section YE (x even) **"; } if (i == 3) { pad2->cd(); sezione = "** Section YO (x odd) **"; } if (i == 0) { pad3->cd(); sezione = "** Section XE (y odd) **"; } if (i == 1) { pad4->cd(); sezione = "** Section XO (y even) **"; } t->SetTextFont(32); t->SetTextColor(1); t->SetTextSize(0.05); t->SetTextAlign(12); t->DrawLatex(33.,97.,sezione); t->SetTextSize(0.05); for (Int_t j = 0; j < 23; j++){ if ( ver[i][j] || lver[i][j] ) { t->SetTextColor(50); if (j == 0) { errore.str(""); errore << "* DSP ack error: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,30.,errore.str().c_str()); check = true; } if (j == 1) { errore.str(""); errore << "* Temperature alarm: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,74.,errore.str().c_str()); check = true; } if (j == 2) { // printf("boh! lalarm[%i] = %i \n",i,lalarm[i]); if ( lalarm[i] ){ errore.str(""); errore << "* Latch up: " << ver[i][j]; errore << " (" << lalarm[i]; errore << " NOT rec!): "; t->DrawLatex(2.,70.,errore.str().c_str()); check = true; } else { errore.str(""); errore << "* Latch up alarm: " << ver[i][j]; errore << " time(s). View(s): "; t->DrawLatex(2.,70.,errore.str().c_str()); check = true; }; errore.str(""); errore << "(" << lver[i][j]; errore << " view(s))"; t->DrawLatex(27.,65.,errore.str().c_str()); Float_t inc=1.; for (Int_t e = 0; e < 11; e++){ if ( pshit[i][e] ) { errore.str(""); errore << " " << (e+1); errore << " "; t->DrawLatex(52.+inc,69.,errore.str().c_str()); Int_t numer = pshit[i][e]; errore.str(""); errore << " (" << numer; errore << ") "; t->SetTextAngle(-70); t->DrawLatex(55.+inc,67.,errore.str().c_str()); t->SetTextAngle(0); inc += 4.; } } } if (j == 3) { errore.str(""); errore << "RAW mode: " << ver[i][j]; errore << " time(s)"; t->SetTextColor(32); t->DrawLatex(2.,90.,errore.str().c_str()); } if (j == 4) { errore.str(""); errore << "* CMD length error: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,66.,errore.str().c_str()); check = true; } if (j == 5) { errore.str(""); errore << "* Execution error: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,62.,errore.str().c_str()); check = true; } if (j == 6) { errore.str(""); errore << "* CRC error (st. word): " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,58.,errore.str().c_str()); check = true; } if (j == 7) { errore.str(""); errore << "View or command not recognized: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,54.,errore.str().c_str()); check = true; } if (j == 8) { errore.str(""); errore << "Missing section: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,50.,errore.str().c_str()); if ( ver[i][j] > 3 ) check = true; } if (j == 10) { errore.str(""); errore << "Calevnum jump: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,42.,errore.str().c_str()); check = true; } if (j == 11) { errore.str(""); errore << "CRC error (data): " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,38.,errore.str().c_str()); if ( ver[i][j] > 10 ) check = true; } if (j == 12) { errore.str(""); errore << "Length problems in RAW mode: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,34.,errore.str().c_str()); check = true; } if (j == 13) { errore.str(""); errore << "Length problems in COMPRESS mode: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,34.,errore.str().c_str()); check = true; } if (j == 14) { errore.str(""); errore << "Length problems in FULL mode: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,26.,errore.str().c_str()); check = true; } if (j == 15) { errore.str(""); errore << "Acquisition mode problems: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,22.,errore.str().c_str()); check = true; } if (j == 16) { errore.str(""); errore << "COMPRESS mode: " << ver[i][j]; errore << " time(s)"; t->SetTextColor(32); t->DrawLatex(2.,86.,errore.str().c_str()); } if (j == 17) { errore.str(""); errore << "FULL mode: " << ver[i][j]; errore << " time(s)"; t->SetTextColor(32); t->DrawLatex(2.,82.,errore.str().c_str()); } if (j == 18) { errore.str(""); errore << "Problems with coding: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,18.,errore.str().c_str()); check = true; } if (j == 19) { errore.str(""); errore << "Pedestal checksum wrong: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,14.,errore.str().c_str()); check = true; } if (j == 20) { errore.str(""); errore << "Thresholds checksum wrong: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,10.,errore.str().c_str()); check = true; } if (j == 21) { errore.str(""); errore << "WARNING! DEXYC < 0: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,6.,errore.str().c_str()); check = true; } if (j == 22) { errore.str(""); errore << "Packet length is zero (YODA input error), skipped: " << ver[i][j]; errore << " time(s)"; t->DrawLatex(2.,3.,errore.str().c_str()); check = true; }; }; }; t->SetTextColor(50); if ( fulldiff[i] !=0 ) { check = true; errore.str(""); errore << "Full mode, differences between RAW and COMPRESS mode: " << fulldiff[i]; errore << " time(s)"; t->DrawLatex(2.,46.,errore.str().c_str()); }; if ( cpre[i] !=0 ) { errore.str(""); errore << "Preamplifier fully transmitted: " << cpre[i]; errore << " time(s)"; t->SetTextColor(32); t->DrawLatex(2.,78.,errore.str().c_str()); }; }; rapporto->cd(); t->SetTextFont(32); t->SetTextColor(1); t->SetTextSize(0.05); t->SetTextAlign(12); t->DrawLatex(22.,95.,"Calorimeter quick look: "); //printf("check %i \n",check); if ( check ){ t->SetTextColor(50); t->DrawLatex(60.,95.," WARNING, CHECK! "); printf("cdexyc %f cdexy %f ctshit %f cbase %f \n",cdexyc,cdexy,ctshit,cbase); } else { t->SetTextColor(32); t->DrawLatex(60.,95.," OK! "); }; // // const string fil = (const char*)filename; // Int_t posiz = fil.find("dw_"); // if ( posiz == -1 ) posiz = fil.find("DW_"); // Int_t posiz2 = posiz+13; // TString file2; // stringcopy(file2,filename,posiz,posiz2); // const string fil = gSystem->BaseName(filename.Data()); Int_t posiz = fil.find(".root"); // TString file2; if ( posiz == -1 ){ file2 = gSystem->BaseName(filename.Data()); } else { Int_t posiz2 = 0; stringcopy(file2,gSystem->BaseName(filename.Data()),posiz2,posiz); }; const char *figrec = file2; // const char *outdir = outDir; stringstream figsave; stringstream figsave1; stringstream figsave2; const char *format = saveas; if ( !strcmp(format,"ps") ) { figsave.str(""); figsave << outdir << "/" ; figsave << figrec << "_CaloQLOOK."; figsave << format << "("; rapporto->Print(figsave.str().c_str(),"Landscape"); figsave1.str(""); figsave1 << outdir << "/" ; figsave1 << figrec << "_CaloQLOOK."; figsave1 << format; figura->Print(figsave1.str().c_str(),"Landscape"); figsave2.str(""); figsave2 << outdir << "/" ; figsave2 << figrec << "_CaloQLOOK."; if ( printfigure3 ){ figsave2 << format; figura2->Print(figsave2.str().c_str(),"Landscape"); figsave2.str(""); figsave2 << outdir << "/" ; figsave2 << figrec << "_CaloQLOOK."; figsave2 << format << ")"; figura3->Print(figsave2.str().c_str(),"Landscape"); } else { figsave2 << format << ")";; figura2->Print(figsave2.str().c_str(),"Landscape"); }; } else { figsave.str(""); figsave << outdir << "/" ; figsave << figrec << "_CaloQLOOK1."; figsave << format; figura->SaveAs(figsave.str().c_str()); figsave.str(""); figsave << outdir << "/" ; figsave << figrec << "_CaloQLOOK2."; figsave << format; figura2->SaveAs(figsave.str().c_str()); if ( printfigure3 ){ figsave.str(""); figsave << outdir << "/" ; figsave << figrec << "_CaloQLOOK2bis."; figsave << format; figura3->SaveAs(figsave.str().c_str()); }; figsave.str(""); figsave << outdir << "/" ; figsave << figrec << "_CaloQLOOK3."; figsave << format; rapporto->SaveAs(figsave.str().c_str()); }; }