EXP1803

#include <string.h>

#include <iostream.h>

#include <fstream.h>

#include "TSystem.h"

#include "TFile.h"

#include "TTree.h"

//#include "/media/users_NAS/Muzalevsky/exp1803/go4/TNeEvent.h"

using namespace std;

void fill1mm(){ 

        string         input_file_name_cs = "/media/users_NAS/Muzalevsky/exp1803/calSi/si_20_03_0002.root";                

        Int_t         max_channel = 300;

        TTree*                 t;

        UShort_t    NeEvent_SQ20[16];

        UShort_t    NeEvent_SQX_L[32];

        TBranch    *b_NeEvent_SQ20;

        TBranch    *b_NeEvent_SQX_L;

  //TNeEvent *revent = new TNeEvent();

        Float_t parSQ1[32], parSQ2[32],par201[16],par202[16];

  TString line1;

  ifstream myfile1;

  Int_t count=0;

  myfile1.open("/home/muzalevsky/AculUti/exp1803/cali1mm/calDL/SQX_L.cal");

  while (! myfile1.eof() ){

                line1.ReadLine(myfile1);

    //cout << line1.Data() << endl;

   // if(count<3) continue;

                if(line1.IsNull()) break;

    sscanf(line1.Data(),"%g %g", parSQ1+count,parSQ2+count);

        count++;

  }

  cout << " pars for 1mm " << endl;

  for(Int_t i=0;i<32;i++){

    cout << parSQ1[i] << " " << parSQ2[i] << endl;   

  }

  ifstream myfile2;

  TString line2;

  count=0;

  myfile2.open("/home/muzalevsky/AculUti/macroexp18/si20cal/SQ20_58.cal");

  while (! myfile2.eof() ){

                line2.ReadLine(myfile2);

    cout << line2.Data() << endl;

   // if(count<3) continue;

                if(line2.IsNull()) break;

    sscanf(line2.Data(),"%g %g", par201+count,par202+count);

    count++;

  }

  cout << " pars for 20 mkm " << endl;

  for(Int_t i=0;i<16;i++){

    cout << par201[i] << " " << par202[i] << endl;   

  }

        f = new TFile(input_file_name_cs.c_str());

        f->GetObject("AnalysisxTree",t);

        t->SetMakeClass(1);

        t->SetBranchAddress("NeEvent.SQX_L[32]", NeEvent_SQX_L, &b_NeEvent_SQX_L);

  t->SetBranchAddress("NeEvent.SQ20[16]", NeEvent_SQ20, &b_NeEvent_SQ20);

        const Long64_t nentries1 = t->GetEntries();

  cout << nentries1 << endl;

        Float_t SQX_L[32],SQ20[16];

        TFile *fw = new TFile("cal_20.root", "RECREATE");

        TTree *tw = new TTree("cal", "20mk spec in energies");

        tw->Branch("SQX_L",&SQX_L,"SQX_L[32]/F");

        tw->Branch("SQ20",&SQ20,"SQ20[16]/F");

  Int_t maxE = nentries1;

  cout<<">>> filling TREE up to "<<maxE<< " event"<<endl;

        for (Long64_t jentry=0; jentry<maxE;jentry++) {

                t->GetEntry(jentry);

                for(Int_t i=0; i<32; i++){

      SQX_L[i]=0.;// обнуление

      if(i<16) {

        SQ20[i]=0.;

      }

      if(NeEvent_SQX_L[i]>50.){ //cut the pedestal, calculating

        SQX_L[i] = NeEvent_SQX_L[i]*parSQ2[i] + parSQ1[i];

      } 

      if(i<14){

        if(NeEvent_SQ20[i]>200.) { //cut the pedestal, calculating

          SQ20[i] = NeEvent_SQ20[i]*par202[i] + par201[i];

        }

      }

          } // i<32

                tw->Fill();                        

        }//entries

        //htest->Draw();

        tw->Write();

        fw->Close();

}