Photospp::PhotosBranch Class Reference

#include <PhotosBranch.h>

Public Member Functions
	PhotosBranch (PhotosParticle *p)
PhotosParticle *	getDecayingParticle ()
vector< PhotosParticle * >	getMothers ()
vector< PhotosParticle * >	getDaughters ()
vector< PhotosParticle * >	getParticles ()
int	getSuppressionStatus ()
int	getForcingStatus ()
bool	checkMomentumConservation ()
void	process ()
	PhotosBranch (PhotosParticle *p)
PhotosParticle *	getDecayingParticle ()
vector< PhotosParticle * >	getMothers ()
vector< PhotosParticle * >	getDaughters ()
vector< PhotosParticle * >	getParticles ()
int	getSuppressionStatus ()
int	getForcingStatus ()
bool	checkMomentumConservation ()
void	process ()

Static Public Member Functions
static vector< PhotosBranch * >	createBranches (vector< PhotosParticle * > particles)
static vector< PhotosBranch * >	createBranches (vector< PhotosParticle * > particles)

Private Member Functions
int	checkSuppressionLevel ()
int	checkForcingLevel ()
int	checkList (bool forceOrSuppress)
int	checkSuppressionLevel ()
int	checkForcingLevel ()
int	checkList (bool forceOrSuppress)

Private Attributes
int	suppression
int	forcing
PhotosParticle *	particle
vector< PhotosParticle * >	mothers
vector< PhotosParticle * >	daughters

Detailed Description

Definition at line 23 of file PhotosBranch.h.

Constructor & Destructor Documentation

PhotosBranch::PhotosBranch

(

PhotosParticle *

p

)

Create branch out of decaying particle

Definition at line 15 of file PhotosBranch.cxx.

References checkForcingLevel(), checkMomentumConservation(), checkSuppressionLevel(), daughters, Photospp::Log::Debug(), Photospp::PhotosParticle::findProductionMothers(), forcing, Photospp::PhotosParticle::getDaughters(), mothers, p, particle, suppression, and Photospp::Log::Warning().

Referenced by createBranches().

{
        daughters = p->getDaughters();

        //Suppress if somehow got stable particle
        if(daughters.size()==0)
        {
                Log::Debug(1)<<"Stable particle."<<endl;
                suppression = 1;
                forcing     = 0;
                particle    = NULL;
                return;
        }
        else if(daughters.at(0)->getMothers().size()==1)
        {
                // Regular case - one mother
                Log::Debug(1)<<"Regular case."<<endl;
                particle  = p;
                mothers   = p->findProductionMothers();
        }
        else
        {
                // Advanced case - branch with multiple mothers - no mid-particle
                Log::Debug(1)<<"Advanced case."<<endl;
                particle  = NULL;
                mothers   = daughters.at(0)->getMothers();
        }
        forcing = checkForcingLevel();
        if(!forcing) suppression = checkSuppressionLevel();
        else         suppression = 0;

        // Even if forced or passed suppression check, we still have to check few things
        if(!suppression)
        {
                // Check momentum conservation
                suppression=!checkMomentumConservation();
                if(suppression) Log::Warning()<<"Branching ignored due to 4-momentum non conservation"<<endl;

                // Check if advanced case has only one daughter
                if(!particle && daughters.size()==1) suppression=-1;

                // If any of special cases is true, we're not forcing this branch
                if(suppression) forcing=0;
        }
}

Photospp::PhotosBranch::PhotosBranch

(

PhotosParticle *

p

)

Create branch out of decaying particle

Member Function Documentation

int Photospp::PhotosBranch::checkForcingLevel ( )

inlineprivate

Checks if branching is forced by PHOTOS.

Definition at line 61 of file PhotosBranch.h.

References checkList().

Referenced by PhotosBranch().

{ return checkList(true);  }

int Photospp::PhotosBranch::checkForcingLevel ( )

inlineprivate

Checks if branching is forced by PHOTOS.

Definition at line 61 of file 3.52/src/photosCInterfaces/PhotosBranch.h.

References checkList().

{ return checkList(true);  }

int PhotosBranch::checkList ( bool  forceOrSuppress )

private

Algorithm used for suppression/forcing check

Definition at line 206 of file PhotosBranch.cxx.

References daughters, Photospp::Photos::forceBremList, Photospp::PhotosParticle::getPdgID(), mothers, particle, and Photospp::Photos::supBremList.

Referenced by checkForcingLevel(), and checkSuppressionLevel().

{
        vector< vector<int>* > *list = (forceOrSuppress) ? Photos::forceBremList : Photos::supBremList;
        if(!list) return 0;

        // Can't check without pdgid
        int motherID;
        if(particle) motherID = particle->getPdgID();
        else
        {
                if(mothers.size()==0) return 0;
                motherID = mothers.at(0)->getPdgID();
        }

        // Create list of daughters
        vector<int> dID;
        for(int j=0;j<(int)daughters.size();j++) dID.push_back(daughters[j]->getPdgID());

        vector< vector<int> *> &patternList = *list;

        // Check if the mother and list of daughters matches any of the declared patterns
        for(int j=0; j<(int)patternList.size();j++)
        {
                // Skip patterns that don't have our mother
                if(motherID!=(*patternList[j])[0]) continue;

                // Compare decay daughters with pattern - max O(n*m)
                vector<int> &pattern = *patternList[j];
                bool fullMatch=true;
                for(int k = 1; k<(int)pattern.size()-1; k++)
                {
                        bool oneMatch=false;
                        for(int l=0;l<(int)dID.size(); l++)
                                if(pattern[k]==dID[l]) { oneMatch=true; break; }
                        if(!oneMatch) { fullMatch=false; break; }
                }
                // Check if the matching pattern is set for consecutive suppression
                /*
                  Currently minimal matching is used.
                  If e.g. 25 -> -15 is suppressed, then 25 -> 15,-15 and 25 -> 15,-15,22 etc. is suppressed too
                  For exact matching (then suppress 25 -> 15,-15 ; 25 -> 15,-15,22 etc. must be done separately) uncoment line ...:
                */
                // if(pattern.size()<=2 || (fullMatch && dID.size()==pattern.size()-2) )
                // ...and comment out line:
                if(pattern.size()<=2 || fullMatch)
                        return (pattern.back()==1) ? 2 : 1;
        }
        return 0;
}

int Photospp::PhotosBranch::checkList ( bool  forceOrSuppress )

private

Algorithm used for suppression/forcing check

bool Photospp::PhotosBranch::checkMomentumConservation

(

)

Checks momentum conservation of decaying particle. If it does not exist, checks momentum of first mother passed to photos

bool PhotosBranch::checkMomentumConservation

(

)

Checks momentum conservation of decaying particle. If it does not exist, checks momentum of first mother passed to photos

Definition at line 85 of file PhotosBranch.cxx.

References Photospp::PhotosParticle::checkMomentumConservation(), mothers, and particle.

Referenced by PhotosBranch(), and process().

{
        if(particle)           return particle->checkMomentumConservation();
        if(mothers.size()>0)   return mothers.at(0)->checkMomentumConservation();
        return true;
}

int Photospp::PhotosBranch::checkSuppressionLevel ( )

inlineprivate

Checks if branching is suppressed by PHOTOS.

Definition at line 58 of file PhotosBranch.h.

References checkList().

Referenced by PhotosBranch().

{ return checkList(false); }

int Photospp::PhotosBranch::checkSuppressionLevel ( )

inlineprivate

Checks if branching is suppressed by PHOTOS.

Definition at line 58 of file 3.52/src/photosCInterfaces/PhotosBranch.h.

References checkList().

{ return checkList(false); }

static vector<PhotosBranch *> Photospp::PhotosBranch::createBranches ( vector< PhotosParticle * >  particles )

static

Create branches from particles list

vector< PhotosBranch * > PhotosBranch::createBranches ( vector< PhotosParticle * >  particles )

static

Create branches from particles list

Definition at line 92 of file PhotosBranch.cxx.

References b, Photospp::Log::Debug(), Photospp::Photos::forceBremList, forcing, Photospp::PhotosParticle::getBarcode(), getDecayingParticle(), Photospp::PhotosParticle::getDecayTree(), getForcingStatus(), getMothers(), Photospp::PhotosParticle::getPdgID(), getSuppressionStatus(), i, Photospp::Photos::isSuppressed, m, mothers, p, particle, PhotosBranch(), suppression, and tree.

Referenced by Photospp::PhotosEvent::process(), and Photospp::Photos::processBranch().

{
        Log::Debug(700)<<"PhotosBranch::createBranches - filtering started"<<endl;
        list<PhotosParticle *> list(particles.begin(),particles.end());
        vector<PhotosBranch *> branches;

        // First - add all forced decays
        if(Photos::forceBremList)
        {
                std::list<PhotosParticle *>::iterator it;
                for(it=list.begin();it!=list.end();it++)
                {
                        PhotosBranch *branch = new PhotosBranch(*it);
                        int forcing = branch->getForcingStatus();
                        if(forcing)
                        {
                                Log::Debug(701)<<" Forced: "<<(*it)->getPdgID()<<" (barcode: "<<(*it)->getBarcode()<<") with forcing status= "<<forcing<<endl;
                                branches.push_back(branch);
                                it = list.erase(it);
                                --it;
                                // If forcing consecutive decays
                                if(forcing==2)
                                {
                                        PhotosParticle *p = branch->getDecayingParticle();
                                        if(!p){
                                                if(branch->getMothers().size()>0) p = branch->getMothers().at(0);
                                                }else {continue;}
                                        vector<PhotosParticle *> tree = p->getDecayTree();
                                        //Add branches for all particles from the list - max O(n*m)
                                        std::list<PhotosParticle *>::iterator it2;
                                        for(it2=list.begin();it2!=list.end();it2++)
                                        {
                                                for(int i=0;i<(int)tree.size();i++)
                                                {
                                                        if(tree.at(i)->getBarcode()==(*it2)->getBarcode())
                                                        {
                                                                PhotosBranch *b = new PhotosBranch(*it2);
                                                                branches.push_back(b);
                                                                // If we were to delete our next particle in line
                                                                if(it==it2) --it;
                                                                it2 = list.erase(it2);
                                                                --it2;
                                                                break;
                                                        }
                                                }
                                        }
                                }
                        }
                        else delete branch;
                }
        }
        // Quit if we're suppressing everything
        if(Photos::isSuppressed) return branches;
        // Now - check if remaining decays are suppressed
        while(!list.empty())
        {
                PhotosParticle *particle = list.front();
                list.pop_front();
                if(!particle) continue;

                PhotosBranch *branch = new PhotosBranch(particle);
                int suppression = branch->getSuppressionStatus();
                if(!suppression) branches.push_back(branch);
                else
                {
                        Log::Debug(702)<<"  Suppressed: "<<particle->getPdgID()<<" (barcode: "<<particle->getBarcode()<<") with suppression status= "<<suppression<<endl;
                        //If suppressing consecutive decays
                        if(suppression==2)
                        {
                                PhotosParticle *p = branch->getDecayingParticle();
                                if(!p){
                                        if(branch->getMothers().size()>0) p = branch->getMothers().at(0);
                                        }else {continue;}
                                vector<PhotosParticle *> tree = p->getDecayTree();
                                //Remove all particles from the list - max O(n*m)
                                std::list<PhotosParticle *>::iterator it;
                                for(it=list.begin();it!=list.end();it++)
                                {
                                        for(int i=0;i<(int)tree.size();i++)
                                        {
                                                if(tree.at(i)->getBarcode()==(*it)->getBarcode())
                                                {
                                                        it = list.erase(it);
                                                        --it;
                                                        break;
                                                }
                                        }
                                }
                        }
                        delete branch;
                        continue;
                }

                //In case we don't have mid-particle erase rest of the mothers from list
                if(!branch->getDecayingParticle())
                {
                        vector<PhotosParticle *> mothers = branch->getMothers();
                        for(int i=0;i<(int)mothers.size();i++)
                        {
                                PhotosParticle *m = mothers.at(i);
                                if(m->getBarcode()==particle->getBarcode()) continue;
                                std::list<PhotosParticle *>::iterator it;
                                for(it=list.begin();it!=list.end();it++)
                                        if(m->getBarcode()==(*it)->getBarcode())
                                        {
                                                it = list.erase(it);
                                                break;
                                        }
                        }
                }
        }
        return branches;
}

vector<PhotosParticle *> Photospp::PhotosBranch::getDaughters ( )

inline

Get list of daughters

Definition at line 36 of file PhotosBranch.h.

References daughters.

Referenced by Photospp::PH_HEPEVT_Interface::set().

{ return daughters; }

vector<PhotosParticle *> Photospp::PhotosBranch::getDaughters ( )

inline

Get list of daughters

Definition at line 36 of file 3.52/src/photosCInterfaces/PhotosBranch.h.

References daughters.

{ return daughters; }

PhotosParticle* Photospp::PhotosBranch::getDecayingParticle ( )

inline

Return decaying particle. NULL if branching does not have mid-particle

Definition at line 30 of file 3.52/src/photosCInterfaces/PhotosBranch.h.

References particle.

{ return particle;  }

PhotosParticle* Photospp::PhotosBranch::getDecayingParticle ( )

inline

Return decaying particle. NULL if branching does not have mid-particle

Definition at line 30 of file PhotosBranch.h.

References particle.

Referenced by createBranches(), and Photospp::PH_HEPEVT_Interface::set().

{ return particle;  }

int Photospp::PhotosBranch::getForcingStatus ( )

inline

Check if branch is forced

Definition at line 45 of file PhotosBranch.h.

References forcing.

Referenced by createBranches().

{ return forcing; }

int Photospp::PhotosBranch::getForcingStatus ( )

inline

Check if branch is forced

Definition at line 45 of file 3.52/src/photosCInterfaces/PhotosBranch.h.

References forcing.

{ return forcing; }

vector<PhotosParticle *> Photospp::PhotosBranch::getMothers ( )

inline

Get list of mothers

Definition at line 33 of file PhotosBranch.h.

References mothers.

Referenced by createBranches(), and Photospp::PH_HEPEVT_Interface::set().

{ return mothers;   }

vector<PhotosParticle *> Photospp::PhotosBranch::getMothers ( )

inline

Get list of mothers

Definition at line 33 of file 3.52/src/photosCInterfaces/PhotosBranch.h.

References mothers.

{ return mothers;   }

vector<PhotosParticle *> Photospp::PhotosBranch::getParticles

(

)

Get list of all particles used by branch

vector< PhotosParticle * > PhotosBranch::getParticles

(

)

Get list of all particles used by branch

Definition at line 77 of file PhotosBranch.cxx.

References daughters, mothers, and particle.

{
        vector<PhotosParticle *> ret = mothers;
        if(particle) ret.push_back(particle);
        ret.insert(ret.end(),daughters.begin(),daughters.end());
        return ret;
}

int Photospp::PhotosBranch::getSuppressionStatus ( )

inline

Check if branch is suppressed

Definition at line 42 of file 3.52/src/photosCInterfaces/PhotosBranch.h.

References suppression.

{ return suppression; }

int Photospp::PhotosBranch::getSuppressionStatus ( )

inline

Check if branch is suppressed

Definition at line 42 of file PhotosBranch.h.

References suppression.

Referenced by createBranches(), and Photospp::Photos::processParticle().

{ return suppression; }

void PhotosBranch::process

(

)

Process single branch

Definition at line 61 of file PhotosBranch.cxx.

References checkMomentumConservation(), Photospp::PH_HEPEVT_Interface::complete(), Photospp::Log::Debug(), Photospp::PH_HEPEVT_Interface::get(), Photospp::PhotosParticle::getBarcode(), mothers, particle, photos_make_c_(), Photospp::PH_HEPEVT_Interface::prepare(), and Photospp::PH_HEPEVT_Interface::set().

Referenced by Photospp::Photos::processParticle().

{
        Log::Debug(703)<<"   Processing barcode: "<<( (particle) ? particle->getBarcode() : ( (mothers.size()) ? mothers.at(0)->getBarcode() : -1) )<<endl;
        /*
        cout<<"Particles send to photos (with barcodes in brackets):"<<endl;
        vector<PhotosParticle *> get = getParticles();
        for(int i=0;i<(int)get.size();i++) cout<<"ID: "<<get.at(i)->getPdgID()<<" ("<<get.at(i)->getBarcode()<<"), "; cout<<endl;
        */
        int index = PH_HEPEVT_Interface::set(this);
        PH_HEPEVT_Interface::prepare();
        photos_make_c_(&index);
        PH_HEPEVT_Interface::complete();
        PH_HEPEVT_Interface::get();
        checkMomentumConservation();
}

void Photospp::PhotosBranch::process

(

)

Process single branch

Member Data Documentation

vector< PhotosParticle * > PhotosBranch::daughters

private

List of daughters

Definition at line 75 of file PhotosBranch.h.

Referenced by checkList(), getDaughters(), getParticles(), and PhotosBranch().

int PhotosBranch::forcing

private

State of branching forcing

Definition at line 69 of file PhotosBranch.h.

Referenced by createBranches(), getForcingStatus(), and PhotosBranch().

vector< PhotosParticle * > PhotosBranch::mothers

private

List of mothers

Definition at line 73 of file PhotosBranch.h.

Referenced by checkList(), checkMomentumConservation(), createBranches(), getMothers(), getParticles(), PhotosBranch(), and process().

PhotosParticle * PhotosBranch::particle

private

Decaying particle

Definition at line 71 of file PhotosBranch.h.

Referenced by checkList(), checkMomentumConservation(), createBranches(), getDecayingParticle(), getParticles(), PhotosBranch(), and process().

int PhotosBranch::suppression

private

State of branching suppression

Definition at line 67 of file PhotosBranch.h.

Referenced by createBranches(), getSuppressionStatus(), and PhotosBranch().

---

The documentation for this class was generated from the following files:

• PhotosBranch.h
• 3.52/src/photosCInterfaces/PhotosBranch.h
• PhotosBranch.cxx