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;
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;
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);
//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)