dtotalEdepPercent.C

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file dtotalEdepPercent.C

void EdepPercentSmall (const char *calo="HO", const char *part="e")
{
  char infile[100];
  // int pz[9] = {1,2,4,8,12,16,32,40,50};
  // double x[9];
 char hname [100];
 sprintf(hname, "%s%sd",calo,part/*,pz[i]*/);
 TH1F *h1 = new TH1F(hname, hname, 100, 0, 10);
  for (int i=0; i<3; i++)
{
  //x[i] = pz[i];
  sprintf(infile, "/sphenix/user/jpinkenburg/ShowerSize/ntuple/%s_eta0_4GeV.root",part/*, pz[i]*/);
              char hname [100];
              // sprintf(hname, "%s%s12GeV",calo,part,pz[i]);
              //TH1F *h1 = new TH1F(hname, "", 100, 0, 1);
              TFile *f = TFile::Open(infile);
              gROOT->cd();
              TNtuple *nt = (TNtuple *)f->Get("de");
              char calid [100];
              sprintf(calid, "ID==%s", i); 
              nt->Project(hname,"edep",calid);
              // char cal [100];
              // sprintf(cal,"%s",calo);
              /*
              if (calo == "HO")
                {
                  nt->Project(hname,"edep","ID==2");
                }
              else if(calo == "HI")
                { 
                  nt->Project(hname,"edep","ID==1");
                } 
              else if(calo == "E")
                { 
                  nt->Project(hname,"edep","ID==0");
                }
              else if (calo == "MAG") 
                {
                  nt->Project(hname,"(MAG)/(EA+ES+HIA+HIS+HOA+HOS+BH+MAG)");
                }
              else if (calo == "BH")
                {
          nt->Project(hname,"(BH)/(EA+ES+HIA+HIS+HOA+HOS+BH+MAG)");
                }
              else 
                {
                  cout << "ERROR" << endl;
                  }*/
              f->Close();
              h1->Draw();
              char fname [100];
              sprintf(fname, "%sShowerPercent.root",part);
              TFile *fout = TFile::Open(fname,"UPDATE");
              h1->Write();
              fout->Write();
              cout << h1->Integral() << " "  <<  hname << " " /*<< pz[i] <<*/ endl;
              //fout->Write();
              fout->Close();
  }
}
void EdepSmall(const char *calo="HOe", const char *part = "e")
{
   int pz[9] = {1,2,4,8,12,16,32,40,50};
   double x[9];
   double y[9];
   double xer[9] = {0};
   double yer[9];
   char hname[100];
   char fname[100];
   sprintf(fname,"%spercent.root",part);
   char area[100];
   sprintf(area,"%s",calo);
  TFile *fin = TFile::Open(fname);
  gROOT->cd();
   for (int i=0; i<9; i++)
     {
       x[i] = pz[i];
       sprintf(hname,"%s%s%dd",calo,part,pz[i]);
TH1 *h1 = (TH1F *) fin->Get(hname);
     y[i] = h1->GetMean(1);
   yer[i] = h1->GetRMS(1);
     }
   gr = new TGraphErrors(9,x,y,xer,yer);
   gr->SetTitle(area);
   gr->Draw("ALP");
   gr->SetMarkerStyle(20);
   gr->SetMarkerColor(2);
   fin->Close();
}