fBackground(0),
fEnergyDistribution(0),
fEventSelectionEfficiency(0),
+ fSharingEfficiency(0),
fCorner1(4.2231, 26.6638),
fCorner2(1.8357, 27.9500),
fEnergyPath("$ALICE_ROOT/PWG2/FORWARD/corrections/EnergyDistribution"),
//____________________________________________________________________
char* AliFMDAnaParameters::GetPath(const char* species) {
- char* path ;
+ char* path = "";
if(species == fgkBackgroundID)
path = Form("%s/%s_%d_%d_%d_%d_%d_%d.root",
TH2F* hMult = 0;
TH2F* hHits = 0;
- TH2F* hHitsNoCuts = 0;
+ // TH2F* hHitsNoCuts = 0;
Int_t nVtxbins = pars->GetNvtxBins();
for(Int_t det =1; det<=3;det++)
hMultTotal->Add(hMultInput);
hMultTotal->Divide(hBg);//,"B");
- /*for(Int_t i = 1; i<=hTmp->GetNbinsX();i++) {
- for(Int_t j = 1; j<=hTmp->GetNbinsY();j++) {
- Float_t mult = hTmp->GetBinContent(i,j);
- if(mult == 0) continue;
- Float_t correction = hBg->GetBinContent(i,j);
-
- Float_t multcor = mult;
- if(correction != 0)
- multcor = multcor/correction;
- else
- std::cout<<"Warning! No correction for bin "<<i<<" , "<<j<<std::endl;
-
- hTmp->SetBinContent(i,j,multcor);
- }
- }*/
+
+ //sharing efficiency correction ?
+
}
}
if(fStandalone) {
fVtxCut(10),
fIsInit(kFALSE),
fIsGenerated(),
- fPrimEvents()
+ fPrimEvents(),
+ fEvents(),
+ fPrimdNdeta()
{
}
//_____________________________________________________________________
const char* AliFMDDndeta::GetAnalysisName(Analysis what, UShort_t det, Char_t ring, Int_t vtxbin) {
- const char* name;
+ const char* name = "";
switch(what) {
case kHits :
//_____________________________________________________________________
const char* AliFMDDndeta::GetPrimName(Analysis what, UShort_t det, Char_t ring, Int_t vtxbin) {
- const char* name;
+ const char* name = "";
switch(what) {
case kHits :
nNonZero++;
}
Int_t nBinsOld = fNbinsToCut;
- if(det == 2 && ringChar =='I') {
- fNbinsToCut = 0;
+ if(det == 2 && ringChar =='I') {
+ fNbinsToCut = 1;
}
TH1F* hRingMult = (TH1F*)fMultList.FindObject(Form("hRingMult_FMD%d%c",det,ringChar));
nNonZeroInData = 0;
TH1F* sumMultHist = (TH1F*)fMultList.FindObject(Form("hMult_vtxbin%d",v));
+
for(Int_t i =1;i<=sumMultHist->GetNbinsX();i++) {
if(multhistproj->GetBinContent(i) != 0 ) {
nNonZeroInData++;
if(nPrimevents > 0)
hPrim->Scale(1/nPrimevents);
- TFile testin("/home/canute/ALICE/FMDanalysis/GridAnalysis/fmd_dNdeta_hits.root","READ");
+ /*
+ TFile testin("/home/canute/ALICE/FMDanalysis/GridAnalysis/fmd_dNdeta_hits.root","READ");
// TFile testin("/home/canute/ALICE/FMDanalysis/productionData/fmd_dNdeta_hits.root","READ");
TH1F* hcorrect = (TH1F*)testin.Get("hReldif");
hcorrect->SetName("djarnis");
sumMult->Divide(hcorrect);
//delete hcorrect;
-
+ */
sumMult->Rebin(rebin);
sumMult->Scale(1/rebin);
//hPrim->GetYaxis()->SetRangeUser(0,20);
hPrim->SetStats(0000);
hPrim->GetXaxis()->SetTitle("#eta");
-
+ sumMult->SetTitle("");
+ sumMult->SetStats(kFALSE);
sumMult->DrawCopy("PE");
- hPrim->DrawCopy("E3same");
- //primMult->DrawCopy("E3same");
+ if(what != kHits)
+ hPrim->DrawCopy("E3same");
+ if(what == kHits)
+ primMult->DrawCopy("E3same");
hPrim->GetXaxis()->SetTitle("#eta");
hPrim->SetMarkerStyle(8);
-
- hPrim->DrawCopy("PE3same");
+ if(what != kHits)
+ hPrim->DrawCopy("PE3same");
sumMult->DrawCopy("PEsame");
std::cout<<"FMD 1 "<<sumMult->Integral(sumMult->FindBin(3),sumMult->FindBin(5),"width")<<std::endl;
/*TFile* itsfile = TFile::Open(itsfilename);
TLegend* leg = new TLegend(0.35,0.2,0.65,0.35,"");
- leg->AddEntry(hPrim,"Primary","pf");
+ if(what != kHits)
+ leg->AddEntry(hPrim,"Primary","pf");
+ else
+ leg->AddEntry(primMult,"Hits","pf");
//leg->AddEntry(primMult,"Analysed prim","f");
leg->AddEntry(sumMult,"Analysis","p");
leg->Draw();