if( ptK <0.300) continue;
- Float_t code = particle->GetPdgCode();
+ Int_t code = particle->GetPdgCode();
Int_t mcProcess=-1011;
if ((code==321)||(code==-321)){
for (Int_t k=firstD;k<=lastD;k++) {
if ( k > 0 ) {
TParticle* daughter1=stack->Particle(k); // 27/8
- Float_t dcode = daughter1->GetPdgCode();
+ Int_t dcode = daughter1->GetPdgCode();
mcProcess=daughter1->GetUniqueID();
if (mcProcess==4) {
TH2F *fZvYv; // two dime of Z vrs Y of vtx main
TH2F *fXvYv; // two dime of X vrs Y of main tracks vtx main
TH1F *fPtPrKink; // pt of Primary PDG kaons inside the selected ones by the kink topology
- TF1 *f1;
- TF1 *f2;
+ TF1 *f1; // upper limit curve for the decay K->mu
+ TF1 *f2; // upper limit curve for the decay pi->mu
TList *fListOfHistos; // list of histos
AliAnalysisKinkESDMC(const AliAnalysisKinkESDMC&); // not implemented
private:
TList *fList; // List
- AliResonanceKink *fKinkResonance;
+ AliResonanceKink *fKinkResonance; // resonance object configured externaly
AliAnalysisTaskKinkResonance(const AliAnalysisTaskKinkResonance&); // not implemented
AliAnalysisTaskKinkResonance& operator=(const AliAnalysisTaskKinkResonance&); // not implemented
continue;
}
- if(particle->GetNDaughters()>2) continue;
-
if(TMath::Abs(particle->GetPdgCode())==fresonancePDGcode) {
+
+ if(particle->GetNDaughters()>2) continue;
Int_t firstD=particle->GetFirstDaughter();
Int_t lastD=particle->GetLastDaughter();
TH2D *fSimEtaPtKink; // Eta pT Sim one kaon kink spectrum
TH1D *fhdr ; // radial impact
TH1D *fhdz ; // z impact
- TF1 *f1;
- TF1 *f2;
+ TF1 *f1; //upper limit curve for the decay K->mu
+ TF1 *f2; //upper limit curve for the decay pi->mu
TString fAnalysisType;//"ESD" or "MC"
TH1D *fvtxz ; // vtx z component
Int_t fNbins; // bins
Float_t fLowX;// lowx
Float_t fHighX; // high x
- Int_t fdaughter1pdg;
- Int_t fdaughter2pdg;
- Int_t fresonancePDGcode;
- Float_t fMaxNSigmaToVertex;
- Double_t fMinPtTrackCut;
- Double_t fMaxDCAxy;
- Double_t fMaxDCAzaxis;
- Int_t fMinTPCclusters;
- Double_t fMaxChi2PerTPCcluster;
- Double_t fMaxCov0;
- Double_t fMaxCov2;
- Double_t fMaxCov5;
- Double_t fMaxCov9;
- Double_t fMaxCov14;
- TH2D *fInvmassPt;
- TH2D *fInvmassPtTrue;
- TH2D *fMCInvmassPt;
- TH2D *fMCInvmassPtTrue;
+ Int_t fdaughter1pdg; // pdg code of the resonance's first daughter
+ Int_t fdaughter2pdg; // pdg code of the resonance's second daughter
+ Int_t fresonancePDGcode; // pdg code of the resonance
+ Float_t fMaxNSigmaToVertex; // standard cut to select primary tracks
+ Double_t fMinPtTrackCut; // lower pt cut for the tracks
+ Double_t fMaxDCAxy; // impact parameter in the xy plane
+ Double_t fMaxDCAzaxis; // impact parameter in the z axis
+ Int_t fMinTPCclusters; // standard cut for the TPC clusters
+ Double_t fMaxChi2PerTPCcluster; // standard cut for the chi2 of the TPC clusters
+ Double_t fMaxCov0; // standard cut
+ Double_t fMaxCov2; // standard cut
+ Double_t fMaxCov5; // standard cut
+ Double_t fMaxCov9; // standard cut
+ Double_t fMaxCov14; // standard cut
+ TH2D *fInvmassPt; // 2D histo for invariant mass calculation in pt bins (all pairs, ESD)
+ TH2D *fInvmassPtTrue; // 2D histo for invariant mass calculation in pt bins (true pairs, ESD)
+ TH2D *fMCInvmassPt; // 2D histo for invariant mass calculation in pt bins (all pairs, MC)
+ TH2D *fMCInvmassPtTrue; // 2D histo for invariant mass calculation in pt bins (true pairs, MC)
// Bool_t fTPCrefitFlag;
AliResonanceKink(const AliResonanceKink&); // not implemented
//Uncomment the following block if the Like Sign is made of K+ kink + positive track
- Int_t IndexKink=trackpos->GetKinkIndex(0);
+ Int_t indexKink=trackpos->GetKinkIndex(0);
Int_t kaonKinkFlag=0;
- if(IndexKink<0){
+ if(indexKink<0){
- AliESDkink *kink=fESD->GetKink(TMath::Abs(IndexKink)-1);
+ AliESDkink *kink=fESD->GetKink(TMath::Abs(indexKink)-1);
const TVector3 motherMfromKink(kink->GetMotherP());
const TVector3 daughterMKink(kink->GetDaughterP());
Float_t qt=kink->GetQt();
private:
AliESDEvent *fESD; //! ESD object
TList *fListOfHistos; //! List
- TF1 *f1;
- TF1 *f2;
+ TF1 *f1; //upper limit curve for the decay K->mu
+ TF1 *f2; //upper limit curve for the decay pi->mu
TH1D *fPosKaonLikeSign; //! negative spectrum
TH2D *fLikeSignInvmassPt; //! negative spectrum