- fhnTrackPt[ij] = new THnSparseF(Form("fhnTrackPt%s",cAdd.Data()),";track number;p_{T};cent;#tracks;RP;trigger",nBinsSparse3,nBins3,xmin3,xmax3);
- fhnTrackPt[ij]->SetBinEdges(1,xPt3);
- fHistList->Add(fhnTrackPt[ij]);
- delete [] xPt3;
-
- // Track QA bins track nr, pTrack, cent, eta, phi bins 5.4 M
- const Int_t nBinsSparse4 = 6;
- const Int_t nBins4[nBinsSparse4] = { 2, 50, 10, 20, 180,2};
- const Double_t xmin4[nBinsSparse4] = { -0.5, 0, 0, -1.0, 0.,-1.5};
- const Double_t xmax4[nBinsSparse4] = { 1.5,150, 100, 1.0,2.*TMath::Pi(),1.5};
-
- // change the binning ot the pT axis:
- Double_t *xPt4 = new Double_t[nBins4[1]+1];
- xPt4[0] = 0.;
- for(int i = 1; i<=nBins4[1];i++){
- if(xPt4[i-1]<2)xPt4[i] = xPt4[i-1] + 0.1;
- else if(xPt4[i-1]<10)xPt4[i] = xPt4[i-1] + 0.5;
- else if(xPt4[i-1]<20)xPt4[i] = xPt4[i-1] + 1.;
- else if(xPt4[i-1]<30)xPt4[i] = xPt4[i-1] + 2.5;
- else xPt4[i] = xPt4[i-1] + 5.;
- }
- fhnTrackPtQA[ij] = new THnSparseF(Form("fhnTrackPtQA%s",cAdd.Data()),";track number;p_{T};cent;#eta;#phi;sign",nBinsSparse4,nBins4,xmin4,xmax4);
- fhnTrackPtQA[ij]->SetBinEdges(1,xPt4);
- fHistList->Add(fhnTrackPtQA[ij]);
- delete [] xPt4;
-
- for(int i = 0;i <= kMaxJets;++i){
- fh1PtIn[ij][i] = new TH1F(Form("fh1Pt%sIn_j%d",cAdd.Data(),i),Form("%s p_T input ;p_{T}",cAdd.Data()),nBinPt,binLimitsPt);
- fHistList->Add(fh1PtIn[ij][i]);
-
-
- if(!fh1TmpRho)fh1TmpRho = new TH1F("fh1TmpRho","tmp histo for jet shape",40,0.,2);
- fh2LTrackPtJetPt[ij][i] = new TH2F(Form("fh2LTrackPtJetPt%s_j%d",cAdd.Data(),i),
- Form("pt of leadin track within a jet vs jet %s;p_{T,lead in jet};p_{T.jet};",
- cAdd.Data()),
- 200,0.,200.,nBinPt,binLimitsPt);
- fHistList->Add(fh2LTrackPtJetPt[ij][i]);
- }
+ fhnTrackPt[ij] = new THnSparseF(Form("fhnTrackPt%s",cAdd.Data()),";track number;p_{T};cent;#tracks;RP;trigger",nBinsSparse3,nBins3,xmin3,xmax3);
+ fhnTrackPt[ij]->SetBinEdges(1,xPt3);
+ fHistList->Add(fhnTrackPt[ij]);
+ delete [] xPt3;
+
+ // Track QA bins track nr, pTrack, cent, eta, phi bins 5.4 M
+ const Int_t nBinsSparse4 = 6;
+ const Int_t nBins4[nBinsSparse4] = { 2, 50, 10, 20, 180,2};
+ const Double_t xmin4[nBinsSparse4] = { -0.5, 0, 0, -1.0, 0.,-1.5};
+ const Double_t xmax4[nBinsSparse4] = { 1.5,150, 100, 1.0,2.*TMath::Pi(),1.5};
+
+ // change the binning ot the pT axis:
+ Double_t *xPt4 = new Double_t[nBins4[1]+1];
+ xPt4[0] = 0.;
+ for(int i = 1; i<=nBins4[1];i++){
+ if(xPt4[i-1]<2)xPt4[i] = xPt4[i-1] + 0.1;
+ else if(xPt4[i-1]<10)xPt4[i] = xPt4[i-1] + 0.5;
+ else if(xPt4[i-1]<20)xPt4[i] = xPt4[i-1] + 1.;
+ else if(xPt4[i-1]<30)xPt4[i] = xPt4[i-1] + 2.5;
+ else xPt4[i] = xPt4[i-1] + 5.;
+ }
+ fhnTrackPtQA[ij] = new THnSparseF(Form("fhnTrackPtQA%s",cAdd.Data()),";track number;p_{T};cent;#eta;#phi;sign",nBinsSparse4,nBins4,xmin4,xmax4);
+ fhnTrackPtQA[ij]->SetBinEdges(1,xPt4);
+ fHistList->Add(fhnTrackPtQA[ij]);
+ delete [] xPt4;
+
+ for(int i = 0;i <= kMaxJets;++i){
+ fh1PtIn[ij][i] = new TH1F(Form("fh1Pt%sIn_j%d",cAdd.Data(),i),Form("%s p_T input ;p_{T}",cAdd.Data()),nBinPt,binLimitsPt);
+ fHistList->Add(fh1PtIn[ij][i]);