188 RhoCandList neutralCands,chargedCands, plusCands,minusCands;
192 RhoCandList xiCands,nonOvCands,dsCands,ds0Cands,ppCands;
193 std::map<Int_t,Float_t > mapp;
203 for (Int_t i1=0; i1<
fMicroArray->GetEntriesFast(); i1++){
206 TLorentzVector l=tc.P4();
210 chargedCands.
Add(&tc);
231 ppiCands.
Combine(minusCands,plusCands);
235 for (
int la=0;la<ppiCands.
GetLength();la++)
237 Lamb->Fill((ppiCands.
Get(la)->
P4()).M());
241 xiCands.
Combine(ppiCands,minusCands);
245 while (t4=itX.Next())
274 while (t2=iterP.Next())
278 TLorentzVector l=t2->
P4();
280 TVector3 pim = l.Vect();
300 for (
int k=0;k<
npi;k++)
320 while (t3=itP.Next())
324 TLorentzVector v4=t3->
P4();
326 TVector3 v3 = v4.Vect();
329 ppi2->Fill(v3.Mag());
384 TLorentzVector vpim=pi->
P4();
385 TVector3 pi3v = vpim.Vect();
388 if(pi3v.Mag()<0.2)
hvtx2[0]->Fill(pi3v.Mag(),0.);
395 for (ii=0;ii<dsi-1;ii++)
401 TLorentzVector vpim=pi->
P4();
402 TVector3 pi3v = vpim.Vect();
404 for (jj=ii+1;jj<dsi;jj++)
409 TLorentzVector vpp=pp->
P4();
410 TVector3 pp3v = vpp.Vect();
421 if (MotherId==-1)Motherpdg = mc->
GetPdgCode();
438 cout<<
" Motherpdg "<<Motherpdg<<
" "<<mc->
GetPdgCode()<<endl;
452 if(pi3v.Mag()>pp3v.Mag())
455 hvtx2[0]->Fill(pi3v.Mag(),pp3v.Mag());
457 if((pi3v.Mag()>0.12&&pi3v.Mag()<0.14)&&(pp3v.Mag()>0.065&&pp3v.Mag()<0.08)) {
462 hvtx2[1]->Fill(pi3v.Mag(),pp3v.Mag());
474 if(Motherpdg==1020040110||Motherpdg==1010050110){
478 if(Motherpdg==1020030090||(Motherpdg==1010040090 && MotherId!=-1)){
489 if((pi3v.Mag()>0.112&&pi3v.Mag()<0.126)&&(pp3v.Mag()>0.09&&pp3v.Mag()<0.103)){
490 hvtx2[2]->Fill(pi3v.Mag(),pp3v.Mag());
498 if((pi3v.Mag()>0.128&&pi3v.Mag()<0.147)&&(pp3v.Mag()>0.0898&&pp3v.Mag()<0.109)){
499 hvtx2[9]->Fill(pi3v.Mag(),pp3v.Mag());
513 if(Motherpdg==1020040100||Motherpdg==1010050100){
520 if((pi3v.Mag()>0.097&&pi3v.Mag()<0.106)&&(pp3v.Mag()>0.094&&pp3v.Mag()<0.103)) {
524 hvtx2[3]->Fill(pi3v.Mag(),pp3v.Mag());
538 if(Motherpdg==1020040120||Motherpdg==1010050120){
546 if((pi3v.Mag()>0.128&&pi3v.Mag()<0.147)&&(pp3v.Mag()>0.110&&pp3v.Mag()<0.124)) {
551 hvtx2[4]->Fill(pi3v.Mag(),pp3v.Mag());
562 if(Motherpdg==1020020060||Motherpdg==1010030060){
570 if((pi3v.Mag()>0.128&&pi3v.Mag()<0.147)&&(pp3v.Mag()>0.124&&pp3v.Mag()<0.143)) {
575 hvtx2[5]->Fill(pi3v.Mag(),pp3v.Mag());
594 if(pi3v.Mag()<pp3v.Mag()) {
597 hvtx2[0]->Fill(pp3v.Mag(),pi3v.Mag());
600 if(Motherpdg==1020040110||Motherpdg==1010050110){
606 if((pp3v.Mag()>0.12&&pp3v.Mag()<0.14)&&(pi3v.Mag()>0.065&&pi3v.Mag()<0.08)) {
609 hvtx2[1]->Fill(pp3v.Mag(),pi3v.Mag());
620 if(Motherpdg==1020030090||(Motherpdg==1010040090&& MotherId!=-1)){
625 if((pp3v.Mag()>0.112&&pp3v.Mag()<0.126)&&(pi3v.Mag()>0.09&&pi3v.Mag()<0.103))
629 hvtx2[2]->Fill(pp3v.Mag(),pi3v.Mag());
640 if((pp3v.Mag()>0.128&&pp3v.Mag()<0.147)&&(pi3v.Mag()>0.0898&&pi3v.Mag()<0.109))
642 hvtx2[9]->Fill(pp3v.Mag(),pi3v.Mag());
653 if(Motherpdg==1020040100||Motherpdg==1010050100){
659 if((pp3v.Mag()>0.097&&pp3v.Mag()<0.106)&&(pi3v.Mag()>0.094&&pi3v.Mag()<0.103)) {
662 hvtx2[3]->Fill(pp3v.Mag(),pi3v.Mag());
672 if(Motherpdg==1020040120||Motherpdg==1010050120){
678 if((pp3v.Mag()>0.128&&pp3v.Mag()<0.147)&&(pi3v.Mag()>0.110&&pi3v.Mag()<0.124)) {
681 hvtx2[4]->Fill(pp3v.Mag(),pi3v.Mag());
693 if(Motherpdg==1020020060||Motherpdg==1010030060){
701 if((pp3v.Mag()>0.128&&pp3v.Mag()<0.147)&&(pi3v.Mag()>0.124&&pi3v.Mag()<0.143)) {
706 hvtx2[5]->Fill(pp3v.Mag(),pi3v.Mag());
void Add(const RhoCandidate *c)
RhoMinusParticleSelector * minusSel
RhoPlusParticleSelector * plusSel
RhoSimpleProtonSelector * pSel
PndPidCandidate * GetRecoCandidate() const
RhoMassParticleSelector * LambMSel
void SetEnergySpectra(int event, int cluster)
void Combine(RhoCandList &l1, RhoCandList &l2)
TClonesArray * fMicroArray
void Select(RhoParticleSelectorBase *pidmgr)
TLorentzVector P4() const
static RhoFactory * Instance()
void SetTotESpectra(int clus)
TVector3 GetStartVertex() const
Int_t GetMotherID() const
RhoSimplePionSelector * piSel
RhoCandidate * Get(Int_t)