slope.cxx

Go to the documentation of this file.

// simply add a slope to the calibrations constants to solve the positive/negative 
// unsymmetry problem 
//1=0.0; 2=-0.05; 3=-0.1; 4=-0.15
// high and low corresponding to momentum
//const double k_high = 0.07;
//const double k_low = 0.07;
const double k_high = 0.0;
const double k_low = 0.0;
const double b= 1.0;
const bool b_trunc = true; 
 
void slope()
{
    TFile *f=new TFile("eangle.root");
    ostringstream strout;
    strout<<"ddgcalib";
    TTree *tree = f->Get(strout.str().c_str()); 
    double dedx[100],eangle[100], fitmean[100], fitmeanerr[100];
    double pos_fitmean[100], pos_fitmeanerr[100], neg_fitmean[100], neg_fitmeanerr[100];
    //    tree->SetBranchAddress("gain",dedx);
    if(b_trunc) tree->SetBranchAddress("truncgain",dedx);
    else tree->SetBranchAddress("gain",dedx);
    tree->SetBranchAddress("eangle",eangle);
    tree->SetBranchAddress("fitmean",fitmean);
    tree->SetBranchAddress("fitmeanerr",fitmeanerr);
    tree->SetBranchAddress("pos_fitmean",pos_fitmean);
    tree->SetBranchAddress("pos_fitmeanerr",pos_fitmeanerr);
    tree->SetBranchAddress("neg_fitmean",neg_fitmean);
    tree->SetBranchAddress("neg_fitmeanerr",neg_fitmeanerr);
 
    cout<<"Entries: "<<tree->GetEntries()<<endl; 
    tree->GetEntry(0);
 
    int entry[1] = {100};
    double denangle[100];
    // for(int i=0; i<100; i++) denangle[i] = (dedx[i])*(k*(eangle[i])+b); // correct by a slope
    // next we try sub-region correction with a slope 
    for(int i=0; i<100; i++){
      if(fabs(eangle[i])<=0.04) denangle[i] = (dedx[i])*(k_high*(eangle[i])+b);
      //      cout << "eangle: " << eangle[i] << endl;
      else denangle[i] = (dedx[i])*(k_low*(eangle[i])+b); // correct by a slope
    }
    TFile *f_save = new TFile("eangle_new.root","recreate");
    TTree* entra = new TTree("entracalib", "entracalib");
    entra->Branch("1denangle_entry", entry, "entry[1]/I");
    entra->Branch("1denangle", denangle, "denangle[100]/D");
    entra->Fill();
    entra->Write();
    f_save->Close();
    f->Close();
}