#include <FP420Test.h>

Inheritance diagram for FP420Test:

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

virtual	~FP420Test ()

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)
	This routine will be called when the appropriate signal arrives. More...

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

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

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

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

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

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

void	update (const EndOfEvent *evt)
	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

FP420NumberingScheme *	theFP420NumberingScheme

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 162 of file FP420Test.h.

Constructor & Destructor Documentation

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

Definition at line 64 of file FP420Test.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity.

                                             {
   //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:
 
   theFP420NumberingScheme = new FP420NumberingScheme();
   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 ( )

virtual

Definition at line 135 of file FP420Test.cc.

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

                       {
   //  delete UserNtuples;
   delete theFP420NumberingScheme;
 
   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::cout << std::endl; 
   }
 
   std::cout<<"FP420Test: End of process"<<std::endl;
 
 
 
 }

Member Function Documentation

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

private

Definition at line 923 of file FP420Test.cc.

References i, cuy::ii, and testEve_cfg::level.

                                                                          {
 
   //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 != 0) 
         name[ii] = pv->GetName();
       else
         name[ii] = "Unknown";
       copyno[ii] = touch->GetReplicaNumber(i);
     }
   }
 }

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

private

Definition at line 901 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 911 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 * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 405 of file FP420Test.cc.

References gather_cfg::cout.

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

void FP420Test::update ( const BeginOfRun * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 412 of file FP420Test.cc.

References gather_cfg::cout.

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

void FP420Test::update ( const EndOfRun * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfRun * >.

Definition at line 419 of file FP420Test.cc.

419 {;}

void FP420Test::update ( const BeginOfEvent * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 424 of file FP420Test.cc.

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

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

void FP420Test::update ( const BeginOfTrack * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 433 of file FP420Test.cc.

References gather_cfg::cout, and PFRecoTauDiscriminationAgainstMuon2_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 * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 560 of file FP420Test.cc.

References gather_cfg::cout, and PFRecoTauDiscriminationAgainstMuon2_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 == 0) {
      std::cout << "FP420Test on aStep: No trk info !!!! abort " << std::endl;
      //     throw Genexception("FP420Test:FP420Test on aStep: cannot get trkInfo");
    } 
   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();
   G4ThreeVector   vert_pos       = theTrack->GetVertexPosition(); // vertex ,where this track was created
   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(); 
   G4ThreeVector     preposition   = preStepPoint->GetPosition();    
   G4ThreeVector     prelocalpoint = theTrack->GetTouchable()->GetHistory()->
                                            GetTopTransform().TransformPoint(preposition);
   G4VPhysicalVolume* currentPV     = preStepPoint->GetPhysicalVolume();
   G4String         prename       = currentPV->GetName();
 
 const G4VTouchable*  pre_touch    = preStepPoint->GetTouchable();
      int          pre_levels   = detLevels(pre_touch);
 
 //     G4String      pre_name1    = detName(pre_touch, pre_levels, 1);
 //     G4String      pre_name2    = detName(pre_touch, pre_levels, 2);
 //     G4String      pre_name3    = detName(pre_touch, pre_levels, 3);
         G4String name1[20]; int copyno1[20];
       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 * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack * >.

Definition at line 450 of file FP420Test.cc.

References gather_cfg::cout, phi, and PFRecoTauDiscriminationAgainstMuon2_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.);
          }
 
      /*
       float th_tr     = lastpo.theta();
       float eta_tr    = -log(tan(th_tr/2));
       float phi_tr    = lastpo.phi();
       if (phi_tr < 0.) phi_tr += twopi;
 */
 
 
   }
 
 }

void FP420Test::update ( const EndOfEvent * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 942 of file FP420Test.cc.

References gather_cfg::cout, eta(), FP420G4Hit::getEnergyLoss(), FP420G4Hit::getEntry(), FP420G4Hit::getEntryLocalP(), FP420G4Hit::getExitLocalP(), FP420G4Hit::getTrackID(), FP420G4Hit::getUnitID(), i, cuy::ii, cmsHarvester::index, j, log, npart, ntfp420_evt, FP420NumberingScheme::packFP420Index(), phi, pi, funct::tan(), FP420NumberingScheme::unpackFP420Index(), PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity, and ecaldqm::zside().

                                              {
   // ==========================================================================
   
   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=0;
 
 
   // 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 == 0)
     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==0) thePrim=avertex->GetPrimary(trackID);
 
        if (thePrim!=0) {
      // 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 != 0) {
 //      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
     varia = 0;
     //                      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 = theFP420NumberingScheme->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 = theFP420NumberingScheme->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 = theFP420NumberingScheme->FP420NumberingScheme::unpackLayerIndex(rn00, zside);// 1,2
        if(justlayer<1||justlayer>2) {
          std::cout << "FP420Test:WRONG  justlayer= " << justlayer << std::endl; 
        }
        int copyinlayer = theFP420NumberingScheme->FP420NumberingScheme::unpackCopyIndex(rn00, zside);// 1,2,3
        if(copyinlayer<1||copyinlayer>3) {
          std::cout << "FP420Test:WRONG  copyinlayer= " << copyinlayer << std::endl; 
        }
        int orientation = theFP420NumberingScheme->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=theFP420NumberingScheme->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 221 of file FP420Test.h.

G4double FP420Test::entot0

private

Definition at line 197 of file FP420Test.h.

std::string FP420Test::fDataLabel

private

Definition at line 235 of file FP420Test.h.

std::string FP420Test::fOutputFile

private

Definition at line 236 of file FP420Test.h.

Float_t FP420Test::fp420eventarray[1]

private

Definition at line 229 of file FP420Test.h.

TNtuple* FP420Test::fp420eventntuple

private

Definition at line 230 of file FP420Test.h.

TFile FP420Test::fp420OutputFile

private

Definition at line 231 of file FP420Test.h.

std::string FP420Test::fRecreateFile

private

Definition at line 237 of file FP420Test.h.

double FP420Test::gapBlade

private

Definition at line 225 of file FP420Test.h.

int FP420Test::iev

private

Definition at line 195 of file FP420Test.h.

int FP420Test::itrk

private

Definition at line 196 of file FP420Test.h.

G4ThreeVector FP420Test::lastpo

private

Definition at line 216 of file FP420Test.h.

int FP420Test::lastTrackID

private

Definition at line 208 of file FP420Test.h.

int FP420Test::numofpart

private

Definition at line 214 of file FP420Test.h.

int FP420Test::pn0

private

Definition at line 221 of file FP420Test.h.

double FP420Test::rinCalo

private

Definition at line 207 of file FP420Test.h.

int FP420Test::rn0

private

Definition at line 221 of file FP420Test.h.

int FP420Test::rn00

private

Definition at line 222 of file FP420Test.h.

int FP420Test::sn0

private

Definition at line 221 of file FP420Test.h.

G4double FP420Test::SumEnerDeposit

private

Definition at line 212 of file FP420Test.h.

G4double FP420Test::SumStepc

private

Definition at line 212 of file FP420Test.h.

G4double FP420Test::SumStepl

private

Definition at line 212 of file FP420Test.h.

FP420NumberingScheme* FP420Test::theFP420NumberingScheme

private

Definition at line 192 of file FP420Test.h.

Fp420AnalysisHistManager* FP420Test::TheHistManager

private

Definition at line 234 of file FP420Test.h.

G4double FP420Test::tracklength0

private

Definition at line 197 of file FP420Test.h.

int FP420Test::verbosity

private

Definition at line 209 of file FP420Test.h.

int FP420Test::whichevent

private

Definition at line 232 of file FP420Test.h.

double FP420Test::z1

private

Definition at line 220 of file FP420Test.h.

double FP420Test::z2

private

Definition at line 220 of file FP420Test.h.

double FP420Test::z3

private

Definition at line 220 of file FP420Test.h.

double FP420Test::z4

private

Definition at line 220 of file FP420Test.h.

double FP420Test::z420

private

Definition at line 223 of file FP420Test.h.

double FP420Test::zBlade

private

Definition at line 225 of file FP420Test.h.

double FP420Test::ZBoundDet

private

Definition at line 224 of file FP420Test.h.

double FP420Test::zD2

private

Definition at line 220 of file FP420Test.h.

double FP420Test::zD3

private

Definition at line 220 of file FP420Test.h.

double FP420Test::ZGapLDet

private

Definition at line 224 of file FP420Test.h.

double FP420Test::zinCalo

private

Definition at line 207 of file FP420Test.h.

double FP420Test::zinibeg

private

Definition at line 224 of file FP420Test.h.

double FP420Test::ZSiDet

private

Definition at line 223 of file FP420Test.h.

double FP420Test::ZSiPlane

private

Definition at line 224 of file FP420Test.h.

double FP420Test::ZSiStep

private

Definition at line 224 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