#include "BTJuly2006.h"

using namespace std;

ClassImp(BTJuly2006)

//============================================================================//

BTJuly2006::BTJuly2006()
{
  NCH=4;
  NF=0;
  datfile=NULL;
  datobj=NULL;
}

//============================================================================//

BTJuly2006::BTJuly2006(string config,string rootfileRW)
{
  NCH=4;
  CONFIG=config;
  ROOTFILE=rootfileRW;
  //sprintf(path,"/media/usbdisk/data/bcm/CERNbt_July2006/");
  sprintf(path,"/afs/f9.ijs.si/home/irena/Delo/phd_analiza/Beam-Tests/CERNbt_july2006/2deff/");
  Int_t read_tracks=1;  
  Char_t pedestal[300]="calib_250806_1420.txt";

  Int_t Polarity[4];
  Float_t low[4],up[4],width_thr[4];
  Float_t noise_time[4],sig_dt[4];



  /************** JULY 2006 ***************************************/
  if(CONFIG=="config1-T11"){
    RUN.push_back(43);
    NINOthr.push_back(75);
    RUN.push_back(44);
    NINOthr.push_back(101.8);
    RUN.push_back(45);
    NINOthr.push_back(64.9);
    RUN.push_back(46);
    NINOthr.push_back(54.6);
    RUN.push_back(47);
    NINOthr.push_back(85.3);
    RUN.push_back(48);
    NINOthr.push_back(85.3);
    RUN.push_back(83);
    NINOthr.push_back(50.3);
    RUN.push_back(84);
    NINOthr.push_back(39.9);
    RUN.push_back(85);
    NINOthr.push_back(29.9);
    RUN.push_back(88);
    NINOthr.push_back(29.9);
    ADC_ChInput[0]="NINO chA, F402 at -1000V";
    ADC_ChInput[1]="NINO chB, F403 at -1000V";
    ADC_ChInput[2]="F404 at +1000V";
    ADC_ChInput[3]="F405 at +1000V";
    for(Int_t ich=0;ich<NCH;ich++){
      Polarity[ich]=1;
      low[ich]=500;
      noise_time[ich]=550;}
    DistZ[0]=226.0, DistZ[1]=187.0;
    DistZ[2]=148.0, DistZ[3]=109.0;
    for(Int_t ii=0;ii<(RUN.size());ii++){
      Pos_x1.push_back(std::vector<Float_t>());
      Pos_y1.push_back(std::vector<Float_t>());
      if(RUN[ii]<83){
	Pos_x1.at(ii).push_back(0), Pos_y1.at(ii).push_back(0);
	Pos_x1.at(ii).push_back(0), Pos_y1.at(ii).push_back(0);
	Pos_x1.at(ii).push_back(7.2), Pos_y1.at(ii).push_back(-20.8);
	Pos_x1.at(ii).push_back(8), Pos_y1.at(ii).push_back(-21.2);}
      else{
	Pos_x1.at(ii).push_back(10), Pos_y1.at(ii).push_back(-22.1);
	Pos_x1.at(ii).push_back(8.1), Pos_y1.at(ii).push_back(-22.4);
	Pos_x1.at(ii).push_back(7.2), Pos_y1.at(ii).push_back(-20.8);
	Pos_x1.at(ii).push_back(8), Pos_y1.at(ii).push_back(-21.2);}
    }

  }
  else if(CONFIG=="config2-T11"){
    RUN.push_back(49);
    NINOthr.push_back(59.9);
    RUN.push_back(50);
    NINOthr.push_back(59.9);
    RUN.push_back(51);
    NINOthr.push_back(100.2);
    RUN.push_back(52);
    NINOthr.push_back(100.2);
    RUN.push_back(54);
    NINOthr.push_back(100.2);
    RUN.push_back(55);
    NINOthr.push_back(100.2);
    RUN.push_back(57);
    NINOthr.push_back(100.2);
    RUN.push_back(59);
    NINOthr.push_back(150.3);
    RUN.push_back(61);
    NINOthr.push_back(128.5);
    RUN.push_back(62);
    NINOthr.push_back(128.5);
    RUN.push_back(68);
    NINOthr.push_back(50.6);
    RUN.push_back(76);
    NINOthr.push_back(63.6);
    RUN.push_back(77);
    NINOthr.push_back(63.6);
    RUN.push_back(78);
    NINOthr.push_back(63.6);
    ADC_ChInput[0]="F402 at -1000V";
    ADC_ChInput[1]="F403 at -1000V";
    ADC_ChInput[2]="NINO chA, F404 at +1000V";
    ADC_ChInput[3]="NINO chB, F405 at +1000V";
    for(Int_t ich=0;ich<NCH;ich++){
      Polarity[ich]=1;
      low[ich]=500;
      noise_time[ich]=550;}
    DistZ[0]=226.0, DistZ[1]=187.0;
    DistZ[2]=148.0, DistZ[3]=109.0;
    for(Int_t ii=0;ii<RUN.size();ii++){
      Pos_x1.push_back(std::vector<Float_t>());
      Pos_y1.push_back(std::vector<Float_t>());
      Pos_x1.at(ii).push_back(9.9), Pos_y1.at(ii).push_back(-22.2);
      Pos_x1.at(ii).push_back(8.2), Pos_y1.at(ii).push_back(-22.3);
      Pos_x1.at(ii).push_back(7.25), Pos_y1.at(ii).push_back(-21.);
      Pos_x1.at(ii).push_back(8.25), Pos_y1.at(ii).push_back(-21.);
    }
    for(Int_t ii=0;ii<RUN.size();ii++){
      for(Int_t ich=0;ich<NCH;ich++){
	width_thr[ich]=675;
	if((RUN[ii]<55)&&(ich==2)){
	  width_thr[ich]=475;
	}}
    }
  }
  else if(CONFIG=="config1-T9"){ //nino ch1, fot ch3 ok; NINo ch0 zlo sumi...
    RUN.push_back(105);   //samo analog, ch3 ok
    NINOthr.push_back(55);
    RUN.push_back(117);   
    NINOthr.push_back(60.8);
    //     RUN.push_back(118);  manjka track
    //     NINOthr.push_back(60.8); 
    //     RUN.push_back(138);  ni cel run syncd s trigg.
    //     NINOthr.push_back(49.7);
    RUN.push_back(145);
    NINOthr.push_back(49.7);
    RUN.push_back(154);
    NINOthr.push_back(20);
    RUN.push_back(155);
    NINOthr.push_back(50);
    //     RUN.push_back(156); manjka
    //     NINOthr.push_back(30);
    //     RUN.push_back(162);  bad run
    //     NINOthr.push_back(30);
    //     RUN.push_back(177);  //bad run
    //     NINOthr.push_back(30);
    RUN.push_back(178);
    NINOthr.push_back(70.2);
    RUN.push_back(180);
    NINOthr.push_back(70.2);
    RUN.push_back(181);
    NINOthr.push_back(9.7);
    ADC_ChInput[0]="NINO chA, F402 at -1000V";
    ADC_ChInput[1]="NINO chB, F403 at -1000V";
    ADC_ChInput[2]="F404 at +1000V";
    ADC_ChInput[3]="F405 at +1000V";
    for(Int_t ich=0;ich<NCH;ich++){
      Polarity[ich]=1;
      low[ich]=500;
      noise_time[ich]=550;}
    DistZ[0]=238.3, DistZ[1]=199.3;
    DistZ[2]=160.3, DistZ[3]=121.3;
    for(Int_t ii=0;ii<RUN.size();ii++){
      Pos_x1.push_back(std::vector<Float_t>());
      Pos_y1.push_back(std::vector<Float_t>());
      if(RUN[ii]==105){
	Pos_x1.at(ii).push_back(0), Pos_y1.at(ii).push_back(0);
	Pos_x1.at(ii).push_back(0), Pos_y1.at(ii).push_back(0);
	Pos_x1.at(ii).push_back(0), Pos_y1.at(ii).push_back(0);
	Pos_x1.at(ii).push_back(8.25), Pos_y1.at(ii).push_back(-21.);}
      else{
	Pos_x1.at(ii).push_back(20.), Pos_y1.at(ii).push_back(-19.5);
	Pos_x1.at(ii).push_back(18.25), Pos_y1.at(ii).push_back(-20.25);
	Pos_x1.at(ii).push_back(0), Pos_y1.at(ii).push_back(0);
	Pos_x1.at(ii).push_back(18), Pos_y1.at(ii).push_back(-19.1);}
    }
  }
  else if(CONFIG=="config2-T9"){
    RUN.push_back(184);
    NINOthr.push_back(30);
//     RUN.push_back(187);
//     NINOthr.push_back(20.0);
//     RUN.push_back(188);
//     NINOthr.push_back(20.0); //manjka
    RUN.push_back(189);
    NINOthr.push_back(10.2);
    RUN.push_back(190);
    NINOthr.push_back(40);
    RUN.push_back(191);
    NINOthr.push_back(50);
    RUN.push_back(192);
    NINOthr.push_back(60);
    ADC_ChInput[0]="F402 at -1000V";
    ADC_ChInput[1]="F403 at -1000V";
    ADC_ChInput[2]="NINO chA, F404 at +1000V";
    ADC_ChInput[3]="NINO chB, F405 at +1000V";
    for(Int_t ich=0;ich<NCH;ich++){
      Polarity[ich]=1;
      low[ich]=500;
      noise_time[ich]=550;}
    DistZ[0]=238.3, DistZ[1]=199.3;
    DistZ[2]=160.3, DistZ[3]=121.3;
    for(Int_t ii=0;ii<RUN.size();ii++){
      Pos_x1.push_back(std::vector<Float_t>());
      Pos_y1.push_back(std::vector<Float_t>());
      Pos_x1.at(ii).push_back(20.), Pos_y1.at(ii).push_back(-19.5);
      Pos_x1.at(ii).push_back(18.1), Pos_y1.at(ii).push_back(-20.1);
      Pos_x1.at(ii).push_back(17.2), Pos_y1.at(ii).push_back(-18.5);
      Pos_x1.at(ii).push_back(18), Pos_y1.at(ii).push_back(-18.9);
    }
  }
  else if(CONFIG=="config3-T9"){
//     RUN.push_back(194); manjka
//     NINOthr.push_back(40);
    RUN.push_back(195);
    NINOthr.push_back(59.8);
    RUN.push_back(196);
    NINOthr.push_back(59.8);
    RUN.push_back(197);
    NINOthr.push_back(59.8);
    RUN.push_back(198);
    NINOthr.push_back(50);
    RUN.push_back(199);
    NINOthr.push_back(50);
    RUN.push_back(200);
    NINOthr.push_back(50);
//     RUN.push_back(201); manjka
//     NINOthr.push_back(60);
//     RUN.push_back(202);
//     NINOthr.push_back(60);
    RUN.push_back(204);
    NINOthr.push_back(70);  //na html pise 60???????????
    RUN.push_back(205);
    NINOthr.push_back(29.9);
    RUN.push_back(206);
    NINOthr.push_back(19.9);
//     RUN.push_back(210); manjka
//     NINOthr.push_back(19.9);
    ADC_ChInput[0]="F406 at +1000V";
    ADC_ChInput[1]="F407 at -1000V";
    ADC_ChInput[2]="NINO chA, F409 at -1000V";
    ADC_ChInput[3]="NINO chB, F413 at +1000V";
    for(Int_t ich=0;ich<NCH;ich++){
      Polarity[ich]=1;
      low[ich]=800;
      noise_time[ich]=850;}
    DistZ[0]=247.3, DistZ[1]=208.3;
    DistZ[2]=169.3, DistZ[3]=130.3;
    for(Int_t ii=0;ii<RUN.size();ii++){
      Pos_x1.push_back(std::vector<Float_t>());
      Pos_y1.push_back(std::vector<Float_t>());
      if(RUN[ii]==195){
	Pos_x1.at(ii).push_back(18.65), Pos_y1.at(ii).push_back(-18.3);
	Pos_x1.at(ii).push_back(18.25), Pos_y1.at(ii).push_back(-18.3);
	Pos_x1.at(ii).push_back(18.9), Pos_y1.at(ii).push_back(-18.5);
	Pos_x1.at(ii).push_back(19.2), Pos_y1.at(ii).push_back(-18.);}
      else if((RUN[ii]==196) || (RUN[ii]==197)|| (RUN[ii]==198)){
	Pos_x1.at(ii).push_back(18.25), Pos_y1.at(ii).push_back(-22); //-21.5);
	Pos_x1.at(ii).push_back(18.25), Pos_y1.at(ii).push_back(-22);//-21.5);
	Pos_x1.at(ii).push_back(18.75), Pos_y1.at(ii).push_back(-22); //Pos_x1.at(ii).push_back(19), Pos_y1.at(ii).push_back(-21.5);
	Pos_x1.at(ii).push_back(18.75), Pos_y1.at(ii).push_back(-21.6);} //Pos_x1.at(ii).push_back(19), Pos_y1.at(ii).push_back(-21.5);}
      else{
	Pos_x1.at(ii).push_back(18.25), Pos_y1.at(ii).push_back(-16.5);
	Pos_x1.at(ii).push_back(18.25), Pos_y1.at(ii).push_back(-16.5);
	Pos_x1.at(ii).push_back(19), Pos_y1.at(ii).push_back(-16.7);
	Pos_x1.at(ii).push_back(19), Pos_y1.at(ii).push_back(-16.5);}
    }
  }
  else if(CONFIG=="config4-T9"){
    RUN.push_back(213);
    NINOthr.push_back(19.9);
    RUN.push_back(214);
    NINOthr.push_back(29.9);
    RUN.push_back(215);
    NINOthr.push_back(40);
    RUN.push_back(216);
    NINOthr.push_back(50);
    RUN.push_back(219);
    NINOthr.push_back(60.2);
    RUN.push_back(220);
    NINOthr.push_back(80);
    RUN.push_back(221);
    NINOthr.push_back(80);
    RUN.push_back(222);
    NINOthr.push_back(109.8);
    RUN.push_back(223);
    NINOthr.push_back(150);
    RUN.push_back(224);
    NINOthr.push_back(200);
    ADC_ChInput[0]="NINO chA, F406 at +1000V";
    ADC_ChInput[1]="NINO chB, F407 at -1000V";
    ADC_ChInput[2]="F409 at -1000V";
    ADC_ChInput[3]="F413 at +1000V";
    for(Int_t ich=0;ich<NCH;ich++){
      Polarity[ich]=1;
      low[ich]=500;
      noise_time[ich]=550;}
    DistZ[0]=247.3, DistZ[1]=208.3;
    DistZ[2]=169.3, DistZ[3]=130.3;
    for(Int_t ii=0;ii<RUN.size();ii++){
      Pos_x1.push_back(std::vector<Float_t>());
      Pos_y1.push_back(std::vector<Float_t>());
      if(RUN[ii]==213){
	Pos_x1.at(ii).push_back(18.4), Pos_y1.at(ii).push_back(-16.7);
	Pos_x1.at(ii).push_back(18.4), Pos_y1.at(ii).push_back(-16.7);
	Pos_x1.at(ii).push_back(19), Pos_y1.at(ii).push_back(-16.75);
	Pos_x1.at(ii).push_back(0), Pos_y1.at(ii).push_back(0);}
      else{
	Pos_x1.at(ii).push_back(18.4), Pos_y1.at(ii).push_back(-16.6);
	Pos_x1.at(ii).push_back(18.4), Pos_y1.at(ii).push_back(-16.6);
	Pos_x1.at(ii).push_back(19), Pos_y1.at(ii).push_back(-16.75);
	Pos_x1.at(ii).push_back(19), Pos_y1.at(ii).push_back(-16.3);}
    }
  }
  else if(CONFIG=="config5-T9"){
//     RUN.push_back(225); prazen..
//     NINOthr.push_back(20.0); 
    RUN.push_back(233);
    NINOthr.push_back(20.9);
    RUN.push_back(234);
    NINOthr.push_back(30.4);
    RUN.push_back(235);
    NINOthr.push_back(60);
    ADC_ChInput[0]="NINO chA, F408 at +1000V";
    ADC_ChInput[1]="NINO chB, F403 at -1000V";
    ADC_ChInput[2]="F404 at -1000V";
    ADC_ChInput[3]="F405 at +1000V";
    for(Int_t ich=0;ich<NCH;ich++){
      Polarity[ich]=1;
      low[ich]=800;
      noise_time[ich]=850;}
    DistZ[0]=283.3, DistZ[1]=199.3;
    DistZ[2]=160.3, DistZ[3]=121.3;
    for(Int_t ii=0;ii<RUN.size();ii++){
      Pos_x1.push_back(std::vector<Float_t>());
      Pos_y1.push_back(std::vector<Float_t>());
      Pos_x1.at(ii).push_back(16.25), Pos_y1.at(ii).push_back(-17.3);
      Pos_x1.at(ii).push_back(17), Pos_y1.at(ii).push_back(-17.5);
      Pos_x1.at(ii).push_back(17.3), Pos_y1.at(ii).push_back(-16.5);
      Pos_x1.at(ii).push_back(18.7), Pos_y1.at(ii).push_back(-16.6);
    }
  }
  else if(CONFIG=="config6-T9"){
    RUN.push_back(236);
    NINOthr.push_back(80.0);
    ADC_ChInput[0]="F408 at +1000V";
    ADC_ChInput[1]="F403 at -1000V";
    ADC_ChInput[2]="NINO chA, F404 at -1000V";
    ADC_ChInput[3]="NINO chB, F405 at +1000V";
    for(Int_t ich=0;ich<NCH;ich++){
      Polarity[ich]=1;
      low[ich]=500;
      noise_time[ich]=550;}
    DistZ[0]=283.3, DistZ[1]=199.3;
    DistZ[2]=160.3, DistZ[3]=121.3;
    for(Int_t ii=0;ii<RUN.size();ii++){
      Pos_x1.push_back(std::vector<Float_t>());
      Pos_y1.push_back(std::vector<Float_t>());
      Pos_x1.at(ii).push_back(16.2), Pos_y1.at(ii).push_back(-17.1);
      Pos_x1.at(ii).push_back(17), Pos_y1.at(ii).push_back(-17.6);
      Pos_x1.at(ii).push_back(17.2), Pos_y1.at(ii).push_back(-16.5);
      Pos_x1.at(ii).push_back(18.8), Pos_y1.at(ii).push_back(-16.5);
    }
  }

  NF=RUN.size();
  Pos_dx=2.25;
  Pos_dy=5.25;


  for(Int_t ii=0;ii<NF;ii++)
    StartEvent.push_back(-1);

  Float_t dBL=50;
  for(Int_t ich=0;ich<NCH;ich++){
    up[ich]=low[ich]+dBL;
    Int_t pos=ADC_ChInput[ich].find("NINO");
    if(pos<ADC_ChInput[ich].size()){
      NINOch.push_back(ich);
      cout<<"NINO ch="<<ich<<endl;
      sig_dt[ich]=20.0;
      if(CONFIG!="config2-T11")
	width_thr[ich]=675;
    }
    else{
      FOTECch.push_back(ich);
      sig_dt[ich]=2.0;
      width_thr[ich]=200;
      cout<<"Fotec ch="<<ich<<endl;
    }
  }

  datobj=new CERNbt*[NF];
  datfile=new TFile*[NF];
  Char_t *ftmp;
  Char_t *rtxt;
  if(ROOTFILE=="WRITE"){
    for(Int_t ifl=0;ifl<NF;ifl++){
      rtxt=new Char_t[1000];
      sprintf(rtxt,"%sroot-files/run%d.root",path,RUN[ifl]);
      datfile[ifl]=new TFile(rtxt,"recreate");
      delete []rtxt;
      rtxt=new Char_t[200];   //obj's key
      sprintf(rtxt,"CERNjuly2006_Obj_RUN%d",RUN[ifl]);
      datobj[ifl]=new CERNbt(path,pedestal,RUN[ifl],ADC_ChInput,read_tracks,2,ftmp);
      for(Int_t kk=0;kk<NINOch.size();kk++){
	Int_t ich=NINOch[kk];
	datobj[ifl]->BLcorr[ich]=1;
      }
      for(Int_t kk=0;kk<FOTECch.size();kk++){
	Int_t ich=FOTECch[kk];
	datobj[ifl]->BLcorr[ich]=0;
      }	
      for(Int_t ich=0;ich<NCH;ich++){
	datobj[ifl]->pol[ich]=Polarity[ich];
	datobj[ifl]->BLup[ich]=up[ich];
	datobj[ifl]->BLdown[ich]=low[ich];
	cout<<"ich="<<ich<<"  blcorr="<<datobj[ifl]->BLcorr[ich]<<"  thr_witdh="<<width_thr[ich]<<" sig_dt="<<sig_dt[ich]<<endl;}
      datobj[ifl]->Fill_Tree(sig_dt,noise_time,width_thr);
      datobj[ifl]->Write(rtxt);
      delete []rtxt;
    }}
  else if(ROOTFILE=="READ"){
    for(Int_t ifl=0;ifl<NF;ifl++){
      rtxt=new Char_t[1000];
      sprintf(rtxt,"%sroot-files/run%d.root",path,RUN[ifl]);
      datfile[ifl]=new TFile(rtxt,"read");
      delete []rtxt;
      rtxt=new Char_t[200];   //obj's key
      sprintf(rtxt,"CERNjuly2006_Obj_RUN%d",RUN[ifl]);
      datobj[ifl]=(CERNbt*)datfile[ifl]->Get(rtxt);
      delete []rtxt;}
  }
  else{
    cout<<"kaj z obj???"<<endl;
    exit(0);}
  for(Int_t ifl=0;ifl<NF;ifl++){
    cout<<"RUN="<<RUN[ifl]<<":"<<endl;
    cout<<"\tfotec CHs:";
    for(Int_t kk=0;kk<FOTECch.size();kk++)
      cout<<FOTECch[kk]<<",";
    cout<<"\tNINO CHs:";
    for(Int_t kk=0;kk<NINOch.size();kk++)
      cout<<NINOch[kk]<<",";
    cout<<endl;
    for(Int_t ich=0;ich<NCH;ich++){
      cout<<"\tich="<<ich<<"  blcorr="<<datobj[0]->BLcorr[ich]<<"  thr_witdh="<<width_thr[ich];
      cout<<" sig_dt="<<sig_dt[ich]<<" noise_t="<<noise_time[ich]<<" BL=["<<datobj[ifl]->BLdown[ich]<<",";
      cout<<datobj[ifl]->BLup[ich]<<"]"<<endl;
    }
  }

}

//============================================================================//

BTJuly2006::~BTJuly2006()
{
//   if(datobj!=NULL){
//     for(Int_t ii=0;ii<NF;ii++){
//       //cout<<"ii="<<ii<<endl;
//       delete datobj[ii];}
//     delete []datobj;
//   }
  if(datfile!=NULL){
    for(Int_t ii=0;ii<NF;ii++){
      //cout<<"ii="<<ii<<endl;
      delete datfile[ii];}
    delete []datfile;   
  }
}

//============================================================================//

TH1F **BTJuly2006::Get_FotecNoiseSpectra(string cut_on_tracks)
//cut_on_tracks= "YES" or "NO"
{
  Float_t xmax=200.0,xmin=-200.0;
  Int_t nbins=101;
  Float_t dx=(xmax-xmin)/(nbins-1.0);
  TH1F *htmp=new TH1F("htmp","htmp",nbins,xmin,xmax);

  TH1F **hst=new TH1F*[FOTECch.size()];
  Char_t *proj,*CUT;
  for(Int_t ich=0;ich<FOTECch.size();ich++){
    Int_t jj=FOTECch[ich];
    proj=new Char_t[200];
    sprintf(proj,"ch%d, %s",jj,ADC_ChInput[jj].c_str());
    hst[ich]=new TH1F(proj,proj,nbins,xmin,xmax);
    delete []proj;
    for(Int_t ii=0;ii<NF;ii++){
      proj=new Char_t[50];
      sprintf(proj,"Noise%d",jj+1);
      //cut
      if(cut_on_tracks=="YES"){
	CUT=new Char_t[5000];
	sprintf(CUT,"((XY2+XYslope2*%f)>=%f)",DistZ[ich],Pos_y1[ii][ich]);            //y1
	sprintf(CUT,"%s && ((XY2+XYslope2*%f)<=%f)",CUT,DistZ[ich],Pos_y1[ii][ich]+Pos_dy);   //y2
	sprintf(CUT,"%s && ((XY1+XYslope1*%f)>=%f)",CUT,DistZ[ich],Pos_x1[ii][ich]);          //x1
	sprintf(CUT,"%s && ((XY1+XYslope1*%f)<=%f)",CUT,DistZ[ich],Pos_x1[ii][ich]+Pos_dx);  //x2
	//sprintf(CUT,"%s && (Event>%d)",CUT,StartEvent[ii]); 
	datobj[ii]->WFdat->Project(htmp->GetName(),proj,CUT);
	cout<<"run="<<RUN[ii]<<"  ch="<<ich<<": "<<ADC_ChInput[ich].c_str()<<endl;
	cout<<"\tcut = "<<CUT<<endl;
	delete []CUT;
      }
      else if(cut_on_tracks=="NO"){
	CUT=new Char_t[500];
	sprintf(CUT,"(Event>%d)",StartEvent[ii]); 
	datobj[ii]->WFdat->Project(htmp->GetName(),proj,CUT);
	delete []CUT;
      }
      else{
	cout<<"do cut="<<cut_on_tracks<<"??? not definded"<<endl;
	exit(0);}
      hst[ich]->Add(htmp);
      htmp->Reset();
    }
  }

  delete htmp;
  return hst;
}

//============================================================================//

TH1F **BTJuly2006::Get_FotecSignalSpectra(string cut_on_tracks)
//cut_on_tracks= "YES" or "NO"
{
  Float_t xmax=1900.0,xmin=-200.0;
  Int_t nbins=211;//96; //191;
  Float_t dx=(xmax-xmin)/(nbins-1.0);
  TH1F *htmp=new TH1F("htmp","htmp",nbins,xmin,xmax);

  TH1F **hst=new TH1F*[FOTECch.size()];
  Char_t *proj,*CUT;
  for(Int_t ich=0;ich<FOTECch.size();ich++){
    Int_t jj=FOTECch[ich];
    proj=new Char_t[200];
    sprintf(proj,"ch%d, %s",ich,ADC_ChInput[jj].c_str());
    hst[ich]=new TH1F(proj,proj,nbins,xmin,xmax);
    delete []proj;
    for(Int_t ii=0;ii<NF;ii++){
      proj=new Char_t[50];
      //projection
      if(datobj[ii]->pol[jj]>0)	
	sprintf(proj,"Ampl%d",jj+1);
      else if(datobj[ii]->pol[jj]<0) 
	sprintf(proj,"-Ampl%d",jj+1);
      else{
	cout<<"polarity="<<datobj[ii]->pol[jj]<<"  ch="<<jj<<endl;
	exit(0);}
      //cut
      if(cut_on_tracks=="YES"){
	CUT=new Char_t[5000];
	sprintf(CUT,"((XY2+XYslope2*%f)>=%f)",DistZ[ich],Pos_y1[ii][jj]);            //y1
	sprintf(CUT,"%s && ((XY2+XYslope2*%f)<=%f)",CUT,DistZ[ich],Pos_y1[ii][jj]+Pos_dy);   //y2
	sprintf(CUT,"%s && ((XY1+XYslope1*%f)>=%f)",CUT,DistZ[jj],Pos_x1[ii][jj]);          //x1
	sprintf(CUT,"%s && ((XY1+XYslope1*%f)<=%f)",CUT,DistZ[jj],Pos_x1[ii][jj]+Pos_dx);  //x2
	//sprintf(CUT,"%s && (Event>%d)",CUT,StartEvent[ii]); 
	datobj[ii]->WFdat->Project(htmp->GetName(),proj,CUT);
	cout<<"run="<<RUN[ii]<<"  ch="<<jj<<": "<<ADC_ChInput[jj].c_str()<<endl;
	cout<<"\tcut = "<<CUT<<endl;
	delete []CUT;
      }
      else if(cut_on_tracks=="NO"){
	CUT=new Char_t[500];
	sprintf(CUT,"(Event>%d)",StartEvent[ii]); 
	datobj[ii]->WFdat->Project(htmp->GetName(),proj,CUT);
	delete []CUT;
      }
      else{
	cout<<"do cut="<<cut_on_tracks<<"??? not definded"<<endl;
	exit(0);}
      hst[ich]->Add(htmp);
      htmp->Reset();
    }
  }

  delete htmp;
  return hst;
}

//============================================================================//

TH1F **BTJuly2006::Get_NINOSignalSpectra(string cut_on_tracks)
//cut_on_tracks= "YES" or "NO"
{
  Float_t xmax=1700.0,xmin=-200.0;
  Int_t nbins=110;//95;
  Float_t dx=(xmax-xmin)/(nbins-1.0);
  TH1F *htmp=new TH1F("htmp","htmp",nbins,xmin,xmax);

  TH1F **hst=new TH1F*[NINOch.size()];
  Char_t *proj,*CUT;
  for(Int_t ich=0;ich<NINOch.size();ich++){
    Int_t jj=NINOch[ich];
    proj=new Char_t[200];
    sprintf(proj,"ch%d, %s",jj,ADC_ChInput[jj].c_str());
    hst[ich]=new TH1F(proj,proj,nbins,xmin,xmax);
    delete []proj;
    for(Int_t ii=0;ii<NF;ii++){
      proj=new Char_t[50];
      //projection
      if(datobj[ii]->pol[jj]>0)	
	sprintf(proj,"Ampl%d",jj+1);
      else if(datobj[ii]->pol[jj]<0) 
	sprintf(proj,"-Ampl%d",jj+1);
      else{
	cout<<"polarity="<<datobj[ii]->pol[jj]<<"  ch="<<jj<<endl;
	exit(0);}
      //cut
      if(cut_on_tracks=="YES"){
	CUT=new Char_t[5000];
	sprintf(CUT,"((XY2+XYslope2*%f)>=%f)",DistZ[jj],Pos_y1[ii][jj]);            //y1
	sprintf(CUT,"%s && ((XY2+XYslope2*%f)<=%f)",CUT,DistZ[jj],Pos_y1[ii][jj]+Pos_dy);   //y2
	sprintf(CUT,"%s && ((XY1+XYslope1*%f)>=%f)",CUT,DistZ[jj],Pos_x1[ii][jj]);          //x1
	sprintf(CUT,"%s && ((XY1+XYslope1*%f)<=%f)",CUT,DistZ[jj],Pos_x1[ii][jj]+Pos_dx);  //x2
	sprintf(CUT,"%s && (Event>%d)",CUT,StartEvent[ii]); 
	datobj[ii]->WFdat->Project(htmp->GetName(),proj,CUT);
	cout<<"run="<<RUN[ii]<<"  ch="<<jj<<": "<<ADC_ChInput[jj].c_str()<<endl;
	cout<<"\tcut = "<<CUT<<endl;
	delete []CUT;
      }
      else if(cut_on_tracks=="NO"){
	CUT=new Char_t[500];
	sprintf(CUT,"(Event>%d)",StartEvent[ii]); 
	datobj[ii]->WFdat->Project(htmp->GetName(),proj,CUT);
	delete []CUT;
      }
      else{
	cout<<"do cut="<<cut_on_tracks<<"??? not definded"<<endl;
	exit(0);}
      hst[ich]->Add(htmp);
      htmp->Reset();
    }
  }

  delete htmp;
  return hst;
}

//============================================================================//

TGraphErrors **BTJuly2006::Get_NINOefficiency(string cut_on_tracks,Int_t drw)
//cut_on_tracks= "YES" or "NO"
{
  TGraphErrors **gtmp=new TGraphErrors*[NINOch.size()];
  for(Int_t ii=0;ii<NINOch.size();ii++){
    gtmp[ii]=new TGraphErrors();
    gtmp[ii]->GetHistogram()->SetTitle(ADC_ChInput[NINOch[ii]].c_str());
    gtmp[ii]->GetHistogram()->SetXTitle("NINO threshold [mV]");
    gtmp[ii]->GetHistogram()->SetYTitle("Efficiency");
    gtmp[ii]->GetHistogram()->SetBins(100,0,200);
  }

  Float_t NINOcut1=900,NINOcut2=90;
  Float_t AMPLcut[4]={0,0,0,0};
  for(Int_t ich=0;ich<NINOch.size();ich++){
    AMPLcut[ich]=NINOcut1;
    Int_t jj=NINOch[ich];
    Int_t pos1=ADC_ChInput[jj].find("opto");
    if(pos1<ADC_ChInput[jj].size())
      AMPLcut[ich]=NINOcut2; 
  }

  TH1F *hst=new TH1F("hstmp","hstmp",1801,-200.5,1600.5);
  Char_t *proj,*CUT;
  for(Int_t ii=0;ii<NF;ii++){
    for(Int_t ich=0;ich<NINOch.size();ich++){
      proj=new Char_t[10];
      Int_t jj=NINOch[ich];
      //projection
      if(datobj[ii]->pol[jj]>0)	
	sprintf(proj,"Ampl%d",jj+1);
      else if(datobj[ii]->pol[jj]<0) 
	sprintf(proj,"-Ampl%d",jj+1);
      else{
	cout<<"polarity="<<datobj[ii]->pol[jj]<<"  ch="<<jj<<endl;
	exit(0);}
      //cut
      if(cut_on_tracks=="YES"){
	CUT=new Char_t[1000];
	sprintf(CUT,"((XY2+XYslope2*%f)>=%f)",DistZ[jj],Pos_y1[ii][jj]);            //y1
	sprintf(CUT,"%s && ((XY2+XYslope2*%f)<=%f)",CUT,DistZ[jj],Pos_y1[ii][jj]+Pos_dy);   //y2
	sprintf(CUT,"%s && ((XY1+XYslope1*%f)>=%f)",CUT,DistZ[jj],Pos_x1[ii][jj]);          //x1
	sprintf(CUT,"%s && ((XY1+XYslope1*%f)<=%f)",CUT,DistZ[jj],Pos_x1[ii][jj]+Pos_dx);  //x2
	sprintf(CUT,"%s && (Event>%d) && (BLrms%d<10)",CUT,StartEvent[ii],jj+1); 
	datobj[ii]->WFdat->Project(hst->GetName(),proj,CUT);
	if(drw==1){
	  cout<<"run="<<RUN[ii]<<"  ch="<<jj<<": "<<ADC_ChInput[jj].c_str()<<endl;
	  cout<<"\tcut = "<<CUT<<endl;}
	delete []CUT;
      }
      else if(cut_on_tracks=="NO"){
	CUT=new Char_t[500];
	sprintf(CUT,"(Event>%d)",StartEvent[ii]); 
	datobj[ii]->WFdat->Project(hst->GetName(),proj,CUT);
	delete []CUT;
      }
      else{
	cout<<"do cut="<<cut_on_tracks<<"??? not definded"<<endl;
	exit(0);}
      Double_t all=hst->GetEntries();
      hst->SetAxisRange(AMPLcut[ich],2000);
      Double_t yes=hst->Integral();
      if(drw) cout<<"\tyes="<<yes<<" all="<<all<<" cut="<<AMPLcut[ich]<<endl;
      Int_t ipt=gtmp[ich]->GetN();
      if((all!=0) && (yes!=0)){
	if((CONFIG=="config3-T9") &&(RUN[ii]==195)){
	  //ne risi
	}
	else{
	  gtmp[ich]->SetPoint(ipt,NINOthr[ii],(yes/all));
	  Float_t err=TMath::Sqrt(yes*(1.0-yes/all))/all;
	  gtmp[ich]->SetPointError(ipt,0,err);
	  if(drw) 
	    cout<<" eff="<<yes/all<<" +/- "<<err<<" THR="<<NINOthr[ii]<<endl;
	}
      }
      delete []proj;
      if(drw==1){
	hst->SetAxisRange(-200,5000);
	hst->Draw();
	gPad->Update();
	cin.get();}
      hst->Reset();
    }}
  delete hst;

  return gtmp;
}

//============================================================================//

TH2F ***BTJuly2006::Get_HitMap()
{
  Float_t xx1=-5,xx2=35;
  Float_t yy1=-35,yy2=5;
  Int_t nbx=100,nby=100;

  Float_t NINOcut1=900,FOTECcut=100,NINOcut2=90;

  Char_t *htit,*proj;
  TH2F ***hmap=new TH2F**[NF+1];
  for(Int_t ii=0;ii<NF;ii++){
    hmap[ii]=new TH2F*[NCH];
    for(Int_t ich=0;ich<NCH;ich++){
      Float_t amplcut=NINOcut1;
      Int_t pos1=ADC_ChInput[ich].find("opto");
      if(pos1<ADC_ChInput[ich].size())
	amplcut=NINOcut2; 
      else{
	for(Int_t kch=0;kch<FOTECch.size();kch++){
	  if(ich==FOTECch[kch]){
	    amplcut=FOTECcut;
	    break;
	  }}
      }
      htit=new Char_t[100];
      sprintf(htit,"Hit Map, CH%d, %s, run%d, signal cut %.1f",ich,ADC_ChInput[ich].c_str(),RUN[ii],amplcut);
      hmap[ii][ich]=new TH2F(htit,htit,nbx,xx1,xx2,nby,yy1,yy2);
      hmap[ii][ich]->SetXTitle("X [mm]");
      hmap[ii][ich]->SetYTitle("Y [mm]");
      delete []htit;
      htit=new Char_t[500];
      if(datobj[ii]->pol[ich]>0)
	sprintf(htit,"(Ampl%d>%.1f)",ich+1,amplcut);
      else if(datobj[ii]->pol[ich]<0)
	sprintf(htit,"((-Ampl%d)>%.1f)",ich+1,amplcut);
      proj=new Char_t[1000];
      sprintf(proj,"(XY2+XYslope2*%f):(XY1+XYslope1*%f)",DistZ[ich],DistZ[ich]);
      cout<<"proj="<<proj<<" cut="<<htit<<endl;
      datobj[ii]->WFdat->Project(hmap[ii][ich]->GetName(),proj,htit);
      delete []htit;
      delete []proj;
    }
    htit=new Char_t[100];
    sprintf(htit,"Hit Map All, run%d",RUN[ii]);
    hmap[ii][4]=new TH2F(htit,htit,nbx,xx1,xx2,nby,yy1,yy2);
    delete []htit;
    proj=new Char_t[1000];
    sprintf(proj,"(XY2+XYslope2*%f):(XY1+XYslope1*%f)",0,0);
    datobj[ii]->WFdat->Project(hmap[ii][4]->GetName(),proj);
    delete []proj;
  }
  
  return hmap;
}
//============================================================================//

TH2F **BTJuly2006::Get_FotecHitMap(Float_t *amplcut)
{
  Float_t xx1=-5,xx2=35;
  Float_t yy1=-35,yy2=5;
  Int_t nbx=100,nby=100;
  TH2F *htm=new TH2F("htit","htit",nbx,xx1,xx2,nby,yy1,yy2);

  Char_t *htit,*proj;
  TH2F **hmap=new TH2F*[FOTECch.size()];
  for(Int_t jj=0;jj<FOTECch.size();jj++){
    Int_t ich=FOTECch[jj];
    htit=new Char_t[500];
    sprintf(htit,"Hit Map, CH%d, %s, signal cut %.1f",ich,ADC_ChInput[ich].c_str(),amplcut[jj]);
    hmap[jj]=new TH2F(htit,htit,nbx,xx1,xx2,nby,yy1,yy2);
    hmap[jj]->SetXTitle("X [mm]");
    hmap[jj]->SetYTitle("Y [mm]");
    delete []htit;
    for(Int_t ii=0;ii<NF;ii++){
      htit=new Char_t[500];
      if(datobj[ii]->pol[ich]>0)
	sprintf(htit,"(Ampl%d>%.1f)",ich+1,amplcut[jj]);
      else if(datobj[ii]->pol[ich]<0)
	sprintf(htit,"((-Ampl%d)>%.1f)",ich+1,amplcut[jj]);
      proj=new Char_t[1000];
      sprintf(proj,"(XY2+XYslope2*%f):(XY1+XYslope1*%f)",DistZ[ich],DistZ[ich]);
      datobj[ii]->WFdat->Project(htm->GetName(),proj,htit);
      delete []htit;
      delete []proj;
      hmap[jj]->Add(htm);
      htm->Reset();
    }
  }
  delete htm;  
  return hmap;
}

//============================================================================//

TH2F **BTJuly2006::Get_NINOHitMap(Float_t *amplcut)
{
  Float_t xx1=-5,xx2=35;
  Float_t yy1=-35,yy2=5;
  Int_t nbx=100,nby=100;

  TH2F *htm=new TH2F("htit","htit",nbx,xx1,xx2,nby,yy1,yy2);

  Char_t *htit,*proj;
  TH2F **hmap=new TH2F*[NINOch.size()];
  for(Int_t jj=0;jj<NINOch.size();jj++){
    Int_t ich=NINOch[jj];
    htit=new Char_t[500];
    sprintf(htit,"Hit Map, CH%d, %s, signal cut %.1f",ich,ADC_ChInput[ich].c_str(),amplcut[jj]);
    hmap[jj]=new TH2F(htit,htit,nbx,xx1,xx2,nby,yy1,yy2);
    hmap[jj]->SetXTitle("X [mm]");
    hmap[jj]->SetYTitle("Y [mm]");
    delete []htit;
    for(Int_t ii=0;ii<NF;ii++){
      htit=new Char_t[500];
      if(datobj[ii]->pol[ich]>0)
	sprintf(htit,"(Ampl%d>%.1f) && (Event>%d)",ich+1,amplcut[jj],StartEvent[ii]);
      else if(datobj[ii]->pol[ich]<0)
	sprintf(htit,"((-Ampl%d)>%.1f) && (Event>%d)",ich+1,amplcut[jj],StartEvent[ii]);
      proj=new Char_t[1000];
      sprintf(proj,"(XY2+XYslope2*%f):(XY1+XYslope1*%f)",DistZ[ich],DistZ[ich]);
      datobj[ii]->WFdat->Project(htm->GetName(),proj,htit);
      delete []htit;
      delete []proj;
      hmap[jj]->Add(htm);
      htm->Reset();
    }
  }
  delete htm;  
  return hmap;
}

//============================================================================//

TH2F **BTJuly2006::ScanZ_HitMap()
{
  TCanvas *ca=new TCanvas("ca","ca",800,800);
  ca->Divide(2,2);
  Float_t MinRMSx[4],MinRMSy[4];
  Int_t MinRMSxINDX[4],MinRMSyINDX[4];
  for(Int_t ich=0;ich<NCH;ich++){
    MinRMSx[ich]=1000;
    MinRMSy[ich]=1000;
    MinRMSxINDX[ich]=-1;
    MinRMSyINDX[ich]=-1;}

  Float_t xx1=-5,xx2=35;
  Float_t yy1=-35,yy2=5;
  Int_t nbx=100,nby=100;

  Float_t NINOcut1=900,FOTECcut=100,NINOcut2=90;
  TH2F *htmp=new TH2F("htht","htht",nbx,xx1,xx2,nby,yy1,yy2);

  Char_t *htit,*proj;
  TH2F **hmap=new TH2F*[NCH];
  Int_t nz=50;
  Float_t dz=10;

  for(Int_t ich=0;ich<NCH;ich++)
    DistZ[ich]=ich*39.0-dz; //BATmiddle+ich*39.0-dz;  //zacnem na zadnjem BATu (najblizji bcm)

  for(Int_t iz=0;iz<nz;iz++){
    for(Int_t ich=0;ich<NCH;ich++){
      DistZ[ich]=DistZ[ich]+dz;
      Float_t amplcut=NINOcut1;
      Int_t pos1=ADC_ChInput[ich].find("opto");
      if(pos1<ADC_ChInput[ich].size())
	amplcut=NINOcut2; 
      else{
	for(Int_t kch=0;kch<FOTECch.size();kch++){
	  if(ich==FOTECch[kch]){
	    amplcut=FOTECcut;
	    break;
	  }}
      }
      if(iz==0){
	htit=new Char_t[100];
	sprintf(htit,"Hit Map, CH%d, %s, signal cut %.1f",ich,ADC_ChInput[ich].c_str(),amplcut);
	hmap[ich]=new TH2F(htit,htit,nbx,xx1,xx2,nby,yy1,yy2);
	hmap[ich]->SetXTitle("X [mm]");
	hmap[ich]->SetYTitle("Y [mm]");
	delete []htit;}
      hmap[ich]->Reset();
      for(Int_t ifl=0;ifl<NF;ifl++){
	htit=new Char_t[50];
	if(datobj[ifl]->pol[ich]>0)
	  sprintf(htit,"Ampl%d>%.1f",ich+1,amplcut);
	else if(datobj[ifl]->pol[ich]<0)
	  sprintf(htit,"(-Ampl%d)>%.1f",ich+1,amplcut);
	proj=new Char_t[1000];
	sprintf(proj,"(XY2+XYslope2*%f):(XY1+XYslope1*%f)",DistZ[ich],DistZ[ich]);
	datobj[ifl]->WFdat->Project(htmp->GetName(),proj,htit);
	hmap[ich]->Add(htmp);
	htmp->Reset();
	delete []htit;
	delete []proj;
      }
      cout<<"ch="<<ich<<" iz="<<iz<<" z="<<DistZ[ich]<<" Xrms="<<hmap[ich]->GetRMS()<<" Yrms="<<hmap[ich]->GetRMS(2)<<endl;
      if(hmap[ich]->GetRMS()<MinRMSx[ich]){
	MinRMSxINDX[ich]=iz;
	MinRMSx[ich]=hmap[ich]->GetRMS();}
      if(hmap[ich]->GetRMS(2)<MinRMSy[ich]){
	MinRMSyINDX[ich]=iz;
	MinRMSy[ich]=hmap[ich]->GetRMS(2);}
      ca->cd(ich+1);
      hmap[ich]->SetMaximum(200);
      hmap[ich]->Draw("colz");
    }
    gPad->Update();
    cout<<endl;
    //cin.get();
  }
  delete htmp;

  for(Int_t ich=0;ich<NCH;ich++){
    cout<<"ch="<<ich<<":"<<endl;
    cout<<"\tminX indx="<<MinRMSxINDX[ich]<<" rmsx="<<MinRMSx[ich]<<"  z-batmiddle="<<ich*39.0+MinRMSxINDX[ich]*dz<<endl;
    cout<<"\tminY indx="<<MinRMSyINDX[ich]<<" rmsY="<<MinRMSy[ich]<<"  z-batmiddle="<<ich*39.0+MinRMSyINDX[ich]*dz<<endl;}
  return hmap;
}
//============================================================================//

TH2F **BTJuly2006::Get_Fotec3Defficiency(Float_t *amplcut)
{
  Float_t xx1=-25,xx2=-14;
  Float_t yy1=8,yy2=19;
  Float_t dxx=0.1,dyy=0.1; 
//   if(CONFIG=="config4-T9"){
//     xx1=12,xx2=27;
//     yy1=-21,yy2=-6;
//   }
//   xx1=8,xx2=28;
//   yy1=-25,yy2=-5;
  xx1=0,xx2=28;
  yy1=-30,yy2=-2;

  if(CONFIG=="config2-T9"){
    yy1=-24,yy2=-11;
    xx1=12,xx2=25;
  }
  if(CONFIG=="config4-T9"){
    yy1=-21,yy2=-7;
    xx1=13,xx2=27;
  }


  Int_t nbx=TMath::Nint((xx2-xx1)/dxx)+1;
  Int_t nby=TMath::Nint((yy2-yy1)/dyy)+1;

  Char_t *txt;
  TH2F **hmap=new TH2F*[FOTECch.size()];
  for(Int_t jj=0;jj<FOTECch.size();jj++){
    Int_t ich=FOTECch[jj];
    txt=new Char_t[1000];
    sprintf(txt,"Efficiency, ch%d %s, amplitude cut at %.0f",ich,ADC_ChInput[ich].c_str(),amplcut[jj]);
    hmap[jj]=new TH2F(txt,txt,nbx,xx1-dxx/2.0,xx2+dxx/2.0,nby,yy1-dyy/2.0,yy2+dyy/2.0);
    delete []txt;
  }

  Pos_dx=1.0;
  Pos_dy=1.0;
  Float_t startX=xx1-Pos_dx/2.0;
  Float_t startY=yy1-Pos_dy/2.0;
  for(Int_t kk=0;kk<NF;kk++){
    for(Int_t ich=0;ich<NCH;ich++){
      Pos_x1[kk][ich]=startX;
      Pos_y1[kk][ich]=startY;
    }}

  for(Int_t iy=0;iy<nby;iy++){
    //ystep
    for(Int_t kk=0;kk<NF;kk++){
      for(Int_t ich=0;ich<NCH;ich++){
	Pos_y1[kk][ich]=startY+iy*dyy;
	  }}
    for(Int_t ix=0;ix<nbx;ix++){
      //xstep
      for(Int_t kk=0;kk<NF;kk++){
	for(Int_t ich=0;ich<NCH;ich++){
	  Pos_x1[kk][ich]=startX+ix*dxx;
	}}
      for(Int_t kc=0;kc<FOTECch.size();kc++){
	Int_t ich=FOTECch[kc];
	Double_t all=Project_Ampl(ich,-100000);
	//cin.get();
	Double_t yes=Project_Ampl(ich,amplcut[kc]);
	//cin.get();
	Double_t eff=0;
	if(all>0)
	  eff=yes/all;
 	Double_t xx=xx1+ix*dxx;
 	Double_t yy=yy1+iy*dyy;
	hmap[kc]->SetBinContent(ix+1,iy+1,eff); //(xx,yy,eff);
	cout<<"CUT="<<amplcut[kc]<<"  ix="<<ix<<"(nbx) iy="<<iy<<"(nby)  all="<<all<<" eff="<<eff<<" posx="<<Pos_x1[0][0]<<" posy="<<Pos_y1[0][0]<<" x="<<xx<<"  y="<<yy<<endl;
// 	if(yy>11.5){
// 	  hmap[kc]->Draw("colz");
// 	  gPad->Update();
// 	  cin.get();
// 	}
      }
    }
  }
  return hmap;
}


//============================================================================//

TH2F ***BTJuly2006::Get_Fotec3Defficiency_v2(Float_t *amplcut,Float_t BoxWidth, Int_t &ngraphs)
{
  //number of output histos
  Int_t NGR=1;
  Int_t *RUN_start;
  Int_t *RUN_end;

  //dimensions of the histograms
  Float_t dxx=0.1,dyy=0.1; 
  Float_t  xx1=0,xx2=28;
  Float_t  yy1=-30,yy2=-2;

  if(CONFIG=="config2-T11"){
    xx1=2.0,xx2=16.0;
    yy1=-27.5;13.5;}
  if(CONFIG=="config2-T9"){
    yy1=-25,yy2=-11;
    xx1=12,xx2=26;}
  if(CONFIG=="config4-T9"){
    yy1=-21,yy2=-7;
    xx1=13,xx2=27;}
  if(CONFIG=="config5-T9"){
    yy1=-21,yy2=-7;
    xx1=13,xx2=27;}
  if(CONFIG=="config6-T9"){
    yy1=-22,yy2=-8;
    xx1=12,xx2=26;}


  // which runs to collect  
  if(CONFIG=="config1-T11"){
    NGR=2;
    RUN_start=new Int_t[NGR];
    RUN_end=new Int_t[NGR];
    RUN_start[0]=RUN[0];
    RUN_end[0]=46;  //,47 mal premaknen...48??
    RUN_start[1]=48;
    RUN_end[1]=RUN[NF-1];
    xx1=2.0,xx2=16.0;
    yy1=-25.5;11.5;
  }
  else if(CONFIG=="config1-T9"){
    NGR=1;
    RUN_start=new Int_t[NGR];
    RUN_end=new Int_t[NGR];
    RUN_start[0]=117;
    RUN_end[0]=RUN[NF-1];
    xx1=12.0,xx2=26.0;
    yy1=-24.0,yy2=-10.0;
  }
  else if(CONFIG=="config3-T9"){
    NGR=2;
    RUN_start=new Int_t[NGR];
    RUN_end=new Int_t[NGR];
    RUN_start[0]=196;
    RUN_end[0]=198;
    RUN_start[1]=199;
    RUN_end[1]=206;
    xx1=8.0,xx2=27.0;
    yy1=-26.0,yy2=-7.0;
  }
  else{
    RUN_start=new Int_t[NGR];
    RUN_end=new Int_t[NGR];
    RUN_start[0]=RUN[0];
    RUN_end[0]=RUN[NF-1];}


  Char_t *txt;
  Int_t nbx=TMath::Nint((xx2-xx1)/dxx)+1;
  Int_t nby=TMath::Nint((yy2-yy1)/dyy)+1;
  TH2F ***hmap=new TH2F**[NGR];
  for(Int_t ii=0;ii<NGR;ii++){
    hmap[ii]=new TH2F*[FOTECch.size()];
    for(Int_t jj=0;jj<FOTECch.size();jj++){
      Int_t ich=FOTECch[jj];
      txt=new Char_t[1000];
      sprintf(txt,"Efficiency%d, ch%d %s, amplitude cut at %.0f",ii,ich,ADC_ChInput[ich].c_str(),amplcut[jj]);
      hmap[ii][jj]=new TH2F(txt,txt,nbx,xx1-dxx/2.0,xx2+dxx/2.0,nby,yy1-dyy/2.0,yy2+dyy/2.0);
      delete []txt;
    }}

  //box size
  Pos_dx=BoxWidth;
  Pos_dy=BoxWidth;
  //start box position
  Float_t startX=xx1-Pos_dx/2.0;
  Float_t startY=yy1-Pos_dy/2.0;
  for(Int_t kk=0;kk<NF;kk++){
    for(Int_t ich=0;ich<NCH;ich++){
      Pos_x1[kk][ich]=startX;
      Pos_y1[kk][ich]=startY;
    }}

  //loop over steps
  for(Int_t jgr=0;jgr<NGR;jgr++){
    for(Int_t iy=0;iy<nby;iy++){
      //ystep
      for(Int_t kk=0;kk<NF;kk++){
	for(Int_t ich=0;ich<NCH;ich++){
	  Pos_y1[kk][ich]=startY+iy*dyy;
	}}
      for(Int_t ix=0;ix<nbx;ix++){
	//xstep
	for(Int_t kk=0;kk<NF;kk++){
	  for(Int_t ich=0;ich<NCH;ich++){
	    Pos_x1[kk][ich]=startX+ix*dxx;
	  }}
	for(Int_t kc=0;kc<FOTECch.size();kc++){
	  Int_t ich=FOTECch[kc];
	  Double_t all=Project_Ampl_v2(ich,-100000,RUN_start[jgr],RUN_end[jgr]);
	  Double_t yes=Project_Ampl_v2(ich,amplcut[kc],RUN_start[jgr],RUN_end[jgr]);
	  Double_t eff=0;
	  if(all>0)
	    eff=yes/all;
	  Double_t xx=xx1+ix*dxx;
	  Double_t yy=yy1+iy*dyy;
	  hmap[jgr][kc]->SetBinContent(ix+1,iy+1,eff); 
	  cout<<"CUT="<<amplcut[kc]<<"  jgr="<<jgr<<"  ix="<<ix<<"(nbx) iy="<<iy<<"(nby)  all="<<all<<"  y="<<yes;
	  cout<<" eff="<<eff<<" posx="<<Pos_x1[0][0]<<" posy="<<Pos_y1[0][0]<<" x="<<xx<<"  y="<<yy<<endl;
	}
      }
    }
  }

  delete []RUN_end;
  delete []RUN_start;
  return hmap;
}


//============================================================================//

TTree *BTJuly2006::Get_Fotec3Defficiency_v2_Tree(Float_t *MIP,Float_t BoxWidth)
{
  //steps and lower upper bounds for x&y
//   Float_t xx1=-24,xx2=-14;
//   Float_t yy1=8,yy2=18;
//   Float_t dxx=0.1,dyy=0.1; //??????????
//   if((CONFIG=="config12-H8") || (CONFIG=="config11-H8")){
//     xx1=13,xx2=21;
//     yy1=11,yy2=19;
//   }
  Float_t dxx=0.1,dyy=0.1; 
  Float_t  xx1=0,xx2=28;
  Float_t  yy1=-30,yy2=-2;
  if(CONFIG=="config2-T11"){
    xx1=2.0,xx2=16.0;
    yy1=-27.5;13.5;}
  if(CONFIG=="config2-T9"){
    yy1=-25,yy2=-11;
    xx1=12,xx2=26;}
  if(CONFIG=="config4-T9"){
    yy1=-21,yy2=-7;
    xx1=13,xx2=27;}
  if(CONFIG=="config5-T9"){
    yy1=-21,yy2=-7;
    xx1=13,xx2=27;}
  if(CONFIG=="config6-T9"){
    yy1=-22,yy2=-8;
    xx1=12,xx2=26;}


  Int_t nbx=TMath::Nint((xx2-xx1)/dxx)+1;
  Int_t nby=TMath::Nint((yy2-yy1)/dyy)+1;

  //signal cut
  Int_t Ncut=7;


  EFFdat_t effdat;
  Char_t *txt=new Char_t[100];
  sprintf(txt,"FotecEff_BT_CERNjuly2006_%s",CONFIG.c_str());
  TTree *TreeEffDat=new TTree(txt,txt);
  delete []txt;
  txt=new Char_t[1000];
  sprintf(txt,"RunIndx/F");
  sprintf(txt,"%s:XX",txt);
  sprintf(txt,"%s:YY",txt);
  for(Int_t ich=0;ich<FOTECch.size();ich++)
    sprintf(txt,"%s:All_CH%d",txt,ich+1);
  for(Int_t ich=0;ich<FOTECch.size();ich++){
    for(Int_t icut=0;icut<Ncut;icut++){
      sprintf(txt,"%s:Yes%d_CH%d",txt,icut+1,ich+1);}}
  TreeEffDat->Branch("2deff",&effdat.RunIndx,txt);
  delete []txt;


  Pos_dx=BoxWidth;
  Pos_dy=BoxWidth;
  Float_t startX=xx1-Pos_dx/2.0;
  Float_t startY=yy1-Pos_dy/2.0;
  for(Int_t kk=0;kk<NF;kk++){
    for(Int_t ich=0;ich<NCH;ich++){
      Pos_x1[kk][ich]=startX;
      Pos_y1[kk][ich]=startY;
    }}

  for(Int_t iy=0;iy<nby;iy++){
    //ystep
    for(Int_t kk=0;kk<NF;kk++){
      for(Int_t ich=0;ich<NCH;ich++){
	Pos_y1[kk][ich]=startY+iy*dyy;
	  }}
    for(Int_t ix=0;ix<nbx;ix++){
      //xstep
      for(Int_t kk=0;kk<NF;kk++){
	for(Int_t ich=0;ich<NCH;ich++){
	  Pos_x1[kk][ich]=startX+ix*dxx;
	}}
      Double_t xx=xx1+ix*dxx;
      Double_t yy=yy1+iy*dyy;
      effdat.XX=xx;
      effdat.YY=yy;
      cout<<"x="<<xx<<"  y="<<yy<<endl;
      for(Int_t irun=0;irun<NF;irun++){
	effdat.RunIndx=irun;
	for(Int_t kc=0;kc<FOTECch.size();kc++){
	  Int_t ich=FOTECch[kc];
	  Float_t startcut=MIP[kc]*0.25;
	  Float_t stepcut=MIP[kc]*0.25;
	  Double_t all=Project_Ampl(ich,-100000,irun);
	  effdat.All[kc]=all;	
	  for(Int_t icut=0;icut<Ncut;icut++){
	    Float_t sig_cut=startcut+stepcut*icut;
	    Double_t yes=Project_Ampl(ich,sig_cut,irun);
	    effdat.Yes[kc][icut]=(Float_t)yes;	    
	    //cout<<"CUT="<<sig_cut<<" ich="<<kc<<" run="<<irun<<"  all="<<all<<"  yes="<<yes<<" x="<<xx<<"  y="<<yy<<endl;
	  }}
	TreeEffDat->Fill();
      }
    }
  }
  return TreeEffDat;
}

//============================================================================//

Double_t BTJuly2006::Project_Ampl(Int_t ich,Float_t amplcut,Int_t irun)
{
  Float_t xmax=2000.0,xmin=-200.0;
  Int_t nbins=221;
  Float_t dx=(xmax-xmin)/(nbins-1.0);
  TH1F *hamp=new TH1F("htmp1","htmp1",nbins,xmin-dx/2.0,xmax+dx/2.0);

  Int_t drw=0;

  Char_t *proj,*CUT;
  proj=new Char_t[10];
  //projection
  if(datobj[irun]->pol[ich]>0)	
    sprintf(proj,"Ampl%d",ich+1);
  else if(datobj[irun]->pol[ich]<0) 
    sprintf(proj,"-Ampl%d",ich+1);
  else{
    cout<<"polarity="<<datobj[irun]->pol[ich]<<"  ch="<<ich<<endl;
    exit(0);}
  //cut
  CUT=new Char_t[5000];
  sprintf(CUT,"((XY2+XYslope2*%f)>=%f)",DistZ[ich],Pos_y1[irun][ich]);            //y1
  sprintf(CUT,"%s && ((XY2+XYslope2*%f)<=%f)",CUT,DistZ[ich],Pos_y1[irun][ich]+Pos_dy);   //y2
  sprintf(CUT,"%s && ((XY1+XYslope1*%f)>=%f)",CUT,DistZ[ich],Pos_x1[irun][ich]);          //x1
  sprintf(CUT,"%s && ((XY1+XYslope1*%f)<=%f)",CUT,DistZ[ich],Pos_x1[irun][ich]+Pos_dx);  //x2
  sprintf(CUT,"%s && (Event>%d) && (%s)>%f",CUT,StartEvent[irun],proj,amplcut); 
  datobj[irun]->WFdat->Project(hamp->GetName(),proj,CUT);
  if(drw==1)
    cout<<CUT<<endl;
  delete []CUT;
  delete []proj;

  Double_t fe=hamp->GetEntries();
  if(drw==1){
    cout<<"ich="<<ich<<"  cut="<<amplcut<<endl;
    hamp->Draw();
    gPad->Update();
    cin.get();}

  delete hamp;
  return fe;
}

//============================================================================//

Double_t BTJuly2006::Project_Ampl(Int_t ich,Float_t amplcut)
{
  Float_t xmax=2000.0,xmin=-200.0;
  Int_t nbins=221;
  Float_t dx=(xmax-xmin)/(nbins-1.0);
  TH1F *hamp=new TH1F("htmp1","htmp1",nbins,xmin-dx/2.0,xmax+dx/2.0);
  TH1F *htmp=new TH1F("htmp","htmp",nbins,xmin-dx/2.0,xmax+dx/2.0);

  Char_t *proj,*CUT;
  for(Int_t ii=0;ii<NF;ii++){
    proj=new Char_t[10];
    //projection
    if(datobj[ii]->pol[ich]>0)	
      sprintf(proj,"Ampl%d",ich+1);
    else if(datobj[ii]->pol[ich]<0) 
      sprintf(proj,"-Ampl%d",ich+1);
    else{
      cout<<"polarity="<<datobj[ii]->pol[ich]<<"  ch="<<ich<<endl;
      exit(0);}
    //cut
    CUT=new Char_t[5000];
    sprintf(CUT,"((XY2+XYslope2*%f)>=%f)",DistZ[ich],Pos_y1[ii][ich]);            //y1
    sprintf(CUT,"%s && ((XY2+XYslope2*%f)<=%f)",CUT,DistZ[ich],Pos_y1[ii][ich]+Pos_dy);   //y2
    sprintf(CUT,"%s && ((XY1+XYslope1*%f)>=%f)",CUT,DistZ[ich],Pos_x1[ii][ich]);          //x1
    sprintf(CUT,"%s && ((XY1+XYslope1*%f)<=%f)",CUT,DistZ[ich],Pos_x1[ii][ich]+Pos_dx);  //x2
    sprintf(CUT,"%s && (Event>%d) && (%s)>%f",CUT,StartEvent[ii],proj,amplcut); 
    //if((RUN[ii]!=213) && (ich!=3))
    if((RUN[ii]==213) && (ich==3)){
      Int_t bla=0;
    }
    else
      datobj[ii]->WFdat->Project(htmp->GetName(),proj,CUT);
    delete []CUT;
    delete []proj;
    hamp->Add(htmp);
    htmp->Reset();
  }

  Double_t fe=hamp->GetEntries();
  delete htmp;
  delete hamp;
  return fe;
}
//============================================================================//

Double_t BTJuly2006::Project_Ampl_v2(Int_t ich,Float_t amplcut,Int_t start_run,Int_t end_run)
{
  Float_t xmax=2000.0,xmin=-200.0;
  Int_t nbins=221;
  Float_t dx=(xmax-xmin)/(nbins-1.0);
  TH1F *hamp=new TH1F("htmp1","htmp1",nbins,xmin-dx/2.0,xmax+dx/2.0);
  TH1F *htmp=new TH1F("htmp","htmp",nbins,xmin-dx/2.0,xmax+dx/2.0);

  Char_t *proj,*CUT;
  for(Int_t ii=0;ii<NF;ii++){
    if((RUN[ii]>=start_run) && (RUN[ii]<=end_run)){
      proj=new Char_t[10];
      //projection
      if(datobj[ii]->pol[ich]>0)	
	sprintf(proj,"Ampl%d",ich+1);
      else if(datobj[ii]->pol[ich]<0) 
	sprintf(proj,"-Ampl%d",ich+1);
      else{
	cout<<"polarity="<<datobj[ii]->pol[ich]<<"  ch="<<ich<<endl;
	exit(0);}
      //cut
      CUT=new Char_t[5000];
      sprintf(CUT,"((XY2+XYslope2*%f)>=%f)",DistZ[ich],Pos_y1[ii][ich]);            //y1
      sprintf(CUT,"%s && ((XY2+XYslope2*%f)<=%f)",CUT,DistZ[ich],Pos_y1[ii][ich]+Pos_dy);   //y2
      sprintf(CUT,"%s && ((XY1+XYslope1*%f)>=%f)",CUT,DistZ[ich],Pos_x1[ii][ich]);          //x1
      sprintf(CUT,"%s && ((XY1+XYslope1*%f)<=%f)",CUT,DistZ[ich],Pos_x1[ii][ich]+Pos_dx);  //x2
      sprintf(CUT,"%s && (Event>%d) && (%s)>%f",CUT,StartEvent[ii],proj,amplcut); 
      //if((RUN[ii]!=213) && (ich!=3))
      if((RUN[ii]==213) && (ich==3)){
	Int_t bla=0;
      }
      else
	datobj[ii]->WFdat->Project(htmp->GetName(),proj,CUT);
      delete []CUT;
      delete []proj;
      hamp->Add(htmp);
      htmp->Reset();
    }
  }

  Double_t fe=hamp->GetEntries();
  delete htmp;
  delete hamp;

  return fe;
}

//============================================================================//
TGraph **BTJuly2006::Draw_Box(Float_t x1, Float_t dx,Float_t y1,Float_t dy,
			      Int_t col,Int_t draw=1)
{
  Float_t x2=x1+dx;
  Float_t y2=y1+dy;

  TGraph **box=new TGraph*[4];
  Int_t nnpt=11;
  for(int ii=0;ii<4;ii++) 
    box[ii]=new TGraph(nnpt);
 
  Double_t xx,yy;
  for(int ii=0;ii<nnpt;ii++){
    xx=x1+ii*(x2-x1)/(nnpt-1.0);
    yy=y1;
    box[0]->SetPoint(ii,xx,yy);
    yy=y2;
    box[1]->SetPoint(ii,xx,yy);
    yy=y1+ii*(y2-y1)/(nnpt-1.0);
    xx=x1;
    box[2]->SetPoint(ii,xx,yy);
    xx=x2;
    box[3]->SetPoint(ii,xx,yy);      
  }  
  for(int ii=0;ii<4;ii++){
    box[ii]->SetLineColor(col);
    box[ii]->SetLineWidth(1);
    if(draw==1) box[ii]->Draw("l");
  }
  return box;
}



/*============================================================================*/
/*           PRIVATE FUNCTIONS                                                */
/*============================================================================*/

//======================= END of CLASS =======================================//