Changeset 7252 for trunk/MagicSoft/Mars/datacenter/macros
- Timestamp:
- 08/02/05 13:41:41 (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/datacenter/macros/fillstar.C
r7228 r7252 105 105 } 106 106 107 TProfile *h1 = (TProfile*)arr.FindObjectInCanvas("RingBroadening", "TProfile", "MHMuonPar"); 108 if (!h1) 109 { 110 cout << "WARNING - Reading of RingBroadening failed." << endl; 111 return 0; 112 } 113 114 Float_t numerator = 0; 115 Float_t denominator = 0; 116 Float_t summe = 0; 117 118 for(Int_t i=5; i<15; i++) 119 { 120 numerator += h1->GetBinContent(i)*h1->GetBinEntries(i); 121 denominator += h1->GetBinEntries(i); 122 } 123 124 summe = numerator/denominator; 125 126 Float_t psf = (summe - 0.04816)/0.001294; 107 MHMuonPar *hmuon = (TProfile*)arr.FindObjectInCanvas("MHMuonPar", "MHMuonPar", "MHMuonPar"); 108 if (!hmuon) 109 { 110 cout << "WARNING - Reading of MHMuon failed." << endl; 111 return 0; 112 } 113 114 Float_t psf = (muon->GetMeanWidth() - 0.04816)/0.001294; 127 115 psf = TMath::Nint(psf*10)/10.; 128 116 TString PSF = Form("%5.1f", psf); 129 Int_t num = (int)h1->GetEntries(); 130 131 TProfile *h2 = (TProfile*)arr.FindObjectInCanvas("SizeVsRadius", "TProfile", "MHMuonPar"); 132 if (!h1) 133 { 134 cout << "WARNING - Reading of SizeVsRadius failed." << endl; 135 return 0; 136 } 137 138 numerator = 0; 139 denominator = 0; 140 summe = 0; 141 142 for(Int_t i=5; i<15; i++) 143 { 144 numerator += h2->GetBinContent(i)*h2->GetBinEntries(i); 145 denominator += h2->GetBinEntries(i); 146 } 147 148 summe = numerator/denominator; 117 Int_t num = (int)muon.GetEntries(); 149 118 150 119 Float_t integralmc = -1.43*psf + 1035; 151 Float_t ratiodatamc = ( summe/integralmc)*100;120 Float_t ratiodatamc = (muon->GetMeanSize()/integralmc)*100; 152 121 TString ratio = Form("%5.1f", ratiodatamc); 153 154 122 155 123 TH1 *h = (TH1*)arr.FindObjectInCanvas("Islands", "TH1F", "MHImagePar");
Note:
See TracChangeset
for help on using the changeset viewer.