FQLOOK/macros/FCaloQLOOK.cxx

//
//   Check for possible errors and shows raw distribution of variables coming from the DSP (no calibrations needed) - Emiliano Mocchiutti
//
//   FCaloQLOOK.c      version 1.13  (2006-09-28)
//
//   The only input needed is the path to the directory created by YODA for the data file you want to analyze.
//
//   Changelog:
//
//   1.12 - 1.13 (2006-09-28): Give wrong number of events in compress mode, fixed.
//
//   1.11 - 1.12 (2006-09-22): Check for different number of events, put output in a file and write down errors for each event. Fixed jumpev bug.
//
//   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 <TApplication.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 sighandler(UInt_t pkt_num, UInt_t obt, UInt_t entry, UInt_t section, UInt_t totev, Bool_t &evdone){
  if ( !evdone ){    
    printf("\n\n=================================================\n");
    printf("PSCU-Pkt N. %u - OBT %u ms\n",pkt_num,obt);
    printf("Total events %u\n",totev);
    printf("(ROOT-tree entry %u )\n",entry);
    printf("=================================================\n");
    evdone = true;
  };
  TString sec = "UNKNOWN";
  if ( section == 0 ) sec = "XE";
  if ( section == 1 ) sec = "XO";
  if ( section == 2 ) sec = "YE";
  if ( section == 3 ) sec = "YO";
  printf(" Section %s (DSP n. %i) message: ",sec.Data(),section);
};

void FCaloQLOOK(TString filename, Int_t fromevent=0, Int_t toevent=0, TString outDir="", TString saveas = "png", Bool_t iactive =false, Bool_t w4i=false){
  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;
  //
  TApplication *app = 0;
  if ( iactive ) app = new TApplication("app",0,0);
  //
  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 = false;
  bool isFULL = false;
  bool isRAW = false;
  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");
  printf("\n Start processing: \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;
  UInt_t headc = 0;
  UInt_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;
  UInt_t obt = 0;
  UInt_t minobt[4];
  UInt_t maxobt[4];
  UInt_t swminobt[4];
  UInt_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;
  };
  Bool_t evdone = false;
  //
  while ( i < maxevent+1){
    evdone = false;
    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++;
    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]++;
          sighandler(headc,obt,i,0,nevents,evdone);
          printf(" jump in the counter calev is %i oldcalev was %i \n",calev0,oldcalev0);
        };
        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]++;
          sighandler(headc,obt,i,1,nevents,evdone);
          printf(" jump in the counter calev is %i oldcalev was %i \n",calev1,oldcalev1);
        };
        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]++;
          sighandler(headc,obt,i,2,nevents,evdone);
          printf(" jump in the counter calev is %i oldcalev was %i \n",calevv2,oldcalev2);
        };
        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]++;
          sighandler(headc,obt,i,3,nevents,evdone);
          printf(" jump in the counter calev is %i oldcalev was %i \n",calev3,oldcalev3);
        };
      }
    };
    //  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 = false;
        isFULL = false;
        isRAW = false;
        if ( ce->stwerr[se] & (1 << 16) ) isCOMP = true; 
        if ( ce->stwerr[se] & (1 << 17) ) isFULL = true; 
        if ( ce->stwerr[se] & (1 << 3) ) isRAW = true; 
        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]++ ;         
                sighandler(headc,obt,i,se,nevents,evdone);
                printf(" latch up (data) in DSP mode alldexy %i planebases %i \n",alldexy,planebases);
              };
            }
            //
            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]++ ;
              sighandler(headc,obt,i,se,nevents,evdone);
              printf(" negative ADC values!! ce->dexyc[%i][%i][%i] = %f \n",l,ii,kk,ce->dexyc[l][ii][kk]);
            };
            //
            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]++ ;
                sighandler(headc,obt,i,se,nevents,evdone);
                printf(" latch up (data) in RAW mode alldexy2 %i \n",alldexy2);
              };
            };
            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 = false;
      isFULL = false;
      cestw=0;
      if ( ce->stwerr[k] ) cestw =  ce->stwerr[k] & cmask ;
      if ( ce->stwerr[k] & (1 << 16) && !(cestw & (1 << 3)) ) isCOMP = true; 
      if ( ce->stwerr[k] & (1 << 17) && !(cestw & (1 << 3)) ) isFULL = true; 
      //
      if ( isCOMP ) ver[k][16]++;
      if ( isFULL ) ver[k][17]++;
      cshit += (int)ce->calstriphit[k];
      if ( cestw ){
        if ( cestw & (1 << 0) ){
          ver[k][6]++ ;
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" CRC error (st. word) \n");
        };
        if ( cestw & (1 << 1) ){
          ver[k][5]++ ;
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" execution error \n");
        };
        if ( cestw & (1 << 2) ){
          ver[k][4]++ ;
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" CMD length error \n");
        };
        if ( cestw & (1 << 3) ) ver[k][3]++ ; // raw mode
        if ( cestw & (1 << 4) ){
          ver[k][2]++ ;
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" latch up alarm (st. word)\n");
          lupstw[k]->Fill(obt);
          if ( swfirstobt[k] ) swminobt[k] = obt;
          if ( obt > swmaxobt[k] ) swmaxobt[k] = obt;
        };
        if ( cestw & (1 << 5) ){
          ver[k][1]++ ;
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" temperature alarm  \n");
        };
        if ( cestw & (1 << 6) ){
          ver[k][0]++ ;         
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" DSP ack error \n");
        };
      };
      if (  ce->stwerr[k] == 0 && ce->perror[k] == 0 ) ver[k][8]++ ; 
      if ( ce->perror[k] != 0. ){
        if (ce->perror[k] == 128){
          ver[k][7]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" view or command not recognized \n");
        };
        if (ce->perror[k] == 129){
          ver[k][8]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" missing section \n");
        };
        if (ce->perror[k] == 132){
          ver[k][11]++ ; 
          if ( ver[k][11] < 30 ){
            sighandler(headc,obt,i,k,nevents,evdone);
            printf(" CRC error (data) \n");
          };
          if ( ver[k][11] == 30 ){
            sighandler(headc,obt,i,k,nevents,evdone);
            printf(" CRC error (data) \n");
            printf(" WARNING TOO MANY CRC ERRORS ON DATA FOR THIS SECTION \n");
            printf(" THEY WILL NOT BE DISPLAYED HERE ANYMORE FOR THIS FILE! \n");
          };
        };
        if (ce->perror[k] == 133){
          ver[k][12]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" length problems in RAW mode \n");
        };
        if (ce->perror[k] == 134){
          ver[k][13]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" length problems in COMPRESS mode \n");
        };
        if (ce->perror[k] == 135){
          ver[k][14]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" length problems in FULL mode \n");
        };
        if (ce->perror[k] == 136){
          ver[k][15]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" acquisition mode problems \n");
        };
        if (ce->perror[k] == 139){
          ver[k][18]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" problems with coding\n");
        };
        if (ce->perror[k] == 140){
          ver[k][19]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" pedestal checksum wrong\n");
        };
        if (ce->perror[k] == 141){
          ver[k][20]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" thresholds checksum wrong\n");
        };
        if (ce->perror[k] == 142){
          ver[k][22]++ ; 
          sighandler(headc,obt,i,k,nevents,evdone);
          printf(" packet length is zero (YODA input error), skipped\n");
        };
      };
      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");
  printf("\n end processing. \n\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();
    gPad->SetLogy();
    Dexy->SetXTitle("ADC channels");
    Dexy->SetYTitle("Number of events");
    Dexy->Draw();
    //    h1max = Dexy->GetMaximum()*1.05;
    h1max = Dexy->GetMaximum()*2.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();
    gPad->SetLogy();
    Dexy->SetXTitle("ADC channels");
    Dexy->SetYTitle("Number of events");
    Dexy->Draw();
    //    h1max = Dexy->GetMaximum()*1.05;
    h1max = Dexy->GetMaximum()*2.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();
  //
  Bool_t redevent = true;
  //
  if ( ver[0][3]+ver[0][16]+ver[0][17] == ver[1][3]+ver[1][16]+ver[1][17] && ver[1][3]+ver[1][16]+ver[1][17] == ver[2][3]+ver[2][16]+ver[2][17] && ver[2][3]+ver[2][16]+ver[2][17] == ver[3][3]+ver[3][16]+ver[3][17] ){
    redevent = false; 
  };
  if ( redevent ) check = true;
  //
  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)";
          if ( !redevent ) 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());
          if ( ver[i][j] > 3 ) 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)";
          if ( !redevent ) 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)";
          if ( !redevent ) 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());    
  };
  if ( iactive && w4i ){
    while ( gROOT->GetListOfCanvases()->FindObject(rapporto) || gROOT->GetListOfCanvases()->FindObject(figura3) || gROOT->GetListOfCanvases()->FindObject(figura2) || gROOT->GetListOfCanvases()->FindObject(figura) ){
      gSystem->ProcessEvents();
      gSystem->Sleep(10);
    };
  };
}
