- if(fbJetsRp){
- // cent : nInpTrks : rp bins : rp wrt jet(2x) : probe pT : probe area: deltaPt : rp phi : rho : correction for RP : local rho
- /*
- entries = 1<<0 | 1<<1 | 1<<3 | 1<<4 | 1<<5 | 1<<6 | 1<<12 | 1<<14 | 1<<2 | 1<<22 | 1<<23 | 1<<24 | 1<<25;
- opt = 1<<4 | 1<<5;*/
- // cent : nInpTrks : rp bins : rp wrt jet(2x) : probe pT : deltaPt : rp phi : rho : correction for RP | pT hard bin
- entries = 1<<0 | 1<<1 | 1<<3 | 1<<4 | 1<<5 | 1<<6 | 1<<14 | 1<<2 | 1<<22 | 1<<23 | 1<<26;
- opt = 1<<4 | 1<<5;
- fhnJetsRp = NewTHnSparseF("fhnJetsRp", entries, opt);
- }
-
- // cent : nInpTrks : rp bins: deltaPt : jetPt(2x) : deltaArea : pT hard bin (hr delta pt)
- if(fbJetsDeltaPt){
- entries = 1<<0 | 1<<1 | 1<<3 | 1<<14 | 1<<6 | 1<<7 | 1<<18 | 1<<26;
- opt = 1<<1 | 1<<14 | 1<<6 | 1<<7 | 1<<18;
- fhnJetsDeltaPt = NewTHnSparseF("fhnJetsDeltaPt", entries, opt);
- }
-
- // cent : nInpTrks : rp bins : deltaPt : jetPt(2x) : deltaR : deltaEta : jetEta(2x) : pT hard bin (hr for eta)
- if(fbJetsEta){
- entries = 1<<0 | 1<<1 | 1<<3 | 1<<14 | 1<<6 | 1<<7 | 1<<17 | 1<<15 | 1<<8 | 1<<9 | 1<<26;
- opt = 1<<15 | 1<<8 | 1<<9;
- fhnJetsEta = NewTHnSparseF("fhnJetsEta", entries, opt);
- }
-
- // cent : nInpTrks : rp bins : jetPt(2x) : jetPhi(2x) : deltaPt : deltaPhi : pT hard bin
- if(fbJetsPhi){
- entries = 1<<0 | 1<<1 | 1<<3 | 1<<6 | 1<<7 | 1<<10 | 1<<11 | 1<<14 | 1<<16 | 1<<26;
- opt = 1<<10 | 1<<11;
- fhnJetsPhi = NewTHnSparseF("fhnJetsPhi", entries, opt);
- }
-
- // cent : nInpTrks : rp bins : deltaArea : jetArea(2x) : deltaR : fraction : distance next rec jet : pT next jet : deltaPt : jetPt(2x) : pT hard bin (hr for area)
- if(fbJetsArea){
- entries = 1<<0 | 1<<1 | 1<<3 | 1<<18 | 1<<12 | 1<<13 | 1<<17 | 1<<19 | 1<<20 | 1<<21 | 1<<14 | 1<<6 | 1<<7 | 1<<26;
- opt = 1<<18 | 1<<12 | 1<<13;
- fhnJetsArea = NewTHnSparseF("fhnJetsArea", entries, opt);
- }
-
-
- //before cut
-
- // cent : nInpTrks : rp bins : fraction : jetPt(2x) : jetEta(2x) : jetPhi(2x) (low resolution) (with fraction, eta, phi, pt cuts possible)
- if(fbJetsBeforeCut1){
- entries = 1<<0 | 1<<1 | 1<<3 | 1<<19 | 1<<6 | 1<<7 | 1<<8 | 1<<9 | 1<<10 | 1<<11;
- opt = 0;
- fhnJetsBeforeCut1 = NewTHnSparseF("fhnJetsBeforeCut1", entries, opt);
- }
-
- // cent : nInpTrks : rp bins : deltaPt : jetPt(2x) : deltaR : deltaEta : jetEta(2x) (low resolution)
- if(fbJetsBeforeCut2){
- entries = 1<<0 | 1<<1 | 1<<3 | 1<<14 | 1<<6 | 1<<7 | 1<<17 | 1<<15 | 1<<8 | 1<<9;
- opt = 0;
- fhnJetsBeforeCut2 = NewTHnSparseF("fhnJetsBeforeCut2", entries, opt);
- }
-
- fOutputList->Add(fHistEvtSelection);
- fOutputList->Add(fHistJetSelection);
- fOutputList->Add(fh2JetSelection);
- if(fbEvent) fOutputList->Add(fhnEvent);
- if(fbJetsMismatch1) fOutputList->Add(fhnJetsMismatch1);
- if(fbJetsMismatch2) fOutputList->Add(fhnJetsMismatch2);
- if(fbJetsRp) fOutputList->Add(fhnJetsRp);
- if(fbJetsDeltaPt) fOutputList->Add(fhnJetsDeltaPt);
- if(fbJetsEta) fOutputList->Add(fhnJetsEta);
- if(fbJetsPhi) fOutputList->Add(fhnJetsPhi);
- if(fbJetsArea) fOutputList->Add(fhnJetsArea);
- if(fbJetsBeforeCut1) fOutputList->Add(fhnJetsBeforeCut1);
- if(fbJetsBeforeCut2) fOutputList->Add(fhnJetsBeforeCut2);
-
- // =========== Switch on Sumw2 for all histos ===========
- for (Int_t i=0; i<fOutputList->GetEntries(); ++i) {
- TH1 *h1 = dynamic_cast<TH1*>(fOutputList->At(i));
- if (h1){
- h1->Sumw2();
- continue;
- }
- THnSparse *hn = dynamic_cast<THnSparse*>(fOutputList->At(i));
- if (hn){
- hn->Sumw2();
- }
- }
- TH1::AddDirectory(oldStatus);
-
- PostData(1, fOutputList);
+ if(fbJetsRp){
+ // cent : nInpTrks : rp bins : rp wrt jet(2x) : probe pT : probe area: deltaPt : rp phi : rho : correction for RP : local rho
+ /*
+ entries = 1<<0 | 1<<1 | 1<<3 | 1<<4 | 1<<5 | 1<<6 | 1<<12 | 1<<14 | 1<<2 | 1<<22 | 1<<23 | 1<<24 | 1<<25;
+ opt = 1<<4 | 1<<5;*/
+ // cent : nInpTrks : rp bins : rp wrt jet(2x) : probe pT : deltaPt : rp phi : rho : correction for RP | pT hard bin
+ entries = 1<<0 | 1<<1 | 1<<3 | 1<<4 | 1<<5 | 1<<6 | 1<<14 | 1<<2 | 1<<22 | 1<<23 | 1<<26;
+ opt = 1<<4 | 1<<5;
+ fhnJetsRp = NewTHnSparseF("fhnJetsRp", entries, opt);
+ }
+
+ // cent : nInpTrks : rp bins: deltaPt : jetPt(2x) : deltaArea : pT hard bin (hr delta pt)
+ if(fbJetsDeltaPt){
+ entries = 1<<0 | 1<<1 | 1<<3 | 1<<14 | 1<<6 | 1<<7 | 1<<18 | 1<<26;
+ opt = 1<<1 | 1<<14 | 1<<6 | 1<<7 | 1<<18;
+ fhnJetsDeltaPt = NewTHnSparseF("fhnJetsDeltaPt", entries, opt);
+ }
+
+ // cent : nInpTrks : rp bins : deltaPt : jetPt(2x) : deltaR : deltaEta : jetEta(2x) : pT hard bin (hr for eta)
+ if(fbJetsEta){
+ entries = 1<<0 | 1<<1 | 1<<3 | 1<<14 | 1<<6 | 1<<7 | 1<<17 | 1<<15 | 1<<8 | 1<<9 | 1<<26;
+ opt = 1<<15 | 1<<8 | 1<<9;
+ fhnJetsEta = NewTHnSparseF("fhnJetsEta", entries, opt);
+ }
+
+ // cent : nInpTrks : rp bins : jetPt(2x) : jetPhi(2x) : deltaPt : deltaPhi : pT hard bin
+ if(fbJetsPhi){
+ entries = 1<<0 | 1<<1 | 1<<3 | 1<<6 | 1<<7 | 1<<10 | 1<<11 | 1<<14 | 1<<16 | 1<<26;
+ opt = 1<<10 | 1<<11;
+ fhnJetsPhi = NewTHnSparseF("fhnJetsPhi", entries, opt);
+ }
+
+ // cent : nInpTrks : rp bins : deltaArea : jetArea(2x) : deltaR : fraction : distance next rec jet : pT next jet : deltaPt : jetPt(2x) : pT hard bin (hr for area)
+ if(fbJetsArea){
+ entries = 1<<0 | 1<<1 | 1<<3 | 1<<18 | 1<<12 | 1<<13 | 1<<17 | 1<<19 | 1<<20 | 1<<21 | 1<<14 | 1<<6 | 1<<7 | 1<<26;
+ opt = 1<<18 | 1<<12 | 1<<13;
+ fhnJetsArea = NewTHnSparseF("fhnJetsArea", entries, opt);
+ }
+
+ // cent : nInpTrks : jetPt(3x) : deltaPt : delta : jetArea(3x) : fraction(2x) : deltaR(1x) : pT hard bin
+ if(fbJets3Branches){
+ entries = 1<<0 | 1<<1 | 1<<6 | 1<<7 | 1<<27 | 1<<14 | 1<<28 | 1<<12 | 1<<13 | 1<<29 | 1<<19 | 1<<30 | 1<<17 | 1<<26;
+ opt = 1<<6 | 1<<7 | 1<<27 | 1<<14 | 1<<28;
+ fhnJets3Branches = NewTHnSparseF("fhnJets3Branches", entries, opt);
+ }
+
+
+
+
+ //before cut
+
+ // cent : nInpTrks : rp bins : fraction : jetPt(2x) : jetEta(2x) : jetPhi(2x) (low resolution) (with fraction, eta, phi, pt cuts possible)
+ if(fbJetsBeforeCut1){
+ entries = 1<<0 | 1<<1 | 1<<3 | 1<<19 | 1<<6 | 1<<7 | 1<<8 | 1<<9 | 1<<10 | 1<<11;
+ opt = 0;
+ fhnJetsBeforeCut1 = NewTHnSparseF("fhnJetsBeforeCut1", entries, opt);
+ }
+
+ // cent : nInpTrks : rp bins : deltaPt : jetPt(2x) : deltaR : deltaEta : jetEta(2x) (low resolution)
+ if(fbJetsBeforeCut2){
+ entries = 1<<0 | 1<<1 | 1<<3 | 1<<14 | 1<<6 | 1<<7 | 1<<17 | 1<<15 | 1<<8 | 1<<9;
+ opt = 0;
+ fhnJetsBeforeCut2 = NewTHnSparseF("fhnJetsBeforeCut2", entries, opt);
+ }
+
+ fOutputList->Add(fHistEvtSelection);
+ fOutputList->Add(fHistJetSelection);
+ fOutputList->Add(fh2JetSelection);
+ if(fbEvent) fOutputList->Add(fhnEvent);
+ if(fbJetsMismatch1) fOutputList->Add(fhnJetsMismatch1);
+ if(fbJetsMismatch2) fOutputList->Add(fhnJetsMismatch2);
+ if(fbJetsRp) fOutputList->Add(fhnJetsRp);
+ if(fbJetsDeltaPt) fOutputList->Add(fhnJetsDeltaPt);
+ if(fbJetsEta) fOutputList->Add(fhnJetsEta);
+ if(fbJetsPhi) fOutputList->Add(fhnJetsPhi);
+ if(fbJetsArea) fOutputList->Add(fhnJetsArea);
+ if(fbJets3Branches) fOutputList->Add(fhnJets3Branches);
+ if(fbJetsBeforeCut1) fOutputList->Add(fhnJetsBeforeCut1);
+ if(fbJetsBeforeCut2) fOutputList->Add(fhnJetsBeforeCut2);
+
+ // =========== Switch on Sumw2 for all histos ===========
+ for (Int_t i=0; i<fOutputList->GetEntries(); ++i) {
+ TH1 *h1 = dynamic_cast<TH1*>(fOutputList->At(i));
+ if (h1){
+ h1->Sumw2();
+ continue;
+ }
+ THnSparse *hn = dynamic_cast<THnSparse*>(fOutputList->At(i));
+ if (hn){
+ hn->Sumw2();
+ }
+ }
+ TH1::AddDirectory(oldStatus);
+
+ PostData(1, fOutputList);