thailand2017/tut_ana_comb.C File Reference

#include "auxtut.C"

Go to the source code of this file.

Functions
void	tut_ana_comb (int nevts=0, TString prefix="signal")

Function Documentation

void tut_ana_comb	(	int	nevts = 0,
		TString	prefix = "signal"
	)

Definition at line 17 of file thailand2017/tut_ana_comb.C.

References fRun, PndAnalysis::GetEntries(), PndAnalysis::GetEvent(), i, initrun(), out, plotmyhistos(), PndAnalysis::PndAnalysis(), RemoveGeoManager(), TString, and writemyhistos().

{
        // *** some variables
        int i=0,j=0, k=0, l=0;
        gStyle->SetOptFit(1011);
        
        // *** Initialize FairRunAna with defaults      
        TString OutFile="out_dummy.root";  
        FairRunAna* fRun = initrun(prefix, OutFile);
        fRun->Init();
        
        // *** create an output file for all histograms
        TFile *out = TFile::Open(prefix+"_ana_comb.root","RECREATE");
        
        // #### EXERCISE:  create some histograms
        //TH1F *hjpsiMass = ...
        
        // ------------------------------
        // Now the analysis stuff comes...
        // ------------------------------
        
        // *** the data reader object
        PndAnalysis* theAnalysis = new PndAnalysis();
        if (nevts==0) nevts= theAnalysis->GetEntries();
        
        // #### EXERCISE: ---> Declare mu+, mu-, pi+, pi-, jpsi and pbarp RhoCandLists 

        // RhoCandList muplus...
        
        
        // *** Mass selector for the jpsi cands
        RhoMassParticleSelector *jpsiMassSel=new RhoMassParticleSelector("jpsi",3.096,1.0);
        
        // ***
        // the event loop
        // ***
        while (theAnalysis->GetEvent() && i++<nevts)
        {
                if ((i%100)==0) cout<<"evt " << i << endl;
                                
                // #### EXERCISE: ---> Fill mu+, mu-, pi+ and pi- lists 

                // theAnalysis->FillList(muplus,  "MuonAllPlus");
                // ....
                
                // #### EXERCISE:  Do combinatorics with the generic lists to form J/psi -> mu+ mu- and pbarp -> J/psi pi+ pi-
                
                //jpsi.Combine(...)
                //pbarp.Combine(...)
                
                // #### EXERCISE:  Fill histograms with masses of J/psi and ppbar candidates
                // fillM(jpsi, hjpsiMass);
                // ...
                
        }
        
        // *** change to directory where histograms are created
        out->cd();

        // *** plot all histos
        plotmyhistos();
        
        // *** write out all the histos to file
        int nhist = writemyhistos();
        cout<<"Writing "<<nhist<<" histograms to file"<<endl;           
        out->Save();    
        
        // *** temporaty fix to avoid error on macro exit
        RemoveGeoManager();
}