#include <AddFlowByShifting.h>

Inheritance diagram for AddFlowByShifting:

Collaboration diagram for AddFlowByShifting:

Public Member Functions
	AddFlowByShifting (const std::string &name, ISvcLocator *pSvcLocator)

StatusCode	initialize ()

StatusCode	execute ()

StatusCode	finalize ()

virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...

virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static double	vn_func (double x, void *params)

static double	vn_func_derivative (double x, void *params)

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
CLHEP::HepRandomEngine *	getRandomEngine (const std::string &streamName, const EventContext &ctx) const

double	SetParentToRanPhi (HepMC::GenParticlePtr parent, CLHEP::HepRandomEngine *rndmEngine)

double	AddFlowToParent (HepMC::GenParticlePtr parent, const HijingEventParams *hijing_pars)

void	MoveDescendantsToParent (HepMC::GenParticlePtr parent, double phishift)

void	jjia_minbias_new (double b, double eta, double pt)

void	jjia_minbias_new_v2only (double b, double eta, double pt)

void	fixed_vn (double b, double eta, double pt)

void	fixed_v2 (double b, double eta, double pt)

void	jjia_minbias_old (double b, double eta, double pt)

void	ao_test (double b, double eta, double pt)

void	custom_vn (double b, double eta, double pt)

void	p_Pb_cent_eta_indep (double b, double eta, double pt)

void	Set_EbE_Fluctuation_Multipliers (HepMC::GenVertexPtr mainvtx, float b, CLHEP::HepRandomEngine *rndmEngine)

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
void(AddFlowByShifting::*	m_flow_function )(double b, double eta, double pt)

TGraph *	m_graph_fluc {}

ServiceHandle< IAthRNGSvc >	m_rndmSvc {this, "RndmSvc", "AthRNGSvc"}

StringProperty	m_inkey {this, "McTruthKey", "GEN_EVENT"}

StringProperty	m_outkey {this, "McFlowKey", "FLOW_EVENT"}

IntegerProperty	m_ranphi_sw {this, "RandomizePhi", 0}

StringProperty	m_flow_function_name {this, "FlowFunctionName", "jjia_minbias_new"}

StringProperty	m_flow_implementation {this, "FlowImplementation", "exact"}

int	m_flow_implementation_type {0}

BooleanProperty	m_flow_fluctuations {this, "FlowFluctuations", false}

IntegerProperty	m_floweta_sw {this, "FlowEtaSwitch", 0 }

FloatProperty	m_flow_maxeta {this, "FlowMaxEtaCut", 10.0}

FloatProperty	m_flow_mineta {this, "FlowMinEtaCut", 0.f}

IntegerProperty	m_flowpt_sw {this, "FlowPtSwitch", 0}

FloatProperty	m_flow_maxpt {this, "FlowMaxPtCut", 1000000.f}

FloatProperty	m_flow_minpt {this, "FlowMinPtCut", 0.f}

IntegerProperty	m_flowb_sw {this, "FlowBSwitch", 0}

FloatProperty	m_custom_v1 {this, "custom_v1", 0.f}

FloatProperty	m_custom_v2 {this, "custom_v2", 0.f}

FloatProperty	m_custom_v3 {this, "custom_v3", 0.f}

FloatProperty	m_custom_v4 {this, "custom_v4", 0.f}

FloatProperty	m_custom_v5 {this, "custom_v5", 0.f}

FloatProperty	m_custom_v6 {this, "custom_v6", 0.f}

int	m_particles_processed {0}

float	m_psi_n [6]

float	m_v_n [6]

float	m_EbE_Multiplier_vn [6]

DataObjIDColl	m_extendedExtraObjects

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 39 of file AddFlowByShifting.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ AddFlowByShifting()

AddFlowByShifting::AddFlowByShifting	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 65 of file AddFlowByShifting.cxx.

                                                                                     :
   AthAlgorithm(name, pSvcLocator)
 {
   m_flow_function= NULL;
   for(int ihar = 0; ihar< 6; ihar++){
     m_psi_n[ihar] =0.0;
     m_v_n  [ihar] =0.0;
     m_EbE_Multiplier_vn[ihar]=1.0;
   }
 }

Member Function Documentation

◆ AddFlowToParent()

double AddFlowByShifting::AddFlowToParent	(	HepMC::GenParticlePtr	parent,
		const HijingEventParams *	hijing_pars
	)

private

Definition at line 358 of file AddFlowByShifting.cxx.

 {
   m_particles_processed++;
   CLHEP::HepLorentzVector momentum(parent->momentum().px(),
                                    parent->momentum().py(),
                                    parent->momentum().pz(),
                                    parent->momentum().e());
   double pt    = parent->momentum().perp();
   double eta   = parent->momentum().pseudoRapidity();
   double phi_0 = parent->momentum().phi();
  
   int error_=0;
   if(pt    !=pt)    {ATH_MSG_ERROR("ERROR pt  of track  is not defined");error_=1;} //true if pt==nan
   if(eta   !=eta)   {ATH_MSG_ERROR("ERROR eta of track  is not defined");error_=1;}
   if(phi_0 !=phi_0) {ATH_MSG_ERROR("ERROR phi of track  is not defined");error_=1;}
   if(error_==1){
     ATH_MSG_ERROR("Original Particle Momentum(px,py,pz,e,m)="<<parent->momentum().px()<<"  "
                   <<parent->momentum().py()<<"  "
                   <<parent->momentum().pz()<<"  "
                   <<parent->momentum().e() <<"  "
                   <<parent->momentum().m() <<"  ");
   }
  
  
  
  
   //Call the appropriate function to set the vn values
   for(int ihar = 0; ihar< 6; ihar++){m_v_n  [ihar]=0.0;} //reset the vn for this particle
   double b = hijing_pars->get_b();
   (*this.*m_flow_function)(b,eta,pt);//get the vn for this particle
   if(m_flow_fluctuations){//add EbE fluctuations
     for(int ihar = 0; ihar< 6; ihar++){
       m_v_n[ihar] *= m_EbE_Multiplier_vn[ihar];
       if(m_v_n[ihar]>0.5) {ATH_MSG_WARNING(" Vn Too large  "<<ihar+1<<"  "<<m_EbE_Multiplier_vn[ihar]<<"  "<<m_v_n[ihar]);m_v_n[ihar]=0.5;}
     }
   }
  
   double phishift=0;
  
   // Old fashioned rotation(approximate)
   if (m_flow_implementation_type==0){
     float phi=phi_0;
     phishift=  -2*( m_v_n[0]*sin(1*(phi-m_psi_n[0]))/1.0 + m_v_n[1]*sin(2*(phi-m_psi_n[1]))/2.0 +
                     m_v_n[2]*sin(3*(phi-m_psi_n[2]))/3.0 + m_v_n[3]*sin(4*(phi-m_psi_n[3]))/4.0 +
                     m_v_n[4]*sin(5*(phi-m_psi_n[4]))/5.0 + m_v_n[5]*sin(6*(phi-m_psi_n[5]))/6.0 );
  
   }
  
   // New fashioned rotation(exact)
   else if (m_flow_implementation_type==1){
     // Thread-safe according to https://www.gnu.org/software/gsl/doc/html/roots.html
     const gsl_root_fsolver_type *T ATLAS_THREAD_SAFE = gsl_root_fsolver_brent;
     gsl_root_fsolver *s = gsl_root_fsolver_alloc (T);
     double x_lo=-2*M_PI,x_hi=2*M_PI;
     float params[13];
     for(int ipar=0;ipar<13;ipar++) {params[ipar]=0;}
     gsl_function F;
     F.function = &(AddFlowByShifting::vn_func);
     F.params =&params;
     gsl_root_fsolver_set (s, &F, x_lo, x_hi);
     int iter=0;
     params[0]=phi_0;
     params[1]=m_v_n  [0]; params[2]=m_v_n  [1]; params[3]=m_v_n  [2]; params[4 ]=m_v_n  [3]; params[5 ]=m_v_n  [4];  params[6 ]=m_v_n  [5];
     params[7]=m_psi_n[0]; params[8]=m_psi_n[1]; params[9]=m_psi_n[2]; params[10]=m_psi_n[3]; params[11]=m_psi_n[4];  params[12]=m_psi_n[5];
     int status;
     double phi=phi_0;
     do
       {
         iter++;
         status = gsl_root_fsolver_iterate (s);
         phi = gsl_root_fsolver_root (s);
         x_lo = gsl_root_fsolver_x_lower (s);
         x_hi = gsl_root_fsolver_x_upper (s);
         status = gsl_root_test_interval (x_lo, x_hi,0, 0.00001);
       }
     while (status == GSL_CONTINUE && iter < 1000);
     gsl_root_fsolver_free (s);
  
     if (iter>=1000) return 0;
  
     phishift = phi-phi_0;
   }
  
   if(std::abs(phishift) > 1e-7) {
     momentum.rotateZ(phishift*Gaudi::Units::rad);
     parent->set_momentum( HepMC::FourVector(momentum.px(),momentum.py(),momentum.pz(),momentum.e()) );
   }
   ATH_MSG_DEBUG( "Parent particle: V1 = " << m_v_n[0] <<
                  " V2 = " << m_v_n[1] <<
                  " V3 = " << m_v_n[2] <<
                  " V4 = " << m_v_n[3] <<
                  " V5 = " << m_v_n[4] <<
                  " V6 = " << m_v_n[5] <<
                  " Phi shift = " << phishift <<
                  " Phi shifted = " << momentum.phi() );
  
   return phishift;
 }

◆ ao_test()

void AddFlowByShifting::ao_test	(	double	b,
		double	eta,
		double	pt
	)

private

Definition at line 542 of file AddFlowByShifting.cxx.

 {
   pt/=1000;
   if(pt>2) pt = 2; // flat max at pt > 2
   m_v_n[1] = 0.02 * b * pt;
   m_v_n[0] = 0;
 }

◆ custom_vn()

void AddFlowByShifting::custom_vn	(	double	b,
		double	eta,
		double	pt
	)

private

Definition at line 551 of file AddFlowByShifting.cxx.

 {
   m_v_n[0]=m_custom_v1;
   m_v_n[1]=m_custom_v2;
   m_v_n[2]=m_custom_v3;
   m_v_n[3]=m_custom_v4;
   m_v_n[4]=m_custom_v5;
   m_v_n[5]=m_custom_v6;
 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	hndl,
		const SG::VarHandleKeyArrayType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	hndl,
		const SG::VarHandleKeyType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	hndl,
		const SG::VarHandleType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
   }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                      {
     typedef typename SG::HandleClassifier<T>::type htype;
     return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
   }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

StatusCode AddFlowByShifting::execute ( )

Definition at line 159 of file AddFlowByShifting.cxx.

                                       {
   ATH_MSG_INFO(">>> AddFlowByShifting from execute");
  
   const EventContext& ctx = Gaudi::Hive::currentContext();
   CLHEP::HepRandomEngine *rndmEngine = getRandomEngine("FLOW", ctx);
   // Get hijing event parameters
   const HijingEventParams *hijing_pars;
   if( evtStore()->retrieve(hijing_pars, "Hijing_event_params").isFailure() ) {
     ATH_MSG_ERROR("Could not retrieve Hijing_event_params");
     return StatusCode::FAILURE;
   }
   ATH_MSG_INFO("Event parameters: B = " << hijing_pars->get_b()<<
                "  BPhi = " << hijing_pars->get_bphi());
  
  
   // FIXME: changing data in the event store
   HijingEventParams *hijing_pars_nc = const_cast<HijingEventParams*> (hijing_pars);
  
  
   // Read Data from Transient Store
   const McEventCollection* mcCollptr;
   if ( evtStore()->retrieve(mcCollptr, m_inkey).isFailure() ) {
     ATH_MSG_ERROR("Could not retrieve truth McEventCollection");
     return StatusCode::FAILURE;
   }
  
  
   // Loop over all events in original McEventCollection and
   // Copy to a new (modifiable) collection
   McEventCollection::const_iterator citr;
   McEventCollection* mcFlowCollptr = new  McEventCollection();
   for (citr = mcCollptr->begin(); citr!=mcCollptr->end(); ++citr) {
     mcFlowCollptr->push_back(new HepMC::GenEvent(*(*citr)));
   }
  
  
   //Geneate the event-plane angles (some of them may or may not be used later on)
   //Store the angles into the hijing event parameters
   for(int ihar=0;ihar<6;ihar++){
     m_psi_n[ihar] =(CLHEP::RandFlat::shoot(rndmEngine)-0.5)*2*M_PI / (ihar+1);   //Principal value must be within -PI/n to PI/n
     hijing_pars_nc->set_psi(ihar+1,m_psi_n[ihar]);
   }
   m_psi_n[1]=hijing_pars->get_bphi()                   ;//the psi2 plane is aligned with the impact parameter
   m_psi_n[1]=std::atan2(std::sin(2*m_psi_n[1]),std::cos(2*m_psi_n[1]))/2.0;//ensure that Psi2 is within [-PI/2,PI/2]
   hijing_pars_nc->set_psi(2,m_psi_n[1]);
   ATH_MSG_DEBUG(" Psi2 for event : "<<(*hijing_pars).get_psi(2));
  
  
   // Add flow by phi angle shifting
   McEventCollection::iterator itr;
   for (itr = mcFlowCollptr->begin(); itr!=mcFlowCollptr->end(); ++itr) {
     ATH_MSG_DEBUG("Next event in the bag ...");
  
  
  
 #ifdef HEPMC3
     auto mainvtx=(*itr)->vertices().front();
     if(m_flow_fluctuations) Set_EbE_Fluctuation_Multipliers(mainvtx,hijing_pars->get_b(),rndmEngine);
     int particles_in_event = (*itr)->particles().size();
     m_particles_processed = 0;
     for ( auto parent: mainvtx->particles_out())
 #else
     auto mainvtx=*((*itr)->vertices_begin());
     if(m_flow_fluctuations) Set_EbE_Fluctuation_Multipliers(mainvtx,hijing_pars->get_b(),rndmEngine);
       int particles_in_event = (*itr)->particles_size();
     m_particles_processed = 0;
     for ( auto parent: *mainvtx)
 #endif
       {
         // Process particles from main vertex
         CLHEP::HepLorentzVector momentum(parent->momentum().px(),
                                          parent->momentum().py(),
                                          parent->momentum().pz(),
                                          parent->momentum().e());
         ATH_MSG_DEBUG("Parent particle: " << parent        <<
                       " Eta = "           << momentum.pseudoRapidity()<<
                       " Phi = "           << momentum.phi()            );
  
         //skip particle if eta is outside implementation range
         if(m_floweta_sw){
           float eta=std::abs(momentum.pseudoRapidity());
           if (eta<m_flow_mineta || eta> m_flow_maxeta) continue;
         }
  
         //skip particle if pT is outside implementation range
         if(m_flowpt_sw){
           float pT=momentum.perp();
           if (pT<m_flow_minpt || pT> m_flow_maxpt) continue;
         }
  
         // Randomize phi if explicitely requested
         if(m_ranphi_sw) {
           double phishift = SetParentToRanPhi(parent, rndmEngine);
           MoveDescendantsToParent(parent, phishift)  ;// adjust decsandants to parent position
         }
  
         // Add flow to particles from main vertex
         double phishift = AddFlowToParent(parent, hijing_pars);
         MoveDescendantsToParent(std::move(parent), phishift);// adjust decsandants to parent position
       }
  
     // correct for double counting
     if(m_ranphi_sw) m_particles_processed /= 2;
     // correct for incoming particles
     ATH_MSG_INFO( " Particles in event: " << particles_in_event <<
                   " Processed for flow: " << m_particles_processed+2);
   }
  
   if(evtStore()->record(mcFlowCollptr, m_outkey).isFailure()){
     ATH_MSG_ERROR("Could not record flow McEventCollection");
     return StatusCode::FAILURE;
   }
   return StatusCode::SUCCESS;
 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

 {
   // If we didn't find any symlinks to add, just return the collection
   // from the base class.  Otherwise, return the extended collection.
   if (!m_extendedExtraObjects.empty()) {
     return m_extendedExtraObjects;
   }
   return Algorithm::extraOutputDeps();
 }

◆ finalize()

StatusCode AddFlowByShifting::finalize ( )

Definition at line 275 of file AddFlowByShifting.cxx.

                                        {
   ATH_MSG_INFO(">>> AddFlowByShifting from finalize <<<");
   // End of finalization step
   return StatusCode::SUCCESS;
 }

◆ fixed_v2()

void AddFlowByShifting::fixed_v2	(	double	b,
		double	eta,
		double	pt
	)

private

Definition at line 523 of file AddFlowByShifting.cxx.

 {
   m_v_n[0]=0.0000; m_v_n[1]=0.0500;
   m_v_n[2]=0.0000; m_v_n[3]=0.0000;
   m_v_n[4]=0.0000; m_v_n[5]=0.0000;
 }

◆ fixed_vn()

void AddFlowByShifting::fixed_vn	(	double	b,
		double	eta,
		double	pt
	)

private

Definition at line 514 of file AddFlowByShifting.cxx.

 {
   m_v_n[0]=0.0000; m_v_n[1]=0.0500;
   m_v_n[2]=0.0280; m_v_n[3]=0.0130;
   m_v_n[4]=0.0045; m_v_n[5]=0.0015;
 }

◆ getRandomEngine()

CLHEP::HepRandomEngine * AddFlowByShifting::getRandomEngine	(	const std::string &	streamName,
		const EventContext &	ctx
	)		const

private

Definition at line 149 of file AddFlowByShifting.cxx.

 {
   ATHRNG::RNGWrapper* rngWrapper = m_rndmSvc->getEngine(this, streamName);
   std::string rngName = name()+streamName;
   rngWrapper->setSeed( rngName, ctx );
   return rngWrapper->getEngine(ctx);
 }

◆ initialize()

StatusCode AddFlowByShifting::initialize ( )

Definition at line 77 of file AddFlowByShifting.cxx.

                                         {
   ATH_MSG_INFO(">>> AddFlowByShifting from Initialize <<<");
  
   ATH_CHECK(m_rndmSvc.retrieve());
  
   ATH_MSG_INFO("**********Settings for Afterburner************");
   ATH_MSG_INFO("McTruthKey         : " << m_inkey              );
   ATH_MSG_INFO("McFlowKey          : " << m_outkey             );
  
   ATH_MSG_INFO("FlowFunctionName   : " << m_flow_function_name );
   ATH_MSG_INFO("FlowInplementation : " << m_flow_implementation);
   ATH_MSG_INFO("FlowFluctuations   : " << m_flow_fluctuations  );
  
   ATH_MSG_INFO("RandomizePhi       : " << m_ranphi_sw          );
  
   ATH_MSG_INFO("FlowEtaSwitch      : " << m_floweta_sw         );
   ATH_MSG_INFO("FlowMinEtaCut      : " << m_flow_mineta        );
   ATH_MSG_INFO("FlowMaxEtaCut      : " << m_flow_maxeta        );
  
   ATH_MSG_INFO("FlowPtSwitch       : " << m_flowpt_sw          );
   ATH_MSG_INFO("FlowMinPtCut       : " << m_flow_minpt         );
   ATH_MSG_INFO("FlowMaxPtCut       : " << m_flow_maxpt         );
  
   ATH_MSG_INFO("FlowV1             : " << m_custom_v1          );
   ATH_MSG_INFO("FlowV2             : " << m_custom_v2          );
   ATH_MSG_INFO("FlowV3             : " << m_custom_v3          );
   ATH_MSG_INFO("FlowV4             : " << m_custom_v4          );
   ATH_MSG_INFO("FlowV5             : " << m_custom_v5          );
   ATH_MSG_INFO("FlowV6             : " << m_custom_v6          );
   ATH_MSG_INFO("FlowBSwitch        : " << m_flowb_sw           );
   ATH_MSG_INFO("********************************r*************");
  
  
   // Select the flow-implementing function based of the function-choice variable
   m_flow_function=&AddFlowByShifting::jjia_minbias_new;
   if     (m_flow_function_name=="jjia_minbias_new"       ) m_flow_function=&AddFlowByShifting::jjia_minbias_new;
   else if(m_flow_function_name=="jjia_minbias_new_v2only") m_flow_function=&AddFlowByShifting::jjia_minbias_new_v2only;
   else if(m_flow_function_name=="fixed_vn"               ) m_flow_function=&AddFlowByShifting::fixed_vn;
   else if(m_flow_function_name=="fixed_v2"               ) m_flow_function=&AddFlowByShifting::fixed_v2;
   else if(m_flow_function_name=="jjia_minbias_old"       ) m_flow_function=&AddFlowByShifting::jjia_minbias_old;
   else if(m_flow_function_name=="ao_test"                ) m_flow_function=&AddFlowByShifting::ao_test;
   else if(m_flow_function_name=="custom"                 ) m_flow_function=&AddFlowByShifting::custom_vn;
   else if(m_flow_function_name=="p_Pb_cent_eta_indep"    ) m_flow_function=&AddFlowByShifting::p_Pb_cent_eta_indep;
  
   m_flow_implementation_type=1;
   if(m_flow_implementation=="approximate") m_flow_implementation_type=0;
   if(m_flow_implementation=="exact"      ) m_flow_implementation_type=1;
  
  
  
   //TGraph storing the v2_RP/delta Vs b_imp values to be used in implementing the EbyE fluctuations
   //the values below are b_imp-low,b_imp-high, delta/v2_RP for different centralities
   //underflow and overflow bins are added for smooth extrapolation
   //The delta/v2_RP values are taken from Fig15 of EbE vn paper (arXiv:1305.2942)
   //                 <0  ,  0-1 ,  1-2 , 2-3  ,  3-4 ,  4-5 , 5-10 , 10-15, 15-20, 20-25, 25-30,
   float b_lo[21]={  -1.00, 0.000, 1.483, 2.098, 2.569, 2.966, 3.317, 4.687, 5.739, 6.627, 7.409,
                     8.117, 8.767, 9.373, 9.943,10.479,10.991,11.479,11.947,15.00 ,100.0};
   //                30-35, 35-40, 40-45, 45-50, 50-55, 55-60, 60-65, 65-70, 70-
   float b_hi[21]={ -1.00 , 1.483, 2.098, 2.569, 2.966, 3.317, 4.687, 5.739, 6.627, 7.409, 8.117,//bimp_high
                    8.767, 9.373, 9.943,10.479,10.991,11.479,11.947,12.399,15.00 ,100.0};
   float val [21]={  5.600, 5.600, 5.600,1.175 ,0.8253,0.7209,0.5324,0.4431,0.3984,0.3844,0.3847,
                     0.3935,0.4106,0.4310,0.4574,0.4674,0.4873,0.4796,0.4856,0.5130,0.5130};
   float bimp_vals[21];
   for(int i=0;i<21;i++) bimp_vals[i]=(b_lo[i]+b_hi[i])/2.0;
   if(m_flow_fluctuations) m_graph_fluc=new TGraph(21,bimp_vals,val);
  
  
   // Initialization terminated
   return StatusCode::SUCCESS;
 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ jjia_minbias_new()

void AddFlowByShifting::jjia_minbias_new	(	double	b,
		double	eta,
		double	pt
	)

private

Definition at line 460 of file AddFlowByShifting.cxx.

 {
   pt=pt/1000.0; //convert to GeV
  
   float a1,a2,a3,a4;
   a1=0.4397*std::exp(-(b-4.526)*(b-4.526)/72.0) + 0.636;
   a2=1.916/(b+2) +0.1;
   a3=4.79*0.0001*(b-0.621)*(b-10.172)*(b-23)+1.2;   // this is >0 for b>0
   a4=0.135*std::exp(-0.5*(b-10.855)*(b-10.855)/4.607/4.607) +0.0120;
  
   float temp1 = std::pow(pt    , a1) / (1+std::exp( (pt-3.0)/a3));
   float temp2 = std::pow(pt+0.1,-a2) / (1+std::exp(-(pt-4.5)/a3));
   float temp3 =  0.01           / (1+std::exp(-(pt-4.5)/a3));
  
   m_v_n[1] = ( a4*(temp1+temp2) + temp3 )* std::exp(-0.5* eta*eta /6.27/6.27) ;
  
   float fb=0.97 +1.06*std::exp(-0.5*b*b/3.2/3.2);
   m_v_n[2]=std::pow(fb*std::sqrt(m_v_n[1]),3);
  
   float gb= 1.096 +1.36 *std::exp(-0.5*b*b/3.0/3.0);
   gb=gb*sqrt(m_v_n[1]);
   m_v_n[3]=pow(gb,4);
   m_v_n[4]=pow(gb,5);
   m_v_n[5]=pow(gb,6);
   m_v_n[0]=0;
 }

◆ jjia_minbias_new_v2only()

void AddFlowByShifting::jjia_minbias_new_v2only	(	double	b,
		double	eta,
		double	pt
	)

private

Definition at line 489 of file AddFlowByShifting.cxx.

 {
   pt=pt/1000.0; //convert to GeV
  
   float a1,a2,a3,a4;
   a1=0.4397*std::exp(-(b-4.526)*(b-4.526)/72.0) + 0.636;
   a2=1.916/(b+2) +0.1;
   a3=4.79*0.0001*(b-0.621)*(b-10.172)*(b-23)+1.2;   // this is >0 for b>0
   a4=0.135*std::exp(-0.5*(b-10.855)*(b-10.855)/4.607/4.607) +0.0120;
  
   float temp1 = std::pow(pt    , a1) / (1+std::exp( (pt-3.0)/a3));
   float temp2 = std::pow(pt+0.1,-a2) / (1+std::exp(-(pt-4.5)/a3));
   float temp3 =  0.01           / (1+std::exp(-(pt-4.5)/a3));
  
   m_v_n[1] = ( a4*(temp1+temp2) + temp3 )* std::exp(-0.5* eta*eta /6.27/6.27) ;
  
   m_v_n[0]=0;
   m_v_n[2]=0;
   m_v_n[3]=0;
   m_v_n[4]=0;
   m_v_n[5]=0;
 }

◆ jjia_minbias_old()

void AddFlowByShifting::jjia_minbias_old	(	double	b,
		double	eta,
		double	pt
	)

private

Definition at line 531 of file AddFlowByShifting.cxx.

 {
   m_v_n[0] = 0;
   m_v_n[1] = 0.03968 * b
     * (1 - 2.1/(1 + std::exp(1.357*(pt/1000))))
     * std::exp(-(eta*eta)/(2*6.37*6.37));
   m_v_n[2]=0.0000; m_v_n[3]=0.0000;
   m_v_n[4]=0.0000; m_v_n[5]=0.0000;
 }

◆ MoveDescendantsToParent()

void AddFlowByShifting::MoveDescendantsToParent	(	HepMC::GenParticlePtr	parent,
		double	phishift
	)

private

Definition at line 309 of file AddFlowByShifting.cxx.

 {
   // Move the branch of descendant vertices and particles
   // by phishift to parent particle position
   auto endvtx = parent->end_vertex();
   if ( endvtx ) {
     ATH_MSG_DEBUG("Processing branch of parent particle "<< parent);
  
     // now rotate descendant vertices
 #ifdef HEPMC3
     for (HepMC::GenVertexPtr descvtx:  HepMC::descendant_vertices(endvtx)) {
 #else
       for ( HepMC::GenVertex::vertex_iterator
               descvtxit = endvtx->vertices_begin(HepMC::descendants);
             descvtxit != endvtx->vertices_end(HepMC::descendants);
             ++descvtxit) {
         auto descvtx = (*descvtxit);
 #endif
  
           ATH_MSG_DEBUG("Processing vertex " << descvtx);
  
           // rotate vertex
           if(std::abs(phishift) > 1e-7) {
             CLHEP::HepLorentzVector position(descvtx->position().x(),
                                              descvtx->position().y(),
                                              descvtx->position().z(),
                                              descvtx->position().t());
             position.rotateZ(phishift*Gaudi::Units::rad);
             descvtx->set_position(HepMC::FourVector( position.x(),position.y(),position.z(),position.t()) );
           }
  
           // now rotate their associated particles
 #ifdef HEPMC3
           for (auto descpart: descvtx->particles_out())
 #else
             for (auto descpart: *descvtx)
 #endif
               {
                 ATH_MSG_DEBUG("Descendant particle: " << descpart <<
                               " Eta = "               << descpart->momentum().pseudoRapidity() <<
                               " Phi = "               << descpart->momentum().phi() );
               }
         }
   }
   return;
 }

◆ msg() [1/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( const MSG::Level lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

bool AthCommonMsg< Algorithm >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ p_Pb_cent_eta_indep()

void AddFlowByShifting::p_Pb_cent_eta_indep	(	double	b,
		double	eta,
		double	pt
	)

private

Definition at line 562 of file AddFlowByShifting.cxx.

 {
   pt=pt/1000.0; //convert to GeV
  
   float an_val[4][3];
  
   an_val[0][0] = 0.1149;
   an_val[0][1] = 1.181;
   an_val[0][2] = 0.3767;
  
   an_val[1][0] = 0.0498;
   an_val[1][1] = 1.688;
   an_val[1][2] = 0.5046;
  
   an_val[2][0] = 0.02095;
   an_val[2][1] = 2.196;
   an_val[2][2] = 0.6259;
  
   an_val[3][0] = 0.00682*0.5;//added in 0.5 factor by hand
   an_val[3][1] = 4.938;
   an_val[3][2] = 1.237;
  
   m_v_n[0]=0;
   m_v_n[1]=an_val[0][0]*std::pow(pt,an_val[0][1])*std::exp(-an_val[0][2]*pt);
   m_v_n[2]=an_val[1][0]*std::pow(pt,an_val[1][1])*std::exp(-an_val[1][2]*pt);
   m_v_n[3]=an_val[2][0]*std::pow(pt,an_val[2][1])*std::exp(-an_val[2][2]*pt);
   m_v_n[4]=an_val[3][0]*std::pow(pt,an_val[3][1])*std::exp(-an_val[3][2]*pt);
   m_v_n[5]=0;
 }

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

   {
     h.renounce();
     PBASE::renounce (h);
   }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ Set_EbE_Fluctuation_Multipliers()

void AddFlowByShifting::Set_EbE_Fluctuation_Multipliers	(	HepMC::GenVertexPtr	mainvtx,
		float	b,
		CLHEP::HepRandomEngine *	rndmEngine
	)

private

Definition at line 593 of file AddFlowByShifting.cxx.

                                                                                                                           {
   int Total_Multiplicity=0;
   double EbE_Vn[6];
   for(int ihar=0;ihar<6;ihar++){m_EbE_Multiplier_vn[ihar]=1.0;EbE_Vn[ihar]=0.0;}
  
   for(auto  parent: *mainvtx) {
     float eta= parent->momentum().pseudoRapidity();
     float pT = parent->momentum().perp();
  
     for(int ihar = 0; ihar< 6; ihar++){m_v_n  [ihar]=0.0;}
     (*this.*m_flow_function)(b,eta,pT);
     for(int ihar = 0; ihar< 6; ihar++){EbE_Vn[ihar] += m_v_n  [ihar];}
     Total_Multiplicity++;
   }
  
  
   for(int ihar = 0; ihar< 6; ihar++){m_v_n  [ihar]=0.0;}//keep vn as zero before we return from function
   if(Total_Multiplicity<=0) return;
  
  
   for(int ihar=0;ihar<6;ihar++){
     EbE_Vn[ihar]/=Total_Multiplicity;
     float vn_rp=0,delta=0;//BG parameterizations
  
     //in the following we assume that the vn in the event is sqrt(<vn^2>)
     //This is because the vn(pT) were tuned to 2PC/EP measurements
     //then : vn_evt=vn_rp^2 + 2*(delta^2)
     //which is used together with the "alpha" to get the "vn_rp" and "delta"
     if(ihar==0) continue; //No EbE fluctuation for v1
     else if(ihar==1) {    //v2
       float alpha=1.0/m_graph_fluc->Eval(b);// ratio of V2_RP over delta (from Fig 15 of EbE vn paper)
       delta=EbE_Vn[ihar]/sqrt(2.0+alpha*alpha);
       vn_rp=alpha*delta;
     }
     else if(ihar>=2) { //v3-v6
       vn_rp =0;
       delta=EbE_Vn[ihar]/sqrt(2.0);
     }
     if(EbE_Vn[ihar]==0) {ATH_MSG_WARNING("Zero EbeV"<<ihar+1); continue;}
     float X=CLHEP::RandGaussQ::shoot(rndmEngine,vn_rp,delta);
     float Y=CLHEP::RandGaussQ::shoot(rndmEngine,0.0  ,delta);
     m_EbE_Multiplier_vn[ihar]=sqrt(X*X+ Y*Y)/EbE_Vn[ihar];
     ATH_MSG_INFO("EbE_Multiplier_v"<<ihar+1<<"="<<m_EbE_Multiplier_vn[ihar]);
   }
 }

◆ SetParentToRanPhi()

double AddFlowByShifting::SetParentToRanPhi	(	HepMC::GenParticlePtr	parent,
		CLHEP::HepRandomEngine *	rndmEngine
	)

private

Definition at line 282 of file AddFlowByShifting.cxx.

 {
   // Set particle to random phi
   // Return phi shift
   m_particles_processed++;
  
   double phi, phishift;
   CLHEP::HepLorentzVector momentum(parent->momentum().px(),
                                    parent->momentum().py(),
                                    parent->momentum().pz(),
                                    parent->momentum().e());
   phi = momentum.phi();
  
   double rannum = CLHEP::RandFlat::shoot(rndmEngine);
   double ranphi = (rannum-0.5)*2*M_PI;
   phishift = ranphi - phi;
  
   momentum.setPhi(ranphi*Gaudi::Units::rad);
   parent->set_momentum(  HepMC::FourVector(momentum.px(),momentum.py(),momentum.pz(),momentum.e()) );
  
   ATH_MSG_INFO("Parent phi randomized = " << momentum.phi());
  
   return phishift;
 }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

                                        {
   StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
  
   if (sc.isFailure()) {
     return sc;
   }
   ServiceHandle<ICondSvc> cs("CondSvc",name());
   for (auto h : outputHandles()) {
     if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
       // do this inside the loop so we don't create the CondSvc until needed
       if ( cs.retrieve().isFailure() ) {
         ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
         return StatusCode::SUCCESS;
       }
       if (cs->regHandle(this,*h).isFailure()) {
         sc = StatusCode::FAILURE;
         ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                       << " with CondSvc");
       }
     }
   }
   return sc;
 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                                {
     // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
     //         << "  size: " << m_vhka.size() << endmsg;
     for (auto &a : m_vhka) {
       std::vector<SG::VarHandleKey*> keys = a->keys();
       for (auto k : keys) {
         k->setOwner(this);
       }
     }
   }

◆ vn_func()

double AddFlowByShifting::vn_func	(	double	x,
		void *	params
	)

static

Definition at line 40 of file AddFlowByShifting.cxx.

 {
   float *par_float = (float*) params;
   float phi_0  = par_float[0];
   float *vn    = par_float+1;
   float *psi_n = vn+6;
   double val=x   +2*(   vn[0]*sin(1*(x-psi_n[0]))/1.0 + vn[1]*sin(2*(x-psi_n[1]))/2.0 +
                         vn[2]*sin(3*(x-psi_n[2]))/3.0 + vn[3]*sin(4*(x-psi_n[3]))/4.0 +
                         vn[4]*sin(5*(x-psi_n[4]))/5.0 + vn[5]*sin(6*(x-psi_n[5]))/6.0 );
   return val-phi_0;
 }

◆ vn_func_derivative()

double AddFlowByShifting::vn_func_derivative	(	double	x,
		void *	params
	)

static

Definition at line 53 of file AddFlowByShifting.cxx.

 {
   float *par_float = (float*) params;
   float *vn    = par_float+1;
   float *psi_n = vn+6;
   double val=1   +2*(   vn[0]*cos(1*(x-psi_n[0]))/1.0 + vn[1]*cos(2*(x-psi_n[1]))/2.0 +
                         vn[2]*cos(3*(x-psi_n[2]))/3.0 + vn[3]*cos(4*(x-psi_n[3]))/4.0 +
                         vn[4]*cos(5*(x-psi_n[4]))/5.0 + vn[5]*cos(6*(x-psi_n[5]))/6.0 );
   return val;
 }

Member Data Documentation

◆ m_custom_v1

FloatProperty AddFlowByShifting::m_custom_v1 {this, "custom_v1", 0.f}

private

Definition at line 97 of file AddFlowByShifting.h.

◆ m_custom_v2

FloatProperty AddFlowByShifting::m_custom_v2 {this, "custom_v2", 0.f}

private

Definition at line 98 of file AddFlowByShifting.h.

◆ m_custom_v3

FloatProperty AddFlowByShifting::m_custom_v3 {this, "custom_v3", 0.f}

private

Definition at line 99 of file AddFlowByShifting.h.

◆ m_custom_v4

FloatProperty AddFlowByShifting::m_custom_v4 {this, "custom_v4", 0.f}

private

Definition at line 100 of file AddFlowByShifting.h.

◆ m_custom_v5

FloatProperty AddFlowByShifting::m_custom_v5 {this, "custom_v5", 0.f}

private

Definition at line 101 of file AddFlowByShifting.h.

◆ m_custom_v6

FloatProperty AddFlowByShifting::m_custom_v6 {this, "custom_v6", 0.f}

private

Definition at line 102 of file AddFlowByShifting.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_EbE_Multiplier_vn

float AddFlowByShifting::m_EbE_Multiplier_vn[6]

private

Definition at line 108 of file AddFlowByShifting.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_flow_fluctuations

BooleanProperty AddFlowByShifting::m_flow_fluctuations {this, "FlowFluctuations", false}

private

Definition at line 86 of file AddFlowByShifting.h.

◆ m_flow_function

void(AddFlowByShifting::* AddFlowByShifting::m_flow_function) (double b, double eta, double pt)

private

Definition at line 61 of file AddFlowByShifting.h.

◆ m_flow_function_name

StringProperty AddFlowByShifting::m_flow_function_name {this, "FlowFunctionName", "jjia_minbias_new"}

private

Definition at line 83 of file AddFlowByShifting.h.

◆ m_flow_implementation

StringProperty AddFlowByShifting::m_flow_implementation {this, "FlowImplementation", "exact"}

private

Definition at line 84 of file AddFlowByShifting.h.

◆ m_flow_implementation_type

int AddFlowByShifting::m_flow_implementation_type {0}

private

Definition at line 85 of file AddFlowByShifting.h.

◆ m_flow_maxeta

FloatProperty AddFlowByShifting::m_flow_maxeta {this, "FlowMaxEtaCut", 10.0}

private

Definition at line 89 of file AddFlowByShifting.h.

◆ m_flow_maxpt

FloatProperty AddFlowByShifting::m_flow_maxpt {this, "FlowMaxPtCut", 1000000.f}

private

Definition at line 93 of file AddFlowByShifting.h.

◆ m_flow_mineta

FloatProperty AddFlowByShifting::m_flow_mineta {this, "FlowMinEtaCut", 0.f}

private

Definition at line 90 of file AddFlowByShifting.h.

◆ m_flow_minpt

FloatProperty AddFlowByShifting::m_flow_minpt {this, "FlowMinPtCut", 0.f}

private

Definition at line 94 of file AddFlowByShifting.h.

◆ m_flowb_sw

IntegerProperty AddFlowByShifting::m_flowb_sw {this, "FlowBSwitch", 0}

private

Definition at line 96 of file AddFlowByShifting.h.

◆ m_floweta_sw

IntegerProperty AddFlowByShifting::m_floweta_sw {this, "FlowEtaSwitch", 0 }

private

Definition at line 88 of file AddFlowByShifting.h.

◆ m_flowpt_sw

IntegerProperty AddFlowByShifting::m_flowpt_sw {this, "FlowPtSwitch", 0}

private

Definition at line 92 of file AddFlowByShifting.h.

◆ m_graph_fluc

TGraph* AddFlowByShifting::m_graph_fluc {}

private

Definition at line 71 of file AddFlowByShifting.h.

◆ m_inkey

StringProperty AddFlowByShifting::m_inkey {this, "McTruthKey", "GEN_EVENT"}

private

Definition at line 78 of file AddFlowByShifting.h.

◆ m_outkey

StringProperty AddFlowByShifting::m_outkey {this, "McFlowKey", "FLOW_EVENT"}

private

Definition at line 79 of file AddFlowByShifting.h.

◆ m_particles_processed

int AddFlowByShifting::m_particles_processed {0}

private

Definition at line 104 of file AddFlowByShifting.h.

◆ m_psi_n

float AddFlowByShifting::m_psi_n[6]

private

Definition at line 107 of file AddFlowByShifting.h.

◆ m_ranphi_sw

IntegerProperty AddFlowByShifting::m_ranphi_sw {this, "RandomizePhi", 0}

private

Definition at line 81 of file AddFlowByShifting.h.

◆ m_rndmSvc

ServiceHandle<IAthRNGSvc> AddFlowByShifting::m_rndmSvc {this, "RndmSvc", "AthRNGSvc"}

private

Definition at line 75 of file AddFlowByShifting.h.

◆ m_v_n

float AddFlowByShifting::m_v_n[6]

private

Definition at line 107 of file AddFlowByShifting.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ AddFlowByShifting()

Member Function Documentation

◆ AddFlowToParent()

◆ ao_test()

◆ custom_vn()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ finalize()

◆ fixed_v2()

◆ fixed_vn()

◆ getRandomEngine()

◆ initialize()

◆ inputHandles()

◆ jjia_minbias_new()

◆ jjia_minbias_new_v2only()

◆ jjia_minbias_old()

◆ MoveDescendantsToParent()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ p_Pb_cent_eta_indep()

◆ renounce()

◆ renounceArray()

◆ Set_EbE_Fluctuation_Multipliers()

◆ SetParentToRanPhi()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

◆ vn_func()

◆ vn_func_derivative()

Member Data Documentation

◆ m_custom_v1

◆ m_custom_v2

◆ m_custom_v3

◆ m_custom_v4

◆ m_custom_v5

◆ m_custom_v6

◆ m_detStore

◆ m_EbE_Multiplier_vn

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_flow_fluctuations

◆ m_flow_function

◆ m_flow_function_name

◆ m_flow_implementation

◆ m_flow_implementation_type

◆ m_flow_maxeta

◆ m_flow_maxpt

◆ m_flow_mineta

◆ m_flow_minpt

◆ m_flowb_sw

◆ m_floweta_sw

◆ m_flowpt_sw

◆ m_graph_fluc

◆ m_inkey