AliJetModelBaseTask("AliJetEmbeddingFromGenTask"),
fGen(0),
fMassless(kFALSE),
+ fChargedOnly(kFALSE),
fHistPt(0),
+ fHistEtaPhi(0),
fHistTrials(0),
fHistXsection(0),
fHistPtHard(0)
AliJetModelBaseTask(name,drawqa),
fGen(0),
fMassless(kFALSE),
+ fChargedOnly(kFALSE),
fHistPt(0),
+ fHistEtaPhi(0),
fHistTrials(0),
fHistXsection(0),
fHistPtHard(0)
fHistPt = new TH1F("fHistpt","fHistPt;#it{p}_{T};N",100,0.,100.);
fOutput->Add(fHistPt);
+ fHistEtaPhi = new TH2F("fHistEtapHI","fHistEtaPhi;#eta;#varphi",100,-3.,3.,100.,0.,TMath::TwoPi());
+ fOutput->Add(fHistEtaPhi);
+
fHistTrials = new TH1F("fHistTrials", "fHistTrials", 1, 0, 1);
fHistTrials->GetYaxis()->SetTitle("trials");
fOutput->Add(fHistTrials);
if (!pdg)
continue;
Int_t c = (Int_t)(TMath::Abs(pdg->Charge()));
- if (c==0)
- continue;
+ if (fChargedOnly && c==0) continue;
Double_t pt = part->Pt();
Double_t eta = part->Eta();
Double_t phi = part->Phi();
Double_t mass = part->GetMass();
if(fMassless) mass = 0.;
fHistPt->Fill(pt);
+ fHistEtaPhi->Fill(eta,phi);
AddTrack(pt, eta, phi,0,0,0,0,0,0,c,mass);
}