//________________________________________________________________________
void AliAnalysisTaskAODCentralityMaker::UserExec(Option_t */*option*/)
{
+// User Exec
AliVEvent* event = InputEvent();
AliESDEvent* esd = dynamic_cast<AliESDEvent*>(event);
if (!esd) {
private:
- AliAODCentrality *fAODCentrality;
+
AliAnalysisTaskAODCentralityMaker(const AliAnalysisTaskAODCentralityMaker &source);
AliAnalysisTaskAODCentralityMaker& operator=(const AliAnalysisTaskAODCentralityMaker& source);
-
+ AliAODCentrality *fAODCentrality; // AOD centrality pointer
TString fDeltaAODFileName; // Name of output file
AliAODHeader* fAODHeader; // Header for replaction
\r
//______________________________________________________________________________\r
void AliAnalysisTaskESDfilter::PrintMCInfo(AliStack *pStack,Int_t label){\r
+// Print MC info\r
if(!pStack)return;\r
label = TMath::Abs(label);\r
TParticle *part = pStack->Particle(label);\r
fFile2->Close();
}
//________________________________________________________________________
-Int_t AliCentralitySelectionTask::SetupRun(AliESDEvent* esd)
+Int_t AliCentralitySelectionTask::SetupRun(AliESDEvent* const esd)
{
// Setup files for run
}
//________________________________________________________________________
-Bool_t AliCentralitySelectionTask::IsOutlierV0MSPD(Float_t spd, Float_t v0, Int_t cent)
+Bool_t AliCentralitySelectionTask::IsOutlierV0MSPD(Float_t spd, Float_t v0, Int_t cent) const
{
+// Clean outliers
Float_t val= -0.143789 + 0.288874 * v0;
Float_t SPDsigma[100]={231.483, 189.446, 183.359, 179.923, 174.229, 170.309, 165.021,
160.84, 159.33, 154.453, 151.644, 148.337, 145.215, 142.353,
}
//________________________________________________________________________
-Bool_t AliCentralitySelectionTask::IsOutlierV0MTPC(Int_t tracks, Float_t v0, Int_t cent)
+Bool_t AliCentralitySelectionTask::IsOutlierV0MTPC(Int_t tracks, Float_t v0, Int_t cent) const
{
+// Clean outliers
Float_t val = -0.540691 + 0.128358 * v0;
Float_t TPCsigma[100]={106.439, 89.2834, 86.7568, 85.3641, 83.379, 81.6093, 79.3189,
78.0616, 77.2167, 75.0021, 73.9957, 72.0926, 71.0442, 69.8395,
}
//________________________________________________________________________
-Bool_t AliCentralitySelectionTask::IsOutlierV0MZDC(Float_t zdc, Float_t v0)
+Bool_t AliCentralitySelectionTask::IsOutlierV0MZDC(Float_t zdc, Float_t v0) const
{
+// Clean outliers
Float_t val1 = 6350. - 0.26 * v0;
Float_t val2 = 5580.;
if ((zdc > val1) || (zdc > val2))
}
//________________________________________________________________________
-Bool_t AliCentralitySelectionTask::IsOutlierV0MZDCECal(Float_t /*zdc*/, Float_t /*v0*/)
+Bool_t AliCentralitySelectionTask::IsOutlierV0MZDCECal(Float_t /*zdc*/, Float_t /*v0*/) const
{
+// Clean outliers
return kFALSE;
}
//________________________________________________________________________
-Float_t AliCentralitySelectionTask::MyGetScaleFactor(Int_t runnumber, Int_t flag)
+Float_t AliCentralitySelectionTask::MyGetScaleFactor(Int_t runnumber, Int_t flag) const
{
+// Get scaling factor
if (! (runnumber >= fLowRunN && runnumber <=fHighRunN)) {
cout << "MyGetScaleFactor error in run number range " << runnumber << endl;
return 0.0;
}
//________________________________________________________________________
-Float_t AliCentralitySelectionTask::MyGetScaleFactorMC(Int_t runnumber)
+Float_t AliCentralitySelectionTask::MyGetScaleFactorMC(Int_t runnumber) const
{
+// Get MC scaling factor
if (! (runnumber >= fLowRunN && runnumber <=fHighRunN)) {
cout << "MyGetScaleFactor error in run number range " << runnumber << endl;
return 0.0;
//________________________________________________________________________
void AliCentralitySelectionTask::MyInitScaleFactor ()
{
+// Initialize the scaling factors
for (int i=0; i<(fHighRunN-fLowRunN); i++) V0MScaleFactor[i] = 0.0;
for (int i=0; i<(fHighRunN-fLowRunN); i++) SPDScaleFactor[i] = 0.0;
for (int i=0; i<(fHighRunN-fLowRunN); i++) TPCScaleFactor[i] = 0.0;
//________________________________________________________________________
void AliCentralitySelectionTask::MyInitScaleFactorMC()
{
+// Initialize the MC scaling factors
for (int i=0; i<(fHighRunN-fLowRunN); i++) V0MScaleFactorMC[i] = 0.0;
// scale factors determined from <V0 charge> on a run-by-run basis
V0MScaleFactorMC[0] = 0.75108;
void ReadCentralityHistos2(TString filename);
private:
- Int_t SetupRun(AliESDEvent* esd);
- Bool_t IsOutlierV0MSPD(Float_t spd, Float_t v0, Int_t cent);
- Bool_t IsOutlierV0MTPC(Int_t tracks, Float_t v0, Int_t cent);
- Bool_t IsOutlierV0MZDC(Float_t zdc, Float_t v0);
- Bool_t IsOutlierV0MZDCECal(Float_t zdc, Float_t v0);
- Float_t MyGetScaleFactor(Int_t runnumber, Int_t flag);
+ Int_t SetupRun(AliESDEvent* const esd);
+ Bool_t IsOutlierV0MSPD(Float_t spd, Float_t v0, Int_t cent) const;
+ Bool_t IsOutlierV0MTPC(Int_t tracks, Float_t v0, Int_t cent) const;
+ Bool_t IsOutlierV0MZDC(Float_t zdc, Float_t v0) const;
+ Bool_t IsOutlierV0MZDCECal(Float_t zdc, Float_t v0) const;
+ Float_t MyGetScaleFactor(Int_t runnumber, Int_t flag) const;
void MyInitScaleFactor();
- Float_t MyGetScaleFactorMC(Int_t runnumber);
+ Float_t MyGetScaleFactorMC(Int_t runnumber) const;
void MyInitScaleFactorMC();
Int_t fDebug; // Debug flag
}
//____________________________________________________________________
-Float_t AliESDtrackCuts::GetSigmaToVertex(AliESDtrack* esdTrack)
+Float_t AliESDtrackCuts::GetSigmaToVertex(AliESDtrack* const esdTrack)
{
// Calculates the number of sigma to the vertex.
}
//____________________________________________________________________
-Int_t AliESDtrackCuts::CountAcceptedTracks(AliESDEvent* esd)
+Int_t AliESDtrackCuts::CountAcceptedTracks(AliESDEvent* const esd)
{
//
// returns an the number of tracks that pass the cuts
Bool_t AcceptTrack(AliESDtrack* esdTrack);
TObjArray* GetAcceptedTracks(AliESDEvent* esd, Bool_t bTPC = kFALSE);
- Int_t CountAcceptedTracks(AliESDEvent* esd);
+ Int_t CountAcceptedTracks(AliESDEvent* const esd);
static Int_t GetReferenceMultiplicity(AliESDEvent* esd, Bool_t tpcOnly);
Bool_t GetAcceptKinkDaughters() const { return fCutAcceptKinkDaughters;}
Bool_t GetAcceptSharedTPCClusters() const {return fCutAcceptSharedTPCClusters;}
Float_t GetMaxFractionSharedTPCClusters() const {return fCutMaxFractionSharedTPCClusters;}
- void GetMaxCovDiagonalElements(Float_t& c1, Float_t& c2, Float_t& c3, Float_t& c4, Float_t& c5)
+ void GetMaxCovDiagonalElements(Float_t& c1, Float_t& c2, Float_t& c3, Float_t& c4, Float_t& c5) const
{c1 = fCutMaxC11; c2 = fCutMaxC22; c3 = fCutMaxC33; c4 = fCutMaxC44; c5 = fCutMaxC55;}
Float_t GetMaxRel1PtUncertainty() const { return fCutMaxRel1PtUncertainty;}
Float_t GetMaxNsigmaToVertex() const { return fCutNsigmaToVertex;}
void SaveHistograms(const Char_t* dir = 0);
void DrawHistograms();
- static Float_t GetSigmaToVertex(AliESDtrack* esdTrack);
+ static Float_t GetSigmaToVertex(AliESDtrack* const esdTrack);
static void EnableNeededBranches(TTree* tree);
//____________________________________________________________________
Bool_t
AliESDv0Cuts::IsSelected(TList* listObj) {
+// Selection cuts
if(listObj->GetSize()!=4) return kFALSE;
AliESDv0 *esdV0 = (AliESDv0*)listObj->At(0);
AliESDtrack *trackPos = (AliESDtrack*)listObj->At(1);
//____________________________________________________________________
Bool_t
-AliESDv0Cuts::AcceptV0(AliESDv0* esdV0, AliESDtrack* trackPos, AliESDtrack* trackNeg, const AliESDVertex* esdVertex) {
+AliESDv0Cuts::AcceptV0(AliESDv0* const esdV0, AliESDtrack* const trackPos, AliESDtrack* const trackNeg, const AliESDVertex* esdVertex) {
//
// figure out if the v0s survives all the v0 cuts defined
//
}
//____________________________________________________________________
-TObjArray* AliESDv0Cuts::GetAcceptedV0s(AliESD* esd)
+TObjArray* AliESDv0Cuts::GetAcceptedV0s(const AliESD* esd)
{
//
// returns an array of all v0s that pass the cuts
}
//____________________________________________________________________
-Int_t AliESDv0Cuts::CountAcceptedV0s(AliESD* esd)
+Int_t AliESDv0Cuts::CountAcceptedV0s(const AliESD* esd)
{
//
// returns an the number of v0s that pass the cuts
}
//____________________________________________________________________
-TObjArray* AliESDv0Cuts::GetAcceptedV0s(AliESDEvent* esd)
+TObjArray* AliESDv0Cuts::GetAcceptedV0s(const AliESDEvent* esd)
{
//
// returns an array of all v0s that pass the cuts
}
//____________________________________________________________________
-Int_t AliESDv0Cuts::CountAcceptedV0s(AliESDEvent* esd)
+Int_t AliESDv0Cuts::CountAcceptedV0s(const AliESDEvent* esd)
{
//
// returns an the number of v0s that pass the cuts
Bool_t IsSelected(TObject* /*obj*/) {return kTRUE;}
Bool_t IsSelected(TList* listObj);
- Bool_t IsSelected(TObject* obj1, TObject* obj2, TObject* obj3, TObject* obj4)
+ Bool_t IsSelected(TObject* const obj1, TObject* const obj2, TObject* const obj3, TObject* const obj4)
{return AcceptV0((AliESDv0*) obj1, (AliESDtrack*) obj2, (AliESDtrack*) obj3, (const AliESDVertex*) obj4);}
- Bool_t AcceptV0(AliESDv0* esdV0, AliESDtrack* trackPos, AliESDtrack* trackNeg, const AliESDVertex* esdVertex);
- TObjArray* GetAcceptedV0s(AliESD* esd);
- Int_t CountAcceptedV0s(AliESD* esd);
- TObjArray* GetAcceptedV0s(AliESDEvent* esd);
- Int_t CountAcceptedV0s(AliESDEvent* esd);
+ Bool_t AcceptV0(AliESDv0* const esdV0, AliESDtrack* const trackPos, AliESDtrack* const trackNeg, const AliESDVertex* esdVertex);
+ TObjArray* GetAcceptedV0s(const AliESD* esd);
+ Int_t CountAcceptedV0s(const AliESD* esd);
+ TObjArray* GetAcceptedV0s(const AliESDEvent* esd);
+ Int_t CountAcceptedV0s(const AliESDEvent* esd);
virtual Long64_t Merge(TCollection* list);
virtual void Copy(TObject &c) const;
TChain* AliEventPoolLoop::GetNextChain()
{
- //
+ // Get the next chain
if (fNIteration > fMaxIterations) {
return (0);
} else {
AliEventPoolLoop(const AliEventPoolLoop& obj);
AliEventPoolLoop& operator=(const AliEventPoolLoop& other);
protected:
- Int_t fMaxIterations;
- Int_t fNIteration;
+ Int_t fMaxIterations; // Maximum number of iterations
+ Int_t fNIteration; // Number of iterations
TChain* fChainClone; // Clone of the original
ClassDef(AliEventPoolLoop, 0);
};
void AliEventPoolOTF::Init()
{
- //
+ // Initialisation
if (!fGridTags) {
fTagAnalysis->ChainLocalTags(fTagDirectory);
} else {
TChain* AliEventPoolOTF::GetNextChain()
{
- //
+ // Get Next Chain
if (fChain) {
delete fChain;
fChain = 0;
void UserExec(Option_t *option);
void Terminate(Option_t *);
- void SetPriors(Double_t conc[AliPID::kSPECIES]) {for(Int_t i=0; i<AliPID::kSPECIES; i++) fPriors[i]=conc[i];}
- void GetPriors(Double_t conc[AliPID::kSPECIES]) {for(Int_t i=0; i<AliPID::kSPECIES; i++) conc[i]=fPriors[i];}
+ void SetPriors(Double_t conc[AliPID::kSPECIES]) {for(Int_t i=0; i<AliPID::kSPECIES; i++) fPriors[i]=conc[i];}
+ void GetPriors(Double_t conc[AliPID::kSPECIES]) const {for(Int_t i=0; i<AliPID::kSPECIES; i++) conc[i]=fPriors[i];}
void SetNiterations(Int_t nIter) {fNiterMax = nIter;}
TH1I *fHistEventsProcessed; // simple histo for monitoring the number of events processed
TH2D *fCalcPriors; // histo monitoring priors during iterations
- Double_t fPriors[AliPID::kSPECIES];
- Double_t fRecId[AliPID::kSPECIES];
- Double_t fMCId[AliPID::kSPECIES];
+ Double_t fPriors[AliPID::kSPECIES]; // Priors
+ Double_t fRecId[AliPID::kSPECIES]; // Reconstructed Id
+ Double_t fMCId[AliPID::kSPECIES]; // MC Id
Int_t fNiterations; //counter
Int_t fNiterMax; //max number of iterations
} while ((xevent = xml.GetNext(xevent)));
}
-Bool_t AliXMLCollection::GetCollectionSummary(Int_t *aTot, Int_t *aAcc, Int_t *aRejRun, Int_t *aRejLHC, Int_t *aRejDet, Int_t *aRejEvt)
+Bool_t AliXMLCollection::GetCollectionSummary(Int_t *aTot, Int_t *aAcc, Int_t *aRejRun, Int_t *aRejLHC, Int_t *aRejDet, Int_t *aRejEvt) const
{
// Return read list summary
*aTot = fTotalEvents;
TEntryList *GetEventList(const char *filename) const;
TEntryList *GetEntryList(const char *filename) { return GetEventList(filename); }
Bool_t OverlapCollection(TGridCollection * comparator);
- Bool_t GetCollectionSummary(Int_t *aTot, Int_t *aAcc, Int_t *aRejRun, Int_t *aRejLHC, Int_t *aRejDet, Int_t *aRejEvt);
+ Bool_t GetCollectionSummary(Int_t *aTot, Int_t *aAcc, Int_t *aRejRun, Int_t *aRejLHC, Int_t *aRejDet, Int_t *aRejEvt) const;
static AliXMLCollection *Open(const char *localcollectionfile);