if(fSystem==2) nBinscb=21;
Double_t maxncn=nBinscb-0.5;
fHistoNcharmed = new TH2F("hncharmed","",100,0.,maxMult,nBinscb,-0.5,maxncn);
if(fSystem==2) nBinscb=21;
Double_t maxncn=nBinscb-0.5;
fHistoNcharmed = new TH2F("hncharmed","",100,0.,maxMult,nBinscb,-0.5,maxncn);
fHistOriginFeeddown->Sumw2();
fHistOriginFeeddown->SetMinimum(0);
fOutput->Add(fHistOriginFeeddown);
fHistOriginFeeddown->Sumw2();
fHistOriginFeeddown->SetMinimum(0);
fOutput->Add(fHistOriginFeeddown);
+ fHistMotherID=new TH1F("hMotherID","",1000,-1.5,998.5);
+ fHistMotherID->SetMinimum(0);
+ fOutput->Add(fHistMotherID);
Double_t disty=part->Vy()-mcVert->GetY();
Double_t distz=part->Vz()-mcVert->GetZ();
Double_t distToVert=TMath::Sqrt(distx*distx+disty*disty+distz*distz);
Double_t disty=part->Vy()-mcVert->GetY();
Double_t distz=part->Vz()-mcVert->GetZ();
Double_t distToVert=TMath::Sqrt(distx*distx+disty*disty+distz*distz);