From: cnattras Date: Fri, 4 Jul 2014 19:07:21 +0000 (-0400) Subject: commenting out use of Oysteins ttree so code will run without it since it isnt used... X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=da90ef374ab54d88521e948d4823602cf6d1aa90;p=u%2Fmrichter%2FAliRoot.git commenting out use of Oysteins ttree so code will run without it since it isnt used anyways --- diff --git a/PWGLF/totEt/macros/emEt/calculateCorrections.C b/PWGLF/totEt/macros/emEt/calculateCorrections.C index 1cf4fa84425..1cbb2319b09 100644 --- 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)