- h->Fill(h->GetBinCenter(bin),1);
-}
-
-
-//____________________________________________________________________
-
-
-void AliAnalysisTaskJetChem::FillPIDhisto(TH1F* hist, Int_t pdg, Float_t weight){
-
- // convert pdg code to Geant ID and fill corresponding bin
-
- Int_t fGID = fpdgdb->ConvertPdgToGeant3(pdg);
-
- //cout<<" pdg "<<pdg<<" fGID "<<fGID<<endl;
-
- if(TMath::Abs(pdg) == 311) fGID = kPDGpm311Bin;
- if(pdg == 333) fGID = kPDG333Bin;
- if(TMath::Abs(pdg) == 313) fGID = kPDGpm313Bin;
- if(pdg == 323) fGID = kPDGp323Bin;
- if(pdg == -323) fGID = kPDGm323Bin;
-
- if(TMath::Abs(pdg)==12 || TMath::Abs(pdg)==14 || TMath::Abs(pdg)==16) fGID = kPDGNeutrinoBin;
-
- if(TMath::Abs(pdg)==4122) fGID = kPDGCharmedBaryonBin;
-
- if(1<=TMath::Abs(pdg) && TMath::Abs(pdg)<=6) fGID = kPDGQuarkBin;
-
- if(TMath::Abs(pdg)==1103 || TMath::Abs(pdg)==2101 || TMath::Abs(pdg)==2103 || TMath::Abs(pdg)==2203 ||
- TMath::Abs(pdg)==3101 || TMath::Abs(pdg)==3103 || TMath::Abs(pdg)==3201 || TMath::Abs(pdg)==3203 ||
- TMath::Abs(pdg)==3303 || TMath::Abs(pdg)==4101 || TMath::Abs(pdg)==4103 || TMath::Abs(pdg)==4201 ||
- TMath::Abs(pdg)==4203 || TMath::Abs(pdg)==4301 || TMath::Abs(pdg)==4303 || TMath::Abs(pdg)==4403 ||
- TMath::Abs(pdg)==5101 || TMath::Abs(pdg)==5103 || TMath::Abs(pdg)==5201 || TMath::Abs(pdg)==5203 ||
- TMath::Abs(pdg)==5301 || TMath::Abs(pdg)==5303 || TMath::Abs(pdg)==5401 || TMath::Abs(pdg)==5403 ||
- TMath::Abs(pdg)==5503) fGID = kPDGDiQuarkBin;
-
-
- hist->Fill(fGID,weight);
-
- if(fGID>hist->GetBinLowEdge(hist->GetNbinsX()+1) ||
- fGID<hist->GetBinLowEdge(1)){
-
- AliError(Form("fGID %d for pdg %d exceeding histo limits ",fGID,pdg));
- }
-
- if(fGID == 0){
- AliInfo(Form("fGID 0 for pdg %d ",pdg));
- }
-
-}
-
-//____________________________________________________________________
-void AliAnalysisTaskJetChem::CreateHistos(){
-
- // make histos
-
- fListOfHistos = new TList();
- fListOfHistos->SetOwner(kTRUE);
-
- fhPrimVertexNCont = new TH1F("hPrimVertexNCont","",52,-2,50);
- fhPrimVertexNCont->SetXTitle("");
- fhPrimVertexNCont->Sumw2();
- fListOfHistos->Add( fhPrimVertexNCont );
-
- fhPrimVertexRho = new TH1F("hPrimVertexRho","",100,0,1);
- fhPrimVertexRho->SetXTitle("");
- fhPrimVertexRho->Sumw2();
- fListOfHistos->Add( fhPrimVertexRho );
-
- fhPrimVertexZ = new TH1F("hPrimVertexZ","",40,-20,20);
- fhPrimVertexZ->SetXTitle("");
- fhPrimVertexZ->Sumw2();
- fListOfHistos->Add( fhPrimVertexZ );
-
- fhNJets = new TH1F("hNJets","",10, 0, 10);
- fhNJets->SetXTitle("# of jets");
- fhNJets->Sumw2();
- fListOfHistos->Add( fhNJets );
-
- fhNJetsMC = new TH1F("hNJetsMC","",10,0,10);
- fhNJetsMC->SetXTitle("# of jets");
- fhNJetsMC->Sumw2();
- fListOfHistos->Add( fhNJetsMC );
-
- fhLeadingEta = new TH1F("hLeadingEta","Leading Jet eta",12,-0.6,0.6);
- fhLeadingEta->SetXTitle("eta");
- fhLeadingEta->SetYTitle("dN/deta");
- fhLeadingEta->Sumw2();
- fListOfHistos->Add(fhLeadingEta);
-
- fhLeadingNTracksVsEta = new TH2F("hLeadingNTracksVsEta","",20,-1.0,1.0,20,0,20);
- fhLeadingNTracksVsEta->SetXTitle("eta");
- fhLeadingNTracksVsEta->SetYTitle("# of tracks");
- fhLeadingNTracksVsEta->Sumw2();
- fListOfHistos->Add(fhLeadingNTracksVsEta);
-
- fhLeadingPtVsEta = new TH2F("hLeadingPtVsEta","",20,-1.0,1.0,50,0,50);
- fhLeadingPtVsEta->SetXTitle("eta");
- fhLeadingPtVsEta->SetYTitle("# of tracks");
- fhLeadingPtVsEta->Sumw2();
- fListOfHistos->Add(fhLeadingPtVsEta);
-
-
- fhLeadingPhi = new TH1F("hLeadingPhi","Leading Jet phi",63,0,6.3);
- fhLeadingPhi->SetXTitle("phi");
- fhLeadingPhi->SetYTitle("dN/dphi");
- fhLeadingPhi->Sumw2();
- fListOfHistos->Add(fhLeadingPhi);
-
-
- fhLeadingPt = new TH1F("hLeadingPt","leading Jet p_{T}",50,0,50);
- fhLeadingPt->SetXTitle("p_{T} (GeV/c)");
- fhLeadingPt->SetYTitle("dN/dp_{T} (1/GeV)");
- fhLeadingPt->Sumw2();
- fListOfHistos->Add( fhLeadingPt );
-
- fhLeadingPtDiffr = new TH1F("hLeadingPtDiffr","leading Jet p_{T}",50,0,50);
- fhLeadingPtDiffr->SetXTitle("P_{T} (GeV/c)");
- fhLeadingPtDiffr->SetYTitle("dN/dp_{T} (1/GeV)");
- fhLeadingPtDiffr->Sumw2();
- fListOfHistos->Add( fhLeadingPtDiffr );
-
- fhLeadingEtaMC = new TH1F("hLeadingEtaMC","Leading Jet eta",12,-0.6,0.6);
- fhLeadingEtaMC->SetXTitle("eta");
- fhLeadingEtaMC->SetYTitle("dN/deta");
- fhLeadingEtaMC->Sumw2();
- fListOfHistos->Add(fhLeadingEtaMC);
-
- fhLeadingPhiMC = new TH1F("hLeadingPhiMC","Leading Jet phi",63,0,6.3);
- fhLeadingPhiMC->SetXTitle("phi");
- fhLeadingPhiMC->SetYTitle("dN/dphi");
- fhLeadingPhiMC->Sumw2();
- fListOfHistos->Add(fhLeadingPhiMC);
-
- fhLeadingPtMC = new TH1F("hLeadingPtMC","Leading Jet p_{T}",50,0,50);
- fhLeadingPtMC->SetXTitle("p_{T} (GeV/c)");
- fhLeadingPtMC->SetYTitle("dN/dp_{T} (1/GeV)");
- fhLeadingPtMC->Sumw2();
- fListOfHistos->Add( fhLeadingPtMC );
-
- fhLeadingPtMCDiffr = new TH1F("hLeadingPtMCDiffr","Leading Jet p_{T}",50,0,50);
- fhLeadingPtMCDiffr->SetXTitle("p_{T} (GeV/c)");
- fhLeadingPtMCDiffr->SetYTitle("dN/dp_{T} (1/GeV)");
- fhLeadingPtMCDiffr->Sumw2();
- fListOfHistos->Add( fhLeadingPtMCDiffr );
-
- fhPhiEtaTracksNoCut = new TH2F("hPhiEtaTracksNoCut","phi vs eta tracks",20,-1.0,1.0,63,0,6.3);
- fhPhiEtaTracksNoCut->SetXTitle("eta");
- fhPhiEtaTracksNoCut->SetYTitle("phi");
- fhPhiEtaTracksNoCut->Sumw2();
- fListOfHistos->Add(fhPhiEtaTracksNoCut);
-
- fhPtTracksNoCut = new TH1F("hPtTracksNoCut","p_{T} tracks",150,0,150);
- fhPtTracksNoCut->SetXTitle("p_{T} (GeV)");
- fhPtTracksNoCut->SetYTitle("dN/dp_{T} (1/GeV)");
- fhPtTracksNoCut->Sumw2();
- fListOfHistos->Add(fhPtTracksNoCut);
-
- fhPhiEtaTracks = new TH2F("hPhiEtaTracks","phi vs eta tracks",20,-1.0,1.0,63,0,6.3);
- fhPhiEtaTracks->SetXTitle("eta");
- fhPhiEtaTracks->SetYTitle("phi");
- fhPhiEtaTracks->Sumw2();
- fListOfHistos->Add(fhPhiEtaTracks);
-
- fhPtTracks = new TH1F("hPtTracks","p_{T} tracks",150,0,150);
- fhPtTracks->SetXTitle("P{T} (GeV)");
- fhPtTracks->SetYTitle("dN/dp_{T} (1/GeV)");
- fhPtTracks->Sumw2();
- fListOfHistos->Add(fhPtTracks);
-
- fhTrackMult = new TH1F("hTrackMult","",150,0,150);
- fhTrackMult->SetXTitle("n Tracks");
- fhTrackMult->SetYTitle("counts");
- fhTrackMult->Sumw2();
- fListOfHistos->Add(fhTrackMult);
-
- fhEtaMCTracks = new TH1F("hEtaMCTracks","eta tracks",30,-1.5,1.5);
- fhEtaMCTracks->SetXTitle("eta");
- fhEtaMCTracks->SetYTitle("dN/deta");
- fhEtaMCTracks->Sumw2();
- fListOfHistos->Add(fhEtaMCTracks);
-
- fhPhiMCTracks = new TH1F("hPhiMCTracks","phi tracks",63,0,6.3);
- fhPhiMCTracks->SetXTitle("phi");
- fhPhiMCTracks->SetYTitle("dN/dphi");
- fhPhiMCTracks->Sumw2();
- fListOfHistos->Add(fhPhiMCTracks);
-
- fhPtMCTracks = new TH1F("hPtMCTracks","p_{T} tracks",50,0,50);
- fhPtMCTracks->SetXTitle("p_{T} (GeV)");
- fhPtMCTracks->SetYTitle("dN/dp_{T} (1/GeV)");
- fhPtMCTracks->Sumw2();
- fListOfHistos->Add(fhPtMCTracks);
-
- fhnTracksVsPtLeading = new TH2F("hnTracksVsPtLeading","",50,0.,50.,20,-0.5,19.5);
- fhnTracksVsPtLeading->SetXTitle("p_{T} (GeV/c)");
- fhnTracksVsPtLeading->SetYTitle("n tracks");
- fhnTracksVsPtLeading->Sumw2();
- fListOfHistos->Add( fhnTracksVsPtLeading );
-
- fhdNdEtaPhiDist = new TH1F("hdNdEtaPhiDist","Charge particle density |#eta|< 1 vs #Delta#phi", 120, 0., 2.*TMath::Pi());
- fhdNdEtaPhiDist->SetXTitle("#Delta#phi");
- fhdNdEtaPhiDist->SetYTitle("dN{ch}/d#etad#phi");
- fhdNdEtaPhiDist->Sumw2();
- fListOfHistos->Add( fhdNdEtaPhiDist );
-
- fhRegionSumPtMaxVsEt = new TH1F("hRegionSumPtMaxVsEt","P{T}^{90, max} vs Leading Jet P{T}",50,0.,50.);
- fhRegionSumPtMaxVsEt->SetXTitle("P{T} (GeV/c)");
- fhRegionSumPtMaxVsEt->Sumw2();
- fListOfHistos->Add( fhRegionSumPtMaxVsEt );
-
- fhRegionMultMaxVsEt = new TH1F("hRegionMultMaxVsEt","N{ch}^{90, max} vs Leading Jet P{T}",50,0.,50.);
- fhRegionMultMaxVsEt->SetXTitle("E (GeV hRegionAveSumPtVsEt/c)");
- fhRegionMultMaxVsEt->Sumw2();
- fListOfHistos->Add( fhRegionMultMaxVsEt );
-
- fhRegionSumPtMinVsEt = new TH1F("hRegionSumPtMinVsEt","P{T}^{90, min} vs Leading Jet P{T}",50,0.,50.);
- fhRegionSumPtMinVsEt->SetXTitle("P{T} (GeV/c)");
- fhRegionSumPtMinVsEt->Sumw2();
- fListOfHistos->Add( fhRegionSumPtMinVsEt );
-
- fhRegionMultMinVsEt = new TH1F("hRegionMultMinVsEt","N{ch}^{90, min} vs Leading Jet P{T}",50,0.,50.);
- fhRegionMultMinVsEt->SetXTitle("E (GeV/c)");
- fhRegionMultMinVsEt->Sumw2();
- fListOfHistos->Add( fhRegionMultMinVsEt );
-
- // V0s
-
- fhNV0s = new TH1F("hNV0s","n V0s",50,0,50);
- fhNV0s->SetXTitle("n V0s");
- fhNV0s->Sumw2();
- fListOfHistos->Add(fhNV0s);
-
- fhV0onFly = new TH1F("hV0onFly","on-the-fly V0",5,0,5);
- fhV0onFly->SetXTitle("is on-the-fly V0");
- fhV0onFly->Sumw2();
- fListOfHistos->Add(fhV0onFly);
-
- fhV0DCADaughters = new TH1F("hV0DCADaughters","V0 DCA daughters",200,0,2.0);
- fhV0DCADaughters->SetXTitle("V0 DCA daughters");
- fhV0DCADaughters->Sumw2();
- fListOfHistos->Add(fhV0DCADaughters);
-
- fhV0Radius = new TH1F("hV0Radius","V0 radius",2500,0,250);
- fhV0Radius->SetXTitle("V0 radius");
- fhV0Radius->Sumw2();
- fListOfHistos->Add(fhV0Radius);
-
- fhV0DCAToVertex = new TH1F("hV0DCAToVertex","",100,0,10);
- fhV0DCAToVertex->SetXTitle("V0 DCA (cm)");
- fhV0DCAToVertex->Sumw2();
- fListOfHistos->Add(fhV0DCAToVertex);
-
- fhV0DCAToVertexK0 = new TH1F("hV0DCAToVertexK0","",100,0,10);
- fhV0DCAToVertexK0->SetXTitle("V0 DCA (cm)");
- fhV0DCAToVertexK0->Sumw2();
- fListOfHistos->Add(fhV0DCAToVertexK0);
-
- fhV0InvMassK0 = new TH1F("hV0InvMassK0","",2000,0,2);
- fhV0InvMassK0->SetXTitle("inv mass (GeV)");
- fhV0InvMassK0->Sumw2();
- fListOfHistos->Add(fhV0InvMassK0);
-
- fhV0PtVsInvMassK0 = new TH2F("hV0PtVsInvMassK0","",100,0.45,0.55,100,0.,50.);
- fhV0PtVsInvMassK0->SetXTitle("inv mass (GeV)");
- fhV0PtVsInvMassK0->SetYTitle("p_{T} (GeV)");
- fhV0PtVsInvMassK0->Sumw2();
- fListOfHistos->Add(fhV0PtVsInvMassK0);
-
- fhV0InvMassK0JetEvt = new TH1F("hV0InvMassK0JetEvt","",2000,0,2);
- fhV0InvMassK0JetEvt->SetXTitle("inv mass (GeV)");
- fhV0InvMassK0JetEvt->Sumw2();
- fListOfHistos->Add(fhV0InvMassK0JetEvt);
-
- fhV0InvMassLambda = new TH1F("hV0InvMassLambda","",2000,0,2);
- fhV0InvMassLambda->SetXTitle("inv mass (GeV)");
- fhV0InvMassLambda->Sumw2();
- fListOfHistos->Add(fhV0InvMassLambda);
-
- fhV0InvMassAntiLambda = new TH1F("hV0InvMassAntiLambda","",2000,0,2);
- fhV0InvMassAntiLambda->SetXTitle("inv mass (GeV)");
- fhV0InvMassAntiLambda->Sumw2();
- fListOfHistos->Add(fhV0InvMassAntiLambda);
-
- fhV0InvMassLambdaJetEvt = new TH1F("hV0InvMassLambdaJetEvt","",2000,0,2);
- fhV0InvMassLambdaJetEvt->SetXTitle("inv mass (GeV)");
- fhV0InvMassLambdaJetEvt->Sumw2();
- fListOfHistos->Add(fhV0InvMassLambdaJetEvt);
-
- fhV0InvMassAntiLambdaJetEvt = new TH1F("hV0InvMassAntiLambdaJetEvt","",2000,0,2);
- fhV0InvMassAntiLambdaJetEvt->SetXTitle("inv mass (GeV)");
- fhV0InvMassAntiLambdaJetEvt->Sumw2();
- fListOfHistos->Add(fhV0InvMassAntiLambdaJetEvt);
-
- fhdROpanK0VsPt = new TH2F("hdROpanK0VsPt","V0 dR vs pt",100,0,10,100,0,10);
- fhdROpanK0VsPt->SetXTitle("opening angle R (rad)");
- fhdROpanK0VsPt->Sumw2();
- fListOfHistos->Add(fhdROpanK0VsPt);
-
- fhdPhiJetV0 = new TH1F("hdPhiJetV0","",640,-3.2,3.2);
- fhdPhiJetV0->SetXTitle("#Delta #phi V0 - jet");
- fhdPhiJetV0->Sumw2();
- fListOfHistos->Add(fhdPhiJetV0);
-
- fhdPhiJetK0 = new TH1F("hdPhiJetK0","",640,-3.2,3.2);
- fhdPhiJetK0->SetXTitle("#Delta #phi V0 - jet");
- fhdPhiJetK0->Sumw2();
- fListOfHistos->Add(fhdPhiJetK0);
-
- fhdRJetK0 = new TH1F("hdRJetK0","dN/dR K0-jet",500,0,5);
- fhdRJetK0->SetXTitle("#Delta R K0 - jet");
- fhdRJetK0->SetYTitle("1/N_{jet} dN/dR");
- fhdRJetK0->Sumw2();
- fListOfHistos->Add(fhdRJetK0);
-
- fhdNdptV0 = new TH1F("hdNdptV0","dN/dpt V0",100,0,10);
- fhdNdptV0->SetXTitle("p_{T} (GeV/c)");
- fhdNdptV0->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptV0->Sumw2();
- fListOfHistos->Add(fhdNdptV0);
-
- fhdNdptK0 = new TH1F("hdNdptK0","",100,0,10);
- fhdNdptK0->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0->Sumw2();
- fListOfHistos->Add(fhdNdptK0);
-
- fhPtVsEtaK0 = new TH2F("hPtVsEtaK0","",20,-1,1,100,0,10);
- fhPtVsEtaK0->SetXTitle("#eta");
- fhPtVsEtaK0->SetYTitle("p_{T} (GeV/c)");
- fhPtVsEtaK0->Sumw2();
- fListOfHistos->Add(fhPtVsEtaK0);
-
- fhV0InvMassK0DCA = new TH1F("hV0InvMassK0DCA","",2000,0,2);
- fhV0InvMassK0DCA->SetXTitle("inv mass (GeV)");
- fhV0InvMassK0DCA->Sumw2();
- fListOfHistos->Add(fhV0InvMassK0DCA);
-
- fhV0InvMassK0DCAdEdx = new TH1F("hV0InvMassK0DCAdEdx","",2000,0,2);
- fhV0InvMassK0DCAdEdx->SetXTitle("inv mass (GeV)");
- fhV0InvMassK0DCAdEdx->Sumw2();
- fListOfHistos->Add(fhV0InvMassK0DCAdEdx);
-
- fhV0PtVsInvMassK0DCAdEdx = new TH2F("hV0PtVsInvMassK0DCAdEdx","",100,0.45,0.55,100,0.,50.);
- fhV0PtVsInvMassK0DCAdEdx->SetXTitle("inv mass (GeV)");
- fhV0PtVsInvMassK0DCAdEdx->SetYTitle("p_{T} (GeV)");
- fhV0PtVsInvMassK0DCAdEdx->Sumw2();
- fListOfHistos->Add(fhV0PtVsInvMassK0DCAdEdx);
-
- fhV0InvMassK0DCAPID = new TH1F("hV0InvMassK0DCAPID","",2000,0,2);
- fhV0InvMassK0DCAPID->SetXTitle("inv mass (GeV)");
- fhV0InvMassK0DCAPID->Sumw2();
- fListOfHistos->Add(fhV0InvMassK0DCAPID);
-
- fhV0InvMassLambdaDCAdEdx = new TH1F("hV0InvMassLambdaDCAdEdx","",2000,0,2);
- fhV0InvMassLambdaDCAdEdx->SetXTitle("inv mass (GeV)");
- fhV0InvMassLambdaDCAdEdx->Sumw2();
- fListOfHistos->Add(fhV0InvMassLambdaDCAdEdx);
-
- fhV0InvMassAntiLambdaDCAdEdx = new TH1F("hV0InvMassAntiLambdaDCAdEdx","",2000,0,2);
- fhV0InvMassAntiLambdaDCAdEdx->SetXTitle("inv mass (GeV)");
- fhV0InvMassAntiLambdaDCAdEdx->Sumw2();
- fListOfHistos->Add(fhV0InvMassAntiLambdaDCAdEdx);
-
- fhdNdptK0DCA = new TH1F("hdNdptK0DCA","",100,0,10);
- fhdNdptK0DCA->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0DCA->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0DCA->Sumw2();
- fListOfHistos->Add(fhdNdptK0DCA);
-
- fhdNdptK0DCAdEdx = new TH1F("hdNdptK0DCAdEdx","",100,0,10);
- fhdNdptK0DCAdEdx->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0DCAdEdx->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0DCAdEdx->Sumw2();
- fListOfHistos->Add(fhdNdptK0DCAdEdx);
-
- fhV0InvMassK0Min = new TH1F("hV0InvMassK0Min","",2000,0,2);
- fhV0InvMassK0Min->SetXTitle("inv mass (GeV)");
- fhV0InvMassK0Min->Sumw2();
- fListOfHistos->Add(fhV0InvMassK0Min);
-
- fhV0InvMassLambdaMin = new TH1F("hV0InvMassLambdaMin","",2000,0,2);
- fhV0InvMassLambdaMin->SetXTitle("inv mass (GeV)");
- fhV0InvMassLambdaMin->Sumw2();
- fListOfHistos->Add(fhV0InvMassLambdaMin);
-
- fhV0InvMassAntiLambdaMin = new TH1F("hV0InvMassAntiLambdaMin","",2000,0,2);
- fhV0InvMassAntiLambdaMin->SetXTitle("inv mass (GeV)");
- fhV0InvMassAntiLambdaMin->Sumw2();
- fListOfHistos->Add(fhV0InvMassAntiLambdaMin);
-
- fhV0InvMassK0Max = new TH1F("hV0InvMassK0Max","",2000,0,2);
- fhV0InvMassK0Max->SetXTitle("inv mass (GeV)");
- fhV0InvMassK0Max->Sumw2();
- fListOfHistos->Add(fhV0InvMassK0Max);
-
- fhV0InvMassLambdaMax = new TH1F("hV0InvMassLambdaMax","",2000,0,2);
- fhV0InvMassLambdaMax->SetXTitle("inv mass (GeV)");
- fhV0InvMassLambdaMax->Sumw2();
- fListOfHistos->Add(fhV0InvMassLambdaMax);
-
- fhV0InvMassAntiLambdaMax = new TH1F("hV0InvMassAntiLambdaMax","",2000,0,2);
- fhV0InvMassAntiLambdaMax->SetXTitle("inv mass (GeV)");
- fhV0InvMassAntiLambdaMax->Sumw2();
- fListOfHistos->Add(fhV0InvMassAntiLambdaMax);
-
- fhV0InvMassK0Jet = new TH1F("hV0InvMassK0Jet","",2000,0,2);
- fhV0InvMassK0Jet->SetXTitle("inv mass (GeV)");
- fhV0InvMassK0Jet->Sumw2();
- fListOfHistos->Add(fhV0InvMassK0Jet);
-
- fhV0InvMassLambdaJet = new TH1F("hV0InvMassLambdaJet","",2000,0,2);
- fhV0InvMassLambdaJet->SetXTitle("inv mass (GeV)");
- fhV0InvMassLambdaJet->Sumw2();
- fListOfHistos->Add(fhV0InvMassLambdaJet);
-
- fhV0InvMassAntiLambdaJet = new TH1F("hV0InvMassAntiLambdaJet","",2000,0,2);
- fhV0InvMassAntiLambdaJet->SetXTitle("inv mass (GeV)");
- fhV0InvMassAntiLambdaJet->Sumw2();
- fListOfHistos->Add(fhV0InvMassAntiLambdaJet);
-
- fhV0InvMassK0Lambda = new TH1F("hV0InvMassK0Lambda","",2000,0,2);
- fhV0InvMassK0Lambda->SetXTitle("inv mass (GeV)");
- fhV0InvMassK0Lambda->Sumw2();
- fListOfHistos->Add(fhV0InvMassK0Lambda);
-
- fhdNdptK0JetEvt = new TH1F("hdNdptK0JetEvt","",100,0,10);
- fhdNdptK0JetEvt->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0JetEvt->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0JetEvt->Sumw2();
- fListOfHistos->Add(fhdNdptK0JetEvt);
-
- fhdNdzK0 = new TH1F("hdNdzK0","",150,0,1.5);
- fhdNdzK0->SetXTitle("z");
- fhdNdzK0->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK0->Sumw2();
- fListOfHistos->Add(fhdNdzK0);
-
- fhdNdzK05to10 = new TH1F("hdNdzK05to10","",150,0,1.5);
- fhdNdzK05to10->SetXTitle("z");
- fhdNdzK05to10->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK05to10->Sumw2();
- fListOfHistos->Add(fhdNdzK05to10);
-
- fhdNdzK010to20 = new TH1F("hdNdzK010to20","",150,0,1.5);
- fhdNdzK010to20->SetXTitle("z");
- fhdNdzK010to20->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK010to20->Sumw2();
- fListOfHistos->Add(fhdNdzK010to20);
-
- fhdNdzK020to30 = new TH1F("hdNdzK020to30","",150,0,1.5);
- fhdNdzK020to30->SetXTitle("z");
- fhdNdzK020to30->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK020to30->Sumw2();
- fListOfHistos->Add(fhdNdzK020to30);
-
- fhdNdzK030to40 = new TH1F("hdNdzK030to40","",150,0,1.5);
- fhdNdzK030to40->SetXTitle("z");
- fhdNdzK030to40->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK030to40->Sumw2();
- fListOfHistos->Add(fhdNdzK030to40);
-
- fhdNdzK040to60 = new TH1F("hdNdzK040to60","",150,0,1.5);
- fhdNdzK040to60->SetXTitle("z");
- fhdNdzK040to60->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK040to60->Sumw2();
- fListOfHistos->Add(fhdNdzK040to60);
-
- fhdNdxiK0 = new TH1F("hdNdxiK0","",100,0,10);
- fhdNdxiK0->SetXTitle("xi");
- fhdNdxiK0->SetYTitle("1/N_{jet} dN/dxi");
- fhdNdxiK0->Sumw2();
- fListOfHistos->Add(fhdNdxiK0);
-
- fhdNdzLambda = new TH1F("hdNdzLambda","",150,0,1.5);
- fhdNdzLambda->SetXTitle("z");
- fhdNdzLambda->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzLambda->Sumw2();
- fListOfHistos->Add(fhdNdzLambda);
-
- fhdNdzAntiLambda = new TH1F("hdNdzAntiLambda","",150,0,1.5);
- fhdNdzAntiLambda->SetXTitle("z");
- fhdNdzAntiLambda->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzAntiLambda->Sumw2();
- fListOfHistos->Add(fhdNdzAntiLambda);
-
- fhdNdzK0Max = new TH1F("hdNdzK0Max","",150,0,1.5);
- fhdNdzK0Max->SetXTitle("z");
- fhdNdzK0Max->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK0Max->Sumw2();
- fListOfHistos->Add(fhdNdzK0Max);
-
- fhdNdxiK0Max = new TH1F("hdNdxiK0Max","",100,0,10);
- fhdNdxiK0Max->SetXTitle("xi");
- fhdNdxiK0Max->SetYTitle("1/N_{jet} dN/dxi");
- fhdNdxiK0Max->Sumw2();
- fListOfHistos->Add(fhdNdxiK0Max);
-
- fhdNdzLambdaMax = new TH1F("hdNdzLambdaMax","",150,0,1.5);
- fhdNdzLambdaMax->SetXTitle("z");
- fhdNdzLambdaMax->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzLambdaMax->Sumw2();
- fListOfHistos->Add(fhdNdzLambdaMax);
-
- fhdNdxiLambdaMax = new TH1F("hdNdxiLambdaMax","",700,0,7);
- fhdNdxiLambdaMax->SetXTitle("xi");
- fhdNdxiLambdaMax->SetYTitle("1/N_{jet} dN/dxi");
- fhdNdxiLambdaMax->Sumw2();
- fListOfHistos->Add(fhdNdxiLambdaMax);
-
- fhdNdptK0Max = new TH1F("hdNdptK0Max","",100,0,10);
- fhdNdptK0Max->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0Max->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0Max->Sumw2();
- fListOfHistos->Add(fhdNdptK0Max);
-
- fhdNdptLambdaMax = new TH1F("hdNdptLambdaMax","",100,0,10);
- fhdNdptLambdaMax->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaMax->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptLambdaMax->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaMax);
-
- fhdNdzK0Min = new TH1F("hdNdzK0Min","",150,0,1.5);
- fhdNdzK0Min->SetXTitle("z");
- fhdNdzK0Min->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK0Min->Sumw2();
- fListOfHistos->Add(fhdNdzK0Min);
-
- fhdNdxiK0Min = new TH1F("hdNdxiK0Min","",100,0,10);
- fhdNdxiK0Min->SetXTitle("xi");
- fhdNdxiK0Min->SetYTitle("1/N_{jet} dN/dxi");
- fhdNdxiK0Min->Sumw2();
- fListOfHistos->Add(fhdNdxiK0Min);
-
- fhdNdzLambdaMin = new TH1F("hdNdzLambdaMin","",150,0,1.5);
- fhdNdzLambdaMin->SetXTitle("z");
- fhdNdzLambdaMin->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzLambdaMin->Sumw2();
- fListOfHistos->Add(fhdNdzLambdaMin);
-
- fhdNdxiLambdaMin = new TH1F("hdNdxiLambdaMin","",700,0,7);
- fhdNdxiLambdaMin->SetXTitle("xi");
- fhdNdxiLambdaMin->SetYTitle("1/N_{jet} dN/dxi");
- fhdNdxiLambdaMin->Sumw2();
- fListOfHistos->Add(fhdNdxiLambdaMin);
-
- fhdNdptK0Min = new TH1F("hdNdptK0Min","",100,0,10);
- fhdNdptK0Min->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0Min->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0Min->Sumw2();
- fListOfHistos->Add(fhdNdptK0Min);
-
- fhdNdptLambdaMin = new TH1F("hdNdptLambdaMin","",100,0,10);
- fhdNdptLambdaMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaMin->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptLambdaMin->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaMin);
-
- fhdNdzK0Jet = new TH1F("hdNdzK0Jet","",150,0,1.5);
- fhdNdzK0Jet->SetXTitle("z");
- fhdNdzK0Jet->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK0Jet->Sumw2();
- fListOfHistos->Add(fhdNdzK0Jet);
-
- fhdNdxiK0Jet = new TH1F("hdNdxiK0Jet","",100,0,10);
- fhdNdxiK0Jet->SetXTitle("xi");
- fhdNdxiK0Jet->SetYTitle("1/N_{jet} dN/dxi");
- fhdNdxiK0Jet->Sumw2();
- fListOfHistos->Add(fhdNdxiK0Jet);
-
- fhdNdzLambdaJet = new TH1F("hdNdzLambdaJet","",150,0,1.5);
- fhdNdzLambdaJet->SetXTitle("z");
- fhdNdzLambdaJet->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzLambdaJet->Sumw2();
- fListOfHistos->Add(fhdNdzLambdaJet);
-
- fhdNdxiLambdaJet = new TH1F("hdNdxiLambdaJet","",700,0,7);
- fhdNdxiLambdaJet->SetXTitle("xi");
- fhdNdxiLambdaJet->SetYTitle("1/N_{jet} dN/dxi");
- fhdNdxiLambdaJet->Sumw2();
- fListOfHistos->Add(fhdNdxiLambdaJet);
-
- fhdNdptK0Jet = new TH1F("hdNdptK0Jet","",100,0,10);
- fhdNdptK0Jet->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0Jet->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0Jet->Sumw2();
- fListOfHistos->Add(fhdNdptK0Jet);
-
- fhdNdptLambdaJet = new TH1F("hdNdptLambdaJet","",100,0,10);
- fhdNdptLambdaJet->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaJet->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptLambdaJet->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaJet);
-
- fhdEdxVsMomV0 = new TH2F("hdEdxVsMomV0","",200,-10,10,200,0,2000);
- fhdEdxVsMomV0->SetXTitle("mom (GeV/c)");
- fhdEdxVsMomV0->SetYTitle("dE/dx");
- fhdEdxVsMomV0->Sumw2();
- fListOfHistos->Add(fhdEdxVsMomV0);
-
- fhdEdxVsMomV0pidEdx = new TH2F("hdEdxVsMomV0pidEdx","",200,-10,10,200,0,2000);
- fhdEdxVsMomV0pidEdx->SetXTitle("mom (GeV/c)");
- fhdEdxVsMomV0pidEdx->SetYTitle("dE/dx");
- fhdEdxVsMomV0pidEdx->Sumw2();
- fListOfHistos->Add(fhdEdxVsMomV0pidEdx);
-
- fhdEdxVsMomV0piPID = new TH2F("hdEdxVsMomV0piPID","",200,-10,10,200,0,2000);
- fhdEdxVsMomV0piPID->SetXTitle("mom (GeV/c)");
- fhdEdxVsMomV0piPID->SetYTitle("dE/dx");
- fhdEdxVsMomV0piPID->Sumw2();
- fListOfHistos->Add(fhdEdxVsMomV0piPID);
-
- fhdPhiJetK0MC = new TH1F("hdPhiJetK0MC","",640,-3.2,3.2);
- fhdPhiJetK0MC->SetXTitle("#Delta #phi K0 - jet");
- fhdPhiJetK0MC->SetYTitle("1/N_{jet} dN/dphi");
- fhdPhiJetK0MC->Sumw2();
- fListOfHistos->Add(fhdPhiJetK0MC);
-
- fhdRJetK0MC = new TH1F("hdRJetK0MC","dN/R K0-jet",500,0,5);
- fhdRJetK0MC->SetXTitle("#Delta R K0 - jet");
- fhdRJetK0MC->SetYTitle("1/N_{jet} dN/dR");
- fhdRJetK0MC->Sumw2();
- fListOfHistos->Add(fhdRJetK0MC);
-
- fhdRV0MC = new TH1F("hdRV0MC","",500,0.,1.);
- fhdRV0MC->SetXTitle("#Delta R");
- fhdRV0MC->SetYTitle("");
- fhdRV0MC->Sumw2();
- fListOfHistos->Add(fhdRV0MC);
-
- fhdNdptchPiMCMax = new TH1F("hdNdptchPiMCMax","",100,0,10);
- fhdNdptchPiMCMax->SetXTitle("p_{T} (GeV/c)");
- fhdNdptchPiMCMax->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptchPiMCMax->Sumw2();
- fListOfHistos->Add(fhdNdptchPiMCMax);
-
- fhdNdptK0MCMax = new TH1F("hdNdptK0MCMax","",100,0,10);
- fhdNdptK0MCMax->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0MCMax->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0MCMax->Sumw2();
- fListOfHistos->Add(fhdNdptK0MCMax);
-
- fhdNdptchKMCMax = new TH1F("hdNdptchKMCMax","",100,0,10);
- fhdNdptchKMCMax->SetXTitle("p_{T} (GeV/c)");
- fhdNdptchKMCMax->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptchKMCMax->Sumw2();
- fListOfHistos->Add(fhdNdptchKMCMax);
-
- fhdNdptpMCMax = new TH1F("hdNdptpMCMax","",100,0,10);
- fhdNdptpMCMax->SetXTitle("p_{T} (GeV/c)");
- fhdNdptpMCMax->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptpMCMax->Sumw2();
- fListOfHistos->Add(fhdNdptpMCMax);
-
- fhdNdptpBarMCMax = new TH1F("hdNdptpBarMCMax","",100,0,10);
- fhdNdptpBarMCMax->SetXTitle("p_{T} (GeV/c)");
- fhdNdptpBarMCMax->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptpBarMCMax->Sumw2();
- fListOfHistos->Add(fhdNdptpBarMCMax);
-
- fhdNdptLambdaMCMax = new TH1F("hdNdptLambdaMCMax","",100,0,10);
- fhdNdptLambdaMCMax->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaMCMax->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptLambdaMCMax->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaMCMax);
-
- fhdNdptLambdaBarMCMax = new TH1F("hdNdptLambdaBarMCMax","",100,0,10);
- fhdNdptLambdaBarMCMax->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaBarMCMax->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptLambdaBarMCMax->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaBarMCMax);
-
-
- fhdNdptchPiMCMin = new TH1F("hdNdptchPiMCMin","",100,0,10);
- fhdNdptchPiMCMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptchPiMCMin->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptchPiMCMin->Sumw2();
- fListOfHistos->Add(fhdNdptchPiMCMin);
-
- fhdNdptK0MCMin = new TH1F("hdNdptK0MCMin","",100,0,10);
- fhdNdptK0MCMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0MCMin->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0MCMin->Sumw2();
- fListOfHistos->Add(fhdNdptK0MCMin);
-
- fhdNdptchKMCMin = new TH1F("hdNdptchKMCMin","",100,0,10);
- fhdNdptchKMCMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptchKMCMin->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptchKMCMin->Sumw2();
- fListOfHistos->Add(fhdNdptchKMCMin);
-
- fhdNdptpMCMin = new TH1F("hdNdptpMCMin","",100,0,10);
- fhdNdptpMCMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptpMCMin->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptpMCMin->Sumw2();
- fListOfHistos->Add(fhdNdptpMCMin);
-
- fhdNdptpBarMCMin = new TH1F("hdNdptpBarMCMin","",100,0,10);
- fhdNdptpBarMCMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptpBarMCMin->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptpBarMCMin->Sumw2();
- fListOfHistos->Add(fhdNdptpBarMCMin);
-
- fhdNdptLambdaMCMin = new TH1F("hdNdptLambdaMCMin","",100,0,10);
- fhdNdptLambdaMCMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaMCMin->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptLambdaMCMin->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaMCMin);
-
- fhdNdptLambdaBarMCMin = new TH1F("hdNdptLambdaBarMCMin","",100,0,10);
- fhdNdptLambdaBarMCMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaBarMCMin->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptLambdaBarMCMin->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaBarMCMin);
-
- fhdNdptOmegaMCMin = new TH1F("hdNdptOmegaMCMin","",100,0,10);;
- fhdNdptOmegaMCMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptOmegaMCMin->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptOmegaMCMin->Sumw2();
- fListOfHistos->Add(fhdNdptOmegaMCMin);
-
- fhdNdptOmegaBarMCMin = new TH1F("hdNdptOmegaBarMCMin","",100,0,10);;
- fhdNdptOmegaBarMCMin->SetXTitle("p_{T} (GeV/c)");
- fhdNdptOmegaBarMCMin->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptOmegaBarMCMin->Sumw2();
- fListOfHistos->Add(fhdNdptOmegaBarMCMin);
-
- fhdNdptchPiMCJet = new TH1F("hdNdptchPiMCJet","",100,0,10);
- fhdNdptchPiMCJet->SetXTitle("p_{T} (GeV/c)");
- fhdNdptchPiMCJet->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptchPiMCJet->Sumw2();
- fListOfHistos->Add(fhdNdptchPiMCJet);
-
- fhdNdptK0MCJet = new TH1F("hdNdptK0MCJet","",100,0,10);
- fhdNdptK0MCJet->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0MCJet->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptK0MCJet->Sumw2();
- fListOfHistos->Add(fhdNdptK0MCJet);
-
- fhdNdptchKMCJet = new TH1F("hdNdptchKMCJet","",100,0,10);
- fhdNdptchKMCJet->SetXTitle("p_{T} (GeV/c)");
- fhdNdptchKMCJet->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptchKMCJet->Sumw2();
- fListOfHistos->Add(fhdNdptchKMCJet);
-
- fhdNdptpMCJet = new TH1F("hdNdptpMCJet","",100,0,10);
- fhdNdptpMCJet->SetXTitle("p_{T} (GeV/c)");
- fhdNdptpMCJet->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptpMCJet->Sumw2();
- fListOfHistos->Add(fhdNdptpMCJet);
-
- fhdNdptpBarMCJet = new TH1F("hdNdptpBarMCJet","",100,0,10);
- fhdNdptpBarMCJet->SetXTitle("p_{T} (GeV/c)");
- fhdNdptpBarMCJet->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptpBarMCJet->Sumw2();
- fListOfHistos->Add(fhdNdptpBarMCJet);
-
- fhdNdptLambdaMCJet = new TH1F("hdNdptLambdaMCJet","",100,0,10);
- fhdNdptLambdaMCJet->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaMCJet->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptLambdaMCJet->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaMCJet);
-
- fhdNdptLambdaBarMCJet = new TH1F("hdNdptLambdaBarMCJet","",100,0,10);
- fhdNdptLambdaBarMCJet->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaBarMCJet->SetYTitle("1/N_{jet} dN/dpt (1/GeV)");
- fhdNdptLambdaBarMCJet->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaBarMCJet);
-
- fhPIDMC = CreatePIDhisto("hPIDMC");
- fhPIDMC->Sumw2();
- fListOfHistos->Add(fhPIDMC);
-
- fhPIDMC_quarkEv = CreatePIDhisto("hPIDMC_quarkEv");
- fhPIDMC_quarkEv->Sumw2();
- fListOfHistos->Add(fhPIDMC_quarkEv);
-
- fhPIDMC_gluonEv = CreatePIDhisto("hPIDMC_gluonEv");
- fhPIDMC_gluonEv->Sumw2();
- fListOfHistos->Add(fhPIDMC_gluonEv);
-
- fhPIDMCAll = CreatePIDhisto("hPIDMCAll");
- fhPIDMCAll->Sumw2();
- fListOfHistos->Add(fhPIDMCAll);
-
- fhPIDMCMin = CreatePIDhisto("hPIDMCMin");
- fhPIDMCMin->Sumw2();
- fListOfHistos->Add(fhPIDMCMin);
-
- fhPIDMCJet = CreatePIDhisto("hPIDMCJet");
- fhPIDMCJet->Sumw2();
- fListOfHistos->Add(fhPIDMCJet);
-
- fhPIDMCMotherK0 = CreatePIDhisto("hPIDMCMotherK0");
- fhPIDMCMotherK0->Sumw2();
- fListOfHistos->Add(fhPIDMCMotherK0);
-
- fhPIDMCGrandMotherK0 = CreatePIDhisto("hPIDMCGrandMotherK0");
- fhPIDMCGrandMotherK0->Sumw2();
- fListOfHistos->Add(fhPIDMCGrandMotherK0);
-
- fhPIDMCMotherChK = CreatePIDhisto("hPIDMCMotherChK");
- fhPIDMCMotherChK->Sumw2();
- fListOfHistos->Add(fhPIDMCMotherChK);
-
- fhPIDMCMotherK0Trans = CreatePIDhisto("hPIDMCMotherK0Trans");
- fhPIDMCMotherK0Trans->Sumw2();
- fListOfHistos->Add(fhPIDMCMotherK0Trans);
-
- fhPIDMCGrandMotherK0Trans = CreatePIDhisto("hPIDMCGrandMotherK0Trans");
- fhPIDMCGrandMotherK0Trans->Sumw2();
- fListOfHistos->Add(fhPIDMCGrandMotherK0Trans);
-
- fhPIDMCMotherChKTrans = CreatePIDhisto("hPIDMCMotherChKTrans");
- fhPIDMCMotherChKTrans->Sumw2();
- fListOfHistos->Add(fhPIDMCMotherChKTrans);
-
- fhdNdptgammaMC = new TH1F("hdNdptgammaMC","",100,0,10);
- fhdNdptgammaMC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptgammaMC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptgammaMC->Sumw2();
- fListOfHistos->Add(fhdNdptgammaMC);
-
- fhdNdptchPiMC = new TH1F("hdNdptchPiMC","",100,0,10);;
- fhdNdptchPiMC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptchPiMC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptchPiMC->Sumw2();
- fListOfHistos->Add(fhdNdptchPiMC);
-
- fhdNdptpi0MC = new TH1F("hdNdptpi0MC","",100,0,10);;
- fhdNdptpi0MC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptpi0MC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptpi0MC->Sumw2();
- fListOfHistos->Add(fhdNdptpi0MC);
-
- fhdNdptK0MC = new TH1F("hdNdptK0MC","",100,0,10);;
- fhdNdptK0MC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0MC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptK0MC->Sumw2();
- fListOfHistos->Add(fhdNdptK0MC);
-
- fhdNdptchKMC = new TH1F("hdNdptchKMC","",100,0,10);;
- fhdNdptchKMC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptchKMC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptchKMC->Sumw2();
- fListOfHistos->Add(fhdNdptchKMC);
-
- fhdNdptpMC = new TH1F("hdNdptpMC","",100,0,10);;
- fhdNdptpMC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptpMC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptpMC->Sumw2();
- fListOfHistos->Add(fhdNdptpMC);
-
- fhdNdptpBarMC = new TH1F("hdNdptpBarMC","",100,0,10);;
- fhdNdptpBarMC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptpBarMC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptpBarMC->Sumw2();
- fListOfHistos->Add(fhdNdptpBarMC);
-
- fhdNdptLambdaMC = new TH1F("hdNdptLambdaMC","",100,0,10);;
- fhdNdptLambdaMC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaMC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptLambdaMC->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaMC);
-
- fhdNdptLambdaBarMC = new TH1F("hdNdptLambdaBarMC","",100,0,10);;
- fhdNdptLambdaBarMC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptLambdaBarMC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptLambdaBarMC->Sumw2();
- fListOfHistos->Add(fhdNdptLambdaBarMC);
-
- fhdNdptOmegaMC = new TH1F("hdNdptOmegaMC","",100,0,10);;
- fhdNdptOmegaMC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptOmegaMC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptOmegaMC->Sumw2();
- fListOfHistos->Add(fhdNdptOmegaMC);
-
- fhdNdptOmegaBarMC = new TH1F("hdNdptOmegaBarMC","",100,0,10);;
- fhdNdptOmegaBarMC->SetXTitle("p_{T} (GeV/c)");
- fhdNdptOmegaBarMC->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptOmegaBarMC->Sumw2();
- fListOfHistos->Add(fhdNdptOmegaBarMC);
-
- fhdNdxiMC = new TH1F("hdNdxiMC","",100,0,10);
- fhdNdxiMC->SetXTitle("#xi");
- fhdNdxiMC->Sumw2();
- fListOfHistos->Add(fhdNdxiMC);
-
- fhdNdxiK0MC = new TH1F("hdNdxiK0MC","",100,0,10);
- fhdNdxiK0MC->SetXTitle("#xi");
- fhdNdxiK0MC->Sumw2();
- fListOfHistos->Add(fhdNdxiK0MC);
-
- fhdNdxiK0MCJet = new TH1F("hdNdxiK0MCJet","",100,0,10);
- fhdNdxiK0MCJet->SetXTitle("#xi");
- fhdNdxiK0MCJet->Sumw2();
- fListOfHistos->Add(fhdNdxiK0MCJet);
-
- fhdNdzK0MC = new TH1F("hdNdzK0MC","",150,0,1.5);
- fhdNdzK0MC->SetXTitle("z");
- fhdNdzK0MC->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK0MC->Sumw2();
- fListOfHistos->Add(fhdNdzK0MC);
-
- fhdNdzK0MCJet = new TH1F("hdNdzK0MCJet","",150,0,1.5);
- fhdNdzK0MCJet->SetXTitle("z");
- fhdNdzK0MCJet->SetYTitle("1/N_{jet} dN/dz");
- fhdNdzK0MCJet->Sumw2();
- fListOfHistos->Add(fhdNdzK0MCJet);
-
- fhdNdptK0MCJetEvt = new TH1F("hdNdptK0MCJetEvt","",100,0,10);;
- fhdNdptK0MCJetEvt->SetXTitle("p_{T} (GeV/c)");
- fhdNdptK0MCJetEvt->SetYTitle("1/N_{event} dN/dpt (1/GeV)");
- fhdNdptK0MCJetEvt->Sumw2();
- fListOfHistos->Add(fhdNdptK0MCJetEvt);
-
- fhnJetsAODvsMC = new TH2F("hnJetsAODvsMC","",20,0,20,20,0.,20.);
- fhnJetsAODvsMC->SetXTitle("n jets MC");
- fhnJetsAODvsMC->SetXTitle("n jets AOD");
- fhnJetsAODvsMC->Sumw2();
- fListOfHistos->Add(fhnJetsAODvsMC);
-
- fhLeadingPtAODvsMC = new TH2F("hLeadingPtAODvsMC","",20,0,20,20,0.,20.);
- fhLeadingPtAODvsMC->SetXTitle("p_{T} MC (GeV/c)");
- fhLeadingPtAODvsMC->SetYTitle("p_{T} AOD (GeV/c)");
- fhLeadingPtAODvsMC->Sumw2();
- fListOfHistos->Add(fhLeadingPtAODvsMC);
-
- fhLeadingEtaAODvsMC = new TH2F("hLeadingEtaAODvsMC","",20,-1,1,20,-1.,1.);
- fhLeadingEtaAODvsMC->SetXTitle("#eta MC");
- fhLeadingEtaAODvsMC->SetYTitle("#eta AOD");
- fhLeadingEtaAODvsMC->Sumw2();
- fListOfHistos->Add(fhLeadingEtaAODvsMC);
-
- fhLeadingPhiAODvsMC = new TH2F("hLeadingPhiAODvsMC","",63,0,6.3,63,0.,6.3);
- fhLeadingPhiAODvsMC->SetXTitle("#phi MC");
- fhLeadingPhiAODvsMC->SetYTitle("#phi AOD");
- fhLeadingPhiAODvsMC->Sumw2();
- fListOfHistos->Add(fhLeadingPhiAODvsMC);
-
-
- fhnTracksLeadingAODvsMC = new TH2F("hnTracksLeadingAODvsMC","",20,0.,20,20,-0.5,19.5);
- fhnTracksLeadingAODvsMC->SetXTitle("nTracks MC");
- fhnTracksLeadingAODvsMC->SetYTitle("nTracks AOD");
- fhnTracksLeadingAODvsMC->Sumw2();
- fListOfHistos->Add(fhnTracksLeadingAODvsMC);
-
- fhLeadingdRAODMC = new TH1F("hLeadingdRAODMC","",40,0.,4);
- fhLeadingdRAODMC->SetXTitle("#Delta R");
- fhLeadingdRAODMC->Sumw2();
- fListOfHistos->Add(fhLeadingdRAODMC);
-
-
- fhLeadingPtAODvsMCdRcut = new TH2F("hLeadingPtAODvsMCdRcut","",40,0,20,40,0.,20.);
- fhLeadingPtAODvsMCdRcut->SetXTitle("p_{T} MC (GeV/c)");
- fhLeadingPtAODvsMCdRcut->SetYTitle("p_{T} AOD (GeV/c)");
- fhLeadingPtAODvsMCdRcut->Sumw2();
- fListOfHistos->Add(fhLeadingPtAODvsMCdRcut);
-
-
- fhdnTracksVsdPtLeadingAODMC = new TH2F("hdnTracksVsdPtLeadingAODMC","",40,-10.,10,20,-10.5,9.5);
- fhdnTracksVsdPtLeadingAODMC->SetXTitle("#Delta Pt AOD-MC (GeV/c)");
- fhdnTracksVsdPtLeadingAODMC->SetYTitle("#Delta N AOD-MC");
- fhdnTracksVsdPtLeadingAODMC->Sumw2();
- fListOfHistos->Add(fhdnTracksVsdPtLeadingAODMC);
-
- fhnTracksJetVsPtAOD = new TH2F("hnTracksJetVsPtAOD","",50,0,50,20,-0.5,19.5);
- fhnTracksJetVsPtAOD->SetXTitle("p_{T} (GeV/c)");
- fhnTracksJetVsPtAOD->SetYTitle("nTracks");
- fhnTracksJetVsPtAOD->Sumw2();
- fListOfHistos->Add(fhnTracksJetVsPtAOD);
-
- fhnTracksJetVsPtAODquarkEv = new TH2F("hnTracksJetVsPtAODquarkEv","",50,0,50,20,-0.5,19.5);
- fhnTracksJetVsPtAODquarkEv->SetXTitle("p_{T} (GeV/c)");
- fhnTracksJetVsPtAODquarkEv->SetYTitle("nTracks");
- fhnTracksJetVsPtAODquarkEv->Sumw2();
- fListOfHistos->Add(fhnTracksJetVsPtAODquarkEv);
-
- fhRadiusJetVsPtAOD = new TH2F("hRadiusJetVsPtAOD","",50,0,50,10,0.,1);
- fhRadiusJetVsPtAOD->SetXTitle("p_{T} (GeV/c)");
- fhRadiusJetVsPtAOD->SetYTitle("radius");
- fhRadiusJetVsPtAOD->Sumw2();
- fListOfHistos->Add(fhRadiusJetVsPtAOD);
-
- fhnTracksJetVsPtMC = new TH2F("hnTracksJetVsPtMC","",50,0,50,20,-0.5,19.5);
- fhnTracksJetVsPtMC->SetXTitle("p_{T} (GeV/c)");
- fhnTracksJetVsPtMC->SetYTitle("nTracks");
- fhnTracksJetVsPtMC->Sumw2();
- fListOfHistos->Add(fhnTracksJetVsPtMC);
-
- fhnTracksJetVsPtMCquarkEv = new TH2F("hnTracksJetVsPtMCquarkEv","",50,0,50,20,-0.5,19.5);
- fhnTracksJetVsPtMCquarkEv->SetXTitle("p_{T} (GeV/c)");
- fhnTracksJetVsPtMCquarkEv->SetYTitle("nTracks");
- fhnTracksJetVsPtMCquarkEv->Sumw2();
- fListOfHistos->Add(fhnTracksJetVsPtMCquarkEv);
-
- fhRadiusJetVsPtMC = new TH2F("hRadiusJetVsPtMC","",50,0,50,10,0.,1);
- fhRadiusJetVsPtMC->SetXTitle("p_{T} (GeV/c)");
- fhRadiusJetVsPtMC->SetYTitle("radius");
- fhRadiusJetVsPtMC->Sumw2();
- fListOfHistos->Add(fhRadiusJetVsPtMC);
-
- fhnTracksJetVsPtMCK0 = new TH2F("hnTracksJetVsPtMCK0","",50,0,50,20,-0.5,19.5);
- fhnTracksJetVsPtMCK0->SetXTitle("p_{T} (GeV/c)");
- fhnTracksJetVsPtMCK0->SetYTitle("nTracks");
- fhnTracksJetVsPtMCK0->Sumw2();
- fListOfHistos->Add(fhnTracksJetVsPtMCK0);
-
-
- fhnTracksJetVsPtMCK0quarkEv = new TH2F("hnTracksJetVsPtMCK0quarkEv","",50,0,50,20,-0.5,19.5);
- fhnTracksJetVsPtMCK0quarkEv->SetXTitle("p_{T} (GeV/c)");
- fhnTracksJetVsPtMCK0quarkEv->SetYTitle("nTracks");
- fhnTracksJetVsPtMCK0quarkEv->Sumw2();
- fListOfHistos->Add(fhnTracksJetVsPtMCK0quarkEv);
-
- fhRadiusJetVsPtMCK0 = new TH2F("hRadiusJetVsPtMCK0","",50,0,50,10,0.,1);
- fhRadiusJetVsPtMCK0->SetXTitle("p_{T} (GeV/c)");
- fhRadiusJetVsPtMCK0->SetYTitle("radius");
- fhRadiusJetVsPtMCK0->Sumw2();
- fListOfHistos->Add(fhRadiusJetVsPtMCK0);
-
- fhnTracksJetVsPtAODK0 = new TH2F("hnTracksJetVsPtAODK0","",50,0,50,20,-0.5,19.5);
- fhnTracksJetVsPtAODK0->SetXTitle("p_{T} (GeV/c)");
- fhnTracksJetVsPtAODK0->SetYTitle("nTracks AODK0");
- fhnTracksJetVsPtAODK0->Sumw2();
- fListOfHistos->Add(fhnTracksJetVsPtAODK0);
-
- fhnTracksJetVsPtAODK0quarkEv = new TH2F("hnTracksJetVsPtAODK0quarkEv","",50,0,50,20,-0.5,19.5);
- fhnTracksJetVsPtAODK0quarkEv->SetXTitle("p_{T} (GeV/c)");
- fhnTracksJetVsPtAODK0quarkEv->SetYTitle("nTracks AODK0quarkEv");
- fhnTracksJetVsPtAODK0quarkEv->Sumw2();
- fListOfHistos->Add(fhnTracksJetVsPtAODK0quarkEv);
-
- fhRadiusJetVsPtAODK0 = new TH2F("hRadiusJetVsPtAODK0","",50,0,50,10,0.,1);
- fhRadiusJetVsPtAODK0->SetXTitle("p_{T} (GeV/c)");
- fhRadiusJetVsPtAODK0->SetYTitle("radius");
- fhRadiusJetVsPtAODK0->Sumw2();
- fListOfHistos->Add(fhRadiusJetVsPtAODK0);
-
- fhnTracksJetVsPtAODpKch = new TH2F("hnTracksJetVsPtAODpKch","",50,0,50,20,-0.5,19.5);
- fhnTracksJetVsPtAODpKch->SetXTitle("p_{T} (GeV/c)");
- fhnTracksJetVsPtAODpKch->SetYTitle("nTracks AODpKch");
- fhnTracksJetVsPtAODpKch->Sumw2();
- fListOfHistos->Add(fhnTracksJetVsPtAODpKch);
-
- fhRadiusJetVsPtAODpKch = new TH2F("hRadiusJetVsPtAODpKch","",50,0,50,20,-0.5,19.5);
- fhRadiusJetVsPtAODpKch->SetXTitle("p_{T} (GeV/c)");
- fhRadiusJetVsPtAODpKch->SetYTitle("Radius AODpKch");
- fhRadiusJetVsPtAODpKch->Sumw2();
- fListOfHistos->Add(fhRadiusJetVsPtAODpKch);
-
- fhPythiaProcess = CreatePythiaIDhisto("hPythiaProcess");
- fListOfHistos->Add(fhPythiaProcess);
-
- fhPythiaProcessK0 = CreatePythiaIDhisto("hPythiaProcessK0");
- fListOfHistos->Add(fhPythiaProcessK0);
-
- fhPythiaProcessKch = CreatePythiaIDhisto("hPythiaProcessKch");
- fListOfHistos->Add(fhPythiaProcessKch);
-
- fhPythiaProcessp = CreatePythiaIDhisto("hPythiaProcessp");
- fListOfHistos->Add(fhPythiaProcessp);
-
- fhPythiaProcesspbar = CreatePythiaIDhisto("hPythiaProcesspbar");
- fListOfHistos->Add(fhPythiaProcesspbar);
-
- fhdNdzJets5to10 = new TH1F("hdNdzJets5to10","",25,0,1.25);
- fhdNdzJets5to10->SetXTitle("z");
- fhdNdzJets5to10->SetYTitle("dN/dz");
- fhdNdzJets5to10->Sumw2();
- fListOfHistos->Add(fhdNdzJets5to10);
-
- fhdNdzJets10to20 = new TH1F("hdNdzJets10to20","",25,0,1.25);
- fhdNdzJets10to20->SetXTitle("z");
- fhdNdzJets10to20->SetYTitle("dN/dz");
- fhdNdzJets10to20->Sumw2();
- fListOfHistos->Add(fhdNdzJets10to20);
-
- fhdNdzJets20to30 = new TH1F("hdNdzJets20to30","",25,0,1.25);
- fhdNdzJets20to30->SetXTitle("z");
- fhdNdzJets20to30->SetYTitle("dN/dz");
- fhdNdzJets20to30->Sumw2();
- fListOfHistos->Add(fhdNdzJets20to30);
-
- fhdNdzJets30to40 = new TH1F("hdNdzJets30to40","",25,0,1.25);
- fhdNdzJets30to40->SetXTitle("z");
- fhdNdzJets30to40->SetYTitle("dN/dz");
- fhdNdzJets30to40->Sumw2();
- fListOfHistos->Add(fhdNdzJets30to40);
-
- fhdNdzJets40to60 = new TH1F("hdNdzJets40to60","",25,0,1.25);
- fhdNdzJets40to60->SetXTitle("z");
- fhdNdzJets40to60->SetYTitle("dN/dz");
- fhdNdzJets40to60->Sumw2();
- fListOfHistos->Add(fhdNdzJets40to60);
-
-
- fhdNdxiJets5to10 = new TH1F("hdNdxiJets5to10","",70,0,7);
- fhdNdxiJets5to10->SetXTitle("z");
- fhdNdxiJets5to10->SetYTitle("dN/dz");
- fhdNdxiJets5to10->Sumw2();
- fListOfHistos->Add(fhdNdxiJets5to10);
-
- fhdNdxiJets10to20 = new TH1F("hdNdxiJets10to20","",70,0,7);
- fhdNdxiJets10to20->SetXTitle("z");
- fhdNdxiJets10to20->SetYTitle("dN/dz");
- fhdNdxiJets10to20->Sumw2();
- fListOfHistos->Add(fhdNdxiJets10to20);
-
- fhdNdxiJets20to30 = new TH1F("hdNdxiJets20to30","",70,0,7);
- fhdNdxiJets20to30->SetXTitle("z");
- fhdNdxiJets20to30->SetYTitle("dN/dz");
- fhdNdxiJets20to30->Sumw2();
- fListOfHistos->Add(fhdNdxiJets20to30);
-
- fhdNdxiJets30to40 = new TH1F("hdNdxiJets30to40","",70,0,7);
- fhdNdxiJets30to40->SetXTitle("z");
- fhdNdxiJets30to40->SetYTitle("dN/dz");
- fhdNdxiJets30to40->Sumw2();
- fListOfHistos->Add(fhdNdxiJets30to40);
-
- fhdNdxiJets40to60 = new TH1F("hdNdxiJets40to60","",70,0,7);
- fhdNdxiJets40to60->SetXTitle("z");
- fhdNdxiJets40to60->SetYTitle("dN/dz");
- fhdNdxiJets40to60->Sumw2();
- fListOfHistos->Add(fhdNdxiJets40to60);
-
- fhdNdptTracksJetPt5to10 = new TH1F("hdNdptTracksJetPt5to10","",250,0,25);
- fhdNdptTracksJetPt5to10->SetXTitle("p_{T} (GeV)");
- fhdNdptTracksJetPt5to10->SetYTitle("dN/dp_{T} 1/GeV");
- fhdNdptTracksJetPt5to10->Sumw2();
- fListOfHistos->Add(fhdNdptTracksJetPt5to10);
-
- fhdNdptTracksJetPt10to20 = new TH1F("hdNdptTracksJetPt10to20","",25,0,25);
- fhdNdptTracksJetPt10to20->SetXTitle("p_{T} (GeV)");
- fhdNdptTracksJetPt10to20->SetYTitle("dN/dp_{T} 1/GeV");
- fhdNdptTracksJetPt10to20->Sumw2();
- fListOfHistos->Add(fhdNdptTracksJetPt10to20);
-
- fhdNdptTracksJetPt20to30 = new TH1F("hdNdptTracksJetPt20to30","",25,0,25);
- fhdNdptTracksJetPt20to30->SetXTitle("p_{T} (GeV)");
- fhdNdptTracksJetPt20to30->SetYTitle("dN/dp_{T} 1/GeV");
- fhdNdptTracksJetPt20to30->Sumw2();
- fListOfHistos->Add(fhdNdptTracksJetPt20to30);
-
- fhdNdptTracksJetPt30to40 = new TH1F("hdNdptTracksJetPt30to40","",25,0,25);
- fhdNdptTracksJetPt30to40->SetXTitle("p_{T} (GeV)");
- fhdNdptTracksJetPt30to40->SetYTitle("dN/dp_{T} 1/GeV");
- fhdNdptTracksJetPt30to40->Sumw2();
- fListOfHistos->Add(fhdNdptTracksJetPt30to40);
-
- fhdNdptTracksJetPt40to60 = new TH1F("hdNdptTracksJetPt40to60","",25,0,25);
- fhdNdptTracksJetPt40to60->SetXTitle("p_{T} (GeV)");
- fhdNdptTracksJetPt40to60->SetYTitle("dN/dp_{T} 1/GeV");
- fhdNdptTracksJetPt40to60->Sumw2();
- fListOfHistos->Add(fhdNdptTracksJetPt40to60);
-
-
- fh1Xsec = new TProfile("h1Xsec","xsec from pyxsec.root",1,0,1);
- fh1Xsec->SetXTitle("<#sigma>");
- fh1Xsec->Sumw2();
- fListOfHistos->Add( fh1Xsec );
-
- fh1Trials = new TH1F("h1Trials","trials from pyxsec.root",1,0,1);
- fh1Trials->SetXTitle("#sum{ntrials}");
- fh1Trials->Sumw2();
- fListOfHistos->Add( fh1Trials );
-
-// fSettingsTree = new TTree("JetChemAnalysisSettings","Analysis Settings");
-// fSettingsTree->Branch("fUseLOConeJets",&fUseLOConeJets,"UseLOConeJets/O");
-// fSettingsTree->Branch("fUseLOConeMCJets",&fUseLOConeMCJets,"UseLOConeMCJets/O");
-// fSettingsTree->Branch("fUsePythiaJets",&fUsePythiaJets,"UsePythiaJets/O");
-// fSettingsTree->Branch("fConeRadius", &fConeRadius,"Rad/D");
-// fSettingsTree->Branch("fTrackPtCutJF", &fTrackPtCutJF,"TrackPt/D");
-// fSettingsTree->Branch("fFilterBitJF",&fFilterBitJF,"FilterBitJF/i");
-// fSettingsTree->Branch("fRequireITSRefitJF",&fRequireITSRefitJF,"RequireITSRefitJF/O");
-// fSettingsTree->Branch("fRejectK0TracksJF",&fRejectK0TracksJF,"RejectK0TracksJF/O");
-// fSettingsTree->Branch("fJetPtCut",&fJetPtCut,"JetPtCut/D");
-// fSettingsTree->Branch("fJetEtaCut",&fJetEtaCut,"JetEtaCut/D");
-// fSettingsTree->Branch("fFilterBit",&fFilterBit,"FilterBit/i");
-// fSettingsTree->Branch("fTrackPtCut",&fTrackPtCut,"TrackPtCut/D");
-// fSettingsTree->Branch("fTrackEtaCut",&fTrackEtaCut,"TrackEtaCut/D");
-// fSettingsTree->Branch("fUseOnFlyV0s",&fUseOnFlyV0s,"UseOnFlyV0s/O");
-// fSettingsTree->Branch("fCutnSigdEdx",&fCutnSigdEdx,"CutnSigdEdx/D");
-// fSettingsTree->Branch("fUseAODMCTracksForUE",&fUseAODMCTracksForUE,"UseAODMCTracksForUE/O");
-// fSettingsTree->Branch("fAreaReg",&fAreaReg,"AreaReg/D");
-// fSettingsTree->Branch("fAvgTrials",&fAvgTrials,"AvgTrials/D");
- // fListOfHistos->Add(fSettingsTree);
-
-}
-
-//____________________________________________________________________
-void AliAnalysisTaskJetChem::Terminate(Option_t */*option*/){
-
- // Terminate analysis
- //
-
- WriteSettings();
-
- // Update pointers reading them from the output slot
- fListOfHistos = dynamic_cast<TList*> (GetOutputData(1));
- if( !fListOfHistos ) {
- AliError("Histogram List is not available");
- return;
- }
-
- fhLeadingPt = (TH1F*)fListOfHistos->FindObject("hLeadingPt");
- fhRegionSumPtMaxVsEt = (TH1F*)fListOfHistos->FindObject("hRegionSumPtMaxVsEt");
- fhRegionSumPtMinVsEt = (TH1F*)fListOfHistos->FindObject("hRegionSumPtMinVsEt");
- fhRegionMultMaxVsEt = (TH1F*)fListOfHistos->FindObject("hRegionMultMaxVsEt");
- fhRegionMultMinVsEt = (TH1F*)fListOfHistos->FindObject("hRegionMultMinVsEt");
-
-
- fhRegionSumPtMaxVsEt->Divide(fhLeadingPt);
- fhRegionSumPtMinVsEt->Divide(fhLeadingPt);
- fhRegionMultMaxVsEt->Divide(fhLeadingPt);
- fhRegionMultMinVsEt->Divide(fhLeadingPt);
-
-
- //Get Normalization
- fh1Xsec = (TProfile*) fListOfHistos->FindObject("fh1Xsec");
- fh1Trials = (TH1F*)fListOfHistos->FindObject("fh1Trials");
-
- //std::cout<<" fh1Xsec "<<fh1Xsec<<" fh1Trials "<<fh1Trials<<std::endl;
-
- if(fh1Xsec && fh1Trials){
-
- Double_t xsec = fh1Xsec->GetBinContent(1);
- Double_t ntrials = fh1Trials->GetBinContent(1);
- Double_t normFactor = xsec/ntrials;
- Printf("xSec %f nTrials %f Norm %f \n",xsec,ntrials,normFactor);
- fhLeadingPt->Scale(normFactor);
- }
-
-
- if(fDebug > 1) AliInfo("End analysis");
-}
-
-// ----------------------------------------------------------------------------
-
-void AliAnalysisTaskJetChem::WriteSettings(){
-
- // fill settings tree
- // not on GRID !
-
- //fSettingsTree->Fill();
-
-}
-
-// ---------------------------------------------------------------------------