GFDstarHistsAnalysisMc

class GFDstarHistsAnalysisMc : public GFDstarHistsAnalysis

    protected:

         TH1* CombinedAccept(const char* variable)
         TH1* CreateAccept(TH1* histGenAcc, TH1* histGen, const char* variable, const char* titlePrefix)
         TH1* CreateHadCorr(TH1* hGen, TH1* hParton, const char* var, const char* titlePrefix)
         TH1* CreatePurity(TH1* hAllRec, TH1* hGenRec, const char* var, Bool_t fagHadCorr)
         TH1* CreateRecEff(const char* var, TH2* hRec2D, TH1* hRecRef, TH1* hGen, const char* titlePrefix)
         TH1* CreateStability(TH1* hGenRec, TH1* hGenNotRec, const char* variable, Bool_t fagHadCor)
         TH1* CreateTrigEff(TH1* histRec, TH1* histTrig, Int_t TE, const char* var, const char* titlePrefix)
         void Draw(void* funcPtr)
         void InitMC(const char* period)
      TArrayD TotalAccept(Double_t sumAcc, Double_t nGen, const char* prefix)
      TArrayD TotalRecEff(TH1* hRec, Double_t numGen, const char* prefix, TH1* hRecRef = 0)
      TArrayD TotalTrigEff(Double_t nTrig, Double_t nRec, Double_t nEffRec, Double_t nFillRec, Int_t TE, const char* prefix)
         TH1* XSection(TH1* h, const char* var, const char* titPrefix)
         TH1* XSectionRecInt(TH1* hNDstar, TH1* hRecEff, TH1* hAccept, TH1* hXSecGen, TH1* hGenNoAccW, TH1* hGen)

    public:

                             GFDstarHistsAnalysisMc(const char* mcFile)
                             GFDstarHistsAnalysisMc(const char* prefix, const char* period)
                             GFDstarHistsAnalysisMc(const GFDstarHistsAnalysisMc&)
                     virtual ~GFDstarHistsAnalysisMc()
              static TClass* Class()
                        TH1* CreateAccept(const char* variable)
                        TH1* CreateAccept(const char* var1, const char* var2, Int_t firstBin2, Int_t lastBin2)
                        TH1* CreateAcceptDsDiJet(const char* var)
                        TH1* CreateAcceptDsJet(const char* var, const char* regio)
                        TH1* CreateDetAcceptance(const char* variable)
                        TH1* CreateDsDiJetPurity(const char* var)
                        TH1* CreateDsDiJetStability(const char* var)
                        TH1* CreateDsJetPurity(const char* var)
                        TH1* CreateDsJetPurityHadCorr(const char* variable)
                        TH1* CreateDsJetStability(const char* variable)
                        TH1* CreateDsJetStabilityHadCorr(const char* var)
                        TH1* CreateHadCorr(const char* var, const char* forwBack)
                        TH1* CreateHadCorrDiJet(const char* var)
                        TH1* CreateRecEff(const char* variable)
                        TH1* CreateRecEff(const char* var1, const char* var1, Int_t firstBin2, Int_t lastBin2)
                        TH1* CreateRecEffDs1Jet(const char* var, const char* forwBack, UInt_t dmRefFlag = 2)
                        TH1* CreateRecEffDsDiJet(const char* var, UInt_t dmRefFlag = 2)
                        TH1* CreateTrigEff(const char* variable, Int_t TE = -1)
                        TH1* CreateTrigEff(const char* var1, const char* var1, Int_t firstBin2, Int_t lastBin2, Int_t TE = -1)
                        TH1* CreateTrigEffDsDiJet(const char* var, Int_t TE = -1)
                        TH1* CreateTrigEffDsJet(const char* variable, const char* region, Int_t TE = -1)
                        TH1* DrawDsDiJetGenVsRec(const char* var = "0")
                        TH1* DrawDsJetGenVsRec(const char* var = "0", Bool_t hadGen = kFALSE)
                        TH1* DrawDsJetQuark(const char* var, const char* dsOrJet, Int_t mode = 0)
                        TH1* DrawDsJetRes(const char* var, const char* rel)
                        TH2* DrawEffCorrWeights(const char* var)
                        void DrawEflow(const char* region = "0")
                        void DrawJetQuark(const char* jet = "jet1")
                        void DrawRecEffAll()
                    Double_t GetBranchRatio() const
                        TH1* GetHistAccept(const char* variable)
                        TH1* GetHistGen(const char* variable)
                        TH1* GetHistGenAccept(const char* variable)
                        TH1* GetHistGenNoAccWeight(const char* variable)
                        TH1* GetHistGenWeight(const char* variable)
                        TH1* GetHistNonTrig(const char* variable)
                        TH1* GetHistTrig(const char* variable, Int_t TE = -1)
             virtual TClass* IsA() const
                virtual void ShowMembers(TMemberInspector& insp, char* parent)
                virtual void Streamer(TBuffer& b)
                        void StreamerNVirtual(TBuffer& b)
                     TArrayD TotalAccept()
                     TArrayD TotalAcceptDsDiJet()
                     TArrayD TotalAcceptDsJet(const char* region)
                     TArrayD TotalHadCorr(const char* forwBack)
                     TArrayD TotalHadCorrDiJet()
                     TArrayD TotalRecEff()
                     TArrayD TotalRecEffDsDiJet(UInt_t dmRefFlag = 2)
                     TArrayD TotalRecEffDsJet(const char* region, UInt_t dmRefFlag = 2)
                     TArrayD TotalTrigEff(Int_t TE = -1)
                     TArrayD TotalTrigEffDsDiJet(Int_t TE = -1)
                     TArrayD TotalTrigEffDsJet(const char* region, Int_t TE = -1)
                        TH1* XSection(const char* variable, Bool_t acceptanceWeighted = kFALSE)
                        TH1* XSection(const char* var1, const char* var2, Int_t firstBin2, Int_t lastBin2)
                        TH1* XSectionDsDiJet(const char* var)
                        TH1* XSectionDsJet(const char* var, const char* forwBack)
                        TH1* XSectionDsJetRec(const char* var, const char* foBack)
                        TH1* XSectionEta(Bool_t acceptanceWeighted = kFALSE)
                        TH1* XSectionPhi(Bool_t acceptanceWeighted = kFALSE)
                        TH1* XSectionPt(Bool_t acceptanceWeighted = kFALSE)
                        TH1* XSectionRec(const char* variable)
                     TArrayD XSectionTot(Bool_t acceptanceWeighted = kFALSE)
                     TArrayD XSectionTotDsDiJet()
                     TArrayD XSectionTotDsJet(const char* forwBack)
                        TH1* XSectionTotHist(Bool_t draw = kFALSE, Bool_t acceptanceWeighted = kFALSE)
                        TH1* XSectionWgammaP(Bool_t acceptanceWeighted = kFALSE)

Data Members

    public:

      static const Double_t kHighestDmMc     up to this value dm bins taken into account for N(D*)
      static const Double_t kBranchRatio     BR(D*->Kpipis) from PDG
      static const Double_t kBranchRatioMc   BR(D*->Kpipis) as used in H1MC
      static const Double_t kXSecCorrFactor  due to wrong c->D* branching

Class Description

 GFDstarHistsAnalysisMc

   Author:      Gero Flucke
   Date:        June 2nd, 2002
   last update: $Date: 2005/11/16 00:39:02 $
   by:          $Author: flucke $

     fMultiDmFitter->SetFitOrNotFlag(1);
  4th digit:    Mean
  3rd :         Sigma
  2nd :         BackgrExp
  1st :         BackgrSqr
  so
     mode = 1011 means: fix Mean, Sigma and BackgrSqr,
                        but fit BackgrExp
                       (and fit BackgrNorm, NDstar)

 nothing to be done

 var: Pt/Eta+Ds/Jet or Deta, Dphi, Dpt
 region: 'Back', 'Forw' or ''
 TE: 19, 31 or -1 (=default: 19&&31)

 var: Pt/Eta + DsJet/OthJet/Ds, DetaDijet, DphiDijet, PtDijetDijet, CosThDijet, xGamDijet,
      DeltaR, DeltaRNoDjet
 TE: 19, 31 or -1 (=default: 19&&31)

pt,eta,phi,wGammaP,...

 var is PtDs, PtJet,EtaJet, EtaDs, Deta, Dphi etc.

 for D* jet

 for D* dijet

 rec eff for D* + 1 jet
 var: Pt/Eta+Ds/Jet (e.g. PtDs), Dpt, Deta, Dphi
 forwBack: 'Forw', 'Back', '' (not NULL)
 if(dmRefFlag != 0): First fit a reference hist with all params free, then fit with fFitMode
 >=2 [default]: D* inclusive
 1          : D* + 1 jet inclusive

 rec eff for D* hists: pt, eta, phi, wGammaP, ...

 hGen:   hist with generator level distribution
 hRec2D: hist with reconstruction level distribution on y-axis and delta(m) on x-axis
 var:    name of variable

 pt, eta, phi, wGammaP, etc. for D*

 var: Pt/Eta+Ds/Jet or Deta, Dphi, Dpt
 region: 'Back', 'Forw' or ''

 var: Pt/Eta + DsJet/OthJet/Ds, Deta, Dphi, PtDijet, CosTh, xGam,
      DeltaR, DeltaRNoDjet

 histGenAcc: generator level, weighted with tagger acceptance
 histgen: generator level

 total cross section for D* + 1 jet
 forwback: Forw, Back, ''

 total cross section for D* dijet

   fHist1D = acceptanceWeighted ? this->GetHistGenAccept("pt") : this->GetHistGen("pt");

   fHist1D = acceptanceWeighted ? this->GetHistGenAccept("eta") :this->GetHistGen("eta");

   fHist1D = acceptanceWeighted ? this->GetHistGenAccept("phi") : this->GetHistGen("phi");

pt,eta,phi,wGammaP,y,nJet,eta/ptJet1/2,xGamma

 hNDstar already corrected for binwidth...!

 Forw, ForwBar, Barrel, Spacal, "" (sum) or NULL: all 3+1 regions;

 jet1, jet2

Dphi Deta xGam PtJet PtDs EtaJet

Dphi Deta xGam PtJet EtaJet

Dphi/Deta/xGam/Pt/Pt2 + Dijet Pt/Eta+Ds/DJetOthJet

 purity as hGenRec/hAllRec, some

Dphi Deta xGam PtJet PtDs EtaJet

xGam Deta Dphi EtaJet PtJet PtJet

Dphi/Deta/xGam/Pt/Pt2 + Dijet Pt/Eta+Ds/DJetOthJet

 variable without 'Res', 'ResRel', if Rel needed: rel = "Rel"
 EtaJetRes EtaJetPt<n>Res PtJetEta<n>ResRel PtJetResRel
 PtDsJetResRel PtDsResRel

 return last drawn hist

 return last drawn hist

 upper case Pt, Eta...

  returning effRec and its error
 cf. log book 13.9.2002

 if dmRefFlag == 2 (default): reference fit from inclusive D*

  returning effRec and its error
 cf. log book 13.9.2002

 D*: 19, 31, -1 (= 19 && 31)

 D* + jet: returning effTrig and its error
 Back, Forw, ''; 19, 31, default: both

 D* dijet: returning effTrig and its error
 Back, Forw, ''; 19, 31, default: both

 returning effTrig and its error:
 nTrig: trig & rec, nRec: all rec, nEffRec: effective number entries for nRec
 nFillRec: hist fill actions for nRec
 if(nEffRec/nRec >= 0.7): do asymm. binomial error GFMath::BinomialError
 then result[2/3] are the asymmetric errors
 result[1]: symmetric easy binomial error

  returning acceptance (and - not yet implemented - its error)

 sumAccept: sum of acceptances of generated events, numGen: num of generated events
 error not implemented!

  D* + 1 jet: returning acceptance (and - not yet implemented - its error)
 region: Forw, Back, ''

  D* dijet: returning acceptance (and - not yet implemented - its error)

Inline Functions

                      Double_t GetBranchRatio() const
                       TClass* Class()
                       TClass* IsA() const
                          void ShowMembers(TMemberInspector& insp, char* parent)
                          void Streamer(TBuffer& b)
                          void StreamerNVirtual(TBuffer& b)
        GFDstarHistsAnalysisMc GFDstarHistsAnalysisMc(const GFDstarHistsAnalysisMc&)