- // Books the histograms where the results of the Resolution analysis are stored
-
- if(fhPhotonEnergy)
- delete fhPhotonEnergy ;
- if(fhPhotonAllEnergy)
- delete fhPhotonAllEnergy ;
- if(fhElectronEnergy)
- delete fhElectronEnergy ;
- if(fhElectronAllEnergy)
- delete fhElectronAllEnergy ;
- if(fhNeutralHadronEnergy)
- delete fhNeutralHadronEnergy ;
- if(fhNeutralEMEnergy)
- delete fhNeutralEMEnergy ;
- if(fhNeutralEMAllEnergy)
- delete fhNeutralEMAllEnergy ;
- if(fhChargedHadronEnergy)
- delete fhChargedHadronEnergy ;
- if(fhPhotonHadronEnergy)
- delete fhPhotonHadronEnergy ;
- if(fhPhotonPosition)
- delete fhPhotonPosition ;
- if(fhPhotonAllPosition)
- delete fhPhotonAllPosition ;
- if(fhElectronPosition)
- delete fhElectronPosition ;
- if(fhElectronAllPosition)
- delete fhElectronAllPosition ;
- if(fhNeutralHadronPosition)
- delete fhNeutralHadronPosition ;
- if(fhNeutralEMPosition)
- delete fhNeutralEMPosition ;
- if(fhNeutralEMAllPosition)
- delete fhNeutralEMAllPosition ;
- if(fhChargedHadronPosition)
- delete fhChargedHadronPosition ;
- if(fhPhotonHadronPosition)
- delete fhPhotonHadronPosition ;
-
- fhPhotonEnergy = new TH2F("hPhotonEnergy", "hPhotonEnergy", 100, 0., 5., 100, 0., 5.);
- fhPhotonAllEnergy = new TH2F("hPhotonAllEnergy", "hPhotonAllEnergy", 100, 0., 5., 100, 0., 5.);
- fhElectronEnergy = new TH2F("hElectronEnergy","hElectronEnergy", 100, 0., 5., 100, 0., 5.);
- fhElectronAllEnergy = new TH2F("hElectronAllEnergy","hElectronAllEnergy", 100, 0., 5., 100, 0., 5.);
- fhNeutralHadronEnergy = new TH2F("hNeutralHadronEnergy", "hNeutralHadronEnergy", 100, 0., 5., 100, 0., 5.);
- fhNeutralEMEnergy = new TH2F("hNeutralEMEnergy", "hNeutralEMEnergy", 100, 0., 5., 100, 0., 5.);
- fhNeutralEMAllEnergy = new TH2F("hNeutralEMAllEnergy", "hNeutralEMAllEnergy", 100, 0., 5., 100, 0., 5.);
- fhChargedHadronEnergy = new TH2F("hChargedHadronEnergy", "hChargedHadronEnergy", 100, 0., 5., 100, 0., 5.);
- fhPhotonHadronEnergy = new TH2F("hPhotonHadronEnergy","hPhotonHadronEnergy", 100, 0., 5., 100, 0., 5.);
- fhPhotonPosition = new TH2F("hPhotonPosition","hPhotonPosition", 20, 0., 5., 100, 0., 5.);
- fhPhotonAllPosition = new TH2F("hPhotonAllPosition","hPhotonAllPosition", 20, 0., 5., 100, 0., 5.);
- fhElectronPosition = new TH2F("hElectronPosition","hElectronPosition", 20, 0., 5., 100, 0., 5.);
- fhElectronAllPosition = new TH2F("hElectronAllPosition","hElectronAllPosition", 20, 0., 5., 100, 0., 5.);
- fhNeutralHadronPosition = new TH2F("hNeutralHadronPosition","hNeutralHadronPosition", 20, 0., 5., 100, 0., 5.);
- fhNeutralEMPosition = new TH2F("hNeutralEMPosition","hNeutralEMPosition", 20, 0., 5., 100, 0., 5.);
- fhNeutralEMAllPosition = new TH2F("hNeutralEMAllPosition","hNeutralEMAllPosition", 20, 0., 5., 100, 0., 5.);
- fhChargedHadronPosition = new TH2F("hChargedHadronPosition","hChargedHadronPosition", 20, 0., 5., 100, 0., 5.);
- fhPhotonHadronPosition = new TH2F("hPhotonHadronPosition","hPhotonHadronPosition", 20, 0., 5., 100, 0., 5.);
- fhPhotonReg = new TH1F("hPhotonReg","hPhotonReg", 20, 0., 5.);
- fhAllReg = new TH1F("hAllReg", "hAllReg", 20, 0., 5.);
- fhNReg = new TH1F("hNReg", "hNReg", 20, 0., 5.);
- fhNBarReg = new TH1F("hNBarReg", "hNBarReg", 20, 0., 5.);
- fhChargedReg= new TH1F("hChargedReg", "hChargedReg", 20, 0., 5.);
-
- fhPhotonEM = new TH1F("hPhotonEM","hPhotonEM", 20, 0., 5.);
- fhAllEM = new TH1F("hAllEM", "hAllEM", 20, 0., 5.);
- fhNEM = new TH1F("hNEM", "hNEM", 20, 0., 5.);
- fhNBarEM = new TH1F("hNBarEM", "hNBarEM", 20, 0., 5.);
- fhChargedEM= new TH1F("hChargedEM", "hChargedEM", 20, 0., 5.);
-
- fhPrimary= new TH1F("hPrimary", "hPrimary", 20, 0., 5.);
-
- fhPhotPhot = new TH1F("hPhotPhot","hPhotPhot", 20, 0., 5.); //Photon registered as photon
- fhPhotElec = new TH1F("hPhotElec","hPhotElec", 20, 0., 5.); //Photon registered as Electron
- fhPhotNeuH = new TH1F("hPhotNeuH","hPhotNeuH", 20, 0., 5.); //Photon registered as Neutral Hadron
- fhPhotNuEM = new TH1F("hPhotNuEM","hPhotNuEM", 20, 0., 5.); //Photon registered as Neutral EM
- fhPhotChHa = new TH1F("hPhotChHa","hPhotChHa", 20, 0., 5.); //Photon registered as Charged Hadron
- fhPhotGaHa = new TH1F("hPhotGaHa","hPhotGaHa", 20, 0., 5.); //Photon registered as Gamma-Hadron
-
+ // Calculates Real and Mixed invariant mass distributions
+ if (fRunLoader == 0x0)
+ {
+ AliError(Form("Error Loading session"));
+ return;
+ }
+
+ AliPHOSLoader* gime = dynamic_cast<AliPHOSLoader*>(fRunLoader->GetLoader("PHOSLoader"));
+ if ( gime == 0 )
+ {
+ AliError(Form("Could not obtain the Loader object !"));
+ return ;
+ }
+
+ gime->LoadRecParticles("READ");
+
+ Int_t nMixedEvents = 4 ; //# of events used for calculation of 'mixed' distribution
+
+
+ //opening file
+ TFile * mfile = new TFile("invmass.root","update");
+
+ //========== Reading /Booking Histograms
+ TH2D * hRealEM = 0 ;
+ hRealEM = (TH2D*) mfile->Get("hRealEM") ;
+ if(hRealEM == 0)
+ hRealEM = new TH2D("hRealEM", "Real for EM particles", 250,0.,1.,40,0.,4.) ;
+ TH2D * hRealPhot = 0 ;
+
+ hRealPhot = (TH2D*)mfile->Get("hRealPhot");
+ if(hRealPhot == 0)
+ hRealPhot = new TH2D("hRealPhot", "Real for kPhoton particles", 250,0.,1.,40,0.,4.) ;
+
+ TH2D * hMixedEM = 0 ;
+ hMixedEM = (TH2D*) mfile->Get("hMixedEM") ;
+ if(hMixedEM == 0)
+ hMixedEM = new TH2D("hMixedEM", "Mixed for EM particles", 250,0.,1.,40,0.,4.) ;
+
+ TH2D * hMixedPhot = 0 ;
+ hMixedPhot = (TH2D*) mfile->Get("hMixedPhot") ;
+ if(hMixedPhot == 0)
+ hMixedPhot = new TH2D("hMixedPhot","Mixed for kPhoton particles",250,0.,1.,40,0.,4.) ;
+