Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
#ifndef _OmegaPiHist_H
#define _OmegaPiHist_H
#include <iostream>
#include <fstream>
#include <string>
#include <vector>
#include <cassert>
#include <boost/shared_ptr.hpp>
#include "TROOT.h"
// #include <TSystem.h>
#include "qft++/topincludes/relativistic-quantum-mechanics.hh"
#include "Examples/MATpbarpToOmegaPi/OmegaPiData.hh"
using OmegaPiData::OmPiEvtData;
class OmegaPiEventList;
class OmegaPiLh;
class TFile;
class TH2F;
class TH1F;
class OmegaPiHist {
public:
// create/copy/destroy:
///Constructor
OmegaPiHist(boost::shared_ptr<const OmegaPiEventList>);
OmegaPiHist(boost::shared_ptr<OmegaPiLh>, OmegaPiData::fitParamVal&);
/** Destructor */
virtual ~OmegaPiHist();
// Getters:
protected:
private:
TFile* _theTFile;
TH1F* _cosOmegaHeliDataHist;
TH1F* _cosOmegaHeliMcHist;
TH1F* _cosOmegaHeliFittedHist;
TH1F* _cosOmegaAccCorHist;
TH1F* _cosPi0FromOmegaDataHeli;
TH1F* _cosPi0FromOmegaMcHeli;
TH1F* _cosPi0FromOmegaFittedHeli;
TH1F* _cosPi0FromOmegaAccCorHeli;
TH1F* _treimanYangDataHist;
TH1F* _treimanYangMcHist;
TH1F* _treimanYangFittedHist;
unsigned _jmax;
unsigned _pbarmom;
void initRootStuff();
void plotCosOmegaHeli(TH1F* theHisto, const OmPiEvtData* theEvtData, double weight);
void plotCosPi0FromOmegaHeli(TH1F* theHisto, const OmPiEvtData* theEvtData, double weight);
void plotTreimanYang(TH1F* theHisto, const OmPiEvtData* theEvtData, double weight);