TH2Poly *Display0()
{
 TH2Poly *Plano0=new TH2Poly(); Plano0->SetTitle("Plano 0");	
 int nchannels=64; int ChList[64],ColList[64],LinList[64];	
 double top=50., bot=-50, left=-50, right=50.;
 double guard_ring=3./10.;//cm
 double borda=1.5/10.;//cm
 double strip=9.6875/10.;//cm
 double height=998./10.;//cm
 //strip+=2.*(borda+guard_ring/2.);//x64=1004.8
 strip=100./64.; //erro de 0.4%... não deve fazer diferença
 for(int ip=0;ip<64;ip++)
   {
     ChList[ip]=ip;
     Plano0->AddBin(left+ip*strip,bot,left+ip*strip+strip,top);
   }
 return Plano0;
}
TH2Poly *DisplayN()
{
 TH2Poly *strangeP=new TH2Poly();
 int nchannels=64; int ChList[64],ColList[64],LinList[64];	
 for(int ip=0;ip<64;ip++)
 {
   ChList[ip]=ip;
    
   if(ChList[ip]==15){ColList[ip]=1;LinList[ip]=1;}
   else if(ChList[ip]==3 ){ColList[ip]=1;LinList[ip]=2;}
   else if(ChList[ip]==16){ColList[ip]=1;LinList[ip]=3;}
   else if(ChList[ip]==4 ){ColList[ip]=1;LinList[ip]=4;}
   else if(ChList[ip]==17){ColList[ip]=1;LinList[ip]=5;}
   else if(ChList[ip]==5 ){ColList[ip]=1;LinList[ip]=6;}
   else if(ChList[ip]==18){ColList[ip]=1;LinList[ip]=7;}
   //similar to 61
   //else if(ChList[ip]==6){ColList[ip]=1;LinList[ip]=7;}
   
   else if(ChList[ip]==19){ColList[ip]=2;LinList[ip]=1;}
   else if(ChList[ip]==7 ){ColList[ip]=2;LinList[ip]=2;}
   else if(ChList[ip]==20){ColList[ip]=2;LinList[ip]=3;}
   else if(ChList[ip]==8 ){ColList[ip]=2;LinList[ip]=4;}
   else if(ChList[ip]==21){ColList[ip]=2;LinList[ip]=5;}
   else if(ChList[ip]==9 ){ColList[ip]=2;LinList[ip]=6;}
   else if(ChList[ip]==22){ColList[ip]=2;LinList[ip]=7;}
   
   else if(ChList[ip]==10){ColList[ip]=3;LinList[ip]=1;}
   else if(ChList[ip]==23){ColList[ip]=3;LinList[ip]=2;}
   else if(ChList[ip]==11){ColList[ip]=3;LinList[ip]=3;}
   else if(ChList[ip]==24){ColList[ip]=3;LinList[ip]=4;}
   else if(ChList[ip]==25){ColList[ip]=3;LinList[ip]=5;}
   else if(ChList[ip]==26){ColList[ip]=3;LinList[ip]=6;}
   else if(ChList[ip]==27){ColList[ip]=3;LinList[ip]=7;}
   
   else if(ChList[ip]==28){ColList[ip]=4;LinList[ip]=1;}
   else if(ChList[ip]==29){ColList[ip]=4;LinList[ip]=2;}
   else if(ChList[ip]==30){ColList[ip]=4;LinList[ip]=3;}
   else if(ChList[ip]==31){ColList[ip]=4;LinList[ip]=4;}
   else if(ChList[ip]==32){ColList[ip]=4;LinList[ip]=5;}
   else if(ChList[ip]==33){ColList[ip]=4;LinList[ip]=6;}
   else if(ChList[ip]==34){ColList[ip]=4;LinList[ip]=7;}
   
   else if(ChList[ip]==35){ColList[ip]=5;LinList[ip]=1;}
   else if(ChList[ip]==36){ColList[ip]=5;LinList[ip]=2;}
   else if(ChList[ip]==37){ColList[ip]=5;LinList[ip]=3;}
   else if(ChList[ip]==38){ColList[ip]=5;LinList[ip]=4;}
   else if(ChList[ip]==39){ColList[ip]=5;LinList[ip]=5;}
   else if(ChList[ip]==52){ColList[ip]=5;LinList[ip]=6;}
   else if(ChList[ip]==40){ColList[ip]=5;LinList[ip]=7;}
   
   else if(ChList[ip]==53){ColList[ip]=6;LinList[ip]=1;}
   else if(ChList[ip]==41){ColList[ip]=6;LinList[ip]=2;}
   else if(ChList[ip]==54){ColList[ip]=6;LinList[ip]=3;}
   else if(ChList[ip]==42){ColList[ip]=6;LinList[ip]=4;}
   else if(ChList[ip]==55){ColList[ip]=6;LinList[ip]=5;}
   else if(ChList[ip]==43){ColList[ip]=6;LinList[ip]=6;}
   else if(ChList[ip]==56){ColList[ip]=6;LinList[ip]=7;}
   
   
   else if(ChList[ip]==44){ColList[ip]=7;LinList[ip]=1;}
   else if(ChList[ip]==57){ColList[ip]=7;LinList[ip]=2;}
   else if(ChList[ip]==45){ColList[ip]=7;LinList[ip]=3;}
   else if(ChList[ip]==58){ColList[ip]=7;LinList[ip]=4;}
   else if(ChList[ip]==46){ColList[ip]=7;LinList[ip]=5;}
   else if(ChList[ip]==59){ColList[ip]=7;LinList[ip]=6;}
   else if(ChList[ip]==47){ColList[ip]=7;LinList[ip]=7;}
   
   ColList[ip]-=4; LinList[ip]-=4;
   
   //	exceptions 
   if(ChList[ip]==60){ColList[ip]=0;  LinList[ip]=-10;}
   else if(ChList[ip]==6 ){ColList[ip]=0;  LinList[ip]=+10;}
   else if(ChList[ip]==12){ColList[ip]=-19;LinList[ip]=0;}
   else if(ChList[ip]==0 ){ColList[ip]=-14;LinList[ip]=0;}
   else if(ChList[ip]==13){ColList[ip]=-13;LinList[ip]=0;}
   else if(ChList[ip]==1 ){ColList[ip]=-12;LinList[ip]=0;}
   else if(ChList[ip]==14){ColList[ip]=-11;LinList[ip]=0;}
   else if(ChList[ip]==2 ){ColList[ip]=-10;LinList[ip]=0;}
   else if(ChList[ip]==62){ColList[ip]=+10;LinList[ip]=0;}
   else if(ChList[ip]==48){ColList[ip]=+10;LinList[ip]=-5;}
   else if(ChList[ip]==49){ColList[ip]=+10;LinList[ip]=+5;}
   else if(ChList[ip]==50){ColList[ip]=+11;LinList[ip]=-7;}
   else if(ChList[ip]==61){ColList[ip]=+11;LinList[ip]=+7;}
   else if(ChList[ip]==63){ColList[ip]=+14;LinList[ip]=0;}
   else if(ChList[ip]==51){ColList[ip]=+19;LinList[ip]=0;}
   
   LinList[ip]=-LinList[ip];
   /*
     Seguir instruções do Luis:
   */
     double guard_ring=3./10.;//cm
     double borda=1./10.;//cm
     double square=38./10.;//cm *7
     double wide_left=261./10.;//cm *1
     double strip_left=12./10.;//cm *5
     double strip_right=12./10.;//cm *2
     double strip_wide=40./10.;//cm *1
     double wide_right=269./10.;//cm *1

     wide_left+=(borda+guard_ring/2.);
     square+=2*(guard_ring/2.+borda);
     strip_left+=2*(guard_ring/2.+borda);
     strip_right+=2*(guard_ring/2.+borda);
     strip_wide+=2*(guard_ring/2.+borda);
     wide_right+=(borda+guard_ring/2.);
     /*
     double tot4=wide_left+5*strip_left+7*square+2*strip_right+strip_wide+wide_right;
     double center1=wide_left+5*strip_left+3.5*square;
     double center3=3.5*square+2*strip_right+strip_wide+wide_right;
     */
     double center=0.;//(center3-center1)/2.;
     double left=-50.;
     double right=50.;
     double top=50.;
     double bot=-50.;
     
     if(abs(LinList[ip])<4 && abs(ColList[ip])<4.)
       strangeP->AddBin((center+ColList[ip]*square)-square/2.,LinList[ip]*square-square/2.,
			(center+ColList[ip]*square)+square/2.,LinList[ip]*square+square/2.);
     
     else if(LinList[ip]==10)
       strangeP->AddBin((center-3.*square)-square/2.,3.*square+square/2.,(center+3.*square)+square/2.,top);	
     else if(LinList[ip]==-10)
       strangeP->AddBin((center-3.*square)-square/2.,-(3.*square+square/2.),(center+3.*square)+square/2.,bot);	

     else if(ChList[ip]==12)
       strangeP->AddBin(left,bot,left+wide_left,top);

     else if(ChList[ip]==0)
       strangeP->AddBin(left+wide_left,bot,left+wide_left+strip_left,top);		
     else if(ChList[ip]==13)
       strangeP->AddBin(left+wide_left+strip_left,bot,left+wide_left+2*strip_left,top);
     else if(ChList[ip]==1)
       strangeP->AddBin(left+wide_left+2*strip_left,bot,left+wide_left+3*strip_left,top);
     else if(ChList[ip]==14)
       strangeP->AddBin(left+wide_left+3*strip_left,bot,left+wide_left+4*strip_left,top);
     else if(ChList[ip]==2)
       strangeP->AddBin(left+wide_left+4*strip_left,bot,left+wide_left+5*strip_left,top);

     else if(ChList[ip]==51)
       strangeP->AddBin(right-wide_right,bot,right,top);
     else if(ChList[ip]==63)
       strangeP->AddBin(right-wide_right-strip_wide,bot,right-wide_right,top);
     else if(ChList[ip]==50)
       strangeP->AddBin(right-wide_right-strip_wide-strip_right,bot,right-wide_right-strip_wide,0);	
     else if(ChList[ip]==61)
       strangeP->AddBin(right-wide_right-strip_wide-strip_right,top,right-wide_right-strip_wide,0);	
   
     else if(ChList[ip]==62)
       strangeP->AddBin(right-wide_right-strip_wide-2*strip_right,top-(top-bot)/3.,
			right-wide_right-strip_wide-strip_right,bot+(top-bot)/3.);	
     else if(ChList[ip]==48)
       strangeP->AddBin(right-wide_right-strip_wide-2*strip_right,top-(top-bot)/3.,
			right-wide_right-strip_wide-strip_right,top);	
     else if(ChList[ip]==49)
       strangeP->AddBin(right-wide_right-strip_wide-2*strip_right,bot+(top-bot)/3.,
			right-wide_right-strip_wide-strip_right,bot);	
 }
 return strangeP;
}
void EventDisplayOnline(int doAll=1, int isave=1)
{
	TFile *f=new TFile("MuTomNow.root");
	TTree *T=(TTree*)f->Get("RpcTree");
	
	//f->ls();

//---  TChain *T=new TChain("RpcTree","RpcTree");
//---  T->Add("MuTomNow.root");



  // T->Add("~/DataCoimbra/RootDay/Day*root");
//  T->Add("Dayly/MuTomShort2t_26Nov21.root");
//  T->Add("Dayly/MuTomShort2t_19-26Nov2021.root");

  gStyle->SetOptStat(0);  gStyle->SetPalette(kCherry); TColor::InvertPalette();
  
  TH2F *DIRMAP=new TH2F("DIRMAP","Dir. (todos os mu#tilde{o}es); O - E; S - N",31,-16.,16.,31,-16.,16.);
  DIRMAP->GetXaxis()->CenterTitle(); DIRMAP->GetYaxis()->CenterTitle();
  
  TH2F *Vertical=new TH2F("Vertical","CorePix Vertical Counts; Coluna; Linha",7,-3.5,3.5,7,-3.5,3.5);
  Vertical->GetXaxis()->CenterTitle(); Vertical->GetYaxis()->CenterTitle();
  TH2F *Direction=new TH2F("Direction","CorePix Direction Map; d Coluna; d Linha",13,-6.5,6.5,13,-6.5,6.5);
  Direction->GetXaxis()->CenterTitle(); Direction->GetYaxis()->CenterTitle();
  
  TH2F *EVertical=new TH2F("Vertical Corr","CorePix Vertical Eff. Counts; Coluna; Linha",7,-3.5,3.5,7,-3.5,3.5);
  EVertical->GetXaxis()->CenterTitle(); EVertical->GetYaxis()->CenterTitle();
  TH2F *EDirection=new TH2F("Direction Corr","CorePix Direction Map; d Coluna; d Linha",13,-6.5,6.5,13,-6.5,6.5);
  EDirection->GetXaxis()->CenterTitle(); EDirection->GetYaxis()->CenterTitle();
  

  TH1F *hMtrig=new TH1F("hMtrig","hMtrig",2000,0,2000);	
	    Direction->GetXaxis()->SetLabelSize(0.05); //Direction->GetXaxis()->SetLabelOffset(-0.02);
	    Direction->GetXaxis()->SetTitleSize(0.05); //Direction->GetXaxis()->SetTitleOffset(0.65);

	    Direction->GetYaxis()->SetLabelSize(0.05); //Direction->GetYaxis()->SetLabelOffset(0.01);
	    Direction->GetYaxis()->SetTitleSize(0.05); //Direction->GetYaxis()->SetTitleOffset(0.70);

  
	    Vertical->GetXaxis()->SetLabelSize(0.05); //Vertical->GetXaxis()->SetLabelOffset(-0.02);
	    Vertical->GetXaxis()->SetTitleSize(0.05); //Vertical->GetXaxis()->SetTitleOffset(0.65);

	    Vertical->GetYaxis()->SetLabelSize(0.05); //Vertical->GetYaxis()->SetLabelOffset(0.01);
	    Vertical->GetYaxis()->SetTitleSize(0.05); //Vertical->GetYaxis()->SetTitleOffset(0.70);


	    EDirection->GetXaxis()->SetLabelSize(0.05); //Direction->GetXaxis()->SetLabelOffset(-0.02);
	    EDirection->GetXaxis()->SetTitleSize(0.05); //Direction->GetXaxis()->SetTitleOffset(0.65);

	    EDirection->GetYaxis()->SetLabelSize(0.05); //Direction->GetYaxis()->SetLabelOffset(0.01);
	    EDirection->GetYaxis()->SetTitleSize(0.05); //Direction->GetYaxis()->SetTitleOffset(0.70);

  
	    EVertical->GetXaxis()->SetLabelSize(0.05); //Vertical->GetXaxis()->SetLabelOffset(-0.02);
	    EVertical->GetXaxis()->SetTitleSize(0.05); //Vertical->GetXaxis()->SetTitleOffset(0.65);

	    EVertical->GetYaxis()->SetLabelSize(0.05); //Vertical->GetYaxis()->SetLabelOffset(0.01);
	    EVertical->GetYaxis()->SetTitleSize(0.05); //Vertical->GetYaxis()->SetTitleOffset(0.70);





  TH2F *DIRMAP0=new TH2F("DIRMAP0","DIRMAP0",31,-16.,16.,31,-16.,16.);
  TH2Poly *Plano0; Plano0=(TH2Poly*)Display0()->Clone("Plano 0"); Plano0->SetTitle(";Plano 0;");
  TH2Poly *strangeP; strangeP=(TH2Poly*)DisplayN()->Clone("strangeP");
  int nchannels=64; int ChList[64],ColList[64],LinList[64];	

  //  strangeP->SetTitleSize(20);
  strangeP->GetXaxis()->SetNdivisions(1); strangeP->GetYaxis()->SetNdivisions(0);
  strangeP->GetXaxis()->SetLabelSize(0);
  strangeP->GetXaxis()->SetTitleSize(0.1);
  strangeP->GetXaxis()->SetTitleOffset(0.275);
  strangeP->GetXaxis()->CenterTitle(1);
  Plano0->GetXaxis()->SetNdivisions(1);  Plano0->GetYaxis()->SetNdivisions(0);
  Plano0->GetXaxis()->SetLabelSize(0);
  Plano0->GetXaxis()->SetTitleSize(0.1);
  Plano0->GetXaxis()->SetTitleOffset(0.275);
  Plano0->GetXaxis()->CenterTitle(1);
  
  //  strangeP->SetOptTitle(0); Plano0->SetOptTitle(0);
  
  TH2Poly *Plano1=(TH2Poly*)strangeP->Clone("Plano1"); Plano1->SetTitle(";Plano -10;");	
  TH2Poly *Plano2=(TH2Poly*)strangeP->Clone("Plano2"); Plano2->SetTitle(";Plano -20;");		
  TH2Poly *Plano3=(TH2Poly*)strangeP->Clone("Plano3"); Plano3->SetTitle(";Plano -30");		
  TH2Poly *P0Empty=(TH2Poly*)Plano0->Clone("P0Empty"); 	

  TH2Poly *DisPl0=(TH2Poly*)Plano0->Clone("Plano0"); Plano0->SetTitle("Plano 0"); 
  TH2Poly *DisPl1=(TH2Poly*)Plano1->Clone("Plano1"); Plano1->SetTitle("Plano 1");	
  TH2Poly *DisPl2=(TH2Poly*)Plano2->Clone("Plano2"); Plano2->SetTitle("Plano 2");	
  TH2Poly *DisPl3=(TH2Poly*)Plano3->Clone("Plano3"); Plano3->SetTitle("Plano 3");	
 
  TH1I *channel=new TH1I("channel","channel",64,0.,64.);
  TH1I *ch0=(TH1I*)channel->Clone("ch0");	
  TH1I *ch1=(TH1I*)channel->Clone("ch1");	
  TH1I *ch2=(TH1I*)channel->Clone("ch2");	
  TH1I *ch3=(TH1I*)channel->Clone("ch3");	

  TH1I *ch0T=(TH1I*)channel->Clone("ch0T");	
  TH1I *ch1T=(TH1I*)channel->Clone("ch1T");	
  TH1I *ch2T=(TH1I*)channel->Clone("ch2T");	
  TH1I *ch3T=(TH1I*)channel->Clone("ch3T");	
  
   TH3F *Disp3D=new TH3F("Disp3D",
			 //			"#acute{u}ltimo mu#tilde{a}o; Coluna; Linha; Plano",
			 //			"Dire#cedil{c}#tilde{a}o do #acute{u}ltimo mu#tilde{a}o; Coluna (Oeste - Este); Linha (Sul - Norte); Plano (Altura)",
		       " ; Oeste - Este ( Coluna ); Sul - Norte ( Linha ); Altura ( Plano )",
			64,-50.,50.,64,-50.,50.,40,-35.,5.);//25,-6.,6.,25,-6.,6.,40,-35.,5.);
   Disp3D->GetXaxis()->SetNdivisions(9);
   Disp3D->GetYaxis()->SetNdivisions(9);
   Disp3D->GetZaxis()->SetNdivisions(9);
   Disp3D->GetXaxis()->CenterTitle(1);
   Disp3D->GetYaxis()->CenterTitle(1);
   //   Disp3D->GetZaxis()->CenterTitle(1);
   //total -50;50 cm; com o centro a -15;15 cm (nas coluna é -3,--,0,--,3)
   float binwidth=100./64.;
   float quadr=4.3;
   
   Disp3D->SetLineColor(kRed); Disp3D->SetFillColor(kRed);
   Disp3D->GetXaxis()->SetTitleOffset(1.5); 
   Disp3D->GetYaxis()->SetTitleOffset(1.5); 
   Disp3D->GetZaxis()->SetTitleOffset(1.5); 
  
   TCanvas *Disp=new TCanvas("Display","Display",1500.,1000.);
   //TCanvas *Disp2=new TCanvas();   


	//ler eficiências dos histogramas
	//e corrigir no mapa vertical e direcoes
	TFile *Eff=new TFile("refEfficiency.root");
	TH2F *xE1=(TH2F*)Eff->Get("E01_SS");
	TH2F *xE2=(TH2F*)Eff->Get("E02_SS");
	TH2F *xE3=(TH2F*)Eff->Get("E03_SS");
	TH2F *xGC=(TH2F*)Eff->Get("GC");


   UInt_t MTrig, Data, first_Data=0, last_Data=0;
	int last_min=0, num_min=0;
   float Hora, first_Hora=0, last_Hora=0;
  int ChMax[4], Linha[4], Coluna[4];//, CColuna[4][64];
  //float Carga[4][64];
  T->SetBranchAddress("Data",&Data);
  T->SetBranchAddress("Hora",&Hora);
  T->SetBranchAddress("ChMax",&ChMax);
  T->SetBranchAddress("Linha",&Linha);
  T->SetBranchAddress("Coluna",&Coluna);
  T->SetBranchAddress("MTrig",&MTrig);
  //T->SetBranchAddress("CColuna",&CColuna);
  //T->SetBranchAddress("Carga",&Carga);

  int tDisplay=0;	
  for(int imu=0;imu<T->GetEntries();imu++)
//  for(int imu=0;imu<100;imu++)
    {
      Disp3D->Reset();    
      last_Hora=Hora;
      last_Data=Data;
      T->GetEntry(imu);

	hMtrig->Fill(MTrig);
	int ichange=-1;
	if(ChMax[1]==10)
		ichange=23;
	if(ChMax[1]==23)
		ichange=10;
	if(ichange!=-1)
		{
			ChMax[1]=ichange;
			if(ichange==23)Linha[1]=2;
			if(ichange==10)Linha[1]=3;
		}
	if(tDisplay>250)
	{
		if(isave==1)cout<<"stop saving images at "<<imu<<" of "<<T->GetEntries()<<endl;
		isave=0;
		//break;//stupid protection
	}
      /*    
	    if(ChMax[3]<999 && 
	    sqrt(pow(Linha[2]/10.-(Linha[1]+Linha[3])/20.,2.)
	    +pow(Coluna[2]/10.-(Coluna[1]+Coluna[3])/20.,2.))<1.){//muon
	    if(ChMax[1]<999 && ChMax[3]<999 && ChMax[2]<999 && ChMax[0]<999){//4 fold
      */
      if(ChMax[1]<999 || ChMax[3]<999){{//almost all 
	//if(ChMax[2]<999 || ChMax[0]<999){//at least 3 fold
	  
	  bool quad1=(ChMax[1]<999 && abs(Coluna[1])<4 && abs(Linha[1])<4);
	  bool quad2=(ChMax[2]<999 && abs(Coluna[2])<4 && abs(Linha[2])<4);
	  bool quad3=(ChMax[3]<999 && abs(Coluna[3])<4 && abs(Linha[3])<4);
	  //---testeJunho if(!quad1 || !quad2 || !quad3)continue;
	  float dirLL=0., dirCC=0.; 
	  float gL0=0,gC0=0,gL3=0,gC3=0;
	  if((quad1 && quad2) || (quad2 && quad3) || (quad1 && quad3)){//direction from central pads
	    if(quad1 && quad3){
	      dirLL=(Linha[1]-Linha[3])/2.; dirCC=(Coluna[1]-Coluna[3])/2.;
	      float gL2=Linha[1]-dirLL;
	      float gC2=Coluna[1]-dirCC;
	      gC0=Coluna[1]+dirCC;
	      gL0=Linha[1]+dirLL;
	      gC3=Coluna[3];
	      gL3=Linha[3];
	      int gib0=P0Empty->FindBin(gC0*quadr,gL0*quadr);
	      int gib2=strangeP->FindBin(gC2*quadr,gL2*quadr);
	      /*
	      if(abs(Coluna[0]+1-gib0)>1.)cout<<"gib0="<<gib0<<" C0+1="<<Coluna[0]+1<<endl;
	      if(ChMax[2]!=gib2-1)cout<<"gib2="<<gib2<<" Ch2+1="<<ChMax[2]+1<<endl;
	      */
	      //--if(abs(Coluna[0]+1-gib0)>1.)continue;
	      //if(ChMax[2]!=gib2-1)continue;
	      if(abs(gL2-Linha[2])>1 || abs(gC2-Coluna[2])>1)continue; //dia 28...
	    }//q1q3
	    else if(quad1 && quad2){
	      //if(ChMax[3]==60 || ChMax[3]==6)continue;
	      dirLL=(Linha[1]-Linha[2]); dirCC=(Coluna[1]-Coluna[2]);
	      gL3=Linha[1]-2.*dirLL;
	      gC3=Coluna[1]-2.*dirCC;
	      gC0=Coluna[1]+dirCC;
	      gL0=Linha[1]+dirLL;
	      //test if in correct channel!!!
	      int gib0=P0Empty->FindBin(gC0*quadr,gL0*quadr);
	      int gib3=strangeP->FindBin(gC3*quadr,gL3*quadr);
	      //--if(abs(Coluna[0]+1-gib0)>1.)continue;//cout<<"gib0="<<gib0<<" C0+1="<<Coluna[0]+1<<endl;
		//continue;
	      if(ChMax[3]!=999 && ChMax[3]!=gib3-1)continue;//cout<<"gib3="<<gib3<<" Ch3+1="<<ChMax[3]+1<<endl;
	    }//q1q2
	    else if(quad3 && quad2){
	      //if(ChMax[1]==60 || ChMax[1]==6)continue;
	      dirLL=(Linha[2]-Linha[3]); dirCC=(Coluna[2]-Coluna[3]);
	      float gL1=Linha[3]+2.*dirLL;
	      float gC1=Coluna[3]+2.*dirCC;
	      gC0=Coluna[3]+3.*dirCC;
	      gL0=Linha[3]+3.*dirLL;
	      gC3=Coluna[3];
	      gL3=Linha[3];
	      //float dist0=abs((Coluna[0]-32)/64.-gC0*quadr);
	      //test if in correct channel!!!
	      int gib1=strangeP->FindBin(gC1*quadr,gL1*quadr);
	      int gib0=P0Empty->FindBin(gC0*quadr,gL0*quadr);
	      //if(abs(Coluna[0]+1-gib0)>1.)continue;
	      //---if(abs(Coluna[0]+1-gib0)>1.)continue;//cout<<"gib0="<<gib0<<" C0+1="<<Coluna[0]+1<<endl;
	      if(ChMax[1]!=999 && ChMax[1]!=gib1-1)continue;//cout<<"gib1="<<gib1<<" Ch1+1="<<ChMax[1]+1<<endl;
	      //if(ChMax[1]!=gib1-1 || dist0>1.)continue;
	    }//q2q3
	    else
	      continue;
	     
	    ch0->Fill(Coluna[0]);
	    ch1->Fill(ChMax[1]);
	    ch2->Fill(ChMax[2]);
	    ch3->Fill(ChMax[3]);
	     
	   //acrescentar eficiencias
	   //nos verticais acrescentar factor geometrico
	  //nas direçoes e tido em conta no fluxo
		float xe1=xE1->GetBinContent(Coluna[1]+4,Linha[1]+4);
		float xe2=xE2->GetBinContent(Coluna[2]+4,Linha[2]+4);
		float xe3=xE3->GetBinContent(Coluna[3]+4,Linha[3]+4);
		float xgc=xGC->GetBinContent(Coluna[1]+4,Linha[1]+4);
		float xeee=xe1*xe2*xe3;

		//if(tDisplay<10)cout<<"xeee="<<xeee<<" gc="<<xgc<<endl;
	    if(quad1 && quad2 && quad3)
		{
		Direction->Fill(Coluna[3]-Coluna[1],Linha[3]-Linha[1]);
		EDirection->Fill(Coluna[3]-Coluna[1],Linha[3]-Linha[1],1./xeee);
		}
	    if(ChMax[3]==ChMax[2] && ChMax[2]==ChMax[1])
		{
		Vertical->Fill(Coluna[1],Linha[1]);
		EVertical->Fill(Coluna[1],Linha[1],xgc/xeee);
		}
	    if(isave==0 || tDisplay>500)continue;	
	    DIRMAP->Fill(dirCC*2.,dirLL*2.);
	    DIRMAP0->Reset(); DIRMAP0->Fill(dirCC*2.,dirLL*2.);
	     
	    double x=0,y=0;
	    for(int il=0;il<64;il++)
	      for(int ic=0;ic<64;ic++)
		{
		  y=-50.+binwidth/2.+il*binwidth;
		  x=-50.+binwidth/2.+ic*binwidth;
		  if(P0Empty->FindBin(x,y)==Coluna[0]+1) Disp3D->Fill(x,y,0.);
		  if(strangeP->FindBin(x,y)==ChMax[1]+1) Disp3D->Fill(x,y,-10.);
		  if(strangeP->FindBin(x,y)==ChMax[2]+1) Disp3D->Fill(x,y,-20.);
		  if(strangeP->FindBin(x,y)==ChMax[3]+1) Disp3D->Fill(x,y,-30.);
		}
	     
	    ch0T->Reset();
	    ch1T->Reset();
	    ch2T->Reset();
	    ch3T->Reset();
	    
/*
	    for(int ik=0;ik<64;ik++)
	      {
		if(Carga[0][ik]>0. && Carga[0][ik]<999.)
		  ch0T->Fill(CColuna[0][ik]);
		if(Carga[1][ik]>0. && Carga[1][ik]<999.)
		  ch1T->Fill(ik);
		if(Carga[2][ik]>0. && Carga[2][ik]<999.)
		  ch2T->Fill(ik);
		if(Carga[3][ik]>0. && Carga[3][ik]<999.)
		  ch3T->Fill(ik);
	      }
*/
	    gStyle->SetOptLogz(0);
	    Disp->Divide(4,1);
	    Disp->cd(1)->SetPad(0.0,0.00,0.2,1.00); Disp->cd(1)->Divide(1,4); 
	    
	    Plano0->ClearBinContents(); 
	    Plano1->ClearBinContents();
	    Plano2->ClearBinContents();
	    Plano3->ClearBinContents();
	    
	    Plano0->SetBinContent(Coluna[0]+1,1); 
	    Plano1->SetBinContent(ChMax[1]+1,1);
	    Plano2->SetBinContent(ChMax[2]+1,1);
	    Plano3->SetBinContent(ChMax[3]+1,1);
	    
	    Plano0->SetMarkerColor(2);	
	    Plano1->SetMarkerColor(2);	
	    Plano2->SetMarkerColor(2);	
	    Plano3->SetMarkerColor(2);	
	    
	    double vmax=0, vmin=1e9; 
	    DisPl0->ClearBinContents();
	    DisPl1->ClearBinContents();
	    DisPl2->ClearBinContents();
	    DisPl3->ClearBinContents();
	    
	    for(int ip=0;ip<64;ip++)
	      {
		DisPl0->SetBinContent(ip+1,ch0T->GetBinContent(ip+1)*10.);
		DisPl1->SetBinContent(ip+1,ch1T->GetBinContent(ip+1)*10.);
		DisPl2->SetBinContent(ip+1,ch2T->GetBinContent(ip+1)*10.);
		DisPl3->SetBinContent(ip+1,ch3T->GetBinContent(ip+1)*10.);
	      }
	    
	    
	    Disp->cd(1)->cd(1)->SetPad(0.00,0.75,0.75,1.00); DisPl0->DrawClone("text,col");
	    Plano0->SetMinimum(0.1); Plano0->SetMaximum(2); Plano0->DrawClone("text,col,same"); 
	    Disp->cd(1)->cd(2)->SetPad(0.00,0.50,0.75,0.75); DisPl1->DrawClone("text,col");
	    Plano1->SetMinimum(0.1); Plano1->SetMaximum(2); Plano1->DrawClone("text,col,same"); 
	    Disp->cd(1)->cd(3)->SetPad(0.00,0.25,0.75,0.50);  DisPl2->DrawClone("text,col");
	    Plano2->SetMinimum(0.1); Plano2->SetMaximum(2); Plano2->DrawClone("text,col,same"); 
	    Disp->cd(1)->cd(4)->SetPad(0.00,0.00,0.75,0.25); DisPl3->DrawClone("text,col");
	    Plano3->SetMinimum(0.1); Plano3->SetMaximum(2); Plano3->DrawClone("text,col,same"); 
	    
	    Disp->cd(3)->SetPad(0.41,0.00,0.82,1.00); 
	    Disp->cd(3)->SetPhi(45.); 
	    Disp3D->DrawClone("box1");
	    int ano,mes,dia;
	    ano=(int)(Data/10000); mes=(int)(Data-ano*10000)/100; dia=(int)(Data-ano*10000-mes*100);
	    int h,m, s; h=(int)(Hora); m=(int)(60.*(Hora-h)); s=(int)(60.*(60.*(Hora-h)-m));
	    /*
	    Hora=Hor+(Min+Seg/60)/60;
	    (Hora-Hor)*60=Min+Seg/60;
	    ((Hora-Hor)*60-Min)*60=Seg;
	    */
	    stringstream das; das<<"Data: "<<dia<<"/"<<mes<<"/20"<<ano;
	    stringstream hos; 
	    if(m>9)
		{
		 if(s>9)
			hos<<" Hora: "<<h<<":"<<m<<":"<<s;
		 else
			hos<<" Hora: "<<h<<":"<<m<<":0"<<s;
		}
	    else
		{
		 if(s>9)
			hos<<" Hora: "<<h<<":0"<<m<<":"<<s;
		 else
			hos<<" Hora: "<<h<<":0"<<m<<":0"<<s;
		}

	    if(first_Data==0)first_Data=Data;
	    if(first_Hora==0)first_Hora=Hora;
	    TText *t1=new TText(0.,0.95,das.str().c_str()); 
	    t1->SetNDC();
	    t1->DrawClone();
	    TText *t2=new TText(0.,0.90,hos.str().c_str()); 
	    t2->SetNDC();
	    t2->DrawClone();





	    
	    Disp->cd(2)->SetPad(0.20,0.22,0.40,0.78); Disp->cd(2)->SetGridy();
	    TH2F *D_zx=(TH2F*)Disp3D->Project3D("zx"); D_zx->SetTitle("");
									 
	    D_zx->GetXaxis()->SetLabelSize(0.06); D_zx->GetXaxis()->SetLabelOffset(-0.02);
	    D_zx->GetXaxis()->SetTitleSize(0.07); D_zx->GetXaxis()->SetTitleOffset(0.65);

	    D_zx->GetYaxis()->SetLabelSize(0.06); D_zx->GetYaxis()->SetLabelOffset(0.01);
	    D_zx->GetYaxis()->SetTitleSize(0.0); D_zx->GetYaxis()->SetTitleOffset(0.70);


									 
	    D_zx->SetMaximum(1); D_zx->DrawClone("box"); 
	    TLine L2(gC0*quadr,0.,gC3*quadr,-30.); L2.DrawClone();
	    //
	    
	    Disp->cd(4)->SetPad(0.82,0.22,1.00,0.78); Disp->cd(4)->SetGridy();
	    TH2F *D_zy=(TH2F*)Disp3D->Project3D("zy"); D_zy->SetTitle("");
	    //	    D_zy->GetXaxis()->SetTitle("Coluna (Sul - Norte)"); 

	    D_zy->GetXaxis()->SetLabelSize(0.06); D_zy->GetXaxis()->SetLabelOffset(-0.02);
	    D_zy->GetXaxis()->SetTitleSize(0.07); D_zy->GetXaxis()->SetTitleOffset(0.65);

	    D_zy->GetYaxis()->SetLabelSize(0.06); D_zy->GetYaxis()->SetLabelOffset(0.01);
	    D_zy->GetYaxis()->SetTitleSize(0.0); D_zy->GetYaxis()->SetTitleOffset(0.70);
	    
	    D_zy->SetMaximum(1); D_zy->DrawClone("box"); 
	    TLine L1(gL0*quadr,0.,gL3*quadr,-30.); L1.DrawClone();


	    //Axis:
	    //SetNdivisions(9)
	    //SetLableSize(0.06); SetLableOffset(0.0) in y; in x  -0.02
	    //SetTitleOffset(0.5); SetTitleSize(0.075)
	    
	    //Title: SetSize(45);
	    //Plano: SetNdivisions(0)


	    //	    cout<<"Time since last="<<(Hora-last_Hora)*3600.<<" s"<<endl;
	    Disp->Update(); 


	    if(doAll!=1)
	      {
		isave=0;
		cout<<"Save (1)? Quit(9)? Or (0) continue..."<<endl;
		cin>>isave;
	      }
		if(isave==1)
		  {
		    stringstream name; name<<"MuTom_Ev"<<tDisplay<<".png";
		    Disp->SaveAs(name.str().c_str());
		  }
		else if(isave==9)
		  return;
	      
	    Disp->Clear();

	     tDisplay++; //cout<<"Event "<<tDisplay<<" ("<<imu<<")"<<endl;
	     //cout<<"Change last_Hora="<<last_Hora<<" to "<<Hora<<endl;
	     num_min++;
	     if(m!=last_min)
	       {
		 last_min=m;
		 cout<<"Rate="<<num_min<<" per minute"<<endl;
		 num_min=0;
	       }
	  }}}//}}
    }
	

	cout<<"Reached end of loop"<<endl;
        TCanvas *Maps=new TCanvas("Maps","Maps",600.,1000.); Maps->Divide(2,2);

/*
	    int ano,mes,dia;
	    ano=(int)(Data/10000); mes=(int)(Data-ano*10000)/100; dia=(int)(Data-ano*10000-mes*100);
	    int h,m, s; h=(int)(Hora); m=(int)(60.*(Hora-h)); s=(int)(60.*(60.*(Hora-h)-m));
	    
	    int fano,fmes,fdia;
	    fano=(int)(first_Data/10000); fmes=(int)(first_Data-fano*10000)/100; fdia=(int)(first_Data-fano*10000-fmes*100);
	    int fh,fm; fh=(int)(first_Hora); fm=(int)(60.*(first_Hora-fh)); 
	    
	stringstream fromto;
	fromto<<"de "<<fdia<<"."<<fmes<<" ; "<<fh<<":"<<fm<<" a "<<dia<<"."<<mes<<" ; "<<hora<<":"<<min;  

	Direction->SetTitle(fromto.str().c_str());
	Vertical->SetTitle(fromto.str().c_str());
*/

	Maps->cd(1); Direction->Draw("colz");
	Maps->cd(2); Vertical->Draw("colz");
	Maps->cd(3); EDirection->Draw("colz");
	Maps->cd(4); EVertical->Draw("colz");



	cout<<"EDO: Write MapsNow from "<<first_Data<<" to "<<last_Data<<endl;
	//cout<<"do we have the maps?"<<endl; f->ls();
	TFile *MapsNow=new TFile("MapsNow.root","recreate");
        TTree *RunNumbers=new TTree("RunNumbers","RunNumbers");
        RunNumbers->Branch("first_Data",&first_Data,"first_Data/I");
        RunNumbers->Branch("last_Data",&last_Data,"last_Data/I");
        RunNumbers->Branch("first_Hora",&first_Hora,"first_Hora/F");
        RunNumbers->Branch("last_Hora",&last_Hora,"last_Hora/F");
	RunNumbers->Fill();
	RunNumbers->Write();
	Direction->Write();
	Vertical->Write();
	EDirection->Write();
	EVertical->Write();
	hMtrig->Write();

       TProfile	*MaxCharge0=(TProfile*)f->Get("MaxCharge0"); MaxCharge0->Write();
       TProfile	*MaxCharge1=(TProfile*)f->Get("MaxCharge1"); MaxCharge1->Write();
       TProfile	*MaxCharge2=(TProfile*)f->Get("MaxCharge2"); MaxCharge2->Write();
       TProfile	*MaxCharge3=(TProfile*)f->Get("MaxCharge3"); MaxCharge3->Write();
//
       TProfile	*OthCharge0=(TProfile*)f->Get("OthCharge0"); OthCharge0->Write();
       TProfile	*OthCharge1=(TProfile*)f->Get("OthCharge1"); OthCharge1->Write();
       TProfile	*OthCharge2=(TProfile*)f->Get("OthCharge2"); OthCharge2->Write();
       TProfile	*OthCharge3=(TProfile*)f->Get("OthCharge3"); OthCharge3->Write();
//
       TProfile	*PedCharge0=(TProfile*)f->Get("PedCharge0"); PedCharge0->Write();
       TProfile	*PedCharge1=(TProfile*)f->Get("PedCharge1"); PedCharge1->Write();
       TProfile	*PedCharge2=(TProfile*)f->Get("PedCharge2"); PedCharge2->Write();
       TProfile	*PedCharge3=(TProfile*)f->Get("PedCharge3"); PedCharge3->Write();

/*
1	ccpp
  KEY: TProfile	OthCharge0;1	ccpp
  KEY: TProfile	PedCharge0;1	ccpp
  KEY: TProfile	MaxCharge1;1	ccpp
  KEY: TProfile	OthCharge1;1	ccpp
  KEY: TProfile	PedCharge1;1	ccpp
  KEY: TProfile	MaxCharge2;1	ccpp
  KEY: TProfile	OthCharge2;1	ccpp
  KEY: TProfile	PedCharge2;1	ccpp
  KEY: TProfile	MaxCharge3;1	ccpp
  KEY: TProfile	OthCharge3;1	ccpp
  KEY: TProfile	PedCharge3;1	ccpp
r
*/
/*
	TH2Poly *MeanPed0=(TH2Poly*)f->Get("MeanPed0"); MeanPed0->Write();
	TH2Poly *MeanPed1=(TH2Poly*)f->Get("MeanPed1"); MeanPed1->Write();
	TH2Poly *MeanPed2=(TH2Poly*)f->Get("MeanPed2"); MeanPed2->Write();
	TH2Poly *MeanPed3=(TH2Poly*)f->Get("MeanPed3"); MeanPed3->Write();
	TH2Poly *SigPed0=(TH2Poly*)f->Get("SigPed0"); SigPed0->Write();
	TH2Poly *SigPed1=(TH2Poly*)f->Get("SigPed1"); SigPed1->Write();
	TH2Poly *SigPed2=(TH2Poly*)f->Get("SigPed2"); SigPed2->Write();
	TH2Poly *SigPed3=(TH2Poly*)f->Get("SigPed3"); SigPed3->Write();
	TH2Poly *ChargeMap0=(TH2Poly*)f->Get("ChargeMap0"); ChargeMap0->Write();
	TH2Poly *ChargeMap1=(TH2Poly*)f->Get("ChargeMap1"); ChargeMap1->Write();
	TH2Poly *ChargeMap2=(TH2Poly*)f->Get("ChargeMap2"); ChargeMap2->Write();
	TH2Poly *ChargeMap3=(TH2Poly*)f->Get("ChargeMap3"); ChargeMap3->Write();
*/
	//MapsNow->ls();
	MapsNow->Close();

  return;
}