CMSSQLoopProcessDiscr Class Reference

#include <CMSSQLoopProcessDiscr.h>

Inheritance diagram for CMSSQLoopProcessDiscr:

Public Member Functions
	CMSSQLoopProcessDiscr (double mass, const G4String &name="SQLooper", G4ProcessType type=fUserDefined)
virtual G4double	GetMeanFreePath (const G4Track &, G4double, G4ForceCondition *)
virtual G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &)
virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
void	SetTimeLimit (G4double)
virtual void	StartTracking (G4Track *aTrack)
virtual	~CMSSQLoopProcessDiscr ()

Protected Attributes
G4ParticleChange *	fParticleChange
double	GenMass

Private Member Functions
	CMSSQLoopProcessDiscr (CMSSQLoopProcessDiscr &)
CMSSQLoopProcessDiscr &	operator= (const CMSSQLoopProcessDiscr &right)

Private Attributes
G4double	globaltimeini
G4ThreeVector	posini

Detailed Description

Definition at line 15 of file CMSSQLoopProcessDiscr.h.

Constructor & Destructor Documentation

CMSSQLoopProcessDiscr::CMSSQLoopProcessDiscr	(	double	mass,
		const G4String &	name = "SQLooper",
		G4ProcessType	type = fUserDefined
	)

Definition at line 10 of file CMSSQLoopProcessDiscr.cc.

References fParticleChange, GenMass, and ResonanceBuilder::mass.

    : G4VDiscreteProcess(name, type) {
  fParticleChange = new G4ParticleChange();
  fParticleChange->ClearDebugFlag();
  GenMass = mass;
}

CMSSQLoopProcessDiscr::~CMSSQLoopProcessDiscr ( )

virtual

Definition at line 17 of file CMSSQLoopProcessDiscr.cc.

References fParticleChange.

{ delete fParticleChange; }

CMSSQLoopProcessDiscr::CMSSQLoopProcessDiscr ( CMSSQLoopProcessDiscr &  )

private

Member Function Documentation

G4double CMSSQLoopProcessDiscr::GetMeanFreePath ( const G4Track &  , G4double  , G4ForceCondition *    )

virtual

Definition at line 65 of file CMSSQLoopProcessDiscr.cc.

{ return DBL_MAX; }

CMSSQLoopProcessDiscr& CMSSQLoopProcessDiscr::operator= ( const CMSSQLoopProcessDiscr &  right )

private

G4VParticleChange * CMSSQLoopProcessDiscr::PostStepDoIt ( const G4Track &  track, const G4Step &  step  )

virtual

Definition at line 19 of file CMSSQLoopProcessDiscr.cc.

References CMSAntiSQ::AntiSQ(), fParticleChange, GenMass, globaltimeini, posini, and HiIsolationCommonParameters_cff::track.

                                                                                               {
  G4Track* mytr = const_cast<G4Track*>(&track);
  mytr->SetPosition(posini);
  if (mytr->GetGlobalTime() / ns > 4990)
    edm::LogWarning("CMSSQLoopProcess::AlongStepDoIt")
        << "going to loose the particle because the GlobalTime is getting close to 5000" << std::endl;

  fParticleChange->Clear();
  fParticleChange->Initialize(track);

  //adding secondary antiS
  fParticleChange->SetNumberOfSecondaries(1);
  G4DynamicParticle* replacementParticle =
      new G4DynamicParticle(CMSAntiSQ::AntiSQ(GenMass), track.GetMomentumDirection(), track.GetKineticEnergy());
  fParticleChange->AddSecondary(replacementParticle, globaltimeini);

  //killing original AntiS
  fParticleChange->ProposeTrackStatus(fStopAndKill);

  // note SL: this way of working makes a very long history of the track,
  // which all get saved recursively in SimTracks. If the cross section
  // is too low such that 10's of thousands of iterations are needed, then
  // this becomes too heavy to swallow writing out this history.
  // So if we ever need very small cross sections, then we really need
  // to change this way of working such that we can throw away all original
  // tracks and only save the one that interacted.

  return fParticleChange;
}

G4double CMSSQLoopProcessDiscr::PostStepGetPhysicalInteractionLength ( const G4Track &  track, G4double  previousStepSize, G4ForceCondition *  condition  )

virtual

Definition at line 49 of file CMSSQLoopProcessDiscr.cc.

References Gflash::intLength, funct::pow(), mathSSE::sqrt(), and HiIsolationCommonParameters_cff::track.

                                                                                                  {
  *condition = NotForced;
  G4double intLength =
      DBL_MAX;  //by default the interaction length is super large, only when outside tracker make it 0 to be sure it will do the reset to the original position
  G4Track* mytr = const_cast<G4Track*>(&track);
  if (sqrt(pow(mytr->GetPosition().rho(), 2)) >
      2.45 *
          centimeter) {  //this is an important cut for the looping: if the radius of the particle is largher than 2.45cm its interaction length becomes 0 which means it will get killed
    // updated from 2.5 to 2.45 so that the Sbar does not start to hit the support of the new inner tracker which was added in 2018
    intLength = 0.0;  //0 interaction length means particle will directly interact.
  }
  return intLength;
}

void CMSSQLoopProcessDiscr::SetTimeLimit

(

G4double

)

void CMSSQLoopProcessDiscr::StartTracking ( G4Track *  aTrack )

virtual

Definition at line 67 of file CMSSQLoopProcessDiscr.cc.

References globaltimeini, and posini.

                                                         {
  posini = aTrack->GetPosition();
  globaltimeini = aTrack->GetGlobalTime();
}

Member Data Documentation

G4ParticleChange* CMSSQLoopProcessDiscr::fParticleChange

protected

Definition at line 34 of file CMSSQLoopProcessDiscr.h.

Referenced by CMSSQLoopProcessDiscr(), PostStepDoIt(), and ~CMSSQLoopProcessDiscr().

double CMSSQLoopProcessDiscr::GenMass

protected

Definition at line 35 of file CMSSQLoopProcessDiscr.h.

Referenced by CMSSQLoopProcessDiscr(), and PostStepDoIt().

G4double CMSSQLoopProcessDiscr::globaltimeini

private

Definition at line 39 of file CMSSQLoopProcessDiscr.h.

Referenced by PostStepDoIt(), and StartTracking().

G4ThreeVector CMSSQLoopProcessDiscr::posini

private

Definition at line 38 of file CMSSQLoopProcessDiscr.h.

Referenced by PostStepDoIt(), and StartTracking().