#include <stdio.h>
#include <UnpackBasic.h>


void UnpackBasic(){
  gSystem->Load("libMbs");
  gSystem->Load("libMbsAPI");
  gSystem->Load("libTCHIEvent");
  gSystem->Load("libTMBALLEvent");
  gSystem->Load("libTLANDEvent");
  gSystem->Load("libTKRATEvent");
  gSystem->Load("libTSTELEvent");
  gSystem->Load("libTRootEvent");
 
  string dir = getenv("VMCWORKDIR");
  string landcaldir = dir + "/CALFILES/LAND/";

  TRootDefine chitree;
  TRootCHIEvent *evt = 0;
  chitree.Open("asyeos.root");
  chitree.GetTree()->SetTitle("AsyEos beam");
  evt = chitree.GetCHIEvt();

  int MBSTSh,MBSTSm,MBSTSs;
  int CHITSh,CHITSm,CHITSs;
  int Trig_Patt;
  int bit_Trig_Patt[16];

  vector<int> pchim;
  Int_t time, de, fastpg, fastgg, slowpg, slowgg, timecsi;
  Float_t x,y;
  TH2F* h2=new  TH2F("fast-slow ntel=212","fs 212",100,0.,4200,100,0.,4200.);
  TH2F* h22=new TH2F("fast-time ntel=212","ft 212",100,0.,4200,100,0.,4200.);

 

  TCHIEvent* chimera_evt = 0;   
  chimera_evt = new TCHIEvent("gsieos.par",2);   
  chimera_evt->ProcessParFile();
  chimera_evt->Get_Pedestal_Data("asyeos.ped");


  TMBALLEvent* mball_evt = 0;
  mball_evt = new TMBALLEvent();   
  mball_evt->LoadMapping("microball_mappingwTDC.txt");
  mball_evt->LoadThresholds("adc_thresh.dat");

  TLANDEvent* land_evt = 0;
  land_evt = new TLANDEvent(landcaldir);   
//  land_evt->PassLANDData(landhistlist);  
  land_evt->DoLANDMapping();//apply mapping!   

  // Krakow array objects. Begining
  Int_t iKrakowEv = 0;
  TKRATRawEvent *krakow_evt = 0;
  TKRATMapping KrattaMapping;
  KrattaMapping.Init( "/home/sebastian/Projects/asyeosroot/asyeosroot/macros/KRATTA_mapping.conf"); 
  KrattaMapping.PrintMapping();


  TSTELEvent* stel_evt = 0;
  stel_evt = new TSTELEvent();
  Int_t word[33];
  int ii;

  int Nruns;
  TString InFile ;
  ifstream inputfile;
  //inputfile.open("runlist_local.dat");
  inputfile.open("runlist.dat");
  inputfile >> Nruns;

 

  MbsFile* myMbs=0;
  Int_t *status=new Int_t;
  Int_t* sebuflenght=new Int_t;
  Short_t* setype=new Short_t;
  Short_t* sesubtype=new Short_t;
  Char_t* sesubcrate=new Char_t;

  Int_t* SubEventDataPtr=new Int_t;

  TH1F * hKrakow_Mult = new TH1F  ("hKrakow_Mult", "Multiplicity", 105, 0, 105);

for(int irun=0;irun<Nruns;irun++){
	inputfile >>InFile;
	myMbs= new MbsFile(InFile,status);
	cout << "$$$$$$$$$$$$$$$$$$$$$$$$$$$$" << endl;
	cout <<irun<<" "<<InFile<<endl;
	cout << "Press any key to start " << InFile << " file processing." << endl;
	getchar();

	if(!(*status)){
		myMbs->PrintFileHeader(); 
		Int_t lastevent=0;
		Int_t i=0;
		Int_t chievt=0;
		Int_t mbsevt=0;
		//   for(i=0;i<100000;i++){
		
		while(!lastevent){
			//    cout << "############################# i="<< i << endl;
			//   getchar();
			if(i%5000==0)cout << "############################# i="<< i << endl;
			lastevent=myMbs->GetEvent(); 
			Int_t nrSubEvts=myMbs->GetNrSubEvent();
 
			for(Int_t j=1;j<=nrSubEvts;j++){
				//    cout << "inside 2 prima "<<(*SubEventDataPtr)<<endl;
				SubEventDataPtr=myMbs->GetSubEvent(j,sebuflenght,setype,sesubtype,sesubcrate); 
				//    cout << "inside: subevent pointing to"<<(SubEventDataPtr)<<" lenght: "<<*sebuflenght<< endl;
				//imp    cout << "inside: nr of se="<<j<<" type="<<*setype<<" subtype="<<*sesubtype<<" subcrate="<<*sesubcrate<<endl;
				//    if(nrSubEvts==1) {

				/// CHIMERA Subvent:
				if((*setype==10) && (*sesubtype==1) && (*SubEventDataPtr==256)){
					//imp     cout << "Found a valid CHIMERA subevent OKKKKKKKKKKKKKKKKKKKKKKKKK"<<endl;
					Int_t ftotfired=chimera_evt->ProcessCurrentEventMBS(SubEventDataPtr); 
					chimera_evt->GetTiMeStamp(&CHITSs,&CHITSm,&CHITSh);
					chievt++;
					//imp     cout <<"CHIMERA TIMESTAMP= "<<CHITSs<<" "<<CHITSm<<" "<<CHITSh<<endl;
					//imp//imp     cout <<"CHIMERA TIMESTAMP= "<<(*(SubEventDataPtr+1) & 0x0000ffff)<<" "<<(*(SubEventDataPtr+2)& 0x0000ffff)<<" "<<(*(SubEventDataPtr+3) & 0x0000ffff)<<endl;
					
					//    cout << "fired= "<<ftotfired<< endl;
					multTOF = chimera_evt->GetTofMultiplicity();
					chimera_evt->GetFired(pchim);
							
					int ntel = pchim.size();
					int num = 0;
					int none;
					for(int k=0; k<ntel; k++) {
						// raw parameters
						int numtel = pchim[k];
						chimera_evt->Get_Data_Telescope(numtel, TOFSi,SI,  LG,LG, &time,   &de);
						chimera_evt->Get_Data_Telescope(numtel, SLOW, FAST,LG,LG, &slowpg, &fastpg);
						chimera_evt->Get_Data_Telescope(numtel, SLOW, FAST,HG,HG, &slowgg, &fastgg);
						chimera_evt->Get_Data_Telescope(numtel, TOF30,NONE,LG,LG, &timecsi,&none);
						//     if(timecsi==0 && time==0)continue;
						if(numtel==212 && slowpg  && fastpg)h2->Fill(slowpg,fastpg);
						if(numtel==212 && timecsi && fastpg)h22->Fill(timecsi,fastpg);
						evt->numtel[num]  = numtel;
						evt->desilpg[num] = de;
						evt->fastpg[num]  = fastpg;
						evt->slowpg[num]  = slowpg;
						evt->fastgg[num]  = fastgg;
						evt->slowgg[num]  = slowgg;
						evt->timesi[num]  = time;
						evt->timecsi[num] = timecsi;		
						num++;
					}
					evt->multi = num;
					
					//    getchar();
					chitree.GetTree()->Fill();
					chimera_evt->Clear_fired_param(0);
					//      cout << "chievt="<<chievt<<endl; 
				}
				/// MBS Subvent:
				if((*setype==10) && (*sesubtype==1) && (*SubEventDataPtr==512))
				{
					//imp     cout << "Found a valid MBS subevent OKKKKKKKKKKKKKKKKKKKKKKKKK"<<endl;
					MBSTSh=*(SubEventDataPtr+1) & 0x0000ffff;
					MBSTSm=*(SubEventDataPtr+2) & 0x0000ffff;
					MBSTSs=*(SubEventDataPtr+3) & 0x0000ffff;
					mbsevt++;
					
					//imp     cout <<"MBS TIMESTAMP= "<<(*(SubEventDataPtr+1) & 0x0000ffff)<<" "<<(*(SubEventDataPtr+2)& 0x0000ffff)<<" "<<(*(SubEventDataPtr+3) & 0x0000ffff)<<endl;
					//     getchar();
					//      cout << "mbsevt="<<mbsevt<<endl; 
					
				}
				/// microBall Subvent:
				if((*setype==36) && (*sesubtype==2700))
				{
					//imp     cout << "Found a valid microBall subevent OKKKKKKKKKKKKKKKKKKKKKKKKK"<<endl;
					mball_evt->ProcessCurrentEventMBS(SubEventDataPtr,*sebuflenght);
					//imp     cout <<"MBALL raw multi="<<mball_evt->GetMBALLMultiplicity()<<endl; 
				}
				/// LAND Subvent:
				if((*setype==36) && (*sesubtype==9494) )
				{
					//imp     cout << "Found a valid LAND subevent OKKKKKKKKKKKKKKKKKKKKKKKKK"<<endl;
					land_evt->ProcessCurrentEventMBS(SubEventDataPtr,*sebuflenght);
					//imp     cout <<"LAND raw multi="<<land_evt->GetLANDMultiplicity()<<endl; 
					//     land_evt->GetLANDHitMaps();
					//     getchar(); 
				} 
				/// ToFWall Subvent:
				if((*setype==32) && (*sesubtype==1130) )
				{
					//imp     cout << "Found a valid ToFWall subevent OKKKKKKKKKKKKKKKKKKKKKKKKK"<<endl;
				} 
				/// SCALER Subvent:
				if((*setype==36) && (*sesubtype==2800) )
				{
					//imp     cout << "Found a valid SCALER subevent OKKKKKKKKKKKKKKKKKKKKKKKKK"<<endl;
					Trig_Patt=*SubEventDataPtr;
					//imp     cout <<"Trigger Pattern= "<<Trig_Patt<<endl;
					for (int k=0;k<16;k++){
						bit_Trig_Patt[k]=0;
					}
					if(Trig_Patt & 0x0001) bit_Trig_Patt[0] =1;
					if(Trig_Patt & 0x0002) bit_Trig_Patt[1] =1;
					if(Trig_Patt & 0x0004) bit_Trig_Patt[2] =1;
					if(Trig_Patt & 0x0008) bit_Trig_Patt[3] =1;
					if(Trig_Patt & 0x0010) bit_Trig_Patt[4] =1;
					if(Trig_Patt & 0x0020) bit_Trig_Patt[5] =1;
					if(Trig_Patt & 0x0040) bit_Trig_Patt[6] =1;
					if(Trig_Patt & 0x0080) bit_Trig_Patt[7] =1;
					if(Trig_Patt & 0x0100) bit_Trig_Patt[8] =1;
					if(Trig_Patt & 0x0200) bit_Trig_Patt[9] =1;
					if(Trig_Patt & 0x0400) bit_Trig_Patt[10]=1;
					for (k=0;k<16;k++){
						//imp      cout << "k= "<<k<<" "<<bit_Trig_Patt[k]<<endl;
					}
				}
				/// CAMAC Subvent:
				if((*setype==34) && (*sesubtype==1) )
				{
					//imp     cout << "Found a valid CAMAC subevent OKKKKKKKKKKKKKKKKKKKKKKKKK"<<endl;
					
					for (int ii=5; ii<13; ii++)
					{
						//		       word = 0xFFFF0313;
						//		       word = (word<<21)>>21;
						word[ii] = (((*(SubEventDataPtr+ii))<<21)>>21);
						//	    cout<<(hex)<<word[i]<<"\t";
					}
					//	   cout<<endl<<"TDC#2: ";
					for (ii=13; ii<21; ii++)
					{
						word[ii] = (((*(SubEventDataPtr+ii))<<21)>>21);
						//		       word = psubevt->Data(i);
						//		       word= (word<<21)>>21;
						//	       cout<<(hex)<<word[i]<<"\t";
					}
					//	   cout<<endl<<"QDC#1: ";
					for (ii=21; ii<33; ii++)
					{
						word[ii] = (((*(SubEventDataPtr+ii))<<20)>>20);
						//	       cout<<(hex)<<word[i]<<"\t";
					}
					stel_evt->X1raw  = word[13];
					stel_evt->X2raw  = word[14];
					stel_evt->Y1raw  = word[15];
					stel_evt->Y2raw  = word[16];
					stel_evt->TS1raw = word[5];
					stel_evt->TS2raw = word[6];
					stel_evt->QS1raw = word[21];
					stel_evt->QS2raw = word[22];
					stel_evt->Calibrate();
					//    stel_evt->FillHistos();
						
				}
				
				/// KRAKOW Subvent:
				if((*setype == 36) && (*sesubtype == 1) ) /// && (iKrakowEv < 400)
				{
					iKrakowEv++;
					
					Int_t sKrkMult = krakow_evt -> ProcessCurrentEventMBS ( SubEventDataPtr, *sebuflenght, KrattaMapping);
					hKrakow_Mult -> Fill ( sKrkMult );
          eKrattaRawEvent->SetId( i );          
          ///cout << eKrattaRawEvent.ToString() << "; id = " << eKrattaRawEvent.GetId() << endl;
          ///cout << eKrattaRawEvent.ToStringPhotodiode() << endl;
				}
			}/// END LOOP OVER SUBEVENTS for(Int_t j=1;j<=nrSubEvts;j++)
			//   cout << "inside 1"<<endl;
			i++;
			//   getchar();
		} /// END LOOP OVER EVENTS while(!lastevent)
		//  cout << "outside"<<endl;
		myMbs->CloseFile();
	} /// END CONDITION if(!(*status))
} /// END LOOP OVER FILES for(int irun=0;irun<Nruns;irun++)
chitree.Close();
inputfile.close();

}