Changeset c61b5ce in git
- Timestamp:
- Dec 10, 2021, 2:11:26 AM (3 years ago)
- Children:
- 3c59f98
- Parents:
- a09b75f
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
Makefile
ra09b75f rc61b5ce 523 523 classes/SortableObject.h 524 524 tmp/classes/DelphesCscClusterFormula.$(ObjSuf): \ 525 526 525 classes/DelphesCscClusterFormula.$(SrcSuf) \ 526 classes/DelphesCscClusterFormula.h 527 527 tmp/classes/DelphesCylindricalFormula.$(ObjSuf): \ 528 528 classes/DelphesCylindricalFormula.$(SrcSuf) \ … … 762 762 modules/CscClusterId.h \ 763 763 classes/DelphesClasses.h \ 764 764 classes/DelphesFactory.h \ 765 765 classes/DelphesCscClusterFormula.h \ 766 766 external/ExRootAnalysis/ExRootClassifier.h \ -
cards/delphes_card_CMS_CSCCluster.tcl
ra09b75f rc61b5ce 5 5 set ExecutionPath { 6 6 ParticlePropagator 7 8 9 7 10 8 ChargedHadronTrackingEfficiency … … 61 59 CutBasedIDEfficiency 62 60 ClusterEfficiency 63 64 61 65 62 TreeWriter … … 414 411 set InputArray Delphes/allParticles 415 412 set OutputArray LLP 416 set DecayRegion 1417 413 # DecayRegion = 0: no cuts on decay region 418 414 # DecayRegion = 1: select LLP that decays in CSC volume 419 415 # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation 416 set DecayRegion 1 420 417 set RequireStatus false 421 418 add PdgCode {1500001} … … 427 424 set InputArray Delphes/allParticles 428 425 set OutputArray LLP 429 set DecayRegion 0430 set RequireStatus false431 432 426 # DecayRegion = 0: no cuts on decay region 433 427 # DecayRegion = 1: select LLP that decays in CSC volume 434 428 # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation 429 set DecayRegion 0 430 set RequireStatus false 435 431 add PdgCode {1500001} 436 432 … … 1034 1030 module TreeWriter TreeWriter { 1035 1031 # add Branch InputArray BranchName BranchClass 1036 #add Branch Delphes/allParticles Particle GenParticle 1037 #add Branch Delphes/stableParticles Particle GenParticle 1038 #add Branch TrackMerger/tracks Track Track 1039 #add Branch Calorimeter/towers Tower Tower 1040 1041 #add Branch HCal/eflowTracks EFlowTrack Track 1042 #add Branch ECal/eflowPhotons EFlowPhoton Tower 1043 #add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower 1044 1045 #add Branch GenJetFinder/jets GenJet Jet 1046 #add Branch GenMissingET/momentum GenMissingET MissingET 1047 1048 1032 # add Branch Delphes/allParticles Particle GenParticle 1049 1033 1050 1034 add Branch UniqueObjectFinder/jets Jet Jet 1051 1035 add Branch UniqueObjectFinder/electrons Electron Electron 1052 #add Branch UniqueObjectFinder/photons Photon Photon1053 1036 add Branch UniqueObjectFinder/muons Muon Muon 1054 1037 1055 #add Branch FatJetFinder/jets FatJet Jet1056 1057 1038 add Branch MissingET/momentum MissingET MissingET 1058 #add Branch ScalarHT/energy ScalarHT ScalarHT1059 1039 add Branch llpFilter/LLP llp CscCluster 1060 1040 add Branch CSCFilter/LLP Cscllp CscCluster 1061 1062 1041 add Branch ClusterEfficiency/cluster CscCluster CscCluster 1063 1042 } -
classes/DelphesClasses.h
ra09b75f rc61b5ce 163 163 Float_t decayT; 164 164 165 Float_t ctau;166 167 165 static CompBase *fgCompare; //! 168 166 const CompBase *GetCompare() const { return fgCompare; } … … 650 648 { 651 649 public: 652 Int_t NHits; //nCSC hits653 650 Float_t Eta; // eta of LLP 654 651 Float_t Phi; // phi of LLP … … 660 657 Float_t Ehad; // had energy of LLP 661 658 Float_t Eem; // em energy of LLP 662 663 659 Float_t pid; // LLP pid 664 665 660 Float_t T; // LLP decay time-photon travel time 666 661 Float_t X; // LLP decay x … … 668 663 Float_t Z; // LLP decay z 669 664 Float_t R; // LLP decay z 670 Float_t beta; 671 Float_t ctau; 665 Float_t beta; // LLP beta 666 Float_t ctau; //LLP ctau 672 667 673 668 -
modules/TreeWriter.cc
ra09b75f rc61b5ce 245 245 entry->decayZ = DecayPosition.Z(); 246 246 entry->decayT = DecayPosition.T()* 1.0E-3 / c_light; 247 float beta = entry->P/momentum.E();248 float gamma = 1./sqrt(1-beta*beta);249 entry->ctau = sqrt(pow(entry->decayX-entry->X,2)+pow(entry->decayY-entry->Y,2)+pow(entry->decayZ-entry->Z,2))/(beta*gamma);// in millimeter250 247 } 251 248 } … … 930 927 entry->Phi = momentum.Phi(); 931 928 932 933 // entry->Eta = position.Eta();934 // entry->Phi = position.Phi();935 936 929 entry->PT = momentum.Pt(); // pt of LLP 937 930 entry->Px = momentum.Px();// px of LLP … … 940 933 entry->E = momentum.E(); // E of LLP 941 934 entry->pid = candidate->PID; // LLP pid 942 entry->Eem = candidate->Eem; // LLP pid943 entry->Ehad = candidate->Ehad; // LLP pid935 entry->Eem = candidate->Eem; // LLP Eem 936 entry->Ehad = candidate->Ehad; // LLP Ehad 944 937 Double_t beta = momentum.P()/momentum.E(); 945 938 Double_t gamma = 1.0/sqrt(1-beta*beta); 946 939 Double_t decayDistance = sqrt(pow(position.X(),2)+pow(position.Y(),2)+pow(position.Z(),2)); // mm 947 940 entry->beta = beta; // LLP pid 948 entry->ctau = decayDistance/(beta * gamma);; // LLP pid 949 950 // entry->T = (position.T()-sqrt(pow(position.X(),2)+pow(position.Y(),2)+pow(position.Z(),2)))* 1.0E-3 / c_light; // LLP decay time-photon travel time 951 941 entry->ctau = decayDistance/(beta * gamma); // LLP travel time in its rest frame 952 942 entry->T = decayDistance*(1./beta-1)* 1.0E-3/c_light*1e9; // ns 953 943 entry->X = position.X(); // LLP decay x 954 944 entry->Y = position.Y(); // LLP decay y 955 945 entry->Z = position.Z(); // LLP decay z 956 // entry->Size = 100;957 946 } 958 947 }
Note:
See TracChangeset
for help on using the changeset viewer.