- 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
+// if(fhAllEnergy)
+// delete fhAllEnergy ;
+// if(fhPhotEnergy)
+// delete fhPhotEnergy ;
+// if(fhEMEnergy)
+// delete fhEMEnergy ;
+// if(fhPPSDEnergy)
+// delete fhPPSDEnergy ;
+
+
+ fhAllEnergy = new TH2F("hAllEnergy", "Energy of any RP with primary photon",100, 0., 5., 100, 0., 5.);
+ fhPhotEnergy = new TH2F("hPhotEnergy", "Energy of kGAMMA with primary photon",100, 0., 5., 100, 0., 5.);
+ fhEMEnergy = new TH2F("hEMEnergy", "Energy of EM with primary photon", 100, 0., 5., 100, 0., 5.);
+ fhPPSDEnergy = new TH2F("hPPSDEnergy", "Energy of PPSD with primary photon", 100, 0., 5., 100, 0., 5.);
+
+// if(fhAllPosition)
+// delete fhAllPosition ;
+// if(fhPhotPosition)
+// delete fhPhotPosition ;
+// if(fhEMPosition)
+// delete fhEMPosition ;
+// if(fhPPSDPosition)
+// delete fhPPSDPosition ;
+
+
+ fhAllPosition = new TH2F("hAllPosition", "Position of any RP with primary photon",100, 0., 5., 100, 0., 5.);
+ fhPhotPosition = new TH2F("hPhotPosition", "Position of kGAMMA with primary photon",100, 0., 5., 100, 0., 5.);
+ fhEMPosition = new TH2F("hEMPosition", "Position of EM with primary photon", 100, 0., 5., 100, 0., 5.);
+ fhPPSDPosition = new TH2F("hPPSDPosition", "Position of PPSD with primary photon", 100, 0., 5., 100, 0., 5.);
+
+// if(fhAllReg)
+// delete fhAllReg ;
+// if(fhPhotReg)
+// delete fhPhotReg ;
+// if(fhNReg)
+// delete fhNReg ;
+// if(fhNBarReg)
+// delete fhNBarReg ;
+// if(fhChargedReg)
+// delete fhChargedReg ;
+
+ fhAllReg = new TH1F("hAllReg", "All primaries registered as photon", 100, 0., 5.);
+ fhPhotReg = new TH1F("hPhotReg", "Photon registered as photon", 100, 0., 5.);
+ fhNReg = new TH1F("hNReg", "N registered as photon", 100, 0., 5.);
+ fhNBarReg = new TH1F("hNBarReg", "NBar registered as photon", 100, 0., 5.);
+ fhChargedReg= new TH1F("hChargedReg", "Charged hadron registered as photon",100, 0., 5.);
+
+// if(fhAllEM)
+// delete fhAllEM ;
+// if(fhPhotEM)
+// delete fhPhotEM ;
+// if(fhNEM)
+// delete fhNEM ;
+// if(fhNBarEM)
+// delete fhNBarEM ;
+// if(fhChargedEM)
+// delete fhChargedEM ;
+
+ fhAllEM = new TH1F("hAllEM", "All primary registered as EM",100, 0., 5.);
+ fhPhotEM = new TH1F("hPhotEM", "Photon registered as EM", 100, 0., 5.);
+ fhNEM = new TH1F("hNEM", "N registered as EM", 100, 0., 5.);
+ fhNBarEM = new TH1F("hNBarEM", "NBar registered as EM", 100, 0., 5.);
+ fhChargedEM= new TH1F("hChargedEM","Charged registered as EM",100, 0., 5.);
+
+// if(fhAllPPSD)
+// delete fhAllPPSD ;
+// if(fhPhotPPSD)
+// delete fhPhotPPSD ;
+// if(fhNPPSD)
+// delete fhNPPSD ;
+// if(fhNBarPPSD)
+// delete fhNBarPPSD ;
+// if(fhChargedPPSD)
+// delete fhChargedPPSD ;
+
+ fhAllPPSD = new TH1F("hAllPPSD", "All primary registered as PPSD",100, 0., 5.);
+ fhPhotPPSD = new TH1F("hPhotPPSD", "Photon registered as PPSD", 100, 0., 5.);
+ fhNPPSD = new TH1F("hNPPSD", "N registered as PPSD", 100, 0., 5.);
+ fhNBarPPSD = new TH1F("hNBarPPSD", "NBar registered as PPSD", 100, 0., 5.);
+ fhChargedPPSD= new TH1F("hChargedPPSD","Charged registered as PPSD",100, 0., 5.);
+
+// if(fhPrimary)
+// delete fhPrimary ;
+ fhPrimary= new TH1F("hPrimary", "hPrimary", 100, 0., 5.);
+
+// if(fhAllRP)
+// delete fhAllRP ;
+// if(fhVeto)
+// delete fhVeto ;
+// if(fhShape)
+// delete fhShape ;
+// if(fhPPSD)
+// delete fhPPSD ;
+
+ fhAllRP = new TH1F("hAllRP","All Reconstructed particles", 100, 0., 5.);
+ fhVeto = new TH1F("hVeto", "All uncharged particles", 100, 0., 5.);
+ fhShape = new TH1F("hShape","All particles with EM shaower",100, 0., 5.);
+ fhPPSD = new TH1F("hPPSD", "All PPSD photon particles", 100, 0., 5.);
+
+
+// if(fhPhotPhot)
+// delete fhPhotPhot ;
+// if(fhPhotElec)
+// delete fhPhotElec ;
+// if(fhPhotNeuH)
+// delete fhPhotNeuH ;
+// if(fhPhotNuEM)
+// delete fhPhotNuEM ;
+// if(fhPhotChHa)
+// delete fhPhotChHa ;
+// if(fhPhotGaHa)
+// delete fhPhotGaHa ;
+
+ fhPhotPhot = new TH1F("hPhotPhot","hPhotPhot", 100, 0., 5.); //Photon registered as photon
+ fhPhotElec = new TH1F("hPhotElec","hPhotElec", 100, 0., 5.); //Photon registered as Electron
+ fhPhotNeuH = new TH1F("hPhotNeuH","hPhotNeuH", 100, 0., 5.); //Photon registered as Neutral Hadron
+ fhPhotNuEM = new TH1F("hPhotNuEM","hPhotNuEM", 100, 0., 5.); //Photon registered as Neutral EM
+ fhPhotChHa = new TH1F("hPhotChHa","hPhotChHa", 100, 0., 5.); //Photon registered as Charged Hadron
+ fhPhotGaHa = new TH1F("hPhotGaHa","hPhotGaHa", 100, 0., 5.); //Photon registered as Gamma-Hadron