StEEmcPool/macros/LCP/compSelection.C

const int mxGr=3;
TGraphErrors * grL[mxGr];

const int  nObs=3;
char *obsL[nObs]={ "A_n", "A_s", "A_d"};


compSelection() {
  gStyle->SetOptStat(0);
  gStyle->SetOptFit(0);

  TString  wrkDir="final/";
  const int  nSel=11;
  char *selL[nSel]={ "defaultB-H/", "default-H/", "maxEta1.4-H/", "highPol-H/", "nPrim5_20-H/","zVert50-H/", "pT1_3-H/", "posCharge-H/","shiftPatt-H/","swapPatt-H/","randPatt-H/" };

  createTgr();

  int isel;
  for (isel=0;isel<nSel;isel++) {
    TString  fname=wrkDir+selL[isel]+"endVer1.hist.root";
    // printf("ii=%d %s %s\n",isel, selL[isel], fname.Data());
    TFile *inpH=new TFile(fname);
    assert(inpH->IsOpen());
    inpH->ls();
    int io;
    for (io=0;io<nObs;io++) { 
      char name[100];
      sprintf(name,"%s*EtaBF",obsL[io]);
      TGraphErrors *gr =(TGraphErrors *) inpH->Get(name);
      assert(gr);
      gr->Print();
      sprintf(name,"avr%s*EtaBF",obsL[io]);
      TF1 *ff=gr->GetFunction(name);
      if(ff==0) continue; // no fit was made
      float val=ff->GetParameter(0);
      float err=ff->GetParError(0);
      printf("pol0= %f +/- %f\n",val,err);
      TGraphErrors * grOut=grL[io];
      int n=grOut->GetN();
      float x=isel+1;
      grOut->SetPoint(n,x,val);
      grOut->SetPointError(n,0,err); 
    }
  }

  grL[0]->Print();
    

  c=new TCanvas("aa","aa"   ,800,800);
  c->Divide(3,1);
  TLine *ln0=new TLine(0,0.,12,0.);
  
  int it; 
  for (it=0;it<mxGr;it++) {
    TGraphErrors * gr=grL[it];
    int n=gr->GetN();
    printf("\nit=%d  name='%s', N=%d\n",it,gr->GetName(),n);
    // gr->Print();
    if(n<=0) continue;
    c->cd(1+it);
    gr->Draw("AP");
    // gr->Fit("pol0");

    TAxis *ax=gr->GetXaxis();
    ax->SetTitle("choice of event selection");
    ln0->Draw();
    ln0->SetLineStyle(3);
  }

  // save Tgraph
  TString fname=wrkDir+"compCuts.hist.root";
  TFile *outH=new TFile(fname,"RECREATE");
  assert(outH->IsOpen());
  printf("save outH -->%s\n", fname.Data());
  outH->cd();
  for (it=0;it<mxGr;it++) {
    grL[it]->Write();
  }
  
  outH->ls();
  outH->Write();
}


//==========================================
//==========================================
void createTgr() {
  char *obsT[nObs]={ "An", "A#Sigma", "A#Delta"};
  // int  ampCol[mxAmp]={kBlue, kYellow, kGreen, kMagenta};

  int io;
  for (io=0;io<nObs;io++) { 
    char name[100];
    sprintf(name,"%s",obsL[io]);
    //printf("ic=%d iam=%d name=%s=\n",ic,iam,name);
    TGraphErrors *gr =new  TGraphErrors;
    gr->SetName(name);
    gr->SetTitle(obsT[io]);
    //gr->SetMarkerColor(ampCol[iam]);
    gr->SetMarkerSize(0.8);
    gr->SetMarkerStyle(21);
    grL[io]=gr;   
  }
}

---