Fork me on GitHub

Changeset c61b5ce in git


Ignore:
Timestamp:
Dec 10, 2021, 2:11:26 AM (3 years ago)
Author:
christinaw97 <christina.wang@…>
Children:
3c59f98
Parents:
a09b75f
Message:

cleaning

Files:
4 edited

Legend:

Unmodified
Added
Removed
  • Makefile

    ra09b75f rc61b5ce  
    523523        classes/SortableObject.h
    524524tmp/classes/DelphesCscClusterFormula.$(ObjSuf): \
    525                 classes/DelphesCscClusterFormula.$(SrcSuf) \
    526                 classes/DelphesCscClusterFormula.h
     525        classes/DelphesCscClusterFormula.$(SrcSuf) \
     526        classes/DelphesCscClusterFormula.h
    527527tmp/classes/DelphesCylindricalFormula.$(ObjSuf): \
    528528        classes/DelphesCylindricalFormula.$(SrcSuf) \
     
    762762        modules/CscClusterId.h \
    763763        classes/DelphesClasses.h \
    764   classes/DelphesFactory.h \
     764        classes/DelphesFactory.h \
    765765        classes/DelphesCscClusterFormula.h \
    766766        external/ExRootAnalysis/ExRootClassifier.h \
  • cards/delphes_card_CMS_CSCCluster.tcl

    ra09b75f rc61b5ce  
    55set ExecutionPath {
    66  ParticlePropagator
    7 
    8 
    97
    108  ChargedHadronTrackingEfficiency
     
    6159  CutBasedIDEfficiency
    6260  ClusterEfficiency
    63 
    6461
    6562  TreeWriter
     
    414411  set InputArray Delphes/allParticles
    415412  set OutputArray LLP
    416   set DecayRegion 1
    417413  # DecayRegion = 0: no cuts on decay region
    418414  # DecayRegion = 1: select LLP that decays in CSC volume
    419415  # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
     416  set DecayRegion 1
    420417  set RequireStatus false
    421418  add PdgCode {1500001}
     
    427424  set InputArray Delphes/allParticles
    428425  set OutputArray LLP
    429   set DecayRegion 0
    430   set RequireStatus false
    431 
    432426  # DecayRegion = 0: no cuts on decay region
    433427  # DecayRegion = 1: select LLP that decays in CSC volume
    434428  # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
     429  set DecayRegion 0
     430  set RequireStatus false
    435431  add PdgCode {1500001}
    436432
     
    10341030module TreeWriter TreeWriter {
    10351031# 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
    10491033
    10501034  add Branch UniqueObjectFinder/jets Jet Jet
    10511035  add Branch UniqueObjectFinder/electrons Electron Electron
    1052   #add Branch UniqueObjectFinder/photons Photon Photon
    10531036  add Branch UniqueObjectFinder/muons Muon Muon
    10541037
    1055   #add Branch FatJetFinder/jets FatJet Jet
    1056 
    10571038  add Branch MissingET/momentum MissingET MissingET
    1058   #add Branch ScalarHT/energy ScalarHT ScalarHT
    10591039  add Branch llpFilter/LLP llp CscCluster
    10601040  add Branch CSCFilter/LLP Cscllp CscCluster
    1061 
    10621041  add Branch ClusterEfficiency/cluster CscCluster CscCluster
    10631042}
  • classes/DelphesClasses.h

    ra09b75f rc61b5ce  
    163163  Float_t decayT;
    164164
    165   Float_t ctau;
    166 
    167165  static CompBase *fgCompare; //!
    168166  const CompBase *GetCompare() const { return fgCompare; }
     
    650648{
    651649public:
    652   Int_t NHits; //nCSC hits
    653650  Float_t Eta; // eta of LLP
    654651  Float_t Phi; // phi of LLP
     
    660657  Float_t Ehad; // had energy of LLP
    661658  Float_t Eem; // em energy of LLP
    662 
    663659  Float_t pid; // LLP pid
    664 
    665660  Float_t T; // LLP decay time-photon travel time
    666661  Float_t X; // LLP decay x
     
    668663  Float_t Z; //  LLP decay z
    669664  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
    672667
    673668
  • modules/TreeWriter.cc

    ra09b75f rc61b5ce  
    245245    entry->decayZ = DecayPosition.Z();
    246246    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 millimeter
    250247  }
    251248}
     
    930927    entry->Phi = momentum.Phi();
    931928
    932 
    933     // entry->Eta = position.Eta();
    934     // entry->Phi = position.Phi();
    935 
    936929    entry->PT = momentum.Pt(); // pt of LLP
    937930    entry->Px = momentum.Px();// px of LLP
     
    940933    entry->E = momentum.E(); // E of LLP
    941934    entry->pid = candidate->PID; // LLP pid
    942     entry->Eem = candidate->Eem; // LLP pid
    943     entry->Ehad = candidate->Ehad; // LLP pid
     935    entry->Eem = candidate->Eem; // LLP Eem
     936    entry->Ehad = candidate->Ehad; // LLP Ehad
    944937    Double_t beta = momentum.P()/momentum.E();
    945938    Double_t gamma = 1.0/sqrt(1-beta*beta);
    946939    Double_t decayDistance = sqrt(pow(position.X(),2)+pow(position.Y(),2)+pow(position.Z(),2)); // mm
    947940    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
    952942    entry->T = decayDistance*(1./beta-1)* 1.0E-3/c_light*1e9; // ns
    953943    entry->X = position.X(); // LLP decay x
    954944    entry->Y = position.Y(); //  LLP decay y
    955945    entry->Z = position.Z(); //  LLP decay z
    956     // entry->Size = 100;
    957946  }
    958947}
Note: See TracChangeset for help on using the changeset viewer.