// Whether to enable low flux specific code
task->SetEnableLowFlux(kFALSE);
- // --- Cuts --------------------------------------------------------
+ // --- Check for MC ------------------------------------------------
// Would like to use dynamic cast but CINT interprets that as a
// static cast - sigh!
- Bool_t mc = (task->IsA()==AliForwardMCMultiplicityTask::Class());
- // Double_t nXi = 2;
- // Bool_t includeSigma = false; //true;
- // Sharing cut
- AliFMDMultCuts cSharingLow;
- Double_t factor = 1.;
- cSharingLow.SetMultCuts(0.3, 0.3, 0.3, 0.3, 0.3);
- // Sharing cut
- AliFMDMultCuts cSharingHigh;
- cSharingHigh.SetMultCuts(-1);
- cSharingHigh.SetNXi(0); // Was 2
- cSharingHigh.SetIncludeSigma(false);
- cSharingHigh.SetMPVFraction(0.6);
- // Density cut
- AliFMDMultCuts cDensity;
- // cDensity.SetMultCuts(0.3, 0.3, 0.3, 0.3, 0.3);
- cDensity.SetMultCuts(-1);
- cDensity.SetMPVFraction(0.6); // Was .7
+ Bool_t mc = (task->IsA()==AliForwardMCMultiplicityTask::Class());
+
+ // --- Cuts --------------------------------------------------------
+ // Note, the absolute lowest signal to consider - ever -
+ // irrespective of MC or real data, is 0.15. Signals below this
+ // will contain remenants of the pedestal (yes, the width of the
+ // pedestal is small, but there are many _many_ channels with only
+ // pedestal value in them, so the absolute number of high-value
+ // pedestal signals is large - despite the low probablity).
+ AliFMDMultCuts cSharingLow(AliFMDMultCuts::kFixed,0.15);
+ AliFMDMultCuts cSharingHigh(AliFMDMultCuts::kMPVFraction,0.8,0.8,0.8,0.8,0.8);
+ AliFMDMultCuts cDensity(AliFMDMultCuts::kMPVFraction,0.7);
// --- Event inspector ---------------------------------------------
// Set the number of SPD tracklets for which we consider the event a
task->GetEventInspector().SetCentralityMethod("V0M");
// --- Sharing filter ----------------------------------------------
- // If the following is uncommented, then the merging of shared
+ // If the following is set to true, then the merging of shared
// signals is disabled completely
- // task->GetSharingFilter().SetDisableMerging(true);
+ task->GetSharingFilter().SetDisableMerging(false);
// Enable use of angle corrected signals in the algorithm
task->GetSharingFilter().SetUseAngleCorrectedSignals(true);
+ // Ignore the ESD information when angle correcting.
+ //
+ // *IMPORTANT*
+ //
+ // This is to counter a known issue with AliESDFMD with ClassDef
+ // version < 4, where the angle correction flag is incorrectly set.
+ // A fix is coming to AliESDFMD to handle it directly in the class.
+ // Only set the flag below to true if you know it to be necessary for
+ // your data set.
+ task->GetSharingFilter().SetIgnoreESDWhenAngleCorrecting(false);
// Disable use of angle corrected signals in the algorithm
task->GetSharingFilter().SetZeroSharedHitsBelowThreshold(false);
// Whether to use simple merging algorithm
task->GetSharingFilter().SetUseSimpleSharing(true);
- // Whether to allow for 3 strip hits
- task->GetSharingFilter().SetAllow3Strips(!mc);
- // Do not cut fixed/hard on multiplicity
- // task->GetSharingFilter().GetHCuts().SetMultCuts(-1);
- // Set the number of xi's (width of landau peak) to stop at
- // task->GetSharingFilter().GetHCuts().SetNXi(nXi);
- // Set whether or not to include sigma in cut
- // task->GetSharingFilter().GetHCuts().SetIncludeSigma(includeSigma);
+ // Whether to allow for 3 strip hits - deprecated
+ task->GetSharingFilter().SetAllow3Strips(false);
// Set upper sharing cut
task->GetSharingFilter().SetHCuts(cSharingHigh);
// Enable use of angle corrected signals in the algorithm
// If this option is enabled, then the summed per-vertex, per-ring
// d2N/detadphi histograms will be stored in the output, as well as
// copies of the secondary maps
- // task->GetHistCollector().SetMakeBGHitMaps(true);
+ task->GetHistCollector().SetMakeBGHitMaps(false);
//
// If this option is enabled, then a 3D histogram will be made for
// each ring, summing dN/deta for each centrality bin.
- // task->GetHistCollector().SetMakeCentralitySums(true);
+ task->GetHistCollector().SetMakeCentralitySums(false);
// --- Eventplane Finder -------------------------------------------
task->GetEventPlaneFinder().SetUsePhiWeights(false);
// output AOD - one for each FMD ring. The branches each contain a
// TH2D object of the (primary) charged particle multiplicity per
// (eta,phi)-bin in that event
- // task->SetStorePerRing(true);
+ task->SetStorePerRing(false);
// --- Set limits on fits the energy -------------------------------
// DO NOT CHANGE THESE UNLESS YOU KNOW WHAT YOU ARE DOING