PndEvtFilterOnInvMassCounts Class Reference

#include <PndEvtFilterOnInvMassCounts.h>

Inheritance diagram for PndEvtFilterOnInvMassCounts:

Public Types
enum	ChargeState {   kNeutral =0, kPlus, kMinus, kCharged,   kAll, kChargeLastElement }
enum	MomState { kMomTotal =0, kMomTrans, kMomZ, kMomLastElement }
enum	GeomState {   kTheta =0, kPhi, kVertexZ, kVertexRho,   kVertexRadius, kGeomLastElement }
enum	LogicOp { kAnd =0, kOr, kLogicOpLastElement }

Public Member Functions
	PndEvtFilterOnInvMassCounts ()
	PndEvtFilterOnInvMassCounts (const char *name, const char *title="PndEvtFilterOnInvMassCounts")
virtual	~PndEvtFilterOnInvMassCounts ()
Bool_t	SetPdgCodesToCombine (Int_t pdgCode1, Int_t pdgCode2, Int_t pdgCode3=kInvalidPdgCode, Int_t pdgCode4=kInvalidPdgCode, Int_t pdgCode5=kInvalidPdgCode)
Bool_t	SetMinMaxInvMass (Double_t min, Double_t max)
Bool_t	SetRhoMassParticleSelector (const char *name, Double_t cv, Double_t w, const char *type)
Bool_t	SetMinMaxCounts (Int_t min, Int_t max)
Bool_t	SetMinCounts (Int_t min)
Bool_t	SetMaxCounts (Int_t max)
Bool_t	EventMatches (Int_t evtNr)
Bool_t	FilterActive ()
Bool_t	Init ()
Bool_t	FillList (RhoCandList &rhoOutList, Int_t inPdgCode, Double_t pdgCodeCharge=kNoChargeSpecified)
Bool_t	SetListOfParticles (TClonesArray *ParticleList)
void	PrintAllTParticleInEvent ()
void	SetVerbose (Int_t verbose=12)
void	ShowAcceptedEventNumbers ()
void	ShowEvtNrsToAccept ()
void	SetTestMode (Int_t *arrayPtr, Int_t nEntries)
Bool_t	TestPassed ()
Bool_t	GetCharge (Int_t inPdgCode, Double_t *pdgCodeCharge)

Protected Attributes
TDatabasePDG *	fdbPdg
TClonesArray *	fParticleList
Int_t	fVerbose
Bool_t	fTestMode
std::set< Int_t >	fAcceptedEventNumbers
std::set< Int_t >	fEvtNrsToAccept
Int_t	fEventNr

Static Protected Attributes
static const Double_t	kNoChargeSpecified = -999.9

Private Member Functions
	ClassDef (PndEvtFilterOnInvMassCounts, 1)

Private Attributes
std::pair< Int_t, Int_t >	fCountsMinMax
RhoMassParticleSelector *	fInvMassSel
std::vector< std::pair< Int_t, Double_t > >	fPdgCodesCharges
Bool_t	fInvMassRangeSet
Bool_t	fPgdCodesSet
Bool_t	fCountRangeSet

Static Private Attributes
static const Int_t	kInvalidPdgCode = 0

Detailed Description

PndEvtFilterOnInvMassCounts.h

Author: Martin Galuska Martin dot J dot Galuska at physik dot uni minus giessen dot de

Purpose: This class is used to filter events right after event generation and before transport through the detector model events will either be accepted (and further simulated) or rejected completely.

This filter class implements such decisions based on multiplicities of invariant mass combinations within a certain mass region

The user has to specify the minimum or maximum multiplicities for the invariant mass combinations and which pdgCodes the particles should correspond to.

only particles which were produced by the event generators are taken into account for the combinations and which have the same charge as the charges corresponding to the input pdgCoded The mass hypotheses of matching particles will be set to the input pdgCodes given by the user

Definition at line 44 of file PndEvtFilterOnInvMassCounts.h.

Member Enumeration Documentation

enum FairEvtFilter::ChargeState

inherited

Enumerator
kNeutral
kPlus
kMinus
kCharged
kAll
kChargeLastElement

Definition at line 39 of file FairEvtFilter.h.

{ kNeutral=0, kPlus, kMinus, kCharged, kAll, kChargeLastElement };

enum FairEvtFilter::GeomState

inherited

Enumerator
kTheta
kPhi
kVertexZ
kVertexRho
kVertexRadius
kGeomLastElement

Definition at line 41 of file FairEvtFilter.h.

{kTheta=0, kPhi, kVertexZ, kVertexRho, kVertexRadius, kGeomLastElement};

enum FairEvtFilter::LogicOp

inherited

Enumerator
kAnd
kOr
kLogicOpLastElement

Definition at line 42 of file FairEvtFilter.h.

{ kAnd=0, kOr, kLogicOpLastElement };

enum FairEvtFilter::MomState

inherited

Enumerator
kMomTotal
kMomTrans
kMomZ
kMomLastElement

Definition at line 40 of file FairEvtFilter.h.

{ kMomTotal=0, kMomTrans, kMomZ, kMomLastElement};

Constructor & Destructor Documentation

PndEvtFilterOnInvMassCounts::PndEvtFilterOnInvMassCounts

(

)

Definition at line 12 of file PndEvtFilterOnInvMassCounts.cxx.

References fInvMassSel.

                                                         : PndEvtFilter(), fInvMassRangeSet(0), fPgdCodesSet(0), fCountRangeSet(0)

{
        fInvMassSel = 0;
}

PndEvtFilterOnInvMassCounts::PndEvtFilterOnInvMassCounts	(	const char *	name,
		const char *	title = "PndEvtFilterOnInvMassCounts"
	)

Definition at line 20 of file PndEvtFilterOnInvMassCounts.cxx.

References fInvMassSel.

                                                                                            : PndEvtFilter(name, title), fInvMassRangeSet(0), fPgdCodesSet(0), fCountRangeSet(0)
{
        fInvMassSel = 0;
}

PndEvtFilterOnInvMassCounts::~PndEvtFilterOnInvMassCounts ( )

virtual

Definition at line 26 of file PndEvtFilterOnInvMassCounts.cxx.

{}

Member Function Documentation

PndEvtFilterOnInvMassCounts::ClassDef ( PndEvtFilterOnInvMassCounts  , 1    )

private

Bool_t PndEvtFilterOnInvMassCounts::EventMatches ( Int_t  evtNr )

virtual

Implements PndEvtFilter.

Definition at line 169 of file PndEvtFilterOnInvMassCounts.cxx.

References RhoCandList::Combine(), FairEvtFilter::fAcceptedEventNumbers, fCountsMinMax, FairEvtFilter::fdbPdg, PndEvtFilter::FillList(), fInvMassSel, FairEvtFilter::fParticleList, fPdgCodesCharges, FairEvtFilter::fVerbose, RhoCandList::GetLength(), i, p1, p2, FairEvtFilter::PrintAllTParticleInEvent(), and RhoCandList::Select().

{
        if (fVerbose > 3){
                std::cout << "\n\n";
                std::cout << "Internal eventNr:" << evtNr << "\n";
                std::cout << "Nr. of simulated particles:" << fParticleList->GetEntries()<<"\n";
                PrintAllTParticleInEvent();
                std::cout << "PndEvtFilterOnInvMassCounts " << this->GetTitle() << ": " << this->GetName() << ": Beginning of EventMatches\n";
        }

        // sanity checks
        if (0==fParticleList){
                std::cout << "\n\n\n PndEvtFilterOnInvMassCounts " << this->GetTitle() << ": " << this->GetName() << ": FATAL ERROR: No particle list! Discard this event.\n\n\n";
                return kFALSE;
        }


        // do not accept events that do not have enough particles
        if ( fParticleList->GetEntriesFast() < (int)fPdgCodesCharges.size() ){
                if (0<fVerbose) std::cout << "\n\n\n PndEvtFilterOnInvMassCounts: Event contains less than " << fPdgCodesCharges.size() << " particles. " << this->GetTitle() << ": " << this->GetName() << " will not accept this event.\n\n\n";
                return kFALSE;
        }
        
        // get rho cand lists of particles and combine them
        RhoCandList p0, p1, p2, p3, p4;
        // also care about charged conjugate 
        RhoCandList ap0, ap1, ap2, ap3, ap4;
        
        std::vector<int> pdg, apdg;
        
        // make list with anti-pdg codes
        for (UInt_t i=0;i<fPdgCodesCharges.size();++i)
        {
                int cpdg = fPdgCodesCharges[i].first;
                pdg.push_back(cpdg);
                
                if (fdbPdg->GetParticle(cpdg)!=0x0 && fdbPdg->GetParticle(cpdg)->AntiParticle()!=0x0) 
                        apdg.push_back(fdbPdg->GetParticle(cpdg)->AntiParticle()->PdgCode()); 
                else 
                        apdg.push_back(cpdg);
        }
        
        // is the cc'd list the same as the original one (= a permutation)?
        bool isperm = is_permutation(pdg.begin(), pdg.end(), apdg.begin());
        
                
        FillList( p0, fPdgCodesCharges[0].first, fPdgCodesCharges[0].second );
        FillList( p1, fPdgCodesCharges[1].first, fPdgCodesCharges[1].second );
        
        // fill anti-particle lists
        if (!isperm)
        {
                FillList( ap0, apdg[0], -fPdgCodesCharges[0].second );
                FillList( ap1, apdg[1], -fPdgCodesCharges[1].second );
        }

        RhoCandList combinedList, combinedAntiList;
        
        switch( fPdgCodesCharges.size() ) {
        case 2:
                combinedList.Combine(p0, p1);
                if (!isperm) combinedAntiList.Combine(ap0, ap1);
                break;
        case 3:
                FillList( p2, fPdgCodesCharges[2].first, fPdgCodesCharges[2].second );
                combinedList.Combine(p0, p1, p2);
                if (!isperm)
                {
                        FillList( ap2, apdg[2], -fPdgCodesCharges[2].second );
                        combinedAntiList.Combine(ap0, ap1, ap2);
                }
                break;
        case 4:
                FillList( p2, fPdgCodesCharges[2].first, fPdgCodesCharges[2].second );
                FillList( p3, fPdgCodesCharges[3].first, fPdgCodesCharges[3].second );
                combinedList.Combine(p0, p1, p2, p3);
                if (!isperm)
                {
                        FillList( ap2, apdg[2], -fPdgCodesCharges[2].second );
                        FillList( ap3, apdg[3], -fPdgCodesCharges[3].second );
                        combinedAntiList.Combine(ap0, ap1, ap2, ap3);
                }
                break;
        case 5:
                FillList( p2, fPdgCodesCharges[2].first, fPdgCodesCharges[2].second );
                FillList( p3, fPdgCodesCharges[3].first, fPdgCodesCharges[3].second );
                FillList( p4, fPdgCodesCharges[4].first, fPdgCodesCharges[4].second );
                combinedList.Combine(p0, p1, p2, p3, p4);
                if (!isperm)
                {
                        FillList( p2, apdg[2], -fPdgCodesCharges[2].second );
                        FillList( p3, apdg[3], -fPdgCodesCharges[3].second );
                        FillList( p4, apdg[4], -fPdgCodesCharges[4].second );
                        combinedAntiList.Combine(ap0, ap1, ap2, ap3, ap4);
                }
                break;
        default:
                std::cout << "FATAL ERROR in PndEvtFilterOnInvMassCounts::EventMatches of " << this->GetTitle() << ": " << this->GetName() << ". The number of pdg codes is not supported. \n";
                return kFALSE;
                break;
        }


        if (fVerbose > 3){
                std::cout << "Looking for inv. masses of " << fPdgCodesCharges.size() << " particles.\n";
                std::cout << p0.GetLength() << " / " << p1.GetLength() << " / " << p2.GetLength() << " / " << p3.GetLength() << " / " << p4.GetLength();
                std::cout << " particles in initial lists\n";
                std::cout << combinedList.GetLength() << " particles in combined list before mass selection\n";
        }

        // select acceptable invariant mass range
        combinedList.Select(fInvMassSel);
        if (!isperm) combinedAntiList.Select(fInvMassSel);
        
        int ncand = combinedList.GetLength() + combinedAntiList.GetLength();
        
        if (fVerbose > 3){
                std::cout << combinedList.GetLength()<<" + "<< combinedAntiList.GetLength()<< " particles in combined list after mass selection\n";
        }

        // check if number of acceptable combinations is within acceptable limits
        if ( ncand < fCountsMinMax.first){
                if (fVerbose >9){
                        std::cout << "Event is not accepted by " << this->GetTitle() << ": " << this->GetName() << " because there are too few matching inv. mass combinations in the event. \n";
                }
                return kFALSE;
        }
        if ( ncand > fCountsMinMax.second ){
                if (fVerbose >9){
                        std::cout << "Event is not accepted by " << this->GetTitle() << ": " << this->GetName() << " because there are too many matching inv. mass combinations in the event. \n";
                }
                return kFALSE;
        }


        // for QA
        fAcceptedEventNumbers.insert(evtNr);
        if (fVerbose >5){
                std::cout << "\n Event matches " << this->GetTitle() << ": " << this->GetName() << "\n\n";
        }
        return kTRUE;
}

Bool_t PndEvtFilter::FillList ( RhoCandList &  rhoOutList, Int_t  inPdgCode, Double_t  pdgCodeCharge = kNoChargeSpecified  )

inherited

Definition at line 31 of file PndEvtFilter.cxx.

References RhoCandList::Add(), RhoCandList::Cleanup(), Double_t, FairEvtFilter::fdbPdg, FairEvtFilter::fParticleList, FairEvtFilter::GetCharge(), FairEvtFilter::kNoChargeSpecified, particle, RhoCandidate::SetMcTruth(), RhoCandidate::SetPos(), and RhoCandidate::SetType().

Referenced by EventMatches().

{
        rhoOutList.Cleanup();


        // Get charge corresponding to pdgCode if it has not been passed as an argument by the user
        if ( kNoChargeSpecified == pdgCodeCharge ){
                if ( kFALSE == GetCharge( inPdgCode, &pdgCodeCharge ) ) {
                        return kFALSE;
                }
        }



        // find all TParticle on stack with a charge(particle) == charge(pdgCode)
        for (Int_t iPart=0; iPart<fParticleList->GetEntries(); ++iPart) {
                TParticle *particle = (TParticle*)fParticleList->At(iPart);

                // get charge for particle
                TParticlePDG*  pdt = particle->GetPDG();
                if (0==pdt) continue; // unknown particle type (KG, 10/2015)

                Double_t pCharge = pdt->Charge()/3.; // TParticlePDG contains charge in units of |e|/3

                if ( pdgCodeCharge != pCharge ){ continue; } // skip all particles with different charge


                //              // mark stable particles
                //              Bool_t isStable = kFALSE;
                //              switch(abs(particle->GetPdgCode())) {
                //              case 22:   isStable = true; break;
                //              case 11:   isStable = true; break;
                //              case 13:   isStable = true; break;
                //              case 211:  isStable = true; break;
                //              case 321:  isStable = true; break;
                //              case 2212: isStable = true; break;
                //              }
                //
                //              if (kFALSE == isStable) { continue; } // skip all particles which are not considered stable


                // store TParticle as RhoCandidate in output TCA and modify 4 momentum to match particle hypothesis
                TVector3 pVertex(particle->Vx(),particle->Vy(),particle->Vz());
                TLorentzVector p4(particle->Px(),particle->Py(),particle->Pz(),particle->Energy());
                // assign mass hypothesis according to pdgCode
                p4.SetVectM(p4.Vect(),fdbPdg->GetParticle(inPdgCode)->Mass());
                // construct RhoCandidate
                RhoCandidate rhoCand(p4,pCharge);
                rhoCand.SetMcTruth(&rhoCand);
                rhoCand.SetPos(pVertex);
                rhoCand.SetType(inPdgCode);
                // save to output
                rhoOutList.Add(&rhoCand);
        }//trackloop


        return kTRUE;
}

Bool_t PndEvtFilterOnInvMassCounts::FilterActive ( )

inlinevirtual

Implements PndEvtFilter.

Definition at line 80 of file PndEvtFilterOnInvMassCounts.h.

References fCountRangeSet, fInvMassRangeSet, and fPgdCodesSet.

{ return ( fInvMassRangeSet && fPgdCodesSet && fCountRangeSet ); };

Bool_t FairEvtFilter::GetCharge ( Int_t  inPdgCode, Double_t *  pdgCodeCharge  )

inherited

Definition at line 62 of file FairEvtFilter.cxx.

References FairEvtFilter::fdbPdg, FairEvtFilter::fVerbose, and FairEvtFilter::kNoChargeSpecified.

Referenced by FairEvtFilterOnSingleParticleCounts::CountCharge(), PndEvtFilter::FillList(), and SetPdgCodesToCombine().

{
        // Try to find the pdg code
        TParticlePDG *ptrToPdg = fdbPdg->GetParticle(inPdgCode);
        if ( 0 == ptrToPdg) {
                // ignore particles with unknown charges
                std::cout << "WARNING from FairEvtFilter::GetCharge Charge of pdgCode " << inPdgCode << " is unknown and will be ignored!\n";
                *pdgCodeCharge = kNoChargeSpecified;
                return kFALSE;
        }
        *pdgCodeCharge = ptrToPdg->Charge()/3.; // TParticlePDG contains charge in units of |e|/3
        if ( fVerbose > 1 ) std::cout << "Found pdgCodeCharge = " << *pdgCodeCharge << " for inPdgCode " << inPdgCode << '\n';
        return kTRUE;
}

Bool_t PndEvtFilter::Init ( )

inlineinherited

Definition at line 54 of file PndEvtFilter.h.

{ return kTRUE;}

void FairEvtFilter::PrintAllTParticleInEvent ( )

inherited

Definition at line 52 of file FairEvtFilter.cxx.

References FairEvtFilter::fParticleList, and particle.

Referenced by EventMatches(), and FairEvtFilterOnSingleParticleCounts::EventMatches().

                                            {
        // sanity checks
        if (0==fParticleList){ return; }
        for (Int_t iPart=0; iPart<fParticleList->GetEntries(); ++iPart) {
                TParticle *particle = (TParticle*)fParticleList->At(iPart);
                particle->Print();
        }
}

Bool_t FairEvtFilter::SetListOfParticles ( TClonesArray *  ParticleList )

inlineinherited

Definition at line 57 of file FairEvtFilter.h.

References FairEvtFilter::fParticleList.

Referenced by FairFilteredPrimaryGenerator::GenerateEvent().

{fParticleList=ParticleList; return kTRUE;};

Bool_t PndEvtFilterOnInvMassCounts::SetMaxCounts ( Int_t  max )

inline

Definition at line 69 of file PndEvtFilterOnInvMassCounts.h.

References SetMinMaxCounts().

{ return SetMinMaxCounts(0, max); };

Bool_t PndEvtFilterOnInvMassCounts::SetMinCounts ( Int_t  min )

inline

Definition at line 68 of file PndEvtFilterOnInvMassCounts.h.

References SetMinMaxCounts().

{ return SetMinMaxCounts(min, 9999); };

Bool_t PndEvtFilterOnInvMassCounts::SetMinMaxCounts	(	Int_t	min,
		Int_t	max
	)

Definition at line 29 of file PndEvtFilterOnInvMassCounts.cxx.

References fCountRangeSet, fCountsMinMax, max(), and min().

Referenced by prod_fsim(), quickfsimana(), SetMaxCounts(), SetMinCounts(), and sim_filter_inv_mass().

{
        if ( kTRUE == fCountRangeSet ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts of " << this->GetTitle() << ": " << this->GetName() << ": \n";
                std::cout << "You are trying to set min and max counts for the events more than once! \n";
                std::cout << "I take the first setting of min = " << fCountsMinMax.first << " max = " << fCountsMinMax.second << " and ignore the later ones! \n\n\n";
                std::cout << "Your current input of min = " << min << " max = " << max << " is ignored! \n\n\n";
                return kFALSE;
        }

        // check if filter can be applied (make sure that min and max make sense)
        if ( min < 0 ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts of " << this->GetTitle() << ": " << this->GetName() << ": Filter could not be added. You are trying to request that your events should have at least a number <= 0 of some particles. That makes no sense. \n\n\n";
                return kFALSE;
        }
        if ( max < min ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts  of " << this->GetTitle() << ": " << this->GetName() << ": Filter could not be added. You are trying to request that your events should have at most a number of some particles which is less than the minimum number that you request. That makes no sense. \n\n\n";
                return kFALSE;
        }

        // set the minimum and maximum number of particle combinations
        fCountsMinMax = std::pair<Int_t, Int_t> (min,max);

        // turn filter on
        fCountRangeSet = kTRUE;

        return kTRUE;
}

Bool_t PndEvtFilterOnInvMassCounts::SetMinMaxInvMass	(	Double_t	min,
		Double_t	max
	)

Definition at line 59 of file PndEvtFilterOnInvMassCounts.cxx.

References Double_t, fInvMassRangeSet, fInvMassSel, max(), and min().

Referenced by prod_fsim(), quickfsimana(), and sim_filter_inv_mass().

{
        if ( kTRUE == fInvMassRangeSet ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts of " << this->GetTitle() << ": " << this->GetName() << ": \n";
                std::cout << "You are trying to set the limits for the invariant mass selection set more than once! \n";
                std::cout << "I take the first setting and ignore the later ones! \n\n\n";
                return kFALSE;
        }

        // check if filter can be applied (make sure that min and max make sense)
        if ( min < 0. ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts " << this->GetTitle() << ": " << this->GetName() << ": Filter could not be added. Min inv. mass too low. \n\n\n";
                return kFALSE;
        }
        if ( max < min ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts " << this->GetTitle() << ": " << this->GetName() << ": Filter could not be added. Max inv. mass lower than min inv. mass. \n\n\n";
                return kFALSE;
        }

        // set up the inv. mass selector
        Double_t mu = (min + max)/2.0;
        Double_t fullWidth = max - min;

        fInvMassSel = new RhoMassParticleSelector("invMassFilter", mu, fullWidth);

        // turn filter on
        fInvMassRangeSet = kTRUE;

        return kTRUE;
}

Bool_t PndEvtFilterOnInvMassCounts::SetPdgCodesToCombine	(	Int_t	pdgCode1,
		Int_t	pdgCode2,
		Int_t	pdgCode3 = kInvalidPdgCode,
		Int_t	pdgCode4 = kInvalidPdgCode,
		Int_t	pdgCode5 = kInvalidPdgCode
	)

Definition at line 119 of file PndEvtFilterOnInvMassCounts.cxx.

References Double_t, fPdgCodesCharges, fPgdCodesSet, FairEvtFilter::GetCharge(), kInvalidPdgCode, and FairEvtFilter::kNoChargeSpecified.

Referenced by prod_fsim(), quickfsimana(), and sim_filter_inv_mass().

                                                                                                                                        {
        if ( kTRUE == fPgdCodesSet ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts of " << this->GetTitle() << ": " << this->GetName() << ": \n";
                std::cout << "You are trying to set the pdgCodes to combine into invariant mass combinations more than once! \n";
                std::cout << "I take the first setting and ignore the later ones! \n\n\n";
                return kFALSE;
        }


        // create temp. vector of all input pdg codes
        std::vector< Int_t > pdgCodes, apdgCodes;
        pdgCodes.push_back(pdgCode1);
        pdgCodes.push_back(pdgCode2);
        pdgCodes.push_back(pdgCode3);
        pdgCodes.push_back(pdgCode4);
        pdgCodes.push_back(pdgCode5);
        
        
//      std::cout <<"[PndEvtFilterOnInvMassCounts]"<<pdgCode1<<" "<<pdgCode2<<" "<<pdgCode3<<" "<<pdgCode4<<" "<<pdgCode5<<std::endl;
        
        // Copy valid entries from pdgCodes to fPdgCodesCharges
        for (UInt_t iPdgCodes = 0; iPdgCodes < pdgCodes.size(); ++iPdgCodes){
                // skip kInvalidPdgCode entries in pdgCodes
                if ( kInvalidPdgCode == pdgCodes[iPdgCodes] ) { continue; }
                // If we have invalid pdg codes warn user
                Double_t pdgCodeCharge = kNoChargeSpecified;
                if ( kFALSE == GetCharge( pdgCodes[iPdgCodes], &pdgCodeCharge ) ) {
                        std::cout << "\n\n\nFATAL ERROR FROM PndEvtFilterOnInvMassCounts::SetPdgCodesToCombine of " << this->GetTitle() << ": " << this->GetName() << " Charge for pdg code " << pdgCodes[iPdgCodes] << " was not found. \n";
                        std::cout << "Filter was not set. \n\n\n”";
                        fPgdCodesSet = kFALSE;
                        return kFALSE;
                }
                // If we get to this point, we have successfully added a pdg code
                fPdgCodesCharges.push_back(std::pair<Int_t, Double_t> (pdgCodes[iPdgCodes], pdgCodeCharge));
        }

        // Check that we had at least two correctly set pdg codes
        if ( 1 < fPdgCodesCharges.size() ){
                fPgdCodesSet = kTRUE;
                return kTRUE;
        } else {
                std::cout << "\n\n\nFATAL ERROR FROM PndEvtFilterOnInvMassCounts::SetPdgCodesToCombine of " << this->GetTitle() << ": " << this->GetName() << " You have to specify at least 2 valid pdg codes. \n";
                std::cout << "Filter was not set. \n\n\n”";
                fPgdCodesSet = kFALSE;
                return kFALSE;
        }
}

Bool_t PndEvtFilterOnInvMassCounts::SetRhoMassParticleSelector	(	const char *	name,
		Double_t	cv,
		Double_t	w,
		const char *	type
	)

Definition at line 91 of file PndEvtFilterOnInvMassCounts.cxx.

References fInvMassRangeSet, and fInvMassSel.

{
        if ( kTRUE == fInvMassRangeSet ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts of " << this->GetTitle() << ": " << this->GetName() << ": \n";
                std::cout << "You are trying to set the limits for the invariant mass selection set more than once! \n";
                std::cout << "I take the first setting and ignore the later ones! \n\n\n";
                return kFALSE;
        }

        // check if filter can be applied (make sure that min and max make sense)
        if ( cv <= 0. ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts " << this->GetTitle() << ": " << this->GetName() << ": Filter could not be added. cv has to be positive. \n\n\n";
                return kFALSE;
        }
        if ( w < 0. ){
                std::cout << "\n\n\n  -WARNING from PndEvtFilterOnInvMassCounts " << this->GetTitle() << ": " << this->GetName() << ": Filter could not be added. w has to be non-negative. \n\n\n";
                return kFALSE;
        }

        fInvMassSel = new RhoMassParticleSelector(name, cv, w, type);

        // turn filter on
        fInvMassRangeSet = kTRUE;

        return kTRUE;
}

void FairEvtFilter::SetTestMode ( Int_t *  arrayPtr, Int_t  nEntries  )

inlineinherited

Definition at line 91 of file FairEvtFilter.h.

References FairEvtFilter::fEvtNrsToAccept, and FairEvtFilter::fTestMode.

                                                         {
                //turns on the test mode with the declared fEvtNrsToAccept
                fTestMode=kTRUE;
                std::set<Int_t> evtNrsToAccept (arrayPtr,arrayPtr+nEntries);
                fEvtNrsToAccept=evtNrsToAccept;
        }

void FairEvtFilter::SetVerbose ( Int_t  verbose = 12 )

inlineinherited

Definition at line 67 of file FairEvtFilter.h.

References FairEvtFilter::fVerbose, and verbose.

Referenced by sim_filter_inv_mass().

                                         {
                if ( verbose >= 0 ){
                        fVerbose = verbose;
                        std::cout << "FairEvtFilter: fVerbose is now set to " << verbose << "\n";
                } else {
                        std::cout << "\n\n\n  -WARNING from FairEvtFilter: verbose must be a positive number! Check your SetVerbose call!\n\n\n";
                }
        }

void FairEvtFilter::ShowAcceptedEventNumbers ( )

inlineinherited

Definition at line 77 of file FairEvtFilter.h.

References FairEvtFilter::fAcceptedEventNumbers.

                                       {
                // for QA
                // shows fAcceptedEventNumbers that is filled after running a simulation
                std::cout << "\n fAcceptedEventNumbers" << " = " << fAcceptedEventNumbers;
        }

void FairEvtFilter::ShowEvtNrsToAccept ( )

inlineinherited

Definition at line 84 of file FairEvtFilter.h.

References FairEvtFilter::fEvtNrsToAccept.

                                 {
                // for QA
                //shows fEvtNrsToAccept that has to be set if you want to run your simulation in test mode
                std::cout << "\n fEvtNrsToAccept" << " = " << fEvtNrsToAccept;
        }

Bool_t FairEvtFilter::TestPassed ( )

inlineinherited

Definition at line 99 of file FairEvtFilter.h.

References FairEvtFilter::fAcceptedEventNumbers, FairEvtFilter::fEvtNrsToAccept, and FairEvtFilter::fTestMode.

                           {
                if(kFALSE==fTestMode){
                        std::cout << "\n\n\n WARNING from FairEvtFilter: Test mode not set.\n\n\n";
                        return kFALSE;
                }
                if(fAcceptedEventNumbers==fEvtNrsToAccept){
                        //std::cout << "\n\n\n FairEvtFilter: Test passed.\n\n\n";
                        return kTRUE;
                }else{
                        //std::cout << "\n\n\n FairEvtFilter: Test failed. Check your SetTestMode calls. \n\n\n";
                        return kFALSE;
                }
        }

Member Data Documentation

std::set<Int_t> FairEvtFilter::fAcceptedEventNumbers

protectedinherited

Definition at line 134 of file FairEvtFilter.h.

Referenced by EventMatches(), FairEvtFilterOnSingleParticleCounts::EventMatches(), FairEvtFilter::FairEvtFilter(), FairEvtFilter::ShowAcceptedEventNumbers(), and FairEvtFilter::TestPassed().

Bool_t PndEvtFilterOnInvMassCounts::fCountRangeSet

private

Definition at line 95 of file PndEvtFilterOnInvMassCounts.h.

Referenced by FilterActive(), and SetMinMaxCounts().

std::pair<Int_t,Int_t> PndEvtFilterOnInvMassCounts::fCountsMinMax

private

Definition at line 80 of file PndEvtFilterOnInvMassCounts.h.

Referenced by EventMatches(), and SetMinMaxCounts().

TDatabasePDG* FairEvtFilter::fdbPdg

protectedinherited

Definition at line 126 of file FairEvtFilter.h.

Referenced by EventMatches(), FairEvtFilter::FairEvtFilter(), PndEvtFilter::FillList(), and FairEvtFilter::GetCharge().

Int_t FairEvtFilter::fEventNr

protectedinherited

Definition at line 136 of file FairEvtFilter.h.

std::set<Int_t> FairEvtFilter::fEvtNrsToAccept

protectedinherited

Definition at line 135 of file FairEvtFilter.h.

Referenced by FairEvtFilter::SetTestMode(), FairEvtFilter::ShowEvtNrsToAccept(), and FairEvtFilter::TestPassed().

Bool_t PndEvtFilterOnInvMassCounts::fInvMassRangeSet

private

Definition at line 93 of file PndEvtFilterOnInvMassCounts.h.

Referenced by FilterActive(), SetMinMaxInvMass(), and SetRhoMassParticleSelector().

RhoMassParticleSelector* PndEvtFilterOnInvMassCounts::fInvMassSel

private

Definition at line 87 of file PndEvtFilterOnInvMassCounts.h.

Referenced by EventMatches(), PndEvtFilterOnInvMassCounts(), SetMinMaxInvMass(), and SetRhoMassParticleSelector().

TClonesArray* FairEvtFilter::fParticleList

protectedinherited

Definition at line 131 of file FairEvtFilter.h.

Referenced by EventMatches(), FairEvtFilterOnSingleParticleCounts::EventMatches(), FairEvtFilter::FairEvtFilter(), PndEvtFilter::FillList(), FairEvtFilter::PrintAllTParticleInEvent(), and FairEvtFilter::SetListOfParticles().

std::vector< std::pair< Int_t, Double_t > > PndEvtFilterOnInvMassCounts::fPdgCodesCharges

private

Definition at line 89 of file PndEvtFilterOnInvMassCounts.h.

Referenced by EventMatches(), and SetPdgCodesToCombine().

Bool_t PndEvtFilterOnInvMassCounts::fPgdCodesSet

private

Definition at line 94 of file PndEvtFilterOnInvMassCounts.h.

Referenced by FilterActive(), and SetPdgCodesToCombine().

Bool_t FairEvtFilter::fTestMode

protectedinherited

Definition at line 133 of file FairEvtFilter.h.

Referenced by FairEvtFilter::SetTestMode(), and FairEvtFilter::TestPassed().

Int_t FairEvtFilter::fVerbose

protectedinherited

Definition at line 132 of file FairEvtFilter.h.

Referenced by FairEvtFilterOnSingleParticleCounts::AcceptChargeCounter(), FairEvtFilterOnSingleParticleCounts::AcceptGeometry(), FairEvtFilterOnSingleParticleCounts::AcceptMomentum(), FairEvtFilterOnSingleParticleCounts::AcceptPdgCounter(), FairEvtFilterOnSingleParticleCounts::AndMinMaxCharge(), FairEvtFilterOnSingleParticleCounts::AndMinMaxGeom(), FairEvtFilterOnSingleParticleCounts::AndMinMaxMom(), FairEvtFilterOnSingleParticleCounts::AndMinMaxPdgCodes(), FairEvtFilterOnSingleParticleCounts::CountCharge(), FairEvtFilterOnSingleParticleCounts::CountPdg(), EventMatches(), FairEvtFilterOnSingleParticleCounts::EventMatches(), FairEvtFilter::GetCharge(), FairEvtFilterOnSingleParticleCounts::InitCounters(), and FairEvtFilter::SetVerbose().

const Int_t PndEvtFilterOnInvMassCounts::kInvalidPdgCode = 0

staticprivate

Definition at line 100 of file PndEvtFilterOnInvMassCounts.h.

Referenced by SetPdgCodesToCombine().

const Double_t FairEvtFilter::kNoChargeSpecified = -999.9

staticprotectedinherited

Definition at line 129 of file FairEvtFilter.h.

Referenced by PndEvtFilter::FillList(), FairEvtFilter::GetCharge(), and SetPdgCodesToCombine().

---

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

• PndEvtFilterOnInvMassCounts.h
• PndEvtFilterOnInvMassCounts.cxx