-//________________________________________________________________________
-void AliAnalysisTaskFragmentationFunction::AliFragFuncDiJetHistos::FillDiJetFF(Int_t jetType, Float_t trackPt, Float_t jetPt, Double_t jetBin, Bool_t incrementJetPt)
-{
- // fill DiJet FF
-
- if(jetType == 0)
- {
- if(incrementJetPt) fh1JetPt->Fill(jetPt);
-
- fh2TrackPtJet->Fill(jetBin, trackPt);
-
- Double_t z = trackPt / jetPt;
- Double_t xi = 0;
- if(z>0) xi = TMath::Log(1/z);
-
- fh2Xi->Fill(jetBin, xi);
- fh2Z->Fill(jetBin, z);
- }
- if(jetType == 1)
- {
- if(incrementJetPt) fh1Jet1Pt->Fill(jetPt);
-
- fh2TrackPtJet1->Fill(jetBin, trackPt);
-
- Double_t z = trackPt / jetPt;
- Double_t xi = 0;
- if(z>0) xi = TMath::Log(1/z);
-
- fh2Xi1->Fill(jetBin, xi);
- fh2Z1->Fill(jetBin, z);
- }
- if(jetType == 2)
- {
- if(incrementJetPt) fh1Jet2Pt->Fill(jetPt);
-
- fh2TrackPtJet2->Fill(jetBin, trackPt);
-
- Double_t z = trackPt / jetPt;
- Double_t xi = 0;
- if(z>0) xi = TMath::Log(1/z);
-
- fh2Xi2->Fill(jetBin, xi);
- fh2Z2->Fill(jetBin, z);
- }
-
-
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskFragmentationFunction::AliFragFuncDiJetHistos::AddToOutput(TList* list)const
-{
- // add histos to list
-
- list->Add(fh1Jet1Pt);
- list->Add(fh1Jet2Pt);
- list->Add(fh1JetPt);
- list->Add(fh2TrackPtJet1);
- list->Add(fh2TrackPtJet2);
- list->Add(fh2TrackPtJet);
- list->Add(fh2Xi1);
- list->Add(fh2Xi2);
- list->Add(fh2Xi);
- list->Add(fh2Z1);
- list->Add(fh2Z2);
- list->Add(fh2Z);
-}
-
-//______________________________________________________________________________________________________
-AliAnalysisTaskFragmentationFunction::AliFragFuncQADiJetHistos::AliFragFuncQADiJetHistos(const char* name, Int_t kindSlices,
- Int_t nInvMass, Float_t invMassMin, Float_t invMassMax,
- Int_t nJetPt, Float_t jetPtMin, Float_t jetPtMax,
- Int_t nDeltaPhi, Float_t deltaPhiMin, Float_t deltaPhiMax,
- Int_t nDeltaEta, Float_t deltaEtaMin, Float_t deltaEtaMax,
- Int_t nDeltaPt, Float_t deltaPtMin, Float_t deltaPtMax,
- Int_t nInBal, Float_t inBalMin, Float_t inBalMax)
- : TObject()
- ,fKindSlices(kindSlices)
- ,fNBinsJetInvMass(nInvMass)
- ,fJetInvMassMin(invMassMin)
- ,fJetInvMassMax(invMassMax)
- ,fNBinsJetPt(nJetPt)
- ,fJetPtMin(jetPtMin)
- ,fJetPtMax(jetPtMax)
- ,fNBinsDeltaPhi(nDeltaPhi)
- ,fDeltaPhiMin(deltaPhiMin)
- ,fDeltaPhiMax(deltaPhiMax)
- ,fNBinsDeltaEta(nDeltaEta)
- ,fDeltaEtaMin(deltaEtaMin)
- ,fDeltaEtaMax(deltaEtaMax)
- ,fNBinsDeltaPt(nDeltaPt)
- ,fDeltaPtMin(deltaPtMin)
- ,fDeltaPtMax(deltaPtMax)
- ,fNBinsInBal(nInBal)
- ,fInBalMin(inBalMin)
- ,fInBalMax(inBalMax)
- ,fh2InvMass(0)
- ,fh2DeltaPhi(0)
- ,fh2DeltaEta(0)
- ,fh2DeltaPt(0)
- ,fh2InBal(0)
- ,fNameQADJ(name)
-{
- // default constructor
-
-}
-
-//______________________________________________________________________________________________________
-AliAnalysisTaskFragmentationFunction::AliFragFuncQADiJetHistos::AliFragFuncQADiJetHistos(const AliFragFuncQADiJetHistos& copy)
- : TObject()
- ,fKindSlices(copy.fKindSlices)
- ,fNBinsJetInvMass(copy.fNBinsJetInvMass)
- ,fJetInvMassMin(copy.fJetInvMassMin)
- ,fJetInvMassMax(copy.fJetInvMassMax)
- ,fNBinsJetPt(copy.fNBinsJetPt)
- ,fJetPtMin(copy.fJetPtMin)
- ,fJetPtMax(copy.fJetPtMax)
- ,fNBinsDeltaPhi(copy.fNBinsDeltaPhi)
- ,fDeltaPhiMin(copy.fDeltaPhiMin)
- ,fDeltaPhiMax(copy.fDeltaPhiMax)
- ,fNBinsDeltaEta(copy.fNBinsDeltaEta)
- ,fDeltaEtaMin(copy.fDeltaEtaMin)
- ,fDeltaEtaMax(copy.fDeltaEtaMax)
- ,fNBinsDeltaPt(copy.fNBinsDeltaPt)
- ,fDeltaPtMin(copy.fDeltaPtMin)
- ,fDeltaPtMax(copy.fDeltaPtMax)
- ,fNBinsInBal(copy.fNBinsInBal)
- ,fInBalMin(copy.fInBalMin)
- ,fInBalMax(copy.fInBalMax)
- ,fh2InvMass(copy.fh2InvMass)
- ,fh2DeltaPhi(copy.fh2DeltaPhi)
- ,fh2DeltaEta(copy.fh2DeltaEta)
- ,fh2DeltaPt(copy.fh2DeltaPt)
- ,fh2InBal(copy.fh2InBal)
- ,fNameQADJ(copy.fNameQADJ)
-{
- // default constructor
-
-}
-
-//_______________________________________________________________________________________________________________________________________________________________
-AliAnalysisTaskFragmentationFunction::AliFragFuncQADiJetHistos& AliAnalysisTaskFragmentationFunction::AliFragFuncQADiJetHistos::operator=(const AliAnalysisTaskFragmentationFunction::AliFragFuncQADiJetHistos& o)
-{
- // assignment
-
- if(this!=&o){
- TObject::operator=(o);
- fKindSlices = o.fKindSlices;
- fNBinsJetInvMass = o.fNBinsJetInvMass;
- fJetInvMassMin = o.fJetInvMassMin;
- fJetInvMassMax = o.fJetInvMassMax;
- fNBinsJetPt = o.fNBinsJetPt;
- fJetPtMin = o.fJetPtMin;
- fJetPtMax = o.fJetPtMax;
- fNBinsDeltaPhi = o.fNBinsDeltaPhi;
- fDeltaPhiMin = o.fDeltaPhiMin;
- fDeltaPhiMax = o.fDeltaPhiMax;
- fNBinsDeltaEta = o.fNBinsDeltaEta;
- fDeltaEtaMin = o.fDeltaEtaMin;
- fDeltaEtaMax = o.fDeltaEtaMax;
- fNBinsDeltaPt = o.fNBinsDeltaPt;
- fDeltaPtMin = o.fDeltaPtMin;
- fDeltaPtMax = o.fDeltaPtMax;
- fNBinsInBal = o.fNBinsInBal;
- fInBalMin = o.fInBalMin;
- fInBalMax = o.fInBalMax;
- fh2InvMass = o.fh2InvMass;
- fh2DeltaPhi = o.fh2DeltaPhi;
- fh2DeltaEta = o.fh2DeltaEta;
- fh2DeltaPt = o.fh2DeltaPt;
- fh2InBal = o.fh2InBal;
- fNameQADJ = o.fNameQADJ;
- }
-
- return *this;
-}
-
-//_________________________________________________________
-AliAnalysisTaskFragmentationFunction::AliFragFuncQADiJetHistos::~AliFragFuncQADiJetHistos()
-{
- // destructor
-
- if(fh2InvMass) delete fh2InvMass;
- if(fh2DeltaPhi) delete fh2DeltaPhi;
- if(fh2DeltaEta) delete fh2DeltaEta;
- if(fh2DeltaPt) delete fh2DeltaPt;
- if(fh2InBal) delete fh2InBal;
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskFragmentationFunction::AliFragFuncQADiJetHistos::DefineQADiJetHistos()
-{
- // define histos
-
- Int_t nBins = 0;
- Double_t min = 0.;
- Double_t max = 0.;
- const char *xaxis = "";
- if(fKindSlices == 1)
- {
- nBins = fNBinsJetInvMass;
- min = fJetInvMassMin;
- max = fJetInvMassMax;
- xaxis = "M_{JJ} [GeV]";
- }
- if(fKindSlices == 2 || fKindSlices == 3)
- {
- nBins = fNBinsJetPt;
- min = fJetPtMin;
- max = fJetPtMax;
- if(fKindSlices == 2) xaxis = "E_{Tmean} [GeV]";
- if(fKindSlices == 3) xaxis ="leading jet p_{T} [GeV/c]";
- }
-
-
- fh2InvMass = new TH2F(Form("fh2DJInvMassPositionCut%s", fNameQADJ.Data()), "",nBins, min, max, fNBinsJetInvMass, fJetInvMassMin, fJetInvMassMax);
- fh2DeltaPhi = new TH2F(Form("fh2DJDeltaPhiPositionCut%s", fNameQADJ.Data()), "",nBins, min, max, fNBinsDeltaPhi, fDeltaPhiMin, fDeltaPhiMax);
- fh2DeltaEta = new TH2F(Form("fh2DJDeltaEtaPositionCut%s", fNameQADJ.Data()), "",nBins, min, max, fNBinsDeltaEta, fDeltaEtaMin, fDeltaEtaMax);
- fh2DeltaPt = new TH2F(Form("fh2DJDeltaPtPositionCut%s", fNameQADJ.Data()), "",nBins, min, max, fNBinsDeltaPt, fDeltaPtMin, fDeltaPtMax);
- fh2InBal = new TH2F(Form("fh2DJInBalPositionCut%s", fNameQADJ.Data()), "",nBins, min, max, fNBinsInBal, fInBalMin, fInBalMax);
-
- AliAnalysisTaskFragmentationFunction::SetProperties(fh2InvMass, xaxis, "Invariant Mass", "Entries");
- AliAnalysisTaskFragmentationFunction::SetProperties(fh2DeltaPhi, xaxis, "#Delta #phi", "Entries");
- AliAnalysisTaskFragmentationFunction::SetProperties(fh2DeltaEta, xaxis, "#Delta #eta", "Entries");
- AliAnalysisTaskFragmentationFunction::SetProperties(fh2DeltaPt, xaxis, "#Delta p_{T}", "Entries");
- AliAnalysisTaskFragmentationFunction::SetProperties(fh2InBal, xaxis, "(p_{T}^{1}-p_{T}^{2})/(p_{T}^{1}+p_{T}^{2})", "Entries");
-
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskFragmentationFunction::AliFragFuncQADiJetHistos::FillDiJetQA(Double_t invMass, Double_t deltaPhi, Double_t deltaEta,Double_t deltaPt, Double_t inbal, Double_t jetBin)
-{
- // fill dijet QA
-
- fh2InvMass->Fill(jetBin, invMass);
- fh2DeltaPhi->Fill(jetBin, deltaPhi);
- fh2DeltaEta->Fill(jetBin, deltaEta);
- fh2DeltaPt->Fill(jetBin, deltaPt);
- fh2InBal->Fill(jetBin, inbal);
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskFragmentationFunction::AliFragFuncQADiJetHistos::AddToOutput(TList* list)const
-{
- // add histos to list
-
- list->Add(fh2InvMass);
- list->Add(fh2DeltaPhi);
- list->Add(fh2DeltaEta);
- list->Add(fh2DeltaPt);
- list->Add(fh2InBal);
-}
-
-//_________________________________________________________________________________
-void AliAnalysisTaskFragmentationFunction::AliFragFuncIntraJetHistos::AddToOutput(TList* list) const
-{
- // add histos to list
-
- list->Add(fh2CosTheta);
- list->Add(fh2PtZ);
- list->Add(fh3ThetaZ);
- list->Add(fh3JtTheta);
- list->Add(fh3JtZ);
-
-}
-
-//_________________________________________________________________________________
-Bool_t AliAnalysisTaskFragmentationFunction::Notify()
-{
- //
- // Implemented Notify() to read the cross sections
- // and number of trials from pyxsec.root
- // (taken from AliAnalysisTaskJetSpectrum2)
- //
- TTree *tree = AliAnalysisManager::GetAnalysisManager()->GetTree();
- Float_t xsection = 0;
- Float_t ftrials = 1;
-
- fAvgTrials = 1;
- if(tree){
- TFile *curfile = tree->GetCurrentFile();
- if (!curfile) {
- Error("Notify","No current file");
- return kFALSE;
- }
- if(!fh1Xsec||!fh1Trials){
- Printf("%s%d No Histogram fh1Xsec",(char*)__FILE__,__LINE__);
- return kFALSE;
- }
- AliAnalysisHelperJetTasks::PythiaInfoFromFile(curfile->GetName(),xsection,ftrials);
- fh1Xsec->Fill("<#sigma>",xsection);
- // construct a poor man average trials
- Float_t nEntries = (Float_t)tree->GetTree()->GetEntries();
- if(ftrials>=nEntries && nEntries>0.)fAvgTrials = ftrials/nEntries;
- }
-
- // Set seed for backg study
- fRandom = new TRandom3();
- fRandom->SetSeed(0);
-
- return kTRUE;
-}
-
-
-
-//__________________________________________________________________
-void AliAnalysisTaskFragmentationFunction::UserCreateOutputObjects()
-{
- // create output objects
-
- if(fDebug > 1) Printf("AliAnalysisTaskFragmentationFunction::UserCreateOutputObjects()");
-
- // create list of tracks and jets
-
- fTracksRec = new TList();
- fTracksRec->SetOwner(kFALSE);
-
- fTracksRecCuts = new TList();
- fTracksRecCuts->SetOwner(kFALSE);
-
- fTracksGen = new TList();
- fTracksGen->SetOwner(kFALSE);
-
- fTracksAODMCCharged = new TList();
- fTracksAODMCCharged->SetOwner(kFALSE);
-
- fTracksAODMCChargedSec = new TList();
- fTracksAODMCChargedSec->SetOwner(kFALSE);
-
- fTracksRecQualityCuts = new TList();
- fTracksRecQualityCuts->SetOwner(kFALSE);
-
- fJetsRec = new TList();
- fJetsRec->SetOwner(kFALSE);
- if(fBranchRecJets.Contains("KT") && fBckgSubMethod) fJetsRec->SetOwner(kTRUE);
-
- fJetsRecCuts = new TList();
- fJetsRecCuts->SetOwner(kFALSE);
- if(fBranchRecJets.Contains("KT") && fBckgSubMethod) fJetsRecCuts->SetOwner(kTRUE);
-
- fJetsGen = new TList();
- fJetsGen->SetOwner(kFALSE);
-
- fJetsRecEff = new TList();
- fJetsRecEff->SetOwner(kFALSE);
-
- fJetsEmbedded = new TList();
- fJetsEmbedded->SetOwner(kFALSE);
-
- // fJetsKine = new TList();
- // fJetsKine->SetOwner(kTRUE); // delete AOD jets using mom from Kine Tree via TList::Clear()
-
- if(fBckgMode &&
- (fBckgType[0]==kBckgClusters || fBckgType[1]==kBckgClusters || fBckgType[2]==kBckgClusters || fBckgType[3]==kBckgClusters || fBckgType[4]==kBckgClusters ||
- fBckgType[0]==kBckgClustersOutLeading || fBckgType[1]==kBckgClustersOutLeading || fBckgType[2]==kBckgClustersOutLeading ||
- fBckgType[3]==kBckgClustersOutLeading || fBckgType[4]==kBckgClustersOutLeading)){
-
- fBckgJetsRec = new TList();
- fBckgJetsRec->SetOwner(kFALSE);
-
- fBckgJetsRecCuts = new TList();
- fBckgJetsRecCuts->SetOwner(kFALSE);
-
- fBckgJetsGen = new TList();
- fBckgJetsGen->SetOwner(kFALSE);
- }
-
- //
- // Create histograms / output container
- //
-
- OpenFile(1);
- fCommonHistList = new TList();
-
- Bool_t oldStatus = TH1::AddDirectoryStatus();
- TH1::AddDirectory(kFALSE);
-
-
- // Histograms
- fh1EvtSelection = new TH1F("fh1EvtSelection", "Event Selection", 6, -0.5, 5.5);
- fh1EvtSelection->GetXaxis()->SetBinLabel(1,"ACCEPTED");
- fh1EvtSelection->GetXaxis()->SetBinLabel(2,"event selection: rejected");
- fh1EvtSelection->GetXaxis()->SetBinLabel(3,"event class: rejected");
- fh1EvtSelection->GetXaxis()->SetBinLabel(4,"vertex Ncontr: rejected");
- fh1EvtSelection->GetXaxis()->SetBinLabel(5,"vertex z: rejected");
- fh1EvtSelection->GetXaxis()->SetBinLabel(6,"vertex type: rejected");
-
- fh1VertexNContributors = new TH1F("fh1VertexNContributors", "Vertex N contributors", 2500,-.5, 2499.5);
- fh1VertexZ = new TH1F("fh1VertexZ", "Vertex z distribution", 30, -15., 15.);
- fh1EvtMult = new TH1F("fh1EvtMult","Event multiplicity, track pT cut > 150 MeV/c, |#eta| < 0.9",120,0.,12000.);
- fh1EvtCent = new TH1F("fh1EvtCent","centrality",100,0.,100.);
- fh2TrackPtVsDCAXY = new TH2F("fh2TrackPtVsDCAXY","",400,-0.02,0.02,100,0.,10.);
- fh2TrackPtVsDCAZ = new TH2F("fh2TrackPtVsDCAZ","",100,-0.050,0.050,100,0.,10.);
-
- fh1Xsec = new TProfile("fh1Xsec","xsec from pyxsec.root",1,0,1);
- fh1Xsec->GetXaxis()->SetBinLabel(1,"<#sigma>");
- fh1Trials = new TH1F("fh1Trials","trials from pyxsec.root",1,0,1);
- fh1Trials->GetXaxis()->SetBinLabel(1,"#sum{ntrials}");
- fh1PtHard = new TH1F("fh1PtHard","PYTHIA Pt hard;p_{T,hard}",350,-.5,349.5);
- fh1PtHardTrials = new TH1F("fh1PtHardTrials","PYTHIA Pt hard weight with trials;p_{T,hard}",350,-.5,349.5);
-
- fh1nRecJetsCuts = new TH1F("fh1nRecJetsCuts","reconstructed jets per event",10,-0.5,9.5);
- fh1nGenJets = new TH1F("fh1nGenJets","generated jets per event",10,-0.5,9.5);
- fh1nRecEffJets = new TH1F("fh1nRecEffJets","reconstruction effiency: jets per event",10,-0.5,9.5);
- fh1nEmbeddedJets = new TH1F("fh1nEmbeddedJets","embedded jets per event",10,-0.5,9.5);
-
- fh2PtRecVsGenPrim = new TH2F("fh2PtRecVsGenPrim","rec vs gen pt",fQATrackNBinsPt,fQATrackPtMin,fQATrackPtMax,fQATrackNBinsPt,fQATrackPtMin,fQATrackPtMax);
- fh2PtRecVsGenSec = new TH2F("fh2PtRecVsGenSec","rec vs gen pt",fQATrackNBinsPt,fQATrackPtMin,fQATrackPtMax,fQATrackNBinsPt,fQATrackPtMin,fQATrackPtMax);
-
-
- // Background
- if(fBckgMode) {
- if(fBckgType[0]==kBckgClusters || fBckgType[1]==kBckgClusters || fBckgType[2]==kBckgClusters || fBckgType[3]==kBckgClusters || fBckgType[4]==kBckgClusters ||
- fBckgType[0]==kBckgClustersOutLeading || fBckgType[1]==kBckgClustersOutLeading || fBckgType[2]==kBckgClustersOutLeading ||
- fBckgType[3]==kBckgClustersOutLeading || fBckgType[4]==kBckgClustersOutLeading){
-
- fh1nRecBckgJetsCuts = new TH1F("fh1nRecBckgJetsCuts","reconstructed background jets per event",10,-0.5,9.5);
- fh1nGenBckgJets = new TH1F("fh1nGenBckgJets","generated background jets per event",10,-0.5,9.5);
- }
-
- if(fBckgType[0]==kBckgPerp || fBckgType[1]==kBckgPerp || fBckgType[2]==kBckgPerp || fBckgType[3]==kBckgPerp || fBckgType[4]==kBckgPerp )
- fh1PerpMult = new TH1F("fh1PerpMult","Background multiplicity - Cone perpendicular to leading jet axis",500,0.,500.);
- if(fBckgType[0]==kBckgASide || fBckgType[1]==kBckgASide || fBckgType[2]==kBckgASide || fBckgType[3]==kBckgASide || fBckgType[4]==kBckgASide)
- fh1ASideMult = new TH1F("fh1ASideMult","Background multiplicity - Cone in the away side of leading jet axis",500,0.,500.);
- if(fBckgType[0]==kBckgASideWindow || fBckgType[1]==kBckgASideWindow || fBckgType[2]==kBckgASideWindow || fBckgType[3]==kBckgASideWindow || fBckgType[4]==kBckgASideWindow)
- fh1ASideWindowMult = new TH1F("fh1ASideWindowMult","Background multiplicity - Cone in the away side of leading jet axis",500,0.,500.);
- if(fBckgType[0]==kBckgPerpWindow || fBckgType[1]==kBckgPerpWindow || fBckgType[2]==kBckgPerpWindow || fBckgType[3]==kBckgPerpWindow || fBckgType[4]==kBckgPerpWindow)
- fh1PerpWindowMult = new TH1F("fh1PerpWindowMult","Background multiplicity - Cone in the perp direction of leading jet axis",500,0.,500.);
- if(fBckgType[0]==kBckgOutLJ || fBckgType[1]==kBckgOutLJ || fBckgType[2]==kBckgOutLJ || fBckgType[3]==kBckgOutLJ || fBckgType[4]==kBckgOutLJ)
- fh1OutLeadingMult = new TH1F("fh1OutLeadingMult","Background multiplicity - Cone outside leading jet",500,0,500.);
- if(fBckgType[0]==kBckgOutLJStat || fBckgType[1]==kBckgOutLJStat || fBckgType[2]==kBckgOutLJStat || fBckgType[3]==kBckgOutLJStat || fBckgType[4]==kBckgOutLJStat)
- fh1OutLeadingStatMult = new TH1F("fh1OutLeadingStatMult","Background multiplicity - Cone outside leading jet",3000,0,3000.);
- if(fBckgType[0]==kBckgOut2J || fBckgType[1]==kBckgOut2J || fBckgType[2]==kBckgOut2J || fBckgType[3]==kBckgOut2J || fBckgType[4]==kBckgOut2J)
- fh1Out2JetsMult = new TH1F("fh1Out2JetsMult","Background multiplicity - Cone outside 2 jets",500,0.,500.);
- if(fBckgType[0]==kBckgOut3J || fBckgType[1]==kBckgOut3J || fBckgType[2]==kBckgOut3J || fBckgType[3]==kBckgOut3J || fBckgType[4]==kBckgOut3J)
- fh1Out3JetsMult = new TH1F("fh1Out3JetsMult","Background multiplicity - Cone outside 3 jets",500,0.,500.);
- if(fBckgType[0]==kBckgClusters || fBckgType[1]==kBckgClusters || fBckgType[2]==kBckgClusters || fBckgType[3]==kBckgClusters || fBckgType[4]==kBckgClusters)
- fh1MedianClustersMult = new TH1F("fh1MedianClustersMult","Background multiplicity - median cluster",500,0.,500.);
- if(fBckgType[0]==kBckgClustersOutLeading || fBckgType[1]==kBckgClustersOutLeading || fBckgType[2]==kBckgClustersOutLeading || fBckgType[3]==kBckgClustersOutLeading
- || fBckgType[4]==kBckgClustersOutLeading)
- fh1OutClustersMult = new TH1F("fh1OutClustersMult","Background multiplicity - clusters outside leading jet",3000,0.,3000.);
- }
-
- // embedding
- if(fBranchEmbeddedJets.Length()){
- fh1FractionPtEmbedded = new TH1F("fh1FractionPtEmbedded","",200,0,2);
- fh1IndexEmbedded = new TH1F("fh1IndexEmbedded","",11,-1,10);
- fh2DeltaPtVsJetPtEmbedded = new TH2F("fh2DeltaPtVsJetPtEmbedded","",250,0,250,200,-100,100);
- fh2DeltaPtVsRecJetPtEmbedded = new TH2F("fh2DeltaPtVsRecJetPtEmbedded","",250,0,250,200,-100,100);
- fh1DeltaREmbedded = new TH1F("fh1DeltaREmbedded","",50,0,0.5);
- fh1nEmbeddedJets = new TH1F("fh1nEmbeddedJets","embedded jets per event",10,-0.5,9.5);
- }
-
- if(fEffMode){
- fh2ptVsDistNN_pt50_rec = new TH2F("fh2ptVsDistNN_pt50_rec","",200,0,0.2,500,0.,100);
- fh2ptVsDistNN_pt50_nonRec = new TH2F("fh2ptVsDistNN_pt50_nonRec","",200,0,0.2,500,0.,100);
- fh2ptVsDistNN_pt10_rec = new TH2F("fh2ptVsDistNN_pt10_rec","",200,0,0.2,500,0.,100);
- fh2ptVsDistNN_pt10_nonRec = new TH2F("fh2ptVsDistNN_pt10_nonRec","",200,0,0.2,500,0.,100);
- }