CurConst_dEdXCal.cxx

Go to the documentation of this file.

 
{
    
  TFile f1("../CurConstRoot/DedxCurConst_data_minrun-maxrun_bossversion.root","recreate"); 
  TTree *curve=new TTree("curve","curve parameters");
  TTree* sigma=new TTree("sigma","sigma parameters");
 
  const int size1=17;
  const int size2=15;
 
  Double_t N_DedxCurve_Parameter[size1];
  Double_t N_DedxSigma_Parameter[size2];
  ifstream filea("../CurConstRoot/DedxCurConst_data_minrun-maxrun_bossversion.dat",ifstream::in); 
  for(int i=0;i<size1;i++)
  {
      filea>>N_DedxCurve_Parameter[i];
      cout<<N_DedxCurve_Parameter[i]<<endl;
  }
  for(int j=0;j<size2;j++)
  {
      filea>>N_DedxSigma_Parameter[j];
      cout<<N_DedxSigma_Parameter[j]<<endl;
  }
 
 
  int CurveSize=sizeof(N_DedxCurve_Parameter)/sizeof(N_DedxCurve_Parameter[0]);
  int SigmaSize=sizeof(N_DedxSigma_Parameter)/sizeof(N_DedxSigma_Parameter[0]);
 
 
 curve->Branch("CurveSize",&CurveSize,"CurveSize/I");
 curve->Branch("curve",N_DedxCurve_Parameter,"N_DedxCurve_Parameter[CurveSize]/D");
 curve->Fill();
 curve->Write();
 sigma->Branch("SigmaSize",&SigmaSize,"SigmaSize/I");
 sigma->Branch("sigma",N_DedxSigma_Parameter,"N_DedxSigma_Parameter[SigmaSize]/D");
 sigma->Fill();
 sigma->Write();
 f1.Close();
 cout<<"end"<<endl;
 }