#include <RandomtXiGunProducer.h>

Inheritance diagram for edm::RandomtXiGunProducer:

Public Member Functions
	RandomtXiGunProducer (const ParameterSet &)

	~RandomtXiGunProducer () override

Public Member Functions inherited from edm::BaseRandomtXiGunProducer
	BaseRandomtXiGunProducer (const ParameterSet &)

void	beginRun (const edm::Run &r, const edm::EventSetup &) override

void	endRun (const edm::Run &r, const edm::EventSetup &) override

void	endRunProduce (edm::Run &r, const edm::EventSetup &) override

	~BaseRandomtXiGunProducer () override

Public Member Functions inherited from edm::one::EDProducer< one::WatchRuns, EndRunProducer >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Protected Attributes
double	fMaxt

double	fMaxXi

double	fMint

double	fMinXi

Protected Attributes inherited from edm::BaseRandomtXiGunProducer
double	fECMS

HepMC::GenEvent *	fEvt

bool	fFireBackward

bool	fFireForward

bool	fLog_t

double	fMaxPhi

double	fMinPhi

std::vector< int >	fPartIDs

ESHandle< HepPDT::ParticleDataTable >	fPDGTable

double	fpEnergy

int	fVerbosity

const HepPDT::ParticleData *	PData

Private Member Functions
HepMC::FourVector	make_particle (double t, double Xi, double phi, int PartID, int direction)

double	Minimum_t (double xi)

void	produce (Event &e, const EventSetup &es) override

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 10 of file RandomtXiGunProducer.h.

Constructor & Destructor Documentation

RandomtXiGunProducer::RandomtXiGunProducer ( const ParameterSet & pset )

Definition at line 26 of file RandomtXiGunProducer.cc.

References edm::BaseRandomtXiGunProducer::fLog_t, fMaxt, fMaxXi, fMint, fMinXi, edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

                                                                       : 
    BaseRandomtXiGunProducer(pset)
 {
 
    ParameterSet defpset ;
    edm::ParameterSet pgun_params = 
       pset.getParameter<edm::ParameterSet>("PGunParameters") ;
   
    fMint = pgun_params.getParameter<double>("Mint");
    fMaxt = pgun_params.getParameter<double>("Maxt");
    fMinXi= pgun_params.getParameter<double>("MinXi");
    fMaxXi= pgun_params.getParameter<double>("MaxXi");
    fLog_t= pgun_params.getUntrackedParameter<bool>("Log_t",false);
   
   produces<HepMCProduct>("unsmeared");
   produces<GenEventInfoProduct>();
 }

RandomtXiGunProducer::~RandomtXiGunProducer ( )

override

Definition at line 44 of file RandomtXiGunProducer.cc.

 {
    // no need to cleanup GenEvent memory - done in HepMCProduct
 }

Member Function Documentation

HepMC::FourVector RandomtXiGunProducer::make_particle	(	double	t,
		double	Xi,
		double	phi,
		int	PartID,
		int	direction
	)

private

Definition at line 152 of file RandomtXiGunProducer.cc.

References funct::cos(), edm::BaseRandomtXiGunProducer::fpEnergy, edm::BaseRandomtXiGunProducer::fVerbosity, ResonanceBuilder::mass, edm::BaseRandomtXiGunProducer::PData, funct::sin(), mathSSE::sqrt(), and theta().

Referenced by produce().

 {
        double mass   = PData->mass().value() ;
        double fpMom  = sqrt(fpEnergy*fpEnergy-mass*mass); // momentum of beam proton
        double sEnergy = (1.0-Xi)*fpEnergy; // energy of scattered particle
        double sMom    = sqrt(sEnergy*sEnergy-mass*mass); // momentum of scattered particle
        double min_t = -2.*(fpMom*sMom-fpEnergy*sEnergy+mass*mass);
        if (t<min_t) t=min_t; // protect against kinemactically forbiden region
        long double theta  = acos((-t/2.- mass*mass + fpEnergy*sEnergy)/(sMom*fpMom)); // use t = -t
 
        if (direction<1) theta = acos(-1.) - theta;
 
        double px     = sMom*cos(phi)*sin(theta);
        double py     = sMom*sin(phi)*sin(theta);
        double pz     = sMom*cos(theta) ;  // the direction is already set by the theta angle
        if (fVerbosity > 0) 
           edm::LogInfo("RandomXiGunProducer") << "-----------------------------------------------------------------------------------------------------\n"
                                               << "Produced a proton with  phi : " << phi << " theta: " << theta << " t: " << t << " Xi: " << Xi << "\n"
                                               << "                         Px : " << px  << " Py : " << py << " Pz : " << pz << "\n"
                                               << "                   direction: " << direction  << "\n"
                                               << "-----------------------------------------------------------------------------------------------------"
                                               << std::endl;
 
        return HepMC::FourVector(px,py,pz,sEnergy) ;
 }

double edm::RandomtXiGunProducer::Minimum_t ( double xi )

inlineprivate

Definition at line 21 of file RandomtXiGunProducer.h.

References edm::BaseRandomtXiGunProducer::fpEnergy, edm::BaseRandomtXiGunProducer::PData, funct::pow(), mathSSE::sqrt(), and protons_cff::xi.

Referenced by produce().

                                 {double partE = fpEnergy*(1.-xi);
                                  double massSQ= pow(PData->mass().value(),2);
                                  double partP = sqrt(partE*partE-massSQ);
                                  return -2.*(sqrt(fpEnergy*fpEnergy-massSQ)*partP-fpEnergy*partE+massSQ);
                                 };

void RandomtXiGunProducer::produce	(	edm::Event &	e,
		const EventSetup &	es
	)

overrideprivate

Definition at line 49 of file RandomtXiGunProducer.cc.

 {
    edm::Service<edm::RandomNumberGenerator> rng;
    CLHEP::HepRandomEngine* engine = &rng->getEngine(e.streamID());
 
    if ( fVerbosity > 0 )
    {
       edm::LogInfo("RandomtXiGunProducer")<< "Begin New Event Generation\n"; 
    }
    // event loop (well, another step in it...)
           
    // no need to clean up GenEvent memory - done in HepMCProduct
    // 
    
    // here re-create fEvt (memory)
    //
    fEvt = new HepMC::GenEvent() ;
    
    // now actualy, cook up the event from PDGTable and gun parameters
    //
    // 1st, primary vertex
    //
    HepMC::GenVertex* Vtx = new HepMC::GenVertex(HepMC::FourVector(0.,0.,0.));
 
    // loop over particles
    //
    int barcode = 1 ;
    for (unsigned int ip=0; ip<fPartIDs.size(); ++ip)
    {
        int PartID = fPartIDs[ip];
 //  t = -2*P*P'*(1-cos(theta)) -> t/(2*P*P')+1=cos(theta)
 // xi = 1 - P'/P  --> P'= (1-xi)*P
 // 
        PData = fPDGTable->particle(HepPDT::ParticleID(abs(PartID)));
        if (!PData) exit(1);
        double t  = 0;
        double Xi = 0;
        double phi=0;
        if (fFireForward) {
           while(true) {
                Xi     = CLHEP::RandFlat::shoot(engine,fMinXi,fMaxXi);
                double min_t = std::max(fMint,Minimum_t(Xi));
                double max_t = fMaxt;
                if (min_t>max_t) {
                   edm::LogInfo("RandomtXiGunProducer") << "WARNING: t limits redefined (unphysical values for given xi).\n";
                   max_t = min_t;
                }
                t      = (fLog_t)?pow(CLHEP::RandFlat::shoot(engine,log10(min_t),log10(max_t)),10):
                                      CLHEP::RandFlat::shoot(engine,min_t,max_t);
                break;
           }
           phi    = CLHEP::RandFlat::shoot(engine,fMinPhi, fMaxPhi) ;
           HepMC::GenParticle* Part = 
               new HepMC::GenParticle(make_particle(t,Xi,phi,PartID,1),PartID,1);
           Part->suggest_barcode( barcode ) ;
           barcode++ ;
           Vtx->add_particle_out(Part);
        }
        if ( fFireBackward) {
           while(true) {
                Xi     = CLHEP::RandFlat::shoot(engine,fMinXi,fMaxXi);
                double min_t = std::max(fMint,Minimum_t(Xi));
                double max_t = fMaxt;
                if (min_t>max_t) {
                   edm::LogInfo("RandomtXiGunProducer") << "WARNING: t limits redefined (unphysical values for given xi)." << endl;
                   max_t = min_t;
                }
                t      = (fLog_t)?pow(CLHEP::RandFlat::shoot(engine,log10(min_t),log10(max_t)),10):
                                      CLHEP::RandFlat::shoot(engine,min_t,max_t);
                break;
           }
           phi    = CLHEP::RandFlat::shoot(engine,fMinPhi, fMaxPhi) ;
       HepMC::GenParticle* Part2 =
          new HepMC::GenParticle(make_particle(t,Xi,phi,PartID,-1),PartID,1);
       Part2->suggest_barcode( barcode ) ;
       barcode++ ;
       Vtx->add_particle_out(Part2) ;
        }
    }
 
    fEvt->add_vertex(Vtx) ;
    fEvt->set_event_number(e.id().event()) ;
    fEvt->set_signal_process_id(20) ; 
         
    if ( fVerbosity > 0 )
    {
       fEvt->print() ;  
    }
 
    std::unique_ptr<HepMCProduct> BProduct(new HepMCProduct()) ;
    BProduct->addHepMCData( fEvt );
    e.put(std::move(BProduct),"unsmeared");
 
    std::unique_ptr<GenEventInfoProduct> genEventInfo(new GenEventInfoProduct(fEvt));
    e.put(std::move(genEventInfo));
     
    if ( fVerbosity > 0 )
    {
       // for testing purpose only
       // fEvt->print() ; // prints empty info after it's made into edm::Event
       edm::LogInfo("RandomtXiGunProducer") << " Event Generation Done \n";
    }
 }