GFDstarAnalyse

class description - source file - inheritance tree

class GFDstarAnalyse : public H1RunDependent

 
    This is an abstract class, constructors will not be documented.
    Look at the header to check for available constructors.

    protected:

                  void CreateHistsBgFinder()
                  void CreateHistsD0Checks()
          virtual void CreateHistsJetMult()
                  void CreateRunDepAnalog(TH1* input, GFHistArray* output, Int_t trigger)
                Bool_t FillHistsDstar(const H1PartDstar* dstar, GFHistArray* dmHists, GFHistArray* dmHists83, GFHistArray* dmHists84, Bool_t requireL4)
                  void FillHistsDstar1Jet(const H1Part* ds, const H1PartJet* jet, GFHistArray* hists, Double_t weight, Double_t dm)
                  void FillHistsDstarDiJet(const H1Part* ds, const H1PartJet* dsJet, const H1PartJet* oJet, Bool_t realDsJet, GFHistArray* hists, Double_t weight, Double_t dm)
                  void FillHistsJetComp(GFHistArray* hists, const H1PartJet* jet, Double_t dm, Double_t weight)
                  void FillHistsTracks(const H1PartDstar* dstar, GFHistArray* toFill, Double_t w)
      virtual Double_t GetCmEnergySqr() const
              Double_t GetSumEt(Int_t threshold = 1) const
      virtual Double_t GetWeight(Int_t trigger, const H1PartDstar* ds) const
       virtual Float_t GetWgammaP(Int_t subTrig) const
       virtual Float_t GetWgammaP33() const
       virtual Float_t GetWgammaP44() const
       virtual Float_t GetWgammaPjb() const
              Double_t GetXgamma(const H1Part* part1, const H1Part* part2, Int_t trig) const
              Double_t GetXgamma(const H1Part* part1, const H1Part* part2, Double_t y) const
              Double_t GetXgammaMass(const H1PartDstar* dstar, const H1PartJet* jet, Double_t yHfs) const
              Double_t GetXgammaMassDiJet(const H1PartDstar* dstar, const H1PartJet* dsJet, const H1PartJet* othJet, Double_t yHfs) const
      virtual Double_t GetXgammaMassDiJetHelp(const H1Part* dstar, const H1PartJet* dsJet, const H1PartJet* othJet) const
      virtual Double_t GetXgammaMassHelp(const H1Part* dstar, const H1PartJet* jet) const
      virtual Double_t GetXpMass(const H1PartDstar* dstar, const H1PartJet* jet) const
      virtual Double_t GetXpMassHelp(const H1Part* dstar, const H1PartJet* jet) const
       virtual Float_t GetY(Int_t subTrig) const
       virtual Float_t GetY33() const
       virtual Float_t GetY44() const
       virtual Float_t GetYjb() const
      virtual Double_t GetZ(const H1PartDstar* ds) const
           TDirectory* MakeHistDir(const char* name, const char* title)
                Bool_t OneL4ClassGiven(const TArrayI& classes) const

    public:

                   virtual ~GFDstarAnalyse()
            virtual Bool_t AreTracksOk(const H1PartDstar* dstar) const
          virtual Double_t Associate(const H1Part* ref1, const H1Part* ref2, H1PartJet*& outCloseTo1, H1PartJet*& outCloseTo2)
                  Double_t CalcSumEt(Double_t thetaCutDegree) const
            static TClass* Class()
              GFHistArray* CreateHists(const char* nameAdd, const char* title, Int_t subTrig, UInt_t dimension, TDirectory* dir = 0)
                      void CreateHistsCentralNotCombinedTrack()
              virtual void CreateHistsD0()
              virtual void CreateHistsDm(GFHistArray*& dmHist, GFHistArray*& dmHist83, GFHistArray*& dmHist84, const TString& nameAdd)
              GFHistArray* CreateHistsDstar1Jet(Int_t forwbackFlag, Int_t dim) const
              virtual void CreateHistsDstarDiJet()
              GFHistArray* CreateHistsDstarDiJetInt(Int_t dim) const
              virtual void CreateHistsDstarJets()
                      void CreateHistsEventClass()
                      void CreateHistsJetComp()
                      void CreateHistsJetProf()
              virtual void CreateHistsL4Eff2()
              virtual void CreateHistsSumEt()
                      void CreateHistsTest()
                      void CreateHistsTrackClus()
                      void CreateHistsTrackErr()
                      void CreateHistsTrackKind()
              virtual void CreateHistsTracks()
                      void CreateHistsTriggers()
                    Bool_t D0CutLoose(const H1PartDstar* ds) const
            virtual Bool_t DstarCutLoose(const H1PartDstar* ds) const
            virtual Bool_t DstarCutLooseNoTrackCut(const H1PartDstar* ds) const
              virtual void EndJob()
              virtual void EndRun(Int_t run)
                      void FillHists(GFHistArray* hists, const H1PartDstar* ds, Int_t subTrig, Double_t weightExtern = -1.)
                      void FillHistsBgFinder(const H1PartDstar* dstar)
                      void FillHistsCentralNotCombinedTrack(const H1PartDstar* dstar)
                      void FillHistsD0(const H1PartDstar* dstar)
                      void FillHistsD0Checks(const H1PartDstar* ds)
                    Bool_t FillHistsDstar(const H1PartDstar* dstar)
                      void FillHistsDstarDiJet(const H1PartDstar* ds)
                      void FillHistsDstarJets(const H1PartDstar* ds)
                      void FillHistsEventClass(const H1PartDstar* dstar)
                      void FillHistsJetComp(const H1PartDstar* ds)
              virtual void FillHistsJetMult(const H1PartDstar* ds, Bool_t accessPartCand)
                      void FillHistsJetProf(const H1PartDstar* ds)
                      void FillHistsL4Eff2(const TObjArray* dstars)
                      void FillHistsSumEt(const H1PartDstar* dstar, Double_t sumEt = 0.)
              virtual void FillHistsTest(const H1PartDstar* ds)
                      void FillHistsTrackClus()
                      void FillHistsTrackErr(const TObjArray* dstars)
                      void FillHistsTrackKind(const H1PartDstar* dstar)
                      void FillHistsTracks(const H1PartDstar* dstar)
                      void FillHistsTriggers(const H1PartDstar* dstar)
              virtual void FillHistsYjbRes(const H1PartDstar* ds)
      static const H1Part* GetCloser(const H1Part* reference, const H1Part* first, const H1Part* second, Float_t* dist = 0) const
                  Double_t GetCutPt1()
                  Double_t GetCutPt2()
                    Bool_t GetDiJets(const H1PartDstar* dstar, const H1PartJet*& dsJet, const H1PartJet*& oJet) const
                      void GetEtRatios(const H1PartJet* jet, Double_t& ratioTrack, Double_t& ratioEmClus, Double_t& ratioHadClus)
              H1EventList* GetEventList()
                     Int_t GetHighestPt(const TObjArray& parts) const
                 TObjArray GetJets(const H1PartDstar* ds, Int_t* dsJet = 0) const
                   Float_t GetNLH(const H1PartDstar* dstar, Int_t part) const
                     Int_t GetNonDstarJets(const H1PartDstar* dstar, TObjArray& jets) const
            const TString* GetPeriod() const
                  Double_t GetRecEnergies(Double_t& forward, Double_t& forwardBar, Double_t& barrel, Double_t& spacal) const
           virtual TClass* IsA() const
                    Bool_t IsDedxOK(const H1PartDstar* dstar, Int_t nDaughter = -1) const
            virtual Bool_t IsL4Found(Int_t hqsel45bit, const H1PartDstar* ds, Int_t hqselMode = -2) const
                    Bool_t IsL4Found(Int_t hqsel45Bit, Int_t hqselMode, const H1BankPtr& yrclYecl) const
            virtual Bool_t IsL4Ref(Int_t subtrig) const
            virtual Bool_t IsS83(const H1PartDstar* ds = 0) const
            virtual Bool_t IsS83NoL4Check(const H1PartDstar* dstar) const
            virtual Bool_t IsS84(const H1PartDstar* ds = 0) const
            virtual Bool_t IsS84NoL4Check(const H1PartDstar* dstar) const
            virtual Bool_t IsSignalRegionDm(const H1PartDstar* ds) const
            virtual Bool_t IsSumEt(const H1PartDstar* dstar) const
            virtual Bool_t IsSumPtKPiOK(const H1PartDstar* ds, Double_t ptCut) const
            virtual Bool_t IsZvtxCut(Float_t cut = 0.) const
                    Bool_t JetContains(const H1PartJet* jet, const H1PartDstar* ds) const
                    Bool_t JetCut(const H1PartJet* jet) const
                      void OpenHat(Bool_t yesNo = kTRUE)
       virtual TDirectory* OpenHistsFile(const char* fileName = "0")
                      void OpenMods(Bool_t yesNo = kTRUE)
                      void OpenOds(Bool_t yesNo = kTRUE)
              virtual void OpenTree()
              virtual void PrepareEvent()
                      void Print(const H1PartDstar* dstar) const
             virtual Int_t RemoveLowPt(TObjArray* parts, Int_t nStay) const
             virtual Int_t RemoveLowPt(TObjArray* parts, Double_t minPt) const
                    Bool_t SetIgnoreDedx(Bool_t ignore)
                    Bool_t SetIgnoreSumEt(Bool_t ignore)
                    Bool_t SetIgnoreSumPtKPi(Bool_t ignore)
                      void SetJetFlag(Int_t flag)
              virtual void ShowMembers(TMemberInspector& insp, char* parent)
              virtual void StartJob()
              virtual void StartRun(Int_t run)
              virtual void Streamer(TBuffer& b)
                      void StreamerNVirtual(TBuffer& b)

Data Members


    protected:

      GFSteerDstarAnalyse* fSteer                         steering object containing cuts and binning
                    Bool_t fIgnoreDedx                    
                    Bool_t fIgnoreSumEt                   
                    Bool_t fIgnoreSumPtKPi                
                  Double_t fDsCutPt1                      used for S84 from 24/9/2002 on
                  Double_t fDsCutPt2                      used for S83 from 24/9/2002 on
                  Double_t fEventCutyMinS83               
                  Double_t fEventCutyMaxS83               
                  Double_t fEventCutyMinS84               
                  Double_t fEventCutyMaxS84               
                   TArrayD fWgammaP33Bins                 
                   TArrayD fWgammaP44Bins                 
                   TArrayD fY33Bins                       
                   TArrayD fY44Bins                       
                   Float_t fXgamProfile                   x_gamma border for profile hists
              GFHistArray* fDmHists                       
              GFHistArray* fDmHistsS83                    
              GFHistArray* fDmHistsS84                    
              GFHistArray* fD0CheckHists                  
              GFHistArray* fD0CheckHistsS83               
              GFHistArray* fD0Hists                       
              GFHistArray* fD0HistsS83                    
              GFHistArray* fD0HistsS84                    
              GFHistArray* fTrackHists                    
              GFHistArray* fTrackHistsS83                 
              GFHistArray* fTrackHistsS84                 
              GFHistArray* fSumEtHistsS83                 
              GFHistArray* fSumEtHistsS84                 
              GFHistArray* fTriggersHists                 
              GFHistArray* fBgFinderHists                 
              GFHistArray* fL4YesHistsS83                 L4 class
              GFHistArray* fL4YesHistsS84                 
              GFHistArray* fL4NoHistsS83                  
              GFHistArray* fL4NoHistsS84                  
              GFHistArray* fL4YesNotDsHistsS83            
              GFHistArray* fL4YesNotDsHistsS84            
              GFHistArray* fL4YesDstarHistsS83            L4 dstar finder in gamma p
              GFHistArray* fL4YesDstarHistsS84            
              GFHistArray* fL4NoDstarHistsS83             
              GFHistArray* fL4NoDstarHistsS84             
              GFHistArray* fEventClassHists               
              GFHistArray* fEventClassHistsS83            
              GFHistArray* fEventClassHistsS84            
              GFHistArray* fOpenCharmFindHists            
              GFHistArray* fTrackKindHists                
              GFHistArray* fCentralNotCombinedTrackHists  
              GFHistArray* fJet1ProfileHists              
              GFHistArray* fJet1ProfileHistsS83           
              GFHistArray* fJet1ProfileHistsS84           
              GFHistArray* fJet2ProfileHists              
              GFHistArray* fJet2ProfileHistsS83           
              GFHistArray* fJet2ProfileHistsS84           
              GFHistArray* fJetCompHists                  
              GFHistArray* fJetCompHistsS83               
              GFHistArray* fJetCompHistsS84               
                 TObjArray fJetProfiler                   
                 TObjArray fJetProfilerS83                
                 TObjArray fJetProfilerS84                
              GFHistArray* fTrackClusHists                
              GFHistArray* fTrackClusHistsForw            
              GFHistArray* fTrackClusHistsForwBar         
              GFHistArray* fTrackClusHistsBar             
              GFHistArray* fTrackClusHistsSpac            
              GFHistArray* fTrackErrHists                 
              GFHistArray* fDstar1JetHists                
              GFHistArray* fDstar1JetHistsS83             
              GFHistArray* fDstar1JetHistsS84             
              GFHistArray* fDstarBack1JetHists            
              GFHistArray* fDstarBack1JetHistsS83         
              GFHistArray* fDstarBack1JetHistsS84         
              GFHistArray* fDstarForw1JetHists            
              GFHistArray* fDstarForw1JetHistsS83         
              GFHistArray* fDstarForw1JetHistsS84         
              GFHistArray* fDstarDiJetHists               
              GFHistArray* fDstarDiJetHistsS83            
              GFHistArray* fDstarDiJetHistsS84            
              GFHistArray* fJetMultHistsS83               
              GFHistArray* fTestHists                     
                 TObjArray fAllHistArrays                 
                    TFile* fHistFile                      
                  TString* fPeriod                        
                    Bool_t fOpenHat                       flags to open chains or not
                    Bool_t fOpenMods                      dito
                    Bool_t fOpenOds                       dito
                     Int_t fJetFlag                       flag to use jets with (0) or without (1) D* as input
              H1EventList* fList                          

    public:

         static const Int_t kDmHistNumBins    
      static const Double_t kDmHistLowEdge    
      static const Double_t kDmHistHighEdge   
         static const Int_t kD0HistNumBins    
       static const Float_t kD0HistLowEdge    
       static const Float_t kD0HistHighEdge   
      static const Double_t kCmEnergySqr      s = 4*E_e*E_p 
         static const Int_t kNumEventClass    
         static const Int_t kL4FinderBit      
         static const Int_t kNumMaxWarn       max number of warnings printed where limit foreseen
         static const Int_t kNumDefaultHists  num hists handled by CreateHists and FillHists!

See also

GFDstarAnalyseData, GFDstarAnalyseMc

Class Description

 GFDstarAnalyse
   Author:      Gero Flucke
   Date:        May 31st, 2002
   last update: $Date: 2006/01/11 11:33:27 $
   by:          $Author: flucke $
~GFDstarAnalyse()
void StartJob()
void StartRun(Int_t run)
void EndRun(Int_t run)
void EndJob() 

void PrepareEvent()



 do things that should manipulate the event and should stay for the whole event
 e.g. calibrate jets!
 nothing done currently...


void CreateHistsDm(GFHistArray*& dmHists, GFHistArray*& dmHistsS83, 
				   GFHistArray*& dmHistsS84, const TString& nameAdd)



 foreseen is nameAdd = "Wc" or ""


void CreateHistsD0()




void CreateHistsTracks()




void CreateHistsSumEt()




void CreateHistsTriggers()




void CreateHistsEventClass()




void CreateHistsTrackKind()




void CreateHistsCentralNotCombinedTrack()




void CreateRunDepAnalog(TH1* input, GFHistArray* output, Int_t trigger)




GFHistArray* CreateHists(const char* nameAdd, 
					      const char* title, 
					      Int_t subTrig, UInt_t dimension, TDirectory* dir)



 create a GFHistArray of hists with identifier 'nameAdd', 'title' and
 for 'subTrig' (if subtrig == -1: subTrig is not put to name,
                                    binning is mixed!)
 'dimension' is one or two so far:
 1: TH1 with x-axis for pt, eta, phi, ptSumEt, wGammaP, y, nJet, ptJet1/2, etaJet1/2, xGamma
 2: x-axis is for delta m, y axes are for pt, eta, phi, ptSumEt, wGammaP, etc.
 if(dir != NULL) dir->cd() before creating hists
 which hists are created has to be in sync with 'FillHists' and ...Mc::FillHistsGen!


void FillHists(GFHistArray* hists, const H1PartDstar* ds,
			       Int_t subTrig, Double_t weightExtern)



 filling histograms in hists with weight by GetWeight(subTrig):
 1st with pt, 2nd with eta, 3rd with phi of 'ds'
 4th with pt(D*)/ sum(E_t)
 5th with wgammap for corresponding 'subTrig'
 (83 or 84, else 5th hist ignored!).
 If a hist is 2D: x-axis will be filled with ds->GetDm().
 if weightExtern == -1. (default), use this->GetWeight(subTrig) as weight, else weightExtern


Bool_t FillHistsDstar(const H1PartDstar* dstar)



 returns true if dstar survives all cuts (only looser pt and eta cuts neccesary)
 (subtrigger don't care for return value!)


Bool_t FillHistsDstar(const H1PartDstar* dstar,
				      GFHistArray* dmHists,
				      GFHistArray* dmHistsS83,
				      GFHistArray* dmHistsS84, Bool_t requireL4)



 returns true if dstar survives all cuts (only looser pt and eta cuts neccesary)
 (subtrigger don't care for return value!)
 first argument may be NULL, if the others also ==> crash!


void FillHistsTracks(const H1PartDstar* dstar)



 all track quantities are filled!


void FillHistsTracks(const H1PartDstar* dstar, 
				     GFHistArray* trackHists, Double_t w)// const



 fill with weight 'w' (except dedx scatters... )


void FillHistsD0(const H1PartDstar* dstar)




void FillHistsTriggers(const H1PartDstar* dstar)



 assuming dstar fulfills cuts that you want
 only pt2 cut for S83 and eta cut for s84!


void CreateHistsBgFinder()




void FillHistsBgFinder(const H1PartDstar* dstar)




void FillHistsEventClass(const H1PartDstar* dstar)



 assuming dstar fulfills cuts that you want
 only pt2 cut for S83 / eta for S84!


void FillHistsSumEt(const H1PartDstar* dstar, Double_t sumEt)



 if sumEt == 0. (default) take it from H1FloatPtr newSumEt("SumEt"), entry [1]


void FillHistsTrackKind(const H1PartDstar* dstar)




void FillHistsCentralNotCombinedTrack(const H1PartDstar* dstar)




void CreateHistsL4Eff2()




void FillHistsL4Eff2(const TObjArray* dstars)




void CreateHistsJetProf()




void FillHistsJetProf(const H1PartDstar* ds)




void CreateHistsTrackClus()




void CreateHistsTrackErr()




void FillHistsTrackErr(const TObjArray *dstars)




void FillHistsTrackClus()




void CreateHistsJetComp()




void FillHistsJetComp(const H1PartDstar* ds)




void FillHistsJetComp(GFHistArray* hists, const H1PartJet* jet, 
				      Double_t dm, Double_t weight)




void CreateHistsDstarJets()




GFHistArray* CreateHistsDstar1Jet(Int_t forwBack, Int_t dim) const



 creating and fill array for D* + 1 jet hists, hist names
 (parallel to FillHistsDstar1Jet)
 - start with 'dstar1Jet',
 - then the variable follows ('', 'Pt/Eta'+'Ds/Jet', 'Dphi', 'Deta', 'Dpt')
 - then the identifier of eta-ordered set follows: 'Forw', 'Back', ''

 which of the last and the chosen binning depends on forwBack flag:
 -1: 'Back', i.e. eta(D*) < eta(jet)
  0: '', i.e. no matter
 +1: 'Forw, i.e. eta(D*) > eta(jet)

 if dim = 2: x-axis is delta m(D*), else 1D hists


void FillHistsDstarJets(const H1PartDstar* ds)



 hists for D* and (other) jet(s)


void FillHistsDstar1Jet(const H1Part* ds, const H1PartJet* jet,
					GFHistArray* hists, Double_t weight, 
					Double_t dm)



 fill hists with D* + 1 jet quantities with the weight
 if(dm>0) assuming 2D hist where x is filled with 'dm'
 (parallel to CreateHistsDstar1Jet)


void CreateHistsDstarDiJet()



 hists for FillHistsDstarDiJets(const H1PartDstar* ds)


GFHistArray* CreateHistsDstarDiJetInt(Int_t dim) const



 creating and fill array for D*-dijet hists, hist names
 (parallel to FillHistsDstarDiJet)
 - start with 'dstarDiJet',
 - then the variable follows ('', 'Pt/Eta'+'Ds/DJet/OthJet', 'Dphi', 'Deta', 'Dpt')

 if dim = 2: x-axis is delta m(D*), else 1D hists


void FillHistsDstarDiJet(const H1PartDstar* ds)



 dijet histos


void FillHistsDstarDiJet(const H1Part* ds, const H1PartJet* dsJet, 
					const H1PartJet* oJet, Bool_t realDsJet, GFHistArray* hists, Double_t weight, 
					Double_t dm)



 fill hists with D* -dijet quantities with the weight
 if(dm>0) assuming 2D hist where x is filled with 'dm'
 (parallel to CreateHistsDstarDiJet)


void CreateHistsD0Checks()




void FillHistsD0Checks(const H1PartDstar* ds)




void CreateHistsJetMult()




void FillHistsJetMult(const H1PartDstar* ds, Bool_t accessPartCand)




void CreateHistsTest()




void FillHistsTest(const H1PartDstar* ds)



   if(!fTestHists) this->CreateHistsTest();
   if(this->DstarCutLoose(ds)){
     const Double_t dm = ds->GetDm();
     if(dm >= fSteer->GetD0DeltaMLow() && dm < fSteer->GetD0DeltaMHigh()){
       const Double_t m = ds->GetMass();
       fTestHists->At(0)->Fill(m);
       if(this->IsS83(ds)) fTestHists->At(1)->Fill(m);
       if(this->IsS84(ds)) fTestHists->At(2)->Fill(m);
     }
   }


TDirectory* OpenHistsFile(const char* fileName)



 opens file, does file->cd() and return previous (!) gDirectory
 please do '->cd()' on return value not to change gDirectory!


TDirectory* MakeHistDir(const char* name, const char* title)



 make new directory in hist output file and cd into it
 the old directory is returned


Bool_t IsSignalRegionDm(const H1PartDstar* ds) const



 true if D* lies in SR of delta m distribution


Bool_t DstarCutLoose(const H1PartDstar* dstar) const



 all cuts common to both taggers (well, not really all...)


Bool_t DstarCutLooseNoTrackCut(const H1PartDstar* dstar) const




Bool_t D0CutLoose(const H1PartDstar* dstar) const




Bool_t AreTracksOk(const H1PartDstar* dstar) const



 all cuts on single track quantities
 cut on radial length (K, pi) and likelihoods (K, pi, pi_s)


Bool_t IsDedxOK(const H1PartDstar* dstar, Int_t nDaughter) const



 check whether dedx is OK:
 nDaughter < 0: check all K, pi, pis (default)
           = 0: check K only
           = 1: check pi only
           = 2: check pi_s only


Bool_t IsSumPtKPiOK(const H1PartDstar* ds, Double_t ptCut) const




Float_t GetNLH(const H1PartDstar* dstar, Int_t part) const



 get normalised likelihood taking into account K, pi and proton for a D* daughter:
 part: 0 = K, 1 = pi, 2 = pis (particle fluxes of pi, K and p assumed to be equal)
 1.09 if no measurement available, FLT_MAX if problems...


void Print(const H1PartDstar* dstar) const




Double_t CalcSumEt(Double_t thetaCutDegree) const




Double_t GetCmEnergySqr() const



 square of center of mass energy (e+p) [i.e. s ~= 4 * E_e * E_p]
 if H1Tree exists first try dstarUser tree for beam energies
 if nothing helps: take hardcoded default


Bool_t IsL4Found(Int_t hqsel45Bit, Int_t hqselMode, 
				 const H1BankPtr& yrclYecl) const



 true if hqsel45Bit is on in
 word hqselMode (mode 0: DIS, 1: very low Q2, 2: tagged, 3: untagged photoproduction,
    -1: any of them, -2: just check L4-class 15, ignore  hqsel45Bit and yrclYecl)
 of given bank yrclYecl (which is supposed to be a YECL [online] or YRCL [offline] bank)


Bool_t IsZvtxCut(Float_t cut) const // 0.: default!



 cutting on abs(zVtx) < cut
 cut == 0 (default) cut is set to fZvtxCut


Bool_t IsSumEt(const H1PartDstar* ds) const




Double_t GetSumEt(Int_t threshold) const



   if(!fOpenUserTree) return 0.;


Float_t GetYjb() const




Float_t GetY(Int_t subTrig) const




Float_t GetWgammaPjb() const




Float_t GetWgammaP(Int_t subTrig) const




Double_t GetZ(const H1PartDstar *ds) const




Bool_t OneL4ClassGiven(const TArrayI& classes) const




Double_t GetXgamma(const H1Part* part1, const H1Part* part2, 
				   Double_t y) const



 x_gamma from pt, eta, y


Double_t GetXgamma(const H1Part* part1, const H1Part* part2, 
				   Int_t trig) const




Double_t GetXgammaMass(const H1PartDstar* dstar,const H1PartJet* jet,
				       Double_t yHfs) const



 xgamma taking into account masses:
 [(E-pz)_{D*_HFS} + sum_{i=1}^{N(jet daught)}(E-pz)_{i}]/2*yHfs*E_{e}


Double_t GetXgammaMassDiJet(const H1PartDstar *dstar,
					    const H1PartJet *dsJet, const H1PartJet *othJet,
					    Double_t yHfs) const



 xgamma taking into account masses:
 (sum_{i=1}^{N(D*-jet daught)}(E-pz)_{i} + sum_{i=1}^{N(other jet daught)}(E-pz)_{i}]) 
   2*yHfs*E_{e}


Double_t GetXpMass(const H1PartDstar* dstar,const H1PartJet* jet) const



 SLOW: looping jet daughters!
 x_p taking into account masses:
 [(E+pz)_{D*_HFS} + sum_{i=1}^{N(jet daught)}(E+pz)_{i}]/2*E_{p}


Int_t GetHighestPt(const TObjArray& parts) const



 parts must contain H1Part, return index of that with highest pt,
 -1 if 'parts' empty


TObjArray GetJets(const H1PartDstar* ds, Int_t* dsJet) const



 get array of jets, input defined by fJetFlag
 output is ordered in eta: starting with negative eta


Int_t GetNonDstarJets(const H1PartDstar* dstar, 
					  TObjArray& jets) const



 array of jets passing the minimal criteria and NOT containing 'dstar'
 will be put into 'jets'
 ordered in eta: most backward first!

 returns the index of the jet, that comes after the D*
 if we put the D* into the list, too:

-1: no jets left
 0: D* more backwards than most backward jet
 1: D* between most backward and next more forward jet
 2: etc.


Bool_t GetDiJets(const H1PartDstar* dstar, const H1PartJet *&dsJet,
				 const H1PartJet *&oJet) const



 dsJet and oJet are set (possibly only one of them!)
 return kFALSE if dstar not in dsJet (i.e. dsJet is closer than oJet to D*)


Bool_t JetCut(const H1PartJet* jet) const//, Bool_t firstJet) const



 true if jet survives cuts, if 'firstJet': require higher pt-cut
 only pt and eta cuts


Bool_t JetContains(const H1PartJet* jet, const H1PartDstar* ds) const



 true if jet contains ds


Double_t GetRecEnergies(Double_t& forward, Double_t& forwardBar,
					Double_t& barrel, Double_t& spacal) const



 fills energie sums in detector region defined by some angle (cf.code)
 returns sum of all regions


void GetEtRatios(const H1PartJet* jet, Double_t& ratioTrack, 
					  Double_t& ratioEmClus, Double_t& ratioHadClus)




const H1Part* GetCloser(const H1Part* reference, const H1Part* first,
					const H1Part* second, Float_t* dist)



 static method
 returns that H1Part from first/second that is closer in eta phi to reference,
 if dist points to any memory it is filled with the distance


Int_t RemoveLowPt(TObjArray* parts, Int_t nStay) const



 remove the lowest pt particles in 'parts' such that 'nStay' are left
 does not change the order
 return number of actually particles left


Int_t RemoveLowPt(TObjArray* parts, Double_t minPt) const



 remove particles in 'parts' that have pt lower than 'minPt'
 does not change the order
 return number of actually particles left


Double_t Associate(const H1Part* ref1, const H1Part* ref2, 
				   H1PartJet*& outCloseTo1, H1PartJet*& outCloseTo2)



 taking outCloseTo1/2 as input, which combination
       (ref1 <-> outCloseTo1 + ref2 <-> outCloseTo2)
 or    (ref1 <-> outCloseTo2 + ref2 <-> outCloseTo1)
 gives the minimal sum of eta-phi distances.
 outCloseTo1/2 are changed that 1 is the one for ref1, analog for2,
 return sum of distances.






Inline Functions




                  void OpenTree()
                  void OpenOds(Bool_t yesNo = kTRUE)
                  void OpenMods(Bool_t yesNo = kTRUE)
                  void OpenHat(Bool_t yesNo = kTRUE)
                  void FillHistsYjbRes(const H1PartDstar* ds)
                Bool_t IsS83(const H1PartDstar* ds = 0) const
                Bool_t IsS84(const H1PartDstar* ds = 0) const
                Bool_t IsS83NoL4Check(const H1PartDstar* dstar) const
                Bool_t IsS84NoL4Check(const H1PartDstar* dstar) const
                Bool_t IsL4Ref(Int_t subtrig) const
                Bool_t IsL4Found(Int_t hqsel45Bit, Int_t hqselMode, const H1BankPtr& yrclYecl) const
              Double_t GetCutPt1()
              Double_t GetCutPt2()
        const TString* GetPeriod() const
                  void SetJetFlag(Int_t flag)
                Bool_t SetIgnoreDedx(Bool_t ignore)
                Bool_t SetIgnoreSumEt(Bool_t ignore)
                Bool_t SetIgnoreSumPtKPi(Bool_t ignore)
          H1EventList* GetEventList()
              Double_t GetWeight(Int_t trigger, const H1PartDstar* ds) const
               Float_t GetWgammaP33() const
               Float_t GetWgammaP44() const
               Float_t GetY33() const
               Float_t GetY44() const
              Double_t GetXgammaMassHelp(const H1Part* dstar, const H1PartJet* jet) const
              Double_t GetXgammaMassDiJetHelp(const H1Part* dstar, const H1PartJet* dsJet, const H1PartJet* othJet) const
              Double_t GetXpMassHelp(const H1Part* dstar, const H1PartJet* jet) const
               TClass* Class()
               TClass* IsA() const
                  void ShowMembers(TMemberInspector& insp, char* parent)
                  void Streamer(TBuffer& b)
                  void StreamerNVirtual(TBuffer& b)







Author: Gero Flucke

Last update: Wed Jan 11 14:06:25 2006









ROOT page - Class index - Top of the page









This page has been automatically generated. If you have any comments or suggestions about the page layout send a mail to ROOT support, or contact the developers with any questions or problems regarding ROOT.