AliAODForwardMult* forward = static_cast<AliAODForwardMult*>(obj);
Float_t cent = forward->GetCentrality();
- //if(cent > 0) {
- // if( cent < 40 || cent >50 ) return;
- // else std::cout<<"selecting event with cent "<<cent<<std::endl;
- //}
-
+ if(cent > 0) {
+ if( cent < fCentLow || cent > fCentHigh ) return;
+ else std::cout<<"selecting event with cent "<<cent<<std::endl;
+ }
+ //std::cout<<fCentLow<<" "<<fCentHigh<<std::endl;
// Get the histogram(s)
TH2D* data = GetHistogram(aod, false);
TH2D* dataMC = GetHistogram(aod, true);
const char* name = GetName();
// Get acceptance normalisation from underflow bins
- TH1D* norm = ProjectX(fSum, Form("norm%s",name), 0, 1, fCorrEmpty, false);
+ TH1D* norm = ProjectX(fSum, Form("norm%s",name), 0, 0, fCorrEmpty, false);
if(fShapeCorr && (fTriggerEff != 1)) {AliInfo("Objects already set for normalization - no action taken"); return; }
- //std::cout<<Form("normalizationHists_%s_%s.root",type.Data(),snn.Data())<<std::endl;
TFile* fin = TFile::Open(Form("$ALICE_ROOT/PWG2/FORWARD/corrections/Normalization/normalizationHists_%s_%s.root",type.Data(),snn.Data()));
if(!fin) AliWarning("no file for normalization");
for(Int_t j = 1; j<=fRebin;j++) {
Int_t bin = (i-1)*fRebin + j;
Double_t c = h->GetBinContent(bin);
-
if (c <= 0) continue;
-
+
if (fCutEdges) {
if (h->GetBinContent(bin+1)<=0 ||
h->GetBinContent(bin-1)<=0) {
*
* @param h Correction
*/
+ void SetShapeCorrection(const TH1* h);
+ /**
+ * Load the normalization data - done automatically if not set from outside
+ *
+ * @param sys system
+ * @param energy energy
+ */
void LoadNormalizationData(UShort_t sys, UShort_t energy);
+ /**
+ * Load the normalization data - done automatically if not set from outside
+ *
+ * @param centlow low cent cut
+ */
+ void SetCentLow(Float_t centlow) { fCentLow = centlow; }
+ /**
+ * Load the normalization data - done automatically if not set from outside
+ *
+ * @param centhigh high cent cut
+ */
+ void SetCentHigh(Float_t centhigh) { fCentHigh = centhigh; }
- void SetShapeCorrection(const TH1* h);
/**
* Destructor
*