#include <CosMuoGenSource.h>

Inheritance diagram for edm::CosMuoGenSource:

Public Member Functions
	CosMuoGenSource (const ParameterSet &, const InputSourceDescription &)

void	endRun (Run &r)

virtual	~CosMuoGenSource ()

Public Member Functions inherited from edm::GeneratedInputSource
	GeneratedInputSource (ParameterSet const &pset, InputSourceDescription const &desc)

virtual	~GeneratedInputSource ()

Public Member Functions inherited from edm::ConfigurableInputSource
	ConfigurableInputSource (ParameterSet const &pset, InputSourceDescription const &desc, bool realData=true)

EventNumber_t	event () const

unsigned int	eventCreationDelay () const

LuminosityBlockNumber_t	luminosityBlock () const

unsigned int	numberEventsInLumi () const

unsigned int	numberEventsInRun () const

unsigned int	numberEventsInThisLumi () const

unsigned int	numberEventsInThisRun () const

TimeValue_t	presentTime () const

RunNumber_t	run () const

unsigned int	timeBetweenEvents () const

virtual	~ConfigurableInputSource ()

Public Member Functions inherited from edm::InputSource
boost::shared_ptr < ActivityRegistry >	actReg () const
	Accessor for Activity Registry. More...

void	closeFile (boost::shared_ptr< FileBlock >, bool cleaningUpAfterException)
	close current file More...

void	doBeginJob ()
	Called by framework at beginning of job. More...

void	doBeginLumi (LuminosityBlockPrincipal &lbp)
	Called by framework at beginning of lumi block. More...

void	doBeginRun (RunPrincipal &rp)
	Called by framework at beginning of run. More...

void	doEndJob ()
	Called by framework at end of job. More...

void	doEndLumi (LuminosityBlockPrincipal &lbp, bool cleaningUpAfterException)
	Called by framework at end of lumi block. More...

void	doEndRun (RunPrincipal &rp, bool cleaningUpAfterException)
	Called by framework at end of run. More...

void	doPostForkReacquireResources (boost::shared_ptr< multicore::MessageReceiverForSource >)

void	doPreForkReleaseResources ()
	Called by the framework before forking the process. More...

ProcessingController::ForwardState	forwardState () const

bool	goToEvent (EventID const &eventID)

	InputSource (ParameterSet const &, InputSourceDescription const &)
	Constructor. More...

void	issueReports (EventID const &eventID)
	issue an event report More...

LuminosityBlockNumber_t	luminosityBlock () const
	Accessor for current luminosity block number. More...

boost::shared_ptr < LuminosityBlockAuxiliary >	luminosityBlockAuxiliary () const
	Called by the framework to merge or insert lumi in principal cache. More...

int	markLumi ()
	Mark lumi as read. More...

int	markRun ()
	Mark run as read. More...

int	maxEvents () const

int	maxLuminosityBlocks () const

ModuleDescription const &	moduleDescription () const
	Accessor for 'module' description. More...

ItemType	nextItemType ()

bool	primary () const
	Accessor for primary input source flag. More...

ProcessConfiguration const &	processConfiguration () const
	Accessor for Process Configuration. More...

std::string const &	processGUID () const
	Accessor for global process identifier. More...

ProcessingMode	processingMode () const
	RunsLumisAndEvents (default), RunsAndLumis, or Runs. More...

boost::shared_ptr < ProductRegistry const >	productRegistry () const
	Accessor for product registry. More...

bool	randomAccess () const

void	readAndCacheLumi (bool merge, HistoryAppender &historyAppender)
	Read next luminosity block. More...

void	readAndCacheRun (bool merge, HistoryAppender &historyAppender)
	Read next run. More...

EventPrincipal *	readEvent (boost::shared_ptr< LuminosityBlockPrincipal > lbCache)

EventPrincipal *	readEvent (EventID const &)
	Read a specific event. More...

boost::shared_ptr< FileBlock >	readFile ()
	Read next file. More...

boost::shared_ptr < LuminosityBlockAuxiliary >	readLuminosityBlockAuxiliary ()
	Read next luminosity block Auxilary. More...

boost::shared_ptr< RunAuxiliary >	readRunAuxiliary ()
	Read next run Auxiliary. More...

ProcessHistoryID const &	reducedProcessHistoryID () const

void	registerProducts ()
	Register any produced products. More...

int	remainingEvents () const

int	remainingLuminosityBlocks () const

void	repeat ()
	Reset the remaining number of events/lumis to the maximum number. More...

ProcessingController::ReverseState	reverseState () const

void	rewind ()
	Begin again at the first event. More...

RunNumber_t	run () const
	Accessor for current run number. More...

boost::shared_ptr< RunAuxiliary >	runAuxiliary () const
	Called by the framework to merge or insert run in principal cache. More...

void	setLuminosityBlockNumber_t (LuminosityBlockNumber_t lb)
	Set the luminosity block ID. More...

void	setRunNumber (RunNumber_t r)
	Set the run number. More...

void	skipEvents (int offset)

Timestamp const &	timestamp () const
	Accessor for the current time, as seen by the input source. More...

virtual	~InputSource ()
	Destructor. More...

Private Member Functions
void	clear ()

virtual bool	produce (Event &e)

Private Attributes
bool	AllMu

double	ClayLayerWidth

bool	cmVerbosity_

CosmicMuonGenerator *	CosMuoGen

double	ELSF

double	extCrossSect

double	extFilterEff

HepMC::GenEvent *	fEvt

double	MaxEn

double	MaxP

double	MaxPh

double	MaxS

double	MaxT

double	MinEn

double	MinP

double	MinP_CMS

double	MinPh

double	MinS

double	MinT

bool	MTCCHalf

bool	MultiMuon

int32_t	MultiMuonFileFirstEvent

std::string	MultiMuonFileName

int32_t	MultiMuonNmin

double	NuPrdAlt

double	PlugVtx

double	PlugVtz

int32_t	RanS

double	RTarget

bool	TIFOnly_constant

bool	TIFOnly_linear

bool	TrackerOnly

double	VarRhoAir

double	VarRhoClay

double	VarRhoPlug

double	VarRhoRock

double	VarRhoWall

double	ZCTarget

double	ZTarget

Additional Inherited Members
Public Types inherited from edm::InputSource
enum	ItemType { IsInvalid, IsStop, IsFile, IsRun, IsLumi, IsEvent, IsRepeat }

enum	ProcessingMode { Runs, RunsAndLumis, RunsLumisAndEvents }

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::ConfigurableInputSource
static void	fillDescription (ParameterSetDescription &desc)

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

static void	fillDescription (ParameterSetDescription &desc)

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::ConfigurableInputSource
void	reallyReadEvent ()

void	setEventNumber (EventNumber_t e)

void	setTime (TimeValue_t t)

Protected Member Functions inherited from edm::InputSource
void	decreaseRemainingEventsBy (int iSkipped)

EventPrincipal *	eventPrincipalCache ()

boost::shared_ptr < LuminosityBlockPrincipal > const	luminosityBlockPrincipal () const

PrincipalCache const &	principalCache () const

PrincipalCache &	principalCache ()

ProductRegistry &	productRegistryUpdate () const

void	reset () const

void	resetLuminosityBlockAuxiliary () const

void	resetRunAuxiliary () const

boost::shared_ptr < RunPrincipal > const	runPrincipal () const

void	setLuminosityBlockAuxiliary (LuminosityBlockAuxiliary *lbp)

void	setRunAuxiliary (RunAuxiliary *rp)

void	setTimestamp (Timestamp const &theTime)
	To set the current time, as seen by the input source. More...

ItemType	state () const

Detailed Description

Definition at line 18 of file CosMuoGenSource.h.

Constructor & Destructor Documentation

edm::CosMuoGenSource::CosMuoGenSource	(	const ParameterSet &	pset,
		const InputSourceDescription &	desc
	)

Definition at line 12 of file CosMuoGenSource.cc.

                                                                                                    :
   GeneratedInputSource(pset, desc ) ,  
   //RanS(pset.getParameter<int>("RanSeed", 123456)), //get seed now from Framework
   MinP(pset.getParameter<double>("MinP")),
   MinP_CMS(pset.getParameter<double>("MinP_CMS")),
   MaxP(pset.getParameter<double>("MaxP")),
   MinT(pset.getParameter<double>("MinTheta")),
   MaxT(pset.getParameter<double>("MaxTheta")),
   MinPh(pset.getParameter<double>("MinPhi")),
   MaxPh(pset.getParameter<double>("MaxPhi")),
   MinS(pset.getParameter<double>("MinT0")),
   MaxS(pset.getParameter<double>("MaxT0")),
   ELSF(pset.getParameter<double>("ElossScaleFactor")),
   RTarget(pset.getParameter<double>("RadiusOfTarget")),
   ZTarget(pset.getParameter<double>("ZDistOfTarget")),
   ZCTarget(pset.getParameter<double>("ZCentrOfTarget")),
   TrackerOnly(pset.getParameter<bool>("TrackerOnly")),
   MultiMuon(pset.getParameter<bool>("MultiMuon")),
   MultiMuonFileName(pset.getParameter<std::string>("MultiMuonFileName")),
   MultiMuonFileFirstEvent(pset.getParameter<int>("MultiMuonFileFirstEvent")),
   MultiMuonNmin(pset.getParameter<int>("MultiMuonNmin")),
   TIFOnly_constant(pset.getParameter<bool>("TIFOnly_constant")),
   TIFOnly_linear(pset.getParameter<bool>("TIFOnly_linear")),
   MTCCHalf(pset.getParameter<bool>("MTCCHalf")),
   PlugVtx(pset.getParameter<double>("PlugVx")),
   PlugVtz(pset.getParameter<double>("PlugVz")),
   VarRhoAir(pset.getParameter<double>("RhoAir")),
   VarRhoWall(pset.getParameter<double>("RhoWall")),
   VarRhoRock(pset.getParameter<double>("RhoRock")),
   VarRhoClay(pset.getParameter<double>("RhoClay")),
   VarRhoPlug(pset.getParameter<double>("RhoPlug")),
   ClayLayerWidth(pset.getParameter<double>("ClayWidth")),
   MinEn(pset.getParameter<double>("MinEnu")),
   MaxEn(pset.getParameter<double>("MaxEnu")),
   NuPrdAlt(pset.getParameter<double>("NuProdAlt")),
   AllMu(pset.getParameter<bool>("AcptAllMu")),
   extCrossSect(pset.getUntrackedParameter<double>("crossSection", -1.)),
   extFilterEff(pset.getUntrackedParameter<double>("filterEfficiency", -1.)),
   cmVerbosity_(pset.getParameter<bool>("Verbosity"))
   {
 
     //if not specified (i.e. negative) then use MinP also for MinP_CMS
     if(MinP_CMS < 0) MinP_CMS = MinP;
 
     //get seed now from Framework
     edm::Service<edm::RandomNumberGenerator> rng;
     RanS = rng->mySeed();
     // set up the generator
     CosMuoGen = new CosmicMuonGenerator();
     CosMuoGen->setNumberOfEvents(numberEventsInRun());
     CosMuoGen->setRanSeed(RanS);
     CosMuoGen->setMinP(MinP);
     CosMuoGen->setMinP_CMS(MinP_CMS);
     CosMuoGen->setMaxP(MaxP);
     CosMuoGen->setMinTheta(MinT);
     CosMuoGen->setMaxTheta(MaxT);
     CosMuoGen->setMinPhi(MinPh);
     CosMuoGen->setMaxPhi(MaxPh);
     CosMuoGen->setMinT0(MinS);
     CosMuoGen->setMaxT0(MaxS);
     CosMuoGen->setElossScaleFactor(ELSF);
     CosMuoGen->setRadiusOfTarget(RTarget);
     CosMuoGen->setZDistOfTarget(ZTarget);
     CosMuoGen->setZCentrOfTarget(ZCTarget);
     CosMuoGen->setTrackerOnly(TrackerOnly);
     CosMuoGen->setMultiMuon(MultiMuon);
     CosMuoGen->setMultiMuonFileName(MultiMuonFileName);
     CosMuoGen->setMultiMuonFileFirstEvent(MultiMuonFileFirstEvent);
     CosMuoGen->setMultiMuonNmin(MultiMuonNmin);
     CosMuoGen->setTIFOnly_constant(TIFOnly_constant);
     CosMuoGen->setTIFOnly_linear(TIFOnly_linear);
     CosMuoGen->setMTCCHalf(MTCCHalf);
     CosMuoGen->setPlugVx(PlugVtx);
     CosMuoGen->setPlugVz(PlugVtz);    
     CosMuoGen->setRhoAir(VarRhoAir);
     CosMuoGen->setRhoWall(VarRhoWall);
     CosMuoGen->setRhoRock(VarRhoRock);
     CosMuoGen->setRhoClay(VarRhoClay);
     CosMuoGen->setRhoPlug(VarRhoPlug);
     CosMuoGen->setClayWidth(ClayLayerWidth);
     CosMuoGen->setMinEnu(MinEn);
     CosMuoGen->setMaxEnu(MaxEn);    
     CosMuoGen->setNuProdAlt(NuPrdAlt);
     CosMuoGen->setAcptAllMu(AllMu);
     CosMuoGen->initialize();
     produces<HepMCProduct>();
     //  fEvt = new HepMC::GenEvent();
     produces<GenRunInfoProduct, edm::InRun>();
   }

edm::CosMuoGenSource::~CosMuoGenSource ( )

virtual

Definition at line 102 of file CosMuoGenSource.cc.

References hitfit::clear().

                                     {
   //CosMuoGen->terminate();
   delete CosMuoGen;
   //  delete fEvt;
   clear();
 }

Member Function Documentation

void edm::CosMuoGenSource::clear ( void )

private

Definition at line 125 of file CosMuoGenSource.cc.

125 {}

void edm::CosMuoGenSource::endRun ( Run & r )

virtual

Reimplemented from edm::ConfigurableInputSource.

Definition at line 110 of file CosMuoGenSource.cc.

References fwrapper::cs, and edm::Run::put().

                                        {
 
   std::auto_ptr<GenRunInfoProduct> genRunInfo(new GenRunInfoProduct());
 
   double cs = CosMuoGen->getRate(); // flux in Hz, not s^-1m^-2
   genRunInfo->setInternalXSec(cs);
   genRunInfo->setExternalXSecLO(extCrossSect);
   genRunInfo->setFilterEfficiency(extFilterEff);
 
   r.put(genRunInfo);
   
   CosMuoGen->terminate();
 }

bool edm::CosMuoGenSource::produce ( Event & e )

privatevirtual

Implements edm::ConfigurableInputSource.

Definition at line 127 of file CosMuoGenSource.cc.

References gather_cfg::cout, Debug, event(), configurableAnalysis::GenParticle, i, edm::Event::put(), and summarizeEdmComparisonLogfiles::success.

 {  
   // generate event
   if (!MultiMuon) {
     CosMuoGen->nextEvent();
   }
   else {
     bool success = CosMuoGen->nextMultiEvent();
     if (!success) return false;
   }
 
   if (Debug) {
     std::cout << "CosMuoGenSource.cc: CosMuoGen->EventWeight=" << CosMuoGen->EventWeight 
           << "  CosMuoGen: Nmuons=" << CosMuoGen->Id_sf.size() << std::endl; 
     std::cout << "CosMuoGen->Id_at=" << CosMuoGen->Id_at
           << "  CosMuoGen->Vx_at=" << CosMuoGen->Vx_at 
           << "  CosMuoGen->Vy_at=" << CosMuoGen->Vy_at
           << "  CosMuoGen->Vz_at=" << CosMuoGen->Vz_at 
           << "  CosMuoGen->T0_at=" << CosMuoGen->T0_at << std::endl;
     std::cout << "  Px=" << CosMuoGen->Px_at
           << "  Py=" << CosMuoGen->Py_at
           << "  Pz=" << CosMuoGen->Pz_at << std::endl;
     for (unsigned int i=0; i<CosMuoGen->Id_sf.size(); ++i) {
       std::cout << "Id_sf[" << i << "]=" << CosMuoGen->Id_sf[i]
         << "  Vx_sf[" << i << "]=" << CosMuoGen->Vx_sf[i]
         << "  Vy_sf=" << CosMuoGen->Vy_sf[i]
         << "  Vz_sf=" << CosMuoGen->Vz_sf[i]
         << "  T0_sf=" << CosMuoGen->T0_sf[i]
         << "  Px_sf=" << CosMuoGen->Px_sf[i]
         << "  Py_sf=" << CosMuoGen->Py_sf[i]
         << "  Pz_sf=" << CosMuoGen->Pz_sf[i] << std::endl;
       std::cout << "phi_sf=" << atan2(CosMuoGen->Px_sf[i],CosMuoGen->Pz_sf[i]) << std::endl;
       std::cout << "Id_ug[" << i << "]=" << CosMuoGen->Id_ug[i] 
         << "  Vx_ug[" << i << "]=" << CosMuoGen->Vx_ug[i] 
         << "  Vy_ug=" << CosMuoGen->Vy_ug[i]
         << "  Vz_ug=" << CosMuoGen->Vz_ug[i]
         << "  T0_ug=" << CosMuoGen->T0_ug[i]
         << "  Px_ug=" << CosMuoGen->Px_ug[i]
         << "  Py_ug=" << CosMuoGen->Py_ug[i]
         << "  Pz_ug=" << CosMuoGen->Pz_ug[i] << std::endl;
       std::cout << "phi_ug=" << atan2(CosMuoGen->Px_ug[i],CosMuoGen->Pz_ug[i]) << std::endl;;
     }
   }
 
 
   fEvt = new HepMC::GenEvent();
   
   HepMC::GenVertex* Vtx_at = new  HepMC::GenVertex(HepMC::FourVector(CosMuoGen->Vx_at, //[mm]
                                  CosMuoGen->Vy_at, //[mm]
                                  CosMuoGen->Vz_at, //[mm]
                                  CosMuoGen->T0_at)); //[mm]
   //cout << "CosMuoGenSource.cc: Vy_at=" << CosMuoGen->Vy_at << endl;
   HepMC::FourVector p_at(CosMuoGen->Px_at,CosMuoGen->Py_at,CosMuoGen->Pz_at,CosMuoGen->E_at);
   HepMC::GenParticle* Part_at =
     new HepMC::GenParticle(p_at,CosMuoGen->Id_at, 3);//Comment mother particle in
   Vtx_at->add_particle_in(Part_at);
 
 
   //loop here in case of multi muon events (else just one iteration)
   for (unsigned int i=0; i<CosMuoGen->Id_sf.size(); ++i) {
 
     HepMC::FourVector p_sf(CosMuoGen->Px_sf[i],CosMuoGen->Py_sf[i],CosMuoGen->Pz_sf[i],CosMuoGen->E_sf[i]);
     HepMC::GenParticle* Part_sf_in =
       new HepMC::GenParticle(p_sf,CosMuoGen->Id_sf[i], 3); //Comment daughter particle
     Vtx_at->add_particle_out(Part_sf_in);
     
     HepMC::GenVertex* Vtx_sf = new HepMC::GenVertex(HepMC::FourVector(CosMuoGen->Vx_sf[i],                             CosMuoGen->Vy_sf[i], CosMuoGen->Vz_sf[i], CosMuoGen->T0_sf[i])); //[mm]
     HepMC::GenParticle* Part_sf_out =
       new HepMC::GenParticle(p_sf,CosMuoGen->Id_sf[i], 3); //Comment daughter particle
     
     Vtx_sf->add_particle_in(Part_sf_in);
     Vtx_sf->add_particle_out(Part_sf_out);
     
     fEvt->add_vertex(Vtx_sf); //one per muon
 
     HepMC::GenVertex* Vtx_ug = new HepMC::GenVertex(HepMC::FourVector(CosMuoGen->Vx_ug[i],                             CosMuoGen->Vy_ug[i], CosMuoGen->Vz_ug[i], CosMuoGen->T0_ug[i])); //[mm]
     
     HepMC::FourVector p_ug(CosMuoGen->Px_ug[i],CosMuoGen->Py_ug[i],CosMuoGen->Pz_ug[i],CosMuoGen->E_ug[i]);
     HepMC::GenParticle* Part_ug =
       new HepMC::GenParticle(p_ug,CosMuoGen->Id_ug[i], 1);//Final state daughter particle
 
     Vtx_ug->add_particle_in(Part_sf_out);
     Vtx_ug->add_particle_out(Part_ug);
 
     fEvt->add_vertex(Vtx_ug); //one per muon
 
   }
 
   fEvt->add_vertex(Vtx_at);
   fEvt->set_signal_process_vertex(Vtx_at);
 
   fEvt->set_event_number(event());
   fEvt->set_signal_process_id(13);
 
   fEvt->weights().push_back( CosMuoGen->EventWeight ); // just one event weight 
   fEvt->weights().push_back( CosMuoGen->Trials ); // int Trials number (unweighted) 
 
 
   if (cmVerbosity_) fEvt->print();
   
   std::auto_ptr<HepMCProduct> CMProduct(new HepMCProduct());
   CMProduct->addHepMCData( fEvt );
   e.put(CMProduct);
 
   return true;
 }