}
//____________________________________________________________________
-TH1F* AliFMDAnaCalibBackgroundCorrection::GetDoubleHitCorrection(Int_t det,
+TH1F* AliFMDAnaCalibBackgroundCorrection::GetDoubleHitCorrection(Int_t det,
Char_t ring) {
- TH1F* hCorrection = (TH1F*)fListOfDoubleHitCorrection.FindObject(Form("hDoubleHitCorrection_FMD%d%c"));
+ TH1F* hCorrection = (TH1F*)fListOfDoubleHitCorrection.FindObject(Form("hDoubleHitCorrection_FMD%d%c",det,ring));
return hCorrection;
}
void AliFMDAnaCalibBackgroundCorrection::SetDoubleHitCorrection(Int_t det,
Char_t ring,
TH1F* hCorrection) {
- hCorrection->SetName(Form("hDoubleHitCorrection_FMD%d%c"));
+ hCorrection->SetName(Form("hDoubleHitCorrection_FMD%d%c",det,ring));
fListOfDoubleHitCorrection.Add(hCorrection);
}
//____________________________________________________________________
: TNamed(o),
fIsInit(o.fIsInit),
fBackground(o.fBackground),
- fEnergyDistribution(o.fEnergyDistribution)
- // fBackgroundArray(o.fBackgroundArray),
- //fEdistArray(o.fEdistArray)
+ fEnergyDistribution(o.fEnergyDistribution),
+ fCorner1(o.fCorner1),
+ fCorner2(o.fCorner2),
+ fEnergyPath(o.fEnergyPath),
+ fBackgroundPath(o.fBackgroundPath)
{}
AliFMDAnaParameters& operator=(const AliFMDAnaParameters&) { return *this; }
virtual ~AliFMDAnaParameters() {}
DefineOutput(0, TList::Class());
}
- fFuncPos = new TF1("funcPos","pol1",0,6);
+ /* fFuncPos = new TF1("funcPos","pol1",0,6);
fFuncPos->SetParameters(0.99925,0.00298301);
fFuncNeg = new TF1("funcNeg","pol1",-6,0);
- fFuncNeg->SetParameters(0.987583,-0.0101022);
+ fFuncNeg->SetParameters(0.987583,-0.0101022);*/
}
// std::cout<<eta<<" "<<etacalc<<std::endl;
//eta = etacalc;
- Float_t m = pars->GetMPV(det,ring,eta);
- Float_t s = pars->GetSigma(det,ring,eta);
+ // Float_t m = pars->GetMPV(det,ring,eta);
+ // Float_t s = pars->GetSigma(det,ring,eta);
//AliFMDParameters* recopars = AliFMDParameters::Instance();
Float_t mult_cut = 0.15;//m-2*s;//0.15;//0.2;//m-3*s;// 0.2;//0.01;//m-2*s;//0.2;
if(ring == 'I')
mult_cut = 0.10;
+
//mult_cut = (5*recopars->GetPedestalWidth(det,ring,sec,strip))/(recopars->GetPulseGain(det,ring,sec,strip)*recopars->GetDACPerMIP());
Float_t nParticles = 0;
if(fESD->GetUniqueID() == kTRUE) {
TH1F* hDoubleHitCorrection = pars->GetDoubleHitCorrection(det,ring);
//if(det == 3)
+
if(hDoubleHitCorrection->GetBinContent(hDoubleHitCorrection->FindBin(eta)) != 0)
correction = correction / hDoubleHitCorrection->GetBinContent(hDoubleHitCorrection->FindBin(eta));
+
//else
// correction = correction / pars->GetDoubleHitCorrection()->GetBinContent(pars->GetDoubleHitCorrection()->FindBin(eta));
}
AliESDVertex* fVertex;
Bool_t fStandalone;
Bool_t fStatus;
- TF1* fFuncPos;
- TF1* fFuncNeg;
+ // TF1* fFuncPos;
+ // TF1* fFuncNeg;
ClassDef(AliFMDAnalysisTaskDensity, 0); // Analysis task for FMD analysis
};
//_____________________________________________________________________
AliFMDAnalysisTaskGenerateBackground::AliFMDAnalysisTaskGenerateBackground():
-AliAnalysisTaskSE()
+AliAnalysisTaskSE(),
+ fListOfHits(),
+ fListOfPrimaries(),
+ fListOfCorrection(),
+ fVertexBins(),
+ fLastTrackByStrip(),
+ fHitsByStrip(),
+ fZvtxCut(10),
+ fNvtxBins(10),
+ fNbinsEta(200),
+ fBackground(0)
{
// Default constructor
}
//_____________________________________________________________________
AliFMDAnalysisTaskGenerateBackground::AliFMDAnalysisTaskGenerateBackground(const char* name):
AliAnalysisTaskSE(name),
+ fListOfHits(),
+ fListOfPrimaries(),
+ fListOfCorrection(),
+ fVertexBins(),
+ fLastTrackByStrip(),
+ fHitsByStrip(),
fZvtxCut(10),
fNvtxBins(10),
- fNbinsEta(200)
+ fNbinsEta(200),
+ fBackground(0)
{
DefineOutput(1, TList::Class());
}
//_____________________________________________________________________
-void AliFMDAnalysisTaskGenerateBackground::ReadFromFile(const Char_t* filename, Bool_t storeInOCDB, Int_t runNo) {
+void AliFMDAnalysisTaskGenerateBackground::ReadFromFile(const Char_t* filename, Bool_t storeInOCDB, Int_t /*runNo*/) {
TFile infile(filename);
TH1F* hVertex = (TH1F*)infile.Get("VertexBins");
fListOfCorrection(),
fVertexBins(o.fVertexBins),
fLastTrackByStrip(o.fLastTrackByStrip),
+ fHitsByStrip(o.fHitsByStrip),
fZvtxCut(o.fZvtxCut),
fNvtxBins(o.fNvtxBins),
fNbinsEta(o.fNbinsEta),
Float_t Enext,
UShort_t det,
Char_t ring,
- UShort_t sec,
- UShort_t strip) {
+ UShort_t /*sec*/,
+ UShort_t /*strip*/) {
AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
Float_t merged_energy = 0;