|
|
|
#include <FP420Test.h>
Public Member Functions | |
| FP420Test (const edm::ParameterSet &p) | |
| virtual | ~FP420Test () |
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 BeginOfEvent *evt) |
| This routine will be called when the appropriate signal arrives. | |
| void | update (const EndOfEvent *evt) |
| This routine will be called when the appropriate signal arrives. | |
| void | update (const EndOfTrack *trk) |
| This routine will be called when the appropriate signal arrives. | |
| void | update (const G4Step *step) |
| This routine will be called when the appropriate signal arrives. | |
| void | update (const BeginOfJob *run) |
| This routine will be called when the appropriate signal arrives. | |
| void | update (const BeginOfRun *run) |
| This routine will be called when the appropriate signal arrives. | |
| void | update (const EndOfRun *run) |
| This routine will be called when the appropriate signal arrives. | |
| void | update (const BeginOfTrack *trk) |
| This routine will be called when the appropriate signal arrives. | |
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 |
Definition at line 162 of file FP420Test.h.
| FP420Test::FP420Test | ( | const edm::ParameterSet & | p | ) |
Definition at line 64 of file FP420Test.cc.
References gather_cfg::cout, edm::ParameterSet::getParameter(), and 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 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;
}
| void FP420Test::detectorLevel | ( | const G4VTouchable * | touch, |
| int & | level, | ||
| int * | copyno, | ||
| G4String * | name | ||
| ) | const [private] |
Definition at line 923 of file FP420Test.cc.
References i, and testEve_cfg::level.
| 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.
| void FP420Test::update | ( | const BeginOfEvent * | ) | [private, virtual] |
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 verbosity.
{
iev = (*evt)()->GetEventID();
if (verbosity > 0) {
std::cout <<"FP420Test:update Event number = " << iev << std::endl;
}
whichevent++;
}
| void FP420Test::update | ( | const EndOfEvent * | ) | [private, virtual] |
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, getHLTprescales::index, j, funct::log(), npart, ntfp420_evt, FP420NumberingScheme::packFP420Index(), phi, pi, funct::tan(), FP420NumberingScheme::unpackFP420Index(), and verbosity.
{
// ==========================================================================
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;
}
}
| void FP420Test::update | ( | const EndOfTrack * | ) | [private, virtual] |
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 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 G4Step * | ) | [private, virtual] |
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 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 BeginOfJob * | ) | [private, virtual] |
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 * | ) | [private, virtual] |
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 * | ) | [private, virtual] |
This routine will be called when the appropriate signal arrives.
Implements Observer< const EndOfRun * >.
Definition at line 419 of file FP420Test.cc.
{;}
| void FP420Test::update | ( | const BeginOfTrack * | ) | [private, virtual] |
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 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.;
}
}
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.
Definition at line 192 of file FP420Test.h.
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.
1.7.3