// coherence check
if (!AliRsnCut::TargetOK(obj1))
{
- AliError(Form("Wrong target. Skipping cut", GetName()));
+ AliError("Wrong target. Skipping cut");
return kTRUE;
}
// coherence check
if (!AliRsnCut::TargetOK(obj1))
{
- AliError(Form("Wrong target. Skipping cut", GetName()));
+ AliError("Wrong target. Skipping cut");
return kTRUE;
}
if (fRef)
{
- AliInfo(Form(".......Px, Py, Pz, Pt (ref) : %d %d %d %d", fP.X(), fP.Y(), fP.Z(), fP.Perp()));
+ AliInfo(Form(".......Px, Py, Pz, Pt (ref) : %f %f %f %f", fP.X(), fP.Y(), fP.Z(), fP.Perp()));
} else AliInfo("....... absent REF");
if (fRefMC)
{
- AliInfo(Form(".......Px, Py, Pz, Pt (ref MC): %d %d %d %d", fP.X(), fP.Y(), fP.Z(), fP.Perp()));
+ AliInfo(Form(".......Px, Py, Pz, Pt (ref MC): %f %f %f %f", fP.X(), fP.Y(), fP.Z(), fP.Perp()));
AliInfo(Form(".......PDG code : %d", fRefMC->Particle()->GetPdgCode()));
AliInfo(Form(".......Mother (label) : %d", fRefMC->Particle()->GetFirstMother()));
AliInfo(Form(".......Mother (PDG code) : %d", fMotherPDG));
// update the various axes using the definitions given in the array of axes here
AliRsnValue *fcnAxis = 0;
- TAxis *axis = 0;
for (Int_t i = 0; i < fSize; i++)
{
fcnAxis = (AliRsnValue*)fAxisList.At(i);
- axis = fHSparse->GetAxis(i);
if (!fcnAxis) {
- axis->Set(1, -1E5, 1E5);
AliError("Empty axis: doing unique bin betweeen -100000 and 100000");
continue;
}
- axis->Set(fcnAxis->GetArray().GetSize() - 1, fcnAxis->GetArray().GetArray());
+ fHSparse->SetBinEdges(i, fcnAxis->GetArray().GetArray());
}
return fHSparse;
}
if (charge != '+' && charge != '-' && charge != '0')
{
- AliError(Form("Character '%c' not recognized as charge sign"));
+ AliError(Form("Character '%c' not recognized as charge sign", charge));
return kFALSE;
}
if (type < 0 && type > (Int_t)AliPID::kSPECIESN)
// to inform about number of events processed
//
- if ((num+1) % fPrintInfoNumber == 0) AliInfo(Form("Events processed %d",num+1));
+ if ((num+1) % fPrintInfoNumber == 0) AliInfo(Form("Events processed %d", (Int_t)num+1));
}
}
AliInfo(Form("=== %s ==================", GetName()));
- AliInfo(Form("Number Of Events Processed : %10d", (Long64_t)hEventInfo->Integral()));
- AliInfo(Form("Number Of Events Accepted : %10d", (Long64_t)hEventInfo->GetBinContent(2)));
- AliInfo(Form("Number Of Events Skipped : %10d", (Long64_t)hEventInfo->GetBinContent(1)));
+ AliInfo(Form("Number Of Events Processed : %10d", (Int_t)hEventInfo->Integral()));
+ AliInfo(Form("Number Of Events Accepted : %10d", (Int_t)hEventInfo->GetBinContent(2)));
+ AliInfo(Form("Number Of Events Skipped : %10d", (Int_t)hEventInfo->GetBinContent(1)));
AliInfo(Form("=== end %s ==============", GetName()));
AliDebug(AliLog::kDebug+2,"->");
}
AliInfo(Form("=== %s ==================",GetName()));
- AliInfo(Form("Number Of Events Processed : %10d",(Long64_t)hEventInfo->Integral()));
- AliInfo(Form("Number Of Events Accepted : %10d",(Long64_t)hEventInfo->GetBinContent(2)));
- AliInfo(Form("Number Of Events Skipped : %10d",(Long64_t)hEventInfo->GetBinContent(1)));
+ AliInfo(Form("Number Of Events Processed : %10d",(Int_t)hEventInfo->Integral()));
+ AliInfo(Form("Number Of Events Accepted : %10d",(Int_t)hEventInfo->GetBinContent(2)));
+ AliInfo(Form("Number Of Events Skipped : %10d",(Int_t)hEventInfo->GetBinContent(1)));
AliInfo(Form("=== end %s ==============",GetName()));
AliDebug(AliLog::kDebug+2, "->");
else fValue = 0;
}
break;
+ case kQInv:
+ {
+ TLorentzVector diff = mother->GetDaughter(0)->P() - mother->GetDaughter(1)->P();
+ fValue = diff.M();
+ }
+ break;
default:
AliWarning("Invalid value type");
return kFALSE;
kPairCosThetaStarMC2,
kAngleToLeading,
kLeadingPt,
+ kQInv,
kEventMult,
kValueTypes
};
// function axes
Double_t ybins[] = {-0.8, -0.7, -0.6, -0.5, 0.5, 0.6, 0.7, 0.8};
- AliRsnValue *axisIM = new AliRsnValue("IM", AliRsnValue::kPairInvMass, 50, 0.9, 1.4);
- AliRsnValue *axisPt = new AliRsnValue("PT", AliRsnValue::kPairPt, 0.0, 20.0, 0.1);
- AliRsnValue *axisY = new AliRsnValue("Y" , AliRsnValue::kPairY, sizeof(ybins)/sizeof(ybins[0]), ybins);
- AliRsnValue *axisEta = new AliRsnValue("Eta", AliRsnValue::kPairEta, sizeof(ybins)/sizeof(ybins[0]), ybins);
+ AliRsnValue *axisIM = new AliRsnValue("IM" , AliRsnValue::kPairInvMass, 50, 0.9, 1.4);
+ AliRsnValue *axisPt = new AliRsnValue("PT" , AliRsnValue::kPairPt, 0.0, 20.0, 0.1);
+ AliRsnValue *axisY = new AliRsnValue("Y" , AliRsnValue::kPairY, sizeof(ybins)/sizeof(ybins[0]), ybins);
+ AliRsnValue *axisQinv = new AliRsnValue("QInv", AliRsnValue::kQInv, 100, 0.0, 10.0);
// functions for TH1-like output
AliRsnFunction *fcnPt = new AliRsnFunction;
fcnPt ->AddAxis(axisIM);
fcnPt ->AddAxis(axisPt);
fcnPt ->AddAxis(axisY);
- fcnPt ->AddAxis(axisEta);
+ fcnPt ->AddAxis(axisQinv);
// add functions to TH1-like output
pairPMhist->AddFunction(fcnPt);
// add values to TNtuple-like output
pairPMntp->AddValue(axisIM);
pairPMntp->AddValue(axisPt);
+ pairPMntp->AddValue(axisY);
+ pairPMntp->AddValue(axisQinv);
// add everything to analysis manager
- task->GetAnalysisManager()->Add(pairPMhist);
- //task->GetAnalysisManager()->Add(pairPMntp);
+ //task->GetAnalysisManager()->Add(pairPMhist);
+ task->GetAnalysisManager()->Add(pairPMntp);
return kTRUE;
}
//
void runLocal
(
- Int_t nReadFiles = 1,
+ Int_t nReadFiles = 0,
Int_t nSkipFiles = 0,
const char *addTaskName = "AddAnalysisTaskRsnTest.C",
- const char *inputSource = "/home/pulvir/analysis/resonances/LHC2010-7TeV-phi/rsn-package/sim.txt",
+ const char *inputSource = "/home/pulvir/analysis/resonances/LHC2010-7TeV-phi/alien+plugin/sim.txt",
const char *dataLabel = "7TeV_pass2_sim_ESD",
const char *outName = "rsn-test"
)