#include <FP420Test.h>

Inheritance diagram for FP420Test:

Public Member Functions
	FP420Test (const edm::ParameterSet &p)

	~FP420Test () override

Public Member Functions inherited from SimWatcher
	SimWatcher ()

virtual	~SimWatcher ()

Public Member Functions inherited from Observer< const BeginOfJob * >
	Observer ()

void	slotForUpdate (const BeginOfJob * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfRun * >
	Observer ()

void	slotForUpdate (const BeginOfRun * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfRun * >
	Observer ()

void	slotForUpdate (const EndOfRun * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfEvent * >
	Observer ()

void	slotForUpdate (const BeginOfEvent * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfTrack * >
	Observer ()

void	slotForUpdate (const BeginOfTrack * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const G4Step * >
	Observer ()

void	slotForUpdate (const G4Step * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfTrack * >
	Observer ()

void	slotForUpdate (const EndOfTrack * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfEvent * >
	Observer ()

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

Private Member Functions
void	detectorLevel (const G4VTouchable , int &, int , G4String *) const

int	detLevels (const G4VTouchable *) const

G4String	detName (const G4VTouchable *, int, int) const

void	update (const BeginOfJob *run) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const EndOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfTrack *trk) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const G4Step *step) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const EndOfTrack *trk) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const EndOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
int	dn0

G4double	entot0

std::string	fDataLabel

std::string	fOutputFile

Float_t	fp420eventarray [1]

TNtuple *	fp420eventntuple

TFile	fp420OutputFile

std::string	fRecreateFile

double	gapBlade

int	iev

int	itrk

G4ThreeVector	lastpo

int	lastTrackID

int	numofpart

int	pn0

double	rinCalo

int	rn0

int	rn00

int	sn0

G4double	SumEnerDeposit

G4double	SumStepc

G4double	SumStepl

Fp420AnalysisHistManager *	TheHistManager

G4double	tracklength0

int	verbosity

int	whichevent

double	z1

double	z2

double	z3

double	z4

double	z420

double	zBlade

double	ZBoundDet

double	zD2

double	zD3

double	ZGapLDet

double	zinCalo

double	zinibeg

double	ZSiDet

double	ZSiPlane

double	ZSiStep

Additional Inherited Members

Detailed Description

Definition at line 153 of file FP420Test.h.

Constructor & Destructor Documentation

FP420Test::FP420Test ( const edm::ParameterSet & p )

Definition at line 51 of file FP420Test.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, HIPAlignmentAlgorithm_cfi::verbosity, FP420Track_cfi::z420, FP420Track_cfi::zD2, and FP420Track_cfi::zD3.

                                              {
   //constructor
   edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("FP420Test");
   verbosity = m_Anal.getParameter<int>("Verbosity");
   //verbosity    = 1;
 
   fDataLabel = m_Anal.getParameter<std::string>("FDataLabel");
   fOutputFile = m_Anal.getParameter<std::string>("FOutputFile");
   fRecreateFile = m_Anal.getParameter<std::string>("FRecreateFile");
 
   if (verbosity > 0) {
     std::cout << "============================================================================" << std::endl;
     std::cout << "FP420Testconstructor :: Initialized as observer" << std::endl;
   }
   // Initialization:
 
   pn0 = 6;
   sn0 = 3;
   rn00 = 7;
 
   z420 = 420000.0;  // mm
   zD2 = 4000.0;     // mm
   zD3 = 8000.0;     // mm
   //
   zBlade = 5.00;
   gapBlade = 1.6;
   double gapSupplane = 1.6;
   ZSiPlane = 2 * zBlade + gapBlade + gapSupplane;
 
   double ZKapton = 0.1;
   ZSiStep = ZSiPlane + ZKapton;
 
   double ZBoundDet = 0.020;
   double ZSiElectr = 0.250;
   double ZCeramDet = 0.500;
   //
   ZSiDet = 0.250;
   ZGapLDet = zBlade / 2 - (ZSiDet + ZSiElectr + ZBoundDet + ZCeramDet / 2);
   //
   //  ZSiStation = 5*(2*(5.+1.6)+0.1)+2*6.+1.0 =  79.5
   double ZSiStation = (pn0 - 1) * (2 * (zBlade + gapBlade) + ZKapton) + 2 * 6. + 0.0;  // =  78.5
   // 11.=e1, 12.=e2 in zzzrectangle.xml
   double eee1 = 11.;
   double eee2 = 12.;
 
   zinibeg = (eee1 - eee2) / 2.;
 
   z1 = zinibeg + (ZSiStation + 10.) / 2 + z420;  // z1 - right after 1st Station
   z2 = z1 + zD2;                                 //z2 - right after middle Station
   z3 = z1 + zD3;                                 //z3 - right after last   Station
   z4 = z1 + 2 * zD3;
   //==================================
 
   fp420eventntuple = new TNtuple("NTfp420event", "NTfp420event", "evt");
 
   whichevent = 0;
 
   //   fDataLabel      = "defaultData";
   //       fOutputFile     = "TheAnlysis.root";
   //       fRecreateFile   = "RECREATE";
 
   TheHistManager = new Fp420AnalysisHistManager(fDataLabel);
 
   if (verbosity > 0) {
     std::cout << "FP420Test constructor :: Initialized Fp420AnalysisHistManager" << std::endl;
   }
 }

FP420Test::~FP420Test ( )

override

Definition at line 119 of file FP420Test.cc.

References gather_cfg::cout, and HIPAlignmentAlgorithm_cfi::verbosity.

                       {
   //  delete UserNtuples;
 
   TFile fp420OutputFile("newntfp420.root", "RECREATE");
   std::cout << "FP420 output root file has been created";
   fp420eventntuple->Write();
   std::cout << ", written";
   fp420OutputFile.Close();
   std::cout << ", closed";
   delete fp420eventntuple;
   std::cout << ", and deleted" << std::endl;
 
   //------->while end
 
   // Write histograms to file
   TheHistManager->WriteToFile(fOutputFile, fRecreateFile);
   if (verbosity > 0) {
     std::cout << std::endl << "FP420Test Destructor  -------->  End of FP420Test : " << std::endl;
   }
 }

Member Function Documentation

void FP420Test::detectorLevel	(	const G4VTouchable *	touch,
		int &	level,
		int *	copyno,
		G4String *	name
	)		const

private

Definition at line 773 of file FP420Test.cc.

References mps_fire::i, cuy::ii, hcalDigis_cfi::level, and MetAnalyzer::pv().

                                                                                                       {
   //Get name and copy numbers
   if (level > 0) {
     for (int ii = 0; ii < level; ii++) {
       int i = level - ii - 1;
       G4VPhysicalVolume* pv = touch->GetVolume(i);
       if (pv != nullptr)
         name[ii] = pv->GetName();
       else
         name[ii] = "Unknown";
       copyno[ii] = touch->GetReplicaNumber(i);
     }
   }
 }

int FP420Test::detLevels ( const G4VTouchable * touch ) const

private

Definition at line 754 of file FP420Test.cc.

                                                         {
   //Return number of levels
   if (touch)
     return ((touch->GetHistoryDepth()) + 1);
   else
     return 0;
 }

G4String FP420Test::detName	(	const G4VTouchable *	touch,
		int	level,
		int	currentlevel
	)		const

private

Definition at line 763 of file FP420Test.cc.

References cuy::ii.

                                                                                         {
   //Go down to current level
   if (level > 0 && level >= currentlevel) {
     int ii = level - currentlevel;
     return touch->GetVolume(ii)->GetName();
   } else {
     return "NotFound";
   }
 }

void FP420Test::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 304 of file FP420Test.cc.

References gather_cfg::cout.

                                             {
   //job
   std::cout << "FP420Test:beggining of job" << std::endl;
   ;
 }

void FP420Test::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 311 of file FP420Test.cc.

References gather_cfg::cout.

                                             {
   //run
 
   std::cout << std::endl << "FP420Test:: Begining of Run" << std::endl;
 }

void FP420Test::update ( const EndOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfRun * >.

Definition at line 317 of file FP420Test.cc.

317 { ; }

void FP420Test::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 320 of file FP420Test.cc.

References gather_cfg::cout, and HIPAlignmentAlgorithm_cfi::verbosity.

                                               {
   iev = (*evt)()->GetEventID();
   if (verbosity > 0) {
     std::cout << "FP420Test:update Event number = " << iev << std::endl;
   }
   whichevent++;
 }

void FP420Test::update ( const BeginOfTrack * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 329 of file FP420Test.cc.

References gather_cfg::cout, and HIPAlignmentAlgorithm_cfi::verbosity.

                                               {
   itrk = (*trk)()->GetTrackID();
   if (verbosity > 1) {
     std::cout << "FP420Test:update BeginOfTrack number = " << itrk << std::endl;
   }
   if (itrk == 1) {
     SumEnerDeposit = 0.;
     numofpart = 0;
     SumStepl = 0.;
     SumStepc = 0.;
     tracklength0 = 0.;
   }
 }

void FP420Test::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 444 of file FP420Test.cc.

References gather_cfg::cout, mps_fire::i, objects.autophobj::particleType, and HIPAlignmentAlgorithm_cfi::verbosity.

                                           {
   // ==========================================================================
 
   if (verbosity > 2) {
     G4int stepnumber = aStep->GetTrack()->GetCurrentStepNumber();
     std::cout << "FP420Test:update Step number = " << stepnumber << std::endl;
   }
   // track on aStep:                                                                                         !
   G4Track* theTrack = aStep->GetTrack();
   TrackInformation* trkInfo = dynamic_cast<TrackInformation*>(theTrack->GetUserInformation());
   if (trkInfo == nullptr) {
     std::cout << "FP420Test on aStep: No trk info !!!! abort " << std::endl;
   }
   G4int id = theTrack->GetTrackID();
   G4String particleType = theTrack->GetDefinition()->GetParticleName();  //   !!!
   G4int parentID = theTrack->GetParentID();                              //   !!!
   G4TrackStatus trackstatus = theTrack->GetTrackStatus();                //   !!!
   G4double tracklength = theTrack->GetTrackLength();                     // Accumulated track length
   G4ThreeVector trackmom = theTrack->GetMomentum();
   G4double entot = theTrack->GetTotalEnergy();  //   !!! deposited on step
   G4int curstepnumber = theTrack->GetCurrentStepNumber();
   // const G4ThreeVector&   vert_pos       = theTrack->GetVertexPosition(); // vertex ,where this track was created
   // const G4ThreeVector&   vert_mom       = theTrack->GetVertexMomentumDirection();
 
   //   double costheta =vert_mom.z()/sqrt(vert_mom.x()*vert_mom.x()+vert_mom.y()*vert_mom.y()+vert_mom.z()*vert_mom.z());
   //   double theta = acos(min(max(costheta,double(-1.)),double(1.)));
   //  float eta = -log(tan(theta/2));
   //   double phi = -1000.;
   //   if (vert_mom.x() != 0) phi = atan2(vert_mom.y(),vert_mom.x());
   //   if (phi < 0.) phi += twopi;
   //   double phigrad = phi*360./twopi;
 
   //G4double       trackvel       = theTrack->GetVelocity();
 
   //std::cout << " trackstatus= " << trackstatus << " entot= " << entot  << std::endl;
 
   // step points:                                                                                         !
   G4double stepl = aStep->GetStepLength();
   G4double EnerDeposit = aStep->GetTotalEnergyDeposit();
   // pointers:                                                                                         !
   //G4VPhysicalVolume*  physvol       = theTrack->GetVolume();
   //G4VPhysicalVolume*  nextphysvol   = theTrack->GetNextVolume();
   //G4Material*       materialtr     = theTrack->GetMaterial();
   //G4Material*       nextmaterialtr = theTrack->GetNextMaterial();
 
   // preStep
   G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
   const G4ThreeVector& preposition = preStepPoint->GetPosition();
   G4ThreeVector prelocalpoint = theTrack->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(preposition);
   G4VPhysicalVolume* currentPV = preStepPoint->GetPhysicalVolume();
   const G4String& prename = currentPV->GetName();
 
   const G4VTouchable* pre_touch = preStepPoint->GetTouchable();
   int pre_levels = detLevels(pre_touch);
 
   G4String name1[20];
   int copyno1[20];
   for (int i = 0; i < 20; ++i) {
     name1[i] = "";
     copyno1[i] = 0;
   }
   if (pre_levels > 0) {
     detectorLevel(pre_touch, pre_levels, copyno1, name1);
   }
 
   //  G4LogicalVolume*   lv            = currentPV->GetLogicalVolume();
   //  G4Material*       mat           = lv->GetMaterial();
   //  std::string prenameVolume;
   //  prenameVolume.assign(prename,0,20);
 
   //   G4double         prebeta          = preStepPoint->GetBeta();
   //   G4double         precharge        = preStepPoint->GetCharge();
   //  G4double          prerad           = mat->GetRadlen();
 
   //  std::cout << " EneryDeposited = " << EnerDeposit << std::endl;
   //  std::cout << " prevolume = "      << prename << std::endl;
   //  std::cout << " preposition = "    << preposition << std::endl;
   /*
   // postStep
   G4StepPoint*      postStepPoint  = aStep->GetPostStepPoint();   
   G4ThreeVector     posposition    = postStepPoint->GetPosition();  
   G4ThreeVector     poslocalpoint  = theTrack->GetTouchable()->GetHistory()->
                                            GetTopTransform().TransformPoint(posposition);
   G4VPhysicalVolume* poscurrentPV      = postStepPoint->GetPhysicalVolume();
   G4String         posname        = poscurrentPV->GetName();
 //  G4LogicalVolume*   poslv             = poscurrentPV->GetLogicalVolume();
 //  G4Material*       posmat            = poslv->GetMaterial();
 //  std::string posnameVolume;
 //  posnameVolume.assign(posname,0,20);
 
 #ifdef ddebug
      std::cout << "============posStep: information:============" << std::endl;
      std::cout << " posposition = "    << posposition
           << " poslocalpoint = "  << poslocalpoint
           << " posvolume = "      << posname
        //          << " posnameVolume = "  << posnameVolume 
           << std::endl;
 
      std::cout << " ==========================================" << std::endl;
 #endif
 
 
 */
 
   //      //
   //      //
 
   // 24.01.2006:
   // tr     :   id    parentID   trackstatus   tracklength   curstepnumber  entot  vert_pos
   // st     :   stepl EnerDeposit
   // prestep:   preposition   prevolume = SBSTm SIDETL SIDETR       name= SISTATION  copy= 1,2,3    name= SIPLANE  copy= 1..5..10
 
   // gen_track:
   //  id=1    parentID=1   trackstatus=0,2   tracklength(accumulated)  curstepnumber   entot  vert_pos
   if (id == 1) {
     // on 1-st step:
     if (curstepnumber == 1) {
       entot0 = entot;
       //UserNtuples->fillg519(entot0,1.);
     }
 
     // on every step:
 
     // for Copper:
     if (prename == "SBST") {
       SumStepc += stepl;
       // =========
     }
     // for ststeel:
     //   if(prename == "SBSTs") {
     if (prename == "SBSTs") {
       SumStepl += stepl;
       // =========
     }
     // =========
     // =========
 
     // exclude last track point if it is in SD (MI was started their)
     if (trackstatus != 2) {
       // for SD:   Si Det.:   SISTATION:SIPLANE:(SIDETL+BOUNDDET        +SIDETR + CERAMDET)
       if (prename == "SIDETL" || prename == "SIDETR") {
         if (prename == "SIDETL") {
           //UserNtuples->fillg569(EnerDeposit,1.);
         }
         if (prename == "SIDETR") {
           //UserNtuples->fillg570(EnerDeposit,1.);
         }
 
         //     double numbcont = 10.*(copyno1[2]-1)+copyno1[3];
 
         // =========
         //     double   xx    = preposition.x();
         //     double   yy    = preposition.y();
         //     double   zz    = preposition.z();
         // =========
         //UserNtuples->fillg580(theta,1.);
         //UserNtuples->fillg07(phigrad,1.);
         //       double xPrimAtZhit = vert_pos.x() + (zz-vert_pos.z())*tan(theta)*cos(phi);
         //       double yPrimAtZhit = vert_pos.y() + (zz-vert_pos.z())*tan(theta)*sin(phi);
         // =========
         //     double  dx = xPrimAtZhit - xx;
         //     double  dy = yPrimAtZhit - yy;
         // =========
 
         //     //UserNtuples->fillp212(numbcont,dx,1.);
         //     //UserNtuples->fillp213(numbcont,dy,1.);
         // =========
 
         // last step of current SD Volume:
         G4String posname = aStep->GetPostStepPoint()->GetPhysicalVolume()->GetName();
         if ((prename == "SIDETL" && posname != "SIDETL") || (prename == "SIDETR" && posname != "SIDETR")) {
           if (name1[2] == "SISTATION") {
             //UserNtuples->fillg539(copyno1[2],1.);
           }
           if (name1[3] == "SIPLANE") {
             //UserNtuples->fillg540(copyno1[3],1.);
           }
 
           if (prename == "SIDETL") {
             //UserNtuples->fillg541(EnerDeposit,1.);
             //UserNtuples->fillg561(numbcont,1.);
             if (copyno1[2] < 2) {
               //UserNtuples->fillg571(dx,1.);
             } else if (copyno1[2] < 3) {
               //UserNtuples->fillg563(dx,1.);
               if (copyno1[3] < 2) {
               } else if (copyno1[3] < 3) {
                 //UserNtuples->fillg572(dx,1.);
               } else if (copyno1[3] < 4) {
                 //UserNtuples->fillg573(dx,1.);
               } else if (copyno1[3] < 5) {
                 //UserNtuples->fillg574(dx,1.);
               } else if (copyno1[3] < 6) {
                 //UserNtuples->fillg575(dx,1.);
               } else if (copyno1[3] < 7) {
                 //UserNtuples->fillg576(dx,1.);
               } else if (copyno1[3] < 8) {
                 //UserNtuples->fillg577(dx,1.);
               } else if (copyno1[3] < 9) {
                 //UserNtuples->fillg578(dx,1.);
               } else if (copyno1[3] < 10) {
                 //UserNtuples->fillg579(dx,1.);
               }
             } else if (copyno1[2] < 4) {
               //UserNtuples->fillg565(dx,1.);
             } else if (copyno1[2] < 5) {
               //UserNtuples->fillg567(dx,1.);
             }
           }
           if (prename == "SIDETR") {
             //UserNtuples->fillg542(EnerDeposit,1.);
             //UserNtuples->fillg562(numbcont,1.);
             if (copyno1[2] < 2) {
               //UserNtuples->fillg581(dy,1.);
             } else if (copyno1[2] < 3) {
               //UserNtuples->fillg564(dy,1.);
               if (copyno1[3] < 2) {
               } else if (copyno1[3] < 3) {
                 //UserNtuples->fillg582(dy,1.);
               } else if (copyno1[3] < 4) {
                 //UserNtuples->fillg583(dy,1.);
               } else if (copyno1[3] < 5) {
                 //UserNtuples->fillg584(dy,1.);
               } else if (copyno1[3] < 6) {
                 //UserNtuples->fillg585(dy,1.);
               } else if (copyno1[3] < 7) {
                 //UserNtuples->fillg586(dy,1.);
               } else if (copyno1[3] < 8) {
                 //UserNtuples->fillg587(dy,1.);
               } else if (copyno1[3] < 9) {
                 //UserNtuples->fillg588(dy,1.);
               } else if (copyno1[3] < 10) {
                 //UserNtuples->fillg589(dy,1.);
               }
             } else if (copyno1[2] < 4) {
               //UserNtuples->fillg566(dy,1.);
             } else if (copyno1[2] < 5) {
               //UserNtuples->fillg568(dy,1.);
             }
           }
         }
       }
       // end of prenames SIDETL // SIDETR
     }
     // end of trackstatus != 2
 
     // deposition of energy on steps along primary track
     //UserNtuples->fillg500(EnerDeposit,1.);
     // collect sum deposited energy on all steps along primary track
     SumEnerDeposit += EnerDeposit;
     // position of step for primary track:
     //UserNtuples->fillg501(preposition.x(),1.);
     //UserNtuples->fillg502(preposition.y(),1.);
     //UserNtuples->fillg503(preposition.z(),1.);
     //UserNtuples->fillg504(preposition.x(),EnerDeposit);
     //UserNtuples->fillg505(preposition.y(),EnerDeposit);
     //UserNtuples->fillg506(preposition.z(),EnerDeposit);
     // 2D step coordinates weighted by energy deposited on step
     //      //UserNtuples->fillp201(preposition.x(),preposition.y(),EnerDeposit);
     //      //UserNtuples->fillp202(preposition.x(),preposition.z(),EnerDeposit);
     //      //UserNtuples->fillp203(preposition.y(),preposition.z(),EnerDeposit);
     //UserNtuples->filld204(preposition.x(),preposition.y(),EnerDeposit);
     //UserNtuples->filld205(preposition.x(),preposition.z(),EnerDeposit);
     //UserNtuples->filld206(preposition.y(),preposition.z(),EnerDeposit);
     //UserNtuples->filld223(preposition.y(),preposition.z(),EnerDeposit);
     // last step of primary track
     if (trackstatus == 2) {
       // primary track length
       //      //UserNtuples->fillg508(tracklength,1.);
       tracklength0 = tracklength;
       // how many steps primary track consist
       //UserNtuples->fillg509(curstepnumber,1.);
       // tot. energy of primary track at the end of trajectory(before disappeare)
       //UserNtuples->fillg510((entot0-entot),1.);
       //UserNtuples->fillg520((entot0-entot),1.);
     }
   }
   // end of primary track !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
   if (parentID == 1 && curstepnumber == 1) {
     // particles deposit their energy along primary track
     numofpart += 1;
     // energy of radiated particle
     //UserNtuples->fillg511(entot,1.);
     // coordinates  of radiated particle
     //UserNtuples->fillg512(vert_pos.x(),1.);
     //UserNtuples->fillg513(vert_pos.y(),1.);
     //UserNtuples->fillg514(vert_pos.z(),1.);
     //UserNtuples->fillg515(vert_pos.x(),entot);
     //UserNtuples->fillg516(vert_pos.y(),entot);
     //UserNtuples->fillg517(vert_pos.z(),entot);
 
     //UserNtuples->filld214(vert_pos.x(),vert_pos.y(),1.);
     //UserNtuples->filld215(vert_pos.x(),vert_pos.z(),1.);
     //UserNtuples->filld216(vert_pos.y(),vert_pos.z(),1.);
     //UserNtuples->filld219(vert_pos.y(),vert_pos.z(),1.);
     //UserNtuples->filld224(vert_pos.y(),vert_pos.z(),1.);
     if (prename == "SBST") {
       //UserNtuples->filld225(vert_pos.y(),vert_pos.z(),1.);
     }
     if (prename == "SBSTs") {
       //UserNtuples->filld226(vert_pos.y(),vert_pos.z(),1.);
     }
   }
 
   // ==========================================================================
 }

void FP420Test::update ( const EndOfTrack * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack * >.

Definition at line 344 of file FP420Test.cc.

References gather_cfg::cout, and HIPAlignmentAlgorithm_cfi::verbosity.

                                             {
   itrk = (*trk)()->GetTrackID();
   if (verbosity > 1) {
     std::cout << "FP420Test:update EndOfTrack number = " << itrk << std::endl;
   }
   if (itrk == 1) {
     G4double tracklength = (*trk)()->GetTrackLength();  // Accumulated track length
 
     TheHistManager->GetHisto("SumEDep")->Fill(SumEnerDeposit);
     TheHistManager->GetHisto("TrackL")->Fill(tracklength);
 
     // direction !!!
     G4ThreeVector vert_mom = (*trk)()->GetVertexMomentumDirection();
     G4ThreeVector vert_pos = (*trk)()->GetVertexPosition();  // vertex ,where this track was created
 
     //    float eta = 0.5 * log( (1.+vert_mom.z()) / (1.-vert_mom.z()) );
     //    float phi = atan2(vert_mom.y(),vert_mom.x());
     //    if (phi < 0.) phi += twopi;
     //    float phigrad = phi*180./pi;
 
     //      float XV = vert_pos.x(); // mm
     //      float YV = vert_pos.y(); // mm
     //UserNtuples->fillg543(phigrad,1.);
     //UserNtuples->fillp203(phigrad,SumStepl,1.);
     //UserNtuples->fillg544(XV,1.);
     //UserNtuples->fillp201(XV,SumStepl,1.);
     //UserNtuples->fillg545(YV,1.);
     //UserNtuples->fillp202(YV,SumStepl,1.);
 
     //UserNtuples->fillg524(eta,1.);
 
     //UserNtuples->fillg534(SumStepl,1.);
     //UserNtuples->fillg535(eta,SumStepl);
     //UserNtuples->fillp207(eta,SumStepl,1.);
     //UserNtuples->filld217(eta,SumStepl,1.);
     //UserNtuples->filld220(phigrad,SumStepl,1.);
     //UserNtuples->filld221(XV,SumStepl,1.);
     //UserNtuples->filld222(YV,SumStepl,1.);
     //UserNtuples->fillg537(SumStepc,1.);
     //UserNtuples->fillg84(SumStepl,1.);
 
     // MI = (multiple interactions):
     if (tracklength < z4) {
       //        //UserNtuples->fillp208(eta,tracklength,1.);
       //UserNtuples->filld218(eta,tracklength,1.);
       //UserNtuples->fillg538(SumStepc,1.);
       //UserNtuples->fillg85(SumStepl,1.);
     }
 
     // last step information
     const G4Step* aStep = (*trk)()->GetStep();
     //   G4int csn = (*trk)()->GetCurrentStepNumber();
     //   G4double sl = (*trk)()->GetStepLength();
     // preStep
     G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
     lastpo = preStepPoint->GetPosition();
 
     // Analysis:
     if (lastpo.z() < z1 && lastpo.perp() < 100.) {
       //UserNtuples->fillg525(eta,1.);
       //UserNtuples->fillg546(XV,1.);
       //UserNtuples->fillg551(YV,1.);
       //UserNtuples->fillg556(phigrad,1.);
     }
     if ((lastpo.z() > z1 && lastpo.z() < z2) && lastpo.perp() < 100.) {
       //UserNtuples->fillg526(eta,1.);
       //UserNtuples->fillg547(XV,1.);
       //UserNtuples->fillg552(YV,1.);
       //UserNtuples->fillg557(phigrad,1.);
     }
     if (lastpo.z() < z2 && lastpo.perp() < 100.) {
       //UserNtuples->fillg527(eta,1.);
       //UserNtuples->fillg548(XV,1.);
       //UserNtuples->fillg553(YV,1.);
       //UserNtuples->fillg558(phigrad,1.);
       //UserNtuples->fillg521(lastpo.x(),1.);
       //UserNtuples->fillg522(lastpo.y(),1.);
       //UserNtuples->fillg523(lastpo.z(),1.);
     }
     if (lastpo.z() < z3 && lastpo.perp() < 100.) {
       //UserNtuples->fillg528(eta,1.);
       //UserNtuples->fillg549(XV,1.);
       //UserNtuples->fillg554(YV,1.);
       //UserNtuples->fillg559(phigrad,1.);
     }
     if (lastpo.z() < z4 && lastpo.perp() < 100.) {
       //UserNtuples->fillg529(eta,1.);
       //UserNtuples->fillg550(XV,1.);
       //UserNtuples->fillg555(YV,1.);
       //UserNtuples->fillg560(phigrad,1.);
       //UserNtuples->fillg531(lastpo.x(),1.);
       //UserNtuples->fillg532(lastpo.y(),1.);
       //UserNtuples->fillg533(lastpo.z(),1.);
     }
   }
 }

void FP420Test::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 790 of file FP420Test.cc.

References gather_cfg::cout, PVValHelper::eta, objects.autophobj::float, FP420G4Hit::getEnergyLoss(), FP420G4Hit::getEntry(), FP420G4Hit::getEntryLocalP(), FP420G4Hit::getExitLocalP(), FP420G4Hit::getTrackID(), FP420G4Hit::getUnitID(), mps_fire::i, cuy::ii, cmsBatch::log, npart, ntfp420_evt, FP420NumberingScheme::packFP420Index(), FP420NumberingScheme::packMYIndex(), pi, FP420NumberingScheme::realzside(), funct::tan(), FP420NumberingScheme::unpackCopyIndex(), FP420NumberingScheme::unpackFP420Index(), FP420NumberingScheme::unpackLayerIndex(), FP420NumberingScheme::unpackOrientation(), HIPAlignmentAlgorithm_cfi::verbosity, FP420Track_cfi::z420, FP420Track_cfi::zD2, FP420Track_cfi::zD3, ecaldqm::zside(), and geometryCSVtoXML::zz.

                                             {
   // ==========================================================================
 
   if (verbosity > 1) {
     iev = (*evt)()->GetEventID();
     std::cout << "FP420Test:update EndOfEvent = " << iev << std::endl;
   }
   // Fill-in ntuple
   fp420eventarray[ntfp420_evt] = (float)whichevent;
 
   //
   int trackID = 0;
   G4PrimaryParticle* thePrim = nullptr;
 
   // prim.vertex:
   G4int nvertex = (*evt)()->GetNumberOfPrimaryVertex();
   if (nvertex != 1)
     std::cout << "FP420Test: My warning: NumberOfPrimaryVertex != 1  -->  = " << nvertex << std::endl;
 
   for (int i = 0; i < nvertex; i++) {
     G4PrimaryVertex* avertex = (*evt)()->GetPrimaryVertex(i);
     if (avertex == nullptr)
       std::cout << "FP420Test  End Of Event ERR: pointer to vertex = 0" << std::endl;
     G4int npart = avertex->GetNumberOfParticle();
     if (npart != 1)
       std::cout << "FP420Test: My warning: NumberOfPrimaryPart != 1  -->  = " << npart << std::endl;
     if (npart == 0)
       std::cout << "FP420Test End Of Event ERR: no NumberOfParticle" << std::endl;
 
     // find just primary track:                                                             track pointer: thePrim
     if (thePrim == nullptr)
       thePrim = avertex->GetPrimary(trackID);
 
     if (thePrim != nullptr) {
       // primary vertex:
       G4double vx = 0., vy = 0., vz = 0.;
       vx = avertex->GetX0();
       vy = avertex->GetY0();
       vz = avertex->GetZ0();
       //UserNtuples->fillh01(vx);
       //UserNtuples->fillh02(vy);
       //UserNtuples->fillh03(vz);
       TheHistManager->GetHisto("VtxX")->Fill(vx);
       TheHistManager->GetHisto("VtxY")->Fill(vy);
       TheHistManager->GetHisto("VtxZ")->Fill(vz);
     }
   }
   // prim.vertex loop end
 
   //=========================== thePrim != 0 ================================================================================
   if (thePrim != nullptr) {
     //      inline G4ParticleDefinition * GetG4code() const
     //      inline G4PrimaryParticle * GetNext() const
     //      inline G4PrimaryParticle * GetDaughter() const
 
     // prim.vertex
     //int ivert = 0 ;
     //G4PrimaryVertex* avertex = (*evt)()->GetPrimaryVertex(ivert);
     //G4double vx=0.,vy=0.,vz=0.;
     //vx = avertex->GetX0();
     //vy = avertex->GetY0();
     //vz = avertex->GetZ0();
 
     //
     // number of secondary particles deposited their energy along primary track
     //UserNtuples->fillg518(numofpart,1.);
     if (lastpo.z() < z4 && lastpo.perp() < 100.) {
       //UserNtuples->fillg536(numofpart,1.);
     }
     //
 
     // direction !!!
     G4ThreeVector mom = thePrim->GetMomentum();
 
     double phi = atan2(mom.y(), mom.x());
     if (phi < 0.)
       phi += twopi;
     double phigrad = phi * 180. / pi;
 
     double th = mom.theta();
     double eta = -log(tan(th / 2));
     // works OK:
     //      double  costheta =mom.z()/sqrt(mom.x()*mom.x()+mom.y()*mom.y()+mom.z()*mom.z());
     //      double theta = acos(min(max(costheta,double(-1.)),double(1.)));
 
     //UserNtuples->fillh04(eta);
     //UserNtuples->fillh05(phigrad);
     //UserNtuples->fillh06(th);
 
     //UserNtuples->fillg13(lastpo.x(),1.);
     //UserNtuples->fillg14(lastpo.y(),1.);
     //UserNtuples->fillg15(lastpo.z(),1.);
 
     TheHistManager->GetHisto("PrimaryEta")->Fill(eta);
     TheHistManager->GetHisto("PrimaryPhigrad")->Fill(phigrad);
     TheHistManager->GetHisto("PrimaryTh")->Fill(th * 180. / pi);
 
     TheHistManager->GetHisto("PrimaryLastpoZ")->Fill(lastpo.z());
     if (lastpo.z() < z4) {
       TheHistManager->GetHisto("PrimaryLastpoX")->Fill(lastpo.x());
       TheHistManager->GetHisto("PrimaryLastpoY")->Fill(lastpo.y());
     }
     if (numofpart > 4) {
       TheHistManager->GetHisto("XLastpoNumofpart")->Fill(lastpo.x());
       TheHistManager->GetHisto("YLastpoNumofpart")->Fill(lastpo.y());
     }
 
     // ==========================================================================
 
     // hit map for FP420
     // ==================================
 
     map<int, float, less<int> > themap;
     map<int, float, less<int> > themap1;
 
     map<int, float, less<int> > themapxy;
     map<int, float, less<int> > themapz;
     // access to the G4 hit collections:  -----> this work OK:
 
     //  edm::LogInfo("FP420Test") << "1";
     G4HCofThisEvent* allHC = (*evt)()->GetHCofThisEvent();
     //  edm::LogInfo("FP420Test") << "2";
     if (verbosity > 0) {
       std::cout << "FP420Test:  accessed all HC" << std::endl;
       ;
     }
     int CAFIid = G4SDManager::GetSDMpointer()->GetCollectionID("FP420SI");
     // edm::LogInfo("FP420Test") << "3";
     // std::cout << " CAFIid = " << CAFIid << std::endl;;
 
     FP420G4HitCollection* theCAFI = (FP420G4HitCollection*)allHC->GetHC(CAFIid);
     //  CaloG4HitCollection* theCAFI = (CaloG4HitCollection*) allHC->GetHC(CAFIid);
     if (verbosity > 0) {
       //std::cout << "FP420Test: theCAFI = " << theCAFI << std::endl;
       std::cout << "FP420Test: theCAFI->entries = " << theCAFI->entries() << std::endl;
     }
     // edm::LogInfo("FP420Test") << "theCAFI->entries="<< theCAFI->entries();
     TheHistManager->GetHisto("NumberOfHits")->Fill(theCAFI->entries());
 
     // access to the G4 hit collections ----> variant 2: give 0 hits
     //  FP420G4HitCollection *   theCAFI;
     //  theCAFI = new FP420G4HitCollection();
     // ==========================================================================
     //   Silicon Hit collection start
     //0) if particle goes into flat beam pipe below detector:
     int varia;  // = 0 -all; =1 - MI; =2 - noMI
     //                      Select MI or noMI over all 3 stations
     // 1)MI:
     //     if particle goes through window into detector:
     // lastpoint of track in lateral dir. outside the detector and inside in z
     // lastpoint of track in lateral dir. outside the detector and inside in z
     // for all except zzzmarta.xml
     //    if(  lastpo.z()<z4  ||  abs(lastpo.x())> 5. || lastpo.y()< 10.2 || lastpo.y()>30.2   ) {
     // for zzzmarta.xml
     //   if(  lastpo.z()<z4  ||  abs(lastpo.x())> 10. || lastpo.y()< 3.2 || lastpo.y()>43.2   ) {
     if (lastpo.z() < z4) {
       //  if(  lastpo.z()<z4 && lastpo.perp()< 100. ) {
       //  if(lastpo.z()<z4  || lastpo.perp()> 45.) {
       //UserNtuples->fillg66(theCAFI->entries(),1.);
       varia = 1;
     } else {
       // 2)   no MI start in detector:
       //UserNtuples->fillg31(numofpart,1.);
       //UserNtuples->fillg67(theCAFI->entries(),1.);
       varia = 2;
     }  // no MI end:
 
     for (int j = 0; j < theCAFI->entries(); j++) {
       FP420G4Hit* aHit = (*theCAFI)[j];
       G4ThreeVector hitPoint = aHit->getEntry();
       double zz = hitPoint.z();
       TheHistManager->GetHisto("zHits")->Fill(zz);
       if (tracklength0 > 8300.)
         TheHistManager->GetHisto("zHitsTrLoLe")->Fill(zz);
     }
     // varia = 0;
     //     if( varia == 0) {
     if (varia == 2) {
       // .............
       // number of hits < 50
       //    if(theCAFI->entries() <50) {
       //    if(theCAFI->entries() > 0) {
       //    if(theCAFI->entries() > -1) {
       // .............
       int nhit11 = 0, nhit12 = 0, nhit13 = 0;
       double totallosenergy = 0.;
       for (int j = 0; j < theCAFI->entries(); j++) {
         FP420G4Hit* aHit = (*theCAFI)[j];
 
         G4ThreeVector hitEntryLocalPoint = aHit->getEntryLocalP();
         G4ThreeVector hitExitLocalPoint = aHit->getExitLocalP();
         G4ThreeVector hitPoint = aHit->getEntry();
         //    double  elmenergy =  aHit->getEM();
         //    double  hadrenergy =  aHit->getHadr();
         //    double incidentEnergyHit  = aHit->getIncidentEnergy();
         int trackIDhit = aHit->getTrackID();
         unsigned int unitID = aHit->getUnitID();
         //    double   timeslice = aHit->getTimeSlice();
         //    int     timesliceID = aHit->getTimeSliceID();
         //    double  depenergy = aHit->getEnergyDeposit();
         //    float   pabs = aHit->getPabs();
         //    float   tof = aHit->getTof();
         double losenergy = aHit->getEnergyLoss();
         //    int   particletype = aHit->getParticleType();
         //    float thetaEntry = aHit->getThetaAtEntry();
         //    float phiEntry = aHit->getPhiAtEntry();
         //    float xEntry = aHit->getX();
         //    float yEntry = aHit->getY();
         //    float zEntry = aHit->getZ();
         //    int  parentID = aHit->getParentId();
         //    float vxget = aHit->getVx();
         //    float vyget = aHit->getVy();
         //    float vzget = aHit->getVz();
 
         //    double th_hit    = hitPoint.theta();
         //    double eta_hit = -log(tan(th_hit/2));
         //    double phi_hit   = hitPoint.phi();
         //    if (phi_hit < 0.) phi_hit += twopi;
         //    double phigrad_hit = phi_hit*180./pi;
         //UserNtuples->fillg60(eta_hit,losenergy);
         //UserNtuples->fillg61(eta_hit,1.);
         //UserNtuples->fillg62(phigrad_hit,losenergy);
         //UserNtuples->fillg63(phigrad_hit,1.);
 
         //    double   xx    = hitPoint.x();
         //    double   yy    = hitPoint.y();
         double zz = hitPoint.z();
 
         TheHistManager->GetHisto("zHitsnoMI")->Fill(zz);
 
         if (verbosity > 2) {
           std::cout << "FP420Test:zHits = " << zz << std::endl;
         }
         //   double   rr    = hitPoint.perp();
         /*
       if(aHit->getTrackID() == 1) {
       emu += aHit->getEnergyDeposit();} else {
       erest += aHit->getEnergyDeposit();}
     */
 
         //collect lost in Si en.of hits in every plane and put it into themap[]
         //UserNtuples->fillg30(losenergy,1.);
         themap[unitID] += losenergy;
         totallosenergy += losenergy;
 
         int det, zside, sector, zmodule;
         //    CaloNumberingPacker::unpackCastorIndex(unitID, det, zside, sector, zmodule);
         FP420NumberingScheme::unpackFP420Index(unitID, det, zside, sector, zmodule);
         int justlayer = FP420NumberingScheme::unpackLayerIndex(rn00, zside);  // 1,2
         if (justlayer < 1 || justlayer > 2) {
           std::cout << "FP420Test:WRONG  justlayer= " << justlayer << std::endl;
         }
         // zside=1,2 ; zmodule=1,10 ; sector=1,3
         //UserNtuples->fillg44(float(sector),1.);
         //UserNtuples->fillg45(float(zmodule),1.);
         //UserNtuples->fillg46(float(zside),1.);
         //      int sScale = 20;
         // intindex is a continues numbering of FP420
         //int zScale = 2; unsigned int intindex = sScale*(sector - 1)+zScale*(zmodule - 1)+zside; //intindex=1-30:X,Y,X,Y,X,Y...
         // int zScale = 10;   unsigned int intindex = sScale*(sector - 1)+zScale*(zside - 1)+zmodule; //intindex=1-30:XXXXXXXXXX,YYYYYYYYYY,...
         //UserNtuples->fillg40(float(intindex),1.);
         //UserNtuples->fillg48(float(intindex),losenergy);
         //
         //=======================================
         //   G4ThreeVector middle = (hitExitLocalPoint+hitEntryLocalPoint)/2.;
         G4ThreeVector middle = (hitExitLocalPoint - hitEntryLocalPoint) / 2.;
         themapz[unitID] = hitPoint.z() + fabs(middle.z());
         if (verbosity > 2) {
           std::cout << "1111111111111111111111111111111111111111111111111111111111111111111111111 " << std::endl;
           std::cout << "FP420Test: det, zside, sector, zmodule = " << det << zside << sector << zmodule << std::endl;
           std::cout << "FP420Test: justlayer = " << justlayer << std::endl;
           std::cout << "FP420Test: hitExitLocalPoint = " << hitExitLocalPoint << std::endl;
           std::cout << "FP420Test: hitEntryLocalPoint = " << hitEntryLocalPoint << std::endl;
           std::cout << "FP420Test:  middle= " << middle << std::endl;
           std::cout << "FP420Test:  hitPoint.z()-419000.= " << hitPoint.z() - 419000. << std::endl;
 
           std::cout << "FP420Test:zHits-419000. = " << themapz[unitID] - 419000. << std::endl;
         }
         //=======================================
         // Y
         if (zside == 1) {
           //UserNtuples->fillg24(losenergy,1.);
           if (losenergy > 0.00003) {
             themap1[unitID] += 1.;
           }
         }
         //X
         if (zside == 2) {
           //UserNtuples->fillg25(losenergy,1.);
           if (losenergy > 0.00005) {
             themap1[unitID] += 1.;
           }
         }
         //     }
         //
         if (sector == 1) {
           nhit11 += 1;
           //UserNtuples->fillg33(rr,1.);
           //UserNtuples->fillg11(yy,1.);
         }
         if (sector == 2) {
           nhit12 += 1;
           //UserNtuples->fillg34(rr,1.);
           //UserNtuples->fillg86(yy,1.);
         }
         if (sector == 3) {
           nhit13 += 1;
           //UserNtuples->fillg35(rr,1.);
           //UserNtuples->fillg87(yy,1.);
         }
         //UserNtuples->fillg10(xx,1.);
         //UserNtuples->fillg12(zz,1.);
         //UserNtuples->fillg32(rr,1.);
 
         // =========
         //    double xPrimAtZhit = vx + (zz-vz)*tan(th)*cos(phi);
         //    double yPrimAtZhit = vy + (zz-vz)*tan(th)*sin(phi);
 
         //       double  dx = xPrimAtZhit - xx;
         //       double  dy = yPrimAtZhit - yy;
 
         //                      Select SD hits
         //    if(rr<120.) {
         //                      Select MI or noMI over all 3 stations
         if (lastpo.z() < z4 && lastpo.perp() < 120.) {
           // MIonly:
           //UserNtuples->fillg16(lastpo.z(),1.);
           //UserNtuples->fillg18(zz,1.);
           //                                                                     Station I
           if (zz < z2) {
             //UserNtuples->fillg54(dx,1.);
             //UserNtuples->fillg55(dy,1.);
           }
           //                                                                     Station II
           if (zz < z3 && zz > z2) {
             //UserNtuples->fillg50(dx,1.);
             //UserNtuples->fillg51(dy,1.);
           }
           //                                                                     Station III
           if (zz < z4 && zz > z3) {
             //UserNtuples->fillg64(dx,1.);
             //UserNtuples->fillg65(dy,1.);
             //UserNtuples->filld209(xx,yy,1.);
           }
         } else {
           // no MIonly:
           //UserNtuples->fillg17(lastpo.z(),1.);
           //UserNtuples->fillg19(zz,1.);
           //UserNtuples->fillg74(incidentEnergyHit,1.);
           //UserNtuples->fillg75(float(trackIDhit),1.);
           //                                                                     Station I
           if (zz < z2) {
             //UserNtuples->fillg56(dx,1.);
             //UserNtuples->fillg57(dy,1.);
             //UserNtuples->fillg20(numofpart,1.);
             //UserNtuples->fillg21(SumEnerDeposit,1.);
             if (zside == 1) {
               //UserNtuples->fillg26(losenergy,1.);
             }
             if (zside == 2) {
               //UserNtuples->fillg76(losenergy,1.);
             }
             if (trackIDhit == 1) {
               //UserNtuples->fillg70(dx,1.);
               //UserNtuples->fillg71(incidentEnergyHit,1.);
               //UserNtuples->fillg79(losenergy,1.);
             } else {
               //UserNtuples->fillg82(dx,1.);
             }
           }
           //                                                                     Station II
           if (zz < z3 && zz > z2) {
             //UserNtuples->fillg52(dx,1.);
             //UserNtuples->fillg53(dy,1.);
             //UserNtuples->fillg22(numofpart,1.);
             //UserNtuples->fillg23(SumEnerDeposit,1.);
             //UserNtuples->fillg80(incidentEnergyHit,1.);
             //UserNtuples->fillg81(float(trackIDhit),1.);
             if (zside == 1) {
               //UserNtuples->fillg27(losenergy,1.);
             }
             if (zside == 2) {
               //UserNtuples->fillg77(losenergy,1.);
             }
             if (trackIDhit == 1) {
               //UserNtuples->fillg72(dx,1.);
               //UserNtuples->fillg73(incidentEnergyHit,1.);
             } else {
               //UserNtuples->fillg83(dx,1.);
             }
           }
           //                                                                     Station III
           if (zz < z4 && zz > z3) {
             //UserNtuples->fillg68(dx,1.);
             //UserNtuples->fillg69(dy,1.);
             //UserNtuples->filld210(xx,yy,1.);
             //UserNtuples->fillg22(numofpart,1.);
             //UserNtuples->fillg23(SumEnerDeposit,1.);
             if (zside == 1) {
               //UserNtuples->fillg28(losenergy,1.);
             }
             if (zside == 2) {
               //UserNtuples->fillg78(losenergy,1.);
             }
           }
         }  // MIonly or noMIonly ENDED
            //    }
 
         //     !!!!!!!!!!!!!
 
       }  // for loop on all hits ENDED  ENDED  ENDED  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
       //     !!!!!!!!!!!!!
 
       //======================================================================================================SUMHIT
       //UserNtuples->fillg29(totallosenergy,1.);
       //UserNtuples->fillg36(nhit11,1.);
       //UserNtuples->fillg37(nhit12,1.);
       //UserNtuples->fillg38(nhit13,1.);
       //======================================================================================================SUMHIT
       //   int rn00=3;//test only with 2 sensors in superlayer, not 4
       //      int rn00=rn0;//always
       if (verbosity > 2) {
         std::cout << "22222222222222222222222222222222222222222222222222222222222222222222222222 " << std::endl;
       }
       int det = 1;
       int allplacesforsensors = 7;
       for (int sector = 1; sector < sn0; sector++) {
         for (int zmodule = 1; zmodule < pn0; zmodule++) {
           for (int zsideinorder = 1; zsideinorder < allplacesforsensors; zsideinorder++) {
             int zside = FP420NumberingScheme::realzside(rn00, zsideinorder);  //1,3,5,2,4,6
             if (verbosity > 2) {
               std::cout << "FP420Test:  sector= " << sector << " zmodule= " << zmodule
                         << " zsideinorder= " << zsideinorder << " zside= " << zside << std::endl;
             }
             if (zside != 0) {
               int justlayer = FP420NumberingScheme::unpackLayerIndex(rn00, zside);  // 1,2
               if (justlayer < 1 || justlayer > 2) {
                 std::cout << "FP420Test:WRONG  justlayer= " << justlayer << std::endl;
               }
               int copyinlayer = FP420NumberingScheme::unpackCopyIndex(rn00, zside);  // 1,2,3
               if (copyinlayer < 1 || copyinlayer > 3) {
                 std::cout << "FP420Test:WRONG  copyinlayer= " << copyinlayer << std::endl;
               }
               int orientation = FP420NumberingScheme::unpackOrientation(rn00, zside);  // Front: = 1; Back: = 2
               if (orientation < 1 || orientation > 2) {
                 std::cout << "FP420Test:WRONG  orientation= " << orientation << std::endl;
               }
 
               // iu is a continues numbering of planes(!)  over two arm FP420 set up
               int detfixed =
                   1;  // use this treatment for each set up arm, hence no sense to do it defferently for +FP420 and -FP420;
               //                                                                    and  ...[ii] massives have prepared in such a way
               unsigned int ii =
                   FP420NumberingScheme::packMYIndex(rn00, pn0, sn0, detfixed, justlayer, sector, zmodule) - 1;
               // ii = 0-19   --> 20 items
               if (verbosity > 2) {
                 std::cout << "FP420Test:  justlayer = " << justlayer << " copyinlayer = " << copyinlayer
                           << " orientation = " << orientation << " ii= " << ii << std::endl;
               }
               double zdiststat = 0.;
               if (sn0 < 4) {
                 if (sector == 2)
                   zdiststat = zD3;
               } else {
                 if (sector == 2)
                   zdiststat = zD2;
                 if (sector == 3)
                   zdiststat = zD3;
               }
               double kplane = -(pn0 - 1) / 2 - 0.5 + (zmodule - 1);  //-3.5 +0...5 = -3.5,-2.5,-1.5,+2.5,+1.5
               double zcurrent = zinibeg + z420 + (ZSiStep - ZSiPlane) / 2 + kplane * ZSiStep + zdiststat;
               //double zcurrent = zinibeg +(ZSiStep-ZSiPlane)/2  + kplane*ZSiStep + (sector-1)*zUnit;
               if (verbosity > 2) {
                 std::cout << "FP420Test:  Leftzcurrent-419000. = " << zcurrent - 419000. << std::endl;
                 std::cout << "FP420Test:  ZGapLDet = " << ZGapLDet << std::endl;
               }
               if (justlayer == 1) {
                 if (orientation == 1)
                   zcurrent += (ZGapLDet + ZSiDet / 2);
                 if (orientation == 2)
                   zcurrent += zBlade - (ZGapLDet + ZSiDet / 2);
               }
               if (justlayer == 2) {
                 if (orientation == 1)
                   zcurrent += (ZGapLDet + ZSiDet / 2) + zBlade + gapBlade;
                 if (orientation == 2)
                   zcurrent += 2 * zBlade + gapBlade - (ZGapLDet + ZSiDet / 2);
               }
               //   .
               //
               if (det == 2)
                 zcurrent = -zcurrent;
               //
               if (verbosity > 2) {
                 std::cout << "FP420Test:  zcurrent-419000. = " << zcurrent - 419000. << std::endl;
               }
               //================================== end of for loops in continuius number iu:
             }  //if(zside!=0
           }    // for superlayer
         }      // for zmodule
       }        // for sector
 
       if (verbosity > 2) {
         std::cout << "----------------------------------------------------------------------------- " << std::endl;
       }
 
       //======================================================================================================CHECK
       if (totallosenergy == 0.0) {
         std::cout << "FP420Test:     number of hits = " << theCAFI->entries() << std::endl;
         for (int j = 0; j < theCAFI->entries(); j++) {
           FP420G4Hit* aHit = (*theCAFI)[j];
           double losenergy = aHit->getEnergyLoss();
           std::cout << " j hits = " << j << "losenergy = " << losenergy << std::endl;
         }
       }
       //======================================================================================================CHECK
 
       //====================================================================================================== HIT  START
 
       //   FIBRE Hit collected analysis
       double totalEnergy = 0.;
       int nhitsX = 0, nhitsY = 0, nsumhit = 0;
       for (int sector = 1; sector < 4; sector++) {
         int nhitsecX = 0, nhitsecY = 0;
         for (int zmodule = 1; zmodule < 11; zmodule++) {
           for (int zside = 1; zside < 3; zside++) {
             int det = 1;
             //      int nhit = 0;
             //  int sScale = 20;
             int index = FP420NumberingScheme::packFP420Index(det, zside, sector, zmodule);
             double theTotalEnergy = themap[index];
             //   X planes
             if (zside < 2) {
               //UserNtuples->fillg47(theTotalEnergy,1.);
               if (theTotalEnergy > 0.00003) {
                 nhitsX += 1;
                 //      nhitsecX += themap1[index];
                 //      nhit=1;
               }
             }
             //   Y planes
             else {
               //UserNtuples->fillg49(theTotalEnergy,1.);
               if (theTotalEnergy > 0.00005) {
                 nhitsY += 1;
                 //      nhitsecY += themap1[index];
                 //      nhit=1;
               }
             }
             // intindex is a continues numbering of FP420
             //        int zScale=2;  unsigned int intindex = sScale*(sector - 1)+zScale*(zmodule - 1)+zside;
             // int zScale=10;       unsigned int intindex = sScale*(sector - 1)+zScale*(zside - 1)+zmodule;
             //UserNtuples->fillg41(float(intindex),theTotalEnergy);
             //UserNtuples->fillg42(float(intindex),1.);
             //UserNtuples->fillp208(float(intindex),float(nhit),1.);
             //UserNtuples->fillp211(float(intindex),float(themap1[index]),1.);
             totalEnergy += themap[index];
           }  // for
         }    // for
         //UserNtuples->fillg39(nhitsecY,1.);
         if (nhitsecX > 10 || nhitsecY > 10) {
           nsumhit += 1;
           //UserNtuples->fillp213(float(sector),float(1.),1.);
         } else {  //UserNtuples->fillp213(float(sector),float(0.),1.);
         }
       }  // for
       //====================================================================================================== HIT  END
 
       //====================================================================================================== HIT  ALL
       //UserNtuples->fillg43(totalEnergy,1.);
       //UserNtuples->fillg58(nhitsX,1.);
       //UserNtuples->fillg59(nhitsY,1.);
       //  if( nsumhit !=0 ) { //UserNtuples->fillp212(vy,float(1.),1.);
       if (nsumhit >= 2) {  //UserNtuples->fillp212(vy,float(1.),1.);
       } else {             //UserNtuples->fillp212(vy,float(0.),1.);
       }
 
       //====================================================================================================== HIT  ALL
 
       //====================================================================================================== number of hits
       // .............
       //    } // number of hits < 50
       // .............
     }  // MI or no MI or all  - end
 
   }  // primary end
      //=========================== thePrim != 0  end   ===
   // ==========================================================================
   if (verbosity > 0) {
     std::cout << "FP420Test:  END OF Event " << (*evt)()->GetEventID() << std::endl;
   }
 }

Member Data Documentation

int FP420Test::dn0

private

Definition at line 205 of file FP420Test.h.

G4double FP420Test::entot0

private

Definition at line 183 of file FP420Test.h.

std::string FP420Test::fDataLabel

private

Definition at line 218 of file FP420Test.h.

std::string FP420Test::fOutputFile

private

Definition at line 219 of file FP420Test.h.

Float_t FP420Test::fp420eventarray[1]

private

Definition at line 212 of file FP420Test.h.

TNtuple* FP420Test::fp420eventntuple

private

Definition at line 213 of file FP420Test.h.

TFile FP420Test::fp420OutputFile

private

Definition at line 214 of file FP420Test.h.

std::string FP420Test::fRecreateFile

private

Definition at line 220 of file FP420Test.h.

double FP420Test::gapBlade

private

Definition at line 209 of file FP420Test.h.

int FP420Test::iev

private

Definition at line 181 of file FP420Test.h.

int FP420Test::itrk

private

Definition at line 182 of file FP420Test.h.

G4ThreeVector FP420Test::lastpo

private

Definition at line 201 of file FP420Test.h.

int FP420Test::lastTrackID

private

Definition at line 193 of file FP420Test.h.

int FP420Test::numofpart

private

Definition at line 199 of file FP420Test.h.

int FP420Test::pn0

private

Definition at line 205 of file FP420Test.h.

double FP420Test::rinCalo

private

Definition at line 192 of file FP420Test.h.

int FP420Test::rn0

private

Definition at line 205 of file FP420Test.h.

int FP420Test::rn00

private

Definition at line 206 of file FP420Test.h.

int FP420Test::sn0

private

Definition at line 205 of file FP420Test.h.

G4double FP420Test::SumEnerDeposit

private

Definition at line 197 of file FP420Test.h.

G4double FP420Test::SumStepc

private

Definition at line 197 of file FP420Test.h.

G4double FP420Test::SumStepl

private

Definition at line 197 of file FP420Test.h.

Fp420AnalysisHistManager* FP420Test::TheHistManager

private

Definition at line 217 of file FP420Test.h.

G4double FP420Test::tracklength0

private

Definition at line 183 of file FP420Test.h.

int FP420Test::verbosity

private

Definition at line 194 of file FP420Test.h.

int FP420Test::whichevent

private

Definition at line 215 of file FP420Test.h.

double FP420Test::z1

private

Definition at line 204 of file FP420Test.h.

double FP420Test::z2

private

Definition at line 204 of file FP420Test.h.

double FP420Test::z3

private

Definition at line 204 of file FP420Test.h.

double FP420Test::z4

private

Definition at line 204 of file FP420Test.h.

double FP420Test::z420

private

Definition at line 207 of file FP420Test.h.

double FP420Test::zBlade

private

Definition at line 209 of file FP420Test.h.

double FP420Test::ZBoundDet

private

Definition at line 208 of file FP420Test.h.

double FP420Test::zD2

private

Definition at line 204 of file FP420Test.h.

double FP420Test::zD3

private

Definition at line 204 of file FP420Test.h.

double FP420Test::ZGapLDet

private

Definition at line 208 of file FP420Test.h.

double FP420Test::zinCalo

private

Definition at line 192 of file FP420Test.h.

double FP420Test::zinibeg

private

Definition at line 208 of file FP420Test.h.

double FP420Test::ZSiDet

private

Definition at line 207 of file FP420Test.h.

double FP420Test::ZSiPlane

private

Definition at line 208 of file FP420Test.h.

double FP420Test::ZSiStep

private

Definition at line 208 of file FP420Test.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation