commenting out use of Oysteins ttree so code will run without it since it isnt used...

diff --git a/PWGLF/totEt/macros/emEt/calculateCorrections.C b/PWGLF/totEt/macros/emEt/calculateCorrections.C
index 1cf4fa84425f17063a3436c7a79c168325de1953..1cbb2319b094fb2cbcff8b338ba06ec223d0499c 100644 (file)
--- a/PWGLF/totEt/macros/emEt/calculateCorrections.C
+++ b/PWGLF/totEt/macros/emEt/calculateCorrections.C
@@ -59,10 +59,10 @@ int calculateCorrections(TString filename="rootFiles/LHC11a10a_bis/Et.ESD.simPbP
   
   TList *l = (TList*)f->Get("out1");
   
-  TTree *primaryTree = (TTree*)l->FindObject(("fPrimaryTree"+detector+"MC").Data());
-  TTree *recTree = (TTree*)l->FindObject(("fEventSummaryTree"+detector+"Rec").Data());
-  TTree *mcTree = (TTree*)l->FindObject(("fEventSummaryTree"+detector+"MC").Data());
-  std::cout << primaryTree << " " << recTree << " " << mcTree << std::endl;
+//   TTree *primaryTree = (TTree*)l->FindObject(("fPrimaryTree"+detector+"MC").Data());
+//   TTree *recTree = (TTree*)l->FindObject(("fEventSummaryTree"+detector+"Rec").Data());
+//   TTree *mcTree = (TTree*)l->FindObject(("fEventSummaryTree"+detector+"MC").Data());
+//   std::cout << primaryTree << " " << recTree << " " << mcTree << std::endl;
  
   Int_t clusterMult = 0;
   Int_t nChargedNotRemoved = 0;
@@ -73,11 +73,11 @@ int calculateCorrections(TString filename="rootFiles/LHC11a10a_bis/Et.ESD.simPbP
   Double_t emEtRec = 0;
   Double_t emEtMc = 0;
 
-  recTree->SetBranchAddress("fNeutralMultiplicity", &clusterMult);
-  mcTree->SetBranchAddress("fChargedNotRemoved", &nChargedNotRemoved);
-  mcTree->SetBranchAddress("fNeutralNotRemoved", &nNeutralNotRemoved);
-  mcTree->SetBranchAddress("fGammaRemoved", &nGammaRemoved);
-  mcTree->SetBranchAddress("fSecondaryNotRemoved", &nSecNotRemoved);
+//   recTree->SetBranchAddress("fNeutralMultiplicity", &clusterMult);
+//   mcTree->SetBranchAddress("fChargedNotRemoved", &nChargedNotRemoved);
+//   mcTree->SetBranchAddress("fNeutralNotRemoved", &nNeutralNotRemoved);
+//   mcTree->SetBranchAddress("fGammaRemoved", &nGammaRemoved);
+//   mcTree->SetBranchAddress("fSecondaryNotRemoved", &nSecNotRemoved);
   
   Float_t maxMult = 99.5;
   Int_t nbins = 100;
@@ -91,16 +91,16 @@ int calculateCorrections(TString filename="rootFiles/LHC11a10a_bis/Et.ESD.simPbP
   TH2I *hGammaVsClusterMult = new TH2I("hGammaVsMult", "hGammaVsMult", nbins, -0.5, maxMult, nbins, -0.5, maxMult);
   TH2I *hSecVsClusterMult = new TH2I("hSecVsMult", "hSecVsMult", nbins, -0.5, maxMult, nbins, -0.5, maxMult);
 
-  Int_t nEvents = mcTree->GetEntriesFast();
-    for(Int_t i = 0; i < nEvents; i++)
-    {
-      mcTree->GetEvent(i);
-      recTree->GetEvent(i);
-      hChargedVsClusterMult->Fill(clusterMult, nChargedNotRemoved);
-      hNeutralVsClusterMult->Fill(clusterMult, nNeutralNotRemoved);
-      hGammaVsClusterMult->Fill(clusterMult, nGammaRemoved);
-      hSecVsClusterMult->Fill(clusterMult, nSecNotRemoved);
-    }
+//   Int_t nEvents = mcTree->GetEntriesFast();
+//     for(Int_t i = 0; i < nEvents; i++)
+//     {
+//       mcTree->GetEvent(i);
+//       recTree->GetEvent(i);
+//       hChargedVsClusterMult->Fill(clusterMult, nChargedNotRemoved);
+//       hNeutralVsClusterMult->Fill(clusterMult, nNeutralNotRemoved);
+//       hGammaVsClusterMult->Fill(clusterMult, nGammaRemoved);
+//       hSecVsClusterMult->Fill(clusterMult, nSecNotRemoved);
+//     }
     
   c1->Divide(2,2);
   c1->cd(1);
@@ -118,8 +118,8 @@ int calculateCorrections(TString filename="rootFiles/LHC11a10a_bis/Et.ESD.simPbP
   //TF1 fitcharged("fitcharged","([0] + [1]*x)*(0.48/([2] + [3]*[2]))", 0, 100);
   TF1 fitcharged("fitcharged","pol2", 0, 100);//fit of number of charged tracks vs detector multiplicity
   //if straight line track matching roughly not dependent on centrality
-//   fitcharged.FixParameter(2, p0);
-//   fitcharged.FixParameter(3, p1);
+  //fitcharged.FixParameter(2, p0);
+   //   fitcharged.FixParameter(3, p1);
   TFitResultPtr chRes = chProf->Fit(&fitcharged,"S");
   TArrayD ch;
   if(!chRes)