getCurConstPar_dEdXCal.cxx

Go to the documentation of this file.

//void getPar(string filename)
{
    //TFile* f=new TFile(filename.c_str());
    TFile* f=new TFile("../CurConstRoot/DedxCurConst_data_minrun-maxrun_bossversion.root");
    TTree *curve = (TTree*) f->Get("curve");
    TTree *sigma = (TTree*) f->Get("sigma");
    
    double Curve[100];
    //double* Curve;
    double Sigma[100];
    int CurveSize,SigmaSize;
    
    curve->SetBranchAddress("CurveSize",&CurveSize);
    curve->SetBranchAddress("curve",Curve);
    sigma->SetBranchAddress("SigmaSize",&SigmaSize);
    sigma->SetBranchAddress("sigma",Sigma);
    
   // cout<<"AS!!!!!"<<endl;
    
    //int n = curve->GetEntries();
    ofstream file("test_data.dat");
    curve->GetEntry(0);
    
    //cout<<CurveSize<<"   sucess!"<<endl;
   // cout<<curve->GetEntries()<<"   sucess!"<<endl;
    //double cur;
    for(int i=0;i<CurveSize;i++)
    {
       // curve->GetEntry(i);
       // file1<<i<<"  "<<CurveSize<<endl;
//        cout<<Curve[i]<<endl;
       // cur = Curve[i];
       file<<i<<"  "<<Curve[i]<<endl; 
    }
    //file1.close();
    
     //cout<<"AS!!!!!"<<endl;
    
    //int m = sigma->GetEntries();
   // ofstream file2("/ihepbatch/besd26/anff/Calib/sigma_jpsi.dat");
    sigma->GetEntry(0);
    for(int i=0;i<SigmaSize;i++)
    {
        //sigma->GetEntry(i);
        file<<i<<"  "<<Sigma[i]<<endl;
       // file2<<SigmaSize<<"  "<<Sigma<<endl;
    }
   //  cout<<"AS!!!!!"<<endl;
    
    file.close();
}